Functional analysis of the Galpha protein, GNA -3, in the development of Neurospora crassa using biochemical and genetic studies

Extracellular signals regulate fungal development and, to sense and respond to these cues, fungi evolved signal transduction pathways similar to those in mammalian systems. In fungi, heterotrimeric G proteins, composed of α, β, and γ subunits, transduce many signals, such as pheromones and nutrients, intracellularly to alter adenylyl cyclase and MAPK cascades activity. ^ Previously, the Gα proteins GNA-1 and GNA-2 were characterized in regulating development in the fungus Neurospora crassa. R. A. Baasiri isolated a third Gα, gna-3, and P. S. Rowley generated Δgna-3 mutants. GNA-3 belongs to a fungal Gα family that regulates cAMP metabolism and virulence. The Δ gna-3 sexual cycle is defective in homozygous crosses, producing inviable spores. Δgna-3 mutants have reduced aerial hyphae formation and derepressed asexual sporulation (conidiation), causing accumulation of asexual spores (conidia). These defects are similar to an adenylyl cyclase mutant, cr-1; cAMP supplementation suppressed Δ gna-3 and cr-1. Inappropriate conidiation and expression of a conidiation gene, con-10, were higher in Δ gna-3 than cr-1 submerged cultures; peptone suppressed conidiation. Adenylyl cyclase activity and expression demonstrated that GNA-3 regulates enzyme levels. ^ A Δgna-1 cr-1 was analyzed with F. D. Ivey to differentiate GNA-1 roles in cAMP-dependent and -independent pathways. Δ gna-1 cr-1 defects were worse than cr-1 and refractory to cAMP, suggesting that GNA-1 is necessary for sensing extracellular CAMP. Submerged culture conidiation was highest in Δgna-1 cr-1, and only high cell density Δgna-1 cultures conidiated, which correlated with con-10 levels. Transcription of a putative heat shock cognate protein was highest in Δgna-1 cr-1. ^ Functional relationships between the three Gαs was analyzed by constructing Δgna-1 Δgna-2 Δ gna-3, Δgna-1 Δgna-3, and Δgna-2 Δgna-3 strains. Δ gna-2 Δgna-3 strains exhibited intensified Δ gna-3 phenotypes; Δgna-1 Δgna-2 Δgna-3 and Δgna-1 Δ gna-3 strains were identical to Δgna-1 cr-1 on plates and were non-responsive to cAMP. The highest levels of conidiation and con-10 were detected in submerged cultures of Δ gna-1 Δgna-2 Δgna-3 and Δgna-1 Δgna-3 mutants, which was partially suppressed by peptone supplementation. Stimulation of adenylyl cyclase is completely deficient in Δgna-1 Δ gna-2 Δgna-3 and Δgna-1 Δ gna-3 strains. Δgna-3 and Δ gna-1 Δgna-3 aerial hyphae and conidiation defects were suppressed by mutation of a PKA regulatory subunit. ^

Structural and functional analysis of p84N5

