Candidate gene identification following linkage: Search for hypertension susceptibility genes

Following up genetic linkage studies to identify the underlying susceptibility gene(s) for complex disease traits is an arduous yet biologically and clinically important task. Complex traits, such as hypertension, are considered polygenic with many genes influencing risk, each with small effects. Chromosome 2 has been consistently identified as a genomic region with genetic linkage evidence suggesting that one or more loci contribute to blood pressure levels and hypertension status. Using combined positional candidate gene methods, the Family Blood Pressure Program has concentrated efforts in investigating this region of chromosome 2 in an effort to identify underlying candidate hypertension susceptibility gene(s). Initial informatics efforts identified the boundaries of the region and the known genes within it. A total of 82 polymorphic sites in eight positional candidate genes were genotyped in a large hypothesis-generating sample consisting of 1640 African Americans, 1339 whites, and 1616 Mexican Americans. To adjust for multiple comparisons, resampling-based false discovery adjustment was applied, extending traditional resampling methods to sibship samples. Following this adjustment for multiple comparisons, SLC4A5, a sodium bicarbonate transporter, was identified as a primary candidate gene for hypertension. Polymorphisms in SLC4A5 were subsequently genotyped and analyzed for validation in two populations of African Americans (N = 461; N = 778) and two of whites (N = 550; N = 967). Again, SNPs within SLC4A5 were significantly associated with blood pressure levels and hypertension status. While not identifying a single causal DNA sequence variation that is significantly associated with blood pressure levels and hypertension status across all samples, the results further implicate SLC4A5 as a candidate hypertension susceptibility gene, validating previous evidence for one or more genes on chromosome 2 that influence hypertension related phenotypes in the population-at-large. The methodology and results reported provide a case study of one approach for following up the results of genetic linkage analyses to identify genes influencing complex traits. ^

Comparative studies of Hox genes in mammalian limb and vertebral pattern formation

Divergence of anterior-posterior (AP) limb pattern and differences in vertebral column morphology are the two main examples of mammalian evolution. The Hox genes (homeobox containing gene) have been implicated in driving evolution of these structures. However, regarding Hox genes, how they contribute to the generation of mammalian morphological diversities, is still unclear. Implementing comparative gene expression and phenotypic rescue studies for different mammalian Hox genes could aid in unraveling this mystery. In the first part of this thesis, the expression pattern of Hoxd13 gene, a key Hox gene in the establishment of the limb AP pattern, was examined in developing limbs of bats and mice. Bat forelimbs exhibit a pronounced asymmetric AP pattern and offer a good model to study the molecular mechanisms that contribute to the variety of mammalian limbs. The data showed that the expression domain of bat Hoxd13 was shifted prior to the asymmetric limb plate expansion, whereas its domain in mice was much more symmetric. This finding reveals a correlation between the divergence of Hoxd13 expression and the AP patterning difference in limb development. The second part of this thesis details a phenotypic rescue approach by human HOXB1-9 transgenes in mice with Hoxb1-9 deletion, The mouse mutants displayed homeosis in cervical and anterior thoracic vertebrae. The human transgenes entirely rescued the mouse mutants, suggesting that these human HOX genes have similar functions to their mouse orthologues in anterior axial skeletal patterning. The anterior expressing human HOXB transgenes such as HOXB1-3 were expressed in the mouse embryonic trunk in a similar manner as their murine orthologues. However, the anterior boundary of human HOXB9 expression domain was more posterior than that of the mouse Hoxb9 by 2-3 somites. These data provide the molecular support for the hypothesis that Hox genes are responsible for maintaining similar anterior axial skeletal architectures cervical and anterior thoracic regions, but different architectures in lumbar and posterior thoracic regions between humans and mice. ^

Insulin-like growth factor binding protein 2 (IGFBP2) promotes glioma development and progression

