Control of chromosomal rearrangements in {\it Escherichia coli\/}

A complete physical map of Escherichia coli K-12 strain MG1655 was constructed by digesting chromosomal DNA with the infrequently cutting restriction enzymes NotI, SfiI and XbaI and separating the fragments by pulsed field gel electrophoresis. The map was used to compare six K-12 strains of E. coli. Although several differences were noted and localized, the map of MG1655 was representative of all the K-12 strains tested. The maps were also used to analyze chromosomal rearrangements in the E. coli strain MG1655. The spontaneous and UV induced frequencies of tandem duplication formation were measured at several loci distributed around the chromosome. The spontaneous duplication frequency varied from 10$\sp{-5}$ to 10$\sp{-3}$ and increased at least ten-fold following mild UV irradiation treatment. Duplications of several regions of the chromosome, including the serA region and the metE region, were mapped using pulsed field gel electrophoresis. Duplications of serA were found to be large, ranging in size from 600 kb to 2100 kb. Several of the duplications isolated at serA were caused by ectopic recombination between IS5 elements and between IS186 elements. Duplications of the metE region, however, were almost exclusively the result of ectopic recombination between ribosomal RNA cistrons. Duplication frequencies were determined at both serA and metE in wild type and mismatch repair mutant strains (mutL, mutS, uvrD and recF). Even though all of the mismatch repair mutations increased duplication frequency of metE, the largest increases were observed in the mutL and mutS strains. Duplication frequency of serA was increased less dramatically by mutations in mismatch repair. Several duplications of metE isolated in a wild type and a mismatch repair mutant were mapped. The results showed that the same repeated sequences were used for duplication formation in the mismatch repair mutant as were used in the wild type strain. Several isolates showed evidence of multiple rearrangements indicating that mismatch repair may play a role in stabilizing the genome by controlling chromosomal rearrangement. ^

Mutagenicity of herpes simplex virus type 1 in a shuttle vector

The shuttle vector plasmid pZ189 was used to find the kinds of mutations that are induced by herpes simplex virus type-1 (HSV-1). In cells infected by HSV-1 the frequency of mutation in supF gene, the mutagenesis marker, was increased over background by from two- to seven-fold, reaching 0.14-0.45%. No increase was induced by infection by vaccinia virus under the same conditions. Mutagenesis was an early event, showing a four-fold increase in mutation frequency at only two hours after infection, and peaking at a seven-fold increase at four hours after infection. DNA sequencing and gel electrophoresis analysis were performed on 105 HSV-1 induced mutants and 65 spontaneous mutants and provided the following information: (1) A change in plasmid size was seen in 54% of HSV-1 related mutants, compared with only 37% of spontaneous mutants. (2) Among point mutations, the predominant type was G:C to A:T transition, which accounted for 51% of point mutations in mutants isolated from cells infected with HSV-1, and 32% of point mutations in spontaneous mutants. (3) Deletions of DNA were seen in HSV-1 related mutants at a frequency of 40%, compared with 29% in spontaneous mutants. The HSV-1 related deletions were about half the length of spontaneous mutants and three contained short filler sequences. (4) Fifteen (15%) of HSV-1 induced mutants revealed the altered restriction patterns on agarose gel electrophoresis analysis and were due either to rearrangements of plasmid DNA, and/or to insertion of sequences derived from chromosomal DNA (seven plasmids). No insertions of DNA from HSV-1 were detected. Among spontaneous mutants, only 5 (7.7%) were rearrangements and none had inserted chromosomal DNA. (5) DNA sequence analysis of seven plasmids with inserted chromosomal DNA revealed that four cases had repetitive DNA sequences integrated and the other three were unidentified sequences from the GenBank database. Three repetitive DNA included $\alpha$ satellite, Alu and KpnI family sequences. The other sequence was identified as tRNA-like component. The observed mutations have implications for the mechanism of malignant transformation of cells by HSV-1. ^

Regulation of Spec genes in aboral ectoderm cells of sea urchin {\it Strongylocentrotus purpuratus}

