Global miRNA expression and correlation with mRNA levels in primary human bone cells.

MicroRNAs (miRNAs) are important post-transcriptional regulators that have recently introduced an additional level of intricacy to our understanding of gene regulation. The aim of this study was to investigate miRNA-mRNA interactions that may be relevant for bone metabolism by assessing correlations and interindividual variability in miRNA levels as well as global correlations between miRNA and mRNA levels in a large cohort of primary human osteoblasts (HOBs) obtained during orthopedic surgery in otherwise healthy individuals. We identified differential expression (DE) of 24 miRNAs, and found 9 miRNAs exhibiting DE between males and females. We identified hsa-miR-29b, hsa-miR-30c2, and hsa-miR-125b and their target genes as important modulators of bone metabolism. Further, we used an integrated analysis of global miRNA-mRNA correlations, mRNA-expression profiling, DE, bioinformatics analysis, and functional studies to identify novel target genes for miRNAs with the potential to regulate osteoblast differentiation and extracellular matrix production. Functional studies by overexpression and knockdown of miRNAs showed that, the differentially expressed miRNAs hsa-miR-29b, hsa-miR-30c2, and hsa-miR-125b target genes highly relevant to bone metabolism, e.g., collagen, type I, α1 (COL1A1), osteonectin (SPARC), Runt-related transcription factor 2 (RUNX2), osteocalcin (BGLAP), and frizzled-related protein (FRZB). These miRNAs orchestrate the activities of key regulators of osteoblast differentiation and extracellular matrix proteins by their convergent action on target genes and pathways to control the skeletal gene expression.

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 = 24 h), and PGE2 (t = 24 h). 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.

Allele dependent silencing of collagen type I using small interfering RNAs targeting 3'UTR Indels - a novel therapeutic approach in osteogenesis imperfecta.

Osteogenesis imperfecta, also known as "brittle bone disease", is a heterogeneous disorder of connective tissue generally caused by dominant mutations in the genes COL1A1 and COL1A2, encoding the α1 and α2 chains of type I (pro)collagen. Symptomatic patients are usually prescribed bisphosphonates, but this treatment is neither curative nor sufficient. A promising field is gene silencing through RNA interference. In this study small interfering RNAs (siRNAs) were designed to target each allele of 3'UTR insertion/deletion polymorphisms (indels) in COL1A1 (rs3840870) and COL1A2 (rs3917). For both indels, the frequency of heterozygous individuals was determined to be approximately 50% in Swedish cohorts of healthy controls as well as in patients with osteogenesis imperfecta. Cultures of primary human bone derived cells were transfected with siRNAs through magnet-assisted transfection. cDNA from transfected cells was sequenced in order to measure targeted allele/non-targeted allele ratios and the overall degree of silencing was assessed by quantitative PCR. Successful allele dependent silencing was observed, with promising results for siRNAs complementary to both the insertion and non-insertion harboring alleles. In COL1A1 cDNA the indel allele ratios were shifted from 1 to 0.09 and 0.19 for the insertion and non-insertion allele respectively while the equivalent resulting ratios for COL1A2 were 0.05 and 0.01. Reductions in mRNA abundance were also demonstrated; in cells treated with siRNAs targeting the COL1A1 alleles the average COL1A1 mRNA levels were reduced 65% and 78% compared to negative control levels and in cells treated with COL1A2 siRNAs the average COL1A2 mRNA levels were decreased 26% and 49% of those observed in the corresponding negative controls. In conclusion, allele dependent silencing of collagen type I utilizing 3'UTR indels common in the general population constitutes a promising mutation independent therapeutic approach for osteogenesis imperfecta.

Genome-wide association meta-analysis of cortical bone mineral density unravels allelic heterogeneity at the RANKL locus and potential pleiotropic effects on bone.

