Tilting the lasso by knowledge-based post-processing

Background It is useful to incorporate biological knowledge on the role of genetic determinants in predicting an outcome. It is, however, not always feasible to fully elicit this information when the number of determinants is large. We present an approach to overcome this difficulty. First, using half of the available data, a shortlist of potentially interesting determinants are generated. Second, binary indications of biological importance are elicited for this much smaller number of determinants. Third, an analysis is carried out on this shortlist using the second half of the data. Results We show through simulations that, compared with adaptive lasso, this approach leads to models containing more biologically relevant variables, while the prediction mean squared error (PMSE) is comparable or even reduced. We also apply our approach to bone mineral density data, and again final models contain more biologically relevant variables and have reduced PMSEs. Conclusion Our method leads to comparable or improved predictive performance, and models with greater face validity and interpretability with feasible incorporation of biological knowledge into predictive models

In silico discovery of blood cell macromolecular associations

Background Physical molecular interactions are the basis of intracellular signalling and gene regulatory networks, and comprehensive, accessible databases are needed for their discovery. Highly correlated transcripts may reflect important functional associations, but identification of such associations from primary data are cumbersome. We have constructed and adapted a user-friendly web application to discover and identify putative macromolecular associations in human peripheral blood based on significant correlations at the transcriptional level. Methods The blood transcriptome was characterized by quantification of 17,328 RNA species, including 341 mature microRNAs in 105 clinically well-characterized postmenopausal women. Intercorrelation of detected transcripts signal levels generated a matrix with > 150 million correlations recognizing the human blood RNA interactome. The correlations with calculated adjusted p-values were made easily accessible by a novel web application. Results We found that significant transcript correlations within the giant matrix reflect experimentally documented interactions involving select ubiquitous blood relevant transcription factors (CREB1, GATA1, and the glucocorticoid receptor (GR, NR3C1)). Their responsive genes recapitulated up to 91% of these as significant correlations, and were replicated in an independent cohort of 1204 individual blood samples from the Framingham Heart Study. Furthermore, experimentally documented mRNAs/miRNA associations were also reproduced in the matrix, and their predicted functional co-expression described. The blood transcript web application is available at http://app.uio.no/med/klinmed/correlation-browser/blood/index.php and works on all commonly used internet browsers. Conclusions Using in silico analyses and a novel web application, we found that correlated blood transcripts across 105 postmenopausal women reflected experimentally proven molecular associations. Furthermore, the associations were reproduced in a much larger and more heterogeneous cohort and should therefore be generally representative. The web application lends itself to be a useful hypothesis generating tool for identification of regulatory mechanisms in complex biological data sets.publishedVersio

Identification of Transcripts with Shared Roles in the Pathogenesis of Postmenopausal Osteoporosis and Cardiovascular Disease

Epidemiological evidence suggests existing comorbidity between postmenopausal osteoporosis (OP) and cardiovascular disease (CVD), but identification of possible shared genes is lacking. The skeletal global transcriptomes were analyzed in trans-iliac bone biopsies (n = 84) from clinically well-characterized postmenopausal women (50 to 86 years) without clinical CVD using microchips and RNA sequencing. One thousand transcripts highly correlated with areal bone mineral density (aBMD) were further analyzed using bioinformatics, and common genes overlapping with CVD and associated biological mechanisms, pathways and functions were identified. Fifty genes (45 mRNAs, 5 miRNAs) were discovered with established roles in oxidative stress, inflammatory response, endothelial function, fibrosis, dyslipidemia and osteoblastogenesis/calcification. These pleiotropic genes with possible CVD comorbidity functions were also present in transcriptomes of microvascular endothelial cells and cardiomyocytes and were differentially expressed between healthy and osteoporotic women with fragility fractures. The results were supported by a genetic pleiotropy-informed conditional False Discovery Rate approach identifying any overlap in single nucleotide polymorphisms (SNPs) within several genes encoding aBMD- and CVD-associated transcripts. The study provides transcriptional and genomic evidence for genes of importance for both BMD regulation and CVD risk in a large collection of postmenopausal bone biopsies. Most of the transcripts identified in the CVD risk categories have no previously recognized roles in OP pathogenesis and provide novel avenues for exploring the mechanistic basis for the biological association between CVD and OP.</p