Overexpression of insulin-like growth factor binding protein 2 (IGFBP2) is associated with progression and poor survival in many types of human cancer (such as prostate, ovarian, adrenocortical, breast, colorectal carcinomas, leukemia, and high-grade gliomas). We therefore hypothesize that IGFBP2 is a key regulator of tumor progression. We tested our hypothesis in gliomas using the somatic gene transfer RCAS-tva mouse model system, which permits the introduction of specific genes into specific, cell lineages, in this case glial cells (RCAS: Replication competent avian sarcomavirus, tv-a: avian RCAS virus receptor). Mice are transgenic and harbor the tv-a receptor under the control of a glial-specific promoter and study genes are cloned into the RCAS vector for post-natal intracranial delivery. For these experiments, the study genes were IGFBP2, platelet-derived growth factor B (PDGFB), K-Ras, Akt, and IIp45 (invasion inhibitory protein 45 kDa; known to bind and block IGFBP2 activity), which were delivered separately and in combination. Our results show that PDGFB signaling leads exclusively to the formation of low-grade (WHO grade II) oligodendrogliomas. PDGFB delivered in combination with IGFBP2 results in the formation of anaplastic oligodendrogliomas (WHO grade III), which are characterized by increased cellularity, vascular proliferation, small regions of necrosis, increased mitotic activity, and increased activation of the Akt pathway. IIp45 injected in combination with PDGFB and IGFBP2 ablates IGFBP2-induced tumor progression, which results in formation of low-grade oligodendrogliomas, and an overall reduction in tumor incidence. K-Ras expression was required to form astrocytomas with either IGFBP2 or Akt, indicating the activation of two separate pathways is necessary for gliomagenesis. In ex vivo experiments, blockade of Akt by an inhibitor led to decreased viability of cells co-expressing IGFBP2 versus PDGFB expression alone. This study provides definitive evidence, for the first time, that: (1) IGFBP2 plays a role in activation of the Akt pathway, (2) IGFBP2 collaborates with K-Ras or PDGFB in the development and progression of two major types of glioma, and (3) IGFBP2-induced tumor progression can be ablated by IIp45 or by specific inhibition of the Akt pathway. ^

Diabetes genes and risk of prostate cancer in the Atherosclerosis Risk in Communities study

Prostate cancer (PrCa) is a leading cause of morbidity and mortality, yet the etiology remains uncertain. Meta-analyses show that PrCa risk is reduced by 16% in men with type 2 diabetes (T2D), but the mechanism is unknown. Recent genome-wide association studies and meta-analyses have found single nucleotide polymorphisms (SNPs) that consistently predict T2D risk. We evaluated associations of incident PrCa with 14 T2D SNPs in the Atherosclerosis Risk in Communities (ARIC) study. From 1987-2000, there were 397 incident PrCa cases ascertained from state or local cancer registries among 6,642 men (1,560 blacks and 5,082 whites) aged 45-64 years at baseline. Genotypes were determined by TaqMan assay. Cox proportional hazards models were used to assess the association between PrCa and increasing number of T2D risk-raising alleles for individual SNPs and for genetic risk scores (GRS) comprised of the number of T2D risk-raising alleles across SNPs. Two-way gene-gene interactions were evaluated with likelihood ratio tests. Using additive genetic models, the T2D risk-raising allele was associated with significantly reduced risk of PrCa for IGF2BP2 rs4402960 (hazard ratio [HR]=0.79; P=0.07 among blacks only), SLC2A2 rs5400 (race-adjusted HR=0.85; P=0.05) and UCP2 rs660339 (race-adjusted HR=0.84; P=0.02), but significantly increased risk of PrCa for CAPN10 rs3792267 (race-adjusted HR=1.20; P=0.05). No other SNPs were associated with PrCa using an additive genetic model. However, at least one copy of the T2D risk-raising allele for TCF7L2 rs7903146 was associated with reduced PrCa risk using a dominant genetic model (race-adjusted HR=0.79; P=0.03). These results imply that the T2D-PrCa association may be partly due to shared genetic variation, but these results should be verified since multiple tests were performed. When the combined, additive effects of these SNPs were tested using a GRS, there was nearly a 10% reduction in risk of PrCa per T2D risk-raising allele (race-adjusted HR=0.92; P=0.02). SNPs in IGF2BP2, KCNJ11 and SLC2A2 were also involved in multiple synergistic gene-gene interactions on a multiplicative scale. In conclusion, it appears that the T2D-PrCa association may be due, in part, to common genetic variation. Further knowledge of T2D gene-PrCa mechanisms may improve understanding of PrCa etiology and may inform PrCa prevention and treatment.^

Identification of diabetic retinopathy genes through a genome-wide association study among Mexican-Americans from Starr County, Texas

