Factors associated with needle sharing among Black male injection drug users in Harris County, Texas

Background. Injection drug users (IDUs) are at increased risk for HIV transmission due to unique risk behaviors, such as sharing needles. In Houston, IDUs account for 18% of all HIV/AIDS cases among Black males. ^ Objectives. This analysis compared demographic, behavioral, and psychosocial characteristics of needle sharing and non-sharing IDUs in a population of Black males in Harris County, Texas. ^ Methods. Data used for this analysis were from the second IDU cycle of the National HIV Behavioral Surveillance System. This dataset included a sample of 288 Black male IDUs. Univariate and multivariate statistical analysis were performed to determine statistically significant associations of needle sharing in this population and to create a functional model to inform local HIV prevention programs. ^ Results. Half of the participants in this analysis shared needles in the past 12 months. Compared to non-sharers, sharers were more likely to be homeless (OR=3.70, p<0.01) or arrested in the past year (OR=2.31, p<0.01), inject cocaine (OR=2.07, p<0.01), report male-to-male sex in the past year (OR=6.97, p<0.01), and to exchange sex for money or drugs. Sharers were less likely than non-sharers to graduate high school (OR=0.36, p<0.01), earn $5,000 or more a year (OR=1.15, p=0.05), get needles from a medical source (OR=0.59, p=0.03), and ever test for HIV (OR=0.17, p<0.01). Sharers were more likely to report depressive symptoms (OR=3.49, p<0.01), lower scores on the family support scale (mean difference 0.41, p=0.01) and decision-making confidence scale (mean difference 0.38, p<0.01), and greater risk-taking (mean difference -0.49, p<0.01) than non-sharers. In a multivariable logistic regression, sharers were less likely to have graduated high school (OR=0.33, p<0.01) and have been tested for HIV (OR=0.12, p<0.01) and were more likely to have been arrested in the past year (OR=2.3, p<0.01), get needles from a street source (OR=3.87, p<0.01), report male-to-male sex (OR=7.01, p<0.01), and have depressive symptoms (OR=2.36, p=0.02) and increased risk-taking (OR=1.78, p=0.01). ^ Conclusions. IDUs that shared needles are different from those that did not, reporting lower socioeconomic status, increased sexual and risk behaviors, increased depressive symptoms and increased risk-taking. These findings suggest that intervention programs that also address these demographic, behavioral, and psychosocial factors may be more successful in decreasing needle sharing among this population.^

THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

AMPA RECEPTOR ALLOSTERISM: MEASUREMENT OF THE CONFORMATIONAL CHANGES IN THE LIGAND BINDING DOMAIN OF A FUNCTIONAL RECEPTOR

Ionotropic glutamate receptors are important excitatory neurotransmitter receptors in the mammalian central nervous system that have been implicated in a number of neuropathologies such as epilepsy, ischemia, and amyotrophic lateral sclerosis. Glutamate binding to an extracellular ligand binding domain initiates a series of structural changes that leads to the formation of a cation selective transmembrane channel, which consequently closes due to desensitization of the receptor. The crystal structures of the AMPA subtype of the glutamate receptor have been particularly useful in providing initial insight into the conformational changes in the ligand binding domain; however, these structures are limited by crystallographic constraint. To gain a clear picture of how agonist binding is coupled to channel activation and desensitization, it is essential to study changes in the ligand binding domain in a dynamic, physiological state. In this dissertation, a technique called Luminescence Resonance Energy Transfer was used to determine the conformational changes associated with activation and desensitization in a functional AMPA receptor (ÄN*-AMPA) that contains the ligand binding domain and transmembrane segments; ÄN*-AMPA has been modified such that fluorophores can be introduced at specific sites to serve as a readout of cleft closure or to establish intersubunit distances. Previous structural studies of cleft closure of the isolated ligand binding domain in conjunction with functional studies of the full receptor suggest that extent of cleft closure correlates with extent of activation. Here, LRET has been used to show that a similar relationship between cleft closure and activation is observed in the “full length” receptor showing that the isolated ligand binding domain is a good model of the domain in the full length receptor for changes within a subunit. Similar LRET investigations were used to study intersubunit distances specifically to probe conformational changes between subunits within a dimer in the tetrameric receptor. These studies show that the dimer interface is coupled in the open state, and decoupled in the desensitized state, similar to the isolated ligand binding domain crystal structure studies. However, we show that the apo state dimer interface is not pre-formed as in the crystal structure, hence suggesting a mechanism for functional transitions within the receptor based on LRET distances obtained.