The Spec genes of the sea urchin Stronylocentrotus purpuratus serves as an excellent model for studying cell type-specific gene expression during early embryogenesis. The Spec1/Spec2 genes encode cytosolic calcium-binding proteins related to the calmodulin/troponin C/myosin light chain superfamily. Members of the Spec gene family are activated shortly after the sixth cleavage as the lineage-specific founder cells giving rise to aboral ectoderm are established, and the accumulation of the Spec mRNAs is limited exclusively to aboral ectoderm cell lineages. In this dissertation, the transcriptional regulation of the Spec genes was studied. Sequence comparisons of the Spec gene 5$\sp\prime$ flanking regions showed that a DNA block of approximately 800 bp from the 3$\sp\prime$ end of the first exon to the 5$\sp\prime$ end of a repetitive DNA element, termed RSR, was highly conserved. In Spec2a, the conserved region was a continuous stretch of DNA, but in Spec1 and Spec2c, DNA insertions interrupt the conserved sequence block and alter the relative placement of the RSR element and other 5$\sp\prime$ flanking DNA. Thus, drastic rearrangements have occurred within the putative control regions of the Spec genes. In vivo expression experiments using the sea urchin embryo gene-transfer system showed that while the 5$\sp\prime$ flanking regions of all three Spec genes conferred proper temporal activation to the reporter CAT gene, only the Spec2a 5$\sp\prime$ flanking region could restrict lacZ gene expression to aboral ectoderm cells. However, the Spec2a conserved region alone was not sufficient to confer proper spatial expression, suggesting that negative spatial elements are also associated with the proper activation of Spec2a. A major positive regulatory region, defined as the RSR enhancer, was identified between base pairs $-$631 and $-$443 on Spec2a. The RSR enhancer was essential for maximal activity and conferred preferential aboral ectoderm expression to a lacZ reporter gene. DNaseI footprinting and band-shift analysis of the RSR enhancer revealed multiple DNA-elements. One of the elements, an A/T-rich sequence called the A/T palindrome was studied in detail. This element binds a single 45-kDa nuclear protein, the A/T palindrome binding protein (A/TBP), whose DNA-binding specificity suggests a possible relationship with the bicoid-class homeodomain proteins. Mutated A/T palindromes are incapable of binding the 45-kDa protein and lower promoter activity by 8-fold. DNA-binding activity for A/TBP is low in unfertilized eggs, increases by the 16-cell stage and continues rising in blastulae. These data suggest that A/TBP plays a major role in the activation of the Spec2a gene in aboral ectoderm cells. ^

The physical map of {\it Escherichia coli\/} K-12 strain MG1655: Construction, comparison with other strains and implications for mechanisms controlling oligonucleotide frequency

Complete NotI, SfiI, XbaI and BlnI cleavage maps of Escherichia coli K-12 strain MG1655 were constructed. Techniques used included: CHEF pulsed field gel electrophoresis; transposon mutagenesis; fragment hybridization to the ordered $\lambda$ library of Kohara et al.; fragment and cosmid hybridization to Southern blots; correlation of fragments and cleavage sites with EcoMap, a sequence-modified version of the genomic restriction map of Kohara et al.; and correlation of cleavage sites with DNA sequence databases. In all, 105 restriction sites were mapped and correlated with the EcoMap coordinate system.^ NotI, SfiI, XbaI and BlnI restriction patterns of five commonly used E. coli K-12 strains were compared to those of MG1655. The variability between strains, some of which are separated by numerous steps of mutagenic treatment, is readily detectable by pulsed-field gel electrophoresis. A model is presented to account for the difference between the strains on the basis of simple insertions, deletions, and in one case an inversion. Insertions and deletions ranged in size from 1 kb to 86 kb. Several of the larger features have previously been characterized and some of the smaller rearrangements can potentially account for previously reported genetic features of these strains.^ Some aspects of the frequency and distribution of NotI, SfiI, XbaI and BlnI cleavage sites were analyzed using a method based on Markov chain theory. Overlaps of Dam and Dcm methylase sites with XbaI and SfiI cleavage sites were examined. The one XbaI-Dam overlap in the database is in accord with the expected frequency of this overlap. The occurrence of certain types of SfiI-Dcm overlaps are overrepresented. Of the four subtypes of SfiI-Dcm overlap, only one has a partial inhibitory effect on the activity of SfiI. Recognition sites for all four enzymes are rarer than expected based on oligonucleotide frequency data, with this effect being much stronger for XbaI and BlnI than for NotI and SfiI. The latter two enzyme sites are rare mainly due to apparent negative selection against GGCC (both) and CGGCCG (NotI). The former two enzyme sites are rare mainly due to effects of the VSP repair system on certain di-tri- and tetranucleotides, most notably CTAG. Models are proposed to explain several of the anomalies of oligonucleotide distribution in E. coli, and the biological significance of the systems that produce these anomalies is discussed. ^

