Published papers

Mutations in the EDA gene cause anhidrotic ectodermal dysplasia (EDA), with lesions in skin appendage formation. To begin to analyze EDA pathways, we have used expression profiling on 15 000-gene mouse cDNA microarrays, comparing adult mouse skin from wild-type, EDA-defective (Tabby) mice, and Tabby mice supplemented with the EDA-A1 isoform, which is sufficient to rescue multiple Tabby phenotypes. Given the sensitivity of the current microarray system, 8500 genes (60%) were estimated to be expressed, including transcription factors and growth-regulatory genes that had not previously been identified in skin; but only 24 (0.16%), one-third of them novel, showed significant differences between wild type and Tabby. An additional eight genes not included in the 15 000 gene set were shown to have expression differences by real-time RT-PCR. Sixteen of 32 affected genes were restored significantly toward wild-type levels in EDA-A1 transgenic Tabby mice. Significant up-regulation in Tabby skin was observed for several dermal matrix genes, including Col1a1, Col1a2, Col3a1 and Sparc. In contrast, down-regulation occurred for the NEMO/NF-kB pathway, already implicated in skin appendage formation, and even more markedly for a second pathway, JNK/c-jun/c-fos and their target genes, that has not previously been clearly associated with skin development. These data are consistent with the regulation of the NF-kB pathway by EDA, and support its involvement in the regulation of the JNK pathway as well.

Many clustering algorithms have been used to analyze microarray gene expression data. Given embryonic stem cell gene expression data, we applied several indices to evaluate the performance of clustering algorithms, including hierarchical clustering, k-means, PAM and SOM. The indices were homogeneity and separation scores, silhouette width, redundant score (based on redundant genes), and WADP (testing the robustness of clustering results after small perturbation). The results showed that the ES cell dataset posed a challenge for cluster analysis in that the clusters generated by different methods were only partially consistent. Using this data set, we were able to evaluate the advantages and weaknesses of algorithms with respect to both internal and external quality measures. This study may provide a guideline on how to select suitable clustering algorithms and it may help raise issues in the extraction of meaningful biological information from microarray expression data.

Mammals achieve gene dosage control by (1) random X-chromosome inactivation in females, (2) parental origin-specific imprinting of selected autosomal genes, and (3) random autosomal inactivation. Genes belonging to the third category of epigenetic phenomenon are just now emerging, with only six identified so far. Here we report three additional genes, Nubp2, Igfals, and Jsap1, that show 50%-methylated CpG sites by Southern blot analyses and primarily monoallelic expression in single-cell allele-specific RT-PCR analysis of bone marrow stromal cells and hepatocytes. Furthermore, we show that, in contrast to X inactivation, alleles can switch between active and inactive states during the formation of daughter cells. These three genes are the first in their category to exist as a tight cluster, in the proximal region of mouse chromosome 17, providing a thus far unique example of a region of autosomal random monoallelic expression.

A comprehensive gene-based map of a genome is a powerful tool for genetic studies and is especially useful for the positional cloning and positional candidate approaches. The availability of gene maps for multiple organisms provides the foundation for detailed conserved-orthology maps showing the correspondence between conserved genomic segments. These maps make it possible to use cross-species information in gene hunts and shed light on the evolutionary forces that shape the genome. Here we report a radiation hybrid map of mouse genes, a combined project of the Whitehead Institute/Massachusetts Institute of Technology Center for Genome Research, the Medical Research Council UK Mouse Genome Centre, and the National Center for Biotechnology Information. The map contains 11,109 genes, screened against the T31 RH panel and positioned relative to a reference map containing 2,280 mouse genetic markers. It includes 3,658 genes homologous to the human genome sequence and provides a framework for overlaying the human genome sequence to the mouse and for sequencing the mouse genome.

Here we report a novel design of linker primer that allows one to differentially amplify long tracts (average 3.0 kb with size ranges of 1-7 kb) or short DNAs (average 1.5 kb with size ranges of 0.5-3 kb) from a complex mixture. The method allows one to generate cDNA libraries enriched for long transcripts without size selection of insert DNAs. One representative library from newborn kidney includes 70% of clones bearing ATG start codons. A comparable library has been generated from 20 mouse blastocysts, containing only ∼40 ng of total RNA. This universal PCR amplification scheme can provide a route to isolate very large cDNAs, even if they are expressed at very low levels.