Normal development and tissue homeostasis requires the carefully orchestrated balance between cell proliferation and cell death. Cell cycle checkpoints control the extent of cell proliferation. Cell death is coordinated through the activation of a cell suicide pathway that results in the morphologically recognizable form of death, apoptosis. Tumorigenesis requires that the balance between these two pathways be disrupted. The tumor suppressor protein Rb has not only been shown to be involved in the enforcement of cell cycle checkpoints, but has also been implicated in playing a role in the regulation of apoptosis. The manner in which Rb enforces cell cycle checkpoints has been well studied; however, its involvement in the regulation of apoptosis is still very unclear. p84N5 is a novel nuclear death domain containing protein that has been shown to interact with the N-terminus of Rb. The fact that it contains a death domain and the fact that it is nuclear localized possibly provides the first known mechanism for apoptotic signaling from the nucleus. The following study tested the hypothesis that the novel exclusively nuclear death domain containing protein p84N5 is an important mediator of programmed cell death and that its apoptotic function is reliant upon its nuclear localization and is regulated by unique functional domains within the p84N5 protein. We identified the p84N5 nuclear localization signal (NLS), eliminated it, and tested the functional significance of nuclear localization by using wild type and mutant sequences fused to EGFP-C1 (Clontech) to create wild type GFPN5 and subsequent mutants. The results of these assays demonstrated exclusive nuclear localization of GFPN5 is required for normal p84N5 induced apoptosis. We further conducted large-scale mutagenesis of the GFPN5 construct to identify a minimal region within p84N5 capable of interacting with Rb. We were able to identify a minimal sequence containing p84N5 amino acids 318 to 464 that was capable of interacting with Rb in co-immunoprecipitation assays. We continued by conducting a structural and functional analysis to identify the region or regions within p84N5 responsible for inducing apoptosis. Point mutations and small-scale deletions within the death domain of p84N5 lessened the effect but did not eliminate p84N5-induced cytotoxicity. Further analysis revealed that the minimal sequence of 318 to 464 of p84N5 was capable of inducing apoptosis to a similar degree as wild-type GFPN5 protein. Since amino acids 318 to 464 of p84N5 are capable of inducing apoptosis and interacting with Rb, we propose possible mechanisms whereby p84N5 may function in a Rb regulated manner. These results demonstrate that p84N5 induced apoptosis is reliant upon its nuclear localization and is regulated by unique functional domains within the p84N5 protein. ^

Cloning and molecular analysis of paraxis: A {\it trans\/}-acting factor required for somite maturation

During vertebrate embryogenesis, cells from the paraxial mesoderm coalesce in a rostral-to-caudal progression to form the somites. Subsequent compartmentalization of the somites yields the sclerotome, myotome and dermatome, which give rise to the axial skeleton, axial musculature, and dermis, respectively. Recently, we cloned a novel basic-Helix-Loop-Helix (bHLH) protein, called scleraxis, which is expressed in the sclerotome, in mesenchymal precursors of bone and cartilage, and in connective tissues. This dissertation focuses on the cloning, expression and functional analysis of a bHLH protein termed paraxis, which is nearly identical to scleraxis within the bHLH region but diverges in both its amino and carboxyl termini. During the process of mouse embryogenesis, paraxis transcripts are first detected at about day 7.5 post coitum within the primitive mesoderm lying posterior to the head and heart primordia. Subsequently, paraxis expression progresses caudally through the paraxial mesoderm, immediately preceding somite formation. Paraxis is expressed at high levels in newly formed somites before the first detectable expression of the myogenic bHLH genes, and as the somite becomes compartmentalized, paraxis becomes downregulated within the myotome.^ To determine the function of paraxis during mammalian embryogenesis, mice were generated with a null mutation in the paraxis locus. Paraxis null mice survived until birth, but exhibited severe foreshortening along the anteroposterior axis due to the absence of vertebrae caudal to the midthoracic region. The phenotype also included axial skeletal defects, particularly shortened bifurcated ribs which were detached from the vertebral column, fused vertebrae and extensive truncation and disorganization caudal to the hindlimbs. Mutant neonates also lacked normal levels of trunk muscle and exhibited defects in the dermis as well as the stratification of the epidermis. Analysis of paraxis -/- mutant embryos has revealed a failure of the somites to both properly epithelialize and compartmentalize, resulting in defects in somite-derived cell lineages. These results suggest that paraxis is an essential component of the genetic pathway regulating somitogenesis. ^

CELL POLARITY REGULATES ORGAN GROWTH THROUGH THE HIPPO PATHWAY

