Sjur Reppe

Genetic sharing with cardiovascular disease risk factors and diabetes reveals novel bone mineral density loci

Bone Mineral Density (BMD) is a highly heritable trait, but genome-wide association studies have identified few genetic risk factors. Epidemiological studies suggest associations between BMD and several traits and diseases, but the nature of the suggestive comorbidity is still unknown. We used a novel genetic pleiotropy-informed conditional False Discovery Rate (FDR) method to identify single nucleotide polymorphisms (SNPs) associated with BMD by leveraging cardiovascular disease (CVD) associated disorders and metabolic traits. By conditioning on SNPs associated with the CVD-related phenotypes, type 1 diabetes, type 2 diabetes, systolic blood pressure, diastolic blood pressure, high density lipoprotein, low density lipoprotein, triglycerides and waist hip ratio, we identified 65 novel independent BMD loci (26 with femoral neck BMD and 47 with lumbar spine BMD) at conditional FDR < 0.01. Many of the loci were confirmed in genetic expression studies. Genes validated at the mRNA levels were characteristic for the osteoblast/osteocyte lineage, Wnt signaling pathway and bone metabolism. The results provide new insight into genetic mechanisms of variability in BMD, and a better understanding of the genetic underpinnings of clinical comorbidity.

Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus

To identify genetic loci influencing bone accrual, we performed a genome-wide association scan for total-body bone mineral density (TB-BMD) variation in 2,660 children of different ethnicities. We discovered variants in 7q31.31 associated with BMD measurements, with the lowest P = 4.1×10−11 observed for rs917727 with minor allele frequency of 0.37. We sought replication for all SNPs located ±500 kb from rs917727 in 11,052 additional individuals from five independent studies including children and adults, together with de novo genotyping of rs3801387 (in perfect linkage disequilibrium (LD) with rs917727) in 1,014 mothers of children from the discovery cohort. The top signal mapping in the surroundings of WNT16 was replicated across studies with a meta-analysis P = 2.6×10−31 and an effect size explaining between 0.6%–1.8% of TB-BMD variance. Conditional analyses on this signal revealed a secondary signal for total body BMD (P = 1.42×10−10) for rs4609139 and mapping to C7orf58. We also examined the genomic region for association with skull BMD to test if the associations were independent of skeletal loading. We identified two signals influencing skull BMD variation, including rs917727 (P = 1.9×10−16) and rs7801723 (P = 8.9×10−28), also mapping to C7orf58 (r2 = 0.50 with rs4609139). Wnt16 knockout (KO) mice with reduced total body BMD and gene expression profiles in human bone biopsies support a role of C7orf58 and WNT16 on the BMD phenotypes observed at the human population level. In summary, we detected two independent signals influencing total body and skull BMD variation in children and adults, thus demonstrating the presence of allelic heterogeneity at the WNT16 locus. One of the skull BMD signals mapping to C7orf58 is mostly driven by children, suggesting temporal determination on peak bone mass acquisition. Our life-course approach postulates that these genetic effects influencing peak bone mass accrual may impact the risk of osteoporosis later in life.

Levetiracetam, Phenytoin, and Valproate Act Differently on Rat Bone Mass, Structure, and Metabolism

Summary:  Purpose: Long‐term treatment with antiepileptic drugs (AEDs) is associated with increased risk of fractures. Phenytoin (PHT) and valproate (VPA) have both been suggested to influence bone health, whereas levetiracetam (LEV) is scarcely studied. The present study compares the effect of these AEDs on bone mass, biomechanical strength, and bone turnover in rats.

Phenotypic Dissection of Bone Mineral Density Reveals Skeletal Site Specificity and Facilitates the Identification of Novel Loci in the Genetic Regulation of Bone Mass Attainment

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (re = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (re = 0.20–0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n∼9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01×10−37), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31×10−14). In addition, we report a novel association between RIN3 (previously associated with Pagets disease) and LL-BMD (rs754388: β = 0.13, SE = 0.02, P = 1.4×10−10). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Eight genes are highly associated with BMD variation in postmenopausal Caucasian women

