Interrogating causal pathways linking genetic variants, small molecule metabolites, and circulating lipids

BACKGROUND

Emerging technologies based on mass spectrometry or nuclear magnetic resonance enable the monitoring of hundreds of small metabolites from tissues or body fluids. Profiling of metabolites can help elucidate causal pathways linking established genetic variants to known disease risk factors such as blood lipid traits.

METHODS

We applied statistical methodology to dissect causal relationships between single nucleotide polymorphisms, metabolite concentrations, and serum lipid traits, focusing on 95 genetic loci reproducibly associated with the four main serum lipids (total-, low-density lipoprotein-, and high-density lipoprotein- cholesterol and triglycerides). The dataset used included 2,973 individuals from two independent population-based cohorts with data for 151 small molecule metabolites and four main serum lipids. Three statistical approaches, namely conditional analysis, Mendelian randomization, and structural equation modeling, were compared to investigate causal relationship at sets of a single nucleotide polymorphism, a metabolite, and a lipid trait associated with one another.

RESULTS

A subset of three lipid-associated loci (FADS1, GCKR, and LPA) have a statistically significant association with at least one main lipid and one metabolite concentration in our data, defining a total of 38 cross-associated sets of a single nucleotide polymorphism, a metabolite and a lipid trait. Structural equation modeling provided sufficient discrimination to indicate that the association of a single nucleotide polymorphism with a lipid trait was mediated through a metabolite at 15 of the 38 sets, and involving variants at the FADS1 and GCKR loci.

CONCLUSIONS

These data provide a framework for evaluating the causal role of components of the metabolome (or other intermediate factors) in mediating the association between established genetic variants and diseases or traits.

Expression of phosphofructokinase in skeletal muscle is influenced by genetic variation and associated with insulin sensitivity

Using an integrative approach in which genetic variation, gene expression, and clinical phenotypes are assessed in relevant tissues may help functionally characterize the contribution of genetics to disease susceptibility. We sought to identify genetic variation influencing skeletal muscle gene expression (expression quantitative trait loci [eQTLs]) as well as expression associated with measures of insulin sensitivity. We investigated associations of 3,799,401 genetic variants in expression of >7,000 genes from three cohorts (n = 104). We identified 287 genes with cis-acting eQTLs (false discovery rate [FDR] <5%; P < 1.96 × 10(-5)) and 49 expression-insulin sensitivity phenotype associations (i.e., fasting insulin, homeostasis model assessment-insulin resistance, and BMI) (FDR <5%; P = 1.34 × 10(-4)). One of these associations, fasting insulin/phosphofructokinase (PFKM), overlaps with an eQTL. Furthermore, the expression of PFKM, a rate-limiting enzyme in glycolysis, was nominally associated with glucose uptake in skeletal muscle (P = 0.026; n = 42) and overexpressed (Bonferroni-corrected P = 0.03) in skeletal muscle of patients with T2D (n = 102) compared with normoglycemic controls (n = 87). The PFKM eQTL (rs4547172; P = 7.69 × 10(-6)) was nominally associated with glucose uptake, glucose oxidation rate, intramuscular triglyceride content, and metabolic flexibility (P = 0.016-0.048; n = 178). We explored eQTL results using published data from genome-wide association studies (DIAGRAM and MAGIC), and a proxy for the PFKM eQTL (rs11168327; r(2) = 0.75) was nominally associated with T2D (DIAGRAM P = 2.7 × 10(-3)). Taken together, our analysis highlights PFKM as a potential regulator of skeletal muscle insulin sensitivity.

Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility

To further understanding of the genetic basis of type 2 diabetes (T2D) susceptibility, we aggregated published meta-analyses of genome-wide association studies (GWAS) including 26,488 cases and 83,964 controls of European, East Asian, South Asian, and Mexican and Mexican American ancestry. We observed significant excess in directional consistency of T2D risk alleles across ancestry groups, even at SNPs demonstrating only weak evidence of association. By following up the strongest signals of association from the trans-ethnic meta-analysis in an additional 21,491 cases and 55,647 controls of European ancestry, we identified seven novel T2D susceptibility loci. Furthermore, we observed considerable improvements in fine-mapping resolution of common variant association signals at several T2D susceptibility loci. These observations highlight the benefits of trans-ethnic GWAS for the discovery and characterisation of complex trait loci, and emphasize an exciting opportunity to extend insight into the genetic architecture and pathogenesis of human diseases across populations of diverse ancestry.

Global analysis of DNA methylation variation in adipose tissue from twins reveals links to disease-associated variants in distal regulatory elements

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h(2)median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.

Gene expression changes with age in skin, adipose tissue, blood and brain

BACKGROUND

Previous studies have demonstrated that gene expression levels change with age. These changes are hypothesized to influence the aging rate of an individual. We analyzed gene expression changes with age in abdominal skin, subcutaneous adipose tissue and lymphoblastoid cell lines in 856 female twins in the age range of 39-85 years. Additionally, we investigated genotypic variants involved in genotype-by-age interactions to understand how the genomic regulation of gene expression alters with age.

RESULTS

Using a linear mixed model, differential expression with age was identified in 1,672 genes in skin and 188 genes in adipose tissue. Only two genes expressed in lymphoblastoid cell lines showed significant changes with age. Genes significantly regulated by age were compared with expression profiles in 10 brain regions from 100 postmortem brains aged 16 to 83 years. We identified only one age-related gene common to the three tissues. There were 12 genes that showed differential expression with age in both skin and brain tissue and three common to adipose and brain tissues.

CONCLUSIONS

Skin showed the most age-related gene expression changes of all the tissues investigated, with many of the genes being previously implicated in fatty acid metabolism, mitochondrial activity, cancer and splicing. A significant proportion of age-related changes in gene expression appear to be tissue-specific with only a few genes sharing an age effect in expression across tissues. More research is needed to improve our understanding of the genetic influences on aging and the relationship with age-related diseases.

The presence of methylation quantitative trait loci indicates a direct genetic influence on the level of DNA methylation in adipose tissue

Genetic variants that associate with DNA methylation at CpG sites (methylation quantitative trait loci, meQTLs) offer a potential biological mechanism of action for disease associated SNPs. We investigated whether meQTLs exist in abdominal subcutaneous adipose tissue (SAT) and if CpG methylation associates with metabolic syndrome (MetSyn) phenotypes. We profiled 27,718 genomic regions in abdominal SAT samples of 38 unrelated individuals using differential methylation hybridization (DMH) together with genotypes at 5,227,243 SNPs and expression of 17,209 mRNA transcripts. Validation and replication of significant meQTLs was pursued in an independent cohort of 181 female twins. We find that, at 5% false discovery rate, methylation levels of 149 DMH regions associate with at least one SNP in a ±500 kilobase cis-region in our primary study. We sought to validate 19 of these in the replication study and find that five of these significantly associate with the corresponding meQTL SNPs from the primary study. We find that none of the 149 meQTL top SNPs is a significant expression quantitative trait locus in our expression data, but we observed association between expression levels of two mRNA transcripts and cis-methylation status. Our results indicate that DNA CpG methylation in abdominal SAT is partly under genetic control. This study provides a starting point for future investigations of DNA methylation in adipose tissue.

