Analysis of BMP signalling through the receptor type IA in embryogenesis using conditional gene inactivation in mice

Although bone morphogenetic proteins (BMPs) were initially identified for their potent bone-inducing activity, their precise roles in processes of endochondral and intramembranous bone formation are far from being clear. Tissue-specific loss-of-function experiments using the BMP receptor type IA (BMPR-IA) are particularly attractive since this receptor is thought to be essential for signaling by the closely related BMPs -2, 4, and 7. To ablate signaling through this receptor during chondrogenesis, we have generated transgenic mice expressing Cre recombinase under the control of the collagen type II (Col2a1) gene regulatory sequences. Mice lacking BMPR-IA function in chondrocytes display a number of skeletal abnormalities, including defects in bones of the chondrocranium, abnormal dorsal vertebral processes, scapulae with severe hypoplasia of dorsal elements, and shortening of the long bones. Alterations in the growth plate of long bones in mutants suggest that BMPR-IA is not required for early steps of the chondrocyte specification, but is rather important in regulation of terminal differentiation. Molecular analysis revealed noticeable downregulation of the Ihh/Ptch signalling pathway, decreased chondrocyte proliferation rate and deregulation of hypertrophy. ^ In order to elucidate the role of BMP signalling in development of the limb and intramembranous ossification, we have used mice expressing Cre recombinase under control of the Prx1 (MHox) regulatory elements (M. Logan, pers comm.). Cre activity was found in those mice in the developing limb bud mesenchyme, as well as in a subset of cranial neural crest cells. Prx1-Cre-induced conditional mutants display prominent defects in distal limb outgrowth, as well as ossification defects in a number of neural crest-derived calvarial bones. Intriguingly, mutant limbs displayed alterations in patterning along all three axes. Molecular analysis revealed ectopic anterior Shh/Ptch signalling pathway activation and expression of some Hox genes. Observed loss of Msx1 and Msx2 expression in the progress zone correlates with downregulation of Cyclin D1 and decreased distal outgrowth. Abnormal ventral localization of Lmx1b-expressing cells along with observed later morphological abnormalities suggest a novel role for BMP signalling in establishment or maintaining of the dorso-ventral polarity in the limb mesoderm. ^

Function and regulation of the Pem homeobox gene in vivo

Pem, a member of the PEPP homeobox family, is expressed in somatic cells in male and female reproductive tissues. In the adult murine testis, Pem is specifically expressed in Sertoli cells, where it is restricted to stages IV–VIII of the seminiferous epithelial cycle. To identify Pem's function in Sertoli cells, transgenic mice were generated that express Pem in Sertoli cells during all stages of the seminiferous epithelial cycle. This resulted in an increase in double-strand DNA breaks in preleptotene spermatocytes and single-strand DNA breaks in elongating spermatids. My results suggest that Pem regulates Sertoli-cell genes that encode secreted or cell-surface proteins that serve to control premeiotic DNA replication, DNA repair, and/or chromatin remodeling in the adjacent germ cells. Three additional transgenic mouse containing varying lengths of the Pem male-specific promoter (Pp) were generated to identify the sequences responsible for regulating Pem expression in the testis and epididymis. My analysis suggests that there are at least two regulatory regions in the Pem Pp. In the testis, region II directs androgen-dependent expression specifically in Sertoli cells whereas region I fine-tunes stage-specific expression by acting as a negative regulator. In the epididymis, region II confers androgen-dependent, developmentally-regulated expression in the caput whereas region I prevents inappropriate expression in the corpus. I also report the identification and characterization of two human PEPP family members related to Pem that I have named hPEPP1 and hPEPP2. The hPEPP1 and hPEPP2 homeodomains are more closely related to PEPP subfamily homeodomains than to any other homeodomain subfamily. Both genes are localized to the specific region of the human X chromosome that shares synteny with the region on the murine X chromosome containing three PEPP homeobox genes, Pem, Psx-1, and Psx-2. hPEPP1 and hPEPP2 mRNA expression is restricted to the testis but is aberrantly expressed in tumor cells of different origins, analogous to the expression pattern of Pem but not of Psx-1 or Psx-2. Unlike all known PEPP members, neither hPEPP1 nor hPEPP2 are expressed in placenta, which suggests that the regulation of the PEPP family has undergone significant alteration since the split between hominids and rodents. ^

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. ^