Due to brain tissue heterogeneity, the molecular genetic profile of any neurotransmitter-specific neuronal subtype is unknown. The purpose of this study was to purify a population of dopamine neurons, construct a cDNA library, and generate an initial gene expression profile and a microarray representative of dopamine neuron transcripts. Ventral mesencephalic dopamine neurons were purified by fluorescent-activated cell sorting from embryonic day 13.5 transgenic mice harboring a 4.5-kb rat tyrosine hydroxylase promoter-lacZ fusion. Nine-hundred sixty dopamine neuron cDNA clones were sequenced and arrayed for use in studies of gene expression changes during methamphetamine neurotoxicity. A neurotoxic dose of methamphetamine produced a greater than twofold up-regulation of the mitochondrial cytochrome c oxidase polypeptide I transcript from adult mouse substantia nigra at 12 h posttreatment. This is the first work to describe a gene expression profile for a neuronal subtype and to identify gene expression changes during methamphetamine neurotoxicity. © 2001 Academic Press.

Fetuin/alpha2-HS glycoprotein (alpha2-HSG) homologs have been identified in several species including rat, sheep, pig, rabbit, guinea pig, cattle, mouse and human. Multiple physiological roles for these homologs have been suggested, including ability to bind to hydroxyapatite crystals and to specifically inhibit the tyrosine kinase (TK) activity of the insulin receptor (IR). In this study we report the identification, cloning, and characterization of the mouse Ahsg gene and its function as an IR-TK inhibitor. Genomic clones derived from a mouse Svj 129 genomic library were sequenced in order to characterize the intron-exon organization of the mouse Ahsg gene, including an 875 bp subclone containing 154 bp upstream from the transcription start site, the first exon, and part of the first intron. A second genomic subclone harboring a 3.45 kb Bgl II fragment contained exons 2, 3 and 4 in addition to two adjacent elements within the first intron-a repetitive element of the B1 family (92 bp) and a 271 bp tract of (T,C)n*(A,G)n. We have mapped mouse Ahsg at 16 cM adjacent to the Diacylglycerol kinase 3 (Dagk3) gene on chromosome 16 by genotyping interspecific backcross panels between C57BL/6J and Mus spretus. The position is syntenic with human chromosome 3q27, where the human AHSG gene resides. Using recombinant mouse alpha2-HSG expressed from a recombinant baculovirus, we demonstrate that mouse alpha2-HSG inhibits insulin-stimulated IR autophosphorylation and IR-TKA in vitro. In addition, mouse alpha2-HSG (25 microg/ml) completely abolishes insulin-induced DNA synthesis in H-35 rat hepatoma cells. Based on the sequence data and functional analysis, we conclude that the mouse Ahsg gene is the true ortholog of the human AHSG gene.

Fundamental questions in developmental biology are: what genes are expressed, where and when they are expressed, what is the level of expression and how are these programs changed by the functional and structural alteration of genes? These questions have been addressed by studying one gene at a time, but a new research field that handles many genes in parallel is emerging. The methodology is at the interface of large-scale genomics approaches and developmental biology. Genomics needs developmental biology because one of the goals of genomics - collection and analysis of all genes in an organism - cannot be completed without working on embryonic tissues in which many genes are uniquely expressed. However, developmental biology needs genomics - the high-throughput approaches of genomics generate information about genes and pathways that can give an integrated view of complex processes. This article discusses these new approaches and their applications to mammalian developmental biology.

Targeted sequencing of the mouse t-complex has started with a 176-kb, gene-rich BAC localized with six PCR-based markers in inversion 2/3 of the highly duplicated region. The sequence contains 11 genes recovered primarily as cDNAs from early embryonic collections, including Igfals (previously placed on chromosome 17), Nubp2 (a fully characterized gene), Jsap1 (a JNK-binding protein), Rsp29 (the mouse homologue of the rat gene), Ndk3 (a nucleoside diphosphate kinase), and six additional putative genes of unknown function. With 50% GC content, 75% of the DNA transcribed, and one gene/16.0 kb (on average), the region may qualify as one of the most gene-dense segments in the mouse genome and provides candidates for dosage-sensitive phenotypes and mouse embryonic lethals mapped to the vicinity.

