The direct effects of tacrolimus and cyclosporin A on isolated human islets: A functional, survival and gene expression study.

The use of immunosuppressive drugs in transplanted patients is associated with the development of diabetes, possibly due to β-cell toxicity. To better understand the mechanisms leading to post-transplant diabetes, we investigated the actions of prolonged exposure of isolated human islets to therapeutical levels of tacrolimus (Tac) or cyclosporin A (CsA). Islets were isolated from the pancreas of multiorgan donors by enzymatic digestion and density gradient centrifugation. Functional, survival and molecular studies were then performed after 4 days of incubation with therapeutical concentrations of Tac or  CsA. Glucose-induced insulin secretion was significantly decreased in Tac, but not in CsA exposed islets, which was associated with a reduction of the amount of insulin granules as shown by electron microscopy. The percentage of apoptotic β-cells was higher in Tac than CsA exposed islets. Microarray experiments followed by Gene Set Enrichment Analysis revealed that gene expression was more markedly affected upon Tac treatment. In conclusion, Tac and CsA affect features of beta-cell differently, with several changes occurring at the molecular level.

Characterization and functional identification of a novel plant extradiol 4,5-dioxygenase involved in betalain pigment biosynthesis in Portulaca Grandiflora

RESUME Les bétalaïnes sont des pigments chromo-alcaloïdes violets et jaunes présents dans les plantes appartenant à l'ordre des Caryophyllales et dans les champignons des genres Amanita et Hygrocybe. Leur courte voie de biosynthèse est élucidée chimiquement depuis de nombreuses années, mais les enzymes impliquées dans cette biosynthèse chez les plantes ne sont toujours pas caractérisées. L'enzyme de la DOPA-dioxygénase d' Amanita muscaria a été identifiée (Girod et Zryd, 1991a), mais de nombreuses tentatives d'isolation d'un homologue chez les plantes ont échoué. Afin d'isoler les gènes spécifiques des bétalaïnes chez les plantes, nous avons construit des banques soustraites d'ADNc à partir d'ARN total de pétales immatures de Portulaca grandiflora (Pg) de génotypes jaunes et blancs, respectivement violets et blancs. Les clones couleur- spécifiques ont été détectés en premier par analyse Northem du RNA de pétales blancs et colorés. Les candidats positifs ont alors été soumis à une analyse de transcription au niveau des tiges colorées, vertes et des feuilles, afin d'établir leur expression spécifique. Deux ARNs messagers complets ont une expression corrélée avec l'accumulation des bétalaïnes dans les tissus. Le premier de ces clones, A.16, code pour une oxydase de l'acyl-Coenzyme A (ACX) putative, mais le domaine de liaison du FAD essentiel pour l'activité d'ACX est absent. Toutes nos tentatives pour démontrer sa fonction ont échoué. Le rôle de cette protéine dans la voie de synthèse des bétalaïnes reste inconnu. Le deuxième de ces clones spécifique aux bétalaïnes, L.6 (isolé par Zaiko, 2000), a été renommé DODA en raison de son homologie avec le domaine LigB (pfam02900) d'une 4,5-dioxygénase extradiol bactérienne. DODA a été identifié in silico comme une dioxygénase extradiol en raison de la conservation stricte, au niveau de sa séquence peptidique, des résidus catalytiques de LigB et de ceux liant le cofacteur fer. Une analyse de transfert Southem a montré que ce gène est unique dans Pg. L'expression transitoire de DODA par transformation biolistique dans des pétales blancs de Pg a produit des taches violettes ou jaunes dans des cellules transformées. Une analyse HPLC de ces taches a démontré leur identité avec les bétalaïnes présentes naturellement dans les pétales violets et jaunes de Pg, confirmant ainsi la complémentation par le gène Pg DODA de l'allèle récessif cc présent dans les pétales blancs de Pg. Des homologues de DODA (DOPA-dioxygénase) ont été identifiés dans de nombreuses espèces de plantes, y compris dans celles sans bétalaïne. L'alignement de ces homologues a permis l'identification d'un motif spécifique aux bétalaïnes à côté d'une histidine catalytique conservée. Ce motif [H-P-(S,A)-(N,D)-x-T-P] remplace le motif [H-N-L-R] conservé dans les plantes sans bétalaïne et le motif [H-N-L-x] présent dans tous les homologues bactériens et archaebactériens. Une modélisation tridimensionnelle préliminaire du site actif de Pg DODA et de son homologue dans la mousse Physcomitrella patens a montré l'importance de ce motif spécifique aux bétalaïnes pour l'accessibilité du substrat au site actif. L'analyse phylogénétique de DODA a confirmé l'évolution séparée de cette protéine chez les plantes à bétalaïnes par comparaison avec celle des plantes sans bétalaïne. Nous avons donc conclu que les bétalaïnes sont apparues par modification de l'affinité pour un substrat d'enzymes similaires à DODA, chez un ancêtre unique des Caryophyllales qui a perdu toute capacité de biosynthèse des anthocyanes. Finalement, Pg DODA n'a aucune similarité avec la protéine DODA d' Amanita muscaria, bien que celle-ci complémente aussi la pigmentation des pétales blancs de Pg. La biosynthèse des bétalaïnes est un exemple remarquable de convergence évolutive biochimique indépendante entre espèces de règnes différents. ABSTRACT Betalains are violet and yellow chromo-alkaloid pigments present in plants belonging to the order Caryophyllales and also in the fungal genera Amanita and Hygrocybe. Their short biosynthetic pathway is chemically well understood since many years, but enzymes involved in the plant pathway are still uncharacterized. The DOPA-dioxygenase from Amanita muscaria was identified (Girod and Zryd, 1991a), but numerous attempts to identify a plant homologue to the corresponding gene, failed. In order to isolate betalain-specific genes in plants, subtractive cDNA libraries were built with total RNA from white and yellow and respectively, violet immature petals from Portulaca grandiflora (Pg) genotypes. Colour-specific clones were first detected by Northern blot analysis using RNA from white and coloured petals. Positive candidates were submitted to further transcription analysis in coloured, green stems and leaves in order to assess their specific expression. Two full-length mRNAs showed a correlated expression with betalain accumulation in tissues. One of them, A.16, encodes a putative acyl-Coenzyme A oxidase (ACX), but missing the FAD binding domain essential for the ACX activity. Thus, all attempts to demonstrate its function failed. The role of this protein in the betalain biosynthesis pathway, if any, is still unknown. The second betalain-specific mRNA, L.6 (isolated by Zaiko, 2000) shows a homology with a LigB domain (pfam02900) from a bacterial extradiol 4,5-dioxygenase. It was then renamed DODA (DOPA-dioxygenase). DODA was identified in silico as a highly conserved extradiol dioxygenase due to the strict conservation of its peptidic sequence with LigB catalytic residues and iron-binding cofactor residues. Southern blot analysis showed that this gene is a single copy-gene in Pg. Transient expression of DODA protein through biolistic transformation of Pg white petals produced violet or yellow spots in individual cells. HPLC analysis of these spots showed an identity with betalain pigments present naturally in yellow and violet Pg petals, thus confirming the complementation of the recessive cc allele present in Pg white petals by Pg DODA gene. DODA homologues were identified in numerous plant species including those without betalain. Alignment of these homologues allowed the identification of a betalain-specific pattern beside a highly conserved catalytic histidine. This [H-P-(S,A)-(N,D)-x-T-P] pattern replaces a [H-N-L-R] pattern strictly conserved in non-betalain plants and a [H-N-L-x] pattern present in all bacterial and archaebacterial homologues. Preliminary three-dimensional modeling of the active site of Pg DODA and its Physcomitrella patens moss homologue revealed the importance of this betalain-specific pattern for the substrate accessibility to the DODA active site. DODA phylogenetic analysis confirmed the separate evolution of this protein in betalain-producing plants. We conclude that betalain pigments appeared in a unique ancestor of the Caryophyllales order in which anthocyanin biosynthetic pathway was impaired, by a modification of enzymes of the DODA family for substrate affinity. The Pg DODA protein has no sequence similarity with Amanita muscaria DODA, despite the fact that they both complement Pg white petals for their pigmentation. Betalain biosynthesis is an interesting example of independent biochemical evolutionary convergence between species from different kingdoms.

