Large-scale comparative sequence analysis of the human and murine Bruton's tyrosine kinase loci reveals conserved regulatory domains

High overexpression of the human alpha-galactosidase A gene driven by its promoter in transgenic mice: implications for the treatment of Fabry disease

BACKGROUND

Human alpha-galactosidase A (alpha-Gal A) is the lysosomal enzyme that cleaves alpha-galactosyl residues from glycoconjugates and is the deficient enzyme in Fabry disease. To date, there have been no studies on the regulation of this "housekeeping" gene.

METHODS

Transgenic mice were established with either 1) a 13.3-kilobase (kb) human genomic fragment that contained 246 bp of 5'- and approximately 2.8 kb of 3'- untranslated sequences, or 2) an "intronless" construct derived from the genomic sequence with the 5' and 3' flanking regions intact. Tissues that expressed high levels of alpha-Gal A activity were examined by light and electron microscopy.

RESULTS

Transgenic mice were generated with 2 and 12 copies of the genomic sequence (Lines 1 and 2) or about 60 copies of the intronless construct (Lines 3 and 4). In mice hemizygous for the genomic transgene (Lines 1 and 2), tissue alpha-Gal A activities were 12 to 155 times higher than those in the respective wild-type tissue, depending on tissue and transgene copy number. Of note, the high overexpression did not alter the cellular or subcellular cytoarchitecture. In contrast, alpha-Gal A activities expressed by mice that carried the intronless construct were only two- to sixfold more than in wild-type tissues in which the genomic transgene was highly expressed.

CONCLUSIONS

The remarkably high levels of alpha-Gal A expression in transgenic mice carrying the intact genomic sequence versus the intronless construct suggested that the genomic sequence contained a strong intronic enhancer element. Identification of this regulatory element or elements may be useful in efforts to overexpress human alpha-Gal A for gene therapy endeavors. In addition, overexpression of human alpha-Gal A did not affect cellular morphology, which indicates that its overexpression in gene therapy endeavors should be safe.

rVista for comparative sequence-based discovery of functional transcription factor binding sites

Identifying transcriptional regulatory elements represents a significant challenge in annotating the genomes of higher vertebrates. We have developed a computational tool, rVista, for high-throughput discovery of cis-regulatory elements that combines clustering of predicted transcription factor binding sites (TFBSs) and the analysis of interspecies sequence conservation to maximize the identification of functional sites. To assess the ability of rVista to discover true positive TFBSs while minimizing the prediction of false positives, we analyzed the distribution of several TFBSs across 1 Mb of the well-annotated cytokine gene cluster (Hs5q31; Mm11). Because a large number of AP-1, NFAT, and GATA-3 sites have been experimentally identified in this interval, we focused our analysis on the distribution of all binding sites specific for these transcription factors. The exploitation of the orthologous human-mouse dataset resulted in the elimination of > 95% of the approximately 58,000 binding sites predicted on analysis of the human sequence alone, whereas it identified 88% of the experimentally verified binding sites in this region.

Complete mutation analysis panel of the 39 human HOX genes

BACKGROUND

The HOX gene family consists of highly conserved transcription factors that specify the identity of the body segments along the anteroposterior axis of the embryo. Because the phenotypes of mice with targeted disruptions of Hox genes resemble some patterns of human malformations, mutations in HOX genes have been expected to be associated with a significant number of human malformations. Thus far, however, mutations have been documented in only three of the 39 human HOX genes (HOXD13, HOXA13, and HOXA11) partly because current knowledge on the complete coding sequence and genome structure is limited to only 20 of the 39 human HOX genes.

METHODS

Taking advantage of the human and mouse draft genome sequences, we attempted to characterize the remaining 19 human HOX genes by bioinformatic analysis including phylogenetic footprinting, the probabilistic prediction method, and comparison of genomic sequences with the complete set of the human anonymous cDNA sequences.

RESULTS

We were able to determine the full coding sequences of 19 HOX genes and their genome structure and successfully designed a complete set of PCR primers to amplify the entire coding region of each of the 39 HOX genes from genomic DNA.

CONCLUSIONS

Our results indicate the usefulness of bioinformatic analysis of the draft genome sequences for clinically oriented research projects. It is hoped that the mutation panel provided here will serve as a launchpad for a new discourse on the genetic basis of human malformations.

Generation and comparative analysis of approximately 3.3 Mb of mouse genomic sequence orthologous to the region of human chromosome 7q11.23 implicated in Williams syndrome

Williams syndrome is a complex developmental disorder that results from the heterozygous deletion of a approximately 1.6-Mb segment of human chromosome 7q11.23. These deletions are mediated by large (approximately 300 kb) duplicated blocks of DNA of near-identical sequence. Previously, we showed that the orthologous region of the mouse genome is devoid of such duplicated segments. Here, we extend our studies to include the generation of approximately 3.3 Mb of genomic sequence from the mouse Williams syndrome region, of which just over 1.4 Mb is finished to high accuracy. Comparative analyses of the mouse and human sequences within and immediately flanking the interval commonly deleted in Williams syndrome have facilitated the identification of nine previously unreported genes, provided detailed sequence-based information regarding 30 genes residing in the region, and revealed a number of potentially interesting conserved noncoding sequences. Finally, to facilitate comparative sequence analysis, we implemented several enhancements to the program, including the addition of links from annotated features within a generated percent-identity plot to specific records in public databases. Taken together, the results reported here provide an important comparative sequence resource that should catalyze additional studies of Williams syndrome, including those that aim to characterize genes within the commonly deleted interval and to develop mouse models of the disorder.

Homology between a 173-kb region from mouse chromosome 10, telomeric to the Ifng locus, and human chromosome 12q15