Defects in apical-basal cell polarity and abnormal expression of cell polarity determinants are linked to human cancer. Loss of polarity is highly correlated with malignancy. In Drosophila, perturbation of apical-basal polarity, including overexpressing the apical determinant Crumbs, can lead to uncontrolled tissue growth. Cells mutant for the basolateral determinant scribble overproliferate and can form neoplastic tumors. Interestingly, scribble mutant clones that arise in wild-type tissues are eliminated and therefore do not manifest their tumorigenic potential. However, the mechanisms by which cell polarity coordinates with growth control pathways in developing organs to achieve appropriate organ size remain obscure. To investigate the function of apical determinants in growth regulation, I investigated the mechanism by which the apical determinant Crumbs affects growth in Drosophila imaginal discs. I found that crumbs gain and loss of function cause overgrowth and induction of Hippo target genes. In addition, Crumbs is required for the proper localization of Expanded, an upstream component of the Hippo pathway. Furthermore, we uncoupled the cell polarity and growth control function of Crb through structure-functional analysis. Taken together, our data identify a role of Crb in growth regulation specifically through modulation of the Hippo pathway. To further explore the role of polarity in growth control, I investigated how cells mutant for basolateral determinants are eliminated by using patches of cells mutant for scribble (scribble mutant clones) as a model system. We found that competitive cell-cell interactions eliminate tumorigenic scribble cells by modulation of the Hippo pathway. The regulation of Hippo signaling is required and sufficient to restrain the tumorous growth of scribble mutant cells. Artificially increasing the relative fitness of scribble mutant cells unleashes their tumorigenic potential. Therefore, we have identified a novel tumor-suppression mechanism that depends on signaling between normal and tumorigenic cells. These data identify evasion of cell competition as a critical step toward malignancy and illustrate a role for wild-type tissue in eliminating abnormal cells and preventing the formation of tumors.

ABERRATIONS OF A PUTATIVE TUMOR SUPPRESSOR GENE SEL1L IN PANCREATIC DUCTAL ADENOCARCINOMA

Introduction: Pancreatic cancer is the fourth leading cause of cancer-related death among males and females in the United States. Sel-1-like (SEL1L) is a putative tumor suppressor gene that is downregulated in a significant proportion of human pancreatic ductal adenocarcinoma (PDAC). It was hypothesized that SEL1L expression could be down-modulated by somatic mutation, loss of heterozygosity (LOH), CpG island hypermethylation and/or aberrantly expressed microRNAs (miRNAs). Material and methods: In 42 PDAC tumors, the SEL1L coding region was amplified using reverse transcription polymerase chain reaction (RT-PCR), and analyzed by agarose gel electrophoresis and sequenced to search for mutations. Using fluorescent fragment analysis, two intragenic microsatellites in the SEL1L gene region were examined to detect LOH in a total of 73 pairs of PDAC tumors and normal-appearing adjacent tissues. Bisulfite DNA sequencing was performed to determine the methylation status of the SEL1L promoter in 41 PDAC tumors and 6 PDAC cell lines. Using real-time quantitative PCR, the expression levels of SEL1L mRNA and 7 aberrantly upregulated miRNAs that potentially target SEL1L were assessed in 42 PDAC tumor and normal pairs. Statistical methods were applied to evaluate the correlation between SEL1L mRNA and the miRNAs. Further the interaction was determined by functional analysis using a molecular biological approach. Results: No mutations were detected in the SEL1L coding region. More than 50% of the samples displayed abnormally alternate or aberrant spliced transcripts of SEL1L. About 14.5% of the tumors displayed LOH at the CAR/CAL microsatellite locus and 10.7% at the RepIN20 microsatellite locus. However, the presence of LOH did not show significant association with SEL1L downregulation. No methylation was observed in the SEL1L promoter. Statistical analysis showed that SEL1L mRNA expression levels significantly and inversely correlated with the expression of hsa-mir-143, hsa-mir-155, and hsa-mir-223. Functional analysis indicated that hsa-mir-155 acted as a suppressor of SEL1L in PL18 and MDAPanc3 PDAC cell lines. Discussion: Evidence from these studies suggested that SEL1L was possibly downregulated by aberrantly upregulated miRNAs in PDAC. Future studies should be directed towards developing a better understanding of the mechanisms for generation of aberrant SEL1L transcripts, and further analysis of miRNAs that may downregulate SEL1L.

Identification of two vitamin D response elements in the rat osteopontin gene