Functional error modeling for uncertainty quantification in hydrogeology

Approximate models (proxies) can be employed to reduce the computational costs of estimating uncertainty. The price to pay is that the approximations introduced by the proxy model can lead to a biased estimation. To avoid this problem and ensure a reliable uncertainty quantification, we propose to combine functional data analysis and machine learning to build error models that allow us to obtain an accurate prediction of the exact response without solving the exact model for all realizations. We build the relationship between proxy and exact model on a learning set of geostatistical realizations for which both exact and approximate solvers are run. Functional principal components analysis (FPCA) is used to investigate the variability in the two sets of curves and reduce the dimensionality of the problem while maximizing the retained information. Once obtained, the error model can be used to predict the exact response of any realization on the basis of the sole proxy response. This methodology is purpose-oriented as the error model is constructed directly for the quantity of interest, rather than for the state of the system. Also, the dimensionality reduction performed by FPCA allows a diagnostic of the quality of the error model to assess the informativeness of the learning set and the fidelity of the proxy to the exact model. The possibility of obtaining a prediction of the exact response for any newly generated realization suggests that the methodology can be effectively used beyond the context of uncertainty quantification, in particular for Bayesian inference and optimization.

Genome-wide meta-analysis identifies six novel loci associated with habitual coffee consumption

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91 462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10-8).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee.