Previous genome-wide association (GWA) studies have identified SNPs associated with areal bone mineral density (aBMD). However, this measure is influenced by several different skeletal parameters, such as periosteal expansion, cortical bone mineral density (BMD(C)) cortical thickness, trabecular number, and trabecular thickness, which may be under distinct biological and genetic control. We have carried out a GWA and replication study of BMD(C), as measured by peripheral quantitative computed tomography (pQCT), a more homogenous and valid measure of actual volumetric bone density. After initial GWA meta-analysis of two cohorts (ALSPAC n = 999, aged ∼15 years and GOOD n = 935, aged ∼19 years), we attempted to replicate the BMD(C) associations that had p<1×10(-5) in an independent sample of ALSPAC children (n = 2803) and in a cohort of elderly men (MrOS Sweden, n = 1052). The rs1021188 SNP (near RANKL) was associated with BMD(C) in all cohorts (overall p = 2×10(-14), n = 5739). Each minor allele was associated with a decrease in BMD(C) of ∼0.14SD. There was also evidence for an interaction between this variant and sex (p = 0.01), with a stronger effect in males than females (at age 15, males -6.77mg/cm(3) per C allele, p = 2×10(-6); females -2.79 mg/cm(3) per C allele, p = 0.004). Furthermore, in a preliminary analysis, the rs1021188 minor C allele was associated with higher circulating levels of sRANKL (p<0.005). We show this variant to be independent from the previously aBMD associated SNP (rs9594738) and possibly from a third variant in the same RANKL region, which demonstrates important allelic heterogeneity at this locus. Associations with skeletal parameters reflecting bone dimensions were either not found or were much less pronounced. This finding implicates RANKL as a locus containing variation associated with volumetric bone density and provides further insight into the mechanism by which the RANK/RANKL/OPG pathway may be involved in skeletal development.

Genetics of tibia bone properties of crossbred commercial laying hens in different housing systems.

Osteoporosis and bone fractures are a severe problem for the welfare of laying hens, with genetics and environment, such as housing system, each making substantial contributions to bone strength. In this work, we performed genetic analyses of bone strength, bone mineral density, and bone composition, as well as body weight, in 860 commercial crossbred laying hens from 2 different companies, kept in either furnished cages or floor pens. We compared bone traits between housing systems and crossbreds and performed a genome-wide association study of bone properties and body weight. As expected, the 2 housing systems produced a large difference in bone strength, with layers housed in floor pens having stronger bones. These differences were accompanied by differences in bone geometry, mineralization, and chemical composition. Genome scans either combining or independently analyzing the 2 housing systems revealed no genome-wide significant loci for bone breaking strength. We detected 3 loci for body weight that were shared between the housing systems on chromosomes 4, 6, and 27 (either genome-wide significant or suggestive) and these coincide with associations for bone length. In summary, we found substantial differences in bone strength, content, and composition between hens kept in floor pens and furnished cages that could be attributed to greater physical activity in pen housing. We found little evidence for large-effect loci for bone strength in commercial crossbred hens, consistent with a highly polygenic architecture for bone strength in the production environment. The lack of consistent genetic associations between housing systems in combination with the differences in bone phenotypes could be due to gene-by-environment interactions with housing system or a lack of power to detect shared associations for bone strength.

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.

Common local founder effects for Wilson's disease and hereditary hemochromatosis; mutation studies of a large family.

UNLABELLED: Wilson's disease (WND) and hereditary hemochromatosis (HH) are two metal loading diseases of copper and iron, respectively, and are both recessively inherited. In central Sweden, where HH is common, 9 Wilson kindred (14 members) were identified. Aims of the study were to test whether nine WND families shared a common origin, a common mutation and if carrying HFE mutations affected their phenotype. RESULTS: The nine families were traced through 13 generations to a common founder origin in the mid-seventeenth century. Despite identity of descent, four different ATP7B mutations appeared with homozygosity in four, with two different mutations, W779X and T977M. There were three compound heterozygotes, W779X/T977M, R1319X/H1069Q and one T977M combined with a new, previously not described mutation, probably of Finnish origin. The founder family also included 26 descendant kindred (55 members) with HH as shown by HFE mutations. This admixture coincided with a migration out of the original parish into hemochromatosis-rich localities. One WND patient had iron overload (serum ferritin 672 µg/l and raised liver enzymes), but lacked HFE mutations. In another family with serious hemochromatosis (two sons dying from bronze diabetes), the coinheritance of congenital spherocytosis was probably the cause rather than an additional effect of WND. CONCLUSIONS: WND though a rare disease may become aggregated like HH in certain areas due to local founder effects. Despite extensive pedigree studies leading back to the local founder family, the authors were unable to find a single defining mutation of the ATP7B gene.

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.

