Immunogenetic characteristics of immunoglobulin E in allergic disease

Patients with allergic diseases produce an excess of allergen-specific IgE, the specific effector molecule that triggers allergic reactions. The provocation for this excess IgE production is still uncertain. Current ideas include oligoclonal expansion of allergen-specific B cells emanating from germinal centres, activation by superantigen of a subset of B cells, or polyclonal B cells class switching to IgE due to an IL-4 predominance. Additionally, genetic elements contribute to a propensity for increased allergen-specific IgE production. The procedure of RT-PCR allows for amplification of infrequent IgE mRNA transcripts from B cells of atopic individuals, and so facilitates examination of expressed Ig cDNA sequences. Better knowledge of the molecular characteristics of IgE produced by patients with allergic diseases would elucidate the immunogenetic basis for elevated allergen-specific IgE levels. The 'immunogenetic footprint' of IgE transcripts may elucidate the origin and activation of IgE-producing B cells in allergic disease. Here we review studies of the immunogenetic features of IgE in allergic diseases, highlighting the major advances and the experimental limitations.

Genome-wide association study identifies five new schizophrenia loci

We examined the role of common genetic variation in schizophrenia in a genome-wide association study of substantial size: a stage 1 discovery sample of 21,856 individuals of European ancestry and a stage 2 replication sample of 29,839 independent subjects. The combined stage 1 and 2 analysis yielded genome-wide significant associations with schizophrenia for seven loci, five of which are new (1p21.3, 2q32.3, 8p23.2, 8q21.3 and 10q24.32-q24.33) and two of which have been previously implicated (6p21.32-p22.1 and 18q21.2). The strongest new finding (P = 1.6 × 10 -11) was with rs1625579 within an intron of a putative primary transcript for MIR137 (microRNA 137), a known regulator of neuronal development. Four other schizophrenia loci achieving genome-wide significance contain predicted targets of MIR137, suggesting MIR137-mediated dysregulation as a previously unknown etiologic mechanism in schizophrenia. In a joint analysis with a bipolar disorder sample (16,374 affected individuals and 14,044 controls), three loci reached genome-wide significance: CACNA1C (rs4765905, P = 7.0 × 10 -9), ANK3 (rs10994359, P = 2.5 × 10 -8) and the ITIH3-ITIH4 region (rs2239547, P = 7.8 × 10 -9).

A promoter-level mammalian expression atlas

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

Association of FOXE1 Polyalanine Repeat Region with Papillary Thyroid Cancer

CONTEXT: Polyalanine tract variations in transcription factors have been identified for a wide spectrum of developmental disorders. The thyroid transcription factor forkhead factor E1 (FOXE1) contains a polymorphic polyalanine tract with 12-22 alanines. Single-nucleotide polymorphisms (SNP) close to this locus are associated with papillary thyroid cancer (PTC), and a strong linkage disequilibrium block extends across this region. OBJECTIVE: The objective of the study was to assess whether the FOXE1 polyalanine repeat region was associated with PTC and to assess the effect of polyalanine repeat region variants on protein expression, DNA binding, and transcriptional function on FOXE1-responsive promoters. DESIGN: This was a case-control study. SETTING: The study was conducted at a tertiary referral hospital. PATIENTS AND METHODS: The FOXE1 polyalanine repeat region and tag SNP were genotyped in 70 PTC, with a replication in a further 92 PTC, and compared with genotypes in 5767 healthy controls (including 5667 samples from the Wellcome Trust Case Control Consortium). In vitro studies were performed to examine the protein expression, DNA binding, and transcriptional function for FOXE1 variants of different polyalanine tract lengths. RESULTS: All the genotyped SNP were in tight linkage disequilibrium, including the FOXE1 polyalanine repeat region. We confirmed the strong association of rs1867277 with PTC (overall P = 1 × 10(-7), odds ratio 1.84, confidence interval 1.31-2.57). rs1867277 was in tight linkage disequilibrium with the FOXE1 polyalanine repeat region (r(2) = 0.95). FOXE1(16Ala) was associated with PTC with an odds ratio of 2.23 (confidence interval 1.42-3.50; P = 0.0005). Functional studies in vitro showed that FOXE1(16Ala) was transcriptionally impaired compared with FOXE1(14Ala), which was not due to differences in protein expression or DNA binding. CONCLUSIONS: We have confirmed the previous association of FOXE1 with PTC. Our data suggest that the coding polyalanine expansion in FOXE1 may be responsible for the observed association between FOXE1 and PTC.

