Loss of the von-Hippel Lindau gene in osteocytes alters B lymphocyte development and immune response

The von-Hippel Lindau protein (VHL) regulates hypoxia-inducible factor (HIF) degradation. Conditional-deletion of VHL in osteoblasts and in hematopoietic progenitors have demonstrated a role for VHL in these cell types, but how changes in bone homeostasis affects immune development are not fully understood. Osteoblasts (OBs) support B cell development. OBs mature into osteocytes (OCYs) and to test if OCYs also support B cell development, we generated OCY-specific (Dmp-Cre) conditional VHL knockout (cKO) mice. cKO mice display dysregulated bone growth, high bone density, smaller bone marrow (BM) cavity volume, and overall decrease in BM cellularity compared to controls. In line with this, the frequencies and numbers of B cell precursors (B220+ IgM−), immature B (B220+ IgMint) and recirculating B (B220+ IgMhigh) cells in the bone marrow are significantly decreased. These data suggest that the change in bone homeostasis adversely affects B cell development in a cell-extrinsic manner. Although the absolute numbers of B cells in the spleens of cKOs were normal, we hypothesized that loss of VHL in OCYs could permanently and adversely affect B cell functional response. To test this, VHL cKO mice were immunized with a T-independent antigen (NP-Ficoll). IgM and IgG3 antigen-specific titers to NP-Ficoll in VHL cKOs were similar to WT mice. Testing of B cell responses to T-dependent antigens (NP-OVA) is ongoing. Taken together, our data suggest that OCY-specific VHL deletion might be helpful to increase bone density, but at the expense of B cell numbers. Further investigation is required to understand the molecular mechanisms by which VHL in OCYs contribute to its putative role as a B cell niche cell.

Shared and Unique Features Distinguishing Follicular T Helper and Regulatory Cells of Peripheral Lymph Node and Peyer's Patches

Follicular helper (TFH) and regulatory (TFR) cells are critical players in managing germinal center (GC) reactions that accomplish effective humoral immune responses. Transcriptome analyses were done comparing gene regulation of TFH and TFR cells isolated from Peyer’s Patches (PP) and immunized peripheral lymph nodes (pLNs) revealing many regulatory patterns common to all follicular cells. However, in contrast to TFH cells, the upregulation or downregulation of many genes was attenuated substantially in pLN TFR cells when compared to those of PP. Additionally, PP but not pLN TFR cells were largely unresponsive to IL2 and expressed Il4 as well as Il21. Together with fundamental differences in gene expression that were found between cells of both compartments this emphasizes specific adaptations of follicular T cell functions to their micro-milieu. Moreover, although GL7 expression distinguishes matured follicular T cells, GL7⁺ as well as GL7⁻ cells are present in the GC.

Genetics of Sost/SOST in sclerosteosis and van Buchem disease animal models

Sclerosteosis and van Buchem disease (VBD) are two rare autosomal recessive disorders that results from osteoblast hyperactivity, in which progressive bone overgrowth leads to very dense bones, distortion of the face, and entrapment of cranial nerves. Sclerosteosis is caused by loss-of-function mutations in the SOST gene which encodes a secreted glycoprotein, sclerostin. VBD is caused by a noncoding deletion that removes a SOST-specific regulatory element in bone. In bone, SOST is expressed predominantly by osteocytes and sclerostin suppresses bone formation by inhibiting the canonical Wnt signaling pathway. Here we describe how human genetics studies in sclerosteosis and VBD patients, in combination with the generation of transgenic and knockout mice, has led to a better understanding of the role of sclerostin in bone metabolism.

Wnt co-receptors Lrp5 and Lrp6 differentially mediate Wnt3a signaling in osteoblasts

Wnt3a is a major regulator of bone metabolism however, very few of its target genes are known in bone. Wnt3a preferentially signals through transmembrane receptors Frizzled and co-receptors Lrp5/6 to activate the canonical signaling pathway. Previous studies have shown that the canonical Wnt co-receptors Lrp5 and Lrp6 also play an essential role in normal postnatal bone homeostasis, yet, very little is known about specific contributions by these co-receptors in Wnt3a-dependent signaling. We used high-throughput sequencing technology to identify target genes regulated by Wnt3a in osteoblasts and to elucidate the role of Lrp5 and Lrp6 in mediating Wnt3a signaling. Our study identified 782 genes regulated by Wnt3a in primary calvarial osteoblasts. Wnt3a up-regulated the expression of several key regulators of osteoblast proliferation/ early stages of differentiation while inhibiting genes expressed in later stages of osteoblastogenesis. We also found that Lrp6 is the key mediator of Wnt3a signaling in osteoblasts and Lrp5 played a less significant role in mediating Wnt3a signaling.

Maternal exposure to an environmentally relevant dose of triclocarban results in perinatal exposure and potential alterations in offspring development in the mouse model

Triclocarban (TCC) is among the top 10 most commonly detected wastewater contaminants in both concentration and frequency. Its presence in water, as well as its propensity to bioaccumulate, has raised numerous questions about potential endocrine and developmental effects. Here, we investigated whether exposure to an environmentally relevant concentration of TCC could result in transfer from mother to offspring in CD-1 mice during gestation and lactation using accelerator mass spectrometry (AMS). ¹⁴C-TCC (100 nM) was administered to dams through drinking water up to gestation day 18, or from birth to post-natal day 10. AMS was used to quantify ¹⁴C-concentrations in offspring and dams after exposure. We demonstrated that TCC does effectively transfer from mother to offspring, both trans-placentally and via lactation. TCC-related compounds were detected in the tissues of offspring with significantly higher concentrations in the brain, heart and fat. In addition to transfer from mother to offspring, exposed offspring were heavier in weight than unexposed controls demonstrating an 11% and 8.5% increase in body weight for females and males, respectively. Quantitative real-time polymerase chain reaction (qPCR) was used to examine changes in gene expression in liver and adipose tissue in exposed offspring. qPCR suggested alterations in genes involved in lipid metabolism in exposed female offspring, which was consistent with the observed increased fat pad weights and hepatic triglycerides. This study represents the first report to quantify the transfer of an environmentally relevant concentration of TCC from mother to offspring in the mouse model and evaluate bio-distribution after exposure using AMS. Our findings suggest that early-life exposure to TCC may interfere with lipid metabolism and could have implications for human health.

Abstract 1977: Tracking cancer colonization in xenografts using ultrasensitive accelerator mass spectrometry methods

Background: The inability to effectively treat metastases is the main reason for the limited progress in reducing the rates of cancer morbidity and mortality. One major drawback is the lack of quantitative assays for assessing the size and tissue prevalence of tumors in newly diagnosed individuals. Current methods for quantifying tumor burden are mainly qualitative and include measuring the gross weight of the affected organ, counting tumors on the surface of the organ, or evaluating a small sample of the organ using histologic sections. These methods are crude measures of tumor burden and size distribution, and in the case of histology, they are time consuming, difficult to process an adequate sample size and non-quantitative.

Methods: Animal models of metastasis have been useful in identifying genes that regulate susceptibility to the development and progression of metastasis and have helped to highlight potential novel targets for drug development. In particular several small animal imaging technologies including magnetic resonance imaging, high frequency ultrasound, and optical imaging have been recently applied to this task. Each of these methods may be useful for specific research projects, based on their unique combination of resolution, image acquisition time, animal throughput, and cost-effectiveness, yet none of these modalities adequately address the need for rapid quantification of tumors across the entire organism, nor do they assess therapeutic effectiveness in eradicating cancer in xenograft models. We have developed an Accelerator Mass Spectrometry (AMS)-based high precision quantitative method for assessing the metastatic potential of primary tumors isolated from newly diagnosed patients.

Results: Our AMS-based methodology to study metastasis uses xenograft cancer cells labeled with 14C-labeled thymidine that are delivered intravenously into NSG mice and allowed to develop metastatic cancer over the course of up to 10 weeks. At the end of the experiment, all vital organs are collected; the DNA is isolated and is examined by AMS for the presence of 14C-signal. The labeling was optimized to achieve sufficient signal such that a tumor derived from a single cell could be detected by AMS, in secondary tumors, in vivo, independent of histological data.

Conclusions: Using this approach we have determined that tissue colonization by tumor cells is a very rare event, where most metastatic tumors are initiated by less than 10 cells delivered into NSG mice. Further optimization of these techniques will allow us to explore the metastatic potential of primary tumors, isolated from biopsies and expanded in Avatar mice.

