Structure-function analysis of human Integrator subunit-4

Structure-function analysis of human Integrator subunit 4 Anupama Sataluri Advisor: Eric. J. Wagner, Ph.D. Uridine-rich small nuclear RNAs (U snRNA) are RNA Polymerase-II (RNAPII) transcripts that are ubiquitously expressed and are known to be essential for gene expression. snRNAs play a key role in mRNA splicing and in histone mRNA expression. Inaccurate snRNA biosynthesis can lead to diseases related to defective splicing and histone mRNA expression. Although the 3′ end formation mechanism and processing machinery of other RNAPII transcripts such as mRNA has been well studied, the mechanism of snRNA 3′ end processing has remained a mystery until the recent discovery of the machinery that mediates this process. In 2005, a complex of 14 subunits (the Integrator complex) associated with RNA Polymerase-II was discovered. The 14subunits were annotated Integrator 1-14 based on their size. The subunits of this complex together were found to facilitate 3′ end processing of snRNA. Identification of the Integrator complex propelled research in the direction of understanding the events of snRNA 3’end processing. Recent studies from our lab confirmed that Integrator subunit (IntS) 9 and 11 together perform the endonucleolytic cleavage of the nascent snRNA 3′ end to generate mature snRNA. However, the role of other members of the Integrator complex remains elusive. Current research in our lab is focused on deciphering the role of each subunit within the Integrator complex This work specifically focuses on elucidating the role of human Integrator subunit 4 (IntS4) and understanding how it facilitates the overall function of the complex. IntS4 has structural similarity with a protein called “Symplekin”, which is part of the mRNA 3’end processing machinery. Symplekin has been thoroughly researched in recent years and structure-function correlation studies in the context of mRNA 3’end processing have reported a scaffold function for Symplekin due to the presence of HEAT repeat motifs in its N-terminus. Based upon the structural similarity between IntS4 and Symplekin, we hypothesized that Integrator subunit 4 may be behaving as a Symplekin-like scaffold molecule that facilitates the interaction between other members of the Integrator Complex. To answer this question, the two important goals of this study were to: 1) identify the region of IntS4, which is important for snRNA 3′ end processing and 2) determine binding partners of IntS4 which promote its function as a scaffold. IntS4 structurally consists of a highly conserved N-terminus with 8 HEAT repeats, followed by a nonconserved C- terminus. A series of siRNA resistant N and C-terminus deletion constructs as well as specific point mutants within its N-terminal HEAT repeats were generated for human IntS4 and, utilizing a snRNA transcriptional readthrough GFP-reporter assay, we tested their ability to rescue misprocessing. This assay revealed a possible scaffold like property of IntS4. To probe IntS4 for interaction partners, we performed co-immunoprecipitation on nuclear extracts of IntS4 expressing stable cell lines and identified IntS3 and IntS5 among other Integrator subunits to be binding partners which facilitate the scaffold like function of hIntS4. These findings have established a critical role for IntS4 in snRNA 3′ end processing, identified that both its N and C termini are essential for its function, and mapped putative interaction domains with other Integrator subunits.

INVESTIGATING THE ROLES OF THE p63 ISOFORMS IN THE MICRORNA BIOGENESIS PATHWAY

