Transcriptional patterns in both host and bacterium underlie a daily rhythm of anatomical and metabolic change in a beneficial symbiosis.

Mechanisms for controlling symbiont populations are critical for maintaining the associations that exist between a host and its microbial partners. We describe here the transcriptional, metabolic, and ultrastructural characteristics of a diel rhythm that occurs in the symbiosis between the squid Euprymna scolopes and the luminous bacterium Vibrio fischeri. The rhythm is driven by the host's expulsion from its light-emitting organ of most of the symbiont population each day at dawn. The transcriptomes of both the host epithelium that supports the symbionts and the symbiont population itself were characterized and compared at four times over this daily cycle. The greatest fluctuation in gene expression of both partners occurred as the day began. Most notable was an up-regulation in the host of >50 cytoskeleton-related genes just before dawn and their subsequent down-regulation within 6 h. Examination of the epithelium by TEM revealed a corresponding restructuring, characterized by effacement and blebbing of its apical surface. After the dawn expulsion, the epithelium reestablished its polarity, and the residual symbionts began growing, repopulating the light organ. Analysis of the symbiont transcriptome suggested that the bacteria respond to the effacement by up-regulating genes associated with anaerobic respiration of glycerol; supporting this finding, lipid analysis of the symbionts' membranes indicated a direct incorporation of host-derived fatty acids. After 12 h, the metabolic signature of the symbiont population shifted to one characteristic of chitin fermentation, which continued until the following dawn. Thus, the persistent maintenance of the squid-vibrio symbiosis is tied to a dynamic diel rhythm that involves both partners.

Effects of colonization, luminescence, and autoinducer on host transcription during development of the squid-vibrio association.

The light-organ symbiosis between the squid Euprymna scolopes and the luminous bacterium Vibrio fischeri offers the opportunity to decipher the hour-by-hour events that occur during the natural colonization of an animal's epithelial surface by its microbial partners. To determine the genetic basis of these events, a glass-slide microarray was used to characterize the light-organ transcriptome of juvenile squid in response to the initiation of symbiosis. Patterns of gene expression were compared between animals not exposed to the symbiont, exposed to the wild-type symbiont, or exposed to a mutant symbiont defective in either of two key characters of this association: bacterial luminescence or autoinducer (AI) production. Hundreds of genes were differentially regulated as a result of symbiosis initiation, and a hierarchy existed in the magnitude of the host's response to three symbiont features: bacterial presence > luminescence > AI production. Putative host receptors for bacterial surface molecules known to induce squid development are up-regulated by symbiont light production, suggesting that bioluminescence plays a key role in preparing the host for bacteria-induced development. Further, because the transcriptional response of tissues exposed to AI in the natural context (i.e., with the symbionts) differed from that to AI alone, the presence of the bacteria potentiates the role of quorum signals in symbiosis. Comparison of these microarray data with those from other symbioses, such as germ-free/conventionalized mice and zebrafish, revealed a set of shared genes that may represent a core set of ancient host responses conserved throughout animal evolution.

Analysis of ESTs from Lutzomyia longipalpis sand flies and their contribution toward understanding the insect-parasite relationship.

An expressed sequence tag library has been generated from a sand fly vector of visceral leishmaniasis, Lutzomyia longipalpis. A normalized cDNA library was constructed from whole adults and 16,608 clones were sequenced from both ends and assembled into 10,203 contigs and singlets. Of these 58% showed significant similarity to known genes from other organisms, <4% were identical to described sand fly genes, and 42% had no match to any database sequence. Our analyses revealed putative proteins involved in the barrier function of the gut (peritrophins, microvillar proteins, glutamine synthase), digestive physiology (secreted and membrane-anchored hydrolytic enzymes), and the immune response (gram-negative binding proteins, thioester proteins, scavenger receptors, galectins, signaling pathway factors, caspases, serpins, and peroxidases). Sequence analysis of this transcriptome dataset has provided new insights into genes that might be associated with the response of the vector to the development of Leishmania.

An annotated cDNA library of juvenile Euprymna scolopes with and without colonization by the symbiont Vibrio fischeri.