This study was supported in part by NIH P41MI03483 and was conducted under the auspices of the USDOE by LLNL (DE-AC52-07NA27344). IM number: LLNL-678306

Citation Format: Nicholas R. Hum, Kelly A. Martin, Michael Malfatti, Kurt Haack, Bruce A. Buchholz, Gabriela G. Loots. Tracking cancer colonization in xenografts using ultrasensitive accelerator mass spectrometry methods [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1977. doi:10.1158/1538-7445.AM2017-1977

Global molecular changes in a tibial compression induced ACL rupture model of post-traumatic osteoarthritis

Joint injury causes post-traumatic osteoarthritis (PTOA). About ∼50% of patients rupturing their anterior cruciate ligament (ACL) will develop PTOA within 1-2 decades of the injury, yet the mechanisms responsible for the development of PTOA after joint injury are not well understood. In this study, we examined whole joint gene expression by RNA sequencing (RNAseq) at 1 day, 1-, 6-, and 12 weeks post injury, in a non-invasive tibial compression (TC) overload mouse model of PTOA that mimics ACL rupture in humans. We identified 1446 genes differentially regulated between injured and contralateral joints. This includes known regulators of osteoarthritis such as MMP3, FN1, and COMP, and several new genes including Suco, Sorcs2, and Medag. We also identified 18 long noncoding RNAs that are differentially expressed in the injured joints. By comparing our data to gene expression data generated using the surgical destabilization of the medial meniscus (DMM) PTOA model, we identified several common genes and shared mechanisms. Our study highlights several differences between these two models and suggests that the TC model may be a more rapidly progressing model of PTOA. This study provides the first account of gene expression changes associated with PTOA development and progression in a TC model. © 2016 The Authors. Journal of Orthopaedic Research Published by Wiley Periodicals, Inc. J Orthop Res 35:474-485, 2017.

Transcriptional control of Sost in bone

Sclerostin is an osteocyte derived negative regulator of bone formation. A highly specific expression pattern and the exclusive bone phenotype have made Sclerostin an attractive target for therapeutic intervention in treating metabolic bone diseases such as osteoporosis and in facilitating fracture repair. Understanding the molecular mechanisms that regulate Sclerostin transcription is of great interest as it may unveil new avenues for therapeutic approaches. Such studies may also elucidate how various signaling pathways intersect to modulate bone metabolism. Here we review the current understanding of the upstream molecular mechanisms that regulate Sost/SOST transcription, in bone.

Wnt Signaling Pathway Inhibitor Sclerostin Inhibits Angiotensin II-Induced Aortic Aneurysm and Atherosclerosis

OBJECTIVE

Sclerostin (SOST) has been identified as an important regulator of bone formation; however, it has not been previously implicated in arterial disease. The aim of this study was to assess the role of SOST in aortic aneurysm (AA) and atherosclerosis using human samples, a mouse model, and in vitro investigations.

APPROACH AND RESULTS

SOST protein was downregulated in human and mouse AA samples compared with controls. Transgenic introduction of human SOST in apolipoprotein E-deficient (ApoE^-/-) mice (SOST^Tg .ApoE^-/-) and administration of recombinant mouse Sost inhibited angiotensin II-induced AA and atherosclerosis. Serum concentrations of several proinflammatory cytokines were significantly reduced in SOST^Tg .ApoE^-/- mice. Compared with controls, the aortas of mice receiving recombinant mouse Sost and SOST^Tg .ApoE^-/- mice showed reduced matrix degradation, reduced elastin breaks, and preserved collagen. Decreased inflammatory cell infiltration and a reduction in the expression of wingless-type mouse mammary virus integration site/β-catenin responsive genes, including matrix metalloproteinase-9, osteoprotegerin, and osteopontin, were observed in the aortas of SOST^Tg .ApoE^-/- mice. SOST expression was downregulated and the wingless-type mouse mammary virus integration site/β-catenin pathway was activated in human AA samples. The cytosine-phosphate-guanine islands in the SOST gene promoter showed significantly higher methylation in human AA samples compared with controls. Incubation of vascular smooth muscle cells with the demethylating agent 5-azacytidine resulted in upregulation of SOST, suggesting that SOST is epigenetically regulated.

CONCLUSIONS

This study identifies that SOST is expressed in the aorta and downregulated in human AA possibly because of epigenetic silencing. Upregulating SOST inhibits AA and atherosclerosis development, with potential important implications for treating these vascular diseases.

Sost, independent of the non-coding enhancer ECR5, is required for bone mechanoadaptation

Sclerostin (Sost) is a negative regulator of bone formation that acts upon the Wnt signaling pathway. Sost is mechanically regulated at both mRNA and protein level such that loading represses and unloading enhances Sost expression, in osteocytes and in circulation. The non-coding evolutionarily conserved enhancer ECR5 has been previously reported as a transcriptional regulatory element required for modulating Sost expression in osteocytes. Here we explored the mechanisms by which ECR5, or several other putative transcriptional enhancers regulate Sost expression, in response to mechanical stimulation. We found that in vivo ulna loading is equally osteoanabolic in wildtype and Sost^-/- mice, although Sost is required for proper distribution of load-induced bone formation to regions of high strain. Using Luciferase reporters carrying the ECR5 non-coding enhancer and heterologous or homologous hSOST promoters, we found that ECR5 is mechanosensitive in vitro and that ECR5-driven Luciferase activity decreases in osteoblasts exposed to oscillatory fluid flow. Yet, ECR5^-/- mice showed similar magnitude of load-induced bone formation and similar periosteal distribution of bone formation to high-strain regions compared to wildtype mice. Further, we found that in contrast to Sost^-/- mice, which are resistant to disuse-induced bone loss, ECR5^-/- mice lose bone upon unloading to a degree similar to wildtype control mice. ECR5 deletion did not abrogate positive effects of unloading on Sost, suggesting that additional transcriptional regulators and regulatory elements contribute to load-induced regulation of Sost.

Abstract 3891: Development of quantitative methods for assessing metastatic potential of human primary tumors

Background: The inability to effectively treat metastases is the main reason for the limited progress in reducing the rates of cancer morbidity and mortality. One major drawback is the lack of quantitative assays for assessing the size and tissue prevalence of tumors in newly diagnosed individuals. Current methods for quantifying tumor burden are mainly qualitative and include measuring the gross weight of the affected organ, counting tumors on the surface of the organ, or evaluating a small sample of the organ using histologic sections. These methods are crude measures of tumor burden and size distribution, and in the case of histology, they are time consuming, difficult to process an adequate sample size and non-quantitative.

Methods: Animal models of metastasis have been useful in identifying genes that regulate susceptibility to the development and progression of metastasis and helped highlight potential novel targets for drug development. In particular several small animal imaging technologies including magnetic resonance imaging, high frequency ultrasound, and optical imaging have been recently applied to this task. Each of these methods may be useful for specific research projects, based on their unique combination of resolution, image acquisition time, animal throughput, and cost-effectiveness, yet none of these modalities adequately address the need for rapid quantification of tumors across the entire organism, nor do they assess therapeutic effectiveness in eradicating cancer in xenograft models. We have developed an Accelerator Mass Spectrometry (AMS)-based high precision quantitative method for assessing the metastatic potential of primary tumors isolated from newly diagnosed patients.

Results: Our Accelerator Mass Spectrometry-based methodology to study metastasis uses xenograft cancer cells labeled with 14C-labeled thymidine that are delivered intravenously into NSG mice and allowed to develop metastatic cancer over the course of up to 10 weeks. At the end of the experiment, all vital organs are collected; the DNA is isolated and is examined by AMS for the presence of 14C-signal. The labeling was optimized to achieve sufficient signal such that a tumor derived from a single cell could be detected by AMS, in secondary tumors, in vivo, independent of histological data.

Conclusions: Using this novel approach we have evaluated the metastatic potential of several prostate cancer cell lines, characterized stem-cell like sublines derived from prostate cancer cell lines [PC3] and examined tissue tropism of cancer sublines derived from kidney and liver metastatic tumors. Further optimization of these techniques will allow us to explore the metastatic potential of primary tumors, isolated from biopsies and expanded in Avatar mice.

Citation Format: Kelly A. Martin, Nicholas Hum, Kurt W. Haack, Bruce A. Buchholz, Gabriela G. Loots. Development of quantitative methods for assessing metastatic potential of human primary tumors. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3891.

Abstract 2106: Transcriptome analysis of patient-derived bladder cancer xenografts identifies genes associated with chemoresistance

Background: Bladder cancer is among the ten most common cancers, with about ∼380,000 new cases and ∼150,000 deaths per year worldwide. Platinum-based combination chemotherapy is commonly used to treat advanced bladder cancer. It has been shown that only ∼50% of the patients with advanced bladder cancer respond to platinum-based therapy.

Methods: We employed a patient-derived bladder cancer xenograft (PDX) platform to characterize the molecular mechanisms contributing to resistance to gemcitabine-cisplatin combination therapy in advanced bladder cancer and to identify novel candidates that can be targeted to treat chemotherapy resistant bladder cancer. Transcriptome profiling of P0 (passage 0) bladder cancer xenograft tumors from 4 PDX lines (2 gemcitabine-cisplatin resistant lines and 2 drug sensitive lines) was performed by RNA-Seq analysis.

Results: PDXs retained the morphology fidelity and shared 92-97% of genetic alterations of parental cancer cells. The RNA-seq data suggested the presence of significant differences between the transcription profiles of drug-sensitive and drug-resistant tumors. We identified 333 genes >2 fold up or down regulated in the drug resistant tumors compared to the drug sensitive tumors. Genes down-regulated in drug resistant tumors include tight junction protein CLDN3 and regulators of G-protein signaling RGS2 and RGS3. Significantly up-regulated genes include metabolic enzymes ALDH2, ALDH3A1, ALDH4A1 and ALDH7A1, transporter proteins ABCA1, SLC1A4, SLC2A5, SLC30A1, SLC39A6, SLC7A5 and SLC9A3, Notch ligand JAG2, Growth hormone receptor GHR and transmembrane glycoprotein GPNMB. Consistent with the change of cell surface proteins such as GHR and GPNMB, the MAPK and the PI3K-AKT pathways were upregulated when PDXs became resistant to cisplatin treatment.

Conclusion: Chemoresistance to gemcitabine and cisplatin is associated with altered expression of several cell surface proteins and upregulation of the downstream signaling pathways. Targeting these cell surface proteins can possibly be harnessed to overcome chemoresistance. GPNMB, a type I transmembrane protein highly up-regulated in the drug resistant tumors, has previously been shown to be over-expressed in various cancers. Targeting GPNMB with an antibody-drug-conjugate, glembatumumab vedotin, has shown promising results in treating several cancers including breast cancer and osteosarcoma. Further studies will elucidate whether targeting GPNMB is an effective strategy for the treatment of chemotherapy resistant bladder cancer.

Citation Format: Aimy Sebastian, Kelly A. Martin, Chong-xian Pan, Ai-hong Ma, Ralph W. deVere White, Gabriela G. Loots. Transcriptome analysis of patient-derived bladder cancer xenografts identifies genes associated with chemoresistance. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2106.

Sostdc1 deficiency accelerates fracture healing by promoting the expansion of periosteal mesenchymal stem cells

Loss of Sostdc1, a growth factor paralogous to Sost, causes the formation of ectopic incisors, fused molars, abnormal hair follicles, and resistance to kidney disease. Sostdc1 is expressed in the periosteum, a source of osteoblasts, fibroblasts and mesenchymal progenitor cells, which are critically important for fracture repair. Here, we investigated the role of Sostdc1 in bone metabolism and fracture repair. Mice lacking Sostdc1 (Sostdc1(-/-)) had a low bone mass phenotype associated with loss of trabecular bone in both lumbar vertebrae and in the appendicular skeleton. In contrast, Sostdc1(-/-) cortical bone measurements revealed larger bones with higher BMD, suggesting that Sostdc1 exerts differential effects on cortical and trabecular bone. Mid-diaphyseal femoral fractures induced in Sostdc1(-/-) mice showed that the periosteal population normally positive for Sostdc1 rapidly expands during periosteal thickening and these cells migrate into the fracture callus at 3days post fracture. Quantitative analysis of mesenchymal stem cell (MSC) and osteoblast populations determined that MSCs express Sostdc1, and that Sostdc1(-/-) 5day calluses harbor >2-fold more MSCs than fractured wildtype controls. Histologically a fraction of Sostdc1-positive cells also expressed nestin and α-smooth muscle actin, suggesting that Sostdc1 marks a population of osteochondral progenitor cells that actively participate in callus formation and bone repair. Elevated numbers of MSCs in D5 calluses resulted in a larger, more vascularized cartilage callus at day 7, and a more rapid turnover of cartilage with significantly more remodeled bone and a thicker cortical shell at 21days post fracture. These data support accelerated or enhanced bone formation/remodeling of the callus in Sostdc1(-/-) mice, suggesting that Sostdc1 may promote and maintain mesenchymal stem cell quiescence in the periosteum.

Sclerostin antibody treatment improves fracture outcomes in a Type I diabetic mouse model

Type 1 diabetes mellitus (T1DM) patients have osteopenia and impaired fracture healing due to decreased osteoblast activity. Further, no adequate treatments are currently available that can restore impaired healing in T1DM; hence a significant need exists to investigate new therapeutics for treatment of orthopedic complications. Sclerostin (SOST), a WNT antagonist, negatively regulates bone formation, and SostAb is a potent bone anabolic agent. To determine whether SOST antibody (SostAb) treatment improves fracture healing in streptozotocin (STZ) induced T1DM mice, we administered SostAb twice weekly for up to 21days post-fracture, and examined bone quality and callus outcomes at 21days and 42days post-fracture (11 and 14weeks of age, respectively). Here we show that SostAb treatment improves bone parameters; these improvements persist after cessation of antibody treatment. Markers of osteoblast differentiation such as Runx2, collagen I, osteocalcin, and DMP1 were reduced, while an abundant number of SP7/osterix-positive early osteoblasts were observed on the bone surface of STZ calluses. These results suggest that STZ calluses have poor osteogenesis resulting from failure of osteoblasts to fully differentiate and produce mineralized matrix, which produces a less mineralized callus. SostAb treatment enhanced fracture healing in both normal and STZ groups, and in STZ+SostAb mice, also reversed the lower mineralization seen in STZ calluses. Micro-CT analysis of calluses revealed improved bone parameters with SostAb treatment, and the mineralized bone was comparable to Controls. Additionally, we found sclerostin levels to be elevated in STZ mice and β-catenin activity to be reduced. Consistent with its function as a WNT antagonist, SostAb treatment enhanced β-catenin activity, but also increased the levels of SOST in the callus and in circulation. Our results indicate that SostAb treatment rescues the impaired osteogenesis seen in the STZ induced T1DM fracture model by facilitating osteoblast differentiation and mineralization of bone.

Parathyroid hormone regulation of hypoxia-inducible factor signaling in osteoblastic cells

Osteoblasts perceive and respond to changes in their pericellular environment, including biophysical signals and oxygen availability, to elicit an anabolic or catabolic response. Parathyroid hormone (PTH) affects each arm of skeletal remodeling, with net anabolic or catabolic effects dependent upon duration of exposure. Similarly, the capacity of osteoblastic cells to perceive pericellular oxygen has a profound effect on skeletal mass and architecture, as mice expressing stable hypoxia-inducible factor (HIF)-1α and -2α demonstrate age-dependent increases in bone volume per tissue volume and osteoblast number. Further, HIF levels and signaling can be influenced in an oxygen-independent manner. Because the cellular mechanisms involved in PTH regulation of the skeleton remain vague, we sought whether PTH could influence HIF-1α expression and HIF-α-driven luciferase activity independently of altered oxygen availability. Using UMR106.01 mature osteoblasts, we observed that 100nM hPTH(1-34) decreased HIF-1α and HIF-responsive luciferase activity in a process involving heat shock protein 90 (Hsp90) and cyclic AMP but not intracellular calcium. Altering activity of the small GTPase RhoA and its effector kinase ROCK altered HIF-α-driven luciferase activity in the absence and presence of PTH. Taken together, these data introduce PTH as a regulator of oxygen-independent HIF-1α levels through a mechanism involving cyclic AMP, Hsp90, and the cytoskeleton.

1α,25-dihydroxyvitamin D3 stimulates human SOST gene expression and sclerostin secretion

Sclerostin, the SOST gene product, is a negative regulator of bone formation and a positive regulator of bone resorption. In this study, treatment of human primary osteoblasts, including cells differentiated to an osteocyte-like stage, with 1α,25-dihydroxyvitaminD3 (1,25D) resulted in the dose-dependent increased expression of SOST mRNA. A similar effect was observed in human trabecular bone samples cultured ex vivo, and in osteocyte-like cultures of differentiated SAOS2 cells. Treatment of SAOS2 cells with 1,25D resulted in the production and secretion of sclerostin protein. In silico analysis of the human SOST gene revealed a single putative DR3-type vitamin D response element (VDRE) at position -6216 bp upstream of the transcription start site (TSS). This sequence was confirmed to have strong VDRE activity by luciferase reporter assays and electrophoretic mobility shift analysis (EMSA). Sequence substitution in the VDR/RXR half-sites abolished VDRE reporter activity and binding of nuclear proteins. A 6.3 kb fragment of the human proximal SOST promoter demonstrated responsiveness to 1,25D. The addition of the evolutionary conserved region 5 (ECR5), a known bone specific enhancer region, ahead of the 6.3 kb fragment increased basal promoter activity but did not increase 1,25D responsiveness. Site-specific mutagenesis abolished the responsiveness of the 6.3 kb promoter to 1,25D. We conclude that 1,25D is a direct regulator of human SOST gene and sclerostin protein expression, extending the pathways of control of sclerostin expression. At least some of this responsiveness is mediated by the identified classical VDRE however the nature of the transcriptional regulation by 1,25D warrants further investigation.

Nanocomposite scaffold for chondrocyte growth and cartilage tissue engineering: effects of carbon nanotube surface functionalization

The goal of this study was to assess the long-term biocompatibility of single-wall carbon nanotubes (SWNTs) for tissue engineering of articular cartilage. We hypothesized that SWNT nanocomposite scaffolds in cartilage tissue engineering can provide an improved molecular-sized substrate for stimulation of chondrocyte growth, as well as structural reinforcement of the scaffold's mechanical properties. The effect of SWNT surface functionalization (-COOH or -PEG) on chondrocyte viability and biochemical matrix deposition was examined in two-dimensional cultures, in three-dimensional (3D) pellet cultures, and in a 3D nanocomposite scaffold consisting of hydrogels+SWNTs. Outcome measures included cell viability, histological and SEM evaluation, GAG biochemical content, compressive and tensile biomechanical properties, and gene expression quantification, including extracellular matrix (ECM) markers aggrecan (Agc), collagen-1 (Col1a1), collagen-2 (Col2a1), collagen-10 (Col10a1), surface adhesion proteins fibronectin (Fn), CD44 antigen (CD44), and tumor marker (Tp53). Our findings indicate that chondrocytes tolerate functionalized SWNTs well, with minimal toxicity of cells in 3D culture systems (pellet and nanocomposite constructs). Both SWNT-PEG and SWNT-COOH groups increased the GAG content in nanocomposites relative to control. The compressive biomechanical properties of cell-laden SWNT-COOH nanocomposites were significantly elevated relative to control. Increases in the tensile modulus and ultimate stress were observed, indicative of a tensile reinforcement of the nanocomposite scaffolds. Surface coating of SWNTs with -COOH also resulted in increased Col2a1 and Fn gene expression throughout the culture in nanocomposite constructs, indicative of increased chondrocyte metabolic activity. In contrast, surface coating of SWNTs with a neutral -PEG moiety had no significant effect on Col2a1 or Fn gene expression, suggesting that the charged nature of the -COOH surface functionalization may promote ECM expression in this culture system. The results of this study indicate that SWNTs exhibit a unique potential for cartilage tissue engineering, where functionalization with bioactive molecules may provide an improved substrate for stimulation of cellular growth and repair.

Evaluation of nanolipoprotein particles (NLPs) as an in vivo delivery platform

Nanoparticles hold great promise for the delivery of therapeutics, yet limitations remain with regards to the use of these nanosystems for efficient long-lasting targeted delivery of therapeutics, including imparting functionality to the platform, in vivo stability, drug entrapment efficiency and toxicity. To begin to address these limitations, we evaluated the functionality, stability, cytotoxicity, toxicity, immunogenicity and in vivo biodistribution of nanolipoprotein particles (NLPs), which are mimetics of naturally occurring high-density lipoproteins (HDLs). We found that a wide range of molecules could be reliably conjugated to the NLP, including proteins, single-stranded DNA, and small molecules. The NLP was also found to be relatively stable in complex biological fluids and displayed no cytotoxicity in vitro at doses as high as 320 µg/ml. In addition, we observed that in vivo administration of the NLP daily for 14 consecutive days did not induce significant weight loss or result in lesions on excised organs. Furthermore, the NLPs did not display overt immunogenicity with respect to antibody generation. Finally, the biodistribution of the NLP in vivo was found to be highly dependent on the route of administration, where intranasal administration resulted in prolonged retention in the lung tissue. Although only a select number of NLP compositions were evaluated, the findings of this study suggest that the NLP platform holds promise for use as both a targeted and non-targeted in vivo delivery vehicle for a range of therapeutics.

The use of nanolipoprotein particles to enhance the immunostimulatory properties of innate immune agonists against lethal influenza challenge

Recent studies have demonstrated that therapies targeting the innate immune system have the potential to provide transient, non-specific protection from a variety of infectious organisms; however, the potential of enhancing the efficacy of such treatments using nano-scale delivery platforms requires more in depth evaluation. As such, we employed a nanolipoprotein (NLP) platform to enhance the efficacy of innate immune agonists. Here, we demonstrate that the synthetic Toll-like receptor (TLR) agonists monophosphoryl lipid A (MPLA) and CpG oligodeoxynucleotides (CpG) can be readily incorporated into NLPs. Conjugation of MPLA and CpG to NLPs (MPLA:NLP and CpG:NLP, respectively) significantly enhanced their immunostimulatory profiles both in vitro and in vivo compared to administration of agonists alone, as evidenced by significant increases in cytokine production, cell surface expression of activation markers, and upregulation of immunoregulatory genes. Importantly, enhancement of cytokine production by agonist conjugation to NLPs was also observed in primary human dendritic cells. Furthermore, BALB/c mice pretreated with CpG:NLP constructs survived a lethal influenza challenge whereas pretreatment with CpG alone had no effect on survival.

Prostate cancer invasion and metastasis: insights from mining genomic data

Prostate cancer (PCa) is the second most commonly diagnosed malignancy in men in the Western world and the second leading cause of cancer-related deaths among men worldwide. Although most cancers have the potential to metastasize under appropriate conditions, PCa favors the skeleton as a primary site of metastasis, suggesting that the bone microenvironment is conducive to its growth. PCa metastasis proceeds through a complex series of molecular events that include angiogenesis at the site of the original tumor, local migration within the primary site, intravasation into the blood stream, survival within the circulation, extravasation of the tumor cells to the target organ and colonization of those cells within the new site. In turn, each one of these steps involves a complicated chain of events that utilize multiple protein-protein interactions, protein signaling cascades and transcriptional changes. Despite the urgent need to improve current biomarkers for diagnosis, prognosis and drug resistance, advances have been slow. Global gene expression methods such as gene microarrays and RNA sequencing enable the study of thousands of genes simultaneously and allow scientists to examine molecular pathways of cancer pathogenesis. In this review, we summarize the current literature that explored high-throughput transcriptome analysis toward the advancement of biomarker discovery for PCa. Novel biomarkers are strongly needed to enable more accurate detection of PCa, improve prediction of tumor aggressiveness and facilitate the discovery of new therapeutic targets for tailored medicine. Promising molecular markers identified from gene expression profiling studies include HPN, CLU1, WT1, WNT5A, AURKA and SPARC.

Interrogating transcriptional regulatory sequences in Tol2-mediated Xenopus transgenics

Identifying gene regulatory elements and their target genes in vertebrates remains a significant challenge. It is now recognized that transcriptional regulatory sequences are critical in orchestrating dynamic controls of tissue-specific gene expression during vertebrate development and in adult tissues, and that these elements can be positioned at great distances in relation to the promoters of the genes they control. While significant progress has been made in mapping DNA binding regions by combining chromatin immunoprecipitation and next generation sequencing, functional validation remains a limiting step in improving our ability to correlate in silico predictions with biological function. We recently developed a computational method that synergistically combines genome-wide gene-expression profiling, vertebrate genome comparisons, and transcription factor binding-site analysis to predict tissue-specific enhancers in the human genome. We applied this method to 270 genes highly expressed in skeletal muscle and predicted 190 putative cis-regulatory modules. Furthermore, we optimized Tol2 transgenic constructs in Xenopus laevis to interrogate 20 of these elements for their ability to function as skeletal muscle-specific transcriptional enhancers during embryonic development. We found 45% of these elements expressed only in the fast muscle fibers that are oriented in highly organized chevrons in the Xenopus laevis tadpole. Transcription factor binding site analysis identified >2 Mef2/MyoD sites within ∼200 bp regions in 6 of the validated enhancers, and systematic mutagenesis of these sites revealed that they are critical for the enhancer function. The data described herein introduces a new reporter system suitable for interrogating tissue-specific cis-regulatory elements which allows monitoring of enhancer activity in real time, throughout early stages of embryonic development, in Xenopus.

Effect of age and cytoskeletal elements on the indentation-dependent mechanical properties of chondrocytes

Articular cartilage chondrocytes are responsible for the synthesis, maintenance, and turnover of the extracellular matrix, metabolic processes that contribute to the mechanical properties of these cells. Here, we systematically evaluated the effect of age and cytoskeletal disruptors on the mechanical properties of chondrocytes as a function of deformation. We quantified the indentation-dependent mechanical properties of chondrocytes isolated from neonatal (1-day), adult (5-year) and geriatric (12-year) bovine knees using atomic force microscopy (AFM). We also measured the contribution of the actin and intermediate filaments to the indentation-dependent mechanical properties of chondrocytes. By integrating AFM with confocal fluorescent microscopy, we monitored cytoskeletal and biomechanical deformation in transgenic cells (GFP-vimentin and mCherry-actin) under compression. We found that the elastic modulus of chondrocytes in all age groups decreased with increased indentation (15-2000 nm). The elastic modulus of adult chondrocytes was significantly greater than neonatal cells at indentations greater than 500 nm. Viscoelastic moduli (instantaneous and equilibrium) were comparable in all age groups examined; however, the intrinsic viscosity was lower in geriatric chondrocytes than neonatal. Disrupting the actin or the intermediate filament structures altered the mechanical properties of chondrocytes by decreasing the elastic modulus and viscoelastic properties, resulting in a dramatic loss of indentation-dependent response with treatment. Actin and vimentin cytoskeletal structures were monitored using confocal fluorescent microscopy in transgenic cells treated with disruptors, and both treatments had a profound disruptive effect on the actin filaments. Here we show that disrupting the structure of intermediate filaments indirectly altered the configuration of the actin cytoskeleton. These findings underscore the importance of the cytoskeletal elements in the overall mechanical response of chondrocytes, indicating that intermediate filament integrity is key to the non-linear elastic properties of chondrocytes. This study improves our understanding of the mechanical properties of articular cartilage at the single cell level.

Long-term administration of AMD3100, an antagonist of SDF-1/CXCR4 signaling, alters fracture repair

Fracture healing involves rapid stem and progenitor cell migration, homing, and differentiation. SDF-1 (CXCL12) is considered a master regulator of CXCR4-positive stem and progenitor cell trafficking to sites of ischemic (hypoxic) injury and regulates their subsequent differentiation into mature reparative cells. In this study, we investigated the role of SDF-1/CXCR4 signaling in fracture healing where vascular disruption results in hypoxia and SDF-1 expression. Mice were injected with AMD3100, a CXCR4 antagonist, or vehicle twice daily until euthanasia with the intent to impair stem cell homing to the fracture site and/or their differentiation. Fracture healing was evaluated using micro-computed tomography, histology, quantitative PCR, and mechanical testing. AMD3100 administration resulted in a significantly reduced hyaline cartilage volume (day 14), callus volume (day 42) and mineralized bone volume (day 42) and reduced expression of genes associated with endochondral ossification including collagen Type 1 alpha 1, collagen Type 2 alpha 1, vascular endothelial growth factor, Annexin A5, nitric oxide synthase 2, and mechanistic target of rapamycin. Our data suggest that the SDF-1/CXCR4 signaling plays a central role in bone healing possibly by regulating the recruitment and/or differentiation of stem and progenitor cells.

Absence of sclerostin adversely affects B-cell survival

Increased osteoblast activity in sclerostin-knockout (Sost(-/-)) mice results in generalized hyperostosis and bones with small bone marrow cavities resulting from hyperactive mineralizing osteoblast populations. Hematopoietic cell fate decisions are dependent on their local microenvironment, which contains osteoblast and stromal cell populations that support both hematopoietic stem cell quiescence and facilitate B-cell development. In this study, we investigated whether high bone mass environments affect B-cell development via the utilization of Sost(-/-) mice, a model of sclerosteosis. We found the bone marrow of Sost(-/-) mice to be specifically depleted of B cells because of elevated apoptosis at all B-cell developmental stages. In contrast, B-cell function in the spleen was normal. Sost expression analysis confirmed that Sost is primarily expressed in osteocytes and is not expressed in any hematopoietic lineage, which indicated that the B-cell defects in Sost(-/-) mice are non-cell autonomous, and this was confirmed by transplantation of wild-type (WT) bone marrow into lethally irradiated Sost(-/-) recipients. WT→Sost(-/-) chimeras displayed a reduction in B cells, whereas reciprocal Sost(-/-) →WT chimeras did not, supporting the idea that the Sost(-/-) bone environment cannot fully support normal B-cell development. Expression of the pre-B-cell growth stimulating factor, Cxcl12, was significantly lower in bone marrow stromal cells of Sost(-/-) mice, whereas the Wnt target genes Lef-1 and Ccnd1 remained unchanged in B cells. Taken together, these results demonstrate a novel role for Sost in the regulation of bone marrow environments that support B cells.

TGF-β regulates sclerostin expression via the ECR5 enhancer

Wnt signaling is critical for skeletal development and homeostasis. Sclerostin (Sost) has emerged as a potent inhibitor of Wnt signaling and, thereby, bone formation. Thus, strategies to reduce sclerostin expression may be used to treat osteoporosis or non-union fractures. Transforming growth factor-beta (TGF-β) elicits various effects upon the skeleton both in vitro and in vivo depending on the duration and timing of administration. In vitro and in vivo studies demonstrate that TGF-β increases osteoprogenitor differentiation but decreases matrix mineralization of committed osteoblasts. Because sclerostin decreases matrix mineralization, this study aimed to examine whether TGF-β achieves such inhibitory effects via transcriptional modulation of Sost. Using the UMR106.01 mature osteoblast cell line, we demonstrated that TGF-βTGF-β(1)-β(2)-β(3) and Activin A increase Sost transcript expression. Pharmacologic inhibition of Alk4/5/7 in vitro and in vivo decreased endogenous Sost expression, and siRNA against Alk4 and Alk5 demonstrated their requirement for endogenous Sost expression. TGF-β(1) targeted the Sost bone enhancer ECR5 and did not affect the transcriptional activity of the endogenous Sost promoter. These results indicate that TGF-β(1) controls Sost transcription in mature osteoblasts, suggesting that sclerostin may mediate the inhibitory effect of TGF-β upon osteoblast differentiation.

Global gene expression analysis of murine limb development

Detailed information about stage-specific changes in gene expression is crucial for understanding the gene regulatory networks underlying development and the various signal transduction pathways contributing to morphogenesis. Here we describe the global gene expression dynamics during early murine limb development, when cartilage, tendons, muscle, joints, vasculature and nerves are specified and the musculoskeletal system of limbs is established. We used whole-genome microarrays to identify genes with differential expression at 5 stages of limb development (E9.5 to 13.5), during fore- and hind-limb patterning. We found that the onset of limb formation is characterized by an up-regulation of transcription factors, which is followed by a massive activation of genes during E10.5 and E11.5 which levels off at later time points. Among the 3520 genes identified as significantly up-regulated in the limb, we find ∼30% to be novel, dramatically expanding the repertoire of candidate genes likely to function in the limb. Hierarchical and stage-specific clustering identified expression profiles that are likely to correlate with functional programs during limb development and further characterization of these transcripts will provide new insights into specific tissue patterning processes. Here, we provide for the first time a comprehensive analysis of developmentally regulated genes during murine limb development, and provide some novel insights into the expression dynamics governing limb morphogenesis.

Prostaglandin E2 signals through PTGER2 to regulate sclerostin expression

The Wnt signaling pathway is a robust regulator of skeletal homeostasis. Gain-of-function mutations promote high bone mass, whereas loss of Lrp5 or Lrp6 co-receptors decrease bone mass. Similarly, mutations in antagonists of Wnt signaling influence skeletal integrity, in an inverse relation to Lrp receptor mutations. Loss of the Wnt antagonist Sclerostin (Sost) produces the generalized skeletal hyperostotic condition of sclerosteosis, which is characterized by increased bone mass and density due to hyperactive osteoblast function. Here we demonstrate that prostaglandin E₂ (PGE₂), a paracrine factor with pleiotropic effects on osteoblasts and osteoclasts, decreases Sclerostin expression in osteoblastic UMR106.01 cells. Decreased Sost expression correlates with increased expression of Wnt/TCF target genes Axin2 and Tcf3. We also show that the suppressive effect of PGE₂ is mediated through a cyclic AMP/PKA pathway. Furthermore, selective agonists for the PGE₂ receptor EP2 mimic the effect of PGE₂ upon Sost, and siRNA reduction in Ptger2 prevents PGE₂-induced Sost repression. These results indicate a functional relationship between prostaglandins and the Wnt/β-catenin signaling pathway in bone.

Hypoxia decreases sclerostin expression and increases Wnt signaling in osteoblasts

Mutations in sclerostin function or expression cause sclerosing bone dysplasias, involving decreased antagonism of Wnt/Lrp5 signaling. Conversely, deletion of the VHL tumor suppressor in osteoblasts, which stabilize HIF-alpha isoforms and thereby enables HIF-alpha/beta-driven gene transcription, increases bone mineral content and cross-sectional area compared to wild-type controls. We examined the influence of cellular hypoxia (1% oxygen) upon sclerostin expression and canonical Wnt signaling. Osteoblasts and osteocytes cultured under hypoxia revealed decreased sclerostin transcript and protein, and increased expression and nuclear localization of activated beta-catenin. Similarly, both hypoxia and the hypoxia mimetic DFO increased beta-catenin gene reporter activity. Hypoxia and its mimetics increased expression of the BMP antagonists gremlin and noggin and decreased Smad-1/5/8 phosphorylation. As a partial explanation for the mechanism of regulation of sclerostin by oxygen, MEF2 reporter assays revealed decreased activity. Modulation of VEGF signaling under normoxia or hypoxia revealed no influence upon Sost transcription. These data suggest that hypoxia inhibits sclerostin expression, through enhanced antagonism of BMP signaling independent of VEGF.

Genetic evidence that SOST inhibits WNT signaling in the limb

SOST is a negative regulator of bone formation, and mutations in human SOST are responsible for sclerosteosis. In addition to high bone mass, sclerosteosis patients occasionally display hand defects, suggesting that SOST may function embryonically. Here we report that overexpression of SOST leads to loss of posterior structures of the zeugopod and autopod by perturbing anterior-posterior and proximal-distal signaling centers in the developing limb. Mutant mice that overexpress SOST in combination with Grem1 and Lrp6 mutations display more severe limb defects than single mutants alone, while Sost(-/-) significantly rescues the Lrp6(-/-) skeletal phenotype, signifying that SOST gain-of-function impairs limb patterning by inhibiting the WNT signaling through LRP5/6.

Human variation in short regions predisposed to deep evolutionary conservation

The landscape of the human genome consists of millions of short islands of conservation that are 100% conserved across multiple vertebrate genomes (termed “bricks”), the majority of which are located in noncoding regions. Several hundred thousand bricks are deeply conserved reaching the genomes of amphibians and fish. Deep phylogenetic conservation of noncoding DNA has been reported to be strongly associated with the presence of gene regulatory elements, introducing bricks as a proxy to the functional noncoding landscape of the human genome. Here, we report a significant overrepresentation of bricks in the promoters of transcription factors and developmental genes, where the high level of phylogenetic conservation correlates with an increase in brick overrepresentation. We also found that the presence of a brick dictates a predisposition to evolutionary constraint, with only 0.7% of the amniota brick central nucleotides being diverged within the primate lineage—an 11-fold reduction in the divergence rate compared with random expectation. Human single-nucleotide polymorphism (SNP) data explains only 3% of primate-specific variation in amniota bricks, thus arguing for a widespread fixation of brick mutations within the primate lineage and prior to human radiation. This variation, in turn, might have been utilized as a driving force for primate- and hominoid-specific adaptation. We also discovered a pronounced deviation from the evolutionary predisposition in the human lineage, with over 20-fold increase in the substitution rate at brick SNP sites over expected values. In addition, contrary to typical brick mutations, brick variation commonly encountered in the human population displays limited, if any, signatures of negative selection as measured by the minor allele frequency and population differentiation (F-statistical measure) measures. These observations argue for the plasticity of gene regulatory mechanisms in vertebrates—with evidence of strong purifying selection acting on the gene regulatory landscape of the human genome, where widespread advantageous mutations in putative regulatory elements are likely utilized in functional diversification and adaptation of species.

Genomic identification of regulatory elements by evolutionary sequence comparison and functional analysis

Despite remarkable recent advances in genomics that have enabled us to identify most of the genes in the human genome, comparable efforts to define transcriptional cis-regulatory elements that control gene expression are lagging behind. The difficulty of this task stems from two equally important problems: our knowledge of how regulatory elements are encoded in genomes remains elementary, and there is a vast genomic search space for regulatory elements, since most of mammalian genomes are noncoding. Comparative genomic approaches are having a remarkable impact on the study of transcriptional regulation in eukaryotes and currently represent the most efficient and reliable methods of predicting noncoding sequences likely to control the patterns of gene expression. By subjecting eukaryotic genomic sequences to computational comparisons and subsequent experimentation, we are inching our way toward a more comprehensive catalog of common regulatory motifs that lie behind fundamental biological processes. We are still far from comprehending how the transcriptional regulatory code is encrypted in the human genome and providing an initial global view of regulatory gene networks, but collectively, the continued development of comparative and experimental approaches will rapidly expand our knowledge of the transcriptional regulome.

Control of the SOST bone enhancer by PTH using MEF2 transcription factors

Expression of the osteocyte-derived bone formation inhibitor sclerostin in adult bone requires a distant enhancer. We show that MEF2 transcription factors control this enhancer and mediate inhibition of sclerostin expression by PTH.

INTRODUCTION

Sclerostin encoded by the SOST gene is a key regulator of bone formation. Lack of SOST expression is the cause for the progressive bone overgrowth disorders sclerosteosis and Van Buchem disease. We have previously identified a distant enhancer within the 52-kb Van Buchem disease deletion downstream of the SOST gene that is essential for its expression in adult bone. Furthermore, we and others have reported that SOST expression is suppressed by PTH. The aim of this study was to identify transcription factors involved in SOST bone enhancer activity and mediating PTH responsiveness.

MATERIALS AND METHODS

Regulation of the SOST enhancer and promoter was studied by luciferase reporter gene assays. Transcription factor binding sites were mapped by footprint analysis and functional mutation analyses using transient transfections of osteoblast-like UMR-106 cells that exhibit endogenous SOST expression. Specific transcription factor binding was predicted by sequence analysis and shown by gel retardation assays and antibody-induced supershifts. Expression of myocyte enhancer factors 2 (MEF2) was detected by in situ hybridization, quantitative RT-PCR (qPCR), and immunohistochemistry. The role of MEF2s in SOST expression was assessed by reporter gene assays and siRNA-mediated RNA knockdown.

RESULTS

PTH completely suppressed the transcriptional activity of the SOST bone enhancer but did not affect the SOST promoter. A MEF2 response element was identified in the bone enhancer. It was essential for transcriptional activation, bound MEF2 transcription factors, and mediated PTH responsiveness. Expression of MEF2s in bone was shown by qPCR, in situ hybridization, and immunohistochemistry. MEF2s and sclerostin co-localized in osteocytes. Enhancer activity was stimulated by MEF2C overexpression and inhibited by co-expression of a dominant negative MEF2C mutant. Finally, siRNA-mediated knockdown of MEF2A, C, and D suppressed endogenous SOST expression in UMR-106 osteoblast-like cells.

CONCLUSIONS

These data strongly suggest that SOST expression in osteocytes of adult bone and its inhibition by PTH is mediated by MEF2A, C, and D transcription factors controlling the SOST bone enhancer. Hence, MEF2s are implicated in the regulation of adult bone mass.

Predicting tissue-specific enhancers in the human genome

Determining how transcriptional regulatory signals are encoded in vertebrate genomes is essential for understanding the origins of multicellular complexity; yet the genetic code of vertebrate gene regulation remains poorly understood. In an attempt to elucidate this code, we synergistically combined genome-wide gene-expression profiling, vertebrate genome comparisons, and transcription factor binding-site analysis to define sequence signatures characteristic of candidate tissue-specific enhancers in the human genome. We applied this strategy to microarray-based gene expression profiles from 79 human tissues and identified 7187 candidate enhancers that defined their flanking gene expression, the majority of which were located outside of known promoters. We cross-validated this method for its ability to de novo predict tissue-specific gene expression and confirmed its reliability in 57 of the 79 available human tissues, with an average precision in enhancer recognition ranging from 32% to 63% and a sensitivity of 47%. We used the sequence signatures identified by this approach to successfully assign tissue-specific predictions to approximately 328,000 human-mouse conserved noncoding elements in the human genome. By overlapping these genome-wide predictions with a data set of enhancers validated in vivo, in transgenic mice, we were able to confirm our results with a 28% sensitivity and 50% precision. These results indicate the power of combining complementary genomic data sets as an initial computational foray into a global view of tissue-specific gene regulation in vertebrates.

Modifying yeast artificial chromosomes to generate Cre/LoxP and FLP/FRT site-specific deletions and inversions

The ability to efficiently and accurately modify genomic DNA through targeted and tissue-specific mutations is an important goal in animal transgenesis. Here we describe how to exploit two systems of homologous recombination, from yeast and bacteria, to engineer yeast artificial chromosomes (YACs) to generate targeted deletions and inversions in vivo, in transgenic animals, and in the presence of DNA-modifying enzymes known as recombinases. Through homologous recombination in yeast, specific recombinogenic sequences are inserted upstream and downstream of a region in the YAC. The sites of integration of these short sequence elements are chosen carefully, such that the YAC is left functionally intact, and this modified transgene represents the wild-type allele. This YAC is subsequently used to generate transgenic animals, which when bred to animals expressing recombinase proteins result in genetic modifications. By expressing recombinase proteins from different tissue-specific promoters, one can mediate site-specific recombination to generate either ubiquitous or tissue-specific deletions or inversion. These modifications can then be carried through the germline or can be studied somatically. A great advantage of this system is the ability to evaluate subtle genetic effects independent of position-effect variegation, and transgene copy number, eliminating the need to examine several independently generated lines of transgenic animals for each genetic variant.

Mulan: multiple-sequence alignment to predict functional elements in genomic sequences

Multiple sequence alignment analysis is a powerful approach for translating the evolutionary selective power into phylogenetic relationships to localize functional coding and noncoding genomic elements. The tool Mulan (http://mulan.dcode.org/) has been designed to effectively perform multiple comparisons of genomic sequences necessary to facilitate bioinformatic-driven biological discoveries. The Mulan network server is capable of comparing both closely and distantly related genomes to identify conserved elements over a broad range of evolutionary time. Several novel algorithms are brought together in this tool: the tba multisequence aligner program used to rapidly identify local sequence conservation and the multiTF program to detect evolutionarily conserved transcription factor binding sites in alignments. Mulan is integrated with the ERC Browser, the UCSC Genome Browser for quick uploads of available sequences and supports two-way communication with the GALA database to overlay GALA functional genome annotation with sequence conservation profiles. Local multiple alignments computed by Mulan ensure reliable representation of short- and large-scale genomic rearrangements in distant organisms. Recently, we have also introduced the ability to handle duplications to permit the reliable reconstruction of evolutionary events that underlie the genome sequence data. Here, we describe the main features of the Mulan tool that include the interactive modification of critical conservation parameters, visualization options, and dynamic access to sequence data from visual graphs for flexible and easy-to-perform analysis of differentially evolving genomic regions.

Array2BIO: from microarray expression data to functional annotation of co-regulated genes

BACKGROUND

There are several isolated tools for partial analysis of microarray expression data. To provide an integrative, easy-to-use and automated toolkit for the analysis of Affymetrix microarray expression data we have developed Array2BIO, an application that couples several analytical methods into a single web based utility.

RESULTS

Array2BIO converts raw intensities into probe expression values, automatically maps those to genes, and subsequently identifies groups of co-expressed genes using two complementary approaches: (1) comparative analysis of signal versus control and (2) clustering analysis of gene expression across different conditions. The identified genes are assigned to functional categories based on Gene Ontology classification and KEGG protein interaction pathways. Array2BIO reliably handles low-expressor genes and provides a set of statistical methods for quantifying expression levels, including Benjamini-Hochberg and Bonferroni multiple testing corrections. An automated interface with the ECR Browser provides evolutionary conservation analysis for the identified gene loci while the interconnection with Crème allows prediction of gene regulatory elements that underlie observed expression patterns.

CONCLUSION

We have developed Array2BIO - a web based tool for rapid comprehensive analysis of Affymetrix microarray expression data, which also allows users to link expression data to Dcode.org comparative genomics tools and integrates a system for translating co-expression data into mechanisms of gene co-regulation. Array2BIO is publicly available at http://array2bio.dcode.org.

Dcode.org anthology of comparative genomic tools

Comparative genomics provides the means to demarcate functional regions in anonymous DNA sequences. The successful application of this method to identifying novel genes is currently shifting to deciphering the non-coding encryption of gene regulation across genomes. To facilitate the practical application of comparative sequence analysis to genetics and genomics, we have developed several analytical and visualization tools for the analysis of arbitrary sequences and whole genomes. These tools include two alignment tools, zPicture and Mulan; a phylogenetic shadowing tool, eShadow for identifying lineage- and species-specific functional elements; two evolutionary conserved transcription factor analysis tools, rVista and multiTF; a tool for extracting cis-regulatory modules governing the expression of co-regulated genes, Creme 2.0; and a dynamic portal to multiple vertebrate and invertebrate genome alignments, the ECR Browser. Here, we briefly describe each one of these tools and provide specific examples on their practical applications. All the tools are publicly available at the http://www.dcode.org/ website.

Strategies for characterising cis-regulatory elements in Xenopus

Understanding the cis-regulatory architecture of metazoan organisms is the greatest challenge facing genome biology today. In vertebrate organisms, distinct sequence elements mediate transcriptional regulation and are scattered throughout the genome, either proximal or distal to promoters. The identification of transcriptional enhancers has proven rather difficult by conventional experimental approaches. In the past decade, the rapid generation of genomic sequences for multiple vertebrate organisms, accompanied by sophisticated comparative tools, has facilitated the identification of non-coding evolutionarily conserved regions that may encode cis-regulatory elements. Validating computational predictions and characterising cis-regulatory elements in vivo, however, has been a major bottleneck, mainly because the most commonly used organism for these experiments has been the mouse, and generating transgenic mice or modifying the mouse genome continues to be a labour-intensive, low-throughput, expensive process. This has led to the use of Xenopus, which holds great promise for high-throughput interrogation of putative cis-regulatory elements. In particular, Xenopus tropicalis may become particularly powerful for elucidating regulatory networks, chiefly because it is amenable to genetic manipulations, and its genome is being sequenced.

Genomic deletion of a long-range bone enhancer misregulates sclerostin in Van Buchem disease

Mutations in distant regulatory elements can have a negative impact on human development and health, yet because of the difficulty of detecting these critical sequences, we predominantly focus on coding sequences for diagnostic purposes. We have undertaken a comparative sequence-based approach to characterize a large noncoding region deleted in patients affected by Van Buchem (VB) disease, a severe sclerosing bone dysplasia. Using BAC recombination and transgenesis, we characterized the expression of human sclerostin (SOST) from normal (SOST(wt)) or Van Buchem (SOST(vbDelta) alleles. Only the SOST(wt) allele faithfully expressed high levels of human SOST in the adult bone and had an impact on bone metabolism, consistent with the model that the VB noncoding deletion removes a SOST-specific regulatory element. By exploiting cross-species sequence comparisons with in vitro and in vivo enhancer assays, we were able to identify a candidate enhancer element that drives human SOST expression in osteoblast-like cell lines in vitro and in the skeletal anlage of the embryonic day 14.5 (E14.5) mouse embryo, and discovered a novel function for sclerostin during limb development. Our approach represents a framework for characterizing distant regulatory elements associated with abnormal human phenotypes.

Evolution and functional classification of vertebrate gene deserts

Large tracts of the human genome, known as gene deserts, are devoid of protein-coding genes. Dichotomy in their level of conservation with chicken separates these regions into two distinct categories, stable and variable. The separation is not caused by differences in rates of neutral evolution but instead appears to be related to different biological functions of stable and variable gene deserts in the human genome. Gene Ontology categories of the adjacent genes are strongly biased toward transcriptional regulation and development for the stable gene deserts, and toward distinctively different functions for the variable gene deserts. Stable gene deserts resist chromosomal rearrangements and appear to harbor multiple distant regulatory elements physically linked to their neighboring genes, with the linearity of conservation invariant throughout vertebrate evolution.

Mulan: multiple-sequence local alignment and visualization for studying function and evolution

Multiple-sequence alignment analysis is a powerful approach for understanding phylogenetic relationships, annotating genes, and detecting functional regulatory elements. With a growing number of partly or fully sequenced vertebrate genomes, effective tools for performing multiple comparisons are required to accurately and efficiently assist biological discoveries. Here we introduce Mulan (http://mulan.dcode.org/), a novel method and a network server for comparing multiple draft and finished-quality sequences to identify functional elements conserved over evolutionary time. Mulan brings together several novel algorithms: the TBA multi-aligner program for rapid identification of local sequence conservation, and the multiTF program for detecting evolutionarily conserved transcription factor binding sites in multiple alignments. In addition, Mulan supports two-way communication with the GALA database; alignments of multiple species dynamically generated in GALA can be viewed in Mulan, and conserved transcription factor binding sites identified with Mulan/multiTF can be integrated and overlaid with extensive genome annotation data using GALA. Local multiple alignments computed by Mulan ensure reliable representation of short- and large-scale genomic rearrangements in distant organisms. Mulan allows for interactive modification of critical conservation parameters to differentially predict conserved regions in comparisons of both closely and distantly related species. We illustrate the uses and applications of the Mulan tool through multispecies comparisons of the GATA3 gene locus and the identification of elements that are conserved in a different way in avians than in other genomes, allowing speculation on the evolution of birds. Source code for the aligners and the aligner-evaluation software can be freely downloaded from http://www.bx.psu.edu/miller_lab/.

Interpreting mammalian evolution using Fugu genome comparisons

Recently, it has been shown that a significant number of evolutionarily conserved human-Fugu noncoding elements function as tissue-specific transcriptional enhancers in vivo, suggesting that distant comparisons are capable of identifying a particular class of regulatory elements. We therefore hypothesized that by juxtaposing human/Fugu and human/mouse conservation patterns we can define conservation criteria for discovering transcriptional regulatory elements specific to mammals. Genome-scale comparisons of noncoding human/Fugu evolutionary conserved elements (ECRs) and their humans/mouse counterparts revealed a particular signature common to human/mouse ECRs (>or=350 bp long, >or=77% identity) that are also conserved in fishes. This newly defined threshold identifies 90% of all human/Fugu noncoding ECRs without the assistance of human-Fugu genome alignments and provides a very efficient filter for identifying functional human/mouse ECRs.

CREME: Cis-Regulatory Module Explorer for the human genome

The binding of transcription factors to specific regulatory sequence elements is a primary mechanism for controlling gene transcription. Eukaryotic genes are often regulated by several transcription factors whose binding sites are tightly clustered and form cis-regulatory modules. In this paper, we present a web server, CREME, for identifying and visualizing cis-regulatory modules in the promoter regions of a given set of potentially co-regulated genes. CREME relies on a database of putative transcription factor binding sites that have been annotated across the human genome using a library of position weight matrices and evolutionary conservation with the mouse and rat genomes. A search algorithm is applied to this data set to identify combinations of transcription factors whose binding sites tend to co-occur in close proximity in the promoter regions of the input gene set. The identified cis-regulatory modules are statistically scored and significant combinations are reported and graphically visualized. Our web server is available at http://creme.dcode.org.

rVISTA 2.0: evolutionary analysis of transcription factor binding sites

Identifying and characterizing the transcription factor binding site (TFBS) patterns of cis-regulatory elements represents a challenge, but holds promise to reveal the regulatory language the genome uses to dictate transcriptional dynamics. Several studies have demonstrated that regulatory modules are under positive selection and, therefore, are often conserved between related species. Using this evolutionary principle, we have created a comparative tool, rVISTA, for analyzing the regulatory potential of noncoding sequences. Our ability to experimentally identify functional noncoding sequences is extremely limited, therefore, rVISTA attempts to fill this great gap in genomic analysis by offering a powerful approach for eliminating TFBSs least likely to be biologically relevant. The rVISTA tool combines TFBS predictions, sequence comparisons and cluster analysis to identify noncoding DNA regions that are evolutionarily conserved and present in a specific configuration within genomic sequences. Here, we present the newly developed version 2.0 of the rVISTA tool, which can process alignments generated by both the zPicture and blastz alignment programs or use pre-computed pairwise alignments of several vertebrate genomes available from the ECR Browser and GALA database. The rVISTA web server is closely interconnected with the TRANSFAC database, allowing users to either search for matrices present in the TRANSFAC library collection or search for user-defined consensus sequences. The rVISTA tool is publicly available at http://rvista.dcode.org/.

ECR Browser: a tool for visualizing and accessing data from comparisons of multiple vertebrate genomes

With an increasing number of vertebrate genomes being sequenced in draft or finished form, unique opportunities for decoding the language of DNA sequence through comparative genome alignments have arisen. However, novel tools and strategies are required to accommodate this large volume of genomic information and to facilitate the transfer of predictions generated by comparative sequence alignment to researchers focused on experimental annotation of genome function. Here, we present the ECR Browser, a tool that provides easy and dynamic access to whole genome alignments of human, mouse, rat and fish sequences. This web-based tool (http://ecrbrowser.dcode.org) provides the starting point for discovery of novel genes, identification of distant gene regulatory elements and prediction of transcription factor binding sites. The genome alignment portal of the ECR Browser also permits fast and automated alignments of any user-submitted sequence to the genome of choice. The interconnection of the ECR Browser with other DNA sequence analysis tools creates a unique portal for studying and exploring vertebrate genomes.

eShadow: a tool for comparing closely related sequences

Primate sequence comparisons are difficult to interpret due to the high degree of sequence similarity shared between such closely related species. Recently, a novel method, phylogenetic shadowing, has been pioneered for predicting functional elements in the human genome through the analysis of multiple primate sequence alignments. We have expanded this theoretical approach to create a computational tool, eShadow, for the identification of elements under selective pressure in multiple sequence alignments of closely related genomes, such as in comparisons of human-to-primate or mouse-to-rat DNA. This tool integrates two different statistical methods and allows for the dynamic visualization of the resulting conservation profile. eShadow also includes a versatile optimization module capable of training the underlying Hidden Markov Model to differentially predict functional sequences. This module grants the tool high flexibility in the analysis of multiple sequence alignments and in comparing sequences with different divergence rates. Here, we describe the eShadow comparative tool and its potential uses for analyzing both multiple nucleotide and protein alignments to predict putative functional elements.

zPicture: dynamic alignment and visualization tool for analyzing conservation profiles

Comparative sequence analysis has evolved as an essential technique for identifying functional coding and noncoding elements conserved throughout evolution. Here, we introduce zPicture, an interactive Web-based sequence alignment and visualization tool for dynamically generating conservation profiles and identifying evolutionarily conserved regions (ECRs). zPicture is highly flexible, because critical parameters can be modified interactively, allowing users to differentially predict ECRs in comparisons of sequences of different phylogenetic distances and evolutionary rates. We demonstrate the application of this module to identify a known regulatory element in the HOXD locus, in which functional ECRs are difficult to discern against the highly conserved genomic background. zPicture also facilitates transcription factor binding-site analysis via the rVista tool portal. We present an example of the HBB complex when zPicture/rVista combination specifically pinpoints to two ECRs containing GATA-1, NF-E2, and TAL1/E47 binding sites that were identified previously as transcriptional enhancers. In addition, zPicture is linked to the UCSC Genome Browser, allowing users to automatically extract sequences and gene annotations for any recorded locus. Finally, we describe how this tool can be efficiently applied to the analysis of nonvertebrate genomes, including those of microbial organisms.