Large-scale association analysis identifies new risk loci for coronary artery disease

Coronary artery disease (CAD) is the commonest cause of death. Here, we report an association analysis in 63,746 CAD cases and 130,681 controls identifying 15 loci reaching genome-wide significance, taking the number of susceptibility loci for CAD to 46, and a further 104 independent variants (r(2) < 0.2) strongly associated with CAD at a 5% false discovery rate (FDR). Together, these variants explain approximately 10.6% of CAD heritability. Of the 46 genome-wide significant lead SNPs, 12 show a significant association with a lipid trait, and 5 show a significant association with blood pressure, but none is significantly associated with diabetes. Network analysis with 233 candidate genes (loci at 10% FDR) generated 5 interaction networks comprising 85% of these putative genes involved in CAD. The four most significant pathways mapping to these networks are linked to lipid metabolism and inflammation, underscoring the causal role of these activities in the genetic etiology of CAD. Our study provides insights into the genetic basis of CAD and identifies key biological pathways.

A comparison of the whole genome approach of MeDIP-seq to the targeted approach of the Infinium HumanMethylation450 BeadChip(®) for methylome profiling

DNA methylation is one of the most studied epigenetic marks in the human genome, with the result that the desire to map the human methylome has driven the development of several methods to map DNA methylation on a genomic scale. Our study presents the first comparison of two of these techniques - the targeted approach of the Infinium HumanMethylation450 BeadChip® with the immunoprecipitation and sequencing-based method, MeDIP-seq. Both methods were initially validated with respect to bisulfite sequencing as the gold standard and then assessed in terms of coverage, resolution and accuracy. The regions of the methylome that can be assayed by both methods and those that can only be assayed by one method were determined and the discovery of differentially methylated regions (DMRs) by both techniques was examined. Our results show that the Infinium HumanMethylation450 BeadChip® and MeDIP-seq show a good positive correlation (Spearman correlation of 0.68) on a genome-wide scale and can both be used successfully to determine differentially methylated loci in RefSeq genes, CpG islands, shores and shelves. MeDIP-seq however, allows a wider interrogation of methylated regions of the human genome, including thousands of non-RefSeq genes and repetitive elements, all of which may be of importance in disease. In our study MeDIP-seq allowed the detection of 15,709 differentially methylated regions, nearly twice as many as the array-based method (8070), which may result in a more comprehensive study of the methylome.

Assessment of gene-by-sex interaction effect on bone mineral density

Sexual dimorphism in various bone phenotypes, including bone mineral density (BMD), is widely observed; however, the extent to which genes explain these sex differences is unclear. To identify variants with different effects by sex, we examined gene-by-sex autosomal interactions genome-wide, and performed expression quantitative trait loci (eQTL) analysis and bioinformatics network analysis. We conducted an autosomal genome-wide meta-analysis of gene-by-sex interaction on lumbar spine (LS) and femoral neck (FN) BMD in 25,353 individuals from 8 cohorts. In a second stage, we followed up the 12 top single-nucleotide polymorphisms (SNPs; p < 1 × 10(-5) ) in an additional set of 24,763 individuals. Gene-by-sex interaction and sex-specific effects were examined in these 12 SNPs. We detected one novel genome-wide significant interaction associated with LS-BMD at the Chr3p26.1-p25.1 locus, near the GRM7 gene (male effect = 0.02 and p = 3.0 × 10(-5) ; female effect = -0.007 and p = 3.3 × 10(-2) ), and 11 suggestive loci associated with either FN- or LS-BMD in discovery cohorts. However, there was no evidence for genome-wide significant (p < 5 × 10(-8) ) gene-by-sex interaction in the joint analysis of discovery and replication cohorts. Despite the large collaborative effort, no genome-wide significant evidence for gene-by-sex interaction was found to influence BMD variation in this screen of autosomal markers. If they exist, gene-by-sex interactions for BMD probably have weak effects, accounting for less than 0.08% of the variation in these traits per implicated SNP. © 2012 American Society for Bone and Mineral Research.

Mapping cis- and trans-regulatory effects across multiple tissues in twins

Sequence-based variation in gene expression is a key driver of disease risk. Common variants regulating expression in cis have been mapped in many expression quantitative trait locus (eQTL) studies, typically in single tissues from unrelated individuals. Here, we present a comprehensive analysis of gene expression across multiple tissues conducted in a large set of mono- and dizygotic twins that allows systematic dissection of genetic (cis and trans) and non-genetic effects on gene expression. Using identity-by-descent estimates, we show that at least 40% of the total heritable cis effect on expression cannot be accounted for by common cis variants, a finding that reveals the contribution of low-frequency and rare regulatory variants with respect to both transcriptional regulation and complex trait susceptibility. We show that a substantial proportion of gene expression heritability is trans to the structural gene, and we identify several replicating trans variants that act predominantly in a tissue-restricted manner and may regulate the transcription of many genes.

Impact of common variation in bone-related genes on type 2 diabetes and related traits

Exploring genetic pleiotropy can provide clues to a mechanism underlying the observed epidemiological association between type 2 diabetes and heightened fracture risk. We examined genetic variants associated with bone mineral density (BMD) for association with type 2 diabetes and glycemic traits in large well-phenotyped and -genotyped consortia. We undertook follow-up analysis in ∼19,000 individuals and assessed gene expression. We queried single nucleotide polymorphisms (SNPs) associated with BMD at levels of genome-wide significance, variants in linkage disequilibrium (r(2) > 0.5), and BMD candidate genes. SNP rs6867040, at the ITGA1 locus, was associated with a 0.0166 mmol/L (0.004) increase in fasting glucose per C allele in the combined analysis. Genetic variants in the ITGA1 locus were associated with its expression in the liver but not in adipose tissue. ITGA1 variants appeared among the top loci associated with type 2 diabetes, fasting insulin, β-cell function by homeostasis model assessment, and 2-h post-oral glucose tolerance test glucose and insulin levels. ITGA1 has demonstrated genetic pleiotropy in prior studies, and its suggested role in liver fibrosis, insulin secretion, and bone healing lends credence to its contribution to both osteoporosis and type 2 diabetes. These findings further underscore the link between skeletal and glucose metabolism and highlight a locus to direct future investigations.

Epigenome-wide scans identify differentially methylated regions for age and age-related phenotypes in a healthy ageing population