Genetic predictors of bladder cancer with an emphasis on methylation-related genes

Although tobacco exposure remains the prevailing risk factor for bladder cancer (BC), only a small percentage of exposed individuals develop cancer, suggesting that tobacco-related carcinogenesis is modulated by genetic susceptibility and possibly by DNA methylation-related events. Methylation patterns established by DNA methyltransferases (DNMTs) are influenced by dietary folate and genetic polymorphisms in the methylene-tetrahydrofolate reductase gene (MTHFR). Therefore, we hypothesized that DNA methylation-related genes, such as DNMT3B and MTHFR, might modulate BC risk. ^ In a study of 514 Caucasian BC cases and 498 healthy Caucasian controls examining the DNMT3B C46359T polymorphism, CC genotype was found to be a risk factor in women (Odds Ratio (OR) = 1.79), but not in men. This risk was further increased among women who were never smokers, consumed low dietary folate, and had adverse variants of MTHFR. In addition, higher DNMT3B expression among smokers was a risk factor (OR = 4.27) and correlated with genetic variants of the DNMT3B C46359T polymorphism, providing salient evidence for the risk associated with the CC variant. This suggests that the DNMT3B CC variant may confer a predisposition toward aberrant de novo methylation of CpG islands in critical tumor suppressor genes. ^ The convergence of alterations in DNMT3B, associated with promoter methylation, and reduced dietary folate consumption, accompanying global hypomethylation and genetic instability, may act synergistically to promote bladder carcinogenesis, especially in women. The results of this study unveiled new gender-specific paradigms of BC risk for women and demonstrated that this risk can be modified by folate consumption as well as polymorphisms in the folate pathway. ^

DNA incisions stimulate genetic instability of triplet repeat sequences related to human hereditary neurological diseases

The molecular mechanisms responsible for the expansion and deletion of trinucleotide repeat sequences (TRS) are the focus of our studies. Several hereditary neurological diseases including Huntington's disease, myotonic dystrophy, and fragile X syndrome are associated with the instability of TRS. Using the well defined and controllable model system of Escherichia coli, the influences of three types of DNA incisions on genetic instability of CTG•CAG repeats were studied: DNA double-strand breaks (DSB), single-strand nicks, and single-strand gaps. The DNA incisions were generated in pUC19 derivatives by in vitro cleavage with restriction endonucleases. The cleaved DNA was then transformed into E. coli parental and mutant strains. Double-strand breaks induced deletions throughout the TRS region in an orientation dependent manner relative to the origin of replication. The extent of instability was enhanced by the repeat length and sequence (CTG•CAG vs. CGG•CCG). Mutations in recA and recBC increased deletions, mutations in recF stabilized the TRS, whereas mutations in ruvA had no effect. DSB were repaired by intramolecular recombination, versus an intermolecular gene conversion or crossover mechanism. 30 nt gaps formed a distinct 30 nt deletion product, whereas single strand nicks and gaps of 15 nts did not induce expansions or deletions. Formation of this deletion product required the CTG•CAG repeats to be present in the single-stranded region and was stimulated by E. coli DNA ligase, but was not dependent upon the RecFOR pathway. Models are presented to explain the DSB induced instabilities and formation of the 30 nucleotide deletion product. In addition to the in vitro creation of DSBs, several attempts to generate this incision in vivo with the use of EcoR I restriction modification systems were conducted. ^

Breast cancer risk and heterocyclic amines and effect modification by NAT gene polymorphisms