MicroRNAs play roles in various biological processes like development, tumorigenesis, metastasis and pluripotency. My thesis work has demonstrated roles for p63, a p53 family member, in the upstream regulation of microRNA biogenesis. The p63 gene has a complex gene structure and has multiple isoforms. The TAp63 isoforms contain an acidic transcription activation domain. The ΔNp63 isoforms, lack the TA domain, but have a proline rich region critical for gene transactivation. To understand the functions of these isoforms, the Flores lab generated TAp63 and ΔNp63 conditional knock out mice. Using these mice and tissues and cells from these mice we have found that TAp63 transcriptionally regulates Dicer while ΔNp63 transcriptionally regulates DGCR8. TAp63 -/- mice are highly tumor prone. These mice develop metastatic mammary adenocarcinomas, squamous cell carcinomas, and lung adenocarcinomas to distant sites including the liver, lungs, and brain. I found that TAp63 suppresses metastasis by transcriptionally activating Dicer. TAp63 and Dicer levels were very low or lost in high grade human tumors like mammary adenocarcinomas, squamous cell carcinomas, and lung adenocarcinomas. Expression of Dicer in these tumor cell lines reduced their invasiveness. Using ΔNp63 -/- mice, I found that ΔNp63 transcriptionally activates DGCR8, resulting in a miRNA profile that is critical to reprogram cells to pluripotency. Analysis of epidermal cells derived from ΔNp63 -/- mice revealed that these cells expressed markers of pluripotency, including Sox2, Oct 4 and Nanog; however, genome-wide analysis revealed a novel profile of genes that are common between ΔNp63 -/- epidermal cells and embryonic stem cells. I also found that mouse cells depleted of ΔNp63 form chimeric mice and teratomas in SCID mice, demonstrating that ΔNp63 deficient cells are pluripotent. Further, I found that restoration of DGCR8 in ΔNp63 -/- epidermal cells reduces their pluripotency and induces terminal differentiation. I also demonstrated that iMS (induced multipotent stem) cells could be generated using human keratinocytes by knockdown of ∆Np63 or DGCR8. Taken together, my work has placed p63 and its isoforms at a critical node in controlling miRNA biogenesis.

Targeted Deletion Of Functionally Validated Enhancers Defines Their Role in Mouse Limb Development

Transcriptional enhancers are genomic DNA sequences that contain clustered transcription factor (TF) binding sites. When combinations of TFs bind to enhancer sequences they act together with basal transcriptional machinery to regulate the timing, location and quantity of gene transcription. Elucidating the genetic mechanisms responsible for differential gene expression, including the role of enhancers, during embryological and postnatal development is essential to an understanding of evolutionary processes and disease etiology. Numerous methods are in use to identify and characterize enhancers. Several high-throughput methods generate large datasets of enhancer sequences with putative roles in embryonic development. However, few enhancers have been deleted from the genome to determine their roles in the development of specific structures, such as the limb. Manipulation of enhancers at their endogenous loci, such as the deletion of such elements, leads to a better understanding of the regulatory interactions, rules and complexities that contribute to faithful and variant gene transcription – the molecular genetic substrate of evolution and disease. To understand the endogenous roles of two distinct enhancers known to be active in the mouse embryo limb bud we deleted them from the mouse genome. I hypothesized that deletion of these enhancers would lead to aberrant limb development. The enhancers were selected because of their association with p300, a protein associated with active transcription, and because the human enhancer sequences drive distinct lacZ expression patterns in limb buds of embryonic day (E) 11.5 transgenic mice. To confirm that the orthologous mouse enhancers, mouse 280 and 1442 (M280 and M1442, respectively), regulate expression in the developing limb we generated stable transgenic lines, and examined lacZ expression. In M280-lacZ mice, expression was detected in E11.5 fore- and hindlimbs in a region that corresponds to digits II-IV. M1442-lacZ mice exhibited lacZ expression in posterior and anterior margins of the fore- and hindlimbs that overlapped with digits I and V and several wrist bones. We generated mice lacking the M280 and M1442 enhancers by gene targeting. Intercrosses between M280 -/+ and M1442 -/+, respectively, generated M280 and M1442 null mice, which are born at expected Mendelian ratios and manifest no gross limb malformations. Quantitative real-time PCR of mutant E11.5 limb buds indicated that significant changes in transcriptional output of enhancer-proximal genes accompanied the deletion of both M280 and M1442. In neonatal null mice we observed that all limb bones are present in their expected positions, an observation also confirmed by histology of E18.5 distal limbs. Fine-scale measurement of E18.5 digit bone lengths found no differences between mutant and control embryos. Furthermore, when the developmental progression of cartilaginous elements was analyzed in M280 and M1442 embryos from E13.5-E15.5, transient development defects were not detected. These results demonstrate that M280 and M1442 are not required for mouse limb development. Though M280 is not required for embryonic limb development it is required for the development and/or maintenance of body size – adult M280 mice are significantly smaller than control littermates. These studies highlight the importance of experiments that manipulate enhancers in situ to understand their contribution to development.