An R2R3 MYB transcription factor determines red petal colour in an Actinidia (kiwifruit) hybrid population

Background Red colour in kiwifruit results from the presence of anthocyanin pigments. Their expression, however, is complex, and varies among genotypes, species, tissues and environments. An understanding of the biosynthesis, physiology and genetics of the anthocyanins involved, and the control of their expression in different tissues, is required. A complex, the MBW complex, consisting of R2R3-MYB and bHLH transcription factors together with a WD-repeat protein, activates anthocyanin 3-O-galactosyltransferase (F3GT1) to produce anthocyanins. We examined the expression and genetic control of anthocyanins in flowers of Actinidia hybrid families segregating for red and white petal colour. Results Four inter-related backcross families between Actinidia chinensis Planch. var. chinensis and Actinidia eriantha Benth. were identified that segregated 1:1 for red or white petal colour. Flower pigments consisted of five known anthocyanins (two delphinidin-based and three cyanidin-based) and three unknowns. Intensity and hue differed in red petals from pale pink to deep magenta, and while intensity of colour increased with total concentration of anthocyanin, no association was found between any particular anthocyanin data and hue. Real time qPCR demonstrated that an R2R3 MYB, MYB110a, was expressed at significant levels in red-petalled progeny, but not in individuals with white petals. A microsatellite marker was developed that identified alleles that segregated with red petal colour, but not with ovary, stamen filament, or fruit flesh colour in these families. The marker mapped to chromosome 10 in Actinidia. The white petal phenotype was complemented by syringing Agrobacterium tumefaciens carrying Actinidia 35S::MYB110a into the petal tissue. Red pigments developed in white petals both with, and without, co-transformation with Actinidia bHLH partners. MYB110a was shown to directly activate Actinidia F3GT1 in transient assays. Conclusions The transcription factor, MYB110a, regulates anthocyanin production in petals in this hybrid population, but not in other flower tissues or mature fruit. The identification of delphinidin-based anthocyanins in these flowers provides candidates for colour enhancement in novel fruits.

Membrane associated transporter protein gene (SLC45A2) and the genetic basis of normal human pigmentation variation

