Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: A Mendelian randomisation analysis

Background To investigate potential cardiovascular and other effects of long-term pharmacological interleukin 1 (IL-1) inhibition, we studied genetic variants that produce inhibition of IL-1, a master regulator of inflammation. Methods We created a genetic score combining the effects of alleles of two common variants (rs6743376 and rs1542176) that are located upstream of IL1RN, the gene encoding the IL-1 receptor antagonist (IL-1Ra; an endogenous inhibitor of both IL-1α and IL-1β); both alleles increase soluble IL-1Ra protein concentration. We compared effects on inflammation biomarkers of this genetic score with those of anakinra, the recombinant form of IL-1Ra, which has previously been studied in randomised trials of rheumatoid arthritis and other inflammatory disorders. In primary analyses, we investigated the score in relation to rheumatoid arthritis and four cardiometabolic diseases (type 2 diabetes, coronary heart disease, ischaemic stroke, and abdominal aortic aneurysm; 453 411 total participants). In exploratory analyses, we studied the relation of the score to many disease traits and to 24 other disorders of proposed relevance to IL-1 signalling (746 171 total participants). Findings For each IL1RN minor allele inherited, serum concentrations of IL-1Ra increased by 0·22 SD (95% CI 0·18–0·25; 12·5%; p=9·3 × 10−33), concentrations of interleukin 6 decreased by 0·02 SD (−0·04 to −0·01; −1·7%; p=3·5 × 10−3), and concentrations of C-reactive protein decreased by 0·03 SD (−0·04 to −0·02; −3·4%; p=7·7 × 10−14). We noted the effects of the genetic score on these inflammation biomarkers to be directionally concordant with those of anakinra. The allele count of the genetic score had roughly log-linear, dose-dependent associations with both IL-1Ra concentration and risk of coronary heart disease. For people who carried four IL-1Ra-raising alleles, the odds ratio for coronary heart disease was 1·15 (1·08–1·22; p=1·8 × 10−6) compared with people who carried no IL-1Ra-raising alleles; the per-allele odds ratio for coronary heart disease was 1·03 (1·02–1·04; p=3·9 × 10−10). Per-allele odds ratios were 0·97 (0·95–0·99; p=9·9 × 10−4) for rheumatoid arthritis, 0·99 (0·97–1·01; p=0·47) for type 2 diabetes, 1·00 (0·98–1·02; p=0·92) for ischaemic stroke, and 1·08 (1·04–1·12; p=1·8 × 10−5) for abdominal aortic aneurysm. In exploratory analyses, we observed per-allele increases in concentrations of proatherogenic lipids, including LDL-cholesterol, but no clear evidence of association for blood pressure, glycaemic traits, or any of the 24 other disorders studied. Modelling suggested that the observed increase in LDL-cholesterol could account for about a third of the association observed between the genetic score and increased coronary risk. Interpretation Human genetic data suggest that long-term dual IL-1α/β inhibition could increase cardiovascular risk and, conversely, reduce the risk of development of rheumatoid arthritis. The cardiovascular risk might, in part, be mediated through an increase in proatherogenic lipid concentrations. Funding UK Medical Research Council, British Heart Foundation, UK National Institute for Health Research, National Institute for Health Research Cambridge Biomedical Research Centre, European Research Council, and European Commission Framework Programme 7.

A new regulatory variant in the interleukin-6 receptor gene associates with asthma risk

The main genetic determinant of soluble interleukin 6 receptor (sIL-6R) levels is the missense variant rs2228145 that maps to the cleavage site of IL-6R. For each Ala allele, sIL-6R serum levels increase by ∼20 ng ml -1 and asthma risk by 1.09-fold. However, this variant does not explain the total heritability for sIL-6R levels. Additional independent variants in IL6R may therefore contribute to variation in sIL-6R levels and influence asthma risk. We imputed 471 variants in IL6R and tested these for association with sIL-6R serum levels in 360 individuals. An intronic variant (rs12083537) was associated with sIL-6R levels independently of rs4129267 (P=0.0005), a proxy single-nucleotide polymorphism for rs2228145. A significant and consistent association for rs12083537 was observed in a replication panel of 354 individuals (P=0.033). Each rs12083537:A allele increased sIL-6R serum levels by 2.4 ng ml -1. Analysis of mRNA levels in two cohorts did not identify significant associations between rs12083537 and IL6R transcription levels. On the other hand, results from 16 705 asthmatics and 30 809 controls showed that the rs12083537:A allele increased asthma risk by 1.04-fold (P=0.0419). Genetic risk scores based on IL6R regulatory variants may prove useful in explaining variation in clinical response to tocilizumab, an anti-IL-6R monoclonal antibody.