Age-related changes in DNA methylation have been implicated in cellular senescence and longevity, yet the causes and functional consequences of these variants remain unclear. To elucidate the role of age-related epigenetic changes in healthy ageing and potential longevity, we tested for association between whole-blood DNA methylation patterns in 172 female twins aged 32 to 80 with age and age-related phenotypes. Twin-based DNA methylation levels at 26,690 CpG-sites showed evidence for mean genome-wide heritability of 18%, which was supported by the identification of 1,537 CpG-sites with methylation QTLs in cis at FDR 5%. We performed genome-wide analyses to discover differentially methylated regions (DMRs) for sixteen age-related phenotypes (ap-DMRs) and chronological age (a-DMRs). Epigenome-wide association scans (EWAS) identified age-related phenotype DMRs (ap-DMRs) associated with LDL (STAT5A), lung function (WT1), and maternal longevity (ARL4A, TBX20). In contrast, EWAS for chronological age identified hundreds of predominantly hyper-methylated age DMRs (490 a-DMRs at FDR 5%), of which only one (TBX20) was also associated with an age-related phenotype. Therefore, the majority of age-related changes in DNA methylation are not associated with phenotypic measures of healthy ageing in later life. We replicated a large proportion of a-DMRs in a sample of 44 younger adult MZ twins aged 20 to 61, suggesting that a-DMRs may initiate at an earlier age. We next explored potential genetic and environmental mechanisms underlying a-DMRs and ap-DMRs. Genome-wide overlap across cis-meQTLs, genotype-phenotype associations, and EWAS ap-DMRs identified CpG-sites that had cis-meQTLs with evidence for genotype–phenotype association, where the CpG-site was also an ap-DMR for the same phenotype. Monozygotic twin methylation difference analyses identified one potential environmentally-mediated ap-DMR associated with total cholesterol and LDL (CSMD1). Our results suggest that in a small set of genes DNA methylation may be a candidate mechanism of mediating not only environmental, but also genetic effects on age-related phenotypes.

Epigenetic patterns vary during healthy ageing and development. Age-related DNA methylation changes have been implicated in cellular senescence and longevity, yet the causes and functional consequences of these variants remain unclear. To understand the biological mechanisms involved in potential longevity and rate of healthy ageing, we performed genome-wide association of epigenetic and genetic variation with both chronological age and age-related phenotypes. We identified hundreds of DNA methylation variants significantly associated with age and replicated these in an independent sample of young adult twins. Only a small proportion of these variants were also associated with age-related phenotypes. Therefore, the majority of age-related epigenetic changes do not contribute to rate of healthy ageing at later stages in life. Our results suggest that age-related changes in methylation occur throughout an individual's lifespan and that a proportion of these may be initiated from an early age. Intriguingly, a fraction of the age differentially methylated regions also associated with genetic variants in our sample, suggesting that DNA methylation may be a candidate mechanism of mediating not only environmental but also genetic effects on age-related phenotypes.

Human metabolic individuality in biomedical and pharmaceutical research

Genome-wide association studies (GWAS) have identified many risk loci for complex diseases, but effect sizes are typically small and information on the underlying biological processes is often lacking. Associations with metabolic traits as functional intermediates can overcome these problems and potentially inform individualized therapy. Here we report a comprehensive analysis of genotype-dependent metabolic phenotypes using a GWAS with non-targeted metabolomics. We identified 37 genetic loci associated with blood metabolite concentrations, of which 25 show effect sizes that are unusually high for GWAS and account for 10-60% differences in metabolite levels per allele copy. Our associations provide new functional insights for many disease-related associations that have been reported in previous studies, including those for cardiovascular and kidney disorders, type 2 diabetes, cancer, gout, venous thromboembolism and Crohn's disease. The study advances our knowledge of the genetic basis of metabolic individuality in humans and generates many new hypotheses for biomedical and pharmaceutical research.

The architecture of gene regulatory variation across multiple human tissues: the MuTHER study

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.

Regulation of gene expression is a fundamental cellular process determining a large proportion of the phenotypic variance. Previous studies have identified genetic loci influencing gene expression levels (eQTLs), but the complexity of their tissue-specific properties has not yet been well-characterized. In this study, we perform cis-eQTL analysis in a unique matched co-twin design for three human tissues derived simultaneously from the same set of individuals. The study design allows validation of the substantial discoveries we make in each tissue. We explore in depth the tissue-dependent features of regulatory variants and estimate the proportions of shared and specific effects. We use continuous measures of eQTL sharing to circumvent the statistical power limitations of comparing direct overlap of eQTLs in multiple tissues. In this framework, we demonstrate that 30% of eQTLs are shared among tissues, while 29% are exclusively tissue-specific. Furthermore, we show that the fold change in expression between eQTL genotypic classes differs between tissues. Even among shared eQTLs, we report a substantial proportion (10%–20%) of significant tissue differences in magnitude of these effects. The complexities we highlight here are essential for understanding the impact of regulatory variants on complex traits.

Global analysis of the impact of environmental perturbation on cis-regulation of gene expression

Genetic variants altering cis-regulation of normal gene expression (cis-eQTLs) have been extensively mapped in human cells and tissues, but the extent by which controlled, environmental perturbation influences cis-eQTLs is unclear. We carried out large-scale induction experiments using primary human bone cells derived from unrelated donors of Swedish origin treated with 18 different stimuli (7 treatments and 2 controls, each assessed at 2 time points). The treatments with the largest impact on the transcriptome, verified on two independent expression arrays, included BMP-2 (t = 2h), dexamethasone (DEX) (t = 24h), and PGE₂ (t = 24h). Using these treatments and control, we performed expression profiling for 18,144 RefSeq transcripts on biological replicates of the complete study cohort of 113 individuals (n_total = 782) and combined it with genome-wide SNP-genotyping data in order to map treatment-specific cis-eQTLs (defined as SNPs located within the gene ±250 kb). We found that 93% of cis-eQTLs at 1% FDR were observed in at least one additional treatment, and in fact, on average, only 1.4% of the cis-eQTLs were considered as treatment-specific at high confidence. The relative invariability of cis-regulation following perturbation was reiterated independently by genome-wide allelic expression tests where only a small proportion of variance could be attributed to treatment. Treatment-specific cis-regulatory effects were, however, 2- to 6-fold more abundant among differently expressed genes upon treatment. We further followed-up and validated the DEX–specific cis-regulation of the MYO6 and TNC loci and found top cis-regulatory variants located 180 kb and 250 kb upstream of the transcription start sites, respectively. Our results suggest that, as opposed to tissue-specificity of cis-eQTLs, the interactions between cellular environment and cis-variants are relatively rare (∼1.5%), but that detection of such specific interactions can be achieved by a combination of functional genomic approaches as described here.

Population variation in normal gene expression has been convincingly shown to be under strong genetic control where the main genetic variants are located within close proximity to the gene itself (so called cis-acting). However, the extent to which controlled, environmental stimuli influences cis-regulation of gene expression is unclear. Here, we combine different functional genomic approaches and examine the role of common genetic variants on induced gene expression in a population panel of primary human cells derived from ∼100 unrelated donors treated under multiple conditions. Using these approaches, we find that the interaction between cellular environment and cis-variants are relatively rare, with only a small proportion of the identified genetic variants being specific to treatment. However, although treatment-specific genetic regulation of gene expression seems to be infrequent, we prove its existence by thorough validation of treatment-specific effects of the glucocorticoid-specific regulation of TNC expression. Taken together, these findings indicate that the regulatory landscape within a cell is very stable but, by combining functional genomic tools gene-environmental interactions of clinical importance, can be detected and possibly used as biomarkers in future pharmacogenomic studies.