Real-World Effectiveness of Anti-Resorptive Treatment in Patients With Incident Fragility Fractures-The STORM Cohort-A Swedish Retrospective Observational Study.

Results from real-world evidence (RWE) from the largest healthcare region in Sweden show low uptake of antiresorptive (AR) treatment, but beneficial effect in those receiving treatment, especially for the composite outcome of hip fracture or death. For RWE studies, Sweden is unique, with virtually complete coverage of electronic medical records (EMRs) and both regional and national registries, in a universal publicly funded healthcare system. To our knowledge, there is no previous RWE study evaluating the efficacy of AR treatment compared to no AR treatment after fragility fracture, including data on parenteral treatments administered in hospital settings. The Stockholm Real World Management (STORM) study cohort was established in the healthcare region of Stockholm to retrospectively assess the effectiveness of AR treatment after first fragility fracture using the regional EMR system for both hospital and primary care. Between 2012 and 2018, we identified 69,577 fragility fracture episodes among 59,078 patients, men and women, 50 years and older. Of those, 21,141 patients met inclusion and exclusion criteria (eligible cohort). From these, the final matched study cohort comprised 9840 fragility fractures (cases receiving AR treatment [n = 1640] and controls not receiving AR treatment [n = 8200]). Propensity scores were estimated using logistic regression models with AR treatment as outcome and confounders as independent variables followed by analysis using Cox proportional hazard models. Real world evidence from Sweden's largest healthcare region, comprising a quarter of the Swedish population, show that only 10% of patients receive AR treatment within 1 year after a fragility fracture. Factors associated with not receiving treatment include having a diagnosis of cardiovascular disease. In those treated, AR have positive effects particularly on the composite of fracture and death (any fracture/death and hip fracture/death) in individuals matched for all major confounders. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

COL1 C-propeptide cleavage site mutations cause high bone mass osteogenesis imperfecta.

Osteogenesis imperfecta (OI) is most often caused by mutations in the type I procollagen genes (COL1A1/COL1A2). We identified two children with substitutions in the type I procollagen C-propeptide cleavage site, which disrupt a unique processing step in collagen maturation and define a novel phenotype within OI. The patients have mild OI caused by mutations in COL1A1 (Patient 1: p.Asp1219Asn) or COL1A2 (Patient 2: p.Ala1119Thr), respectively. Patient 1 L1-L4 DXA Z-score was +3.9 and pQCT vBMD was+3.1; Patient 2 had L1-L4 DXA Z-score of 0.0 and pQCT vBMD of -1.8. Patient BMD contrasts with radiographic osteopenia and histomorphometry without osteosclerosis. Mutant procollagen processing is impaired in pericellular and in vitro assays. Patient dermal collagen fibrils have irregular borders. Incorporation of pC-collagen into matrix leads to increased bone mineralization. FTIR imaging confirms elevated mineral/matrix ratios in both patients, along with increased collagen maturation in trabecular bone, compared to normal or OI controls. Bone mineralization density distribution revealed a marked shift toward increased mineralization density for both patients. Patient 1 has areas of higher and lower bone mineralization than controls; Patient 2's bone matrix has a mineral content exceeding even classical OI bone. These patients define a new phenotype of high BMD OI and demonstrate that procollagen C-propeptide cleavage is crucial to normal bone mineralization.

Hypervariable intronic region in NCX1 is enriched in short insertion-deletion polymorphisms and showed association with cardiovascular traits.