Mapping the regulatory sequences controlling 93 breast cancer-associated miRNA genes leads to the identification of two functional promoters of the Hsa-mir-200b cluster, methylation of which is associated with metastasis or hormone receptor status in advanced breast cancer

MicroRNAs (miRNAs) are small non-coding RNAs of 20 nt in length that are capable of modulating gene expression post-transcriptionally. Although miRNAs have been implicated in cancer, including breast cancer, the regulation of miRNA transcription and the role of defects in this process in cancer is not well understood. In this study we have mapped the promoters of 93 breast cancer-associated miRNAs, and then looked for associations between DNA methylation of 15 of these promoters and miRNA expression in breast cancer cells. The miRNA promoters with clearest association between DNA methylation and expression included a previously described and a novel promoter of the Hsa-mir-200b cluster. The novel promoter of the Hsa-mir-200b cluster, denoted P2, is located 2 kb upstream of the 5′ stemloop and maps within a CpG island. P2 has comparable promoter activity to the previously reported promoter (P1), and is able to drive the expression of miR-200b in its endogenous genomic context. DNA methylation of both P1 and P2 was inversely associated with miR-200b expression in eight out of nine breast cancer cell lines, and in vitro methylation of both promoters repressed their activity in reporter assays. In clinical samples, P1 and P2 were differentially methylated with methylation inversely associated with miR-200b expression. P1 was hypermethylated in metastatic lymph nodes compared with matched primary breast tumours whereas P2 hypermethylation was associated with loss of either oestrogen receptor or progesterone receptor. Hypomethylation of P2 was associated with gain of HER2 and androgen receptor expression. These data suggest an association between miR-200b regulation and breast cancer subtype and a potential use of DNA methylation of miRNA promoters as a component of a suite of breast cancer biomarkers.

Association between the interleukin 23 receptor and ankylosing spondylitis is confirmed by a new UK case-control study and meta-analysis of published series

Objectives. It has been shown previously that IL-23R variants are associated with AS. We conducted an extended analysis in the UK population and a meta-analysis with the previously published studies, in order to refine these IL-23R associations with AS. Methods. The UK case-control study included 730 new cases and 1331 healthy controls. In the extended study, the 730 cases were combined with 1088 published cases. Allelic associations were analysed using contingency tables. In the meta-analysis, 3482 cases and 3150 controls from four different published studies and the new UK cases were combined. DerSimonian-Laird test was used to calculate random effects pooled odds ratios (ORs). Results. In the UK case-control study with new cases, four of the eight SNPs showed significant associations, whereas in the extended UK study, seven of the eight IL-23R SNPs showed significant associations (P < 0.05) with AS, maximal with rs11209032 (P < 10-5, OR 1.3), when cases with IBD and/or psoriasis were excluded. The meta-analysis showed significant associations with all eight SNPs; the strongest associations were again seen not only with rs11209032 (P = 4.06 × 10-9, OR ∼1.2) but also with rs11209026 (P < 10-10, OR ∼0.6). Conclusions. IL-23R polymorphisms are clearly associated with AS, but the primary causal association(s) is(are) still not established. These polymorphisms could contribute either increased or decreased susceptibility to AS; functional studies will be required for their full evaluation. Additionally, observed stronger associations with SNPs rs11209026 and rs11465804 upon exclusion of IBD and/or psoriasis cases may represent an independent association with AS. © The Author 2009. Published by Oxford University Press on behalf of the British Society for Rheumatology. All rights reserved.

Estrogen receptor α regulates area-adjusted bone mineral content in late pubertal girls