BACKGROUND: Biologists are becoming increasingly aware that the interaction of animals, including humans, with their coevolved bacterial partners is essential for health. This growing awareness has been a driving force for the development of models for the study of beneficial animal-bacterial interactions. In the squid-vibrio model, symbiotic Vibrio fischeri induce dramatic developmental changes in the light organ of host Euprymna scolopes over the first hours to days of their partnership. We report here the creation of a juvenile light-organ specific EST database. RESULTS: We generated eleven cDNA libraries from the light organ of E. scolopes at developmentally significant time points with and without colonization by V. fischeri. Single pass 3' sequencing efforts generated 42,564 expressed sequence tags (ESTs) of which 35,421 passed our quality criteria and were then clustered via the UIcluster program into 13,962 nonredundant sequences. The cDNA clones representing these nonredundant sequences were sequenced from the 5' end of the vector and 58% of these resulting sequences overlapped significantly with the associated 3' sequence to generate 8,067 contigs with an average sequence length of 1,065 bp. All sequences were annotated with BLASTX (E-value < -03) and Gene Ontology (GO). CONCLUSION: Both the number of ESTs generated from each library and GO categorizations are reflective of the activity state of the light organ during these early stages of symbiosis. Future analyses of the sequences identified in these libraries promise to provide valuable information not only about pathways involved in colonization and early development of the squid light organ, but also about pathways conserved in response to bacterial colonization across the animal kingdom.

Identifying components of the NF-kappaB pathway in the beneficial Euprymna scolopes-Vibrio fischeri light organ symbiosis.

The Toll/NF-kappaB pathway is a common, evolutionarily conserved innate immune pathway that modulates the responses of animal cells to microbe-associated molecular patterns (MAMPs). Because MAMPs have been implicated as critical elements in the signaling of symbiont-induced development, an expressed sequence tag library from the juvenile light organ of Euprymna scolopes was used to identify members of the Toll/NF-kappaB pathway. Full-length transcripts were identified by using 5' and 3' RACE PCR. Seven transcripts critical for MAMP-induced triggering of the Toll/NF-kappaB phosphorylation cascade have been identified, including receptors, signal transducers, and a transcription factor. Further investigations should elucidate the role of the Toll/NF-kappaB pathway in the initiation of the beneficial symbiosis between E. scolopes and Vibrio fischeri.

Comprehensive in silico functional specification of mouse retina transcripts.

BACKGROUND: The retina is a well-defined portion of the central nervous system (CNS) that has been used as a model for CNS development and function studies. The full specification of transcripts in an individual tissue or cell type, like retina, can greatly aid the understanding of the control of cell differentiation and cell function. In this study, we have integrated computational bioinformatics and microarray experimental approaches to classify the tissue specificity and developmental distribution of mouse retina transcripts. RESULTS: We have classified a set of retina-specific genes using sequence-based screening integrated with computational and retina tissue-specific microarray approaches. 33,737 non-redundant sequences were identified as retina transcript clusters (RTCs) from more than 81,000 mouse retina ESTs. We estimate that about 19,000 to 20,000 genes might express in mouse retina from embryonic to adult stages. 39.1% of the RTCs are not covered by 60,770 RIKEN full-length cDNAs. Through comparison with 2 million mouse ESTs, spectra of neural, retinal, late-generated retinal, and photoreceptor -enriched RTCs have been generated. More than 70% of these RTCs have data from biological experiments confirming their tissue-specific expression pattern. The highest-grade retina-enriched pool covered almost all the known genes encoding proteins involved in photo-transduction. CONCLUSION: This study provides a comprehensive mouse retina transcript profile for further gene discovery in retina and suggests that tissue-specific transcripts contribute substantially to the whole transcriptome.

Differential expression of genes within the cochlea as defined by a custom mouse inner ear microarray.

Microarray analyses have contributed greatly to the rapid understanding of functional genomics through the identification of gene networks as well as gene discovery. To facilitate functional genomics of the inner ear, we have developed a mouse inner-ear-pertinent custom microarray chip (CMA-IE1). Nonredundant cDNA clones were obtained from two cDNA library resources: the RIKEN subtracted inner ear set and the NIH organ of Corti library. At least 2000 cDNAs unique to the inner ear were present on the chip. Comparisons were performed to examine the relative expression levels of these unique cDNAs within the organ of Corti, lateral wall, and spiral ganglion. Total RNA samples were obtained from the three cochlear-dissected fractions from adult CF-1 mice. The total RNA was linearly amplified, and a dendrimer-based system was utilized to enhance the hybridization signal. Differentially expressed genes were verified by comparison to known gene expression patterns in the cochlea or by correlation with genes and gene families deduced to be present in the three tissue types. Approximately 22-25% of the genes on the array had significant levels of expression. A number of differentially expressed genes were detected in each tissue fraction. These included genes with known functional roles, hypothetical genes, and various unknown or uncharacterized genes. Four of the differentially expressed genes found in the organ of Corti are linked to deafness loci. None of these are hypothetical or unknown genes.

The status, quality, and expansion of the NIH full-length cDNA project: the Mammalian Gene Collection (MGC).