Genome -wide and locus -specific instability

The discovery of expanded simple repeated sequences causing or associated with human disease has lead to a new area of research involved in the elucidation of how the expanded repeat causes disease and how the repeat becomes unstable. ^ To study the genetic basis of the (CTG)n repeat instability in the DMPK gene in myotonic dystrophy (DM1) patients, somatic cell hybrids were constructed between the lymphocytes of DM1 patients and a variety of Chinese hamster ovary (CHO) cell DNA repair gene deficient mutants. By using small pool PCR (SP-PCR), the instability of the (CTG)n can be quantitated for both the frequency and sizes of length change mutations. ^ Additional SP-PCR analysis on 2/11 subclones generated from this original hybrid showed a marked increase in large repeat deletions, ∼50%. A bimodal distribution of repeats was seen around the progenitor allele and at a large deleted product (within the normal range) with no intermediate products present. ^ To determine if the repair capacity of the CHO cell led to a mutator phenotype in the hamster and hybrid clones, SP-PCR was also done on 3 hamster microsatellites in a variety of hamster cell backgrounds. No variant alleles were seen in over 2500 genome equivalents screened. ^ Human-hamster hybrids have long been shown to be chromosomally unstable, yet information about the stability of repeated sequences was not known. To test if repeat instability was associated with either intact or non-intact human chromosomes, more than 300 microsatellite repeats on 13 human chromosomes (intact and non-intact) were analyzed in eight hybrid cells. No variants were seen between the hybrid and patient alleles in the hybrids. ^ To identify whether DM1 patients have a previously undetected level of genome wide instability or if the instability is truly locus specific, SP-PCR was done on 6 human microsatellites within the patient used to make the hybrid cells. No variants were seen in over 1000 genomes screened. ^ These studies show that the somatic cell hybrid approach is a genetically stable system that allows for the determination of factors that could lead to changes in microsatellite instability. It also shows that there is something inherent about the DM1 expanded (CTG)n repeat that it is solely targeted by, as of yet, and unknown mechanism that causes the repeat to be unstable. (Abstract shortened by UMI.)^

Activation of c -Met contributes to growth and metastasis of human colorectal carcinoma

c-Met is the protein tyrosine kinase receptor for hepatocyte growth factor/scatter factor (HGF/SF) and mediates several normal cellular functions including proliferation, survival, and migration. Overexpression of c-Met correlates with progression and metastasis of human colorectal carcinoma (CRC). The goals of this study were to determine if overexpression of c-Met directly contributes to tumorigenicity and liver metastatic potential of colon cancer, and what are the critical pathways regulated by c-Met in this process. The studies used two colon tumor cell lines, KM12SM and KM20, which express high levels of constitutively active c-Met and are highly metastatic in nude mice. To examine the effects of c-Met overexpression, subclones of theses lines with reduced c-Met expression were obtained following transfection with a c-Met specific targeting ribozyme. Reduction of c-Met in KM12SM cells abolished liver metastases when cells were injected intrasplenically in an experimental metastasis assay. However, c-Met downregulation in theses clones was unstable. Three stable KM20 clones with a 25–35% reduction in c-Met protein levels but 60–90% reduction in basal c-Met autophosphorylation and kinase activity were obtained. While HGF increased c-Met kinase activity in the clones with reduced c-Met, the activity was less than that observed in parental or control transfected cells. Correlating with the reduction in c-Met kinase activity, subclones with reduced c-Met expression had significantly reduced in vitro growth rates, soft-agar colony forming abilities, and increased apoptosis. HGF/SF treatment did not affect anchorage-dependent growth or soft-agar colony forming abilities. Further, c-Met downregulation significantly impaired the ability of HGF/SF to induce migration. To examine the effects of reduced c-Met on tumor formation, parental and c-Met reduced KM20 cells were grown subcutaneously and intrahepatically in nude mice. c-Met downregulation delayed, but did not abolish growth at the subcutaneous site. When these cells were injected intrahepatically, both tumor incidences and size were significantly reduced. To further understand the molecular basis of c-Met in promoting tumor growth, the activation of several signaling intermediates that have been implicated in c-Met mediated growth, survival and migration were compared between KM20 parental cells and subclones with reduced c-Met expression levels. The expression and activity (as determined by phosphorylation) of AKT and Erk1/2 were unaltered. In contrast, Src kinase activity, as measured by immune complex kinase assay, was reduced 2–5 fold following c-Met downregulation. As Src has been implicated in growth, survival and migration, Src activation in c-Met overexpressing lines is likely contributing to the tumorigenic and metastatic capabilities of colon tumor cell lines that overexpress c-Met. Collectively, these results suggest that c-Met overexpression plays a causal role in the development of CRC liver metastases, and that c-Src and c-Met inhibitors may be of potential therapeutic benefit for late-stage colon cancer. ^