Context: Whether the action of estrogen in skeletal development depends on estrogen receptor α as encoded by the ESR1 gene is unknown. Objectives: The aim of this study was to establish whether the gain in area-adjusted bone mineral content (ABMC) in girls occurs in late puberty and to examine whether the magnitude of this gain is related to ESR1 polymorphisms. Design: We conducted a cross-sectional analysis. Setting: The study involved the Avon Longitudinal Study of Parents and Children (ALSPAC), a population-based prospective study. Participants: Participants included 3097 11-yr-olds with DNA samples, dual x-ray absorptiometry measurements, and pubertal stage information. Outcomes: Outcome measures included separate prespecified analyses in boys and girls of the relationship between ABMC derived from total body dual x-ray absorptiometry scans and Tanner stage and of the interaction between ABMC, Tanner stage, and ESR1 polymorphisms. Results: Total body less head and spinal ABMC were higher in girls in Tanner stages 4 and 5, compared with those in Tanner stages 1, 2, and 3. In contrast, height increased throughout puberty. No differences were observed in ABMC according to Tanner stage in boys. For rs2234693 (PvuII) and rs9340799 (XbaI) polymorphisms, differences in spinal ABMC in late puberty were 2-fold greater in girls who were homozygous for the C and G alleles, respectively (P = 0.001). For rs7757956, the difference in total body less head ABMC in late puberty was 50% less in individuals homozygous or heterozygous for the A allele (P = 0.006). Conclusions: Gains in ABMC in late pubertal girls are strongly associated with ESR1 polymorphisms, suggesting that estrogen contributes to this process via an estrogen receptor α-dependent pathway.

Effect of an estrogen receptor-α intron 4 polymorphism on fat mass in 11-year-old children

Context: in the ESR1 gene encoding estrogen receptor (ER)-α may be associated with fat mass in adults. Objectives: The objective of the study was to establish whether ESR1 polymorphisms influence fat mass in childhood. Design: This was a cross-sectional analysis after genotyping of rs9340799, rs2234693, and rs7757956 ESR1 polymorphisms. Setting: The Avon Longitudinal Study of Parents and Children (ALSPAC) was a population-based prospective study. Participants: Participants included 3097 11-yr-old children with results for ESR1 genotyping, puberty measures, and dual-energy x-ray absorptiometry results. Outcomes: Relationships between ESR1 polymorphisms and indices of body composition were measured. Results: The rs7757956 polymorphism was associated with fat mass (P = 0.002). Total body fat mass (adjusted for height) was reduced by 6% in children with TA/AA genotypes, and risk of being overweight (≥85th centile of fat mass) was decreased by 20%. This genetic effect appeared to interact with puberty in girls (P = 0.05 for interaction): in those with the TT genotype, total body fat mass (adjusted for height) was 18% higher in Tanner stages 3-5 vs. stages 1-2; the equivalent difference was 7% in those with TA/AA genotypes. Furthermore, the risk of being overweight was 36% lower in girls with TA/AA genotypes in Tanner stages 3-5, but no reduction was seen in those in stages 1-2. Neither rs9340799 nor rs2234693 polymorphisms were associated with body composition measures. Conclusions: Fat mass in 11-yr-old children was related to the rs7757956 ESR1 polymorphism. This association was strongest in girls in more advanced puberty, in whom the risk of being overweight was reduced by 36% in those with the TA/AA genotype.

A recurrent mutation in the BMP type I receptor ACVR1 causes inherited and sporadic fibrodysplasia ossificans progressiva

Fibrodysplasia ossificans progressiva (FOP) is a rare autosomal dominant disorder of skeletal malformations and progressive extraskeletal ossification. We mapped FOP to chromosome 2q23-24 by linkage analysis and identified an identical heterozygous mutation (617G→A; R206H) in the glycine-serine (GS) activation domain of ACVR1, a BMP type I receptor, in all affected individuals examined. Protein modeling predicts destabilization of the GS domain, consistent with constitutive activation of ACVR1 as the underlying cause of the ectopic chondrogenesis, osteogenesis and joint fusions seen in FOP.

Genetic control of bone density and turnover: Role of the collagen 1α1, estrogen receptor, and vitamin D receptor genes