Identification of Transcripts with Shared Roles in the Pathogenesis of Postmenopausal Osteoporosis and Cardiovascular Disease

Epidemiological evidence suggests existing comorbidity between postmenopausal osteoporosis (OP) and cardiovascular disease (CVD), but identification of possible shared genes is lacking. The skeletal global transcriptomes were analyzed in trans-iliac bone biopsies (n = 84) from clinically well-characterized postmenopausal women (50 to 86 years) without clinical CVD using microchips and RNA sequencing. One thousand transcripts highly correlated with areal bone mineral density (aBMD) were further analyzed using bioinformatics, and common genes overlapping with CVD and associated biological mechanisms, pathways and functions were identified. Fifty genes (45 mRNAs, 5 miRNAs) were discovered with established roles in oxidative stress, inflammatory response, endothelial function, fibrosis, dyslipidemia and osteoblastogenesis/calcification. These pleiotropic genes with possible CVD comorbidity functions were also present in transcriptomes of microvascular endothelial cells and cardiomyocytes and were differentially expressed between healthy and osteoporotic women with fragility fractures. The results were supported by a genetic pleiotropy-informed conditional False Discovery Rate approach identifying any overlap in single nucleotide polymorphisms (SNPs) within several genes encoding aBMD- and CVD-associated transcripts. The study provides transcriptional and genomic evidence for genes of importance for both BMD regulation and CVD risk in a large collection of postmenopausal bone biopsies. Most of the transcripts identified in the CVD risk categories have no previously recognized roles in OP pathogenesis and provide novel avenues for exploring the mechanistic basis for the biological association between CVD and OP.</p

Genetic polymorphism of miR-196a-2 is associated with bone mineral density (BMD)

MicroRNAs (miRNAs) are small non-coding RNA molecules that post-transcriptionally regulate the translation of messenger RNAs. Given the crucial role of miRNAs in gene expression, genetic variants within miRNA-related sequences may affect miRNA function and contribute to disease risk. Osteoporosis is characterized by reduced bone mass, and bone mineral density (BMD) is a major diagnostic proxy to assess osteoporosis risk. Here, we aimed to identify miRNAs that are involved in BMD using data from recent genome-wide association studies (GWAS) on femoral neck, lumbar spine and forearm BMD. Of 242 miRNA-variants available in the GWAS data, we found rs11614913:C > T in the precursor miR-196a-2 to be significantly associated with femoral neck-BMD (p-value = 9.9 × 10-7, β = −0.038) and lumbar spine-BMD (p-value = 3.2 × 10-11, β = −0.061). Furthermore, our sensitivity analyses using the Rotterdam study data showed a sex-specific association of rs11614913 with BMD only in women. Subsequently, we highlighted a number of miR-196a-2 target genes, expressed in bone and associated with BMD, that may mediate the miRNA function in BMD. Collectively, our results suggest that miR-196a-2 may contribute to variations in BMD level. Further biological investigations will give more insights into the mechanisms by which miR-196a-2 control expression of BMD-related genes

Cytokines Explored in Saliva and Tears from Radiated Cancer Patients Correlate with Clinical Manifestations, Influencing Important Immunoregulatory Cellular Pathways

publishedVersio

LINE-1 RNA triggers matrix formation in bone cells via a PKR-mediated inflammatory response

Transposable elements (TEs) are mobile genetic modules of viral derivation that have been co-opted to become modulators of mammalian gene expression. TEs are a major source of endogenous dsRNAs, signaling molecules able to coordinate inflammatory responses in various physiological processes. Here, we provide evidence for a positive involvement of TEs in inflammation-driven bone repair and mineralization. In newly fractured mice bone, we observed an early transient upregulation of repeats occurring concurrently with the initiation of the inflammatory stage. In human bone biopsies, analysis revealed a significant correlation between repeats expression, mechanical stress and bone mineral density. We investigated a potential link between LINE-1 (L1) expression and bone mineralization by delivering a synthetic L1 RNA to osteoporotic patient-derived mesenchymal stem cells and observed a dsRNA-triggered protein kinase (PKR)-mediated stress response that led to strongly increased mineralization. This response was associated with a strong and transient inflammation, accompanied by a global translation attenuation induced by eIF2α phosphorylation. We demonstrated that L1 transfection reshaped the secretory profile of osteoblasts, triggering a paracrine activity that stimulated the mineralization of recipient cells