We sequenced a 173-kb region of mouse chromosome 10, telomeric to the Ifng locus, and compared it with the human homologous sequence located on chromosome 12q15 using various sequence analysis programs. This region has a low density of genes: one gene was detected in the mouse and the human sequences and a second gene was detected only in the human sequence. The mouse gene and its human orthologue, which are expressed in the immune system at a low level, produce a noncoding mRNA. Nonexpressed sequences show a higher degree of conservation than exons in this genomic region. At least three of these conserved sequences are also conserved in a third mammalian species (sheep or cow).

Evolutionarily conserved sequences on human chromosome 21

Comparison of human sequences with the DNA of other mammals is an excellent means of identifying functional elements in the human genome. Here we describe the utility of high-density oligonucleotide arrays as a rapid approach for comparing human sequences with the DNA of multiple species whose sequences are not presently available. High-density arrays representing approximately 22.5 Mb of nonrepetitive human chromosome 21 sequence were synthesized and then hybridized with mouse and dog DNA to identify sequences conserved between humans and mice (human-mouse elements) and between humans and dogs (human-dog elements). Our data show that sequence comparison of multiple species provides a powerful empiric method for identifying actively conserved elements in the human genome. A large fraction of these evolutionarily conserved elements are present in regions on chromosome 21 that do not encode known genes.

cDNA sequence and genomic structure of the rat RET proto-oncogene

The RET proto-oncogene, a member of the Receptor Tyrosine Kinase family, plays a crucial role during the development of the excretory system and the enteric nervous system, as demonstrated by in vivo animal studies and by its involvement in the pathogenesis of several human neurocristopathies like Hirschsprung disease and Multiple Endocrine Neoplasia type 2. Using a multistep RT-PCR approach we have isolated and sequenced the cDNA of the whole rat RET proto-oncogene, reporting the deduced amino acid sequence in comparison with the human and mouse counterparts. Moreover, two different isoforms (RET9 and RET51) have been confirmed in the rat, while a third RET isoform demonstrated in human (RET43) has not resulted to be conserved in this species. Finally, we have determined the genomic structure of the rat RET proto-oncogene comparing the exon-intron boundaries and intron sizes with the known structure of the human homologous gene. Our findings will facilitate the molecular study of appropriate rat models of RET related human diseases.

Analysis of the cat eye syndrome critical region in humans and the region of conserved synteny in mice: a search for candidate genes at or near the human chromosome 22 pericentromere

We have sequenced a 1.1-Mb region of human chromosome 22q containing the dosage-sensitive gene(s) responsible for cat eye syndrome (CES) as well as the 450-kb homologous region on mouse chromosome 6. Fourteen putative genes were identified within or adjacent to the human CES critical region (CESCR), including three known genes (IL-17R, ATP6E, and BID) and nine novel genes, based on EST identity. Two putative genes (CECR3 and CECR9) were identified, in the absence of EST hits, by comparing segments of human and mouse genomic sequence around two solitary amplified exons, thus showing the utility of comparative genomic sequence analysis in identifying transcripts. Of the 14 genes, 10 were confirmed to be present in the mouse genomic sequence in the same order and orientation as in human. Absent from the mouse region of conserved synteny are CECR1, a promising CES candidate gene from the center of the contig, neighboring CECR4, and CECR7 and CECR8, which are located in the gene-poor proximal 400 kb of the contig. This latter proximal region, located approximately 1 Mb from the centromere, shows abundant duplicated gene fragments typical of pericentromeric DNA. The margin of this region also delineates the boundary of conserved synteny between the CESCR and mouse chromosome 6. Because the proximal CESCR appears abundant in duplicated segments and, therefore, is likely to be gene poor, we consider the putative genes identified in the distal CESCR to represent the majority of candidate genes for involvement in CES.

Sequence conservation at human and mouse orthologous common fragile regions, FRA3B/FHIT and Fra14A2/Fhit

It has been suggested that delayed DNA replication underlies fragility at common human fragile sites, but specific sequences responsible for expression of these inducible fragile sites have not been identified. One approach to identify such cis-acting sequences within the large nonexonic regions of fragile sites would be to identify conserved functional elements within orthologous fragile sites by interspecies sequence comparison. This study describes a comparison of orthologous fragile regions, the human FRA3B/FHIT and the murine Fra14A2/Fhit locus. We sequenced over 600 kbp of the mouse Fra14A2, covering the region orthologous to the fragile epicenter of FRA3B, and determined the Fhit deletion break points in a mouse kidney cancer cell line (RENCA). The murine Fra14A2 locus, like the human FRA3B, was characterized by a high AT content. Alignment of the two sequences showed that this fragile region was stable in evolution despite its susceptibility to mitotic recombination on inhibition of DNA replication. There were also several unusual highly conserved regions (HCRs). The positions of predicted matrix attachment regions (MARs), possibly related to replication origins, were not conserved. Of known fragile region landmarks, five cancer cell break points, one viral integration site, and one aphidicolin break cluster were located within or near HCRs. Thus, comparison of orthologous fragile regions has identified highly conserved sequences with possible functional roles in maintenance of fragility.

Characterization of the genomic Xist locus in rodents reveals conservation of overall gene structure and tandem repeats but rapid evolution of unique sequence

The Xist locus plays a central role in the regulation of X chromosome inactivation in mammals, although its exact mode of action remains to be elucidated. Evolutionary studies are important in identifying conserved genomic regions and defining their possible function. Here we report cloning, sequence analysis, and detailed characterization of the Xist gene from four closely related species of common vole (field mouse), Microtus arvalis. Our analysis reveals that there is overall conservation of Xist gene structure both between different vole species and relative to mouse and human Xist/XIST. Within transcribed sequence, there is significant conservation over five short regions of unique sequence and also over Xist-specific tandem repeats. The majority of unique sequences, however, are evolving at an unexpectedly high rate. This is also evident from analysis of flanking sequences, which reveals a very high rate of rearrangement and invasion of dispersed repeats. We discuss these results in the context of Xist gene function and evolution.