BACKGROUND: Conserved non-coding regions (CNR) have been shown to harbor gene expression regulatory elements. Genetic variations in these regions may potentially contribute to complex disease susceptibility. METHODS: We targeted CNRs of cardiovascular disease (CVD) candidate gene, Na(+)-Ca(2+) exchanger (NCX1) with polymorphism screening among CVD patients (n = 46) using DHPLC technology. The flanking region (348 bp) of the 14 bp indel in intron 2 was further genotyped by DGGE assay in two Eastern-European CVD samples: essential hypertension (HYPEST; 470 cases, 652 controls) and coronary artery disease, CAD (CADCZ; 257 cases, controls 413). Genotype-phenotype associations were tested by regression analysis implemented in PLINK. Alignments of primate sequences were performed by ClustalW2. RESULTS: Nine of the identified NCX1 variants were either singletons or targeted by commercial platforms. The 14 bp intronic indel (rs11274804) was represented with substantial frequency in HYPEST (6.82%) and CADCZ (14.58%). Genotyping in Eastern-Europeans (n = 1792) revealed hypervariable nature of this locus, represented by seven alternative alleles. The alignments of human-chimpanzee-macaque sequences showed that the major human variant (allele frequency 90.45%) was actually a human-specific deletion compared to other primates. In humans, this deletion was surrounded by other short (5-43 bp) deletion variants and a duplication (40 bp) polymorphism possessing overlapping breakpoints. This indicates a potential indel hotspot, triggered by the initial deletion in human lineage. An association was detected between the carrier status of 14 bp indel ancestral allele and CAD (P = 0.0016, OR = 2.02; Bonferroni significance level alpha = 0.0045), but not with hypertension. The risk for the CAD development was even higher among the patients additionally diagnosed with metabolic syndrome (P = 0.0014, OR = 2.34). Consistent with the effect on metabolic processes, suggestive evidence for the association with heart rate, serum triglyceride and LDL levels was detected (P = 0.04). CONCLUSIONS: Compared to SNPs targeted by large number of locus-specific and genome-wide assays, considerably less attention has been paid to short indel variants in the human genome. The data of genome dynamics, mutation rate and population genetics of short indels, as well as their impact on gene expressional profile and human disease susceptibility is limited. The characterization of NCX1 intronic hypervariable non-coding region enriched in human-specific indel variants contributes to this gap of knowledge.

The association between Single Nucleotide Polymorphisms of Klotho Gene and Mortality in Elderly Men: The MrOS Sweden Study.

The Klotho (KL) gene is involved in phosphate homeostasis. Polymorphisms in this gene have been reported to be associated with the risk of cardiovascular disease. Here we used computational tools to predict the damage-associated single nucleotide polymorphisms (SNPs) in the human KL gene. We further investigated the association of SNPs in the KL gene and mortality in the Swedish multicenter prospective Osteoporotic Fractures in Men (MrOS) cohort. This study included 2921 men (aged 69-81 years) with mean 4.49 ± 1.03 years follow-up. 18 SNPs in the KL gene were genotyped using Sequenom. These SNPs were identified by in silico tools for the coding and noncoding genome to predict the damaging SNPs. After quality analyses, SNPs were analyzed for mortality risk using two steps approach on logistic regression model screening and then Cox regression model confirmation. Two non-synonymous SNPs rs9536314 and rs9527025 were found to be potentially damaging SNPs that affect KL protein stability and expression. However, these two SNPs were not statistically significantly associated with all-cause mortality (crude Hazard ratio [HR] 1.72, 95% confidence interval [CI] 0.96-3.07 in rs9536314; crude HR 1.82, 95% CI 0.998-3.33 in rs9527025) or cardiovascular mortality (crude HR 1.52, 95% CI 0.56-4.14 in rs9536314; crude HR 1.54, 95% CI 0.55-4.33 in rs9527025) in additive model using Cox regression analysis. In conclusion, these two potentially damaging SNPs (rs9536314 and rs9527025) in the KL gene were not associated with all-cause mortality or cardiovascular mortality in MrOs cohort. Larger scales studies and meta-analysis are needed to confirm the correlation between polymorphisms of the KL gene and mortality.

Heterozygosity for a coding SNP in COL1A2 confers a lower BMD and an increased stroke risk.