Defining the regions of homology between human chromosome 16p12--13 and mouse chromosomes

In this study, the evolutionary relationship between human chromosome 16p12-p13 and mouse chromosomes was investigated by determining the order of marker loci in the region and then identifying the chromosomal locations of the homologous loci in mice. Eighteen genes from human 16 were mapped to fifteen subchromosomal regions by a variety of mapping approaches.^ Thirteen of the genes were mapped in the mouse. Linkage analysis with backcross mice and segregation analysis in a mouse - Chinese Hamster Ovary (CHO) somatic cell hybrid panel informative for different regions of mouse genome were used. The results assigned the thirteen genes to three different mouse chromosomes.^ A group of six genes on mouse 16 was found to be closely linked to Scid. The order of Myh11 and Mrp remains ambiguous since no recombination was detected in backcross analysis. Their relative position in human is also uncertain since they were shown to be very close to each other. For the other mouse loci, an unambiguous gene order could be determined and was found to be identical to that in human. Therefore, they comprise a new conserved linkage group between the two species. The orientation of the group was inverted relative to the centromeres, i.e. the proximal loci in one species become distal in another. The size of the group was estimated to be from 4.4 to 8 Mb and 10 to 32 cM in human. In mouse, it was about 21 cM in the backcross analysis. The two boundaries of the conserved linkage were defined within a 1 Mb range. It is now possible to predict the locations of mouse homologs for some human disease genes based on their locations on human 16p.^ The six human 16p genes that map to MMU7 showed a different gene order in mouse than in human. No recombination was found between Crym and Umod while Crym was distal to D16S79A and proximal to D16S92. The location of Stp and Cdr2 with respect to the above four loci was not determined since they were not mapped in the same set of backcross mice. These genes greatly expanded an existing conserved synteny group between the human 16p12-p13 region and the MMU7. It now consists of eleven loci that span a region of probably more than 10 Mb in human. The gene order derived from this study provided further evidence for chromosomal rearrangements within the conserved synteny. (Abstract shortened by UMI.) ^

Molecular analysis of the human $\gamma$-chain T cell receptor genes

