Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.

Computational identification and functional validation of regulatory motifs in cartilage-expressed genes.

Chondrocyte gene regulation is important for the generation and maintenance of cartilage tissues. Several regulatory factors have been identified that play a role in chondrogenesis, including the positive transacting factors of the SOX family such as SOX9, SOX5, and SOX6, as well as negative transacting factors such as C/EBP and delta EF1. However, a complete understanding of the intricate regulatory network that governs the tissue-specific expression of cartilage genes is not yet available. We have taken a computational approach to identify cis-regulatory, transcription factor (TF) binding motifs in a set of cartilage characteristic genes to better define the transcriptional regulatory networks that regulate chondrogenesis. Our computational methods have identified several TFs, whose binding profiles are available in the TRANSFAC database, as important to chondrogenesis. In addition, a cartilage-specific SOX-binding profile was constructed and used to identify both known, and novel, functional paired SOX-binding motifs in chondrocyte genes. Using DNA pattern-recognition algorithms, we have also identified cis-regulatory elements for unknown TFs. We have validated our computational predictions through mutational analyses in cell transfection experiments. One novel regulatory motif, N1, found at high frequency in the COL2A1 promoter, was found to bind to chondrocyte nuclear proteins. Mutational analyses suggest that this motif binds a repressive factor that regulates basal levels of the COL2A1 promoter.

Eomesodermin, a target gene of Pou4f2, is required for retinal ganglion cell and optic nerve development in the mouse.

The mechanisms regulating retinal ganglion cell (RGC) development are crucial for retinogenesis and for the establishment of normal vision. However, these mechanisms are only vaguely understood. RGCs are the first neuronal lineage to segregate from pluripotent progenitors in the developing retina. As output neurons, RGCs display developmental features very distinct from those of the other retinal cell types. To better understand RGC development, we have previously constructed a gene regulatory network featuring a hierarchical cascade of transcription factors that ultimately controls the expression of downstream effector genes. This has revealed the existence of a Pou domain transcription factor, Pou4f2, that occupies a key node in the RGC gene regulatory network and that is essential for RGC differentiation. However, little is known about the genes that connect upstream regulatory genes, such as Pou4f2 with downstream effector genes responsible for RGC differentiation. The purpose of this study was to characterize the retinal function of eomesodermin (Eomes), a T-box transcription factor with previously unsuspected roles in retinogenesis. We show that Eomes is expressed in developing RGCs and is a mediator of Pou4f2 function. Pou4f2 directly regulates Eomes expression through a cis-regulatory element within a conserved retinal enhancer. Deleting Eomes in the developing retina causes defects reminiscent of those in Pou4f2(-/-) retinas. Moreover, myelin ensheathment in the optic nerves of Eomes(-/-) embryos is severely impaired, suggesting that Eomes regulates this process. We conclude that Eomes is a crucial regulator positioned immediately downstream of Pou4f2 and is required for RGC differentiation and optic nerve development.

Role of acetyl-phosphate in activation of the Rrp2-RpoN-RpoS pathway in Borrelia burgdorferi.

Borrelia burgdorferi, the Lyme disease spirochete, dramatically alters its transcriptome and proteome as it cycles between the arthropod vector and mammalian host. During this enzootic cycle, a novel regulatory network, the Rrp2-RpoN-RpoS pathway (also known as the σ(54)-σ(S) sigma factor cascade), plays a central role in modulating the differential expression of more than 10% of all B. burgdorferi genes, including the major virulence genes ospA and ospC. However, the mechanism(s) by which the upstream activator and response regulator Rrp2 is activated remains unclear. Here, we show that none of the histidine kinases present in the B. burgdorferi genome are required for the activation of Rrp2. Instead, we present biochemical and genetic evidence that supports the hypothesis that activation of the Rrp2-RpoN-RpoS pathway occurs via the small, high-energy, phosphoryl-donor acetyl phosphate (acetyl∼P), the intermediate of the Ack-Pta (acetate kinase-phosphate acetyltransferase) pathway that converts acetate to acetyl-CoA. Supplementation of the growth medium with acetate induced activation of the Rrp2-RpoN-RpoS pathway in a dose-dependent manner. Conversely, the overexpression of Pta virtually abolished acetate-induced activation of this pathway, suggesting that acetate works through acetyl∼P. Overexpression of Pta also greatly inhibited temperature and cell density-induced activation of RpoS and OspC, suggesting that these environmental cues affect the Rrp2-RpoN-RpoS pathway by influencing acetyl∼P. Finally, overexpression of Pta partially reduced infectivity of B. burgdorferi in mice. Taken together, these findings suggest that acetyl∼P is one of the key activating molecule for the activation of the Rrp2-RpoN-RpoS pathway and support the emerging concept that acetyl∼P can serve as a global signal in bacterial pathogenesis.