Low bone mineral density (BMD) is an important risk factor for skeletal fractures which occur in about 40% of women >/=50 years in the western world. We describe the transcriptional changes in 84 trans-iliacal bone biopsies associated with BMD variations in postmenopausal females (50 to 86 years), aiming to identify genetic determinants of bone structure. The women were healthy or having a primary osteopenic or osteoporotic status with or without low energy fractures. The total cohort of 91 unrelated women representing a wide range of BMDs, were consecutively registered and submitted to global gene Affymetrix microarray expression analysis or histomorphometry. Among almost 23,000 expressed transcripts, a set represented by ACSL3 (acyl-CoA synthetase long-chain family member 3), NIPSNAP3B (nipsnap homolog 3B), DLEU2 (Deleted in lymphocytic leukemia, 2), C1ORF61 (Chromosome 1 open reading frame 61), DKK1 (Dickkopf homolog 1), SOST (Sclerostin), ABCA8, (ATP-binding cassette, sub-family A, member 8), and uncharacterized (AFFX-M27830-M-at), was significantly correlated to total hip BMD (5% false discovery rate) explaining 62% of the BMD variation expressed as T-score, 53% when adjusting for the influence of age (Z-score) and 44% when further adjusting for body mass index (BMI). Only SOST was previously associated to BMD, and the majority of the genes have previously not been associated with a bone phenotype. In molecular network analyses, SOST shows a strong, positive correlation with DKK1, both being members of the Wnt signaling pathway. The results provide novel insight in the underlying biology of bone metabolism and osteoporosis which is the ultimate consequence of low BMD.

Osteopenia, decreased bone formation and impaired osteoblast development in Sox4 heterozygous mice

The transcription factor Sox4 is vital for fetal development, as Sox4–/– homozygotes die in utero. Sox4 mRNA is expressed in the early embryonic growth plate and is regulated by parathyroid hormone, but its function in bone modeling/remodeling is unknown. We report that Sox4+/– mice exhibit significantly lower bone mass (by dual-energy X-ray absorptiometry) from an early age, and fail to obtain the peak bone mass of wild-type (WT) animals. Microcomputed tomography (μCT), histomorphometry and biomechanical testing of Sox4+/– bones show reduced trabecular and cortical thickness, growth plate width, ultimate force and stiffness compared with WT. Bone formation rate (BFR) in 3-month-old Sox4+/– mice is 64% lower than in WT. Primary calvarial osteoblasts from Sox4+/– mice demonstrate markedly inhibited proliferation, differentiation and mineralization. In these cultures, osterix (Osx) and osteocalcin (OCN) mRNA expression was reduced, whereas Runx2 mRNA was unaffected. No functional defects were found in osteoclasts. Silencing of Sox4 by siRNA in WT osteoblasts replicated the defects observed in Sox4+/– cells. We demonstrate inhibited formation and altered microarchitecture of bone in Sox4+/– mice versus WT, without apparent defects in bone resorption. Our results implicate the transcription factor Sox4 in regulation of bone formation, by acting upstream of Osx and independent of Runx2.

Identifying common genetic variants in blood pressure due to polygenic pleiotropy with associated phenotypes

Blood pressure is a critical determinant of cardiovascular morbidity and mortality. It is affected by environmental factors, but has a strong heritable component. Despite recent large genome-wide association studies, few genetic risk factors for blood pressure have been identified. Epidemiological studies suggest associations between blood pressure and several diseases and traits, which may partly arise from a shared genetic basis (genetic pleiotropy). Using genome-wide association studies summary statistics and a genetic pleiotropy-informed conditional false discovery rate method, we systematically investigated genetic overlap between systolic blood pressure (SBP) and 12 comorbid traits and diseases. We found significant enrichment of single nucleotide polymorphisms associated with SBP as a function of their association with body mass index, low-density lipoprotein, waist/hip ratio, schizophrenia, bone mineral density, type 1 diabetes mellitus, and celiac disease. In contrast, the magnitude of enrichment due to shared polygenic effects was smaller with the other phenotypes (triglycerides, high-density lipoproteins, type 2 diabetes mellitus, rheumatoid arthritis, and height). Applying the conditional false discovery rate method to the enriched phenotypes, we identified 62 loci associated with SBP (false discovery rate <0.01), including 42 novel loci. The observed polygenic overlap between SBP and several related disorders indicates that the epidemiological associations are not mediated solely via lifestyle factors but also reflect an etiologic relation that warrants further investigation. The new gene loci identified implicate novel genetic mechanisms related to lipid biology and the immune system in SBP.

Efficient secretion of human parathyroid hormone by Saccharomyces cerevisiae