Refractory intraretinal or subretinal fluid in neovascular age-related macular degeneration treated with intravitreal ranizubimab: Functional and Structural Outcome.

PURPOSE: To investigate the visual acuity results of eyes with neovascular age-related macular degeneration and refractory fluid despite monthly treatment with ranibizumab, and to investigate differences between refractory subretinal fluid and intraretinal cystic changes. METHODS: Retrospective chart review of consecutive treatment-refractory neovascular age-related macular degeneration, defined as persistent intraretinal or subretinal fluid despite monthly ranibizumab injections during 12 months or more. Data were evaluated for baseline characteristics, type and location of the refractory fluid, mean visual acuity change, number of injections, and the time point of first complete disappearance of all fluid on spectral domain optical coherence tomography. RESULTS: Seventy-six eyes (74 patients, mean age, 76.8 years) were identified. The mean follow-up was 33.6 months (range, 12-73 months). The mean number of injections was 11.4 in the first year and 27.7 over follow-up. The refractory fluid was located subfoveally in 61.8%. In 27 eyes (35.5%), the fluid resolved after a mean of 21.8 months (range, 13-49 months). Mean visual acuity increased by 9.0, 7.9, and 7.9 letters by Month 12, Month 24, and Month 36, respectively. Subgroup analysis revealed a higher risk for fibrosis (odds ratio, 3.30) or atrophy (odds ratio, 3.34) in patients with refractory cysts as compared with refractory subretinal fluid. Furthermore, refractory cysts showed a higher risk for a 10-letter visual acuity loss (P = 0.018). CONCLUSION: Fluid refractory to monthly treatment with ranibizumab for neovascular age-related macular degeneration still allowed for well-maintained visual improvement, even in subfoveal location. Late fluid resolution may occur. However, refractory cysts were associated with poorer anatomical and functional outcome than subretinal fluid.

Functional importance of cardiac enhancer-associated noncoding RNAs in heart development and disease.

The key information processing units within gene regulatory networks are enhancers. Enhancer activity is associated with the production of tissue-specific noncoding RNAs, yet the existence of such transcripts during cardiac development has not been established. Using an integrated genomic approach, we demonstrate that fetal cardiac enhancers generate long noncoding RNAs (lncRNAs) during cardiac differentiation and morphogenesis. Enhancer expression correlates with the emergence of active enhancer chromatin states, the initiation of RNA polymerase II at enhancer loci and expression of target genes. Orthologous human sequences are also transcribed in fetal human hearts and cardiac progenitor cells. Through a systematic bioinformatic analysis, we identified and characterized, for the first time, a catalog of lncRNAs that are expressed during embryonic stem cell differentiation into cardiomyocytes and associated with active cardiac enhancer sequences. RNA-sequencing demonstrates that many of these transcripts are polyadenylated, multi-exonic long noncoding RNAs. Moreover, knockdown of two enhancer-associated lncRNAs resulted in the specific downregulation of their predicted target genes. Interestingly, the reactivation of the fetal gene program, a hallmark of the stress response in the adult heart, is accompanied by increased expression of fetal cardiac enhancer transcripts. Altogether, these findings demonstrate that the activity of cardiac enhancers and expression of their target genes are associated with the production of enhancer-derived lncRNAs.

High doses of onabotulinumtoxinA in post-stroke spasticity: a retrospective analysis.