Identification of genetic variation influencing apolipoprotein E levels: Follow-up of genome-wide linkage scans in the GENOA study

Apolipoprotein E (ApoE) plays a major role in the metabolism of high density and low density lipoproteins (HDL and LDL). Its common protein isoforms (E2, E3, E4) are risk factors for coronary artery disease (CAD) and explain between 16 to 23% of the inter-individual variation in plasma apoE levels. Linkage analysis has been completed for plasma apoE levels in the GENOA study (Genetic Epidemiology Network of Atherosclerosis). After stratification of the population by lipoprotein levels and body mass index (BMI) to create more homogeneity with regard to biological context for apoE levels, Hispanic families showed significant linkage on chromosome 17q for two strata (LOD=2.93 at 104 cM for a low cholesterol group, LOD=3.04 at 111 cM for a low cholesterol, high HDLC group). Replication of 17q linkage was observed for apoB and apoE levels in the unstratified Hispanic and African-American populations, and for apoE levels in African-American families. Replication of this 17q linkage in different populations and strata provides strong support for the presence of gene(s) in this region with significant roles in the determination of inter-individual variation in plasma apoE levels. Through a positional and functional candidate gene approach, ten genes were identified in the 17q linked region, and 62 polymorphisms in these genes were genotyped in the GENOA families. Association analysis was performed with FBAT, GEE, and variance-component based tests followed by conditional linkage analysis. Association studies with partial coverage of TagSNPs in the gene coding for apolipoprotein H (APOH) were performed, and significant results were found for 2 SNPs (APOH_20951 and APOH_05407) in the Hispanic low cholesterol strata accounting for 3.49% of the inter-individual variation in plasma apoE levels. Among the other candidate genes, we identified a haplotype block in the ACE1 gene that contains two major haplotypes associated with apoE levels as well as total cholesterol, apoB and LDLC levels in the unstratified Hispanic population. Identifying genes responsible for the remaining 60% of inter-individual variation in plasma apoE level, will yield new insights into the understanding of genetic interactions involved in the lipid metabolism, and a more precise understanding of the risk factors leading to CAD. ^

Control of the master virulence regulatory gene atxA in Bacillus anthracis

Transcription of the Bacillus anthracis structural genes for the anthrax toxin proteins and biosynthetic operon for capsule are positively regulated by AtxA, a transcription regulator with unique properties. Consistent with the role of atxA in virulence factor expression, a B. anthracis atxA-null mutant is avirulent in a murine model for anthrax. In batch culture, multiple signals impact atxA transcript levels, and the timing and steady state level of atxA expression is critical for optimal toxin and capsule synthesis. Despite the apparent complex control of atxA transcription, only one trans-acting protein, the transition state regulator AbrB, has been demonstrated to directly interact with the atxA promoter. The AbrB-binding site has been described, but additional cis-acting control sequences have not been defined. Using transcriptional lacZ fusions, electrophoretic mobility shift assays, and Western blot analysis, the cis-acting elements and trans-acting factors involved in regulation of atxA in B. anthracis strains containing either both virulence plasmids, pXO1 and pXO2, or only one plasmid, pXO1, were studied. This work demonstrates that atxA transcription from the major start site P1 is dependent upon a consensus sequence for the housekeeping sigma factor SigA, and an A+T-rich upstream element (UP-element) for RNA polymerase (RNAP). In addition, the data show that a trans-acting protein(s) other than AbrB negatively impacts atxA transcription when it binds specifically to a 9-bp palindrome within atxA promoter sequences located downstream of P1. Mutation of the palindrome prevents binding of the trans-acting protein(s) and results in a corresponding increase in AtxA and anthrax toxin production in a strain- and culture-dependent manner. The identity of the trans-acting repressor protein(s) remains elusive; however, phenotypes associated with mutation of the repressor binding site have revealed that the trans-acting repressor protein(s) indirectly controls B. anthracis development. Mutation of the repressor binding site results in misregulation and overexpression of AtxA in conditions conducive for development, leading to a marked sporulation defect that is both atxA- and pXO2-61-dependent. pXO2-61 is homologous to the sensor domain of sporulation sensor histidine kinases and is proposed to titrate an activating signal away from the sporulation phosphorelay when overexpressed by AtxA. These results indicate that AtxA is not only a master virulence regulator, but also a modulator of proper B. anthracis development. Also demonstrated in this work is the impact of the developmental regulators AbrB, Spo0A, and SigH on atxA expression and anthrax toxin production in a genetically incomplete (pXO1+, pXO2-) and genetically complete (pXO1+, pXO2+) strain background. AtxA and anthrax toxin production resulting from deletion of the developmental regulators are strain-dependent suggesting that factors on pXO2 are involved in control of atxA. The only developmental deletion mutant that resulted in a prominent and consistent strain-independent increase in AtxA protein levels was an abrB-null mutant. As a result of increased AtxA levels, there is early and increased production of anthrax toxins in an abrB-null mutant. In addition, the abrB-null mutant exhibited an increase in virulence in a murine model for anthrax. In contrast, virulence of the atxA promoter mutant was unaffected in a murine model for anthrax despite the production of 5-fold more AtxA than the abrB-null mutant. These results imply that AtxA is not the only factor impacting pathogenesis in an abrB-null mutant. Overall, this work highlights the complex regulatory network that governs expression of atxA and provides an additional role for AtxA in B. anthracis development.