Analysis of the impact of genetic variation on human gene expression

Interindividual variation in gene expression has been convincingly shown to be controlled, in part, by genetic differences. Determining the architecture of genetic variation, the underlying gene expression may allow deeper insight into complex phenotypes, such as differences in disease susceptibility. Mapping genetic variants accounting for expression phenotypes in human cell and tissue panels has rapidly progressed from proof-of-principle experiments to general tools in biomedical discovery. We discuss the general approach and critical considerations for carrying out expression quantitative trait mapping in human tissues.

A genome-wide association study reveals variants in ARL15 that influence adiponectin levels

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P< or =5x10(-8)). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P< or =0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2x10(-19) for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9x10(-8), n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5x10(-6), n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2x10(-3), n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.

Global patterns of cis variation in human cells revealed by high-density allelic expression analysis

Cis-acting variants altering gene expression are a source of phenotypic differences. The cis-acting components of expression variation can be identified through the mapping of differences in allelic expression (AE), which is the measure of relative expression between two allelic transcripts. We generated a map of AE associated SNPs using quantitative measurements of AE on Illumina Human1M BeadChips. In 53 lymphoblastoid cell lines derived from donors of European descent, we identified common cis variants affecting 30% (2935/9751) of the measured RefSeq transcripts at 0.001 permutation significance. The pervasive influence of cis-regulatory variants, which explain 50% of population variation in AE, extend to full-length transcripts and their isoforms as well as to unannotated transcripts. These strong effects facilitate fine mapping of cis-regulatory SNPs, as demonstrated by dissection of heritable control of transcripts in the systemic lupus erythematosus-associated C8orf13-BLK region in chromosome 8. The dense collection of associations will facilitate large-scale isolation of cis-regulatory SNPs.

Twenty bone-mineral-density loci identified by large-scale meta-analysis of genome-wide association studies

Bone mineral density (BMD) is a heritable complex trait used in the clinical diagnosis of osteoporosis and the assessment of fracture risk. We performed meta-analysis of five genome-wide association studies of femoral neck and lumbar spine BMD in 19,195 subjects of Northern European descent. We identified 20 BMD loci that reached genome-wide significance (GWS; P < 5 x 10(-8)), of which 13 map to regions not previously associated with this trait: 1p31.3 (GPR177), 2p21 (SPTBN1), 3p22 (CTNNB1), 4q21.1 (MEPE), 5q14 (MEF2C), 7p14 (STARD3NL), 7q21.3 (FLJ42280), 11p11.2 (LRP4, ARHGAP1, F2), 11p14.1 (DCDC5), 11p15 (SOX6), 16q24 (FOXL1), 17q21 (HDAC5) and 17q12 (CRHR1). The meta-analysis also confirmed at GWS level seven known BMD loci on 1p36 (ZBTB40), 6q25 (ESR1), 8q24 (TNFRSF11B), 11q13.4 (LRP5), 12q13 (SP7), 13q14 (TNFSF11) and 18q21 (TNFRSF11A). The many SNPs associated with BMD map to genes in signaling pathways with relevance to bone metabolism and highlight the complex genetic architecture that underlies osteoporosis and variation in BMD.

Spectrum of mutations in MMACHC, allelic expression, and evidence for genotype-phenotype correlations

Methylmalonic aciduria and homocystinuria, cblC type, is a rare disorder of intracellular vitamin B(12) (cobalamin [Cbl]) metabolism caused by mutations in the MMACHC gene. MMACHC was sequenced from the gDNA of 118 cblC individuals. Eleven novel mutations were identified, as well as 23 mutations that were observed previously. Six sequence variants capture haplotype diversity in individuals across the MMACHC interval. Genotype-phenotype correlations of common mutations were apparent; individuals with c.394C>T tend to present with late-onset disease whereas patients with c.331C>T and c.271dupA tend to present in infancy. Other missense variants were also associated with late- or early-onset disease. Allelic expression analysis was carried out on human cblC fibroblasts compound heterozygous for different combinations of mutations including c.271dupA, c.331C>T, c.394C>T, and c.482G>A. The early-onset c.271dupA mutation was consistently underexpressed when compared to control alleles and the late-onset c.394C>T and c.482G>A mutations. The early-onset c.331C>T mutation was also underexpressed when compared to control alleles and the c.394C>T mutation. Levels of MMACHC mRNA transcript in cell lines homozygous for c.271dupA, c.331C>T, and c.394C>T were assessed using quantitative real-time RT-PCR. Cell lines homozygous for the late onset c.394C>T mutation had significantly higher levels of transcript when compared to cell lines homozygous for the early-onset mutations. Differential or preferential MMACHC transcript levels may provide a clue as to why individuals carrying c.394C>T generally present later in life.

Allele-specific chromatin remodeling in the ZPBP2/GSDMB/ORMDL3 locus associated with the risk of asthma and autoimmune disease

Common SNPs in the chromosome 17q12-q21 region alter the risk for asthma, type 1 diabetes, primary biliary cirrhosis, and Crohn disease. Previous reports by us and others have linked the disease-associated genetic variants with changes in expression of GSDMB and ORMDL3 transcripts in human lymphoblastoid cell lines (LCLs). The variants also alter regulation of other transcripts, and this domain-wide cis-regulatory effect suggests a mechanism involving long-range chromatin interactions. Here, we further dissect the disease-linked haplotype and identify putative causal DNA variants via a combination of genetic and functional analyses. First, high-throughput resequencing of the region and genotyping of potential candidate variants were performed. Next, additional mapping of allelic expression differences in Yoruba HapMap LCLs allowed us to fine-map the basis of the cis-regulatory differences to a handful of candidate functional variants. Functional assays identified allele-specific differences in nucleosome distribution, an allele-specific association with the insulator protein CTCF, as well as a weak promoter activity for rs12936231. Overall, this study shows a common disease allele linked to changes in CTCF binding and nucleosome occupancy leading to altered domain-wide cis-regulation. Finally, a strong association between asthma and cis-regulatory haplotypes was observed in three independent family-based cohorts (p = 1.78 x 10(-8)). This study demonstrates the requirement of multiple parallel allele-specific tools for the investigation of noncoding disease variants and functional fine-mapping of human disease-associated haplotypes.

Tissue effect on genetic control of transcript isoform variation

Current genome-wide association studies (GWAS) are moving towards the use of large cohorts of primary cell lines to study a disease of interest and to assign biological relevance to the genetic signals identified. Here, we use a panel of human osteoblasts (HObs) to carry out a transcriptomic survey, similar to recent studies in lymphoblastoid cell lines (LCLs). The distinct nature of HObs and LCLs is reflected by the preferential grouping of cell type–specific genes within biologically and functionally relevant pathways unique to each tissue type. We performed cis-association analysis with SNP genotypes to identify genetic variations of transcript isoforms, and our analysis indicates that differential expression of transcript isoforms in HObs is also partly controlled by cis-regulatory genetic variants. These isoforms are regulated by genetic variants in both a tissue-specific and tissue-independent fashion, and these associations have been confirmed by RT–PCR validation. Our study suggests that multiple transcript isoforms are often present in both tissues and that genetic control may affect the relative expression of one isoform to another, rather than having an all-or-none effect. Examination of the top SNPs from a GWAS of bone mineral density show overlap with probeset associations observed in this study. The top hit corresponding to the FAM118A gene was tested for association studies in two additional clinical studies, revealing a novel transcript isoform variant. Our approach to examining transcriptome variation in multiple tissue types is useful for detecting the proportion of genetic variation common to different cell types and for the identification of cell-specific isoform variants that may be functionally relevant, an important follow-up step for GWAS.