In our studies we have focused on the issue of variability and diversity of the $\gamma$ (or $\delta)$ chain T cell receptor (TCR) genes by studying cDNA transcripts in peripheral blood mononuclear cells or $\gamma\delta$ TCR+ T cell clones. The significance of these studies lies in the better understanding of the molecular biology of the $\gamma\delta$ T cell receptor as well as in answering the question whether certain molecular forms predominate in $\gamma\delta$ T cells exhibiting specific immunologic functions. We establish that certain $\gamma$-chain TCR genes exhibit particular patterns of rearrangements in cDNA transcripts in normal individuals. V$\gamma$I subgroup were shown to preferentially rearrange to J$\gamma$2C$\gamma$2 gene segments. These preferential VJC rearrangements, may have implications regarding the potential for diversity and polymorphism of the $\gamma$-chain TCR gene. In addition, the preferential association of V$\gamma$I genes with J$\gamma$2C$\gamma$2, which encode a non-disulfide-linked $\gamma\delta$ TCR, suggests that $\gamma$ chains utilizing V$\gamma$I are predominantly expressed as non-disulfide-linked $\gamma\delta$ TCR heterodimers. The implications of this type of expression remain to be determined. We identified two alternative splicing events of the $\gamma$-chain TCR genes occurring in high frequency in all the normal individuals examined. These events may suggest additional mechanisms of regulation and control as well as diversification of $\gamma\delta$ TCR gene expression. The question whether particular forms of $\gamma$ or $\delta$-chain TCR genes are involved in HLA Class I recognition by specific $\gamma\delta$ cytotoxic T cell clones was addressed. Our results indicated that the T cell clones expressed identical $\gamma$ but distinct $\delta$-chains suggesting that the specificity for recognition of HLA-A2 or HLA-A3 may be conferred by the $\delta$-chain TCR. The issue of the degree of diversity and polymorphism of the $\delta$-chain TCR genes in a patient with a primary immunodeficiency (Omenn's syndrome) was addressed. A limited pattern of rearrangements in peripheral blood transcripts was found, suggesting that a limited $\gamma\delta$ TCR repertoire may be expressed in this particular primary immunodeficiency syndrome. Overall, our findings suggest that $\delta$-chain TCR genes exhibit the potential for significant diversity and that there are certain preferential patterns of expression that may be associated with particular immunologic functions. ^

GENETIC EVENTS PREDISPOSING TO WILMS' TUMOR

Wilms' tumor (WT) is a childhood embryonic tumor of the kidney. In some cases, WT has been associated with a chromosome deletion in the region 11p13. The majority of WT cases, however, have normal karyotypes with no discernable deletions or rearrangements of chromosome 11.^ To study the genetic events predisposing to the development of WT, I have used a number of gene markers specific for chromosome 11. Gene probes for human catalase and apolipoprotein A1 were localized to chromosome 11 by in situ hybridization. A number of other probes previously mapped to chromosome 11 were also used. Nine WT patients who were heterozygous for at least one 11p marker were shown to lose heterozygosity in their tumor DNA. Gene dosage experiments demonstrated that two chromosomes 11 were present although loss of heterozygosity had occurred in all but two cases. By using gene probes from the short and long arms of chromosome 11, I discerned that loss of heterozygosity was due to somatic recombination in four cases, chromosome deletion in two cases, and chromosome loss and reduplication or somatic recombination in these cases. Examination of DNAs from the parents of six of these patients indicated that the alleles that were lost in tumor tissues were alleles inherited from the mother. In sporadic WT cases one would expect the loss of alleles to be random. These data suggest that the loss of alleles resulting in the development of WT is not a random event, however, the significance of this is not known. ^

GCN5 maintains genome stability during development

The histone acetyltransferase, GCN5, is essential for survival of mice during embryogenesis. GCN5 null embryos die early during development due to increased apoptosis. We have demonstrated that the increased apoptosis in associated with increased p53 protein levels. Loss of p53 rescues the embryonic apoptosis in the GCN5 null embryos. These results raised the question of what molecular trigger leads to p53 stabilization and cell death in the absence of GCN5. p53 is generally referred to as the gatekeeper of the cell, monitoring cellular responses to DNA damage, genotoxic stress, and other unfavorable conditions in the cell. Therefore, we examined individual cells in wild type and mutant embryos for gross chromosomal aberrations that might trigger a genome integrity checkpoint. Karyotype analysis indicates that approximately 30% of the cells in an E8.5 GCN5 null embryo display chromosomal aberrations, predominantly chromosomal end adhesions and associations. In wild type E8.5 embryos, only 6% of the cells have chromosomal aberrations. Recent data using telomeric FISH demonstrates that cells from GCN5 null embryos have a decreased telomeric signal. Telomere maintenance is essential for maintaining genome integrity. Telomeric defects are associated with loss of chromosomes and chromosomal rearrangements that can lead to detrimental gene fusions involved in many types of cancers. Little is known about the chromatin structures present near the telomeric ends, or whether any of the telomere-associated proteins are subject to post-translational modification such as acetylation. Our results are the first data to demonstrate the involvement of a histone acetyltransferase, GCN5, in maintaining genome integrity through telomere maintenance and/or capping. ^

The roles of UPF3a and UPF3b in nonsense mediated decay

Nonsense-mediated decay (NMD) degrades aberrant transcripts containing premature termination codons (PTCs). The T-cell receptor (TCR) locus undergoes error-prone rearrangements that frequently acquire PTCs. Transcripts harboring PTCs from this locus are downregulated much more than transcripts from non-rearranging genes. Efficient splicing is essential for this robust downregulation. ^ Here I show that TCR NMD is unique in another respect: it is not impaired by RNAi-mediated depletion of the NMD factor UPF3b. This differentiates TCR transcripts from classical NMD (assayed using β-globin or triose phosphate isomerase transcripts), which does depend on UPF3b. Depletion of UPF3a, which encodes a gene related to UPF3b, also had no effect on TCR NMD. Mapping experiments identified TCR sequences that when deleted or mutated caused a switch to UPF3b dependence. Since UPF3b dependence was invariably accompanied by less efficient RNA splicing, this suggests that UPF3b-dependent NMD occurs when transcripts are generated by inefficient splicing. Microarray analysis revealed the existence of many NMD-targeted mRNAs from wild-type genes whose downregulation is impervious to UPF3b depletion. This suggests the existence of an alternative NMD pathway independent of UPF3b that is widely used to downregulate the level of both normal and mutant transcripts. ^ During the course of my studies, I also found that the function of UPF3a is fundamentally distinct from that of UPF3b in several aspects. First, classical NMD failed to be impaired by UPF3a depletion, whereas it was reversed by UPF3b depletion. Second, UPF3a depletion had no effect on NMD elicited by tethered UPF2, whereas UPF3b depletion blocked this response. Thus, UPF3a does not function in classical NMD. Third, UPF3b depletion upregulated the expression of UPF3a, whereas UPF3a depletion had no effect on UPF3b expression. This suggests that a UPF3b-mediated feedback network exists that regulates the UPF3a expression. Lastly, UPF3a depletion but not UPF3b depletion significantly upregulated TCR precursor RNAs. This suggests that UPF3a, not UPF3b, functions in the surveillance of precursor RNAs, which typically contain many PTCs in the introns. Collectively, my data suggests that UPF3a and UPF3b are not functionally redundant, as previously thought, but instead have separable functions. ^

Whole genome analysis of Fusobacterium nucleatum, Rickettsia typhi and Francisella tularensis

The genomes of Fusobacterium nucleatum subspecies polymorphum strain ATCC 10953, Rickettsia typhi strain Wilmington, and Francisella tularensis subspecies holarctica strain OSU18 were sequenced, annotated, and analyzed. Each genome was then compared to the sequenced genomes of closely related bacteria. The genome of F. nucleatum ATCC 10953 was compared to two additional F. nucleatum subspecies, subspecies nucleatum and subspecies vincentii. This analysis revealed substantial evidence of horizontal gene transfer along with considerable genetic diversity within the species of F. nucleatum. R. typhi was compared to R. prowazekii and R. conorii. This analysis uncovered a hotspot for chromosomal rearrangements in the Spotted Fever Group but not the Typhus Group Rickettsia and revealed the close genetic relationship between the Typhus Group rickettsial species. F. tularensis OSU18 was compared to two additional F. tularensis strains. These comparisons uncovered significant chromosomal rearrangements between F. tularensis subspecies due to recombination between insertion sequence elements. ^

Functional differences of beta-tubulin isotypes: A study of microtubule assembly and antimitotic drug resistance

Mammalian cells express 7 β-tubulin isotypes in a tissue specific manner. This has long fueled the speculation that different isotypes carry out different functions. To provide direct evidence for their functional significance, class III, IVa, and VI β-tubulin cDNAs were cloned into a tetracycline regulated expression vector and stably transfected Chinese hamster ovary cell lines expressing different levels of ectopic β-tubulin were compared for effects on microtubule organization, microtubule assembly and sensitivity to antimitotic drugs. It was found that all three isotypes coassembled with endogenous β-tubulin. βVI expression caused distinct microtubule rearrangements including microtubule dissociation from the centrosome and accumulation at the cell periphery; whereas expression of βIII and βVIa caused no observable changes in the interphase microtubule network. Overexpression of all 3 isotypes caused spindle malformation and mitotic defects. Both βIII and βIVa disrupted microtubule assembly in proportion to their abundance and thereby conferred supersensitivity to microtubule depolymerizing drugs. In contrast, βVI stabilized microtubules at low stoichiometry and thus conferred resistance to many microtubule destabilizing drugs but not vinblastine. The 3 isotypes caused differing responses to microtubule stabilizing drugs. Expression of βIII conferred paclitaxel resistance while βVI did not. Low expression of βIVa caused supersensitivity to paclitaxel, whereas higher expression resulted in the loss of supersensitivity. The results suggest that βIVa may possess an enhanced ability to bind paclitaxel that increases sensitivity to the drug and acts substoichiometrically. At high levels of βVIa expression, however, microtubule disruptive effects counteract the assembly promoting pressure exerted by increased paclitaxel binding, and drug supersensitivity is lost. From this study, I concluded that β-tubulin isotypes behave differently from each other in terms of microtubule organization, microtubule assembly and dynamics, and antimitotic drug sensitivity. The isotype composition of cell can impart subtle to dramatic effects on the properties of microtubules leading to potential functional consequences and opening the opportunity to exploit differences in microtubule isotype composition for therapeutic gain. ^

Emerging roles for the double strand break repair protein 53BP1 in transcriptional regulation

Disruption of the mechanisms that regulate cell-cycle checkpoints, DNA repair, and apoptosis results in genomic instability and often leads to the development of cancer. In response to double stranded breaks (DSBs) as induced by ionizing radiation (IR), generated during DNA replication, or through immunoglobulin heavy chain (IgH) rearrangements in T and B cells of lymphoid origin, the protein kinases ATM and ATR are central players that activate signaling pathways leading to DSB repair. p53 binding protein 1 (53BP1) participates in the repair of DNA double stranded breaks (DSBs) where it is recruited to or near sites of DNA damage. In addition to its well established role in DSB repair, multiple lines of evidence implicate 53BP1 in transcription which stem from its initial discovery as a p53 binding protein in a yeast two-hybrid screen. However, the mechanisms behind the role of 53BP1 in these processes are not well understood. ^ 53BP1 possesses several motifs that are likely important for its role in DSB repair including two BRCA1 C-terminal repeats, tandem Tudor domains, and a variety of phosphorylation sites. In addition to these motifs, we identified a glycine and arginine rich region (GAR) upstream of the Tudor domains, a sequence that is oftentimes serves as a site for protein arginine methylation. The focus of this project was to characterize the methylation of 53BP1 and to evaluate how methylation influenced the role of 53BP1 as a tumor suppressor. ^ Using a variety of biochemical techniques, we demonstrated that 53BP1 is methylated by the PRMT1 methyltransferase in vivo. Moreover, GAR methylation occurs on arginine residues in an asymmetric manner. We further show that sequences upstream of the Tudor domains that do not include the GAR stretch are sufficient for 53BP1 oligomerization in vivo. While investigating the role of arginine methylation in 53BP1 function, we discovered that 53BP1 associates with proteins of the general transcription apparatus as well as to other factors implicated in coordinating transcription with chromatin function. Collectively, these data support a role for 53BP1 in regulating transcription and provide insight into the possible mechanisms by which this occurs. ^

Mechanisms by which nonsense codons downregulate T-cell receptor transcripts

Nonsense-mediated mRNA decay (NMD) is a quality control mechanism that degrades aberrant mRNAs harboring premature termination codons (PTCs). Two out of three T-cell receptor β (TCRβ) transcripts carry PTCs as a result of error-prone programmed rearrangements that occur at this locus during lymphocyte maturation. PTCs decrease TCRβ mRNA levels to a much greater extent than mRNAs transcribed from non-rearranging genes. This robust decrease in TCRβ mRNA levels is not a unique characteristic of the T-cell environment or the TCRβ promoter. The simplest explanation for this is that PTC-bearing TCRβ mRNAs elicit a stronger NMD response. An alternative explanation is NMD collaborates with another mechanism to dramatically decrease PTC-bearing TCRβ mRNA levels. ^ In my dissertation, I investigated the molecular mechanism behind the strong decrease in TCRβ mRNA levels triggered by PTCs. To determine the location of this response, I performed mRNA half-life analysis and found that PTCs elicited more rapid TCRβ mRNA decay in the nuclear fraction, not the cytoplasmic fraction. Although decay was restricted to the nuclear fraction, PTC-bearing TCRβ transcript levels were extremely low in the cytoplasm, a phenomenon that I named the nonsense-codon induced partitioning shift (NIPS). I established that NIPS shares several qualities with NMD, including its dependence on translation and NMD factors. Several lines of evidence suggested that NIPS results from PTCs eliciting retention of TCRβ transcripts in the nuclear fraction. This retention, as well as rapid TCRβ mRNA decay, most likely occurs in either the nucleoplasm or the outer nuclear membrane, based on analysis of nuclear and cytoplasmic markers in the highly purified nuclei I used for my studies. To further address the location of decay, I asked whether nuclear or cytoplasmic RNA decay factors mediated the destruction of PTC-bearing mRNAs. My results suggested that a nuclear component of the 3'-to-5' exosome, as well as an endonucleolytic activity, are involved in the destruction of PTC-containing TCRβ mRNAs. Individual endogenous NMD substrates had differential requirements for nuclear and cytoplasmic exonucleases. In summary, my results provide evidence that PTCs trigger multiple mechanisms involving multiple decay factors to remove and regulate mRNAs in mammalian cells. ^

INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

NOVEL INDUCERS OF GLIOTOXIN PRODUCTION IN ASPERGILLUS FUMIGATUS

Secondary metabolites are produced by numerous organisms and can either be benign to humans or harmful. Genes involved in the synthesis and transport of these secondary metabolites are frequently found in gene clusters, which are often located in subtelomeric regions of the chromosome. These clusters are often coordinately regulated, being almost exclusively dependent on transcription factors that are located within the clusters themselves. Secondary metabolites are also regulated by a variety of factors, including nutritional factors, environmental factors and developmental processes. Gliotoxin, which is produced by a variety of Aspergillus species, Trichoderma species, and Penicillium species, exhibits immunosuppressive properties and has therefore been the subject of research for many laboratories. There have been a few proteins shown to regulate the gliotoxin cluster, most notably GliZ, a Zn2Cys6 binuclear finger transcription factor that lies within the cluster, and LaeA, a putative methyltransferase that globally regulates secondary metabolism clusters within numerous fungal organisms, although no study has demonstrated the direct binding of any protein to a promoter region in the gliotoxin cluster. I report here two novel proteins, GipA, a C2H2 transcription factor and GipB, a hybrid sensor kinase, which are involved in regulating the gliotoxin biosynthetic cluster. GipA plays an important role in gliotoxin production, as high-copy expression of gipA induces gliotoxin biosynthesis and loss of gipA reduces gliotoxin biosynthesis by 50%. GipB is also involved in regulating gliotoxin production, as high-copy expression of gipB induces gliotoxin biosynthesis, but only during certain stages of asexual development. Furthermore, loss of gipB reduces gliotoxin biosynthesis by 10%. Based on data obtained from this project, I propose a model for the regulation of gliA, the efflux pump of the gliotoxin cluster, which involves GipB signaling through both GliZ and GipA. I propose that GliZ and GipA are interdependent, as mutation of the GipA DNA binding site in the gliA promoter negatively affects both GliZ-mediated and GipA-mediated induction of gliA. This is further supported by the fact that GliZ cannot fully induce gliA in the absence of GipA and vice versa. This is the first time that anyone has shown evidence of a protein directly binding to the gliotoxin cluster. Even though biosynthetic clusters are often coordinately regulated, my model raises the possibility that gliA is independently regulated, as the layout of the binding site in the gliA promoter is not present upstream of any other genes in the gliotoxin cluster, except for gliZ.