Studies to elucidate the function of vitamin D have demonstrated an important role in regulating bone-related cells, including osteoblasts and osteoclasts. A seemingly paradoxical observation is that 1,25(OH)$\sb2$D$\sb3$, the active metabolite of vitamin D, stimulates bone resorption, yet regulates transcription of genes expressed by osteoblasts. One mechanism that could explain these actions is the upregulation of transcription of osteoblast-specific genes. These gene products could then act as effectors to influence osteoclastic activity. We hypothesized that molecular signals could be deposited directly into the mineralized matrix in the form of noncollagenous proteins, such as osteopontin (OPN). The structure, biosynthesis and localization of OPN suggest that it could function to mediate the molecular "cross talk" between osteoblasts and osteoclasts in response to 1,25(OH)$\sb2$D$\sb3$. To begin to address this hypothesis, elucidation of the molecular mechanisms of action involved in the transactivation of OPN by 1,25(OH)$\sb2$D$\sb3$ is essential.^ In the present study, the rat opn gene was isolated and characterized. Functional analysis by transient transfection of the 5$\sp\prime$ flanking sequences of the rat opn gene fused to the luciferase gene demonstrated that OPN is transcriptionally upregulated by 1,25(OH)$\sb2$D$\sb3$, mediated through two vitamin D response elements (VDRE). Both proximal and distal VDREs are structurally similar (two imperfect direct repeats separated by a 3 nucleotide spacer) and bind protein complexes that include the VDR and retinoid-X receptor (RXR). Isolated VDRE expression constructs produce functional activity of equivalent magnitude of responsiveness to 1,25(OH)$\sb2$D$\sb3$. However, expression constructs containing either VDRE and at least 200 bp of 5$\sp\prime$ and 3$\sp\prime$ flanking sequence demonstrated that the distal VDRE produces an amplitude of response significantly higher than the proximal VDRE. We conclude that the transcriptional upregulation of the opn gene by 1,25(OH)$\sb2$D$\sb3$ involves the transactivation of two VDREs, while maximal responsiveness requires interaction of the VDREs with additional cis-elements contained in the 5$\sp\prime$ sequence. ^

Structure-function analyses of {\it PAX6\/} and the molecular bases of aniridia

PAX6 is a transcription activator that regulates eye development in animals ranging from Drosophila to human. The C-terminal region of PAX6 is proline/serine/threonine-rich (PST) and functions as a potent transactivation domain when attached to a heterologous DNA-binding domain of the yeast transcription factor, GAL4. The PST region comprises 152 amino acids encoded by four exons. The transactivation function of the PST region has not been defined and characterized in detail by in vitro mutagenesis. I dissected the PST domain in two independent systems, a heterologous system using a GAL4 DNA-binding site and the native system of PAX6. In both systems, the results show consistently that all four constituent exons of the PST domain are responsible for the transactivation function. The four exon fragments act cooperatively to stimulate transcription, although none of them can function individually as an independent transactivation domain. Combinations of two or more exon fragments can reconstitute substantial transactivation activity when fused to the DNA-binding domain of GAL4, but they surprisingly do not produce much activity in the context of native PAX6 even though the mutant PAX6 proteins are stable and their DNA-binding function remains unaffected. I conclude that the PAX6 protein contains an unusually large transactivation domain that is evolutionarily conserved to a high degree, and that its full transactivation activity relies on the cooperative action of the four exon fragments.^ Most PAX6 mutations detected in patients with aniridia result in truncations of the protein. Some of the truncation mutations occur in the PST region of PAX6, resulting in mutant proteins that retain their DNA-binding ability but have no significant transactivation activity. It is not clear whether such mutants are true loss-of-function or dominant-negative mutants. I show that these mutants are dominant-negative if they are coexpressed with wild-type PAX6 in cultured cells and that the dominant-negative effects result from enhanced DNA-binding ability of these mutants due to removal of the PST domain. These mutants are able to repress the wild-type PAX6 activity not only at target genes with paired domain binding sites but also at target genes with homeodomain binding sites.^ Mutations in the human PAX6 gene produce various phenotypes, including aniridia, Peters' anomaly, autosomal dominant keratitis, and familial foveal dysplasia. The various phenotypes may arise from different mutations in the same gene. To test this theory, I performed a functional analysis of two missense mutations in the paired domain: the R26G mutation reported in a case of Peters' anomaly, and the I87R mutation identified in a patient with aniridia. While both the R26 and the I87 positions are conserved in the paired boxes of all known PAX genes, X-ray crystallography has shown that only R26 makes contact with DNA. I found that the R26G mutant failed to bind a subset of paired domain binding sites but, surprisingly, bound other sites and successfully transactivated promoters containing those sites. In contrast, the I87R mutant had lost the ability to bind DNA at all tested sites and failed to transactivate promoters. My data support the haploinsufficiency hypothesis of aniridia, and the hypothesis that R26G is a hypomorphic allele. ^