Molecular basis of Corynebacterium diphtheriae virulence and infection in the Caenorhabditis elegans model host

Corynebacterium diphtheriae is the causative agent of cutaneous and pharyngeal diphtheria in humans. While lethality is certainly caused by diphtheria toxin, corynebacterial colonization may primarily require proteinaceous fibers called pili, which mediate adherence to specific tissues. The type strain of C. diphtheriae possesses three distinct pilus structures, namely the SpaA, SpaD, and SpaH-type pili, which are encoded by three distinct pilus gene clusters. The pilus is assembled onto the bacterial peptidoglycan by a specific transpeptidase enzyme called sortase. Although the SpaA pili are shown to be specific for pharyngeal cells in vitro, little is known about functions of the three pili in bacterial pathogenesis. This is mainly due to lack of in vivo models of corynebacterial infection. As an alternative to mouse models as mice do not have functional receptors for diphtheria toxin, in this study I use Caenorhabditis elegans as a model host for C. diphtheriae. A simple C. elegans model would be beneficial in determining the specific role of each pilus-type and the literature suggests that C. elegans infection model can be used to study a variety of bacterial species giving insight into bacterial virulence and host-pathogen interactions. My study examines the hypothesis that pili and toxin are major virulent determinants of C. diphtheriae in the C. elegans model host.

DEVELOPING A METHOD FOR IDENTIFYING AND REDUCING FUNCTIONAL DISCREPANCIES OF INFORMATION SYSTEMS

Currently more than half of Electronic Health Record (EHR) projects fail. Most of these failures are not due to flawed technology, but rather due to the lack of systematic considerations of human issues. Among the barriers for EHR adoption, function mismatching among users, activities, and systems is a major area that has not been systematically addressed from a human-centered perspective. A theoretical framework called Functional Framework was developed for identifying and reducing functional discrepancies among users, activities, and systems. The Functional Framework is composed of three models – the User Model, the Designer Model, and the Activity Model. The User Model was developed by conducting a survey (N = 32) that identified the functions needed and desired from the user’s perspective. The Designer Model was developed by conducting a systemic review of an Electronic Dental Record (EDR) and its functions. The Activity Model was developed using an ethnographic method called shadowing where EDR users (5 dentists, 5 dental assistants, 5 administrative personnel) were followed quietly and observed for their activities. These three models were combined to form a unified model. From the unified model the work domain ontology was developed by asking users to rate the functions (a total of 190 functions) in the unified model along the dimensions of frequency and criticality in a survey. The functional discrepancies, as indicated by the regions of the Venn diagrams formed by the three models, were consistent with the survey results, especially with user satisfaction. The survey for the Functional Framework indicated the preference of one system over the other (R=0.895). The results of this project showed that the Functional Framework provides a systematic method for identifying, evaluating, and reducing functional discrepancies among users, systems, and activities. Limitations and generalizability of the Functional Framework were discussed.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Differential gene expression in multistep tumorigenesis using an {\it in vitro\/} human cell model

Using a human terato-carcinoma cell line, PA-1, the functional role of the oncogenes and tumor suppressor gene involved in the multistep process of carcinogenesis have been analyzed. The expression of AP-2 was strongly correlated with the susceptibility to ras transformation. The differential responsiveness to growth factors between stage 1 ras resistant cells and stage 2 ras susceptible cells was observed, indicating that the ability of stage 2 cells to respond to the mutated ras oncogenes in transformation correlated with the ability to be stimulated by certain growth factors. Using differential screening of cDNA libraries, a number of differentially expressed cDNA clones was isolated. One of those, clone 12, is overexpressed in ras transformed stage 3 cells. The amino acid sequence of clone 12 is almost identical to a mouse LLrep3 gene that was growth-regulated, and 78% similar to a yeast ribosomal protein S4. These results suggest that the S4 gene may be involved in regulation of growth. Clone 9 is expressed in stage 1 ras resistant cells (3.5-kb and 3.0-kb transcripts) but the expression of this clone in stage 2 ras susceptible cells and stage 3 ras-transformed cells is greatly diminished. The expression of this cDNA clone was increased to at least five fold in ras resistant cells and nontumorigenic hybrids treated with retinoic acid but not increased in retinoic acid treated ras susceptible cells, ras transformed cells and the tumorigenic segregants. Partial sequence of this clone showed no homology to the sequences in Genbank. These findings suggest that clone 9 could be a suppressor gene or the genes that are involved in the biochemical pathway of tumor suppression or neurogenic differentiation. The apparent pleiotropic effect of the loss of this suppressor gene function support Harris' proposal that tumor suppressor genes regulate differentiation. The tumor suppressor gene may act as negative regulator of tumor growth by controlling gene expression in differentiation. ^