Genetic factors are known to influence both the peak bone mass and probably the rate of change in bone density. A range of regulatory and structural genes has been proposed to be involved including collagen 1α1 (COL1A1), the estrogen receptor (ER), and the vitamin D receptor (VDR), but the actual genes involved are uncertain. We therefore studied the role of the COL1A1 and VDR loci in control of bone density by linkage in 45 dizygotic twin pairs and 29 nuclear families comprising 120 individuals. The influences on bone density of polymorphisms of COL1A1, VDR, and ER were studied by association both cross-sectionally and longitudinally in 193 elderly postmenopausal women (average age, 69 years) over a mean follow-up time of 6.3 years. Weak linkage of the COL1A1 locus with bone density was observed in both twins and families (p = 0.02 in both data sets), confirming previous observations of linkage of this locus with bone density. Association between the MscI polymorphism of COL1A1 and rate of lumbar spine bone loss was observed with significant gene-environment interaction related to dietary calcium intake (p = 0.0006). In the lowest tertile of dietary calcium intake, carriers of "s" alleles lost more bone than "SS" homozygotes (p = 0.01), whereas the opposite was observed in the highest dietary calcium intake (p = 0.003). Association also was observed between rate of bone loss at both the femoral neck and the lumbar spine and the TaqI VDR polymorphism (p = 0.03). This association was strongest in those in the lowest tertile of calcium intake, also suggesting the presence of gene-environment interaction involving dietary calcium and VDR, influencing bone turnover. No significant association was observed between the PvuII ER polymorphism alone or in combination with VDR or COL1A1 genotypes, with either bone density or its rate of change. These data support the involvement of COL1A1 in determination of bone density and the interaction of both COL1A1 and VDR with calcium intake in regulation of change of bone density over time.

Suggestive linkage of the parathyroid receptor type 1 to osteoporosis

We have investigated the role of 23 candidate genes in the control of bone mineral density (BMD) by linkage studies in families of probands with osteoporosis (lumbar spine [LS] or femoral neck [FN] BMD T score < -2.5) and low BMD relative to an age- and gender-matched cohort (Z score < -2.0). One hundred and fifteen probands (35 male, 80 female) and 499 of their first- or second-degree relatives (223 males and 276 females) were recruited for the study. BMD was measured at the LS and FN using dual-energy X-ray absorptiometry and expressed as age- and gender-matched Z scores corrected for body mass index. The candidate genes studied were the androgen receptor, type I collagen A1 (COLIA1), COLIA2, COLIIA1, vitamin D receptor (VDR), colony-stimulating factor 1, calcium-sensing receptor, epidermal growth factor (EGF), estrogen receptor 1 (ESR1), fibrillin type 1, insulin-like growth factor 1, interleukin-1 alpha (IL-1α), interleukin-4 (IL-4), interleukin-6 (IL-6), interleukin-11 (IL-11), osteopontin, parathyroid hormone (PTH), PTH-related peptide, PTH receptor type 1 (PTHR1), transforming growth factor-beta 1, and tumor necrosis factors alpha and beta. Sixty-four microsatellites lying close to or within these genes were investigated for linkage with BMD. Using the program MapMaker/Sibs there was suggestive evidence of linkage between BMD and PTHR1 (maximum LOD score obtained [MLS] 2.7-3.5). Moderate evidence of linkage was also observed with EGF (MLS 1.8), COLIA1 (MLS 1.7), COLIIA1/VDR (MLS 1.7), ESR1 (MLS 1.4), IL-1α (MLS 1.4), IL-4 (MLS 1.2), and IL-6 (MLS 1.2). Variance components analysis using the program ACT, correcting for proband-wise ascertainment, also showed evidence of linkage (p ≤0.05) at markers close to or within the candidate genes IL- 1α, PTHR1, IL-6, and COLIIA1/VDR. Further studies will be required to confirm these findings, to refine the location of gene responsible for the observed linkage, and to screen the candidate genes targeted at these loci for mutations.

Topological mimicry and epitope duplication in the guanylyl cyclase C receptor