Breast cancer is the most common cancer in women in the United States and is a leading cause of cancer-related deaths (1). Recently, dietary heterocyclic amines (HCAs) have been proposed to be a risk factor for breast cancer (2). This study uses the data collected for a case-control study conducted at the M.D. Anderson Cancer Center to assess the association between breast cancer risk and HCAs {2-amino-1-methyl-6-phenylimidazole [4,5-b] pyridine (PhIP), 2-amino-3,8-dimethylimidazo [4,5-f] quinoxaline (MeIQx), 2-amino-3,4,8-trimethylimidazo [4,5-f] quinoxaline (DiMeIQx) and mutagenicity of HCAs} and to examine if this association is modified by genetic polymorphisms of N-acetyl transferases (NAT1/NAT2). The NAT1/2 genotype was determined using Taqman technology. HCAs were estimated by using a meat preparation questionnaire on meat type, cooking method, and doneness, combined with a quantitative HCA database. Three hundred and fifty patients with breast cancer attending the Diagnostic Radiology Clinic at M. D. Anderson Cancer Center and fulfilling the eligibility criteria were compared to three hundred and fifty patients attending the same clinic for benign breast lesions to answer these questions. Logistic regression models were used to control for known risk factors and showed no statistically significant association between breast cancer versus benign breast cancer lesions and dietary intake of heterocyclic amines. There was no clear difference in their effect after subgroup analyses in different acetylator strata of NAT1/2 and no statistical interactions were found between NAT1/2 genotypes and HCAs, suggesting no effect modification by NAT1/2 acetylator status. These results suggest the need for further research to analyze if these null associations were because of the benign breast lesions sharing the risk factors with breast cancer or any other factors which haven't been explored yet.^

Cis-acting elements and developmental regulators govern toxin gene expression in Bacillus anthracis

Expression of the structural genes for the anthrax toxin proteins is coordinately controlled by host-related signals such as elevated CO2 , and the trans-acting positive regulator, AtxA. No specific binding of AtxA to the toxin gene promoters has been demonstrated and no sequence-based similarities are apparent in the promoter regions of toxin genes. We hypothesized that the toxin genes possess common structural features that are required for positive regulation. To test this hypothesis, I performed an extensive characterization of the toxin gene promoters. I determined the minimal sequences required for atxA-mediated toxin gene expression and compared these sequences for structural similarities. In silico modeling and in vitro experiments indicated significant curvature within these regions. Random mutagenesis revealed that point mutations associated with reduced transcriptional activity, mostly mapped to areas of high curvature. This work enabled the identification of two potential cis-acting elements implicated in AtxA-mediated regulation of the toxin genes. In addition to the growth condition requirements and AtxA, toxin gene expression is under growth phase regulation. The transition state regulator AbrB represses atxA expression to influence toxin synthesis. Here I report that toxin gene expression also requires sigH, a gene encoding the RNA polymerase sigma factor associated with development in B. subtilis. In the well-studied B. subtilis system, σH is part of a feedback control pathway that involves AbrB and the major response regulator of sporulation initiation, Spo0A. My data indicate that in B. anthracis, regulatory relationships exist between these developmental regulators and atxA . Interestingly, during growth in toxin-inducing conditions, sigH and abrB expression deviates from that described for B. subtilis, affecting expression of the atxA gene. These findings, combined with previous observations, suggest that the steady state level of atxA expression is critical for optimal toxin gene transcription. I propose a model whereby, under toxin-inducing conditions, control of toxin gene expression is fine-tuned by the independent effects of the developmental regulators on the expression of atxA . The growth condition-dependent changes in expression of these regulators may be crucial for the correct timing and uninterrupted expression of the toxin genes during infection. ^

Characterization of glucocorticoid receptor (NR3C1) gene polymorphisms and a family-based association study with hypertension

