Prediction of regulatory elements

Expression Profiling and Functional Analysis of Candidate Col10a1 Regulators Identified by the TRAP Program

Hypertrophic chondrocytes and their specific marker, the type X collagen gene (Col10a1), are critical components of endochondral bone formation during skeletal development. We previously found that Runx2 is an indispensable mouse Col10a1 gene regulator and identified many other transcription factors (TFs) that potentially interact with the 150-bp Col10a1 cis-enhancer. However, the roles of these candidate TFs in Col10a1 expression and chondrocyte hypertrophy have not been elucidated. Here, we focus on 32 candidate TFs recently identified by analyzing the 150-bp Col10a1 enhancer using the transcription factor affinity prediction (TRAP) program. We found that 12 TFs (Hoxa3, Lsx, Evx2, Dlx5, S8, Pax2, Egr2, Mef2a, Barhl2, GKlf, Sox17, and Crx) were significantly upregulated and four TFs (Lhx4, Tbx5, Mef2c, and Hb9) were significantly downregulated in hypertrophic MCT cells, which show upregulation of Col10a1 expression. Most of the differential expression pattern of these TFs conformed with the results obtained from ATDC5 cell model and primary mouse chondrocytes. Notably, Tbx5 was downregulated upon Col10a1 upregulation, overexpression of Tbx5 decreased Col10a1 expression, and knock-down of Tbx5 increased Col10a1 expression in hypertrophic chondrocytes, suggesting that Tbx5 is a negative regulator of Col10a1. We further generated a stable Tbx5-overexpressing ATDC5 cell line and ColX-Tbx5 transgenic mice driven by Col10a1-specific enhancers and promoters. Tbx5 overexpression decreased Col10a1 expression in ATDC5 cells cultured as early as day 7 and in limb tissue on post-natal day 1. Slightly weaker alkaline phosphatase staining was also observed in cell culture on day 7 and in limb digits on embryonic day 17.5, indicating mildly delayed ossification. Further characterization of these candidate Col10a1 transcriptional regulators could help identify novel therapeutic targets for skeletal diseases associated with abnormal chondrocyte hypertrophy.

Identification and characterization of the novel Col10a1 regulatory mechanism during chondrocyte hypertrophic differentiation

The majority of human skeleton develops through the endochondral pathway, in which cartilage-forming chondrocytes proliferate and enlarge into hypertrophic chondrocytes that eventually undergo apoptosis and are replaced by bone. Although at a terminal differentiation stage, hypertrophic chondrocytes have been implicated as the principal engine of bone growth. Abnormal chondrocyte hypertrophy has been seen in many skeletal dysplasia and osteoarthritis. Meanwhile, as a specific marker of hypertrophic chondrocytes, the type X collagen gene (COL10A1) is also critical for endochondral bone formation, as mutation and altered COL10A1 expression are often accompanied by abnormal chondrocyte hypertrophy in many skeletal diseases. However, how the type X collagen gene is regulated during chondrocyte hypertrophy has not been fully elucidated. We have recently demonstrated that Runx2 interaction with a 150-bp mouse Col10a1 cis-enhancer is required but not sufficient for its hypertrophic chondrocyte-specific reporter expression in transgenic mice, suggesting requirement of additional Col10a1 regulators. In this study, we report in silico sequence analysis of this 150-bp enhancer and identification of its multiple binding factors, including AP1, MEF2, NFAT, Runx1 and TBX5. Using this enhancer as bait, we performed yeast one-hybrid assay and identified multiple candidate Col10a1-interacting genes, including cyclooxygenase 1 (Cox-1) and Cox-2. We have also performed mass spectrometry analysis and detected EF1-alpha, Fus, GdF7 and Runx3 as components of the specific complex formed by the cis-enhancer and nuclear extracts from hypertrophic MCT (mouse chondrocytes immortalized with large T antigen) cells that express Col10a1 abundantly. Notably, some of the candidate genes are differentially expressed in hypertrophic MCT cells and have been associated with chondrocyte hypertrophy and Runx2, an indispensible Col10a1 regulator. Intriguingly, we detected high-level Cox-2 expression in hypertrophic MCT cells. Electrophoretic mobility shift assay and chromatin immunoprecipitation assays confirmed the interaction between Cox-2 and Col10a1 cis-enhancer, supporting its role as a candidate Col10a1 regulator. Together, our data support a Cox-2-containing, Runx2-centered Col10a1 regulatory mechanism, during chondrocyte hypertrophic differentiation.