A novel human X-linked gene shows placenta-specific expression and has been named PLAC1. The gene maps 65 kb telomeric to HPRT at Xq26 and has been completely sequenced at the cDNA and genomic levels. The mouse orthologue Plac1 maps to the syntenically equivalent region of the mouse X chromosome. In situ hybridization studies with the antisense mRNA during mouse embryogenesis detect Plac1 expression from 7.5 dpc (days postcoitum) to 14.5 dpc in ectoplacental cone, giant cells, and labyrinthine trophoblasts. The putative human and murine PLAC1 proteins are 60% identical and 77% homologous. Both include a signal peptide and a peptide sequence also found in an interaction domain of the ZP3 (zona pellucida 3) protein. These results make PLAC1 a marker for placental development, with a possible role in the establishment of the mother-fetus interface. (C) 2000 Academic Press.

Phemx (Pan hematopoietic expression) is a novel murine gene expressed in developmentally regulated sites of hematopoiesis from early in embryogenesis through adulthood. Phemx is expressed in hematopoietic progenitors and mature cells of the three main hematopoietic lineages. Conceptual translation of the murine Phemx cDNA predicts a 25-kDa polypeptide with four hydrophobic regions and several potential phosphorylation sites, suggestive of a transmembrane protein involved in cell signaling. The PHEMX protein is structurally similar to tetraspanin CD81 (TAPA-1), a transmembrane protein involved in leukocyte activation, adhesion, and proliferation. Phemx maps to the distal region of chromosome 7, a segment of the mouse genome that contains a cluster of genes that exhibit genomic imprinting. However, imprinting analysis of Phemx at the whole organ level shows that it is biallelically expressed, suggesting that mechanisms leading to monoallelic expression are not imposed at this locus. The human PHEMX ortholog is specifically expressed in hematopoietic organs and tissues and, in contrast to murine Phemx, undergoes alternative splicing. The unique mode and range of Phemx expression suggest that it plays a role in hematopoietic cell function. (C) 2000 Academic Press.

cDNA microarray technology has been increasingly used to monitor global gene expression patterns in various tissues and cell types. However, applications to mammalian development have been hampered by the lack of appropriate cDNA collections, particularly for early developmental stages. To overcome this problem, a PCR-based cDNA library construction method was used to derive 52,374 expressed sequence tags from pre- and peri-implantation embryos, embryonic day (E) 12.5 female gonad/mesonephros, and newborn ovary. From these cDNA collections, a microarray representing 15,264 unique genes (78% novel and 22% known) was assembled. In initial applications, the divergence of placental and embryonic gene expression profiles was assessed. At stage E12.5 of development, based on triplicate experiments, 720 genes (6.5%) displayed statistically significant differences in expression between placenta and embryo. Among 289 more highly expressed in placenta, 61 placenta-specific genes encoded, for example, a novel prolactin-like protein. The number of genes highly expressed (and frequently specific) for placenta has thereby been increased 5-fold over the total previously reported, illustrating the potential of the microarrays for tissue-specific gene discovery and analysis of mammalian developmental programs.

Little is known about gene action in the preimplantation events that initiate mammalian development. Based on cDNA collections made from each stage from egg to blastocyst, 25438 3'-ESTs were derived, and represent 9718 genes, half of them novel. Thus, a considerable fraction of mammalian genes is dedicated to embryonic expression. This study reveals profound changes in gene expression that include the transient induction of transcripts at each stage. These results raise the possibility that development is driven by the action of a series of stage-specific expressed genes. The new genes, 798 of them placed on the mouse genetic map, provide entry points for analyses of human and mouse developmental disorders.