Transcriptional regulation of the human {\it PAX6\/} gene

PAX6, a member of the paired-type homeobox gene family, is expressed in a partially and temporally restricted pattern in the developing central nervous system, and its mutation is responsible for human aniridia (AN) and mouse small eye (Sey). The objective of this study was to characterize the PAX6 gene regulation at the transcriptional level, and thereby gain a better understanding of the molecular basis of the dynamic expression pattern and the diversified function of the human PAX6 gene.^ Initially, we examined the transcriptional regulation of the PAX6 gene by transient transfection assays and identified multiple cis-regulatory elements that function differently in different cell lines. The transcriptional initiation site was identified by RNase protection and primer extension assays. Examination of the genomic DNA sequence indicated that the PAX6 promoter has a TATA like-box (ATATTTT) at $-$26 bp, and two CCAAT-boxes are located at positions $-$70 and $-$100 bp. A 38 bp ply (CA) sequence was located 992 bp upstream from the initiation site. Transient transfection assays in glioblastoma cells and leukemia cells indicate that a 92 bp region was required for basal level PAX6 promoter activity. Gel retardation assays showed that this 92 bp sequence can form four DNA-protein complexes which can be specifically competed by a 31-mer oligonucleotide containing a PAX6 TATA-like sequence or an adenovirus TATA box. The activation of the promoter is positively correlated with the expression of PAX6 transcripts in cells tested.^ Based on the results obtained from the in vitro transfection assays, we did further dissection assay and functional analysis in both cell-culture and transgenic mice. We found that a 5 kb upstream promoter sequence is required for the tissue specific expression in the forebrain region which is consistent with that of the endogenous PAX6 gene. A 267 bp cell-type specific repressor located within the 5 kb fragment was identified and shown to direct forebrain specific expression. The cell-type specific repressor element has been narrowed to a 30 bp region which contains a consensus E-box by in vitro transfection assays. The third regulatory element identified was contained in a 162 bp sequence (+167 to +328) which functions as a midbrain repressor, and it appeared to be required for establishing the normal expression pattern of the PAX6 gene. Finally, a highly conserved 216 bp sequence identified in intron 4 exhibited as a spinal cord specific enhancer. And this 216 bp cis-regulatory element can be used as a marker to trace the differentiation and migration of progenitor cells in the developing spinal cord. These studies show that the concerted action of multiple cis-acting regulatory elements located upstream and downstream of the transcription initiation site determines the tissue specific expression of PAX6 gene. ^

Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

IDENTIFYING GENETIC VARIANTS AND CHARACTERIZING THEIR ROLE IN CLUBFOOT