The SasS -SasR two-component signal transduction system required for early Myxococcus xanthus multicellular *development

Initiation of Myxococcus xanthus multicellular development requires both nutrient limitation and high cell density. The extracellular signal, A signal, which consists of a set of amino acids at specific concentrations, serves as a cell density signal in M. xanthus early development. A reporter gene, designated 4521, that requires both starvation and A signal for developmental expression was used to identify mutations in the signal transduction pathways. A group of point mutations located in the chromosomal sasB locus that bypasses both requirements was previously isolated. One of these point mutations, sasB7, was mapped to the sasS gene, which is predicted to encode a transmembrane histidine protein kinase required for normal development. SasS is a positive regulator of 4521 and a candidate A signal sensor. This dissertation continues the characterization of the sasB locus, focusing on the sasR gene and the functional relationship of SasS and SasR. ^ The sasR gene is located 2.2-kb downstream of sasS. It is predicted to encode an NtrC-like response regulator, which belongs to the family of sigma54 transcriptional activators. SasR is a positive regulator of 4521 gene and is required for normal development. The sasR mutant displays phenotypes similar to that of sasS mutant. Both SasS and SasR are required for the A-signal-dependent 4521 expression. Genetic epistasis analysis indicates that SasR functions downstream of SasS. Biochemical studies show that SasS has autokinase activity, and phosphorylated SasS is able to transfer its phosphate to SasR. We propose that SasS and SasR form a two-component signal transduction system in the A signal transduction pathway. ^ To search for the genes regulated by SasS and SasR, expression patterns of a group of developmental genes were compared in wild-type and sasS null mutant backgrounds. SasS and SasR were found to positively regulate sasN and 4521. The sasN gene was previously identified as a negative regulator of 4521, located at about 170-bp downstream of sasR. It is required for normal fruiting body development. Based on the above data, a regulatory network consisting of sasS, sasR, sasN, and 4521 is hypothesized, and the interactions of the components in this network can now be further studied. ^

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

REST maintains self-renewal and pluripotency of embryonic stem cells.

The neuronal repressor REST (RE1-silencing transcription factor; also called NRSF) is expressed at high levels in mouse embryonic stem (ES) cells, but its role in these cells is unclear. Here we show that REST maintains self-renewal and pluripotency in mouse ES cells through suppression of the microRNA miR-21. We found that, as with known self-renewal markers, the level of REST expression is much higher in self-renewing mouse ES cells than in differentiating mouse ES (embryoid body, EB) cells. Heterozygous deletion of Rest (Rest+/-) and its short-interfering-RNA-mediated knockdown in mouse ES cells cause a loss of self-renewal-even when these cells are grown under self-renewal conditions-and lead to the expression of markers specific for multiple lineages. Conversely, exogenously added REST maintains self-renewal in mouse EB cells. Furthermore, Rest+/- mouse ES cells cultured under self-renewal conditions express substantially reduced levels of several self-renewal regulators, including Oct4 (also called Pou5f1), Nanog, Sox2 and c-Myc, and exogenously added REST in mouse EB cells maintains the self-renewal phenotypes and expression of these self-renewal regulators. We also show that in mouse ES cells, REST is bound to the gene chromatin of a set of miRNAs that potentially target self-renewal genes. Whereas mouse ES cells and mouse EB cells containing exogenously added REST express lower levels of these miRNAs, EB cells, Rest+/- ES cells and ES cells treated with short interfering RNA targeting Rest express higher levels of these miRNAs. At least one of these REST-regulated miRNAs, miR-21, specifically suppresses the self-renewal of mouse ES cells, corresponding to the decreased expression of Oct4, Nanog, Sox2 and c-Myc. Thus, REST is a newly discovered element of the interconnected regulatory network that maintains the self-renewal and pluripotency of mouse ES cells.