The National Institutes of Health's Mammalian Gene Collection (MGC) project was designed to generate and sequence a publicly accessible cDNA resource containing a complete open reading frame (ORF) for every human and mouse gene. The project initially used a random strategy to select clones from a large number of cDNA libraries from diverse tissues. Candidate clones were chosen based on 5'-EST sequences, and then fully sequenced to high accuracy and analyzed by algorithms developed for this project. Currently, more than 11,000 human and 10,000 mouse genes are represented in MGC by at least one clone with a full ORF. The random selection approach is now reaching a saturation point, and a transition to protocols targeted at the missing transcripts is now required to complete the mouse and human collections. Comparison of the sequence of the MGC clones to reference genome sequences reveals that most cDNA clones are of very high sequence quality, although it is likely that some cDNAs may carry missense variants as a consequence of experimental artifact, such as PCR, cloning, or reverse transcriptase errors. Recently, a rat cDNA component was added to the project, and ongoing frog (Xenopus) and zebrafish (Danio) cDNA projects were expanded to take advantage of the high-throughput MGC pipeline.

PLET1 (C11orf34), a highly expressed and processed novel gene in pig and mouse placenta, is transcribed but poorly spliced in human.

Sequencing of porcine cDNAs identified a novel EST with high frequency in placenta tissue. Full-length PLET1 (placenta-expressed transcript 1, also called C11orf34) matched a mouse cDNA and many bovine and mouse ESTs but no human transcripts or ESTs. However, the porcine cDNA matched several putative exons within a human genomic DNA fragment on chromosome 11. This human locus is in a region of conserved synteny with pig chromosome 9, to which the porcine gene was subsequently mapped. RNA blot hybridization showed that this gene had high expression in porcine and mouse conceptus and throughout placenta development. In situ hybridization using mouse placenta showed PLET1 expression in trophoblast cells of the labyrinth, as well as in spongiotrophoblast and glycogen trophoblast cells. However, no expression of PLET1 was detected by RNA blot analysis of human placenta, although RT-PCR analysis detected very small amounts of partially spliced RNA that were significantly less abundant than the RNA levels in mouse placenta. Donor and acceptor splicing site sequences in the exons of the human gene are poorly conserved and may be the cause of inefficient splicing found specifically in human tissue. Our data correct GenomeScan annotation of this region of the human genome and describe functional gene discovery in mammals not recognized in human EST projects.

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.

Saci-1, -2, and -3 and Perere, four novel retrotransposons with high transcriptional activities from the human parasite Schistosoma mansoni.

Using the data set of 180,000 expressed sequence tags (ESTs) of the blood fluke Schistosoma mansoni generated recently by our group, we identified three novel long-terminal-repeat (LTR)- and one novel non-LTR-expressed retrotransposon, named Saci-1, -2, and -3 and Perere, respectively. Full-length sequences were reconstructed from ESTs and have deduced open reading frames (ORFs) with several uncorrupted features, characterizing them as possible active retrotransposons of different known transposon families. Alignment of reconstructed sequences to available preliminary genome sequence data confirmed the overall structure of the transposons. The frequency of sequenced transposon transcripts in cercariae was 14% of all transcripts from that stage, twofold higher than that in schistosomula and three- to fourfold higher than that in adults, eggs, miracidia, and germ balls. We show by Southern blot analysis, by EST annotation and tallying, and by counting transposon tags from a Serial Analysis of Gene Expression library, that the four novel retrotransposons exhibit a 10- to 30-fold lower copy number in the genome and a 4- to 200-fold-higher transcriptional rate per copy than the four previously described S. mansoni retrotransposons [corrected]. Such differences lead us to hypothesize that there are two different populations of retrotransposons in S. mansoni genome, occupying different niches in its ecology. Examples of retrotransposon fragment inserts were found into the 5' and 3' untranslated regions of four different S. mansoni target gene transcripts. The data presented here suggest a role for these elements in the dynamics of this complex human parasite genome.

Migration of the plastid genome to the nucleus in a peridinin dinoflagellate.

Dinoflagellate algae are important primary producers and of significant ecological and economic impact because of their ability to form "red tides". They are also models for evolutionary research because of an unparalleled ability to capture photosynthetic organelles (plastids) through endosymbiosis. The nature and extent of the plastid genome in the dominant perdinin-containing dinoflagellates remain, however, two of the most intriguing issues in plastid evolution. The plastid genome in these taxa is reduced to single-gene minicircles encoding an incomplete (until now 15) set of plastid proteins. The location of the remaining photosynthetic genes is unknown. We generated a data set of 6,480 unique expressed sequence tags (ESTs) from the toxic dinoflagellate Alexandrium tamarense (for details, see the Experimental Procedures in the Supplemental Data) to find the missing plastid genes and to understand the impact of endosymbiosis on genome evolution. Here we identify 48 of the non-minicircle-encoded photosynthetic genes in the nuclear genome of A. tamarense, accounting for the majority of the photosystem. Fifteen genes that are always found on the plastid genome of other algae and plants have been transferred to the nucleus in A. tamarense. The plastid-targeted genes have red and green algal origins. These results highlight the unique position of dinoflagellates as the champions of plastid gene transfer to the nucleus among photosynthetic eukaryotes.

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