Clubfoot is a common, complex birth defect affecting 4,000 newborns in the United States and 135,000 world-wide each year. The clubfoot deformity is characterized by inward and rigid downward displacement of one or both feet, along with persistent calf muscle hypoplasia. Despite strong evidence for a genetic liability, there is a limited understanding of the genetic and environmental factors contributing to the etiology of clubfoot. The studies described in this dissertation were performed to identify variants and/or genes associated with clubfoot. Genome-wide linkage scan performed on ten multiplex clubfoot families identified seven new chromosomal regions that provide new areas to search for clubfoot genes. Troponin C (TNNC2) the strongest candidate gene, located in 20q12-q13.11, is involved in muscle contraction. Exon sequencing of TNNC2 did not identify any novel coding variants. Interrogation of fifteen muscle contraction genes found strong associations with SNPs located in potential regulatory regions of TPM1 (rs4075583 and rs3805965), TPM2 (rs2025126 and rs2145925) and TNNC2 (rs383112 and rs437122). In previous studies, a strong association was found with rs3801776 located in the basal promoter of HOXA9, a gene also involved in muscle development and patterning. Altogether, this data suggests that SNPs located in potential regulatory regions of genes involved in muscle development and function could alter transcription factor binding leading to changes in gene expression. Functional analysis of 3801776/HOXA9, rs2025126/TPM2 and rs2145925/TPM2 showed altered protein binding, which significantly influenced promoter activity. Although the ancestral allele (G) of rs4075583/TPM1 creates a DNA-protein complex, it did not affect TPM1 promoter activity. However and importantly, in the context of a haplotype, rs4075583/G significantly decreased TPM1 promoter activity. These results suggest dysregulation of multiple skeletal muscle genes, TPM1, TPM2, TNNC2 and HOXA9, working in concert may contribute to clubfoot. However, specific allelic combinations involving these four regulatory SNPs did not confer a significantly higher risk for clubfoot. Other combinations of these variants are being evaluated. Moreover, these variants may interact with yet to be discovered variants in other genes to confer a higher clubfoot risk. Collectively, we show novel evidence for the role of skeletal muscle genes in clubfoot indicating that there are multiple genetic factors contributing to this complex birth defect.

{\it cis\/} -acting elements and {\it trans\/} -acting factors that control serum amyloid A gene expression during inflammation

To understand how the serum amyloid A (SAA) genes are regulated, the cis-acting elements and trans-acting factors involved in the regulation of mouse SAA3 and rat SAA1 genes expression during inflammation were analyzed.^ To identify DNA sequences involved in the liver-specific expression of the mouse SAA3 gene, the 5$\sp\prime$ flanking region of this gene was analyzed by transient transfection studies. Results suggest that C/EBP, a liver-enriched transcription factor, plays an important role for the enhanced expression of the mouse SAA3 gene in hepatocytes.^ Transfection studies of the regulation of the expression of rat SAA1 gene indicated that a 322 bp fragment ($-$304 to +18) of the gene contains sufficient information for cytokine-induced expression of the reporter gene in a liver cell-specific manner. Further functional analysis of the 5$\sp\prime$ flanking region of the rat SAA1 gene demonstrated that a 65 bp DNA fragment ($-$138/$-$73) can confer cytokine-inducibility onto a heterologous promoter both in liver and nonliver cells. DNase I footprint and gel retardation assays identified five putative cis-regulatory elements within the 5$\sp\prime$ flanking region of the gene: one inducible element, a NF$\kappa$B binding site and four constitutive elements. Two constitutive elements, footprint regions I and III, were identified as C/EBP binding sites with region III having over a 10-fold higher affinity for C/EBP binding than region I. Functional analysis of the cis-elements indicated that C/EBP(I) and C/EBP(III) confer liver cell-specific activation onto a heterologous promoter, while sequences corresponding to the NF$\kappa$B element and C/EBP(I) impart cytokine responsiveness onto the heterologous promoter. These results suggest that C/EBP(I) possesses two functions: liver-specific activation and cytokine responsiveness. The identification of two cytokine responsive elements (NF$\kappa$B and C/EBP(I)), and two liver-specific elements (C/EBP(I) and C/EBP(III)) implies that multiple cis-acting elements are involved in the regulation of the expression of the rat SAA1 gene. The tissue-specific and cytokine-induced expression of rat SAA1 gene is likely the result of the interactions of these cis-acting elements with their cognate trans-acting factors as well as the interplay between the different cis-acting elements and their binding factors. (Abstract shortened with permission of author.) ^

A novel DNA damage inducible inhibitor of p53

