Generation of a high-density rat EST map

Detection of low-abundant novel transcripts in mouse hematopoietic stem cells

Gene expression profiles of hematopoietic stem cells (HSCs) provide clues for understanding molecular mechanisms of HSC behavior, including self-renewal and differentiation. We took advantage of serial analysis of gene expression (SAGE) to identify medium- and low-abundant transcripts expressed in HSCs/hematopoietic progenitor cells (HPCs). Among a total of 31,380 unique transcripts, 17,326 (55%) correspond to known genes and, 14,054 (45%) are low-copy transcripts that have no matches to currently known genes. Among the former class, 3,899 (23%) were alternatively spliced transcripts and 3,754 (22%) represent anti-sense transcripts from known genes. Mapping of the SAGE tags to the mouse genome showed that differences in gene expression exist among chromosomes. In addition, comparison of the HSCs/HPCs SAGE data to that of myeloid progenitor cells revealed that massive genetic reprogramming occurs in hematopoietic cell differentiation. Our results demonstrate a previously unrecognized complexity of gene expression in HSCs/HPCs, and indicate the need for further efforts to fully identify and characterize the transcripts expressed in this cell type.

Annotation of expressed sequence tags for the East African cichlid fish Astatotilapia burtoni and evolutionary analyses of cichlid ORFs

Background

The cichlid fishes in general, and the exceptionally diverse East African haplochromine cichlids in particular, are famous examples of adaptive radiation and explosive speciation. Here we report the collection and annotation of more than 12,000 expressed sequence tags (ESTs) generated from three different cDNA libraries obtained from the East African haplochromine cichlid species Astatotilapia burtoni and Metriaclima zebra.

Results

We first annotated more than 12,000 newly generated cichlid ESTs using the Gene Ontology classification system. For evolutionary analyses, we combined these ESTs with all available sequence data for haplochromine cichlids, which resulted in a total of more than 45,000 ESTs. The ESTs represent a broad range of molecular functions and biological processes. We compared the haplochromine ESTs to sequence data from those available for other fish model systems such as pufferfish (Takifugu rubripes and Tetraodon nigroviridis), trout, and zebrafish. We characterized genes that show a faster or slower rate of base substitutions in haplochromine cichlids compared to other fish species, as this is indicative of a relaxed or reinforced selection regime. Four of these genes showed the signature of positive selection as revealed by calculating K_a/K_s ratios.

Conclusion

About 22% of the surveyed ESTs were found to have cichlid specific rate differences suggesting that these genes might play a role in lineage specific characteristics of cichlids. We also conclude that the four genes with a K_a/K_s ratio greater than one appear as good candidate genes for further work on the genetic basis of evolutionary success of haplochromine cichlid fishes.

Characterization of a new gene WX2 in Toxoplasma gondii

Using hybridization techniques, we prepared the monoclonal antibody (Mab) 7C3-C3 against Toxoplasma gondii. The protection tests showed that the protein (Mab7C3-C3) inhibited the invasion and proliferation of T. gondii RH strain in HeLa cells. The passive transfer test indicated that the antibody significantly prolonged the survival time of the challenged mice. It was also shown that the antibody could be used for the detection of the circulating antigen of T. gondii. After immunoscreening the T. gondii tachyzoite cDNA library with Mab7C3-C3, a new gene wx2 of T. gondii was obtained. Immunofluorescence analysis showed that the WX2 protein was located on the membrane of the parasite. Nucleotide sequence comparison showed 28% identity to the calcium channel alpha-1E unit and shared with the surface antigen related sequence in some conservative residues. However, no match was found in protein databases. Therefore, it was an unknown gene in T. gondii encoding a functional protein on the membrane of T. gondii. Because it has been shown to have a partial protective effect against T. gondii infection and is released as a circulating antigen, it could be a candidate molecule for vaccine or a novel target for new drugs.

Nicotine causes age-dependent changes in gene expression in the adolescent female rat brain

Humans often start smoking during adolescence. Recent results suggest that rodents may also be particularly vulnerable to nicotine dependence during adolescence. We examined the effect of chronic nicotine exposure on gene expression profiles during adolescence in female rats, who were dosed with nicotine (and control animals were dosed with saline) via subcutaneously implanted osmotic minipumps. Brain samples were collected at four ages: before puberty (postnatal day 25), at about the time of puberty in females (postnatal day 35), and after puberty (postnatal days 45 and 55). The expression of 7931 genes in three brain areas was measured using DNA microarrays. Quantitative RT-PCR was also employed to confirm the expression patterns of selected genes. We used a novel clustering technique (principal cluster analysis) to classify 162 nicotine-regulated genes into five clusters, of which only one (cluster A) showed similar patterns of gene expression across all three brain areas (ventral striatum, prefrontal cortex, and hippocampus). Three clusters of genes (A, B, and C) showed dramatic peaks in their nicotine responses at the same age (p35). The other two clusters (D1 and D2) showed smaller peaks and/or valleys in their nicotine responses at p35 and p45. Thus, the age of maximal gene expression response to nicotine in female rats corresponds approximately to the age of maximal behavioral response and the age of puberty.

