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.

Identification of a coordinate regulator of interleukins 4, 13, and 5 by cross-species sequence comparisons

Long-range regulatory elements are difficult to discover experimentally; however, they tend to be conserved among mammals, suggesting that cross-species sequence comparisons should identify them. To search for regulatory sequences, we examined about 1 megabase of orthologous human and mouse sequences for conserved noncoding elements with greater than or equal to 70% identity over at least 100 base pairs. Ninety noncoding sequences meeting these criteria were discovered, and the analysis of 15 of these elements found that about 70% were conserved across mammals. Characterization of the largest element in yeast artificial chromosome transgenic mice revealed it to be a coordinate regulator of three genes, interleukin-4, interleukin-13, and interleukin-5, spread over 120 kilobases.

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.

Hox cluster genomics in the horn shark, Heterodontus francisci

Reconstructing the evolutionary history of Hox cluster origins will lead to insights into the developmental and evolutionary significance of Hox gene clusters in vertebrate phylogeny and to their role in the origins of various vertebrate body plans. We have isolated two Hox clusters from the horn shark, Heterodontus francisci. These have been sequenced and compared with one another and with other chordate Hox clusters. The results show that one of the horn shark clusters (HoxM) is orthologous to the mammalian HoxA cluster and shows a structural similarity to the amphioxus cluster, whereas the other shark cluster (HoxN) is orthologous to the mammalian HoxD cluster based on cluster organization and a comparison with noncoding and Hox gene-coding sequences. The persistence of an identifiable HoxA cluster over an 800-million-year divergence time demonstrates that the Hox gene clusters are highly integrated and structured genetic entities. The data presented herein identify many noncoding sequence motifs conserved over 800 million years that may function as genetic control motifs essential to the developmental process.

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.

Winnowing sequences from a database search

In database searches for sequence similarity, matches to a distinct sequence region (e.g., protein domain) are frequently obscured by numerous matches to another region of the same sequence. In order to cope with this problem, algorithms are developed to discard redundant matches. One model for this problem begins with a list of intervals, each with an associated score; each interval gives the range of positions in the query sequence that align to a database sequence, and the score is that of the alignment. If interval I is contained in interval J, and I's score is less than J's, then I is said to be dominated by J. The problem is then to identify each interval that is dominated by at least K other intervals, where K is a given level of "tolerable redundancy." An algorithm is developed to solve the problem in O(N log N) time and O(N*) space, where N is the number of intervals and N* is a precisely defined value that never exceeds N and is frequently much smaller. This criterion for discarding database hits has been implemented in the Blast program, as illustrated herein with examples. Several variations and extensions of this approach are also described.

A greedy algorithm for aligning DNA sequences

For aligning DNA sequences that differ only by sequencing errors, or by equivalent errors from other sources, a greedy algorithm can be much faster than traditional dynamic programming approaches and yet produce an alignment that is guaranteed to be theoretically optimal. We introduce a new greedy alignment algorithm with particularly good performance and show that it computes the same alignment as does a certain dynamic programming algorithm, while executing over 10 times faster on appropriate data. An implementation of this algorithm is currently used in a program that assembles the UniGene database at the National Center for Biotechnology Information.

Post-processing long pairwise alignments

MOTIVATION

The local alignment problem for two sequences requires determining similar regions, one from each sequence, and aligning those regions. For alignments computed by dynamic programming, current approaches for selecting similar regions may have potential flaws. For instance, the criterion of Smith and Waterman can lead to inclusion of an arbitrarily poor internal segment. Other approaches can generate an alignment scoring less than some of its internal segments.

RESULTS

We develop an algorithm that decomposes a long alignment into sub-alignments that avoid these potential imperfections. Our algorithm runs in time proportional to the original alignment's length. Practical applications to alignments of genomic DNA sequences are described.

Characterization of the human and mouse unconventional myosin XV genes responsible for hereditary deafness DFNB3 and shaker 2