Identification of 5′-deoxyribose phosphate lyase activity in human DNA polymerase γ and its role in mitochondrial base excision repair in vitro

Mitochondria have been proposed to possess base excision repair processes to correct oxidative damage to the mitochondrial genome. As the only DNA polymerase (pol) present in mitochondria, pol γ is necessarily implicated in such processes. Therefore, we tested the ability of the catalytic subunit of human pol γ to participate in uracil-provoked base excision repair reconstituted in vitro with purified components. Subsequent to actions of uracil-DNA glycosylase and apurinic/apyrimidinic endonuclease, human pol γ was able to fill a single nucleotide gap in the presence of a 5′ terminal deoxyribose phosphate (dRP) flap. We report here that the catalytic subunit of human pol γ catalyzes release of the dRP residue from incised apurinic/apyrimidinic sites to produce a substrate for DNA ligase. The heat sensitivity of this activity suggests the dRP lyase function requires a three-dimensional protein structure. The dRP lyase activity does not require divalent metal ions, and the ability to trap covalent enzyme-DNA complexes with NaBH4 strongly implicates a Schiff base intermediate in a β-elimination reaction mechanism.

Natural genetic exchange between Haemophilus and Neisseria: Intergeneric transfer of chromosomal genes between major human pathogens

Members of the bacterial families Haemophilus and Neisseria, important human pathogens that commonly colonize the nasopharynx, are naturally competent for DNA uptake from their environment. In each genus this process is discriminant in favor of its own and against foreign DNA through sequence specificity of DNA receptors. The Haemophilus DNA uptake apparatus binds a 29-bp oligonucleotide domain containing a highly conserved 9-bp core sequence, whereas the neisserial apparatus binds a 10-bp motif. Each motif (“uptake sequence”, US) is highly over-represented in the chromosome of the corresponding genus, particularly concentrated with core sequences in inverted pairs forming gene terminators. Two Haemophilus core USs were unexpectedly found forming the terminator of sodC in Neisseria meningitidis (meningococcus), and sequence analysis strongly suggests that this virulence gene, located next to IS1106, arose through horizontal transfer from Haemophilus. By using USs as search strings in a computer-based analysis of genome sequence, it was established that while USs of the “wrong” genus do not occur commonly in Neisseria or Haemophilus, where they do they are highly likely to flag domains of chromosomal DNA that have been transferred from Haemophilus. Three independent domains of Haemophilus-like DNA were found in the meningococcal chromosome, associated respectively with the virulence gene sodC, the bio gene cluster, and an unidentified orf. This report identifies intergenerically transferred DNA and its source in bacteria, and further identifies transformation with heterologous chromosomal DNA as a way of establishing potentially important chromosomal mosaicism in these pathogenic bacteria.