Functional and motor neuronal analysis of defensive siphon responses in {\it Aplysia californica}

One of the central goals of neuroscience research is to determine how networks of neurons control and modify behavior. One of the most influential model systems for this kind of analysis is the siphon and gill withdrawal reflex of the marine mollusc A. californica. In response to tactile stimulation, the siphon displays 3 different responses: (1) a posterior pointing and leveling (flaring) of the siphon in response to tail stimulation (the siphon T response), (2) constriction and anterior pointing to head stimulation (the siphon H response) and (3) constriction and withdrawal between the animal's parapodia (the siphon S response). The siphon S response is pseudoconditioned by a noxious tail stimulus to resemble the siphon T response. Behavioral and combined behavioral/intracellular studies were conducted to determine the motor neuronal control of these behaviors and to search for mechanisms of siphon response transformation following pseudoconditioning. The present studies have found that the flaring component of pseudoconditioned siphon S responses occurs during mantle pumping (MP) triggered by noxious tail stimulation. Siphon stimulation also triggers MP, as recorded in neurons of the Interneuron II pattern generator which commands MP. The 4 LF$\rm\sb{SB}$ siphon motor neurons (SMNs) were found necessary and sufficient for the siphon T response, while SMNs RD$\rm\sb S$ and LD$\rm\sb{S1}$ were found necessary and sufficient for the siphon H response. Following pseudoconditioning, there is an increase in the number of evoked spikes to the test stimulus for the LF$\rm\sb{SB}$ cells and a decreased number for RD$\rm\sb S.$ Siphon flaring occurring during the pseudoconditioned response correlates with increased LF$\rm\sb{SB}$ activity during triggered MP cycles. This suggests that psuedoconditioning is in part due to reconfiguration of the motor outputs of the Interneuron II network. These results suggest that these defensive responses are controlled and patterned by a well-defined, finite set of motor neurons and interneurons (Interneuron II) that are dedicated to specific behavioral functions, but also have parallel distributed properties. ^

Quantitative MRI of spinal cord injury in a rat model: Correlative studies

Serial quantitative and correlative studies of experimental spinal cord injury (SCI) in rats were conducted using three-dimensional magnetic resonance imaging (MRI). Correlative measures included morphological histopathology, neurobehavioral measures of functional deficit, and biochemical assays for N-acetyl-aspartate (NAA), lactate, pyruvate, and ATP. A spinal cord injury device was characterized and provided a reproducible injury severity. Injuries were moderate and consistent to within $\pm$20% (standard deviation). For MRI, a three-dimensional implementation of the single spin-echo FATE (Fast optimum angle, short TE) pulse sequence was used for rapid acquisition, with a 128 x 128 x 32 (x,y,z) matrix size and a 0.21 x 0.21 x 1.5 mm resolution. These serial studies revealed a bimodal characteristic in the evolution in MRI pathology with time. Early and late phases of SCI pathology were clearly visualized in $T\sb2$-weighted MRI, and these corresponded to specific histopathological changes in the spinal cord. Centralized hypointense MRI regions correlated with evidence of hemorrhagic and necrotic tissue, while surrounding hyperintense regions represented edema or myelomalacia. Unexpectedly, $T\sb2$-weighted MRI pathology contrast at 24 hours after injury appeared to subside before peaking at 72 hours after injury. This change is likely attributable to ongoing secondary injury processes, which may alter local $T\sb2$ values or reduce the natural anisotropy of the spinal cord. MRI, functional, and histological measures all indicated that 72 hours after injury was the temporal maximum for quantitative measures of spinal cord pathology. Thereafter, significant improvement was seen only in neurobehavioral scores. Significant correlations were found between quantitated MRI pathology and histopathology. Also, NAA and lactate levels correlated with behavioral measures of the level of function deficit. Asymmetric (rostral/caudal) changes in NAA and lactate due to injury indicate that rostral and caudal segments from the injury site are affected differently by the injury. These studies indicate that volumetric quantitation of MRI pathology from $T\sb2$-weighted images may play an important role in early prediction of neurologic deficit and spinal cord pathology. The loss of $T\sb2$ contrast at 24 hours suggests MR may be able to detect certain delayed mechanisms of secondary injury which are not resolved by histopathology or other radiological modalities. Furthermore, in vivo proton magnetic resonance spectroscopy (MRS) studies of SCI may provide a valuable addition source of information about changes in regional spinal cord lactate and NAA levels, which are indicative of local metabolic and pathological changes. ^