p53 is required for the maintenance of the genomic stability of cells. Mutations in the p53 tumor-suppressor gene occur in more than 50% of human cancers of diverse types. In addition, 70% of families with Li-Fraumeni syndrome have a germline mutation in p53, predisposing these individuals to multiple forms of cancer. In response to DNA damage, p53 becomes stabilized and activated. However the exact mechanism by which DNA damage signals the stabilization and activation of p53 still remains elusive. The biochemical activity of p53 that is required for tumor suppression, and presumably the cellular response to DNA damage, involves the ability of the protein to bind to specific DNA sequences and to function as a transcription factor. For the downstream targets, p53 transactivates many genes involved in growth arrest, apoptosis and DNA repair such as p21, Bax and GADD45, respectively. An open question in the field is how cells can determine the downstream effects of p53. ^ We hypothesize that, through its associated proteins, p53 can differentially transactivate its target genes, which determine its downstream effect. Additionally, p53 interacting proteins may be involved in signaling for the stabilization and activation of p53. Therefore, a key aspect to understanding p53 function is the identification and analysis of proteins that interact with it. We have employed the Sos recruitment system (SRS), a cytoplasmic yeast two-hybrid screen to identify p53 interacting proteins. The SRS is based on the ability of Sos to activate Ras when it becomes localized to the plasma membrane. The system takes advantage of an S. cerevisiae strain, cdc25-2 temperature sensitive mutant, harboring a mutation in Sos. In this strain, fusion proteins containing a truncated Sos will only localize to the membrane by protein-protein interaction, which allows growth at non-permissive temperature. This system allows the use of intact transcriptional activators such as p53. ^ To date, using a modified SRS library screen to identify p53 interacting proteins, I have identified p53 (known to interact with itself) and a novel p53-interacting protein (PIP). PIP is a specific p53 interacting protein in the SRS. The interaction of p53 and PIP was further confirmed by performing in vitro and in vivo binding assays. In the in vivo binding study, the interaction can only be detected in the presence of ionizing radiation suggesting that this interaction might be involved in DNA-damage induced p53-signalling pathway. After screening cDNA and genomic libraries, a full-length PIP-cDNA clone ( ∼ 3kb) was obtained which encodes a protein of 429 amino acids with calculated molecular weight of 46 kDa. The results of genebank search indicated that the PIP is an unidentified gene and contains a conserved ring-finger domain, which is present in a diverse family of regulatory proteins involved in different aspects of cellular function. Northern blot analysis revealed that the size of its messenge is approximately 3 kb preferentially expressed in brain, heart, liver and kidney. The PIP protein is mainly located in the cytoplasm as determined by the cellular localization of a green fluorescence fusion protein. Preliminary functional analysis revealed that PIP downregulated the transactivation activity of p53 on both p21 and mdm2 promoters. Thus, PIP may be a novel negative regulator of p53 subsequent to DNA damage. ^

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

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

NONLINEAR TIME SERIES/SYSTEM ANALYSIS (STATE-SPACE, DYNAMICS, INTERVENTION, RESILIENCE, KALMAN)

A discussion of nonlinear dynamics, demonstrated by the familiar automobile, is followed by the development of a systematic method of analysis of a possibly nonlinear time series using difference equations in the general state-space format. This format allows recursive state-dependent parameter estimation after each observation thereby revealing the dynamics inherent in the system in combination with random external perturbations.^ The one-step ahead prediction errors at each time period, transformed to have constant variance, and the estimated parametric sequences provide the information to (1) formally test whether time series observations y(,t) are some linear function of random errors (ELEM)(,s), for some t and s, or whether the series would more appropriately be described by a nonlinear model such as bilinear, exponential, threshold, etc., (2) formally test whether a statistically significant change has occurred in structure/level either historically or as it occurs, (3) forecast nonlinear system with a new and innovative (but very old numerical) technique utilizing rational functions to extrapolate individual parameters as smooth functions of time which are then combined to obtain the forecast of y and (4) suggest a measure of resilience, i.e. how much perturbation a structure/level can tolerate, whether internal or external to the system, and remain statistically unchanged. Although similar to one-step control, this provides a less rigid way to think about changes affecting social systems.^ Applications consisting of the analysis of some familiar and some simulated series demonstrate the procedure. Empirical results suggest that this state-space or modified augmented Kalman filter may provide interesting ways to identify particular kinds of nonlinearities as they occur in structural change via the state trajectory.^ A computational flow-chart detailing computations and software input and output is provided in the body of the text. IBM Advanced BASIC program listings to accomplish most of the analysis are provided in the appendix. ^