This work is concerned with the genetic basis of normal human pigmentation variation. Specifically, the role of polymorphisms within the solute carrier family 45 member 2 (SLC45A2 or membrane associated transporter protein; MATP) gene were investigated with respect to variation in hair, skin and eye colour ― both between and within populations. SLC45A2 is an important regulator of melanin production and mutations in the gene underly the most recently identified form of oculocutaneous albinism. There is evidence to suggest that non-synonymous polymorphisms in SLC45A2 are associated with normal pigmentation variation between populations. Therefore, the underlying hypothesis of this thesis is that polymorphisms in SLC45A2 will alter the function or regulation of the protein, thereby altering the important role it plays in melanogenesis and providing a mechanism for normal pigmentation variation. In order to investigate the role that SLC45A2 polymorphisms play in human pigmentation variation, a DNA database was established which collected pigmentation phenotypic information and blood samples of more than 700 individuals. This database was used as the foundation for two association studies outlined in this thesis, the first of which involved genotyping two previously-described non-synonymous polymorphisms, p.Glu272Lys and p.Phe374Leu, in four different population groups. For both polymorphisms, allele frequencies were significantly different between population groups and the 272Lys and 374Leu alleles were strongly associated with black hair, brown eyes and olive skin colour in Caucasians. This was the first report to show that SLC45A2 polymorphisms were associated with normal human intra-population pigmentation variation. The second association study involved genotyping several SLC45A2 promoter polymorphisms to determine if they also played a role in pigmentation variation. Firstly, the transcription start site (TSS), and hence putative proximal promoter region, was identified using 5' RNA ligase mediated rapid amplification of cDNA ends (RLM-RACE). Two alternate TSSs were identified and the putative promoter region was screened for novel polymorphisms using denaturing high performance liquid chromatography (dHPLC). A novel duplication (c.–1176_–1174dupAAT) was identified along with other previously described single nucleotide polymorphisms (c.–1721C>G and c.–1169G>A). Strong linkage disequilibrium ensured that all three polymorphisms were associated with skin colour such that the –1721G, +dup and –1169A alleles were associated with olive skin in Caucasians. No linkage disequilibrium was observed between the promoter and coding region polymorphisms, suggesting independent effects. The association analyses were complemented with functional data, showing that the –1721G, +dup and –1169A alleles significantly decreased SLC45A2 transcriptional activity. Based on in silico bioinformatic analysis that showed these alleles remove a microphthalmia-associated transcription factor (MITF) binding site, and that MITF is a known regulator of SLC45A2 (Baxter and Pavan, 2002; Du and Fisher, 2002), it was postulated that SLC45A2 promoter polymorphisms could contribute to the regulation of pigmentation by altering MITF binding affinity. Further characterisation of the SLC45A2 promoter was carried out using luciferase reporter assays to determine the transcriptional activity of different regions of the promoter. Five constructs were designed of increasing length and their promoter activity evaluated. Constitutive promoter activity was observed within the first ~200 bp and promoter activity increased as the construct size increased. The functional impact of the –1721G, +dup and –1169A alleles, which removed a MITF consensus binding site, were assessed using electrophoretic mobility shift assays (EMSA) and expression analysis of genotyped melanoblast and melanocyte cell lines. EMSA results confirmed that the promoter polymorphisms affected DNA-protein binding. Interestingly, however, the protein/s involved were not MITF, or at least MITF was not the protein directly binding to the DNA. In an effort to more thoroughly characterise the functional consequences of SLC45A2 promoter polymorphisms, the mRNA expression levels of SLC45A2 and MITF were determined in melanocyte/melanoblast cell lines. Based on SLC45A2’s role in processing and trafficking TYRP1 from the trans-Golgi network to stage 2 melanosmes, the mRNA expression of TYRP1 was also investigated. Expression results suggested a coordinated expression of pigmentation genes. This thesis has substantially contributed to the field of pigmentation by showing that SLC45A2 polymorphisms not only show allele frequency differences between population groups, but also contribute to normal pigmentation variation within a Caucasian population. In addition, promoter polymorphisms have been shown to have functional consequences for SLC45A2 transcription and the expression of other pigmentation genes. Combined, the data presented in this work supports the notion that SLC45A2 is an important contributor to normal pigmentation variation and should be the target of further research to elucidate its role in determining pigmentation phenotypes. Understanding SLC45A2’s function may lead to the development of therapeutic interventions for oculocutaneous albinism and other disorders of pigmentation. It may also help in our understanding of skin cancer susceptibility and evolutionary adaptation to different UV environments, and contribute to the forensic application of pigmentation phenotype prediction.

The genetic basis of human height : the role of estrogen