Expanding the functional human mitochondrial DNA database by the establishment of primate xenomitochondrial cybrids

The nuclear and mitochondrial genomes coevolve to optimize approximately 100 different interactions necessary for an efficient ATP-generating system. This coevolution led to a species-specific compatibility between these genomes. We introduced mitochondrial DNA (mtDNA) from different primates into mtDNA-less human cells and selected for growth of cells with a functional oxidative phosphorylation system. mtDNA from common chimpanzee, pigmy chimpanzee, and gorilla were able to restore oxidative phosphorylation in the context of a human nuclear background, whereas mtDNA from orangutan, and species representative of Old-World monkeys, New-World monkeys, and lemurs were not. Oxygen consumption, a sensitive index of respiratory function, showed that mtDNA from chimpanzee, pigmy chimpanzee, and gorilla replaced the human mtDNA and restored respiration to essentially normal levels. Mitochondrial protein synthesis was also unaltered in successful “xenomitochondrial cybrids.” The abrupt failure of mtDNA from primate species that diverged from humans as recently as 8–18 million years ago to functionally replace human mtDNA suggests the presence of one or a few mutations affecting critical nuclear–mitochondrial genome interactions between these species. These cellular systems provide a demonstration of intergenus mtDNA transfer, expand more than 20-fold the number of mtDNA polymorphisms that can be analyzed in a human nuclear background, and provide a novel model for the study of nuclear–mitochondrial interactions.

ETS target genes: Identification of Egr1 as a target by RNA differential display and whole genome PCR techniques

ETS transcription factors play important roles in hematopoiesis, angiogenesis, and organogenesis during murine development. The ETS genes also have a role in neoplasia, for example in Ewing’s sarcomas and retrovirally induced cancers. The ETS genes encode transcription factors that bind to specific DNA sequences and activate transcription of various cellular and viral genes. To isolate novel ETS target genes, we used two approaches. In the first approach, we isolated genes by the RNA differential display technique. Previously, we have shown that the overexpression of ETS1 and ETS2 genes effects transformation of NIH 3T3 cells and specific transformants produce high levels of the ETS proteins. To isolate ETS1 and ETS2 responsive genes in these transformed cells, we prepared RNA from ETS1, ETS2 transformants, and normal NIH 3T3 cell lines and converted it into cDNA. This cDNA was amplified by PCR and displayed on sequencing gels. The differentially displayed bands were subcloned into plasmid vectors. By Northern blot analysis, several clones showed differential patterns of mRNA expression in the NIH 3T3-, ETS1-, and ETS2-expressing cell lines. Sixteen clones were analyzed by DNA sequence analysis, and 13 of them appeared to be unique because their DNA sequences did not match with any of the known genes present in the gene bank. Three known genes were found to be identical to the CArG box binding factor, phospholipase A2-activating protein, and early growth response 1 (Egr1) genes. In the second approach, to isolate ETS target promoters directly, we performed ETS1 binding with MboI-cleaved genomic DNA in the presence of a specific mAb followed by whole genome PCR. The immune complex-bound ETS binding sites containing DNA fragments were amplified and subcloned into pBluescript and subjected to DNA sequence and computer analysis. We found that, of a large number of clones isolated, 43 represented unique sequences not previously identified. Three clones turned out to contain regulatory sequences derived from human serglycin, preproapolipoprotein C II, and Egr1 genes. The ETS binding sites derived from these three regulatory sequences showed specific binding with recombinant ETS proteins. Of interest, Egr1 was identified by both of these techniques, suggesting strongly that it is indeed an ETS target gene.

CD68+ cells of monocyte/macrophage lineage in the environment of AIDS-associated and classic-sporadic Kaposi sarcoma are singly or doubly infected with human herpesviruses 7 and 6B