We retrospectively evaluated the efficacy and safety of high doses of onabotulinumtoxinA (from 600 to 800 units) in 26 patients affected by upper and/or lower limb post-stroke spasticity. They were assessed before, 30 and 90 days after treatment. We observed a significant muscle tone reduction and a significant functional improvement (assessed with the Disability Assessment Scale). No adverse events were reported. In our retrospective analysis the treatment with high doses of onabotulinumtoxinA showed to be effective and safe.

Embolic Strokes of Undetermined Source in the Athens Stroke Registry: an Outcome Analysis.

BACKGROUND AND PURPOSE: Information about outcomes in Embolic Stroke of Undetermined Source (ESUS) patients is unavailable. This study provides a detailed analysis of outcomes of a large ESUS population. METHODS: Data set was derived from the Athens Stroke Registry. ESUS was defined according to the Cryptogenic Stroke/ESUS International Working Group criteria. End points were mortality, stroke recurrence, functional outcome, and a composite cardiovascular end point comprising recurrent stroke, myocardial infarction, aortic aneurysm rupture, systemic embolism, or sudden cardiac death. We performed Kaplan-Meier analyses to estimate cumulative probabilities of outcomes by stroke type and Cox-regression to investigate whether stroke type was outcome predictor. RESULTS: 2731 patients were followed-up for a mean of 30.5±24.1months. There were 73 (26.5%) deaths, 60 (21.8%) recurrences, and 78 (28.4%) composite cardiovascular end points in the 275 ESUS patients. The cumulative probability of survival in ESUS was 65.6% (95% confidence intervals [CI], 58.9%-72.2%), significantly higher compared with cardioembolic stroke (38.8%, 95% CI, 34.9%-42.7%). The cumulative probability of stroke recurrence in ESUS was 29.0% (95% CI, 22.3%-35.7%), similar to cardioembolic strokes (26.8%, 95% CI, 22.1%-31.5%), but significantly higher compared with all types of noncardioembolic stroke. One hundred seventy-two (62.5%) ESUS patients had favorable functional outcome compared with 280 (32.2%) in cardioembolic and 303 (60.9%) in large-artery atherosclerotic. ESUS patients had similar risk of composite cardiovascular end point as all other stroke types, with the exception of lacunar strokes, which had significantly lower risk (adjusted hazard ratio, 0.70 [95% CI, 0.52-0.94]). CONCLUSIONS: Long-term mortality risk in ESUS is lower compared with cardioembolic strokes, despite similar rates of recurrence and composite cardiovascular end point. Recurrent stroke risk is higher in ESUS than in noncardioembolic strokes.

Multiplicative interaction of functional inflammasome genetic variants in determining the risk of gout.

INTRODUCTION: The acute gout flare results from a localised self-limiting innate immune response to monosodium urate (MSU) crystals deposited in joints in hyperuricaemic individuals. Activation of the caspase recruitment domain-containing protein 8 (CARD8) NOD-like receptor pyrin-containing 3 (NLRP3) inflammasome by MSU crystals and production of mature interleukin-1β (IL-1β) is central to acute gouty arthritis. However very little is known about genetic control of the innate immune response involved in acute gouty arthritis. Therefore our aim was to test functional single nucleotide polymorphism (SNP) variants in the toll-like receptor (TLR)-inflammasome-IL-1β axis for association with gout. METHODS: 1,494 gout cases of European and 863 gout cases of New Zealand (NZ) Polynesian (Māori and Pacific Island) ancestry were included. Gout was diagnosed by the 1977 ARA gout classification criteria. There were 1,030 Polynesian controls and 10,942 European controls including from the publicly-available Atherosclerosis Risk in Communities (ARIC) and Framingham Heart (FHS) studies. The ten SNPs were either genotyped by Sequenom MassArray or by Affymetrix SNP array or imputed in the ARIC and FHS datasets. Allelic association was done by logistic regression adjusting by age and sex with European and Polynesian data combined by meta-analysis. Sample sets were pooled for multiplicative interaction analysis, which was also adjusted by sample set. RESULTS: Eleven SNPs were tested in the TLR2, CD14, IL1B, CARD8, NLRP3, MYD88, P2RX7, DAPK1 and TNXIP genes. Nominally significant (P < 0.05) associations with gout were detected at CARD8 rs2043211 (OR = 1.12, P = 0.007), IL1B rs1143623 (OR = 1.10, P = 0.020) and CD14 rs2569190 (OR = 1.08; P = 0.036). There was significant multiplicative interaction between CARD8 and IL1B (P = 0.005), with the IL1B risk genotype amplifying the risk effect of CARD8. CONCLUSION: There is evidence for association of gout with functional variants in CARD8, IL1B and CD14. The gout-associated allele of IL1B increases expression of IL-1β - the multiplicative interaction with CARD8 would be consistent with a synergy of greater inflammasome activity (resulting from reduced CARD8) combined with higher levels of pre-IL-1β expression leading to increased production of mature IL-1β in gout.