To identify genetic susceptibility loci for severe diabetic retinopathy, 286 Mexican-Americans with type 2 diabetes from Starr County, Texas completed detailed physical and ophthalmologic examinations including fundus photography for diabetic retinopathy grading. 103 individuals with moderate-to-severe non-proliferative diabetic retinopathy or proliferative diabetic retinopathy were defined as cases for this study. DNA samples extracted from study subjects were genotyped using the Affymetrix GeneChip® Human Mapping 100K Set, which includes 116,204 single nucleotide polymorphisms (SNPs) across the whole genome. Single-marker allelic tests and 2- to 8-SNP sliding-window Haplotype Trend Regression implemented in HelixTreeTM were first performed with these direct genotypes to identify genes/regions contributing to the risk of severe diabetic retinopathy. An additional 1,885,781 HapMap Phase II SNPs were imputed from the direct genotypes to expand the genomic coverage for a more detailed exploration of genetic susceptibility to diabetic retinopathy. The average estimated allelic dosage and imputed genotypes with the highest posterior probabilities were subsequently analyzed for associations using logistic regression and Fisher's Exact allelic tests, respectively. To move beyond these SNP-based approaches, 104,572 directly genotyped and 333,375 well-imputed SNPs were used to construct genetic distance matrices based on 262 retinopathy candidate genes and their 112 related biological pathways. Multivariate distance matrix regression was then used to test hypotheses with genes and pathways as the units of inference in the context of susceptibility to diabetic retinopathy. This study provides a framework for genome-wide association analyses, and implicated several genes involved in the regulation of oxidative stress, inflammatory processes, histidine metabolism, and pancreatic cancer pathways associated with severe diabetic retinopathy. Many of these loci have not previously been implicated in either diabetic retinopathy or diabetes. In summary, CDC73, IL12RB2, and SULF1 had the best evidence as candidates to influence diabetic retinopathy, possibly through novel biological mechanisms related to VEGF-mediated signaling pathway or inflammatory processes. While this study uncovered some genes for diabetic retinopathy, a comprehensive picture of the genetic architecture of diabetic retinopathy has not yet been achieved. Once fully understood, the genetics and biology of diabetic retinopathy will contribute to better strategies for diagnosis, treatment and prevention of this disease.^

Identification and characterization of mutations in SMC contractile genes involved in thoracic aortic aneurysms and dissections

Aortic aneurysms and dissections are the 15th most common cause of death in the United States. Genetic factors contribute to the pathogenesis of thoracic aortic aneurysms and dissections (TAAD). Currently, six loci and four genes have been identified for familial TAAD. Notably, mutations in smooth muscle cell (SMC) contractile genes, ACTA2 and MYH11, are responsible for 15% of familial TAAD, suggesting that proper SMC contraction is important for normal aorta function. Therefore, we hypothesize that mutations in other genes encoding SMC contractile proteins also cause familial TAAD. ^ To test this hypothesis, we used a candidate gene approach to identify causative mutations in SMC contractile genes for familial TAAD. Sequencing DNA in 80 TAAD patients from unrelated families, we identified putative mutations in eight contractile genes. We chose myosin light chain kinase (MLCK ) S1759P for further study for the following reasons: (1) Serine 1759 is conserved between vertebrates and invertebrates. (2) S1759P is predicted to be functionally deleterious by bioinformatics. (3) Low blood pressure is observed in SMC-selective MLCK-deficient mice. ^ In the presence of Ca2+/Calmodulin (CaM), MLCK containing CaM binding and kinase domains are activated to phosphorylate myosin light chain, thereby initiate SMC contraction. The CaM binding sequence of MLCK forms an α-helix structure required for CaM binding. MLCK Serine 1759 is located within the CaM binding domain. S1759P is predicted to decrease the α-helix composition in the CaM binding domain. Hence, we hypothesize that MLCK mutations cause TAAD through disturbing CaM binding and MLCK activity. ^ We further sequenced MLCK in DNA samples from additional 86 probands with familial TAAD. Two more mutations, MLCK A1754T and R1480Stop, were identified, supporting that MLCK mutations cause familial TAAD. ^ To define whether MLCK mutations disrupted CaM binding and MLCK activity, we performed co-immunoprecipitation and kinase assays. Decreased CaM binding and kinase activity was detected in A1754T and S1759P. Moreover, R1480Stop is predicted to truncate kinase and CaM binding domains. We conclude that MLCK mutations disrupt CaM binding and MLCK activity. ^ Collectively, our study is first to show mutations in genes regulating SMC contraction cause TAAD. This finding further highlights the importance of SMC contraction in maintaining aorta function. ^