Long-range comparison of human and mouse SCL loci: localized regions of sensitivity to restriction endonucleases correspond precisely with peaks of conserved noncoding sequences

Long-range comparative sequence analysis provides a powerful strategy for identifying conserved regulatory elements. The stem cell leukemia (SCL) gene encodes a bHLH transcription factor with a pivotal role in hemopoiesis and vasculogenesis, and it displays a highly conserved expression pattern. We present here a detailed sequence comparison of 193 kb of the human SCL locus to 234 kb of the mouse SCL locus. Four new genes have been identified together with an ancient mitochondrial insertion in the human locus. The SCL gene is flanked upstream by the SIL gene and downstream by the MAP17 gene in both species, but the gene order is not collinear downstream from MAP17. To facilitate rapid identification of candidate regulatory elements, we have developed a new sequence analysis tool (SynPlot) that automates the graphical display of large-scale sequence alignments. Unlike existing programs, SynPlot can display the locus features of more than one sequence, thereby indicating the position of homology peaks relative to the structure of all sequences in the alignment. In addition, high-resolution analysis of the chromatin structure of the mouse SCL gene permitted the accurate positioning of localized zones accessible to restriction endonucleases. Zones known to be associated with functional regulatory regions were found to correspond precisely with peaks of human/mouse homology, thus demonstrating that long-range human/mouse sequence comparisons allow accurate prediction of the extent of accessible DNA associated with active regulatory regions.

Human and mouse alpha-synuclein genes: comparative genomic sequence analysis and identification of a novel gene regulatory element

The human alpha-synuclein gene (SNCA) encodes a presynaptic nerve terminal protein that was originally identified as a precursor of the non-beta-amyloid component of Alzheimer's disease plaques. More recently, mutations in SNCA have been identified in some cases of familial Parkinson's disease, presenting numerous new areas of investigation for this important disease. Molecular studies would benefit from detailed information about the long-range sequence context of SNCA. To that end, we have established the complete genomic sequence of the chromosomal regions containing the human and mouse alpha-synuclein genes, with the objective of using the resulting sequence information to identify conserved regions of biological importance through comparative sequence analysis. These efforts have yielded approximately 146 and approximately 119 kb of high-accuracy human and mouse genomic sequence, respectively, revealing the precise genetic architecture of the alpha-synuclein gene in both species. A simple repeat element upstream of SNCA/Snca has been identified and shown to be necessary for normal expression in transient transfection assays using a luciferase reporter construct. Together, these studies provide valuable data that should facilitate more detailed analysis of this medically important gene.

Comparative sequence analysis of a region on human chromosome 13q14, frequently deleted in B-cell chronic lymphocytic leukemia, and its homologous region on mouse chromosome 14

Previous studies have indicated the presence of a putative tumor suppressor gene on human chromosome 13q14, commonly deleted in patients with B-cell chronic lymphocytic leukemia (B-CLL). We have recently identified a minimally deleted region encompassing parts of two adjacent genes, termed LEU1 and LEU2 (leukemia-associated genes 1 and 2), and several additional transcripts. In addition, 50 kb centromeric to this region we have identified another gene, LEU5/RFP2. To elucidate further the complex genomic organization of this region, we have identified, mapped, and sequenced the homologous region in the mouse. Fluorescence in situ hybridization analysis demonstrated that the region maps to mouse chromosome 14. The overall organization and gene order in this region were found to be highly conserved in the mouse. Sequence comparison between the human deletion hotspot region and its homologous mouse region revealed a high degree of sequence conservation with an overall score of 74%. However, our data also show that in terms of transcribed sequences, only two of those, human LEU2 and LEU5/RFP2, are clearly conserved, strengthening the case for these genes as putative candidate B-CLL tumor suppressor genes.

Active conservation of noncoding sequences revealed by three-way species comparisons

Human and mouse genomic sequence comparisons are being increasingly used to search for evolutionarily conserved gene regulatory elements. Large-scale human-mouse DNA comparison studies have discovered numerous conserved noncoding sequences of which only a fraction has been functionally investigated A question therefore remains as to whether most of these noncoding sequences are conserved because of functional constraints or are the result of a lack of divergence time.

Genomic sequence analysis of the mouse Naip gene array

A mouse locus called Lgn1 determines differences in macrophage permissiveness for the intracellular replication of Legionella pneumophila. The only regional candidate genes for this phenotype difference lie within a cluster of closely linked paralogs of the Neuronal Apoptosis Inhibitory Protein (Naip) gene. Previous genetic and physical mapping of the Lgn1 phenotype narrowed it to an interval containing only Naip2 and Naip5, suggesting that there is not complete functional overlap among the mouse Naip loci. In order to gather more information about polymorphisms among the Naip genes of the 129 mouse haplotype, we have determined the genomic sequence of a substantial portion of the 129 Naip gene array. We have constructed an evolutionary model for the expansion of the Naip gene array from a single progenitor Naip gene. This model predicts the presence of two distinct families of Naip paralogs: Naip1/2/3 and Naip4/5/6/7. Unlike the divergences among all the other Naip paralogs, the splits among Naip4, Naip5, Naip6, and Naip7 occurred relatively recently. The high degree of sequence conservation within the Naip4/5/6/7 family increases the likelihood of functional overlap among these genes.

Alfresco--a workbench for comparative genomic sequence analysis

Comparative analysis of genomic sequences provides a powerful tool for identifying regions of potential biologic function; by comparing corresponding regions of genomes from suitable species, protein coding or regulatory regions can be identified by their homology. This requires the use of several specific types of computational analysis tools. Many programs exist for these types of analysis; not many exist for overall view/control of the results, which is necessary for large-scale genomic sequence analysis. Using Java, we have developed a new visualization tool that allows effective comparative genome sequence analysis. The program handles a pair of sequences from putatively homologous regions in different species. Results from various different existing external analysis programs, such as database searching, gene prediction, repeat masking, and alignment programs, are visualized and used to find corresponding functional sequence domains in the two sequences. The user interacts with the program through a graphic display of the genome regions, in which an independently scrollable and zoomable symbolic representation of the sequences is shown. As an example, the analysis of two unannotated orthologous genomic sequences from human and mouse containing parts of the UTY locus is presented.