Earlier studies have shown that Kaposi sarcomas contain cells infected with human herpesvirus (HHV) 6B, and in current studies we report that both AIDS-associated and classic-sporadic Kaposi sarcoma contain HHV-7 genome sequences detectable by PCR. To determine the distribution of HHV-7-infected cells relative to those infected with HHV-6, sections from paraffin-embedded tissues were allowed to react with antibodies to HHV-7 virion tegument phosphoprotein pp85 and to HHV-6B protein p101. The antibodies are specific for HHV-7 and HHV-6B, respectively, and they retained reactivity for antigens contained in formalin-fixed, paraffin-embedded tissue samples. We report that (i) HHV-7 pp85 was present in 9 of 32 AIDS-associated Kaposi sarcomas, and in 1 of 7 classical-sporadic HIV-negative Kaposi sarcomas; (ii) HHV-7 pp85 was detected primarily in cells bearing the CD68 marker characteristic of the monocyte/macrophage lineage present in or surrounding the Kaposi sarcoma lesions; and (iii) in a number of Kaposi sarcoma specimens, tumor-associated CD68+ monocytes/macrophages expressed simultaneously antigens from both HHV-7 and HHV-6B, and therefore appeared to be doubly infected with the two viruses. CD68+ monocytes/macrophages infected with HHV-7 were readily detectable in Kaposi sarcoma, but virtually absent from other normal or pathological tissues that harbor macrophages. Because all of the available data indicate that HHV-7 infects CD4+ T lymphocytes, these results suggest that the environment of the Kaposi sarcoma (i) attracts circulating peripheral lymphocytes and monocytes, triggers the replication of latent viruses, and thereby increases the local concentration of viruses, (ii) renders CD68+ monocytes/macrophages susceptible to infection with HHV-7, and (iii) the combination of both events enables double infections of cells with both HHV-6B and HHV-7.

Octamer-based genome scanning distinguishes a unique subpopulation of Escherichia coli O157:H7 strains in cattle

Multilocus-genotyping methods have shown that Escherichia coli O157:H7 is a geographically disseminated clone. However, high-resolution methods such as pulse-field gel electrophoresis demonstrate significant genomic diversity among different isolates. To assess the genetic relationship of human and bovine isolates of E. coli O157:H7 in detail, we have developed an octamer-based genome-scanning methodology, which compares the distance between over-represented, strand-biased octamers that occur in the genome. Comparison of octamer-based genome-scanning products derived from >1 megabase of the genome demonstrated the existence of two distinct lineages of E. coli O157:H7 that are disseminated within the United States. Human and bovine isolates are nonrandomly distributed among the lineages, suggesting that one of these lineages may be less virulent for humans or may not be efficiently transmitted to humans from bovine sources. Restriction fragment length polymorphism analysis with lambdoid phage genomes indicates that phage-mediated events are associated with divergence of the lineages, thereby providing one explanation for the degree of diversity that is observed among E. coli O157:H7 by other molecular-fingerprinting methods.

The chitinase PfCHT1 from the human malaria parasite Plasmodium falciparum lacks proenzyme and chitin-binding domains and displays unique substrate preferences

Within hours after the ingestion of a blood meal, the mosquito midgut epithelium synthesizes a chitinous sac, the peritrophic matrix. Plasmodium ookinetes traverse the peritrophic matrix while escaping the mosquito midgut. Chitinases (EC 3.2.1.14) are critical for parasite invasion of the midgut: the presence of the chitinase inhibitor, allosamidin, in an infectious blood meal prevents oocyst development. A chitinase gene, PgCHT1, recently has been identified in the avian malaria parasite P. gallinaceum. We used the sequence of PgCHT1 to identify a P. falciparum chitinase gene, PfCHT1, in the P. falciparum genome database. PfCHT1 differs from PgCHT1 in that the P. falciparum gene lacks proenzyme and chitin-binding domains. PfCHT1 was expressed as an active recombinant enzyme in Escherichia coli. PfCHT1 shares with PgCHT1 a substrate preference unique to Plasmodium chitinases: the enzymes cleave tri- and tetramers of GlcNAc from penta- and hexameric oligomers and are unable to cleave smaller native chitin oligosaccharides. The pH activity profile of PfCHT1 and its IC50 (40 nM) to allosamidin are distinct from endochitinase activities secreted by P. gallinaceum ookinetes. Homology modeling predicts that PgCHT1 has a novel pocket in the catalytic active site that PfCHT1 lacks, which may explain the differential sensitivity of PfCHT1 and PgCHT1 to allosamidin. PfCHT1 may be the ortholog of a second, as yet unidentified, chitinase gene of P. gallinaceum. These results may allow us to develop novel strategies of blocking human malaria transmission based on interfering with P. falciparum chitinase.