Measuring individuals' response quality in self-administered psychological tests : an introduction to Gendre's functional method

The functional method is a new test theory using a new scoring method that assumes complexity in test structure, and thus takes into account every correlation between factors and items. The main specificity of the functional method is to model test scores by multiple regression instead of estimating them by using simplistic sums of points. In order to proceed, the functional method requires the creation of hyperspherical measurement space, in which item responses are expressed by their correlation with orthogonal factors. This method has three main qualities. First, measures are expressed in the absolute metric of correlations; therefore, items, scales and persons are expressed in the same measurement space using the same single metric. Second, factors are systematically orthogonal and without errors, which is optimal in order to predict other outcomes. Such predictions can be performed to estimate how one would answer to other tests, or even to model one's response strategy if it was perfectly coherent. Third, the functional method provides measures of individuals' response validity (i.e., control indices). Herein, we propose a standard procedure in order to identify whether test results are interpretable and to exclude invalid results caused by various response biases based on control indices.

Gamma Knife Surgery versus Linac for Brain Metastasis: Which Technique for Which Patient? Analysis of an Historical Cohort of 63 Consecutive Patients Harboring 130 Lesions

Objectives: We present the retrospective analysis of a single-institution experience for radiosurgery (RS) in brain metastasis (BM) with Gamma Knife (GK) and Linac. Methods: From July 2010 to July 2012, 28 patients (with 83 lesions) had RS with GK and 35 patients (with 47 lesions) with Linac. The primary outcome was the local progression-free survival (LPFS). The secondary outcome was the overall survival (OS). Apart a standard statistical analysis, we included a Cox regression model with shared frailty, to modulate the within-patient correlation (preliminary evaluation showed a significant frailty effect, meaning that the correlation within patient could be ignored). Results: The mean follow-up period was 11.7 months (median 7.9, 1.7-22.7) for GK and 18.1 (median 17, 7.5-28.7) for Linac. The median number of lesions per patient was 2.5 (1-9) in GK compared with 1 (1-3) in Linac. There were more radioresistant lesions (melanoma) and more lesions located in functional areas for the GK group. The median dose was 24 Gy (GK) compared with 20 Gy (Linac). The LPFS actuarial rate was as follows: for GK at 3, 6, 9, 12, and 17 months: 96.96, 96.96, 96.96, 88.1, and 81.5%, and remained stable till 32 months; for Linac at 3, 6, 12, 17, 24, and 33 months, it was 91.5, 91.5, 91.5, 79.9, 55.5, and 17.1%, respectively (p = 0.03, chi-square test). After the Cox regression analysis with shared frailty, the p-value was not statistically significant between groups. The median overall survival was 9.7 months for GK and 23.6 months for Linac group. Uni- and multivariate analysis showed a lower GPA score and noncontrolled systemic status were associated with lower OS. Cox regression analysis adjusting for these two parameters showed comparable OS rate. Conclusions: In this comparative report between GK and Linac, preliminary analysis showed that more difficult cases are treated by GK, with patients harboring more lesions, radioresistant tumors, and highly functional located. The groups look, in this sense, very heterogeneous at baseline. After a Cox frailty model, the LPFS rates seemed very similar (p < 0.05). The OS was similar, after adjusting for systemic status and GPA score (p < 0.05). The technical reasons for choosing GK instead of Linac were the anatomical location related to highly functional areas, histology, technical limitations of Linac movements, especially lower posterior fossa locations, or closeness of multiple lesions to highly functional areas optimal dosimetry with Linac

Functional studies of a germ-line polymorphism at codon 47 within the p53 gene.

