In vitro hematopoietic and endothelial potential of flk-1−/− embryonic stem cells and embryos

Mice deficient in the Flk-1 receptor tyrosine kinase are known to die in utero because of defective vascular and hematopoietic development. Here, we show that flk-1−/− embryonic stem cells are nevertheless able to differentiate into hematopoietic and endothelial cells in vitro, although they give rise to a greatly reduced number of blast colonies, a measure of hemangioblast potential. Furthermore, normal numbers of hematopoietic progenitors are found in 7.5-day postcoitum flk-1−/− embryos, even though 8.5-day postcoitum flk-1−/− embryos are known to be deficient in such cells. Our results suggest that hematopoietic/endothelial progenitors arise independently of Flk-1, but that their subsequent migration and expansion require a Flk-1-mediated signal.

MSF (MLL septin-like fusion), a fusion partner gene of MLL, in a therapy-related acute myeloid leukemia with a t(11;17)(q23;q25)

MLL (ALL1, Htrx, HRX), which is located on chromosome band 11q23, frequently is rearranged in patients with therapy-related acute myeloid leukemia who previously were treated with DNA topoisomerase II inhibitors. In this study, we have identified a fusion partner of MLL in a 10-year-old female who developed therapy-related acute myeloid leukemia 17 months after treatment for Hodgkin’s disease. Leukemia cells of this patient had a t(11;17)(q23;q25), which involved MLL as demonstrated by Southern blot analysis. The partner gene was cloned from cDNA of the leukemia cells by use of a combination of adapter reverse transcriptase–PCR, rapid amplification of 5′ cDNA ends, and blast database analysis to identify expressed sequence tags. The full-length cDNA of 2.8 kb was found to be an additional member of the septin family, therefore it was named MSF (MLL septin-like fusion). Members of the septin family conserve the GTP binding domain, localize in the cytoplasm, and interact with cytoskeletal filaments. A major 4-kb transcript of MSF was expressed ubiquitously; a 1.7-kb transcript was found in most tissues. An additional 3-kb transcript was found only in hematopoietic tissues. By amplification with MLL exon 5 forward primer and reverse primers in MSF, the appropriately sized products were obtained. MSF is highly homologous to hCDCrel-1, which is a partner gene of MLL in leukemias with a t(11;22)(q23;q11.2). Further analysis of MSF may help to delineate the function of MLL partner genes in leukemia, particularly in therapy-related leukemia.

Enrichment for murine keratinocyte stem cells based on cell surface phenotype

The identification and physical isolation of epithelial stem cells is critical to our understanding of their growth regulation during homeostasis, wound healing, and carcinogenesis. These stem cells remain poorly characterized because of the absence of specific molecular markers that permit us to distinguish them from their progeny, the transit amplifying (TA) cells, which have a more restricted proliferative potential. Cell kinetic analyses have permitted the identification of murine keratinocyte stem cells (KSCs) as slowly cycling cells that retain [3H]thymidine ([3H]Tdr) label, termed label-retaining cells (LRCs), whereas TA cells are visualized as rapidly cycling cells after a single pulse of [3H]Tdr, termed pulse-labeled cells (PLCs). Here, we report on the successful separation of KSCs from TA cells through the combined use of in vivo cell kinetic analysis and fluorescence-activated cell sorting. Specifically, we demonstrate that murine dorsal keratinocytes characterized by their high levels of α6 integrin and low to undetectable expression of the transferrin receptor (CD71) termed α6briCD71dim cells, are enriched for epithelial stem cells because they represent a minor (≈8%) and quiescent subpopulation of small blast-like cells, with a high nuclear:cytoplasmic ratio, containing ≈70% of label-retaining cells, the latter being a well documented characteristic of stem cells. Conversely, TA cells could be enriched in a phenotypically distinct subpopulation termed α6briCD71bri, representing the majority (≈60%) of basal keratinocytes that are actively cycling, and importantly contain ≈70% of [3H]Tdr pulse-labeled cells. Importantly, immunostaining of dorsal skin revealed the presence of CD71dim cells in the hair follicle bulge region, a well documented location for KSCs.

Y chromosomal fertility factors kl-2 and kl-3 of Drosophila melanogaster encode dynein heavy chain polypeptides