Hypertension is usually defined as having values of systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg. Hypertension is one of the main adverse effects of glucocorticoid on the cardiovascular system. Glucocorticoids are essential hormones, secreted from adrenal glands in circadian fashion. Glucocorticoid's effect on blood pressure is conveyed by the glucocorticoid receptor (NR3C1), an omnipresent nuclear transcription factor. Although polymorphisms in this gene have long been implicated to be a causal factor for cardiovascular diseases such as hypertension, no study has yet thoroughly interrogated the gene's polymorphisms for their effect on blood pressure levels. Therefore, I have first resequenced ∼30 kb of the gene, encompassing all exons, promoter regions, 5'/3' UTRs as well as at least 1.5 kb of the gene's flanking regions from 114 chromosome 5 monosomic cell lines, comprised of three major American ethnic groups—European American, African American and Mexican American. I observed 115 polymorphisms and 14 common molecularly phased haplotypes. A subset of markers was chosen for genotyping study populations of GENOA (Genetic Epidemiology Network of Atherosclerosis; 1022 non-Hispanic whites, 1228 African Americans and 954 Mexican Americans). Since these study populations include sibships, the family-based association test was performed on 4 blood pressure-related quantitative variables—pulse, systolic blood pressure, diastolic blood pressure and mean arterial pressure. Using these analyses, multiple correlated SNPs are significantly protective against high systolic blood pressure in non-Hispanic whites, which includes rsb198, a SNP formerly associated with beneficial body compositions. Haplotype association analysis also supports this finding and all p-values remained significant after permutation tests. I therefore conclude that multiple correlated SNPs on the gene may confer protection against high blood pressure in non-Hispanic whites. ^

The molecular basis for beta-lactamase gene expression in B. anthracis, B. cereus and B. thuringiensis

The susceptibility of most Bacillus anthracis strains to β-lactam antibiotics is intriguing considering that the B. anthracis genome harbors two β-lactamase genes, bla1 and bla2, and closely-related species, Bacillus cereus and Bacillus thuringiensis, typically produce β-lactamases. This work demonstrates that B. anthracis bla expression is affected by two genes, sigP and rsp, predicted to encode an extracytoplasmic function sigma factor and an antisigma factor, respectively. Deletion of the sigP/rsp locus abolished bla expression in a penicillin-resistant clinical isolate and had no effect on bla expression in a prototypical penicillin-susceptible strain. Complementation with sigP/rsp from the penicillin-resistant strain, but not the penicillin-susceptible strain, conferred β-lactamase activity upon both mutants. These results are attributed to a nucleotide deletion near the 5' end of rsp in the penicillin-resistant strain that is predicted to result in a nonfunctional protein. B. cereus and B. thuringiensis sigP and rsp homologues are required for inducible penicillin resistance in those species. Expression of the B. cereus or B. thuringiensis sigP and rsp genes in a B. anthracis sigP/rsp-null mutant confers resistance to β-lactam antibiotics, suggesting that while B. anthracis contains the genes necessary for sensing β-lactam antibiotics, the B. anthracis sigP/rsp gene products are insufficient for bla induction. ^ Because alternative sigma factors recognize unique promoter sequence, direct targets can be elucidated by comparing transcriptional profiling results with an in silico search using the sigma factor binding sequence. Potential σP -10 and -35 promoter elements were identified upstream from bla1 bla2 and sigP. Results obtained from searching the B. anthracis genome with the conserved sequences were evaluated against transcriptional profiling results comparing B. anthracis 32 and an isogenic sigP/rsp -null strain. Results from these analyses indicate that while the absence of the sigP gene significantly affects the transcript levels of 16 genes, only bla1, bla2 and sigP are directly regulated by σP. The genomes of B. cereus and B. thuringiensis strains were also analyzed for the potential σP binding elements. The sequence was located upstream from the sigP and bla genes, and previously unidentified genes predicted to encode a penicillin-binding protein (PBP) and a D-alanyl-D-alanine carboxypeptidase, indicating that the σ P regulon in these species responds to cell-wall stress caused by β-lactam antibiotics. ^ β-lactam antibiotics prevent attachment of new peptidoglycan to the cell wall by blocking the active site of PBPs. A B. cereus and B. thuringiensis pbp-encoding gene located near bla1 contains a potential σP recognition sequence upstream from the annotated translational start. Deletion of this gene abolished β-lactam resistance in both strains. Mutations in the active site of the PBP were detrimental to β-lactam resistance in B. cereus, but not B. thuringiensis, indicating that the transpeptidase activity is only important in B. cereus. I also found that transcript levels of the PBP-encoding gene are not significantly affected by the presence of β-lactam antibiotic. Based on these data I hypothesize that the gene product acts a sensor of β-lactam antibiotic. ^