Guanylyl cyclase C (GCC) is the receptor for the gastrointestinal hormones, guanylin, and uroguanylin, in addition to the bacterial heat-stable enterotoxins, which are one of the major causes of watery diarrhea the world over. GCC is expressed in intestinal cells, colorectal tumor tissue and tumors originating from metastasis of the colorectal carcinoma. We have earlier generated a monoclonal antibody to human GCC, GCC:B10, which was useful for the immunohistochemical localization of the receptor in the rat intestine (Nandi A et al., 1997, J Cell Biochem 66:500-511), and identified its epitope to a 63-amino acid stretch in the intracellular domain of GCC. In view of the potential that this antibody has for the identification of colorectal tumors, we have characterized the epitope for GCC:B10 in this study. Overlapping peptide synthesis indicated that the epitope was contained in the sequence HIPPENIFPLE. This sequence was unique to GCC, and despite a short stretch of homology with serum amyloid protein and pertussis toxin, no cross reactivity was detected. The core epitope was delineated using a random hexameric phage display library, and two categories of sequences were identified, containing either a single, or two adjacent proline residues. No sequence identified by phage display was identical to the epitope present in GCC, indicating that phage sequences represented mimotopes of the native epitope. Alignment of these sequences with HIPPENIFPLE suggested duplication of the recognition motif, which was confirmed by peptide synthesis. These studies allowed us not only to define the requirements of epitope recognition by GCC:B10 monoclonal antibody, but also to describe a novel means of epitope recognition involving topological mimicry and probable duplication of the cognate epitope in the native guanylyl cyclase C receptor sequence.

Immunization of male bonnet monkeys (M-radiata) with a recombinant FSH receptor preparation affects testicular function and fertility

Immunization of proven fertile adult male monkeys (n = 3) with a recombinant FSH receptor protein preparation (oFSHR-P) (representing amino acids 1-134 of the extracellular domain of the receptor Mr similar to 15KDa) resulted in production of receptor blocking antibodies. The ability of the antibody to bind a particulate FSH receptor preparation and receptors in intact granulosa cells was markedly (by 30-80%) inhibited by FSH. Serum T levels and LH receptor function following immunization remained unchanged. The immunized monkeys showed a 50% reduction (p<0.001) in transformation of spermatogonia(2C) to primary spermatocytes (4C) as determined by flow cytometry and the 4C:2C ratio showed a correlative change (R 0.81, p<0.0007) with reduction in fertility index (sperm counts X motility score). Breeding studies indicated that monkeys became infertile between 242-368 days of immunization when the fertility index was in the range of 123+/-76 to 354+/-42 (compared to a value of 1602+/-384 on day 0). As the effects observed ate near identical to that seen following immunization with FSH it is suggestive that oFSHR-P can substitute for FSH in the development of a contraceptive vaccine.

Immunochemical approach to the mapping of an assembled epitope of human chorionic gonadotropin: Proximity of CTP-alpha to the receptor binding region of the beta-subunit

A single-step solid-phase RIA (SS-SPRIA) developed in our laboratory using hybridoma culture supernatants has been utilised for the quantitation of epitope-paratope interactions. Using SS-SPRIA as a quantitative tool for the assessment of epitope stability, it was found that several assembled epitopes of human chorionic gonadotropin (hCG) are differentially stable to proteolysis and chemical modification. Based on these observations an approach has been developed for identifying the amino acid residues constituting an epitopic region. This approach has now been used to map an assembled epitope at/near the receptor binding region of the hormone. The mapped site forms a part of the seat belt region and the cystine knot region (C34-C38-C88-C90-H106). The carboxy terminal region of the alpha-subunit forms a part of the epitope indicating its proximity to the receptor binding region. These results are in agreement with the reported receptor binding region identified through other approaches and the X-ray crystal structure of hCG.

Immunization of male rabbits with sheep luteal receptor to LH results in production of antibodies exhibiting hormone-agonistic and -antagonistic activities