Mutations in myosin XV are responsible for congenital profound deafness DFNB3 in humans and deafness and vestibular defects in shaker 2 mice. By combining direct cDNA analyses with a comparison of 95.2 kb of genomic DNA sequence from human chromosome 17p11.2 and 88.4 kb from the homologous region on mouse chromosome 11, we have determined the genomic and mRNA structures of the human (MYO15) and mouse (Myo15) myosin XV genes. Our results indicate that full-length myosin XV transcripts contain 66 exons, are >12 kb in length, and encode 365-kDa proteins that are unique among myosins in possessing very long approximately 1200-aa N-terminal extensions preceding their conserved motor domains. The tail regions of the myosin XV proteins contain two MyTH4 domains, two regions with similarity to the membrane attachment FERM domain, and a putative SH3 domain. Northern and dot blot analyses revealed that myosin XV is expressed in the pituitary gland in both humans and mice. Myosin XV transcripts were also observed by in situ hybridization within areas corresponding to the sensory epithelia of the cochlea and vestibular systems in the developing mouse inner ear. Immunostaining of adult mouse organ of Corti revealed that myosin XV protein is concentrated within the cuticular plate and stereocilia of cochlear sensory hair cells. These results indicate a likely role for myosin XV in the formation or maintenance of the unique actin-rich structures of inner ear sensory hair cells.

Comparison of five methods for finding conserved sequences in multiple alignments of gene regulatory regions

Conserved segments in DNA or protein sequences are strong candidates for functional elements and thus appropriate methods for computing them need to be developed and compared. We describe five methods and computer programs for finding highly conserved blocks within previously computed multiple alignments, primarily for DNA sequences. Two of the methods are already in common use; these are based on good column agreement and high information content. Three additional methods find blocks with minimal evolutionary change, blocks that differ in at most k positions per row from a known center sequence and blocks that differ in at most k positions per row from a center sequence that is unknown a priori. The center sequence in the latter two methods is a way to model potential binding sites for known or unknown proteins in DNA sequences. The efficacy of each method was evaluated by analysis of three extensively analyzed regulatory regions in mammalian beta-globin gene clusters and the control region of bacterial arabinose operons. Although all five methods have quite different theoretical underpinnings, they produce rather similar results on these data sets when their parameters are adjusted to best approximate the experimental data. The optimal parameters for the method based on information content varied little for different regulatory regions of the beta-globin gene cluster and hence may be extrapolated to many other regulatory regions. The programs based on maximum allowed mismatches per row have simple parameters whose values can be chosen a priori and thus they may be more useful than the other methods when calibration against known functional sites is not available.

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.

The gene mutated in bare patches and striated mice encodes a novel 3beta-hydroxysteroid dehydrogenase

X-linked dominant disorders that are exclusively lethal prenatally in hemizygous males have been described in human and mouse. None of the genes responsible has been isolated in either species. The bare patches (Bpa) and striated (Str) mouse mutations were originally identified in female offspring of X-irradiated males. Subsequently, additional independent alleles were described. We have previously mapped these X-linked dominant, male-lethal mutations to an overlapping region of 600 kb that is homologous to human Xq28 (ref. 4) and identified several candidate genes in this interval. Here we report mutations in one of these genes, Nsdhl, encoding an NAD(P)H steroid dehydrogenase-like protein, in two independent Bpa and three independent Str alleles. Quantitative analysis of sterols from tissues of affected Bpa mice support a role for Nsdhl in cholesterol biosynthesis. Our results demonstrate that Bpa and Str are allelic mutations and identify the first mammalian locus associated with an X-linked dominant, male-lethal phenotype. They also expand the spectrum of phenotypes associated with abnormalities of cholesterol metabolism.

Analysis of 9-fluorenylmethoxycarbonyl (FMOC) loading of solid-phase synthesis resins by gas chromatography