The 'NK-triggering-antigen regulator' (Nktar) gene is a locus identified in the C57BL/6 genome which regulates the ability of unlabeled activated Con A blasts to compete for recognition of labeled syngeneic Con A blasts by BALB/c NK cells. Linkage analysis on Con A blasts from (BALB/c x CByB6F1) N2 backcross progeny for (1) relative level of competitive inhibition of BALB/c NK lysis of syngeneic Con A blasts and (2) genotypes at polymorphic microsatellite markers distributed throughout the mouse genome mapped the Nktar gene locus to a 5-cM region of chromosome 10 containing the interferon- γ receptor (Ifngr) gene locus. N2 Con A blasts exhibited an inverse relationship between (a) their cell surface density of IFN-γR molecules detected by FACS with monoclonal anti-CD119 and (b) their cold target inhibition of BALB/c NK self-reactivity. Con A blasts from Ifngr(-1-) knockout mice showed a relatively high level of inhibition of BALB/c NK self- lysis and a relatively low level of class I MHC, which were both reversed by transient transfection with the Ifngr gene. Sequencing studies showed that Balb/c Ifngr encodes a Gly69 whereas C57BL/6 Ifngr encodes Glu69 due to a difference at nucleotide 284. Sequencing studies of N2 progeny demonstrated 100% concordance between their Nktar inhibitory phenotype and their Ifngr genotype. These findings demonstrate that the Nktar and Ifngr loci are identical. They further indicate that polymorphisms related to the Ifngr locus and affecting expression of cell surface IFN-γR may underlie genetic differences in the availability of NK-triggering antigens (NKTAgs) to recognition by self-reactive BALB/c NK cells by differentially affecting the ability of IFN-γ/R molecules to mediate up-regulation of NKTAg-masking class I molecules. (C) 2000 Academic Press.

The ability to quantitatively measure the expression of all genes in a given tissue or cell with a single assay is an exciting promise of gene-expression profiling technology. An in situ-synthesized 60-mer oligonucleotide microarray designed to detect transcripts from all mouse genes was validated, as well as a set of exogenous RNA controls derived from the yeast genome (made freely available without restriction), which allow quantitative estimation of absolute endogenous transcript abundance.

Interferon regulatory factor-8 (IRF-8)/interferon consensus sequence-binding protein (ICSBP) is a transcription factor that controls myeloid-cell development. Microarray gene expression analysis of Irf-8 -/- myeloid progenitor cells expressing an IRF-8/estrogen receptor chimera (which differentiate into macrophages after addition of estradiol) was used to identify 69 genes altered by IRF-8 during early differentiation (62 up-regulated and 7 down-regulated). Among them, 4 lysosomal/endosomal enzyme-related genes (cystatin C, cathepsin C, lysozyme, and prosaposin) did not require de novo protein synthesis for induction, suggesting that they were direct targets of IRF-8. We developed a reporter assay system employing a self-inactivating retrovirus and analyzed the cystatin C and cathepsin C promoters. We found that a unique cis element mediates IRF-8-induced activation of both promoters. Similar elements were also found in other IRF-8 target genes with a consensus sequence (GAAANN[N]GGAA) comprising a core IRF-binding motif and an Ets-binding motif; this sequence is similar but distinct from the previously reported Ets/IRF composite element. Chromatin immunoprecipitation assays demonstrated that IRF-8 and the PU.1 Ets transcription factor bind to this element in vivo. Collectively, these data indicate that IRF-8 stimulates transcription of target genes through a novel cis element to specify macrophage differentiation. © 2005 by The American Society of Hematology.

To build a mouse gene index with the most comprehensive coverage of alternative transcription/splicing (ATS), we developed an algorithm and a fully automated computational pipeline for transcript assembly from expressed sequences aligned to the genome. We identified 191,946 genomic loci, which included 27,497 protein-coding genes and 11,906 additional gene candidates (e.g., nonprotein-coding, but multiexon). Comparison of the resulting gene index with TIGR, UniGene, DoTS, and ESTGenes databases revealed that it had a greater number of transcripts, a greater average number of exons and introns with proper splicing sites per gene, and longer ORFs. The 27,497 protein-coding genes had 77,138 transcripts, i.e., 2.8 transcripts per gene on average. Close examination of transcripts led to a combinatorial table of 23 types of ATS units, only nine of which were previously described, i.e., 14 types of alternative splicing, seven types of alternative starts, and two types of alternative termination. The 47%, 18%, and 14% of 20,323 multiexon protein-coding genes with proper splice sites had alternative splicings, alternative starts, and alternative terminations, respectively. The gene index with the comprehensive ATS will provide a useful platform for analyzing the nature and mechanism of ATS, as well as for designing the accurate exon-based DNA microarrays. ©2005 by Cold Spring Harbor Laboratory Press.