Divergence of the SigB regulon and pathogenesis of the Bacillus cereus sensu lato group

Background

The Bacillus cereus sensu lato group currently includes seven species (B. cereus, B. anthracis, B. mycoides, B. pseudomycoides, B. thuringiensis, B. weihenstephanensis and B. cytotoxicus) that recent phylogenetic and phylogenomic analyses suggest are likely a single species, despite their varied phenotypes. Although horizontal gene transfer and insertion-deletion events are clearly important for promoting divergence among these genomes, recent studies have demonstrated that a major basis for phenotypic diversity in these organisms may be differential regulation of the highly similar gene content shared by these organisms. To explore this hypothesis, we used an in silico approach to evaluate the relationship of pathogenic potential and the divergence of the SigB-dependent general stress response within the B. cereus sensu lato group, since SigB has been demonstrated to support pathogenesis in Bacillus, Listeria and Staphylococcus species.

Results

During the divergence of these organisms from a common “SigB-less” ancestor, the placement of SigB promoters at varied locations in the B. cereus sensu lato genomes predict alternative structures for the SigB regulon in different organisms. Predicted promoter changes suggesting differential transcriptional control of a common gene pool predominate over evidence of indels or horizontal gene transfer for explaining SigB regulon divergence.

Conclusions

Four lineages of the SigB regulon have arisen that encompass different gene contents and suggest different strategies for supporting pathogenesis. This is consistent with the hypothesis that divergence within the B. cereus sensu lato group rests in part on alternative strategies for regulation of a common gene pool.

Identifying functional single nucleotide polymorphisms in the human CArGome

Regulatory SNPs (rSNPs) reside primarily within the nonprotein coding genome and are thought to disturb normal patterns of gene expression by altering DNA binding of transcription factors. Nevertheless, despite the explosive rise in SNP association studies, there is little information as to the function of rSNPs in human disease. Serum response factor (SRF) is a widely expressed DNA-binding transcription factor that has variable affinity to at least 1,216 permutations of a 10 bp transcription factor binding site (TFBS) known as the CArG box. We developed a robust in silico bioinformatics screening method to evaluate sequences around RefSeq genes for conserved CArG boxes. Utilizing a predetermined phastCons threshold score, we identified 8,252 strand-specific CArGs within an 8 kb window around the transcription start site of 5,213 genes, including all previously defined SRF target genes. We then interrogated this CArG dataset for the presence of previously annotated common polymorphisms. We found a total of 118 unique CArG boxes harboring a SNP within the 10 bp CArG sequence and 1,130 CArG boxes with SNPs located just outside the CArG element. Gel shift and luciferase reporter assays validated SRF binding and functional activity of several new CArG boxes. Importantly, SNPs within or just outside the CArG box often resulted in altered SRF binding and activity. Collectively, these findings demonstrate a powerful approach to computationally define rSNPs in the human CArGome and provide a foundation for similar analyses of other TFBS. Such information may find utility in genetic association studies of human disease where little insight is known regarding the functionality of rSNPs.

Exploring the landscape of the genome

Genome browsers are powerful tools for biologists - offering fundamental information on genes, regulatory elements, genomic variants, genome structure, and evolution. The comprehensive range of information presented in tools such as the UCSC genome browser and Ensembl enables integrated queries of data that are otherwise reserved to the most skilled computational biologists. However, for the non-specialist user, the juxtaposition of so many different forms of data in one small space can be an information overload. Getting the most out of these tools requires some understanding of the key concepts and caveats of genome visualization and annotation. Genome analysis can be carried out at different levels of detail - at a macro level; it improves understanding of issues like genome structure and species evolution. While at a micro level, genome annotation can help to describe the full complexity of gene regulation, variation, and transcript diversity. Once demystified, it is clear that genome browsers are more than the sum of their parts - they are the most comprehensive portals available for browsing and analysis of biological data.