A cDNA encoding mature human parathyroid hormone (hPTH) was expressed in Saccharomyces cerevisiae, after fusion to the prepro region of yeast mating factor alpha (MF alpha). Radioimmunoassay showed high levels of hPTH immunoreactive material in the growth medium (up to 10 micrograms/ml). More than 95% of the immunoreactive material was found extracellularly as multiple forms of hormone peptides. Three internal cleavage sites were identified in the hPTH molecule. The major cleavage site, after a pair of basic amino acids (aa) (Arg25Lys26 decreases Lys27), resembles that recognized by the KEX2 gene product on which the MF alpha expression-secretion system depends. The use of a protease-deficient yeast strain and the addition of high concentrations of aa to the growth medium, however, not only changed the peptide pattern, but also resulted in a significant increase in the yield of intact hPTH (1-84) (more than 20% of the total amount of immunoreactive material). The secreted hPTH (1-84) migrates like a hPTH standard in two different gel-electrophoretic systems, co-elutes with standard hPTH on reverse-phase high-performance liquid chromatography, reacts with two hPTH antibodies raised against different parts of the peptide, has a correct N-terminal aa sequence, and has full biological activity in a hormone-sensitive osteoblast adenylate cyclase assay.

Methylation of bone SOST, its mRNA, and serum sclerostin levels correlate strongly with fracture risk in postmenopausal women

Inhibition of sclerostin, a glycoprotein secreted by osteocytes, offers a new therapeutic paradigm for treatment of osteoporosis (OP) through its critical role as Wnt/catenin signaling regulator. This study describes the epigenetic regulation of SOST expression in bone biopsies of postmenopausal women. We correlated serum sclerostin to bone mineral density (BMD), fractures, and bone remodeling parameters, and related these findings to epigenetic and genetic disease mechanisms. Serum sclerostin and bone remodeling biomarkers were measured in two postmenopausal groups: healthy (BMD T‐score > –1) and established OP (BMD T‐score < –2.5, with at least one low‐energy fracture). Bone specimens were used to analyze SOST mRNAs, single nucleotide polymorphisms (SNPs), and DNA methylation changes. The SOST gene promoter region showed increased CpG methylation in OP patients (n = 4) compared to age and body mass index (BMI) balanced controls (n = 4) (80.5% versus 63.2%, p = 0.0001) with replication in independent cohorts (n = 27 and n = 36, respectively). Serum sclerostin and bone SOST mRNA expression correlated positively with age‐adjusted and BMI‐adjusted total hip BMD (r = 0.47 and r = 0.43, respectively; both p < 0.0005), and inversely to serum bone turnover markers. Five SNPs, one of which replicates in an independent population‐based genomewide association study (GWAS), showed association with serum sclerostin or SOST mRNA levels under an additive model (p = 0.0016 to 0.0079). Genetic and epigenetic changes in SOST influence its bone mRNA expression and serum sclerostin levels in postmenopausal women. The observations suggest that increased SOST promoter methylation seen in OP is a compensatory counteracting mechanism, which lowers serum sclerostin concentrations and reduces inhibition of Wnt signaling in an attempt to promote bone formation.

Molecular disease map of bone characterizing the postmenopausal osteoporosis phenotype

Genome‐wide gene expressions in bone biopsies from patients with postmenopausal osteoporosis and healthy controls were profiled, to identify osteoporosis candidate genes. All osteoporotic patients (n = 27) in an unbiased cohort of Norwegian women presented with bone mineral density (BMD) T‐scores of less than −2.5 SD and one or more confirmed low‐energy fracture(s). A validation group (n = 18) had clinical and laboratory parameters intermediate to the control (n = 39) and osteoporosis groups. RNA from iliac crest bone biopsies were analyzed by Affymetrix microarrays and real‐time reverse‐transcriptase polymerase chain reaction (RT‐PCR). Differentially expressed genes in osteoporosis versus control groups were identified using the Bayesian ANOVA for microarrays (BAMarray) method, whereas the R‐package Limma (Linear Models for Microarray Data) was used to determine whether these transcripts were explained by disease, age, body mass index (BMI), or combinations thereof. Laboratory tests showed normal ranges for the cohort. A total of 609 transcripts were differentially expressed in osteoporotic patients relative to controls; 256 transcripts were confirmed for disease when controlling for age or BMI. Most of the osteoporosis susceptibility genes (80%) also were confirmed to be regulated in the same direction in the validation group. Furthermore, 217 of 256 transcripts were correlated with BMD (adjusted for age and BMI) at various skeletal sites (|r| > 0.2, p < .05). Among the most distinctly expressed genes were Wnt antagonists DKK1 and SOST, the transcription factor SOX4, and the bone matrix proteins MMP13 and MEPE, all reduced in osteoporosis versus control groups. Our results identify potential osteoporosis susceptibility candidate genes adjusted for confounding factors (ie, age and BMI) with or without a significant correlation with BMD.