Genetic variation plays an important role in osteoporosis and a prime candidate gene is Collagen alpha2(I) (COL1A2). A coding polymorphism (rs42524) in COL1A2 has previously been associated with intracranial aneurysms. Here the effects of this polymorphism have been studied in relation to bone mineral density (BMD) and prevalences of stroke and myocardial infarction (MI). rs42524 was genotyped in elderly men (n = 2004) from the Swedish MrOS cohort. Genotypes were analysed for association to BMD and certain health parameters. Significant associations (overall P < 0.05), were observed between rs42524 genotype and BMD at several skeletal sites. Surprisingly, the heterozygote genotype class exhibited lower BMD than either homozygote group. When subjects were classified as heterozygotes or homozygotes, the heterozygous genotype was found to confer a lower BMD at total hip, femoral neck and trochanter Furthermore, the heterozygote genotype had an increased risk of stroke and MI, with population Attributable Risks being 0.12 and 0.08, respectively.

Expanding the phenotypic spectrum of osteogenesis imperfecta type V including heterotopic ossification of muscle origins and attachments.

BACKGROUND: Osteogenesis imperfecta (OI) is a clinical and genetic heterogeneous group of connective tissue disorders, characterized by bone fragility and a propensity to fracture. METHODS: In this report we describe the clinical phenotype of two patients, a 28-year-old woman and her mother (54 years old), both with a history of short stature and multiple fractures. RESULTS: Exome sequencing revealed the recurring IFITM5:c.-14 C>T variant causing OI type V. Both patients had several fractures during childhood. CT-scan and scintigraphy showed ossification of the origin and attachment of muscles and hypertrophic callus formation. CONCLUSION: Ossification of the origin and attachment of muscles seems to be part of the phenotype in patients with OI type V.

Gaucher Disease in Bone: From Pathophysiology to Practice.

Gaucher disease (GD) is a rare, genetic lysosomal disorder leading to lipid accumulation and dysfunction in multiple organs. Involvement of the skeleton is one of the most prevalent aspects of GD and a major cause of pain, disability, and reduced quality of life. Uniform recommendations for contemporary evaluation and management are needed. To develop practical clinical recommendations, an international group of experienced physicians conducted a comprehensive review of 20 years' of the literature, defining terms according to pathophysiological understanding and pointing out best practice and unmet needs related to the skeletal features of this disorder. Abnormalities of bone modeling, reduced bone density, bone infarction, and plasma cell dyscrasias accompany the displacement of healthy adipocytes in adult marrow. Exposure to excess bioactive glycosphingolipids appears to affect hematopoiesis and the balance of osteoblast and osteoclast numbers and activity. Imbalance between bone formation and breakdown induces disordered trabecular and cortical bone modeling, cortical bone thinning, fragility fractures, and osteolytic lesions. Regular assessment of bone mineral density, marrow infiltration, the axial skeleton and searching for potential malignancy are recommended. MRI is valuable for monitoring skeletal involvement: It provides semiquantitative assessment of marrow infiltration and the degree of bone infarction. When MRI is not available, monitoring of painful acute bone crises and osteonecrosis by plain X-ray has limited value. In adult patients, we recommend DXA of the lumbar spine and left and right hips, with careful protocols designed to exclude focal disease; serial follow-up should be done using the same standardized instrument. Skeletal health may be improved by common measures, including adequate calcium and vitamin D and management of pain and orthopedic complications. Prompt initiation of specific therapy for GD is crucial to optimizing outcomes and preventing irreversible skeletal complications. Investing in safe, clinically useful, and better predictive methods for determining bone integrity and fracture risk remains a need. © 2019 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals Inc.

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.

The genetic architecture of a female sexual ornament.

Understanding the evolution of sexual ornaments, and particularly that of female sexual ornaments, is an enduring challenge in evolutionary biology. Key to this challenge are establishing the relationship between ornament expression and female reproductive investment, and determining the genetic basis underpinning such relationship. Advances in genomics provide unprecedented opportunities to study the genetic architecture of sexual ornaments in model species. Here, we present a quantitative trait locus (QTL) analysis of a female sexual ornament, the comb of the fowl, Gallus gallus, using a large-scale intercross between red junglefowl and a domestic line, selected for egg production. First, we demonstrate that female somatic investment in comb reflects female reproductive investment. Despite a trade-off between reproductive and skeletal investment mediated by the mobilization of skeletal minerals for egg production, females with proportionally large combs also had relatively high skeletal investment. Second, we identify a major QTL for bisexual expression of comb mass and several QTL specific to female comb mass. Importantly, QTL for comb mass were nonrandomly clustered with QTL for female reproductive and skeletal investment on chromosomes one and three. Together, these results shed light onto the physiological and genetic architecture of a female ornament.