Cloning and expression of mouse deafness dystonia peptide 1 cDNA

Complementary DNA of mouse deafness dystonia peptide 1 (DDP1) was isolated from adipocyte cDNA library and expressed in mammalian cells. The sequence shares homology of 92 and 97% on the nucleic acid and the amino acid levels with human DDP1. In comparison to mouse Bruton's tyrosine kinase (Btk) locus, the coding region spans 2 exons and the splice point is the same as human DDP1. Northern blot analysis suggests that mouse DDP1 expresses ubiquitously. In vitro transcription/translation study showed that the cDNA of mouse DDP1 codes about 11 kDa peptide. DDP1 tagged with FLAG localized in mitochondria and cytoplasm of COS7 cells. P19 embryonal carcinoma cells transfected with anti sense DDP1 cDNA were frequently dead after subculture and all the survivals expressed endogenous DDP1 mRNA. Therefore, mouse DDP1 may play an important role to survive in contrast to Tim8p, a yeast homologue, which was unnecessary in yeast.

Human and mouse gene structure: comparative analysis and application to exon prediction

We describe a novel analytical approach to gene recognition based on cross-species comparison. We first undertook a comparison of orthologous genomic loci from human and mouse, studying the extent of similarity in the number, size and sequence of exons and introns. We then developed an approach for recognizing genes within such orthologous regions by first aligning the regions using an iterative global alignment system and then identifying genes based on conservation of exonic features at aligned positions in both species. The alignment and gene recognition are performed by new programs called and, respectively. performed well at exact identification of coding exons in 117 orthologous pairs tested.

Mutations of the human BTK gene coding for bruton tyrosine kinase in X-linked agammaglobulinemia

X-linked agammaglobulinemia (XLA) is an immunodeficiency caused by mutations in the gene coding for Bruton agammaglobulinemia tyrosine kinase (BTK). A database (BTKbase) of BTK mutations lists 544 mutation entries from 471 unrelated families showing 341 unique molecular events. In addition to mutations, a number of variants or polymorphisms have been found. Mutations in all the five domains of BTK cause the disease, the single most common event being missense mutations. Most mutations lead to truncation of the enzyme. The mutations appear almost uniformly throughout the molecule. About one-third of point mutations affect CpG sites, which usually code for arginine residues. The putative structural implications of all the missense mutations are provided in the database. BTKbase is available at http://www.uta.fi/imt/bioinfo.

Comparative genome sequence analysis of the Bpa/Str region in mouse and Man

The progress of human and mouse genome sequencing programs presages the possibility of systematic cross-species comparison of the two genomes as a powerful tool for gene and regulatory element identification. As the opportunities to perform comparative sequence analysis emerge, it is important to develop parameters for such analyses and to examine the outcomes of cross-species comparison. Our analysis used gene prediction and a database search of 430 kb of genomic sequence covering the Bpa/Str region of the mouse X chromosome, and 745 kb of genomic sequence from the homologous human X chromosome region. We identified 11 genes in mouse and 13 genes and two pseudogenes in human. In addition, we compared the mouse and human sequences using pairwise alignment and searches for evolutionary conserved regions (ECRs) exceeding a defined threshold of sequence identity. This approach aided the identification of at least four further putative conserved genes in the region. Comparative sequencing revealed that this region is a mosaic in evolutionary terms, with considerably more rearrangement between the two species than realized previously from comparative mapping studies. Surprisingly, this region showed an extremely high LINE and low SINE content, low G+C content, and yet a relatively high gene density, in contrast to the low gene density usually associated with such regions.

Comparative genomic sequencing identifies novel tissue-specific enhancers and sequence elements for methylation-sensitive factors implicated in Igf2/H19 imprinting

A differentially methylated region (DMR) and endoderm-specific enhancers, located upstream and downstream of the mouse H19 gene, respectively, are known to be essential for the reciprocal imprinting of Igf2 and H19. To explain the same imprinting patterns in non-endodermal tissues, additional enhancers have been hypothesized. We determined and compared the sequences of human and mouse H19 over 40 kb and identified 10 evolutionarily conserved downstream segments, 2 of which were coincident with the known enhancers. Reporter assays in transgenic mice showed that 5 of the other 8 segments functioned as enhancers in specific mesodermal and/or ectodermal tissues. We also identified a conserved 39-bp element that appeared repeatedly within the DMR and formed complexes with specific nuclear factors. Binding of one of the factors was inhibited when the target sequence contained methylated CpGs. These complexes may contribute to the presumed boundary function of the unmethylated DMR, which is proposed to insulate maternal Igf2 from the enhancers. Our results demonstrate that comparative genomic sequencing is highly efficient in identifying regulatory elements.

Shotgun sample sequence comparisons between mouse and human genomes

A mixed 'clone-by-clone' and 'whole-genome shotgun' strategy will be used to determine the genomic sequence of the mouse. This method will allow a phase of rapid annotation of the contemporaneous human sequence draft, through whole-genome 'sample sequence comparisons'.

PipMaker--a web server for aligning two genomic DNA sequences