A rare germ-line polymorphism in codon 47 of the p53 gene replaces the wild-type proline (CCG) with a serine (TCG). Restriction analysis of 101 human samples revealed the frequency of the rare allele to be 0% (n = 69) in Caucasians and 4.7% (3/64, n = 32) among African-Americans. To investigate the consequence of this amino acid substitution, a cDNA construct (p53 mut47ser) containing the mutation was introduced into a lung adenocarcinoma cell line (Calu-6) that does not express p53. A growth suppression similar to that obtained after introduction of a wild-type p53 cDNA construct was observed, in contrast to the result obtained by introduction of p53 mut143ala. Furthermore, expression of neither p53 mut47ser nor wild-type p53 was tolerated by growing cells. In transient expression assays, both mut47ser and wild-type p53 activated the expression of a reporter gene linked to a p53 binding sequence (PG13-CAT) and inhibited the expression of the luciferase gene under the control of the Rous sarcoma virus promoter (RSVluc). In the same assay, mut143ala did not activate the expression of PG13-CAT and produced only a slight inhibitory effect on RSVluc. These findings indicate that the p53 variant with a serine at codon 47 should be considered as a rare germ-line polymorphism that does not alter the growth-suppression activity of p53.

The difference of disease perception by juvenile idiopathic arthritis patients and their parents : analysis of the JAMAR questionnaire.

BACKGROUND: The JAMAR (Juvenile Arthritis Multidimensional Assessment Report) has been developed to evaluate the perception of the patient and his parents on different items: well-being, pain, functional status, quality of life, disease activity, disease course, side effects of medication, therapeutic compliance and satisfaction with illness outcome. Our aim was to compare disease's perception by JIA patients and their parents. METHODS: We included into the study 100 consecutive patients over 7 years of age. We asked both parent and child to complete the JAMAR questionnaire. For each patient we recorded demographic and disease related data. We examined the level of disagreement between children and parents for the quantitative items of the JAMAR: VAS Pain, VAS Disease Activity, VAS Well Being, Juvenile Arthritis Functional Score, HRQoL. Then we looked for a relation between discordance-rate and demographic and clinical variables. RESULTS: Children and parents' median scores for all five items were similar. Individual dyads agreement was low, with a large amount of pairs (80) discordant for at least one item. We found higher MD VAS and JADAS in more discordant dyads, suggesting that when the disease is more active discordance between child and parent increase. CONCLUSION: The JAMAR questionnaire is an important tool that helps clinicians to detect divergent child and parent's disease perceptions. It is essential that both patients and parents fill the JAMAR questionnaire for a complete clinical and psychosocial evaluation.

Genome-wide enrichment analysis between endometriosis and obesity-related traits reveals novel susceptibility loci.

Endometriosis is a chronic inflammatory condition in women that results in pelvic pain and subfertility, and has been associated with decreased body mass index (BMI). Genetic variants contributing to the heritable component have started to emerge from genome-wide association studies (GWAS), although the majority remain unknown. Unexpectedly, we observed an intergenic locus on 7p15.2 that was genome-wide significantly associated with both endometriosis and fat distribution (waist-to-hip ratio adjusted for BMI; WHRadjBMI) in an independent meta-GWAS of European ancestry individuals. This led us to investigate the potential overlap in genetic variants underlying the aetiology of endometriosis, WHRadjBMI and BMI using GWAS data. Our analyses demonstrated significant enrichment of common variants between fat distribution and endometriosis (P = 3.7 × 10(-3)), which was stronger when we restricted the investigation to more severe (Stage B) cases (P = 4.5 × 10(-4)). However, no genetic enrichment was observed between endometriosis and BMI (P = 0.79). In addition to 7p15.2, we identify four more variants with statistically significant evidence of involvement in both endometriosis and WHRadjBMI (in/near KIFAP3, CAB39L, WNT4, GRB14); two of these, KIFAP3 and CAB39L, are novel associations for both traits. KIFAP3, WNT4 and 7p15.2 are associated with the WNT signalling pathway; formal pathway analysis confirmed a statistically significant (P = 6.41 × 10(-4)) overrepresentation of shared associations in developmental processes/WNT signalling between the two traits. Our results demonstrate an example of potential biological pleiotropy that was hitherto unknown, and represent an opportunity for functional follow-up of loci and further cross-phenotype comparisons to assess how fat distribution and endometriosis pathogenesis research fields can inform each other.