Studies of the expression of DNA repair related genes

This dissertation examines the biological functions and the regulation of expression of DNA ligase I by studying its expression under different conditions.^ The gene expression of DNA ligase I was induced two- to four-fold in S-phase lymphoblastoid cells but was decreased to 15% of control after administration of a DNA damaging agent, 4-nitroquinoline-1-oxide. When cells were induced into differentiation, the expression level of DNA ligase I was decreased to less than 15% of that of the control cells. When the gene of DNA ligase I was examined for tissue specific expression in adult rats, high levels of DNA ligase I mRNA were observed in testis (8-fold), intermediate levels in ovary and brain (4-fold), and low levels were found in intestine, spleen, and liver (1- to 2-fold).^ In confluent cells of normal skin fibroblasts, UV irradiation induced the gene expression of DNA ligase I at 24 and 48 h. The induction of DNA ligase I gene expression requires active p53 protein. Introducing a vector containing the wild type p53 protein in the cells caused an induction of the DNA ligase I protein 24 h after the treatment.^ Our results indicate that, in addition to the regulation by phosphorylation/dephosphorylation, cellular DNA ligase I activity can be regulated at the gene transcription level, and the p53 tumor suppresser is one of the transcription factors for the DNA ligase I gene. Also, our results suggest that DNA ligase I is involved in DNA repair as well as in DNA replication.^ Also, as an early attempt to clone the human homolog of the yeast CDC9 gene which has been shown to be involved in DNA replication, DNA repair, and DNA recombination, we have identified a human gene with mRNA of 1.7 kb. This dissertation studies the gene regulation and the possible biological functions of this new human gene by examining its expression at different stages of the cell cycle, during cell differentiation, and in cellular response to DNA damage.^ The new gene that we recently identified from human cells is highly expressed in brain and reproductive organs (BRE). This BRE gene encodes an mRNA of 1.7-1.9 kb, with an open reading frame of 1,149 bp, and gives rise to a deduced polypeptide of 383 amino acid residues. No extensive homology was found between BRE and sequences from the EMBL-Gene Banks. BRE showed tissue-specific expression in adult rats. The steady state mRNA levels were high in testis (5-6 fold), ovary and brain (3-4 fold) compared to the spleen level, but low in intestine and liver (1-2 fold). The expression of this gene is responsive to DNA damage and/or retinoic acid (RA) treatment. Treatment of fibroblast cells with UV irradiation and 4-nitroquinoline-1-oxide caused more than 90% and 50% decreases in BRE mRNA, respectively. Similar decreases in BRE expression were observed after treatment of the brain glioma cell line U-251 and the promyelocytic cell line HL-60 with retinoic acid. (Abstract shortened by UMI). ^

Longitudinal and logistic analysis of endocrine hormone profiles, physical maturation, and candidate genes related to pubertal events in children