Functional studies of homeobox gene {\it Cart1\/} during mouse development

Cart1 is a paired-class homeobox-containing gene that is expressed in head mesenchyme, branchial arches, limb buds, and various cartilages during embryogenesis. To understand the role of Cart1 during mammalian development, I generated Cart1-mutant mice by gene targeting in mouse embryonic stem cells. Cart1-homozygous mutants were born alive but all died soon after birth. Most had acrania (absence of the cranial vault) and meroanencephaly (absence of part of the brain). In situ hybridization studies showed that Cart1 is expressed specifically in forebrain mesenchyme but not in midbrain or hindbrain mesenchyme nor in the neural tube. Developmental studies revealed a transient deficiency of forebrain mesenchyme cells due to apoptosis associated with a delay in neural tube closure in that region. Subsequently, the forebrain region became filled with mesenchyme and closed, however, the midbrain neural tube region never initiated closure and remained open. These results suggest that Cart1 is required for the survival of forebrain mesenchyme and that its absence disrupts cranial neural tube morphogenesis by blocking the initiation of closure in the midbrain region, and this ultimately leads to the generation of lethal craniofacial defects. Prenatal treatment of Cart1 homozygous mutants with folic acid suppressed the development of the acrania/meroanencephaly phenotype. Thus, Cart1 mutant mice provide a novel animal model for understanding the cellular, molecular, and genetic etiology of neural tube defects and for the development of prenatal therapeutic protocols using folic acid. ^

Functional analysis of MRF4 during mouse embryogenesis

MRF4 is one of four skeletal muscle specific regulatory genes, (the other three genes being MyoD, myf5, and myogenin), each of which has the unique ability to orchestrate an entire program of muscle-specific transcription when introduced into diverse cell types. These findings have led to the notion that these factors function as master regulators of muscle cell fate. Analysis of mice lacking MyoD, myf5, and myogenin have further defined their roles in the commitment and differentiation of myotomal progenitor cells. Current data strongly supports the model that MyoD and myf5 share functional redundancy in determining the muscle cell lineage, while myogenin acts downstream of MyoD and myf5, to initiate myoblast differentiation. Unlike other myogenic bHLH genes, MRF4 is expressed predominantly in the adult, suggesting that it may function to regulate adult muscle maturation and maintenance. To test this hypothesis and to eventually incorporate MRF4 into a general model for muscle specification, differentiation, maturation and maintenance, I deleted the MRF4 gene. MRF4-null mice are viable and fertile, however, they show mild rib anomalies. In addition, the expression of myogenin is dramatically upregulated only in the adult, suggesting that myogenin may compensate for the loss of MRF4 in the adult, and MRF4 may normally suppress the expression of myogenin after birth. MRF4 is also required during muscle regeneration after injury.^ To determine the degree of genetic redundancy between MRF4-myogenin; and MRF4-MyoD, I crossed the MRF4-null mice with MyoD- and myogenin-null mice respectively. There are no additional muscle phenotypes in double-null progeny from a MRF4 and myogenin cross, suggesting that the existence of residual fibers in myogenin-null mice is not due to the presence of MRF4. MRF4 expression also cannot account for the ability of myogenin-null myoblasts to differentiate in vitro. However, the combination of the MRF4-null mutation with the myogenin-null mutation results in a novel rib phenotype. This result suggests that MRF4 modifies the myogenin-null rib phenotype, and MRF4 and myogenin play redundant roles in rib development.^ MRF4 also shares dosage effects with MyoD during mouse development. (MyoD+/$-$;MRF4$-$/$-$)mice are fertile and viable, while (MyoD$-$/$-$;MRF4+/$-$) mice die between birth and two weeks after birth, and have a small skeletal structure. The double homozygous mice for MRF4 and MyoD mutations are embryonic lethal and die at around E10.5. These results suggest that MRF4 and MyoD share overlapping functions during mouse embryogenesis. ^