Antibodies to LH/chorionic gonadotrophin receptor (LH/CG-R; molecular weight 67 000), isolated in a homogenous state (established by SDS-PAGE and ligand blotting) from sheep luteal membrane using human CG (hCG)-Sepharose affinity chromatography, were raised in three adult male rabbits (R-I, R-II and R-III). Each of the rabbits received 20-30 mu g oi the purified receptor in Freund's complete adjuvant at a time. Primary immunization was followed by booster injection at intervals. Production of receptor antibodies was monitored by (1) determining the dilution of the serum (IgG fraction) that could specifically bind 50% of I-125-LH/CG-R added and (2) analysing sera for any chance in testosterone levels. Following primary immunization and the first booster, all three rabbits exhibited a 2.5- to 6.0-fold increase in serum testosterone over basal levels and this effect was spread over a period of time (similar to 40 days) coinciding with the rise and fall of receptor antibodies. The maximal antibody titre (ED(50)) produced at this time ranged from 1:350 to 1:100 to below detectable limits for R-I, R-II and R-III respectively. Subsequent immunizations followed by the second booster resulted in a substantial increase in antibody titre (ED(50) of 1:5000) in R-I, but this was not accompanied by any change in serum testosterone over preimmune levels, suggesting that with the progress of immunization the character of the antibody produced had also changed. Two pools of antisera from R-I collected 10 days following the booster (at day 70 (bleed I) and day 290 (bleed II)) were used in further experiments. IgG isolated from bleed I but not from bleed II antiserum showed a dose-dependent stimulation of testosterone production by mouse Leydig cells in vitro, thus confirming the in vivo hormone-mimicking activity antibodies generated during the early immunization phase. The IgG fractions from both bleeds were, however, capable of inhibiting (1) I-125-hCG binding to crude sheep luteal membrane (EC(50) of 1:70 and 1:350 for bleed I and II antisera respectively) and (2) ovine LH-stimulated testosterone production by mouse Leydig cells in vitro, indicating the presence oi antagonistic antibodies irrespective of the period of time during which the rabbits were immunized. The: fact that bleed I-stimulated testosterone production could be inhibited in a dose-dependent manner by the addition of IgG from bleed II to the mouse Leydig cell in vitro assay system showed that the agonistic activity is intrinsic to the bleed I antibody. The receptor antibody (bleed II) was also capable of blocking LH action in vivo, as rabbits passively (for 24 h with LH/CG-R antiserum) as well as actively (for 130 days) immunized against LH/CG-R failed to respond to a bolus injection of LH (50 mu g). At no time, however, was the serum testosterone reduced below the basal level. This study clearly shows that, unlike with LH antibody, attempts to achieve an LH deficiency effect in vivo by resorting to immunization with hole LH receptor is difficult, as receptor antibodies exhibit both hormone-mimicking (agonistic) as well as hormone-blocking (antagonistic) activities.

Expression of the extracellular domain of the human heat-stable enterotoxin receptor in Escherichia coli and generation of neutralizing antibodies

The entire extracellular domain of the human heat-stable enterotoxin (ST) receptor as well as a truncated N-terminal domain were cloned as glutathione S-transferase fusion proteins and expressed in Escherichia coli. The recombinant fusion proteins were purified from both the cytosol and the inclusion body fractions by selective detergent extraction followed by glutathione-agarose affinity chromatography. The purified protein, corresponding to the entire extracellular domain, bound the stable toxin peptide with an affinity comparable to that of the native receptor characterized from the human colonic T84 cell line. No binding was observed with the N-terminal truncated fragment of the receptor under similar conditions, Polyclonal antibodies were raised to the entire extracellular domain fusion protein as well as the truncated extracellular domain fusion protein, and the antibodies were purified by affinity chromatography. Addition of the purified antibodies to T84 cells inhibited ST binding and abolished ST-mediated cGMP production, indicating that critical epitopes involved in ligand interaction are present in the N-terminal fragment of the receptor, Purified antibodies recognized a single protein of M(r) 160,000 Da on Western blotting with T84 membranes, corresponding to a size of the native glycosylated receptor in T84 cells. These studies are the first report of the expression, purification, and characterization of any member of the guanylyl cyclase family of receptors in E. coli and show that binding of the toxin to the extracellular domain of the receptor is possible in the absence of any posttranslational modifications such as glycosylation. The recombinant fusion proteins as well as the antibodies that we have generated could serve as useful tools in the identification of critical residues of the extracellular domain involved in ligand interaction.