Children who experience early pubertal development have an increased risk of developing cancer (breast, ovarian, and testicular), osteoporosis, insulin resistance, and obesity as adults. Early pubertal development has been associated with depression, aggressiveness, and increased sexual prowess. Possible explanations for the decline in age of pubertal onset include genetics, exposure to environmental toxins, better nutrition, and a reduction in childhood infections. In this study we (1) evaluated the association between 415 single nucleotide polymorphisms (SNPs) from hormonal pathways and early puberty, defined as menarche prior to age 12 in females and Tanner Stage 2 development prior to age 11 in males, and (2) measured endocrine hormone trajectories (estradiol, testosterone, and DHEAS) in relation to age, race, and Tanner Stage in a cohort of children from Project HeartBeat! At the end of the 4-year study, 193 females had onset of menarche and 121 males had pubertal staging at age 11. African American females had a younger mean age at menarche than Non-Hispanic White females. African American females and males had a lower mean age at each pubertal stage (1-5) than Non-Hispanic White females and males. African American females had higher mean BMI measures at each pubertal stage than Non-Hispanic White females. Of the 415 SNPs evaluated in females, 22 SNPs were associated with early menarche, when adjusted for race ( p<0.05), but none remained significant after adjusting for multiple testing by False Discovery Rate (p<0.00017). In males, 17 SNPs were associated with early pubertal development when adjusted for race (p<0.05), but none remained significant when adjusted for multiple testing (p<0.00017). ^ There were 4955 hormone measurements taken during the 4-year study period from 632 African American and Non-Hispanic White males and females. On average, African American females started and ended the pubertal process at a younger age than Non-Hispanic White females. The mean age of Tanner Stage 2 breast development in African American and Non-Hispanic White females was 9.7 (S.D.=0.8) and 10.2 (S.D.=1.1) years, respectively. There was a significant difference by race in mean age for each pubertal stage, except Tanner Stage 1 for pubic hair development. Both Estradiol and DHEAS levels in females varied significantly with age, but not by race. Estradiol and DHEAS levels increased from Tanner Stage 1 to Tanner Stage 5.^ African American males had a lower mean age at each Tanner Stage of development than Non-Hispanic White males. The mean age of Tanner Stage 2 genital development in African American and Non-Hispanic White males was 10.5 (S.D.=1.1) and 10.8 (S.D.=1.1) years, respectively, but this difference was not significant (p=0.11). Testosterone levels varied significantly with age and race. Non-Hispanic White males had higher levels of testosterone than African American males from Tanner Stage 1-4. Testosterone levels increased for both races from Tanner Stage 1 to Tanner Stage 5. Testosterone levels had the steepest increase from ages 11-15 for both races. DHEAS levels in males varied significantly with age, but not by race. DHEAS levels had the steepest increase from ages 14-17. ^ In conclusion, African American males and females experience pubertal onset at a younger age than Non-Hispanic White males and females, but in this study, we could not find a specific gene that explained the observed variation in age of pubertal onset. Future studies with larger study populations may provide a better understanding of the contribution of genes in early pubertal onset.^

FZD6, MATN2 AND SLC25A32, POSSIBLE CANDIDATE GENES IN NONSYNDROMIC CLEFT LIP AND PALATE

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect with a multifactorial etiology. Despite decades of research, the genetic underpinnings of NSCLP still remain largely unexplained. A genome wide association study (GWAS) of a large NSCLP African American family with seven affected individuals across three generations found evidence for linkage at 8q21.3-24.12 (LOD = 2.98). This region contained three biologically relevant candidate genes: Frizzled-6 (FZD6) (LOD = 2.8), Matrilin-2 (MATN2) (LOD = 2.3), and Solute Carrier Family 25, Member 32 (SLC26A32) (LOD = 1.6). Sequencing of the coding regions and the 5’ and 3’ UTRs of these genes in two affected family members identified a rare intronic variant, rs138557689 (c.-153+432A>C), in FZD6. The rs138557689/C allele segregated with the NSCLP phenotype; in silico analysis predicted and EMSA analysis showed that the 138557689/C allele creates new DNA binding sites. FZD6 is part of the WNT pathway, which is involved in craniofacial development, including midface development and upper lip fusion. Our novel findings suggest that an alteration in FZD6 gene regulation may perturb this tightly controlled biological pathway and in turn contribute to the development of NSCLP in this family. Studies are underway to further define how the rs138557689/C variant affects expression of FZD6.

2-methylthio-N$\sp6$-dimethylallyl modification of adenosine 37 of tRNAs in {\it Escherichia coli\/} K-12