The molecular identity and function of the Drosophila melanogaster Y-linked fertility factors have long eluded researchers. Although the D. melanogaster genome sequence was recently completed, the fertility factors still were not identified, in part because of low cloning efficiency of heterochromatic Y sequences. Here we report a method for iterative blast searching to assemble heterochromatic genes from shotgun assemblies, and we successfully identify kl-2 and kl-3 as 1β- and γ-dynein heavy chains, respectively. Our conclusions are supported by formal genetics with X-Y translocation lines. Reverse transcription–PCR was successful in linking together unmapped sequence fragments from the whole-genome shotgun assembly, although some sequences were missing altogether from the shotgun effort and had to be generated de novo. We also found a previously undescribed Y gene, polycystine-related (PRY). The closest paralogs of kl-2, kl-3, and PRY (and also of kl-5) are autosomal and not X-linked, suggesting that the evolution of the Drosophila Y chromosome has been driven by an accumulation of male-related genes arising de novo from the autosomes.

Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships

Pairwise sequence comparison methods have been assessed using proteins whose relationships are known reliably from their structures and functions, as described in the scop database [Murzin, A. G., Brenner, S. E., Hubbard, T. & Chothia C. (1995) J. Mol. Biol. 247, 536–540]. The evaluation tested the programs blast [Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403–410], wu-blast2 [Altschul, S. F. & Gish, W. (1996) Methods Enzymol. 266, 460–480], fasta [Pearson, W. R. & Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444–2448], and ssearch [Smith, T. F. & Waterman, M. S. (1981) J. Mol. Biol. 147, 195–197] and their scoring schemes. The error rate of all algorithms is greatly reduced by using statistical scores to evaluate matches rather than percentage identity or raw scores. The E-value statistical scores of ssearch and fasta are reliable: the number of false positives found in our tests agrees well with the scores reported. However, the P-values reported by blast and wu-blast2 exaggerate significance by orders of magnitude. ssearch, fasta ktup = 1, and wu-blast2 perform best, and they are capable of detecting almost all relationships between proteins whose sequence identities are >30%. For more distantly related proteins, they do much less well; only one-half of the relationships between proteins with 20–30% identity are found. Because many homologs have low sequence similarity, most distant relationships cannot be detected by any pairwise comparison method; however, those which are identified may be used with confidence.

Lineage commitment in the progeny of murine hematopoietic preprogenitor cells: Influence of thrombopoietin and interleukin 5

Normal mouse marrow cells were stimulated by stem cell factor (SCF) to form dispersed or multicentric blast colonies containing progenitor cells committed to various hematopoietic lineages. Combination of the eosinophil-specific regulator interleukin 5 with SCF increased the frequency of colonies containing eosinophil-committed progenitor cells with multicentric but not dispersed blast colonies. Combination of thrombopoietin with SCF increased the frequency of colonies containing megakaryocyte-committed progenitor cells with both types of blast colony. Neither interleukin 5 nor thrombopoietin significantly altered the number or total cell content of blast colonies or progenitor cell numbers in blast colonies from those stimulated by SCF alone. No correlation was observed between total progenitor cell content and the presence or absence of either eosinophil or megakaryocyte progenitors in either type of blast colony. The data argue against a random process as being responsible for the formation of particular committed progenitor cells or the possibility that lineage-specific regulators merely enhance survival of such committed progenitor cells formed in developing blast colonies.

Gene2EST: a BLAST2 server for searching expressed sequence tag (EST) databases with eukaryotic gene-sized queries