Mechanism of methylation gene reactivation in human cancer cells

Epigenetic silencing of tumor suppressor genes by DNA hypermethylation at promoter regions is a common event in carcinogenesis and tumor progression. Abrogation of methylation and reversal of epigenetic silencing is a very potent way in cancer treatment. However, the reactivation mechanisms are poorly understood. In this study, we first developed a cell line model system named YB5, derived from SW48 cancer cell line, which bears one copy of stably integrated EGFP gene on Chromosome 1p31.1 region. The GFP gene expression is transcriptionally silenced due to the hypermethylated promoter CMV. However, the GFP expression can be restored using demethylating agent 5-aza-2' deoxycytidine (DAC), and detected by FACS and fluorescent microscopy. Using this system, we observed the heterogeneous reactivation induced by DAC treatment. After flow sorting, GFP negative cells exhibited similar level of incomplete demethylation compared to GFP positive cells on repetitive LINE1 element, tumor suppressor genes such as P16, CDH13, and RASSF1a, and CMV promoter as well. However, the local chromatin of CMV-GFP locus altered to an open structure marked by high H3 lysine 9 acetylation and low H3 lysine 27 tri-methylation in GFP positive cells, while the GFP negative cells retained mostly the original repressive marks. Thus, we concluded that DAC induced DNA hypomethylation alone does not directly determine the level of re-expression, and the resetting of the local chromatin structure under hypomethylation environment is required for gene reactivation. Besides, a lentivirus vector-based shRNA screening was performed using the YB5 system. Although it is the rare chance that vector lands in the neighboring region of GFP, we found that the exogenous vector DNA inserted into the upstream region of GFP gene locus led to the promoter demethylation and reactivated the silenced GFP gene. Thus, epigenetic state can be affected by changing of the adjacent nucleic acid sequences. Further, this hypermethylation silenced system was utilized for epigenetic drug screening. We have found that DAC combined with carboplatin would enhance the GFP% yield and increase expression of other tumor suppressor genes than DAC alone, and this synergistic effect may be related to DNA repair process. In summary, these studies reveal that reversing of methylation silencing requires coordinated alterations of DNA methylation, chromatin structure, and local microenvironment. ^

Characterization of cytochrome P450 4F subfamily: Response in traumatic brain injury and gene regulation

CYP4F enzymes metabolize endogenous molecules including arachidonic acid, leukotrienes and prostaglandins. The involvement of these eisosanoids in inflammation has led to the hypothesis that CYP4Fs may modulate inflammatory conditions after traumatic brain injury (TBI). In rat, TBI elicited changes in mRNA expression of CYP4Fs as a function of time in the cerebrum region. These changes in CYP4F mRNA levels inversely correlated with the cerebral leukotriene B4 (LTB4) level following injury at the same time points. TBI also resulted in changes in CYP4F protein expression and localization around the injury site, where CYP4F1 and CYP4F6 immunoreactivity increased in surrounding astrocytes and CYP4F4 immunoreactivity shifted from endothelia of cerebral vessels to astrocytes. The study with rat primary astrocytes indicated that pro-inflammatory cytokines TNFα and IL-1β could affect the transcription of CYP4Fs to a certain degree, whereas the changing pattern in the primary astrocytes appeared to be different from that in the in vivo TBI model.^ In addition, the regulation of CYP4F genes has been an unsolved issue although factors including cytokines and fatty acids appear to affect CYP4Fs expression in multiple models. In this project, HaCaT cells were used as an in vitro cellular model to define signaling pathways involved in the regulation of human CYP4F genes. Retinoic acids inhibited CYP4F11 expression, whereas cytokines TNFα and IL-1β induced transcription of CYP4F11 in HaCaT cells. The induction of CYP4F11 by both cytokines could be blocked by a JNK specific inhibitor, indicating the involvement of the JNK pathway in the up-regulation of CYP4F11. Retinoic acids are known to function in gene regulation through nuclear receptors RARs and RXRs. The RXR agonist LG268 greatly induced transcription of CYP4F11, whereas RAR agonist TTNPB obviously inhibited CYP4F11 transcription, indicating that the down-regulation of CYP4F11 by retinoic acid was mediated by RARs, and that inhibition of CYP4F11 by retinoic acid may also be related to the competition for RXR receptors. Thus, the CYP4F11 gene is regulated by signaling pathways including the RXR pathway and the JNK pathway. In contrast, the regulation mechanism of other CYP4Fs by retinoic acids appears to be different from that of CYP4F11.^

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.^