Development of a deer mouse whole-genome radiation hybrid panel and comparative mapping of Mus chromosome 11 loci

A 5000-rad whole-genome radiation hybrid cell panel (BW5000) was developed for mapping the deer mouse (Peromyscus maniculatus bairdii) genome. The panel consists of 103 cell lines and has an estimated marker retention frequency of 63.9% (range, 28%-88%) based on PCR typing of 30 Type I (coding gene) and 25 Type II (microsatellite) markers. Using the composite Mus map, Type I markers were selected from six Mus chromosomes, 22 of which are on Mus Chr 11. Fifteen of the Mus Chr 11 markers were simultaneously mapped on an interspecific (P. maniculatus x P. polionotus) backcross panel to test the utility of the radiation hybrid panel, create a framework map, and help establish gene order. The radiation hybrids have effectively detected linkage in the deer mouse genome between markers as far apart as 6.7 cM and resolved markers that are, in the Mus genome, as close as 0.2 Mb. Combined results from both panels have indicated a high degree of gene order conservation of the telomeric 64 cM of Mus Chr 11 in the deer mouse genome. The remaining centromeric portion also shows gene order conservation with the deer mouse but as a separate linkage group. This indicates a translocation of that portion of Mus Chr 11 in P. maniculatus and is consistent with rearrangement breakpoints observed between Mus and other mammalian genomes, including rat and human. Furthermore, this separate linkage group is likely to reside in a chromosomal region of inversion polymorphism between P. maniculatus and P. polionotus.

Chromosome bin map of expressed sequence tags in homoeologous group 1 of hexaploid wheat and homoeology with rice and Arabidopsis

A total of 944 expressed sequence tags (ESTs) generated 2212 EST loci mapped to homoeologous group 1 chromosomes in hexaploid wheat (Triticum aestivum L.). EST deletion maps and the consensus map of group 1 chromosomes were constructed to show EST distribution. EST loci were unevenly distributed among chromosomes 1A, 1B, and 1D with 660, 826, and 726, respectively. The number of EST loci was greater on the long arms than on the short arms for all three chromosomes. The distribution of ESTs along chromosome arms was nonrandom with EST clusters occurring in the distal regions of short arms and middle regions of long arms. Duplications of group 1 ESTs in other homoeologous groups occurred at a rate of 35.5%. Seventy-five percent of wheat chromosome 1 ESTs had significant matches with rice sequences (E < or = e(-10)), where large regions of conservation occurred between wheat consensus chromosome 1 and rice chromosome 5 and between the proximal portion of the long arm of wheat consensus chromosome 1 and rice chromosome 10. Only 9.5% of group 1 ESTs showed significant matches to Arabidopsis genome sequences. The results presented are useful for gene mapping and evolutionary and comparative genomics of grasses.

Single nucleotide polymorphisms associated with rat expressed sequences