Distinct Subsets of Noncoding RNAs Are Strongly Associated With BMD and Fracture, Studied in Weight-Bearing and Non–Weight-Bearing Human Bone

We investigated mechanisms resulting in low bone mineral density (BMD) and susceptibility to fracture by comparing noncoding RNAs (ncRNAs) in biopsies of non–weight-bearing (NWB) iliac (n = 84) and weight bearing (WB) femoral (n = 18) postmenopausal bone across BMDs varying from normal (T-score > −1.0) to osteoporotic (T-score ≤ −2.5). Global bone ncRNA concentrations were determined by PCR and microchip analyses. Association with BMD or fracture, adjusted by age and body mass index, were calculated using linear and logistic regression and least absolute shrinkage and selection operator (Lasso) analysis. At 10% false discovery rate (FDR), 75 iliac bone ncRNAs and 94 femoral bone ncRNAs were associated with total hip BMD. Eight of the ncRNAs were common for the two sites, but five of them (miR-484, miR-328-3p, miR-27a-5p, miR-28-3p, and miR-409-3p) correlated positively to BMD in femoral bone, but negatively in iliac bone. Of predicted pathways recognized in bone metabolism, ECM-receptor interaction and prote

Identification of a novel locus on chromosome 2q13, which predisposes to clinical vertebral fractures independently of bone density

OBJECTIVES: To identify genetic determinants of susceptibility to clinical vertebral fractures, which is an important complication of osteoporosis. METHODS: Here we conduct a genome-wide association study in 1553 postmenopausal women with clinical vertebral fractures and 4340 controls, with a two-stage replication involving 1028 cases and 3762 controls. Potentially causal variants were identified using expression quantitative trait loci (eQTL) data from transiliac bone biopsies and bioinformatic studies. RESULTS: A locus tagged by rs10190845 was identified on chromosome 2q13, which was significantly associated with clinical vertebral fracture (P=1.04×10-9) with a large effect size (OR 1.74, 95% CI 1.06 to 2.6). Bioinformatic analysis of this locus identified several potentially functional SNPs that are associated with expression of the positional candidate genes TTL (tubulin tyrosine ligase) and SLC20A1 (solute carrier family 20 member 1). Three other suggestive loci were identified on chromosomes 1p31, 11q12 and 15q11. All these loci were novel and had not previously been associated with bone mineral density or clinical fractures. CONCLUSION: We have identified a novel genetic variant that is associated with clinical vertebral fractures by mechanisms that are independent of BMD. Further studies are now in progress to validate this association and evaluate the underlying mechanism.Funding: ORCADES was supported by the Chief Scientist Office of the Scottish Government (CZB/4/276, CZB/4/710), the Royal Society, the MRC Human Genetics Unit, Arthritis Research UK and the European Union framework programme 6 EUROSPAN project (contract no. LSHG-CT-2006-018947). DNA extractions were performed at the Wellcome Trust Clinical Research Facility in Edinburgh. We would like to acknowledge the invaluable contributions of Lorraine Anderson and the research nurses in Orkney, the administrative team in Edinburgh and the people of Orkney. CABRIO was supported by the Instituto de Salud Carlos III and Fondos FEDER from the EU (PI 11/1092 and PI12/615). The AOGC study was funded by the Australian National Health and Medical Research Council (Project grant 511132). Lothian Birth Cohort 1921 phenotype collection was supported by the UK’s Biotechnology and Biological Sciences Research Council (BBSRC), The Royal Society and The Chief Scientist Office of the Scottish Government. Phenotype collection in the Lothian Birth Cohort 1936 was supported by Age UK (The Disconnected Mind project). Genotyping of the cohorts was funded by the BBSRC. The work was undertaken by the University of Edinburgh Centre for Cognitive Ageing and Cognitive Epidemiology, part of the cross council Lifelong Health and Wellbeing Initiative (MR/K026992/1). Funding from the BBSRC and Medical Research Council (MRC) is gratefully acknowledged. Research work on Slovenian case and control samples was funded by Slovenian Research Agency (project no. P3-0298 and J3-2330). The Danish National Birth Cohort (DNBC) is a result of major grants from the Danish National Research Foundation, the Danish Pharmacists’Fund, the Egmont Foundation, the March of Dimes Birth Defects Foundation, the Augustinus Foundation and the Health Fund of the Danish Health Insurance Societies. The DNBC biobank is a part of the Danish National Biobank resource, which is supported by the Novo Nordisk Foundation. Dr Bjarke Feenstra is supported by an Oak Foundation Fellowship. The Framingham Study was funded by grants from the US National Institute for Arthritis, Musculoskeletal and Skin Diseases and National Institute on Aging (R01 AR 41398 and R01 AR061162; DPK and R01 AR 050066; DK). The Framingham Heart Study of the National Heart, Lung, and Blood Institute of the National Institutes of Health and Boston University School of Medicine were supported by the National Heart, Lung, and Blood Institute’s Framingham Heart Study (N01-HC-25195) and its contract with Affymetrix, Inc. for genotyping services (N02-HL-6-4278). Analyses reflect intellectual input and resource development from the Framingham Heart Study investigators participating in the SNP Health Association Resource (SHARe) project. A portion of this research was conducted using the Linux Cluster for Genetic Analysis (LinGA-II) funded by the Robert Dawson Evans Endowment of the Department of Medicine at Boston University School of Medicine and Boston Medical Center. This research was performed within the Genetic Factors for Osteoporosis (GEFOS) consortium, funded by the European Commission (HEALTH-F2-2008-201865-GEFOS).Acknowledgments: The authors are grateful to the patients and controls from the different centres who agreed to participate in this study. We would like to thank Ms Dilruba Kabir at the Rheumatology and Bone Disease Unit, CGEM-IGMM, Edinburgh, UK; Mr Matt Sims at the MRC Epidemiology Unit, University of Cambridge, UK; Ms Mila Jhamai and Ms Sarah Higgins at the Genetics Laboratory of Erasmus MC, Rotterdam, The Netherlands; Ms Johanna Hadler, Ms Kathryn A Addison and Ms Karena Pryce of the University of Queensland Centre for Clinical Genomics, Brisbane, Australia, for technical support on the genotyping stage; and Mr Marijn Verkerk and Dr Anis Abuseiris at the Genetics Laboratory of Erasmus MC, Rotterdam, for assistance on the data analysis. We would like to acknowledge the invaluable contributions of Lorraine Anderson and the research nurses in Orkney, the administrative team in Edinburgh and the people of Orkney. We would also like to thank Professor Nick Gilbert and Dr Giovanny Rodriguez-Blanco for their comments and advice on the manuscript preparation. This study makes use of data generated by the Wellcome Trust Case Control Consortium. A full list of the investigators who contributed to the generation of the data is available at www.wtccc.org.uk