The transcriptome of any given cell type is a complex program of controlled gene expression underlying its biological function. An additional layer of molecular complexity involving individual genetic variation can modulate the transcriptome within the same tissue type, conferring potential phenotypic differences between individuals at the cellular level. This study highlights common and unique aspects of the transcriptome between the well-characterized lymphoblastoid cell lines from the International HapMap Project and those of a cultured primary cell type, human osteoblasts. We observe that inter-individual genetic variation can regulate transcript isoform expression in tissue-specific and tissue-independent manners, indicating that genetic differences among individuals can alter the transcriptome in one or more tissues, ultimately leading to altered biological functions within the lymphoblasts and/or osteoblasts. Pursuant to this, genome wide association studies on bone mineral density (BMD) have identified a number of significant loci and polymorphisms highly linked to the BMD quantitative phenotype. A small proportion of these polymorphisms overlap with our highly significant SNPs regulating the osteoblast transcriptome, revealing a potential molecular basis for this phenotype at the transcriptional level. This study highlights the importance of examining the differing transcriptomes and cis-regulatory mechanisms governing the biological and functional roles of varied tissue types.

Population genomics in a disease targeted primary cell model

The common genetic variants associated with complex traits typically lie in noncoding DNA and may alter gene regulation in a cell type-specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an expression quantitative trait loci (eQTL) study of primary human osteoblasts (HOb) derived from 95 unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top 2000 BMD-associated SNPs (P < approximately 10(-3)) were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment (threefold) of converging cis-eQTLs as compared to LCLs. The top 10 BMD loci with SNPs showing strong cis-effects on gene expression in HObs (P = 6 x 10(-10) - 7 x 10(-16)) were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All 10 variants were tested in the Swedish MrOS Cohort (n = 3014), providing evidence for two novel BMD loci (SRR and MSH3). These variants were then tested in the Rotterdam Study (n = 2090), yielding converging evidence for BMD association at the 17p13.3 SRR locus (P(combined) = 5.6 x 10(-5)). The cis-regulatory effect was further fine-mapped to the proximal promoter of the SRR gene (rs3744270, r(2) = 0.5, P = 2.6 x 10(-15)). Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.

Type I collagen alpha1 Sp1 polymorphism and the risk of cruciate ligament ruptures or shoulder dislocations

BACKGROUND

Cruciate ligament ruptures and shoulder dislocations are often caused by trauma, but predisposing intrinsic factors might also influence the risk. These injuries are more common in those with a previously injured sibling, an observation that might indicate a genetic predisposition. It is well known that polymorphisms in the collagen I gene are associated not only with osteoporosis and osteoporotic fracture risk, but also with osteoarthritis.

HYPOTHESIS

Because collagen I is abundant in ligaments and tendons, the authors hypothesized that collagen I alpha1 Sp1 polymorphism also was related to the occurrence of cruciate ligament ruptures and shoulder dislocations.

STUDY DESIGN

Case-control study; Level of evidence, 3.

METHODS

A total of 358 patients and 325 randomly selected population-based female controls were included in the study. Of the cases, 233 had a cruciate ligament rupture and 126 had had a shoulder dislocation. Age-adjusted odds ratios (ORs) with 95% confidence intervals (CIs) estimated by unconditional logistic regression were used as measures of association.

RESULTS

Compared with the homozygous SS category, the heterozygous participants displayed a similar risk (OR, 1.06; 95% CI, 0.76-1.49), whereas the ss genotype was underrepresented in the injured population compared with the controls (OR, 0.15; 95% CI, 0.03-0.68). This latter estimate was similar for both cruciate ligament ruptures and shoulder dislocations, and was furthermore not modified by general joint laxity.

CONCLUSION

The authors found a substantially decreased risk of these injuries associated with collagen type I alpha1 Sp1 polymorphism. The study might encourage other investigators to consider further research in the area of genes and soft tissue injuries.

Targeted screening of cis-regulatory variation in human haplotypes

Regulatory cis-acting variants account for a large proportion of gene expression variability in populations. Cis-acting differences can be specifically measured by comparing relative levels of allelic transcripts within a sample. Allelic expression (AE) mapping for cis-regulatory variant discovery has been hindered by the requirements of having informative or heterozygous single nucleotide polymorphisms (SNPs) within genes in order to assign the allelic origin of each transcript. In this study we have developed an approach to systematically screen for heritable cis-variants in common human haplotypes across >1,000 genes. In order to achieve the highest level of information per haplotype studied, we carried out allelic expression measurements by using both intronic and exonic SNPs in primary transcripts. We used a novel RNA pooling strategy in immortalized lymphoblastoid cell lines (LCLs) and primary human osteoblast cell lines (HObs) to allow for high-throughput AE. Screening hits from RNA pools were further validated by performing allelic expression mapping in individual samples. Our results indicate that >10% of expressed genes in human LCLs show genotype-linked AE. In addition, we have validated cis-acting variants in over 20 genes linked with common disease susceptibility in recent genome-wide studies. More generally, our results indicate that RNA pooling coupled with AE read-out by second generation sequencing or by other methods provides a high-throughput tool for cataloging the impact of common noncoding variants in the human genome.

A risk haplotype of STAT4 for systemic lupus erythematosus is over-expressed, correlates with anti-dsDNA and shows additive effects with two risk alleles of IRF5

Systemic lupus erythematosus (SLE) is the prototype autoimmune disease where genes regulated by type I interferon (IFN) are over-expressed and contribute to the disease pathogenesis. Because signal transducer and activator of transcription 4 (STAT4) plays a key role in the type I IFN receptor signaling, we performed a candidate gene study of a comprehensive set of single nucleotide polymorphism (SNPs) in STAT4 in Swedish patients with SLE. We found that 10 out of 53 analyzed SNPs in STAT4 were associated with SLE, with the strongest signal of association (P = 7.1 x 10(-8)) for two perfectly linked SNPs rs10181656 and rs7582694. The risk alleles of these 10 SNPs form a common risk haplotype for SLE (P = 1.7 x 10(-5)). According to conditional logistic regression analysis the SNP rs10181656 or rs7582694 accounts for all of the observed association signal. By quantitative analysis of the allelic expression of STAT4 we found that the risk allele of STAT4 was over-expressed in primary human cells of mesenchymal origin, but not in B-cells, and that the risk allele of STAT4 was over-expressed (P = 8.4 x 10(-5)) in cells carrying the risk haplotype for SLE compared with cells with a non-risk haplotype. The risk allele of the SNP rs7582694 in STAT4 correlated to production of anti-dsDNA (double-stranded DNA) antibodies and displayed a multiplicatively increased, 1.82-fold risk of SLE with two independent risk alleles of the IRF5 (interferon regulatory factor 5) gene.