Single nucleotide polymorphisms (SNPs) are the most common source of genetic variation in populations and are thus most likely to account for the majority of phenotypic and behavioral differences between individuals or strains. Although the rat is extensively studied for the latter, data on naturally occurring polymorphisms are mostly lacking. We have used publicly available sequences consisting of whole-genome shotgun (WGS), expressed sequence tag (EST), and mRNA data as a source for the in silico identification of SNPs in gene-coding regions and have identified a large collection of 33,305 high-quality candidate SNPs. Experimental verification of 471 candidate SNPs using a limited set of rat isolates revealed a confirmation rate of approximately 50%. Although the majority of SNPs were identified between Sprague-Dawley (EST data) and Brown Norway (WGS data) strains, we found that 66% of the verified variations are common among different rat strains. All SNPs were extensively annotated, including chromosomal and genetic map information, and nonsynonymous SNPs were analyzed by SIFT and PolyPhen prediction programs for their potential deleterious effect on protein function. Interestingly, we retrieved three SNPs from the database that result in the introduction of a premature stop codon and that could be confirmed experimentally. Two of these "in silico-identified knockouts" reside in interesting QTL regions. Data are publicly available via a Web interface (http://cascad.niob.knaw.nl), allowing simple and advanced search queries.

High-density rat radiation hybrid maps containing over 24,000 SSLPs, genes, and ESTs provide a direct link to the rat genome sequence

The laboratory rat is a major model organism for systems biology. To complement the cornucopia of physiological and pharmacological data generated in the rat, a large genomic toolset has been developed, culminating in the release of the rat draft genome sequence. The rat draft sequence used a variety of assembly packages, as well as data from the Radiation Hybrid (RH) map of the rat as part of their validation. As part of the Rat Genome Project, we have been building a high-density RH map to facilitate data integration from multiple maps and now to help validate the genome assembly. By incorporating vectors from our lab and several other labs, we have doubled the number of simple sequence length polymorphisms (SSLPs), genes, expressed sequence tags (ESTs), and sequence-tagged sites (STSs) compared to any other genome-wide rat map, a total of 24,437 elements. During the process, we also identified a novel approach for integrating the RH placement results from multiple maps. This new integrated RH map contains approximately 10 RH-mapped elements per Mb on the genome assembly, enabling the RH maps to serve as a scaffold for a variety of data visualization tools.

Integrative genomics: in silico coupling of rat physiology and complex traits with mouse and human data

Integration of the large variety of genome maps from several organisms provides the mechanism by which physiological knowledge obtained in model systems such as the rat can be projected onto the human genome to further the research on human disease. The release of the rat genome sequence provides new information for studies using the rat model and is a key reference against which existing and new rat physiological results can be aligned. Previously, we described comparative maps of the rat, mouse, and human based on EST sequence comparisons combined with radiation hybrid maps. Here, we use new data and introduce the Integrated Genomics Environment, an extensive database of curated and integrated maps, markers, and physiological results. These results are integrated by using VCMapview, a java-based map integration and visualization tool. This unique environment allows researchers to relate results from cytogenetic, genetic, and radiation hybrid studies to the genome sequence and compare regions of interest between human, mouse, and rat. Integrating rat physiology with mouse genetics and clinical results from human by using the respective genomes provides a novel route to capitalize on comparative genomics and the strengths of model organism biology.

A comprehensive nonredundant expressed sequence tag collection for the developing Rattus norvegicus heart

Congenital heart defects affect ∼1,000,000 people in the United States, with 40,000 new births contributing to that number every year. A large percentage of these defects can be attributed to septal defects. We assembled a nonredundant collection of over 12,000 expressed sequence tags (ESTs) from a total of 30,000 ESTs, with the ultimate goal of identifying spatially and/or temporally regulated genes during heart septation. These ESTs were compiled from nonnormalized, normalized, and serially subtracted cDNA libraries derived from two sets of tissue samples. The first includes microdissected rat hearts from embryonic (E) days E13, E15, and E16.5–E18.5 and adult heart. The second includes hearts from embryonic days E17, E19, and E21 and postnatal (P) days P1, P12, P74, and P200. Over 6,000 novel ESTs were identified in the libraries derived from these two sets of tissues, all of which have been contributed to the NCBI rat UniGene collection. It is anticipated that such EST and cDNA clone resources will prove invaluable to gene expression studies aimed at the understanding of the molecular mechanisms underlying heart septation defects.

Integration of the rat recombination and EST maps in the rat genomic sequence and comparative mapping analysis with the mouse genome

Inbred strains of the laboratory rat are widely used for identifying genetic regions involved in the control of complex quantitative phenotypes of biomedical importance. The draft genomic sequence of the rat now provides essential information for annotating rat quantitative trait locus (QTL) maps. Following the survey of unique rat microsatellite (11,585 including 1648 new markers) and EST (10,067) markers currently available, we have incorporated a selection of 7952 rat EST sequences in an improved version of the integrated linkage-radiation hybrid map of the rat containing 2058 microsatellite markers which provided over 10,000 potential anchor points between rat QTL and the genomic sequence of the rat. A total of 996 genetic positions were resolved (avg. spacing 1.77 cM) in a single large intercross and anchored in the rat genomic sequence (avg. spacing 1.62 Mb). Comparative genome maps between rat and mouse were constructed by successful computational alignment of 6108 mapped rat ESTs in the mouse genome. The integration of rat linkage maps in the draft genomic sequence of the rat and that of other species represents an essential step for translating rat QTL intervals into human chromosomal targets.

Integrated and sequence-ordered BAC- and YAC-based physical maps for the rat genome

As part of the effort to sequence the genome of Rattus norvegicus, we constructed a physical map comprised of fingerprinted bacterial artificial chromosome (BAC) clones from the CHORI-230 BAC library. These BAC clones provide approximately 13-fold redundant coverage of the genome and have been assembled into 376 fingerprint contigs. A yeast artificial chromosome (YAC) map was also constructed and aligned with the BAC map via fingerprinted BAC and P1 artificial chromosome clones (PACs) sharing interspersed repetitive sequence markers with the YAC-based physical map. We have annotated 95% of the fingerprint map clones in contigs with coordinates on the version 3.1 rat genome sequence assembly, using BAC-end sequences and in silico mapping methods. These coordinates have allowed anchoring 358 of the 376 fingerprint map contigs onto the sequence assembly. Of these, 324 contigs are anchored to rat genome sequences localized to chromosomes, and 34 contigs are anchored to unlocalized portions of the rat sequence assembly. The remaining 18 contigs, containing 54 clones, still require placement. The fingerprint map is a high-resolution integrative data resource that provides genome-ordered associations among BAC, YAC, and PAC clones and the assembled sequence of the rat genome.

Physical and linkage mapping of mammary-derived expressed sequence tags in cattle

This study describes the physical and linkage mapping of 42 gene-associated markers developed from mammary gland-derived expressed sequence tags to the cattle genome. Of the markers, 25 were placed on the USDA reference linkage map and 37 were positioned on the Roslin 3000-rad radiation hybrid (RH) map, with 20 assignments shared between the maps. Although no novel regions of conserved synteny between the cattle and the human genomes were identified, the coverage was extended for bovine chromosomes 3, 7, 15, and 29 compared with previously published comparative maps between human and bovine genomes. Overall, these data improve the resolution of the human-bovine comparative maps and will assist future efforts to integrate bovine RH and linkage map data.

FTIR spectroscopy demonstrates biochemical differences in mammalian cell cultures at different growth stages

We have observed differences in the infrared spectra of viable fibroblast cells depending on whether the cells were in the exponential (proliferating) or plateau (nonproliferating) phase of growth. The spectral changes were observed even after correcting for cell number and volume, ruling out these trivial explanations. Several of the changes occurred for both transformed and normal cell lines and were greater for the normal cell line. The biochemical basis of the spectral changes was estimated by fitting the cell spectra to a linear superposition of spectra for the major biochemical components of mammalian cells (DNA, RNA, protein, lipids, and glycogen). The ratios of RNA/lipid and protein/lipid increased when the cells were in the exponential phase compared to the plateau phase of growth. The fits of cell spectra to individual biochemical components also demonstrated that DNA is a relatively minor spectroscopic component as would be expected biochemically. Contrary to other reports in the literature, our data demonstrate that determining DNA content or structure using Fourier transform infrared spectroscopy data is difficult due to the relatively small amount of DNA and the overlap of DNA bands with the absorption bands of other biochemical components.

Gene discovery in the apicomplexa as revealed by EST sequencing and assembly of a comparative gene database

Large-scale EST sequencing projects for several important parasites within the phylum Apicomplexa were undertaken for the purpose of gene discovery. Included were several parasites of medical importance (Plasmodium falciparum, Toxoplasma gondii) and others of veterinary importance (Eimeria tenella, Sarcocystis neurona, and Neospora caninum). A total of 55192 ESTs, deposited into dbEST/GenBank, were included in the analyses. The resulting sequences have been clustered into nonredundant gene assemblies and deposited into a relational database that supports a variety of sequence and text searches. This database has been used to compare the gene assemblies using BLAST similarity comparisons to the public protein databases to identify putative genes. Of these new entries, approximately 15%-20% represent putative homologs with a conservative cutoff of p < 10(-9), thus identifying many conserved genes that are likely to share common functions with other well-studied organisms. Gene assemblies were also used to identify strain polymorphisms, examine stage-specific expression, and identify gene families. An interesting class of genes that are confined to members of this phylum and not shared by plants, animals, or fungi, was identified. These genes likely mediate the novel biological features of members of the Apicomplexa and hence offer great potential for biological investigation and as possible therapeutic targets.

[beta]-defensins in lung host defense

Host defenses at the mucosal surface of the airways evolved to present many layers of protection against inhaled microbes. Normally, the intrapulmonary airways are sterile. Airway secretions contain numerous factors with antimicrobial activity that contribute to innate defenses. Many protein and peptide components exert bacteriostatic or bacteriocidal effects against a wide variety of organisms and may act in synergistic or additive combinations. The beta-defensins are a relatively recently described family of peptide antimicrobials that are widely expressed at mucosal surfaces, including airway and submucosal gland epithelia. These small cationic peptides are products of individual genes that exhibit broad-spectrum activity against bacteria, fungi, and some enveloped viruses. Their expression in airway epithelia may be constitutive or inducible by bacterial products or pro-inflammatory cytokines. beta-defensins also act as chemokines for adaptive immune cells, including immature dendritic cells and T cells via the CCR6 receptor, and provide a link between innate and adaptive immunity. Alterations in the function of the beta-defensins may contribute to disease states. Here we review much of the biology of the beta-defensins, including gene discovery, genomic organization, molecular structure, regulation of expression, and function.

The LN54 radiation hybrid map of zebrafish expressed sequences

To increase the density of a gene map of the zebrafish, Danio rerio, we have placed 3119 expressed sequence tags (ESTs) and cDNA sequences on the LN54 radiation hybrid (RH) panel. The ESTs and genes mapped here join 748 SSLp markers and 459 previously mapped genes and ESTs, bringing the total number of markers on the LN54 RH panel to 4226. Addition of these new markers brings the total LN54 map size to 14,372 cR, with 118 kb/cR. The distribution of ESTs according to linkage groups shows relatively little variation (minimum, 73; maximum, 201). This observation, combined with a relatively uniform size for zebrafish chromosomes, as previously indicated by karyotyping, indicates that there are no especially gene-rich or gene-poor chromosomes in this species. We developed an algorithm to provide a semiautomatic method for the selection of additional framework markers for the LN54 map. This algorithm increased the total number of framework markers to 1150 and permitted the mapping of a high percentage of sequences that could not be placed on a previous version of the LN54 map. The increased concentration of expressed sequences on the LN54 map of the zebrafish genome will facilitate the molecular characterization of mutations in this species.

Automated construction of high-density comparative maps between rat, human, and mouse

Animal models have been used primarily as surrogates for humans, having similar disease-based phenotypes. Genomic organization also tends to be conserved between species, leading to the generation of comparative genome maps. The emergence of radiation hybrid (RH) maps, coupled with the large numbers of available Expressed Sequence Tags (ESTs), has revolutionized the way comparative maps can be built. We used publicly available rat, mouse, and human data to identify genes and ESTs with interspecies sequence identity (homology), identified their UniGene relationships, and incorporated their RH map positions to build integrated comparative maps with >2100 homologous UniGenes mapped in more than one species (approximately 6% of all mammalian genes). The generation of these maps is iterative and labor intensive; therefore, we developed a series of computer tools (not described here) based on our algorithm that identifies anchors between species and produces printable and on-line clickable comparative maps that link to a wide variety of useful tools and databases. The maps were constructed using sequence-based comparisons, thus creating "hooks" for further sequence-based annotation of human, mouse, and rat sequences. Currently, this map enables investigators to link the physiology of the rat with the genetics of the mouse and the clinical significance of the human.

A high-resolution consensus linkage map of the rat, integrating radiation hybrid and genetic maps

We have constructed a high-resolution consensus genetic map of the rat in a single large intercross, which integrates 747 framework markers and 687 positions of our whole-genome radiation hybrid (RH) map of the rat. We selected 136 new gene markers from the GenBank database and assigned them either genetically or physically to rat chromosomes to evaluate the accuracy of the integrated linkage-RH maps in the localization of new markers and to enrich existing comparative mapping data. These markers and 631 D-Got- markers, which are physically mapped but still uncharacterized for evidence of polymorphism, were tested for allele variations in a panel of 16 rat strains commonly used in genetic studies. The consensus linkage map constructed in the GK x BN cross now comprises 1620 markers of various origins, defining 840 resolved genetic positions with an average spacing of 2.2 cM between adjacent loci, and includes 407 gene markers. This whole-genome genetic map will contribute to the advancement of genetic studies in the rat by incorporating gene/EST maps, physical mapping information, and sequence data generated in rat and other mammalian species into genetic intervals harboring disease susceptibility loci identified in rat models of human genetic disorders.

Evolution of mammalian genome organization inferred from comparative gene mapping

Comparative genome analyses, including chromosome painting in over 40 diverse mammalian species, ordered gene maps from several representatives of different mammalian and vertebrate orders, and large-scale sequencing of the human and mouse genomes are beginning to provide insight into the rates and patterns of chromosomal evolution on a whole-genome scale, as well as into the forces that have sculpted the genomes of extant mammalian species.