Evolution of sensory complexity recorded in a myxobacterial genome.

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced delta-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell-cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

Diversity of neuronal connexins in the mammalian retina

Many neurons in the mammalian retina are electrically coupled by intercellular channels or gap junctions, which are assembled from a family of proteins called connexins. Numerous studies indicate that gap junctions differ in properties such as conductance and tracer permeability. For example, A-type horizontal cell gap junctions are permeable to Lucifer Yellow, but B-type horizontal cell gap junctions are not. This suggests the two cell types express different connexins. My hypothesis is that multiple neuronal connexins are expressed in the mammalian retina in a cell type specific manner. Immunohistochemical techniques and confocal microscopy were used to localize certain connexins within well-defined neuronal circuits. The results of this study can be summarized as follows: AII amacrine cells, which receive direct input from rod bipolar cells, are well-coupled to neighboring AIIs. In addition, AII amacrine cells also form gap junctions with ON cone bipolar cells. This is a complex heterocellular network. In both rabbit and primate retina, connexin36 occurs at dendritic crossings in the AII matrix as well as between AIIs and ON cone bipolar cells. Coupling in the AII network is thought to reduce noise in the rod pathway while AII/bipolar gap junctions are required for the transmission of rod signals to ON ganglion cells. In the outer plexiform layer, connexin36 forms gap junctions between cones and between rods and cones via cone telodendria. Cone to cone coupling is thought to reduce noise and is partly color selective. Rod to cone coupling forms an alternative rod pathway thought to operate at intermediate light intensity. A-type horizontal cells in the rabbit retina are strongly coupled via massive low resistance gap junctions composed from Cx50. Coupling dramatically extends the receptive field of horizontal cells and the modulation of coupling is thought to change the strength of the feedback signal from horizontal cells to cones. Finally, there are other coupled networks, such as B-type horizontal cells and S1/S2 amacrine cells, which do not use either connexin36 or Cx50. These results confirm the hypothesis that multiple neuronal connexins are expressed in the mammalian retina and these connexins are localized to particular retinal circuits. ^

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

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.

REGULATION OF ALTERNATIVE CARBON METABOLISM IN CANDIDA ALBICANS

Candida albicans is the most important fungal pathogen of humans. Transcript profiling studies show that upon phagocytosis by macrophages, C. albicans undergoes a massive metabolic reorganization activating genes involved in alternative carbon metabolism, including the glyoxylate cycle, β-oxidation and gluconeogenesis. Mutations in key enzymes such as ICL1 (glyoxylate cycle) and FOX2 (fatty acid β-oxidation) revealed that alternative carbon metabolic pathways are required for full virulence in C. albicans. These studies indicate C. albicans uses non-preferred carbon sources allowing its adaptation to microenvironments were nutrients are scarce. It has become apparent that the regulatory networks required for regulation of alternative carbon metabolism in C. albicans are considerably different from the Saccharomyces cerevisiae paradigm and appear more analogous to the Aspergillus nidulans systems. Well-characterized transcription factors in S. cerevisiae have no apparent phenotype or are missing in C. albicans. CTF1 was found to be a single functional homolog of the A. nidulans FarA/FarB proteins, which are transcription factors required for fatty acid utilization. Both FOX2 and ICL1 were found to be part of a large CTF1 regulon. To increase our understanding of how CTF1 regulates its target genes, including whether regulation is direct or indirect, the FOX2 and ICL1 promoter regions were analyzed using a combination of bioinformatics and promoter deletion analysis. To begin characterizing the FOX2 and ICL1 promoters, 5’ rapid amplification of cDNA ends (5’RACE) was used to identify two transcriptional initiation sites in FOX2 and one in ICL1. GFP reporter assays show FOX2 and ICL1 are rapidly expressed in the presence of alternative carbon sources. Both FOX2 and ICL1 harbor the CCTCGG sequence known to be bound by the Far proteins, hence rendering the motif as a putative CTF1 DNA binding element. In this study, the CCTCGG sequence was found to be essential for FOX2 regulation. However, this motif does not appear to be equally important for the regulation of ICL1. This study supports the notion that although C. albicans has diverged from the paradigms of model fungi, C. albicans has made specific adaptations to its transcription-based regulatory network that may contribute to its metabolic flexibility.