Genetic and functional analyses of cell division proteins FtsZ and FtsA in Escherichia coli

In the current model for bacterial cell division, the FtsZ protein forms a ring that marks the division plane, creating a cytoskeletal framework for the subsequent action of other essential division proteins such as FtsA and ZipA. The putative protein complex ultimately generates the division septum. The essential cell division protein FtsZ is a functional and structural homolog of eukaryotic tubulin, and like tubulin, FtsZ hydrolyzes GTP and self-assembles into protein filaments in a strictly GTP-dependent manner. FtsA shares sequence similarity with members of the ATPase superfamily that include actin, but its actual function remains unknown. To test the division model and elucidate functions of the division proteins, this dissertation primarily focuses on the analysis of FtsZ and FtsA in Escherichia coli. ^ By tagging with green fluorescent protein, we first demonstrated that FtsA also exhibits a ring-like structure at the potential division site. The localization of FtsA was dependent on functional FtsZ, suggesting that FtsA is recruited to the septum by the FtsZ ring. In support of this idea, we showed that FtsA and FtsZ directly interact. Using a novel E. coli in situ assay, we found that the FtsA-FtsZ interaction appears to be species-specific, although an interspecies interaction could occur between FtsA and FtsZ proteins from two closely related organisms. In addition, mutagenesis of FtsA revealed that no single domain is solely responsible for its septal localization or interaction with FtsZ. To explore the function of FtsA, we purified FtsA protein and demonstrated that it has ATPase activity. Furthermore, purified FtsA stimulates disassembly of FtsZ polymers in a sedimentation assay but does not affect GTP hydrolysis of FtsZ. This result suggests that in the cell, FtsA may function similarly in regulating dynamic instability of the FtsZ ring during the cell division process. ^ To elucidate the structure-function relationship of FtsZ, we carried out thorough genetic and functional analyses of the mutagenized FtsZ derivatives. Our results indicate that the conserved N-terminal domain of FtsZ is necessary and sufficient for FtsZ self-assembly and localization. Moreover, we discovered a critical role for an extreme C-terminal domain of FtsZ that consists of only 12 residues. Truncated FtsZ derivatives lacking this domain, though able to polymerize and localize, are defective in ring formation in vivo as well as interaction with FtsA and ZipA. Alanine scanning mutagenesis of this region pinpointed at least five residues necessary for the function of FtsZ. Studies of protein levels and protein-protein interactions suggested that these residues may be involved in regulating protein stability and/or FtsZ-FtsA interactions. Interestingly, two of the point mutants exhibited dominant-negative phenotypes. ^ In summary, results from this thesis work have provided additional support for the division machinery model and will contribute to a better understanding of the coordinate functions of FtsA and FtsZ in the cell division process. ^

Gene regulation during placenta development: Murine adenosine deaminase as a model to study placenta specific expression

The formation of the placenta is one of the first and most important developmental events that occur in early mammalian embryogenesis. Even given this importance of the placenta, the academic community has largely ignored studying gene regulation during the development and maturation of the placenta. For this reason, an in-depth study of gene regulation in the trophoblast layer of the placenta using murine Adenosine Deaminase (Ada) as a model system has been undertaken. It has been determined that Ada is highly expressed in the placenta and is critical for embryo development. Dr. Kellems' laboratory has previously described a 1.8 kb fragment of the Ada 5 ′ flanking region that is capable of directing trophoblast specific expression in a transgenic model system. Preliminary studies have demonstrated several critical portions of this fragment that are necessary for the correct tissue specific expression in the placenta. My first specific aim was to elucidate the trans factor binding to one of these sequences, the FP3. Through electromobility shift assays (EMSA), the 30 bp FP3 was narrowed to a 5 bp sequence which computer databases predicted bound to Acute Myeloid Leukemia 1 (AML-1). This was confirmed by supershift analysis. The functional importance of this binding was demonstrated by a transgenic approach. A significant difference in expression of the reporter in the placenta was seen when the 5 bp sequence was mutated. This finding is a novel use for the AML-1 transcription factor which is the DNA binding portion of the heterodimer Core Binding Protein (CBP). The 5′ 240 bp region has also been demonstrated to contain functionally significant sequence. Through EMSA assays and computer predictions, the area has been narrowed to two pertinent regions that are predicted to contain GATA binding motifs. ^