Large-scale association study between two coding LRP5 gene polymorphisms and bone phenotypes and fractures in men

Herein we investigated the association between polymorphisms in the LRP5 gene and bone phenotypes and fractures in three large male cohorts based on the rationale that mutations in LRP5 cause severe bone phenotypes. Results showed an association of the Val667Met SNP with spine BMD in 3,800 young and elderly men.

INTRODUCTION

The low-density lipoprotein receptor-related protein 5 (LRP5)-Wnt signalling system is of importance for regulating osteoblastic activity, which became clear after findings that inactivating mutations in LRP5 cause osteoporosis. The overall aim of this study was to investigate the association between polymorphisms in the LRP5 gene and bone mineral density (BMD) in three large cohorts of young and elderly men.

METHODS

The cohorts used were MrOS Sweden (n = 3014, aged 69-81 years) and MrOs Hong Kong (n = 2000, aged > 65 years) and the Swedish GOOD study (n = 1068, aged 18-20 years). The polymorphisms Val667Met and Ala1330Val were genotyped using a TaqMan assay.

RESULTS

When combining the data from the Swedish cohorts in a meta-analysis (n = 3,800), men carrying the 667Met-allele had 3% lower BMD at lumbar spine compared with non-carriers (p < 0.05). The Val667Met SNP was not polymorphic in the Hong Kong population and thus were not included. There were no associations between the Ala1330Val SNP and bone phenotypes in the study populations. No associations between the LRP5 polymorphisms and self-reported fractures were seen in MrOs Sweden.

CONCLUSIONS

Results from these three large cohorts indicate that the Val667Met polymorphism but not the Ala1330Val contributes to the observed variability in BMD in the Swedish populations.

Systematic assessment of the human osteoblast transcriptome in resting and induced primary cells

Osteoblasts are key players in bone remodeling. The accessibility of human primary osteoblast-like cells (HObs) from bone explants makes them a lucrative model for studying molecular physiology of bone turnover, for discovering novel anabolic therapeutics, and for mesenchymal cell biology in general. Relatively little is known about resting and dynamic expression profiles of HObs, and to date no studies have been conducted to systematically assess the osteoblast transcriptome. The aim of this study was to characterize HObs and investigate signaling cascades and gene networks with genomewide expression profiling in resting and bone morphogenic protein (BMP)-2- and dexamethasone-induced cells. In addition, we compared HOb gene expression with publicly available samples from the Gene Expression Omnibus. Our data show a vast number of genes and networks expressed predominantly in HObs compared with closely related cells such as fibroblasts or chondrocytes. For instance, genes in the insulin-like growth factor (IGF) signaling pathway were enriched in HObs (P = 0.003) and included the binding proteins (IGFBP-1, -2, -5) and IGF-II and its receptor. Another HOb-specific expression pattern included leptin and its receptor (P < 10(-8)). Furthermore, after stimulation of HObs with BMP-2 or dexamethasone, the expression of several interesting genes and pathways was observed. For instance, our data support the role of peripheral leptin signaling in bone cell function. In conclusion, we provide the landscape of tissue-specific and dynamic gene expression in HObs. This resource will allow utilization of osteoblasts as a model to study specific gene networks and gene families related to human bone physiology and diseases.

IL6 and IL1B polymorphisms are associated with fat mass in older men: the MrOS Study Sweden

There is growing evidence that immune functions are linked to the regulation of body fat. Our studies of knockout mice indicate that both endogenous interleukin (IL)-6 and IL-1 can suppress mature-onset obesity. We now investigated whether four common polymorphisms of the IL6 and IL1 systems are associated with the fat mass measured with dual-energy X-ray absorptiometry (DXA) in elderly men (n = 3,014). The study subjects were from the Swedish part of the MrOS multicenter population study and 69-81 years of age. The IL6 -174 G>C (Minor allele frequency (MAF) = 48%) gene promoter polymorphism was associated with the primary outcome total fat mass (P = 0.006) and regional fat masses, but not with lean body mass. The IL1B -31T>C (MAF = 34%) polymorphism was also associated with total fat (P = 0.007) and regional fat masses, but not lean body mass. The IL-1 receptor antagonist (IL-1ra) gene (IL1RN) +2018 T>C (MAF = 27%) polymorphism (in linkage disequilibrium (LD) with a well-studied variable number tandem repeat of 86 base pair (bp)) and IL1B +3953 C>T (MAF = 26%) polymorphism were not associated with total fat mass. In conclusion, the IL-1 and IL-6 systems, shown to suppress mature-onset obesity in experimental animals, contain gene polymorphisms that are associated with fat, but not lean, mass in elderly men.

Fibroblast growth factor-23 is associated with parathyroid hormone and renal function in a population-based cohort of elderly men

OBJECTIVE

Fibroblast growth factor-23 (FGF23) is a circulating factor involved in phosphate (Pi) and vitamin D metabolism. Serum FGF23 is increased at later stages of chronic kidney disease due to chronic hyperphosphatemia and decreased renal clearance. Recent studies also indicate that FGF23 may directly regulate the expression of parathyroid hormone (PTH) in vitro. Therefore, the objective of the current study was to determine the relationship between FGF23, PTH, and other biochemistries in vivo in subjects with no history of renal disease.

DESIGN

Serum biochemistries were measured in a subsample of the population-based Swedish part of the MrOS study. In total, 1000 Caucasian men aged 70-80 years were randomly selected from the population.

METHODS

Intact FGF23, Pi, calcium, albumin, estimated glomerular filtration rate (eGFR, calculated from cystatin C), PTH, and 25(OH)D3 were measured. Association studies were performed using linear univariate and multivariate regression analyses.

RESULTS

The median FGF23 level was 36.6 pg/ml, ranging from 0.63 to 957 pg/ml. There was a significant correlation between log FGF23 and eGFR (r=-0.21; P<0.00001) and log PTH (r=0.13; P<0.001). These variables remained as independent predictors of FGF23 in multivariate analysis. In addition, log PTH (beta=0.082; P<0.05) and eGFR (beta=-0.090; P<0.05) were associated with log FGF23 in subjects with eGFR>60 ml/min. Only eGFR (beta=-0.35; P<0.0001) remained as a predictor of log FGF23 in subjects with eGFR<60 ml/min.

CONCLUSIONS

Serum FGF23 and PTH are associated in vivo, supporting recent findings that FGF23 directly regulates PTH expression in vitro. Additionally, eGFR is associated with FGF23 in subjects with normal or mildly impaired renal function, indicating that GFR may modulate FGF23 levels independent of serum Pi.

The COMT val158met polymorphism is associated with prevalent fractures in Swedish men

INTRODUCTION

Sex steroids are important for growth and maintenance of the skeleton. Catechol-O-methyltransferase (COMT) is an estrogen degrading enzyme. The COMT val158met polymorphism results in a 60-75% difference in enzyme activity between the val (high activity=H) and met (low activity=L) variants. We have previously reported that this polymorphism is associated with bone mineral density (BMD) in young men. The aim of this study was to investigate associations between COMT val158met, BMD and fractures in elderly men.