Systems Biology Approaches to Probe Gene Regulation in Bacteria

Mechanisms that allow pathogens to colonize the host are not the product of isolated genes, but instead emerge from the concerted operation of regulatory networks. Therefore, identifying components and the systemic behavior of networks is necessary to a better understanding of gene regulation and pathogenesis. To this end, I have developed systems biology approaches to study transcriptional and post-transcriptional gene regulation in bacteria, with an emphasis in the human pathogen Mycobacterium tuberculosis (Mtb). First, I developed a network response method to identify parts of the Mtb global transcriptional regulatory network utilized by the pathogen to counteract phagosomal stresses and survive within resting macrophages. As a result, the method unveiled transcriptional regulators and associated regulons utilized by Mtb to establish a successful infection of macrophages throughout the first 14 days of infection. Additionally, this network-based analysis identified the production of Fe-S proteins coupled to lipid metabolism through the alkane hydroxylase complex as a possible strategy employed by Mtb to survive in the host. Second, I developed a network inference method to infer the small non-coding RNA (sRNA) regulatory network in Mtb. The method identifies sRNA-mRNA interactions by integrating a priori knowledge of possible binding sites with structure-driven identification of binding sites. The reconstructed network was useful to predict functional roles for the multitude of sRNAs recently discovered in the pathogen, being that several sRNAs were postulated to be involved in virulence-related processes. Finally, I applied a combined experimental and computational approach to study post-transcriptional repression mediated by small non-coding RNAs in bacteria. Specifically, a probabilistic ranking methodology termed rank-conciliation was developed to infer sRNA-mRNA interactions based on multiple types of data. The method was shown to improve target prediction in Escherichia coli, and therefore is useful to prioritize candidate targets for experimental validation.

Detecting genetic and nutritional lung cancer risk factors related to folate metabolism using Bayesian generalized linear models

Complex diseases, such as cancer, are caused by various genetic and environmental factors, and their interactions. Joint analysis of these factors and their interactions would increase the power to detect risk factors but is statistically. Bayesian generalized linear models using student-t prior distributions on coefficients, is a novel method to simultaneously analyze genetic factors, environmental factors, and interactions. I performed simulation studies using three different disease models and demonstrated that the variable selection performance of Bayesian generalized linear models is comparable to that of Bayesian stochastic search variable selection, an improved method for variable selection when compared to standard methods. I further evaluated the variable selection performance of Bayesian generalized linear models using different numbers of candidate covariates and different sample sizes, and provided a guideline for required sample size to achieve a high power of variable selection using Bayesian generalize linear models, considering different scales of number of candidate covariates. ^ Polymorphisms in folate metabolism genes and nutritional factors have been previously associated with lung cancer risk. In this study, I simultaneously analyzed 115 tag SNPs in folate metabolism genes, 14 nutritional factors, and all possible genetic-nutritional interactions from 1239 lung cancer cases and 1692 controls using Bayesian generalized linear models stratified by never, former, and current smoking status. SNPs in MTRR were significantly associated with lung cancer risk across never, former, and current smokers. In never smokers, three SNPs in TYMS and three gene-nutrient interactions, including an interaction between SHMT1 and vitamin B12, an interaction between MTRR and total fat intake, and an interaction between MTR and alcohol use, were also identified as associated with lung cancer risk. These lung cancer risk factors are worthy of further investigation.^