Height is a complex physical trait that displays strong heritability. Adult height is related to length of the long bones, which is determined by growth at the epiphyseal growth plate. Longitudinal bone growth occurs via the process of endochondral ossification, where bone forms over the differentiating cartilage template at the growth plate. Estrogen plays a major role in regulating longitudinal bone growth and is responsible for inducing the pubertal growth spurt and fusion of the epiphyseal growth plate. However, the mechanism by which estrogen promotes epiphyseal fusion is poorly understood. It has been hypothesised that estrogen functions to regulate growth plate fusion by stimulating chondrocyte apoptosis, angiogenesis and bone cell invasion in the growth plate. Another theory has suggested that estrogen exposure exhausts the proliferative capacity of growth plate chondrocytes, which accelerates the process of chondrocyte senescence, leading to growth plate fusion. The overall objective of this study was to gain a greater understanding of the molecular mechanisms behind estrogen-mediated growth and height attainment by examining gene regulation in chondrocytes and the role of some of these genes in normal height inheritance. With the heritability of height so well established, the initial hypothesis was that genetic variation in candidate genes associated with longitudinal bone growth would be involved in normal adult height variation. The height-related genes FGFR3, CBFA1, ER and CBFA1 were screened for novel polymorphisms using denaturing HPLC and RFLP analysis. In total, 24 polymorphisms were identified. Two SNPs in ER (rs3757323 C>T and rs1801132 G>C) were strongly associated with adult male height and displayed an 8 cm and 9 cm height difference between homozygous genotypes, respectively. The TC haplotype of these SNPs was associated with a 6 cm decrease in height and remarkably, no homozygous carriers of the TC haplotype were identified in tall subjects. No significant associations with height were found for polymorphisms in the FGFR3, CBFA1 or VDR genes. In the epiphyseal growth plate, chondrocyte proliferation, matrix synthesis and chondrocyte hypertrophy are all major contributors to long bone growth. As estrogen plays such a significant role in both growth and final height attainment, another hypothesis of this study was that estrogen exerted its effects in the growth plate by influencing chondrocyte proliferation and mediating the expression of chondrocyte marker genes. The examination of genes regulated by estrogen in chondrocyte-like cells aimed to identify potential regulators of growth plate fusion, which may further elucidate mechanisms involved in the cessation of linear growth. While estrogen did not dramatically alter the proliferation of the SW1353 cell line, gene expression experiments identified several estrogen regulated genes. Sixteen chondrocyte marker genes were examined in response to estrogen concentrations ranging from 10-12 M to 10-8 M over varying time points. Of the genes analysed, IHH, FGFR3, collagen II and collagen X were not readily detectable and PTHrP, GHR, ER, BMP6, SOX9 and TGF1 mRNAs showed no significant response to estrogen treatments. However, the expression of MMP13, CBFA1, BCL-2 and BAX genes were significantly decreased. Interestingly, the majority of estrogen regulated genes in SW1353 cells are expressed in the hypertrophic zone of the growth plate. Estrogen is also known to regulate systemic GH secretion and local GH action. At the molecular level, estrogen functions to inhibit GH action by negatively regulating GH signalling. GH treated SW1353 cells displayed increases in MMP9 mRNA expression (4.4-fold) and MMP13 mRNA expression (64-fold) in SW1353 cells. Increases were also detected in their respective proteins. Treatment with AG490, an established JAK2 inhibitor, blocked the GH mediated stimulation of both MMP9 and MMP13 mRNA expression. The application of estrogen and GH to SW1353 cells attenuated GH-stimulated MMP13 levels, but did not affect MMP9 levels. Investigation of GH signalling revealed that SW1353 cells have high levels of activated JAK2 and exposure to GH, estrogen, AG490 and other signalling inhibitors did not affect JAK2 phosphorylation. Interestingly, AG490 treatment dramatically decreased ERK2 signalling, although GH did stimulate ERK2 phosphorylation above control levels. AG490 also decreased CBFA1 expression, a transcription factor known to activate MMP9 and MMP13. Finally, GH and estrogen treatment increased expression of SOCS3 mRNA, suggesting that SOCS3 may regulate JAK/STAT signalling in SW1353 cells. The modulation of GH-mediated MMP expression by estrogen in SW1353 cells represents a potentially novel mechanism by which estrogen may regulate longitudinal bone growth. However, further investigation is required in order to elucidate the precise mechanisms behind estrogen and GH regulation of MMP13 expression in SW1353 cells. This study has provided additional evidence that estrogen and the ER gene are major factors in the regulation of growth and the determination of adult height. Newly identified polymorphisms in the ER gene not only contribute to our understanding of the genetic basis of human height, but may also be useful in association studies examining other complex traits. This study also identified several estrogen regulated genes and indicated that estrogen modifies the expression of genes which are primarily expressed in the hypertrophic region of the epiphyseal growth plate. Furthermore, synergistic studies incorporating GH and estrogen have revealed the ability of estrogen to attenuate the effects of GH on MMP13 expression, revealing potential pathways by which estrogen may modulate growth plate fusion, longitudinal bone growth and even arthritis.

Evolving noisy oscillatory dynamics in genetic regulatory networks

We introduce a genetic programming (GP) approach for evolving genetic networks that demonstrate desired dynamics when simulated as a discrete stochastic process. Our representation of genetic networks is based on a biochemical reaction model including key elements such as transcription, translation and post-translational modifications. The stochastic, reaction-based GP system is similar but not identical with algorithmic chemistries. We evolved genetic networks with noisy oscillatory dynamics. The results show the practicality of evolving particular dynamics in gene regulatory networks when modelled with intrinsic noise.

A mathematical model for genetic regulation of the lactose operon

In this paper we construct a mathematical model for the genetic regulatory network of the lactose operon. This mathematical model contains transcription and translation of the lactose permease (LacY) and a reporter gene GFP. The probability of transcription of LacY is determined by 14 binding states out of all 50 possible binding states of the lactose operon based on the quasi-steady-state assumption for the binding reactions, while we calculate the probability of transcription for the reporter gene GFP based on 5 binding states out of 19 possible binding states because the binding site O2 is missing for this reporter gene. We have tested different mechanisms for the transport of thio-methylgalactoside (TMG) and the effect of different Hill coefficients on the simulated LacY expression levels. Using this mathematical model we have realized one of the experimental results with different LacY concentrations, which are induced by different concentrations of TMG.