METHODS

Population-based study of Swedish men 75.4, SD 3.2, years of age. Fractures were reported using standardized questionnaires. Fracture and genotype data were available from 2,822 individuals.

RESULTS

Total number of individuals with self-reported fracture was 989 (35.0%). Prevalence of >or=1 fracture was 37.2% in COMT(LL), 35.7% in COMT(HL) and 30.4% in COMT(HH) (p<0.05). Early fractures (<or=50 years of age) were less common in COMT(HH) than in the combined COMT(LL+HL) genotype, OR 0.78 (95% CI 0.63-0.97). No associations were found for late fractures (>50 years of age). The OR for fracture of the non-weight bearing skeleton in COMT(HH) compared with COMT(LL+HL) was 0.74 (95% CI 0.59-0.92). No associations between COMT val158met and BMD were found in this cohort of elderly men.

CONCLUSIONS

The COMT val158met polymorphism is associated with life time fracture prevalence in elderly Swedish men. This association is mainly driven by early fractures (<or=50 years of age).

Vitamin D receptor 3' haplotypes are unequally expressed in primary human bone cells and associated with increased fracture risk: the MrOS Study in Sweden and Hong Kong

The VDR is a prime candidate gene for osteoporosis. Here, we studied three common VDR haplotypes in relation to bone phenotypes in 5014 participants of the global MrOS Study. We also studied the relative expression of the haplotypes in human bone cells. One haplotype was associated with increased fracture risk and differently expressed in primary human bone cells.

INTRODUCTION

Vitamin D plays an essential role in skeletal metabolism by binding to its nuclear steroid receptor, the vitamin D receptor (VDR). The heritability of BMD is well established, and the VDR gene is considered a prime candidate suggested to partially account for genetically controlled BMD variance in the population.

MATERIALS AND METHODS

Here, we reconstructed common haplotypes in the VDR 3' untranslated region (UTR) and studied the association to BMD and risk of vertebral fractures in elderly men from Sweden (n = 3014) and Hong Kong (n = 2000), all participants of the global MrOS Study. To assess any functional implications of the VDR polymorphisms, we studied allele-specific expressions of the different VDR 3' UTR haplotypes in the normal chromosomal context of 70 unrelated human trabecular bone samples. This was performed by quantitative genotyping of coding polymorphisms in RNA samples and in corresponding DNA samples isolated from the bone samples.

RESULTS

Three major haplotypes were reconstructed and in agreement with the previously well-defined baT, BAt, and bAT haplotypes, herein denoted Hap1, Hap2, and Hap3. The Hap1 haplotype was independently associated with increased risk of vertebral fractures in Swedish men (OR, 1.655; 95% CI, 1.146-2.391; p < 0.01) and with lower lumbar spine BMD in elderly men from Sweden (p < 0.01) and Hong Kong (p < 0.05). The VDR gene was also shown to exhibit a 3' UTR haplotype dependent allelic imbalance, indicating that the VDR Hap1 allele was overexpressed in human trabecular bone samples.

CONCLUSIONS

The results indicate that the relatively overexpressed VDR Hap1 haplotype could be considered a risk allele for osteoporosis regardless of ethnicity.

The impact of estradiol on bone mineral density is modulated by the specific estrogen receptor-alpha cofactor retinoblastoma-interacting zinc finger protein-1 insertion/deletion polymorphism

CONTEXT

Estrogens regulate bone mass by binding to the estrogen receptor (ER)-alpha as well as ER-beta. The specific ERalpha cofactor retinoblastoma-interacting zinc finger protein (RIZ)-1 enhances ERalpha function in the presence of estrogen.

OBJECTIVE

The objective of the study was to determine whether a RIZ P704 insertion (+)/deletion (-) (indel) polymorphism modulates the impact of estradiol on bone mineral density (BMD) and study the association between the polymorphism and BMD in elderly subjects.

DESIGN

This was a population-based, prospective, and cross-sectional study, the Swedish MrOS Study, and the Malmö OPRA Study, respectively.

SETTING

The study was conducted at three academic medical centers: Sahlgrenska Academy in Gothenburg, Malmö University Hospital, and Uppsala University Hospital.

PARTICIPANTS

In total, 4058 men and women, aged 69-81 yr, were randomly selected from population registries.

MAIN OUTCOME MEASURES

BMD (grams per square centimeter) was measured at femoral neck, trochanter, lumbar spine, and total body.

RESULTS

The RIZ P704(+/+) genotype was associated with low BMD in both women (femoral neck, P < 0.001; trochanter, P < 0.01; lumbar spine, P < 0.05; total body, P < 0.01) and men (lumbar spine, P < 0.05). However, the association between the polymorphism and BMD was dependent on estradiol status. The positive correlation between serum estradiol and BMD was significantly modulated by the genotype with a stronger correlation in the P704(+/+) group than the P704(-/-) group (r = 0.19 vs. r = 0.08, P < 0.05).

CONCLUSIONS

These large-scale studies of elderly men and women indicate that the ERalpha cofactor RIZ gene has a prominent effect on BMD, and the P704 genotype modulates the impact of estradiol on BMD. Further studies are required to determine whether this polymorphism modulates the estrogenic response to estradiol treatment.

Novel estrogen receptor alpha promoter polymorphism increases ventricular hypertrophic response to hypertension

Given the strong genetic contribution to blood pressure and left ventricular hypertrophy (LVH), and the influence of estrogen on these parameters, we hypothesized that polymorphisms in the estrogen receptor alpha (ERalpha) promoter may influence LVH. Three novel polymorphisms were identified upstream of the ERalpha alternatively spliced exon 1E, within sequence which demonstrated significant promoter activity in vitro. Demonstration of ERalpha E isoform expression in human ventricle by RT-PCR supported a possible functional role for the 1E novel polymorphisms in estrogen signaling in the heart. Indeed, G>A (-721 E) was significantly associated with LVH after controlling for systolic blood pressure and sex in a healthy population (n=74), contributing to 23% of interventricular septum (IVS) width variance (p<0.001) and 9.4% of left ventricular mass index (LVMI) variance (p=0.035). In a separate hypertensive cohort, male carriers of the A allele (n=8) had a 17% increase in IVS (95% CI: 6-28%) and a 19% increase in LVMI (3-34%) compared to GG homozygotes (n=84). We conclude that a novel polymorphism in the promoter of a cardiac mRNA splice isoform of ERalpha is associated with LVH.

The positive effect of dietary vitamin D intake on bone mineral density in men is modulated by the polyadenosine repeat polymorphism of the vitamin D receptor

INTRODUCTION

Few studies have considered the dietary influence of vitamin D intake on bone mineral density (BMD). Numerous studies have examined the association between VDR polymorphism and BMD, but no previous study has examined the joint influence of dietary vitamin D intake and VDR polymorphism on BMD.

METHODS