PipMaker (http://bio.cse.psu.edu) is a World-Wide Web site for comparing two long DNA sequences to identify conserved segments and for producing informative, high-resolution displays of the resulting alignments. One display is a percent identity plot (pip), which shows both the position in one sequence and the degree of similarity for each aligning segment between the two sequences in a compact and easily understandable form. Positions along the horizontal axis can be labeled with features such as exons of genes and repetitive elements, and colors can be used to clarify and enhance the display. The web site also provides a plot of the locations of those segments in both species (similar to a dot plot). PipMaker is appropriate for comparing genomic sequences from any two related species, although the types of information that can be inferred (e.g., protein-coding regions and cis-regulatory elements) depend on the level of conservation and the time and divergence rate since the separation of the species. Gene regulatory elements are often detectable as similar, noncoding sequences in species that diverged as much as 100-300 million years ago, such as humans and mice, Caenorhabditis elegans and C. briggsae, or Escherichia coli and Salmonella spp. PipMaker supports analysis of unfinished or "working draft" sequences by permitting one of the two sequences to be in unoriented and unordered contigs.

Comparative genomic sequence analysis of the human and mouse cystic fibrosis transmembrane conductance regulator genes

The identification of the cystic fibrosis transmembrane conductance regulator gene (CFTR) in 1989 represents a landmark accomplishment in human genetics. Since that time, there have been numerous advances in elucidating the function of the encoded protein and the physiological basis of cystic fibrosis. However, numerous areas of cystic fibrosis biology require additional investigation, some of which would be facilitated by information about the long-range sequence context of the CFTR gene. For example, the latter might provide clues about the sequence elements responsible for the temporal and spatial regulation of CFTR expression. We thus sought to establish the sequence of the chromosomal segments encompassing the human CFTR and mouse Cftr genes, with the hope of identifying conserved regions of biologic interest by sequence comparison. Bacterial clone-based physical maps of the relevant human and mouse genomic regions were constructed, and minimally overlapping sets of clones were selected and sequenced, eventually yielding approximately 1.6 Mb and approximately 358 kb of contiguous human and mouse sequence, respectively. These efforts have produced the complete sequence of the approximately 189-kb and approximately 152-kb segments containing the human CFTR and mouse Cftr genes, respectively, as well as significant amounts of flanking DNA. Analyses of the resulting data provide insights about the organization of the CFTR/Cftr genes and potential sequence elements regulating their expression. Furthermore, the generated sequence reveals the precise architecture of genes residing near CFTR/Cftr, including one known gene (WNT2/Wnt2) and two previously unknown genes that immediately flank CFTR/Cftr.

Analysis of vertebrate SCL loci identifies conserved enhancers

The SCL gene encodes a highly conserved bHLH transcription factor with a pivotal role in hemopoiesis and vasculogenesis. We have sequenced and analyzed 320 kb of genomic DNA composing the SCL loci from human, mouse, and chicken. Long-range sequence comparisons demonstrated multiple peaks of human/mouse homology, a subset of which corresponded precisely with known SCL enhancers. Comparisons between mammalian and chicken sequences identified some, but not all, SCL enhancers. Moreover, one peak of human/mouse homology (+23 region), which did not correspond to a known enhancer, showed significant homology to an analogous region of the chicken SCL locus. A transgenic Xenopus reporter assay was established and demonstrated that the +23 region contained a new neural enhancer. This combination of long-range comparative sequence analysis with a high-throughput transgenic bioassay provides a powerful strategy for identifying and characterizing developmentally important enhancers.

Generation and analysis of 25 Mb of genomic DNA from the pufferfish Fugu rubripes by sequence scanning

We have generated and analyzed >50,000 shotgun clones from 1059 Fugu cosmid clones. All sequences have been minimally edited and searched against protein and DNA databases. These data are all displayed on a searchable, publicly available web site at. With an average of 50 reads per cosmid, this is virtually nonredundant sequence skimming, covering 30%-50% of each clone. This essentially random data set covers nearly 25 Mb (>6%) of the Fugu genome and forms the basis of a series of whole genome analyses which address questions regarding gene density and distribution in the Fugu genome and the similarity between Fugu and mammalian genes. The Fugu genome, with eight times less DNA but a similar gene repertoire, is ideally suited to this type of study because most cosmids contain more than one identifiable gene. General features of the genome are also discussed. We have made some estimation of the syntenic relationship between mammals and Fugu and looked at the efficacy of ORF prediction from short, unedited Fugu genomic sequences. Comparative DNA sequence analyses are an essential tool in the functional interpretation of complex vertebrate genomes. This project highlights the utility of using the Fugu genome in this kind of study.

Complete sequence determination of the mouse and human CTLA4 gene loci: cross-species DNA sequence similarity beyond exon borders

CTLA4 (CD152), a receptor for the B7 costimulatory molecules (CD80 and CD86), is considered a fundamental regulator of T-cell activation. In this paper, we present the complete primary structure of the mouse and human CTLA4 gene loci. Sequence comparison between the mouse and the human CTLA4 gene loci revealed a high degree of sequence conservation both for homologous noncoding regions (65-78% identity) and for coding regions (72-98% identity), with an overall score of 71% over the entire length of the two genes. Of the CTLA4 genomic regions aligned, five simple repetitive elements were found in the mouse locus, whereas two simple repetitive sequences were localized on the human locus. RNA blot analysis of mouse and human primary tissues indicated that both CTLA4 and T-cell receptor transcripts were found in most organs with generally higher levels in lymphoid tissues. The conservation of CTLA4 gene patterning raises the possibility that constrained gene evolution of CTLA4 may be linked to conserved transcriptional control of this locus.

Comparative analysis of noncoding regions of 77 orthologous mouse and human gene pairs

A data set of 77 genomic mouse/human gene pairs has been compiled from the EMBL nucleotide database, and their corresponding features determined. This set was used to analyze the degree of conservation of noncoding sequences between mouse and human. A new alignment algorithm was developed to cope with the fact that large parts of noncoding sequences are not alignable in a meaningful way because of genetic drift. This new algorithm, DNA Block Aligner (DBA), finds colinear-conserved blocks that are flanked by nonconserved sequences of varying lengths. The noncoding regions of the data set were aligned with DBA. The proportion of the noncoding regions covered by blocks >60% identical was 36% for upstream regions, 50% for 5' UTRs, 23% for introns, and 56% for 3' UTRs. These blocks of high identity were more or less evenly distributed across the length of the features, except for upstream regions in which the first 100 bp upstream of the transcription start site was covered in up to 70% of the gene pairs. This data set complements earlier sets on the basis of cDNA sequences and will be useful for further comparative studies. [This paper contains supplementary data that can be found at http://www.genome.org [corrected]].

Comparative sequence analysis of the mouse and human Lgn1/SMA interval

Human chromosome 5q11.2-q13.3 and its ortholog on mouse chromosome 13 contain candidate genes for an inherited human neurodegenerative disorder called spinal muscular atrophy (SMA) and for an inherited mouse susceptibility to infection with Legionella pneumophila (Lgn1). These homologous genomic regions also have unusual repetitive organizations that create practical difficulties in mapping and raise interesting issues about the evolutionary origin of the repeats. In an attempt to analyze this region in detail, and as a way to identify additional candidate genes for these diseases, we have determined the sequence of 179 kb of the mouse Lgn1/SMA interval. We have analyzed this sequence using BLAST searches and various exon prediction programs to identify potential genes. Since these methods can generate false-positive exon declarations, our alignments of the mouse sequence with available human orthologous sequence allowed us to discriminate rapidly among this collection of potential coding regions by indicating which regions were well conserved and were more likely to represent actual coding sequence. As a result of our analysis, we accurately mapped two additional genes in the SMA interval that can be tested for involvement in the pathogenesis of SMA. While no new Lgn1 candidates emerged, we have identified new genetic markers that exclude Smn as an Lgn1 candidate. In addition to providing important resources for studying SMA and Lgn1, our data provide further evidence of the value of sequencing the mouse genome as a means to help with the annotation of the human genomic sequence and vice versa.

Chromosomal location effects on gene sequence evolution in mammals

BACKGROUND

Nucleotide substitution rates and G + C content vary considerably among mammalian genes. It has been proposed that the mammalian genome comprises a mosaic of regions - termed isochores - with differing G + C content. The regional variation in gene G + C content might therefore be a reflection of the isochore structure of chromosomes, but the factors influencing the variation of nucleotide substitution rate are still open to question.

RESULTS

To examine whether nucleotide substitution rates and gene G + C content are influenced by the chromosomal location of genes, we compared human and murid (mouse or rat) orthologues known to belong to one of the chromosomal (autosomal) segments conserved between these species. Multiple members of gene families were excluded from the dataset. Sets of neighbouring genes were defined as those lying within 1 centiMorgan (cM) of each other on the mouse genetic map. For both synonymous substitution rates and G + C content at silent sites, neighbouring genes were found to be significantly more similar to each other than sets of genes randomly drawn from the dataset. Moreover, we demonstrated that the regional similarities in G + C content (isochores) and synonymous substitution rate were independent of each other.

CONCLUSIONS

Our results provide the first substantial statistical evidence for the existence of a regional variation in the synonymous substitution rate within the mammalian genome, indicating that different chromosomal regions evolve at different rates. This regional phenomenon which shapes gene evolution could reflect the existence of 'evolutionary rate units' along the chromosome.

Alignment of whole genomes

A new system for aligning whole genome sequences is described. Using an efficient data structure called a suffix tree, the system is able to rapidly align sequences containing millions of nucleotides. Its use is demonstrated on two strains of Mycoplasma tuberculosis, on two less similar species of Mycoplasma bacteria and on two syntenic sequences from human chromosome 12 and mouse chromosome 6. In each case it found an alignment of the input sequences, using between 30 s and 2 min of computation time. From the system output, information on single nucleotide changes, translocations and homologous genes can easily be extracted. Use of the algorithm should facilitate analysis of syntenic chromosomal regions, strain-to-strain comparisons, evolutionary comparisons and genomic duplications.

Isolation and characterization of the mouse Aire gene

Autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED) is a rare autosomal recessive disorder characterized by Addison's disease and/or hypoparathyroidism and/or chronic mucocutaneous candidiasis. Patients may also have other clinical symptoms both within and outside the endocrine system, mainly as a result of autoimmunity against organ-specific autoantigens. The gene for APECED has recently been identified and termed AIRE (for AutoImmune REgulator). APECED is a model of organ-specific autoimmunity and isolation and characterization of the homologous mouse gene, Aire, will provide tools for dissection of the mechanisms underlying this human disorder and defining molecular pathways involved in organ-specific autoimmunity. We have isolated and completely sequenced the mouse Aire gene which is split into 14 exons over 13 kb and encodes a predicted protein of 552 amino acids. The predicted mouse and human AIRE proteins are 71% identical and contain motifs suggestive of a transcriptional regulator. Additional conserved motifs are emerging in the AIRE/Aire proteins including a nuclear localization signal, an "ASS" domain, and a "SAND" domain. The human and mouse AIRE promoters have conserved sites for several thymus-specific transcription factors and others important in hematopoesis, consistent with its expression in rare cells of the thymus medulla, lymph nodes, and fetal liver. We have mapped mouse Aire to mouse chromosome 10 by FISH, to the same region as Pwp2 and Pfkl, confirming synteny to the corresponding region of human chromosome 21.

Comparative genomics: the key to understanding the Human Genome Project

The sequencing of the human genome is well underway. Technology has advanced, such that the total genomic sequence is possible, along with an extensive catalogue of genes via comprehensive cDNA libraries. With the recent completion of the Saccharomyces cerevisiae sequencing project and the imminent completion of that of Caenorhabditis elegans, the most frequently asked question is how much can sequence data alone tell us? The answer is that that a DNA sequence taken in isolation from a single organism reveals very little. The vast majority of DNA in most organisms is noncoding. Protein coding sequences or genes cannot function as isolated units without interaction with noncoding DNA and neighboring genes. This genomic environment is specific to each organism. In order to understand this we need to look at similar genes in different organisms, to determine how function and position has changed over the course of evolution. By understanding evolutionary processes we can gain a greater insight into what makes a gene and the wider processes of genetics and inheritance. Comparative genomics (with model organisms), once the poor relation of the human genome project, is starting to provide the key to unlock the DNA code.

Unusual mutations in Btk: an insertion, a duplication, an inversion, and four large deletions

Mutations in Bruton's tyrosine kinase (Btk) result in the immunodeficiency X-linked agammaglobulinemia (XLA). In a previous study of 101 patients with presumed XLA, we identified seven patients with large genomic alterations in Btk. The recent completion of 100 kb of contiguous DNA sequence at the Btk locus has allowed us to characterize these mutations in detail and to identify four different types of alterations. These alterations included a 253-bp retroposon insertion at position +5 within intron 9, an inversion of greater than 48 kb that disrupted Btk between exons 4 and 5, a 12.9-kb duplication including Btk exons 2 to 5, and four deletions ranging from 2.8 to 38 kb in size. The duplication and three of the deletions resulted from unequal crossovers of Alu repeats. Further, three of the deletions terminated within a repeat-rich cluster spanning 30 kb of sequence 3' of Btk exon 19, suggesting that this region was more susceptible to unequal crossovers than the rest of the Btk gene. These studies describe the first reports of an insertion, an inversion, and a duplication in Btk and demonstrate the utility of large-scale sequencing in the elucidation of disease-causing mutations.

Comparative Sequence of Human and Mouse BAC Clones from the mnd2 Region of Chromosome 2p13

The mnd2 mutation on mouse chromosome 6 produces a progressive neuromuscular disorder. To determine the gene content of the 400-kb mnd2 nonrecombinant region, we sequenced 108 kb of mouse genomic DNA and 92 kb of human genomic sequence from the corresponding region of chromosome 2p13.3. Three genes with the indicated sizes and intergenic distances were identified:D6Mm5e (⩾81 kb)–787 bp–DOK (2 kb)–845 bp–LOR2 (⩾6 kb). D6Mm5e is expressed in many tissues at very low abundance and the predicted 526-residue protein contains no known functional domains. DOK encodes the p62dok rasGAP binding protein involved in signal transduction. LOR2 encodes a novel lysyl oxidase-related protein of 757 amino acid residues. We describe a simple search protocol for identification of conserved internal exons in genomic sequence. Evolutionary conservation proved to be a useful criterion for distinguishing between authentic exons and artifactual products obtained by exon amplification, RT–PCR, and 5′ RACE. Conserved noncoding sequence elements longer than 80 bp with ⩾75% nucleotide sequence identity comprise ∼1% of the genomic sequence in this region. Comparative analysis of this human and mouse genomic DNA sequence was an efficient method for gene identification and is independent of developmental stage or quantitative level of gene expression.

[The sequence data described in this paper have been submitted to the GenBank data library under the following accession numbers: AC003061, mouse BAC clone 245c12; AC003065, human BAC clone h173(E10); AF053368, mouse Lor2 cDNA; AF084363, 108-kb contig from mouse BAC 245c12; AF084364, mouse D6Mm5ecDNA.]

Severe B cell deficiency and disrupted splenic architecture in transgenic mice expressing the E41K mutated form of Bruton's tyrosine kinase

To identify B-cell signaling pathways activated by Bruton's tyrosine kinase (Btk) in vivo, we generated transgenic mice in which Btk expression is driven by the MHC class II Ea gene locus control region. Btk overexpression did not have significant adverse effects on B cell function, and essentially corrected the X-linked immunodeficiency (xid) phenotype in Btk- mice. In contrast, expression of a constitutively activated form of Btk carrying the E41K gain-of-function mutation resulted in a B cell defect that was more severe than xid. The mice showed a marked reduction of the B cell compartment in spleen, lymph nodes, peripheral blood and peritoneal cavity. The levels in the serum of most immunoglobulin subclasses decreased with age, and B cell responses to both T cell-independent type II and T cell-dependent antigens were essentially absent. Expression of the E41K Btk mutant enhanced blast formation of splenic B cells in vitro in response to anti-IgM stimulation. Furthermore, the mice manifested a disorganization of B cell areas and marginal zones in the spleen. Our findings demonstrate that expression of constitutively activated Btk blocks the development of follicular recirculating B cells.

A computer program for aligning a cDNA sequence with a genomic DNA sequence

We address the problem of efficiently aligning a transcribed and spliced DNA sequence with a genomic sequence containing that gene, allowing for introns in the genomic sequence and a relatively small number of sequencing errors. A freely available computer program, described herein, solves the problem for a 100-kb genomic sequence in a few seconds on a workstation.

Mouse models of genetic disease: new approaches, new paradigms

The mouse mutant resource is a valuable tool for gene function studies in the post-genomics era. However, despite a seemingly large catalogue of mouse mutants, it is recognized that we have access to mutations at only a small fraction of the total number of mouse genes. There is a phenotype gap that needs to be narrowed by the implementation of large-scale, systematic mutagenesis programmes in the mouse. Both genotype-driven and phenotype-driven approaches can be employed to recover new mouse mutations. Genotype-driven approaches include large-scale genome-wide mutagenesis by gene trapping in embryonic stem cells. For genotype-driven approaches, the initial focus is on the characterization of the mutational change to the genome. Identification of the mutated gene is relatively trivial, but the genotype-driven route provides little indication of the likely phenotypic outcome of the mutation. In contrast, phenotype-driven approaches employ mutagenesis procedures that emphasize the recovery of novel phenotypes without prior assumptions about the underlying gene or pathway that has been disrupted--although identifying the underlying gene may not be trivial. One phenotype-driven approach includes chemical mutagenesis using N-ethyl-N-nitrosourea (ENU). ENU mutagenesis programmes are increasingly being brought to bear on increasing the breadth and depth of the mouse mutant resource, and in so doing narrowing the phenotype gap.

Genetic basis of abnormal B cell development

A susceptibility gene in the MHC class III region may underlie the defective B-cell differentiation in familial IgA deficiency and common variable immunodeficiency. Mutations in Bruton's tyrosine kinase, immunoglobulin heavy chain and lambda 5/14.1 surrogate light chain loci disrupt B-cell development to cause profound antibody deficiency. Mutational, biochemical and transgenic studies offer insight into the function of these and other 'antibody deficiency genes'.

Disruption and sequence identification of 2,000 genes in mouse embryonic stem cells

The dramatic increase in sequence information in the form of expressed sequence tags (ESTs) and genomic sequence has created a 'gene function gap' with the identification of new genes far outpacing the rate at which their function can be identified. The ability to create mutations in embryonic stem (ES) cells on a large scale by tagged random mutagenesis provides a powerful approach for determining gene function in a mammalian system; this approach is well established in lower organisms. Here we describe a high-throughput mutagenesis method based on gene trapping that allows the automated identification of sequence tags from the mutated genes. This method traps and mutates genes regardless of their expression status in ES cells. To facilitate the study of gene function on a large scale, we are using these techniques to create a library of ES cells called Omnibank, from which sequence-tagged mutations in 2,000 genes are described.

Physical mapping of the mouse genome

This review is intended to provide an overview of techniques and a source of reagents for physical mapping of the mouse genome. It focuses on those applications, methods, or resources unique to the mouse and on the generation of comparative physical maps. The reference list is not comprehensive; rather, recent reviews on each topic and selected representative examples are given.

Transcriptional Regulatory Elements Within the First Intron of Bruton's Tyrosine Kinase

Defects in the gene for Bruton's tyrosine kinase (Btk) result in the disorder X-linked agammaglobulinemia (XLA). Whereas XLA is characterized by a profound defect in B-cell development, Btk is expressed in both the B lymphocyte and myeloid cell lineages. We evaluated a patient with XLA who had reduced amounts of Btk transcript but no abnormalities in his coding sequence. A single base-pair substitution in the first intron of Btk was identified in this patient, suggesting that this region may contain regulatory elements. Using reporter constructs we identified two transcriptional control elements in the first 500 bp of intron 1. A strong positive regulator, active in both pre-B cells and B cells, was identified within the first 43 bp of the intron. Gel-shift assays identified two Sp1 binding sites within this element. The patient's mutation results in an altered binding specificity of the proximal Sp1 binding site. A negative regulator, active in pre-B cells only, was located between base pairs 281 and 491 of the intron. These findings indicate that regulation of Btk transcription is complex and may involve several transcriptional regulatory factors at the different stages of B-cell differentiation.

Transcriptional Regulatory Elements Within the First Intron of Bruton's Tyrosine Kinase

Defects in the gene for Bruton's tyrosine kinase (Btk) result in the disorder X-linked agammaglobulinemia (XLA). Whereas XLA is characterized by a profound defect in B-cell development, Btk is expressed in both the B lymphocyte and myeloid cell lineages. We evaluated a patient with XLA who had reduced amounts of Btk transcript but no abnormalities in his coding sequence. A single base-pair substitution in the first intron of Btk was identified in this patient, suggesting that this region may contain regulatory elements. Using reporter constructs we identified two transcriptional control elements in the first 500 bp of intron 1. A strong positive regulator, active in both pre-B cells and B cells, was identified within the first 43 bp of the intron. Gel-shift assays identified two Sp1 binding sites within this element. The patient's mutation results in an altered binding specificity of the proximal Sp1 binding site. A negative regulator, active in pre-B cells only, was located between base pairs 281 and 491 of the intron. These findings indicate that regulation of Btk transcription is complex and may involve several transcriptional regulatory factors at the different stages of B-cell differentiation.

BTKbase, mutation database for X-linked agammaglobulinemia (XLA)

X-linked agammaglobulinemia (XLA) is an immunodeficiency caused by mutations in the gene coding for Bruton's agammaglobulinemia tyrosine kinase (BTK). A database (BTKbase) of BTK mutations has been compiled and the recent update lists 463 mutation entries from 406 unrelated families showing 303 unique molecular events. In addition to mutations, the database also lists variants or polymorphisms. Each patient is given a unique patient identity number (PIN). Information is included regarding the phenotype including symptoms. Mutations in all the five domains of BTK have been noticed to cause the disease, the most common event being missense mutations. The mutations appear almost uniformly throughout the molecule and frequently affect CpG sites that code for arginine residues. The putative structural implications of all the missense mutations are given in the database. The improved version of the registry having a number of new features is available at http://www. helsinki.fi/science/signal/btkbase.html

Zooming in on the human-mouse comparative map: genome conservation re-examined on a high-resolution scale

Over the past decade, conservation of genetic linkage groups has been shown in mammals and used to great advantage, fueling significant exchanges of gene mapping and functional information especially between the genomes of humans and mice. As human physical maps increase in resolution from chromosome bands to nucleotide sequence, comparative alignments of mouse and human regions have revealed striking similarities and surprising differences between the genomes of these two best-mapped mammalian species. Whereas, at present, very few mouse and human regions have been compared on the physical level, existing studies provide intriguing insights to genome evolution, including the observation of recent duplications and deletions of genes that may play significant roles in defining some of the biological differences between the two species. Although high-resolution conserved marker-based maps are currently available only for human and mouse, a variety of new methods and resources are speeding the development of comparative maps of additional organisms. These advances mark the first step toward establishment of the human genome as a reference map for vertebrate species, providing evolutionary and functional annotation to human sequence and vast new resources for genetic analysis of a variety of commercially, medically, and ecologically important animal models.