Association between prostinogen (KLK15) genetic variants and prostate cancer risk and aggressiveness in Australia and a meta-analysis of GWAS data

Background: Kallikrein 15 (KLK15)/Prostinogen is a plausible candidate for prostate cancer susceptibility. Elevated KLK15 expression has been reported in prostate cancer and it has been described as an unfavorable prognostic marker for the disease. Objectives: We performed a comprehensive analysis of association of variants in the KLK15 gene with prostate cancer risk and aggressiveness by genotyping tagSNPs, as well as putative functional SNPs identified by extensive bioinformatics analysis. Methods and Data Sources: Twelve out of 22 SNPs, selected on the basis of linkage disequilibrium pattern, were analyzed in an Australian sample of 1,011 histologically verified prostate cancer cases and 1,405 ethnically matched controls. Replication was sought from two existing genome wide association studies (GWAS): the Cancer Genetic Markers of Susceptibility (CGEMS) project and a UK GWAS study. Results: Two KLK15 SNPs, rs2659053 and rs3745522, showed evidence of association (p, 0.05) but were not present on the GWAS platforms. KLK15 SNP rs2659056 was found to be associated with prostate cancer aggressiveness and showed evidence of association in a replication cohort of 5,051 patients from the UK, Australia, and the CGEMS dataset of US samples. A highly significant association with Gleason score was observed when the data was combined from these three studies with an Odds Ratio (OR) of 0.85 (95% CI = 0.77-0.93; p = 2.7610 24). The rs2659056 SNP is predicted to alter binding of the RORalpha transcription factor, which has a role in the control of cell growth and differentiation and has been suggested to control the metastatic behavior of prostate cancer cells. Conclusions: Our findings suggest a role for KLK15 genetic variation in the etiology of prostate cancer among men of European ancestry, although further studies in very large sample sets are necessary to confirm effect sizes.

Investigation of Association between PFO complicated by cryptogenic stroke and a common variant of the cardiac transcription factor GATA4

Patent foramen ovale (PFO) is associated with clinical conditions including cryptogenic stroke, migraine and varicose veins. Data from studies in humans and mouse suggest that PFO and the secundum form of atrial septal defect (ASDII) exist in an anatomical continuum of septal dysmorphogenesis with a common genetic basis. Mutations in multiple members of the evolutionarily conserved cardiac transcription factor network, including GATA4, cause or predispose to ASDII and PFO. Here, we assessed whether the most prevalent variant of the GATA4 gene, S377G, was significantly associated with PFO or ASD. Our analysis of world indigenous populations showed that GATA4 S377G was largely Caucasian-specific, and so subjects were restricted to those of Caucasian descent. To select for patients with larger PFO, we limited our analysis to those with cryptogenic stroke in which PFO was a subsequent finding. In an initial study of Australian subjects, we observed a weak association between GATA4 S377G and PFO/Stroke relative to Caucasian controls in whom ASD and PFO had been excluded (OR = 2.16; p = 0.02). However, in a follow up study of German Caucasians no association was found with either PFO or ASD. Analysis of combined Australian and German data confirmed the lack of a significant association. Thus, the common GATA4 variant S377G is likely to be relatively benign in terms of its participation in CHD and PFO/Stroke.

Common variants in the regulative regions of GRIA1 and GRIA3 receptor genes are associated with migraine susceptibility

Background Glutamate is the principal excitatory neurotransmitter in the central nervous system which acts by the activation of either ionotropic (AMPA, NMDA and kainate receptors) or G-protein coupled metabotropic receptors. Glutamate is widely accepted to play a major role in the path physiology of migraine as implicated by data from animal and human studies. Genes involved in synthesis, metabolism and regulation of both glutamate and its receptors could be, therefore, considered as potential candidates for causing/predisposing to migraine when mutated. Methods The association of polymorphic variants of GRIA1-GRIA4 genes which encode for the four subunits (GluR1-GluR4) of the alpha-amino-3- hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) receptor for glutamate was tested in migraineurs with and without aura (MA and MO) and healthy controls. Results Two variants in the regulative regions of GRIA1 (rs2195450) and GRIA3 (rs3761555) genes resulted strongly associated with MA (P = 0.00002 and P = 0.0001, respectively), but not associated with MO, suggesting their role in cortical spreading depression. Whereas the rs548294 variant in GRIA1 gene showed association primarily with MO phenotype, supporting the hypothesis that MA and MO phenotypes could be genetically related. These variants modify binding sites for transcription factors altering the expression of GRIA1 and GRIA3 genes in different conditions. Conclusions This study represents the first genetic evidence of a link between glutamate receptors and migraine.