Expressed sequence tags (ESTs) are randomly sequenced cDNA clones. Currently, nearly 3 million human and 2 million mouse ESTs provide valuable resources that enable researchers to investigate the products of gene expression. The EST databases have proven to be useful tools for detecting homologous genes, for exon mapping, revealing differential splicing, etc. With the increasing availability of large amounts of poorly characterised eukaryotic (notably human) genomic sequence, ESTs have now become a vital tool for gene identification, sometimes yielding the only unambiguous evidence for the existence of a gene expression product. However, BLAST-based Web servers available to the general user have not kept pace with these developments and do not provide appropriate tools for querying EST databases with large highly spliced genes, often spanning 50 000–100 000 bases or more. Here we describe Gene2EST (http://woody.embl-heidelberg.de/gene2est/), a server that brings together a set of tools enabling efficient retrieval of ESTs matching large DNA queries and their subsequent analysis. RepeatMasker is used to mask dispersed repetitive sequences (such as Alu elements) in the query, BLAST2 for searching EST databases and Artemis for graphical display of the findings. Gene2EST combines these components into a Web resource targeted at the researcher who wishes to study one or a few genes to a high level of detail.

RISSC: a novel database for ribosomal 16S–23S RNA genes spacer regions

A novel database, under the acronym RISSC (Ribosomal Intergenic Spacer Sequence Collection), has been created. It compiles more than 1600 entries of edited DNA sequence data from the 16S–23S ribosomal spacers present in most prokaryotes and organelles (e.g. mitochondria and chloroplasts) and is accessible through the Internet (http://ulises.umh.es/RISSC), where systematic searches for specific words can be conducted, as well as BLAST-type sequence searches. Additionally, a characteristic feature of this region, the presence/absence and nature of tRNA genes within the spacer, is included in all the entries, even when not previously indicated in the original database. All these combined features could provide a useful documen­tation tool for studies on evolution, identification, typing and strain characterization, among others.

The EMBL nucleotide sequence database

The EMBL Nucleotide Sequence Database (http://www.ebi.ac.uk/embl/) is maintained at the European Bioinformatics Institute (EBI) in an international collaboration with the DNA Data Bank of Japan (DDBJ) and GenBank at the NCBI (USA). Data is exchanged amongst the collaborating databases on a daily basis. The major contributors to the EMBL database are individual authors and genome project groups. Webin is the preferred web-based submission system for individual submitters, whilst automatic procedures allow incorporation of sequence data from large-scale genome sequencing centres and from the European Patent Office (EPO). Database releases are produced quarterly. Network services allow free access to the most up-to-date data collection via ftp, email and World Wide Web interfaces. EBI’s Sequence Retrieval System (SRS), a network browser for databanks in molecular biology, integrates and links the main nucleotide and protein databases plus many specialized databases. For sequence similarity searching a variety of tools (e.g. Blitz, Fasta, BLAST) are available which allow external users to compare their own sequences against the latest data in the EMBL Nucleotide Sequence Database and SWISS-PROT.

A rapid classification protocol for the CATH Domain Database to support structural genomics

In order to support the structural genomic initiatives, both by rapidly classifying newly determined structures and by suggesting suitable targets for structure determination, we have recently developed several new protocols for classifying structures in the CATH domain database (http://www.biochem.ucl.ac.uk/bsm/cath). These aim to increase the speed of classification of new structures using fast algorithms for structure comparison (GRATH) and to improve the sensitivity in recognising distant structural relatives by incorporating sequence information from relatives in the genomes (DomainFinder). In order to ensure the integrity of the database given the expected increase in data, the CATH Protein Family Database (CATH-PFDB), which currently includes 25 320 structural domains and a further 160 000 sequence relatives has now been installed in a relational ORACLE database. This was essential for developing more rigorous validation procedures and for allowing efficient querying of the database, particularly for genome analysis. The associated Dictionary of Homologous Superfamilies [Bray,J.E., Todd,A.E., Pearl,F.M.G., Thornton,J.M. and Orengo,C.A. (2000) Protein Eng., 13, 153–165], which provides multiple structural alignments and functional information to assist in assigning new relatives, has also been expanded recently and now includes information for 903 homo­logous superfamilies. In order to improve coverage of known structures, preliminary classification levels are now provided for new structures at interim stages in the classification protocol. Since a large proportion of new structures can be rapidly classified using profile-based sequence analysis [e.g. PSI-BLAST: Altschul,S.F., Madden,T.L., Schaffer,A.A., Zhang,J., Zhang,Z., Miller,W. and Lipman,D.J. (1997) Nucleic Acids Res., 25, 3389–3402], this provides preliminary classification for easily recognisable homologues, which in the latest release of CATH (version 1.7) represented nearly three-quarters of the non-identical structures.

Database resources of the National Center for Biotechnology Information

In addition to maintaining the GenBank® nucleic acid sequence database, the National Center for Biotechnology Information (NCBI) provides data analysis and retrieval resources that operate on the data in GenBank and a variety of other biological data made available through NCBI’s Web site. NCBI data retrieval resources include Entrez, PubMed, LocusLink and the Taxonomy Browser. Data analysis resources include BLAST, Electronic PCR, OrfFinder, RefSeq, UniGene, HomoloGene, Database of Single Nucleotide Polymorphisms (dbSNP), Human Genome Sequencing, Human MapViewer, GeneMap’99, Human–Mouse Homology Map, Cancer Chromosome Aberration Project (CCAP), Entrez Genomes, Clusters of Orthologous Groups (COGs) database, Retroviral Genotyping Tools, Cancer Genome Anatomy Project (CGAP), SAGEmap, Gene Expression Omnibus (GEO), Online Mendelian Inheri­tance in Man (OMIM), the Molecular Modeling Database (MMDB) and the Conserved Domain Database (CDD). Augmenting many of the Web applications are custom implementations of the BLAST program optimized to search specialized data sets. All of the resources can be accessed through the NCBI home page at: http://www.ncbi.nlm.nih.gov.

Mendel-GFDb and Mendel-ESTS: databases of plant gene families and ESTs annotated with gene family numbers and gene family names

There is no control over the information provided with sequences when they are deposited in the sequence databases. Consequently mistakes can seed the incorrect annotation of other sequences. Grouping genes into families and applying controlled annotation overcomes the problems of incorrect annotation associated with individual sequences. Two databases (http://www.mendel.ac.uk) were created to apply controlled annotation to plant genes and plant ESTs: Mendel-GFDb is a database of plant protein (gene) families based on gapped-BLAST analysis of all sequences in the SWISS-PROT family of databases. Sequences are aligned (ClustalW) and identical and similar residues shaded. The families are visually curated to ensure that one or more criteria, for example overall relatedness and/or domain similarity relate all sequences within a family. Sequence families are assigned a ‘Gene Family Number’ and a unified description is developed which best describes the family and its members. If authority exists the gene family is assigned a ‘Gene Family Name’. This information is placed in Mendel-GFDb. Mendel-ESTS is primarily a database of plant ESTs, which have been compared to Mendel-GFDb, completely sequenced genomes and domain databases. This approach associated ESTs with individual sequences and the controlled annotation of gene families and protein domains; the information being placed in Mendel-ESTS. The controlled annotation applied to genes and ESTs provides a basis from which a plant transcription database can be developed.

ParAlign: a parallel sequence alignment algorithm for rapid and sensitive database searches

There is a need for faster and more sensitive algorithms for sequence similarity searching in view of the rapidly increasing amounts of genomic sequence data available. Parallel processing capabilities in the form of the single instruction, multiple data (SIMD) technology are now available in common microprocessors and enable a single microprocessor to perform many operations in parallel. The ParAlign algorithm has been specifically designed to take advantage of this technology. The new algorithm initially exploits parallelism to perform a very rapid computation of the exact optimal ungapped alignment score for all diagonals in the alignment matrix. Then, a novel heuristic is employed to compute an approximate score of a gapped alignment by combining the scores of several diagonals. This approximate score is used to select the most interesting database sequences for a subsequent Smith–Waterman alignment, which is also parallelised. The resulting method represents a substantial improvement compared to existing heuristics. The sensitivity and specificity of ParAlign was found to be as good as Smith–Waterman implementations when the same method for computing the statistical significance of the matches was used. In terms of speed, only the significantly less sensitive NCBI BLAST 2 program was found to outperform the new approach. Online searches are available at http://dna.uio.no/search/

PLS1, a gene encoding a tetraspanin-like protein, is required for penetration of rice leaf by the fungal pathogen Magnaporthe grisea

We describe in this study punchless, a nonpathogenic mutant from the rice blast fungus M. grisea, obtained by plasmid-mediated insertional mutagenesis. As do most fungal plant pathogens, M. grisea differentiates an infection structure specialized for host penetration called the appressorium. We show that punchless differentiates appressoria that fail to breach either the leaf epidermis or artificial membranes such as cellophane. Cytological analysis of punchless appressoria shows that they have a cellular structure, turgor, and glycogen content similar to those of wild type before penetration, but that they are unable to differentiate penetration pegs. The inactivated gene, PLS1, encodes a putative integral membrane protein of 225 aa (Pls1p). A functional Pls1p-green fluorescent protein fusion protein was detected only in appressoria and was localized in plasma membranes and vacuoles. Pls1p is structurally related to the tetraspanin family. In animals, these proteins are components of membrane signaling complexes controlling cell differentiation, motility, and adhesion. We conclude that PLS1 controls an appressorial function essential for the penetration of the fungus into host leaves.

A unified statistical framework for sequence comparison and structure comparison

We present an approach for assessing the significance of sequence and structure comparisons by using nearly identical statistical formalisms for both sequence and structure. Doing so involves an all-vs.-all comparison of protein domains [taken here from the Structural Classification of Proteins (scop) database] and then fitting a simple distribution function to the observed scores. By using this distribution, we can attach a statistical significance to each comparison score in the form of a P value, the probability that a better score would occur by chance. As expected, we find that the scores for sequence matching follow an extreme-value distribution. The agreement, moreover, between the P values that we derive from this distribution and those reported by standard programs (e.g., blast and fasta validates our approach. Structure comparison scores also follow an extreme-value distribution when the statistics are expressed in terms of a structural alignment score (essentially the sum of reciprocated distances between aligned atoms minus gap penalties). We find that the traditional metric of structural similarity, the rms deviation in atom positions after fitting aligned atoms, follows a different distribution of scores and does not perform as well as the structural alignment score. Comparison of the sequence and structure statistics for pairs of proteins known to be related distantly shows that structural comparison is able to detect approximately twice as many distant relationships as sequence comparison at the same error rate. The comparison also indicates that there are very few pairs with significant similarity in terms of sequence but not structure whereas many pairs have significant similarity in terms of structure but not sequence.