We therefore conducted a study in 230 men aged 41-76 years of age. BMD was measured with DXA. A second bone scan was performed on average 2.7 years after the first investigation. Dietary habits were assessed by 14 dietary 24-h recall interviews. The polyadenosine (A) VDR genotypes were determined.

RESULTS

Dietary vitamin D intake was associated with BMD at all sites, also after multivariate adjustment. Those in the highest quintile of intake had 9% higher femoral neck BMD (p = 0.004), 6% higher BMD at the lumbar spine (p = 0.06) and 5% higher total body BMD (p = 0.003) compared to men in the lowest quintile of dietary vitamin D intake. However, the positive association between vitamin D intake and BMD was especially apparent among those with the L/L polyadenosine (A) VDR genotype explaining between 10 and 15% of the variability in BMD depending on site (p < 0.004). There was furthermore a trend, in the lumbar spine, of less reduction in BMD with increasing vitamin D intake (p = 0.07) but not at the other sites. Calcium intake conferred no association with BMD.

CONCLUSIONS

Our results indicate that the extent of positive association between dietary vitamin D intake and BMD in men is dependent on VDR polymorphism, a novel conceivable important gene-environmental interaction.

A TA-repeat polymorphism in the gene for the estrogen receptor alpha does not correlate with muscle strength or body composition in young adult Swedish women

OBJECTIVES

There are conflicting data in the literature whether estrogens affect muscle strength. Prospective studies with hormone replacement therapy have not been able to convincingly demonstrate a muscular effect and the putative role of estrogen in the development of lean body mass is not established. Both lean mass and fat mass are known to be under strong genetic control and therefore we have investigated the relation between a TA-repeat in the gene for the estrogen receptor alpha (ERalpha) and muscle strength and body composition.

METHODS

175 healthy Swedish women, aged 20-39 were randomly selected from the population registry and included in the study. Body mass measurements (lean mass, fat mass, body weight and BMI) and muscle strength (quadriceps, hamstring and grip strength) were evaluated. The TA-repeat in the ERalpha gene was amplified by polymerase chain reaction.

RESULTS

Alleles with a TA-repeat length of 16 repeats or shorter were denoted short (e), and repeat length of 17 repeats or longer were denoted long (E). Women homozygous for the short and long genotype were denoted ee (31%) and EE (21%), respectively, while heterozygous individuals were denoted Ee (48%). The frequencies were in Hardy-Weinberg equilibrium. No associations were found between ERalpha genotypes and muscle strength or body composition.

CONCLUSION

The TA-repeat in the human ERalpha gene does not correlate with muscle strength or body mass measurements, indicating that body composition is not as sensitive to genetic variation in this receptor as other target organs for estrogen.

A deletion polymorphism in the RIZ gene, a female sex steroid hormone receptor coactivator, exhibits decreased response to estrogen in vitro and associates with low bone mineral density in young Swedish women

Low bone mineral density (BMD) is a major risk factor for osteoporotic fracture, and the trait is under genetic control by a large number of genes. It is recognized that estrogen plays an important role in the maintenance of bone mass by binding to estrogen receptor alpha (ERalpha). RIZ1 has previously been shown to be a specific ERalpha coactivator and strongly enhances its function both in vivo and in vitro. We performed in vitro studies comparing the abilities of RIZ1 P704 polymorphic variants (homozygous presence, P704+; absence, P704-; heterozygosity P704(+/-) of a proline at position 704) to coactivate the ERalpha and also examined the polymorphism associated to BMD of 343 Swedish women, aged 20-39 yr. The expression vector containing P704- RIZ1 showed an impaired response in coactivating ERalpha in a ligand- and dose-dependent manner compared with P704+ RIZ (P < 0.0001). The genotype frequencies were 19% (P704+), 32% (P704-), and 49% (P704(+/-)) and were in Hardy-Weinberg equilibrium. BMD at the heel was higher in the P704+ genotype group than in the P704(+/-) group (P = 0.02), which was evident also after corrections for fat and lean mass (P = 0.03). We conclude that RIZ1 may be a new candidate gene for involvement in the variation seen in BMD.

Genetic variation in the human vitamin D receptor is associated with muscle strength, fat mass and body weight in Swedish women

OBJECTIVE

Bone mineral density (BMD) is under strong genetic control and a number of candidate genes have been associated with BMD. Both muscle strength and body weight are considered to be important predictors of BMD but far less is known about the genes affecting muscle strength and fat mass. The purpose of this study was to investigate the poly adenosine (A) repeat and the BsmI SNP in the vitamin D receptor (VDR) in relation to muscle strength and body composition in healthy women.

DESIGN

A population-based study of 175 healthy women aged 20-39 years was used.

METHODS

The polymorphic regions in the VDR gene (the poly A repeat and the BsmI SNP) were amplified by PCR. Body mass measurements (fat mass, lean mass, body weight and body mass index) and muscle strength (quadriceps, hamstring and grip strength) were evaluated.

RESULTS

Individuals with shorter poly A repeat, ss and/or absence of the linked BsmI restriction site (BB) have higher hamstring strength (ss vs LL, P=0.02), body weight (ss vs LL, P=0.049) and fat mass (ss vs LL, P=0.04) compared with women with a longer poly A repeat (LL) and/or the presence of the linked BsmI restriction site (bb).

CONCLUSIONS

Genetic variation in the VDR is correlated with muscle strength, fat mass and body weight in premenopausal women. Further functional studies on the poly A microsatellite are needed to elucidate whether this is the functionally relevant locus or if the polymorphism is in linkage disequilibrium with a functional variant in a closely situated gene further downstream of the VDR 3'UTR.

A poly adenosine repeat in the human vitamin D receptor gene is associated with bone mineral density in young Swedish women

Peak bone mass (PBM) and subsequent bone loss are important risk factors for development of osteoporosis later in life, and twin studies have reported strong genetic influence on PBM. The genetic factor influencing PBM is polygenetic, and many genes most likely exert relatively small effects on bone mass. The poly adenosine (A) microsatellite in the 3' untranslated region (UTR) of the VDR gene has been associated with both prostate and breast cancer risk but little is known about the effect of bone mineral density (BMD). In this report the poly A microsatellite and the linked BsmI SNP have been investigated in a population-based cohort of 343 Swedish women, aged 20-39. BMD was measured by dual x-ray absorptiometry at the spine, proximal femur, total body and heel and by quantitative ultrasound at the heel. Correlations were found between VDR genotypes and BMD at lumbar spine L2-L4, (ss versus LL, P = 0.03 and BB versus bb, P = 0.02, respectively), with a similar pattern concerning total hip (ss versus LL, P = 0.12 and BB versus bb, P = 0.16 respectively). After corrections for age, height, fat and lean mass, the VDR BsmI genotype was still associated to BMD at the lumbar spine (BB versus bb, P = 0.03). The polymorphisms were in linkage disequilibrium (Chi-square = 566, P < 0.0001). In conclusion, genetic variation in the VDR is associated with BMD in premenopausal women, and further studies are needed to evaluate a possible functional role of the VDR 3'UTR poly A repeat, a region that has shown to be of important for mRNA stability.