Integrin αvβ3 upregulates integrin-linked kinase expression in human ovarian cancer cells via enhancement of ILK gene transcription

We previously showed that integrin alphavbeta3 overexpression and engagement by its ligand vitronectin increased adhesion, motility, and proliferation of human ovarian cancer cells. In search of differentially regulated genes involved in these tumor biological events, we previously identified the integrin-linked kinase (ILK) to be under control of alphavbeta3. In the present investigation we demonstrated significantly upregulated ILK protein as a function of alphavbeta3 in two ovarian cancer cell lines, OV-MZ-6 and OVCAR-3, and proved co-localization at the surface of alphavbeta3-overexpressing cells adherent to vitronectin. Increase of ILK protein was reflected by enhanced ILK promoter activity, an effect, which we further characterized with regard to transcriptional response elements involved. Abrogation of NF-kappaB/c-rel or p53 binding augmented ILK promoter activity and preserved induction by alphavbeta3. The AP1-mutant exhibited decreased promoter activity but was also still inducible by alphavbeta3. Disruption of the two DNA consensus motifs for Ets proteins led to divergent observations: mutation of the Ets motif at promoter position -462 bp did not significantly alter promoter activity but still allowed response to alphavbeta3. In contrast, disruption of the second Ets motif at position -85 bp did not only lead to slightly diminished promoter activity but also, in that case, abrogated ILK promoter induction by alphavbeta3. Subsequent co-transfection studies with ets-1 in the presence of the second Ets motif led to additional induction of ILK promoter activity. Taken together, these data suggest that ets-1 binding to the second Ets DNA motif strongly contributes to alphavbeta3-mediated ILK upregulation. By increasing ILK as an important integrin-proximal kinase, alphavbeta3 may promote its intracellular signaling and tumor biological processes arising thereof in favor of ovarian cancer metastasis.

The osteogenic transcription factor Runx2 regulates components of the fibroblast growth factor/proteoglycan signaling axis in osteoblasts

Heparan sulfate proteoglycans cooperate with basic fibroblast growth factor (bFGF/FGF2) signaling to control osteoblast growth and differentiation, as well as metabolic functions of osteoblasts. FGF2 signaling modulates the expression and activity of Runt-related transcription factor 2 (Runx2/Cbfa1), a key regulator of osteoblast proliferation and maturation. Here, we have characterized novel Runx2 target genes in osteoprogenitors under conditions that promote growth arrest while not yet permitting sustained phenotypic maturation. Runx2 enhances expression of genes related to proteoglycan-mediated signaling, including FGF receptors (e.g., FGFR2 and FGFR3) and proteoglycans (e.g., syndecans [Sdc1, Sdc2, Sdc3], glypicans [Gpc1], versican [Vcan]). Runx2 increases expression of the glycosyltransferase Exostosin-1 (Ext1) and heparanase, as well as alters the relative expression of N-linked sulfotransferases (Ndst1 = Ndst2 > Ndst3) and enzymes mediating O-linked sulfation of heparan sulfate (Hs2st > Hs6st) or chondroitin sulfate (Cs4st > Cs6st). Runx2 cooperates with FGF2 to induce expression of Sdc4 and the sulfatase Galns, but Runx2 and FGF2 suppress Gpc6, thus suggesting intricate Runx2 and FGF2 dependent changes in proteoglycan utilization. One functional consequence of Runx2 mediated modulations in proteoglycan-related gene expression is a change in the responsiveness of bone markers to FGF2 stimulation. Runx2 and FGF2 synergistically enhance osteopontin expression (>100 fold), while FGF2 blocks Runx2 induction of alkaline phosphatase. Our data suggest that Runx2 and the FGF/proteoglycan axis may form an extracellular matrix (ECM)-related regulatory feed-back loop that controls osteoblast proliferation and execution of the osteogenic program.