Estrogen receptor alpha gene polymorphism and endometrial cancer risk--a case-control study.

BACKGROUND: Estrogen is an established endometrial carcinogen. One of the most important mediators of estrogenic action is the estrogen receptor alpha. We have investigated whether polymorphic variation in the estrogen receptor alpha gene (ESR1) is associated with endometrial cancer risk. METHODS: In 702 cases with invasive endometrial cancer and 1563 controls, we genotyped five markers in ESR1 and used logistic regression models to estimate odds ratios (OR) and 95 percent confidence intervals (CI). RESULTS: We found an association between rs2234670, rs2234693, as well as rs9340799, markers in strong linkage disequilibrium (LD), and endometrial cancer risk. The association with rs9340799 was the strongest, OR 0.75 (CI 0.60-0.93) for heterozygous and OR 0.53 (CI 0.37-0.77) for homozygous rare compared to those homozygous for the most common allele. Haplotype models did not fit better to the data than single marker models. CONCLUSION: We found that intronic variation in ESR1 was associated with endometrial cancer risk.

Allele dependent silencing of COL1A2 using small interfering RNAs.

Osteogenesis imperfecta (OI) is generally caused by a dominant mutation in Collagen I, encoded by the genes COL1A1 and COL1A2. To date there is no satisfactory therapy for OI, but inactivation of the mutant allele through small interfering RNAs (siRNA) is a promising approach, as siRNAs targeting each allele of a polymorphism could be used for allele-specific silencing irrespective of the location of the actual mutations. In this study we examined the allele dependent effects of several tiled siRNAs targeting a region surrounding an exonic COL1A2 T/C polymorphism (rs1800222) in heterozygous primary human bone cells. Relative abundances of COL1A2 alleles were determined by cDNA sequencing and overall COL1A2 abundance was analyzed by quantitative PCR. One of the siRNAs decreased overall COL1A2 abundance by 71% of which 75% was due to silencing of the targeted T-allele. In conclusion, allele-preferential silencing of Collagen type I genes may be a future therapeutic approach for OI.

Haplotypes in the CYP2R1 gene are associated with levels of 25(OH)D and bone mineral density, but not with other markers of bone metabolism (MrOS Sweden).

OBJECTIVE: Polymorphisms in the CYP2R1 gene encoding Vitamin D 25-hydroxylase have been reported to correlate with circulating levels of 25-OH vitamin D3 (25(OH)D). It is unknown whether these variations also affect overall bone metabolism. In order to elucidate the overall associations of polymorphisms in the CYP2R1, we studied haplotype tagging single nucleotide polymorphisms (SNPs) in the gene and serum levels of 25(OH)D, calcium, phosphate, parathyroid hormone (PTH) and fibroblast growth factor-23 (FGF23), as well as bone mineral density (BMD). METHODS: Baseline data on serum parameters and BMD from MrOS Sweden, a prospective population-based cohort study of elderly men (mean age 75 years, range 69-81), were analyzed. Genotyping was performed for eight SNPs covering the CYP2R1 gene in 2868 men with available samples of DNA. Subjects were followed up concerning incidence of fracture during five years. RESULTS: There was a significant genetic association with circulating levels of 25(OH)D (4.6-18.5% difference in mean values between SNP alleles), but there were no correlations with levels of calcium, phosphate, PTH or FGF23 for any genetic variant. No differences were found in fracture incidence between the variants. There was an inverse relationship between lower BMD and concomitant higher 25(OH)D for three of the haplotypes (p < 0.005). CONCLUSIONS: Common variants in the CYP2R1 gene encoding Vitamin D 25-hydroxylase correlate with levels of circulating 25(OH)D but do not otherwise associate with measures of calcium and phosphate homeostasis. Presence of the specific haplotypes may be an indicator of risk for low 25(OH)D levels, and may in addition be correlated to bone mineral density.