Skeletal Site-Related Variation in Human Trabecular Bone Transcriptome and Signaling

BACKGROUND: The skeletal site-specific influence of multiple genes on bone morphology is recognised, but the question as to how these influences may be exerted at the molecular and cellular level has not been explored. METHODOLOGY: To address this question, we have compared global gene expression profiles of human trabecular bone from two different skeletal sites that experience vastly different degrees of mechanical loading, namely biopsies from iliac crest and lumbar spinal lamina. PRINCIPAL FINDINGS: In the lumbar spine, compared to the iliac crest, the majority of the differentially expressed genes showed significantly increased levels of expression; 3406 transcripts were up- whilst 838 were down-regulated. Interestingly, all gene transcripts that have been recently demonstrated to be markers of osteocyte, as well as osteoblast and osteoclast-related genes, were markedly up-regulated in the spine. The transcriptome data is consistent with osteocyte numbers being almost identical at the two anatomical sites, but suggesting a relatively low osteocyte functional activity in the iliac crest. Similarly, osteoblast and osteoclast expression data suggested similar numbers of the cells, but presented with higher activity in the spine than iliac crest. This analysis has also led to the identification of expression of a number of transcripts, previously known and novel, which to our knowledge have never earlier been associated with bone growth and remodelling. CONCLUSIONS AND SIGNIFICANCE: This study provides molecular evidence explaining anatomical and micro-architectural site-related changes in bone cell function, which is predominantly attributable to alteration in cell transcriptional activity. A number of novel signaling molecules in critical pathways, which have been hitherto not known to be expressed in bone cells of mature vertebrates, were identified