The hypermodified, hydrophobic 2-methylthio-N$\sp6$-(dimethylallyl)-adenosine (ms${2{\cdot}6}\atop1$A) residue occurs $3\sp\prime$ to the anticodon in tRNA species that read codons beginning with U. The first step (i$\sp6$A37 formation) of this modification is catalyzed by dimethylallyl diphosphate:tRNA dimethyallyltransferase (EC 2.5.1.8), which is the product of the miaA gene. Subsequent steps were proposed to be catalyzed by MiaB and MiaC enzymes to complete the ms${2{\cdot}6}\atop1$A37 modification. The study of functions of the ms${2{\cdot}6}\atop1$A37 is very important because this modified base is one of the best candidates for a role in global control in response to environmental stress. This dissertation describes the further delineation of functions of the ms${2{\cdot}6}\atop1$A37 modification in E. coli K-12 cells. This work provides significant information on functions of tRNA modifications in E. coli cells to adapt to stressful environmental conditions. Three hypotheses were tested in this work.^ The first hypothesis tested was that non-optimal translation processes cause increased spontaneous mutagenesis by the induction of SOS response in starving cells. To test this hypothesis, I measured spontaneous mutation rates of wild type cells and various mutant strains which are defective in tRNA modification, SOS response, or oxidative damage repair. I found that the miaA mutation acts as a mutator that increased Lac$\sp+$ reversion rates and Trp$\sp+$ reversion frequencies of the wild-type cells in starving conditions. However, the lexA3(Ind)(which abolishes the induction of SOS response) mutation abolished the mutator phenotype of the miaA mutant. The recA430 mutation, not other identified SOS genes, decreased the Lac$\sp+$ reversion to a less extent than that of the lexA3(Ind) mutation. These results suggest that RecA together with another unidentified SOS gene product are responsible for the process.^ The second hypothesis tested was that MiaA protein binds to full-length tRNA$\sp{\rm Phe}$ molecules in form of a protein dimer. To test this hypothesis, three versions of the MiaA protein and seven species of tRNA substrates were purified. Binding studies by gel mobility shift assays, filter binding assays and gel filtration shift assays support the hypothesis that MiaA protein binds to full-length tRNA$\sp{\rm Phe}$ as a protein dimer but as a monomer to the anticodon stem-and-loop. These results were further supported by using steady state enzyme kinetic studies.^ The third hypothesis tested in this work was that the miaB gene in E. coli exists and is clonable. The miaB::Tn10dCm insertion mutation of Salmonella typhimurium was transduced to E. coli K-12 cells by using P$\sb1$ and P$\sb{22}$ bacteriophages. The insertion was confirmed by HPLC analyses of nucleotide profiles of miaB mutants of E. coli. The insertion mutation was cloned and DNA sequences adjacent to the transposon were sequenced. These DNA sequences were 86% identical to the f474 gene at 14.97 min chromosome of E. coli. The f474 gene was then cloned by PCR from the wild-type chromosome of E. coli. The recombinant plasmid complemented the mutant phenotype of the miaB mutant of E. coli. These results support the hypothesis that the miaB gene of E. coli exists and is clonable. In summary, functions of the ms${2{\cdot}6}\atop1$A37 modification in E. coli cells are further delineated in this work in perspectives of adaptation to stressful environmental conditions and protein:tRNA interaction. (Abstract shortened by UMI.) ^

Transcriptional regulation of cell lineage-specific expression of the murine type I collagen genes and of osteoblast differentiation

The aim of my project is to examine the mechanisms of cell lineage-specific transcriptional regulation of the two type I collagen genes by characterizing critical cis-acting elements and trans-acting factors. I hypothesize that the transcription factors that are involved in the cell lineage-specific expression of these genes may have a larger essential role in cell lineage commitment and differentiation. I first examined the proximal promoters of the proα1(I) and the proα2(I) collagen genes for cell type-specific DNA-protein interactions, using in vitro DNaseI and in vivo DMS footprinting. These experiments demonstrated that the cis-acting elements in these promoters are accessible to ubiquitous DNA-binding proteins in fibroblasts that express these genes, but not in other cells that do not express these genes. I speculate that in type I collagen-expressing cells, cell type-specific enhancer elements facilitate binding of ubiquitous proteins to the proximal promoters of these genes. Subsequently, examination of the upstream promoter of the proα(I) collagen gene by transgenic mice experiments delineated a 117 bp sequence (-1656 to -1540 bp) as the minimum element required for osteoblast-specific expression. This 117 bp element contained two segments that appeared to have different functions: (1) the A-segment, which was necessary to obtain osteoblast-specific expression and (2) the C-segment, which was dispensable for osteoblast-specific expression, but was necessary to obtain high-level expression. In experiments to identify trans-acting factors that bind to the 117 bp element, I have demonstrated that the cell lineage-restricted homeodomain proteins, Dlx2, Dlx5 and mHOX, bound to the A-segment and that the ubiquitous transcription factor, Sp1, bound to the C-segment of this element. These results suggested a model where the binding of cell lineage-restricted proteins to the A-segment and of ubiquitous proteins to the C-segment of the 117 bp element of the proα1 (I) collagen gene activated this gene in osteoblasts. These results, combined with additional evidence that Dlx2, Dlx5 and mHOX are probably involved in osteoblast differentiation, support my hypothesis that the transcription factors involved in osteoblast-specific expression of type I collagen genes may have essential role in osteoblast lineage commitment and differentiation. ^