Base-catalyzed cleavage of the 9-fluorenylmethoxycarbonyl (FMOC) group and subsequent analysis by UV spectrophotometry is commonly employed to measure the "loading" of functional groups on solid support. Recent work suggests that 1,8-diazabicyclo[5.4. 0]undec-7-ene (DBU) may be superior to piperidine (the most commonly used base for the cleavage) for quantitative analysis by UV. We have compared deprotection of FMOC-bearing compounds by both DBU and piperidine, and have observed by GC-MS the formation of a dibenzofulvene-piperidine adduct (piperidine deprotection), and the formation of unassociated dibenzofulvene (DBU deprotection). We have further been able to use GC analysis of dibenzofulvene produced in the DBU deprotection mixture in a quantitative analysis of resin loading, which gave results comparable to UV methods. Sample preparation for this method has been automated using the Nautilus 2400 organic synthesizer to reduce the amount of operator time and increase throughput of sample analysis. Copyright 1998 John Wiley & Sons, Inc.

The mouse cornichon gene family

As part of a large scale mouse Expressed Sequence Tag (EST) project to identify molecules involved in the initiation of mammalian development, a homolog of the Drosophila cornichon gene was detected as a mouse maternal transcript present in the two-cell embryo. Cornichon is a multigene family in the mouse: the new gene, Cnih, maps to mouse chromosome 10, another cornichon homolog, Cnil, maps to chromosome 14 and two additional cornichon-related loci, possibly pseudogenes, localize to chromosomes 3 and 10, respectively. Cnih encodes an open reading frame (ORF) of 144 amino acids that is 93% homologous (68% identical) to the Drosophila protein, whereas the ORF of Cnil contains two extra polypeptide regions not found in these other proteins. Transcripts of Cnih are highly abundant in the full grown oocyte and the ovulated unfertilized egg, while Cnil message is only detectable after activation of the embryonic genome at the eight-cell stage. In situ hybridization shows specific localization of Cnih transcripts to ovarian oocytes. The lack of cytoplasmic polyadenylation of the maternally inherited Cnih transcript suggests that Cnih mRNA is translated in the full grown oocyte before, but not after, ovulation. In Drosophila, cornichon is involved in the establishment of both anterior-posterior and dorso-ventral polarity via the epidermal growth factor (EGF)-receptor signaling pathway. Finding Cnih in the mammalian oocyte opens a new perspective on the investigation of EGF-signaling in the oocyte.

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 p62(dok) 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 approximately 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 D6Mm5e cDNA.]

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

The Paced Auditory Serial Addition Task (PASAT): norms for age, education, and ethnicity

The Paced Auditory Serial Addition Task (PASAT) is often used to measure attention, concentration, working memory, and speed of information processing. Using a modified 200-item version of the PASAT with presentation rates of 3.0, 2.4, 2.0, and 1.6 items per second, we analyzed demographic influences on test performance in a large sample (N = 566) of healthy North American adults. We found that age, education, and ethnicity were significant predictors, accounting for nearly 23% of the variance in test performance. We discuss these results in comparison to previous normative studies, and present a formula and tables for computing age, education, and ethnicity-adjusted T scores for performance on the PASAT 200.

A 1996 NASN research award winner. Health services provided by school nurses for students with chronic health problems

Increased numbers of students are entering schools with chronic health problems and in need of nursing and other health services. Some studies have been conducted in recent years describing nursing practice in schools. However, further study of school nursing practice was warranted owing to two major problems. The first was limited financial resources to provide needed nursing services. The second was inadequate documentation of both the health needs of students and current school nursing practice in meeting those needs. The purpose of this descriptive study was to examine the nature and extent of nursing services provided in the schools in the late 1990s and to identify potentially reimbursable nursing services being provided to students with chronic health problems. A survey was mailed to active members of the School Nurse Organization of Minnesota (SNOM). The findings indicated the considerable diversity in the practice setting as well as the large volume of nursing services being delivered in schools. Results will be essential in the acquisition of financial resources, including collaboration with third-party payors for reimbursement of nursing and other health services.