Derepression of human embryonic ζ-globin promoter by a locus-control region sequence

A multiple protein–DNA complex formed at a human α-globin locus-specific regulatory element, HS-40, confers appropriate developmental expression pattern on human embryonic ζ-globin promoter activity in humans and transgenic mice. We show here that introduction of a 1-bp mutation in an NF-E2/AP1 sequence motif converts HS-40 into an erythroid-specific locus-control region. Cis-linkage with this locus-control region, in contrast to the wild-type HS-40, allows erythroid lineage-specific derepression of the silenced human ζ-globin promoter in fetal and adult transgenic mice. Furthermore, ζ-globin promoter activities in adult mice increase in proportion to the number of integrated DNA fragments even at 19 copies/genome. The mutant HS-40 in conjunction with human ζ-globin promoter thus can be used to direct position-independent and copy number-dependent expression of transgenes in adult erythroid cells. The data also supports a model in which competitive DNA binding of different members of the NF-E2/AP1 transcription factor family modulates the developmental stage specificity of an erythroid enhancer. Feasibility to reswitch on embryonic/fetal globin genes through the manipulation of nuclear factor binding at a single regulatory DNA motif is discussed.

Cluster analysis and display of genome-wide expression patterns

A system of cluster analysis for genome-wide expression data from DNA microarray hybridization is described that uses standard statistical algorithms to arrange genes according to similarity in pattern of gene expression. The output is displayed graphically, conveying the clustering and the underlying expression data simultaneously in a form intuitive for biologists. We have found in the budding yeast Saccharomyces cerevisiae that clustering gene expression data groups together efficiently genes of known similar function, and we find a similar tendency in human data. Thus patterns seen in genome-wide expression experiments can be interpreted as indications of the status of cellular processes. Also, coexpression of genes of known function with poorly characterized or novel genes may provide a simple means of gaining leads to the functions of many genes for which information is not available currently.

Genome scan for teratogen-induced clefting susceptibility loci in the mouse: Evidence of both allelic and locus heterogeneity distinguishing cleft lip and cleft palate

Nonsyndromic clefting of the lip and palate in humans has a highly complex etiology, with both multiple genetic loci and exposure to teratogens influencing susceptibility. Previous studies using mouse models have examined only very small portions of the genome. Here we report the findings of a genome-wide search for susceptibility genes for teratogen-induced clefting in the AXB and BXA set of recombinant inbred mouse strains. We compare results obtained using phenytoin (which induces cleft lip) and 6-aminonicotinamide (which induces cleft palate). We use a new statistical approach based on logistic regression suitable for these categorical data to identify several chromosomal regions as possible locations of clefting susceptibility loci, and we review candidate genes located within each region. Because cleft lip and cleft palate do not frequently co-aggregate in human families and because these structures arise semi-independently during development, these disorders are usually considered to be distinct in etiology. Our data, however, implicate several of the same chromosomal regions for both forms of clefting when teratogen-induced. Furthermore, different parental strain alleles are usually associated with clefting of the lip versus that of the palate (i.e., allelic heterogeneity). Because several other chromosomal regions are associated with only one form of clefting, locus heterogeneity also appears to be involved. Our findings in this mouse model suggest several priority areas for evaluation in human epidemiological studies.