Congenital heart defects affect approximately 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.

Large-scale gene discovery in human airway epithelia reveals novel transcripts.

The airway epithelium represents an important barrier between the host and the environment. It is a first site of contact with pathogens, particulates, and other stimuli, and has evolved the means to dynamically respond to these challenges. In an effort to define the transcript profile of airway epithelia, we created and sequenced cDNA libraries from cystic fibrosis (CF) and non-CF epithelia and from human lung tissue. Sequencing of these libraries produced approximately 53,000 3'-expressed sequence tags (3'-ESTs). From these, a nonredundant UniGene set of more than 19,000 sequences was generated. Despite the relatively small contribution of airway epithelia to the total mass of the lung, focused gene discovery in this tissue yielded novel results. The ESTs included several thousand transcripts (6,416) not previously identified from cDNA sequences as expressed in the lung. Among the abundant transcripts were several genes involved in host defense. Most importantly, the set also included 879 3'-ESTs that appear to be novel sequences not previously represented in the National Center for Biotechnology Information UniGene collection. This UniGene set should be useful for studies of pulmonary diseases involving the airway epithelium including cystic fibrosis, respiratory infections and asthma. It also provides a reagent for large-scale expression profiling.

Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni.

Schistosoma mansoni is the primary causative agent of schistosomiasis, which affects 200 million individuals in 74 countries. We generated 163,000 expressed-sequence tags (ESTs) from normalized cDNA libraries from six selected developmental stages of the parasite, resulting in 31,000 assembled sequences and 92% sampling of an estimated 14,000 gene complement. By analyzing automated Gene Ontology assignments, we provide a detailed view of important S. mansoni biological systems, including characterization of metazoa-specific and eukarya-conserved genes. Phylogenetic analysis suggests an early divergence from other metazoa. The data set provides insights into the molecular mechanisms of tissue organization, development, signaling, sexual dimorphism, host interactions and immune evasion and identifies novel proteins to be investigated as vaccine candidates and potential drug targets.

Generation and initial analysis of more than 15,000 full-length human and mouse cDNA sequences.

The National Institutes of Health Mammalian Gene Collection (MGC) Program is a multiinstitutional effort to identify and sequence a cDNA clone containing a complete ORF for each human and mouse gene. ESTs were generated from libraries enriched for full-length cDNAs and analyzed to identify candidate full-ORF clones, which then were sequenced to high accuracy. The MGC has currently sequenced and verified the full ORF for a nonredundant set of >9,000 human and >6,000 mouse genes. Candidate full-ORF clones for an additional 7,800 human and 3,500 mouse genes also have been identified. All MGC sequences and clones are available without restriction through public databases and clone distribution networks (see http:mgc.nci.nih.gov).

Annotated expressed sequence tags and cDNA microarrays for studies of brain and behavior in the honey bee.

To accelerate the molecular analysis of behavior in the honey bee (Apis mellifera), we created expressed sequence tag (EST) and cDNA microarray resources for the bee brain. Over 20,000 cDNA clones were partially sequenced from a normalized (and subsequently subtracted) library generated from adult A. mellifera brains. These sequences were processed to identify 15,311 high-quality ESTs representing 8912 putative transcripts. Putative transcripts were functionally annotated (using the Gene Ontology classification system) based on matching gene sequences in Drosophila melanogaster. The brain ESTs represent a broad range of molecular functions and biological processes, with neurobiological classifications particularly well represented. Roughly half of Drosophila genes currently implicated in synaptic transmission and/or behavior are represented in the Apis EST set. Of Apis sequences with open reading frames of at least 450 bp, 24% are highly diverged with no matches to known protein sequences. Additionally, over 100 Apis transcript sequences conserved with other organisms appear to have been lost from the Drosophila genome. DNA microarrays were fabricated with over 7000 EST cDNA clones putatively representing different transcripts. Using probe derived from single bee brain mRNA, microarrays detected gene expression for 90% of Apis cDNAs two standard deviations greater than exogenous control cDNAs. [The sequence data described in this paper have been submitted to Genbank data library under accession nos. BI502708-BI517278. The sequences are also available at http://titan.biotec.uiuc.edu/bee/honeybee_project.htm.]