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.

Genetic and epigenetic variations contributed by Alu retrotransposition.

BACKGROUND: De novo retrotransposition of Alu elements has been recognized as a major driver for insertion polymorphisms in human populations. In this study, we exploited Alu-anchored bisulfite PCR libraries to identify evolutionarily recent Alu element insertions, and to investigate their genetic and epigenetic variation. RESULTS: A total of 327 putatively recent Alu insertions were identified, altogether represented by 1,762 sequence reads. Nearly all such de novo retrotransposition events (316/327) were novel. Forty-seven out of forty-nine randomly selected events, corresponding to nineteen genomic loci, were sequence-verified. Alu element insertions remained hemizygous in one or more individuals in sixteen of the nineteen genomic loci. The Alu elements were found to be enriched for young Alu families with characteristic sequence features, such as the presence of a longer poly(A) tail. In addition, we documented the occurrence of a duplication of the AT-rich target site in their immediate flanking sequences, a hallmark of retrotransposition. Furthermore, we found the sequence motif (TT/AAAA) that is recognized by the ORF2P protein encoded by LINE-1 in their 5'-flanking regions, consistent with the fact that Alu retrotransposition is facilitated by LINE-1 elements. While most of these Alu elements were heavily methylated, we identified an Alu localized 1.5 kb downstream of TOMM5 that exhibited a completely unmethylated left arm. Interestingly, we observed differential methylation of its immediate 5' and 3' flanking CpG dinucleotides, in concordance with the unmethylated and methylated statuses of its internal 5' and 3' sequences, respectively. Importantly, TOMM5's CpG island and the 3 Alu repeats and 1 MIR element localized upstream of this newly inserted Alu were also found to be unmethylated. Methylation analyses of two additional genomic loci revealed no methylation differences in CpG dinucleotides flanking the Alu insertion sites in the two homologous chromosomes, irrespective of the presence or absence of the insertion. CONCLUSIONS: We anticipate that the combination of methodologies utilized in this study, which included repeat-anchored bisulfite PCR sequencing and the computational analysis pipeline herein reported, will prove invaluable for the generation of genetic and epigenetic variation maps.

Cyanobacterial contribution to algal nuclear genomes is primarily limited to plastid functions.

A single cyanobacterial primary endosymbiosis that occurred approximately 1.5 billion years ago is believed to have given rise to the plastid in the common ancestor of the Plantae or Archaeplastida--the eukaryotic supergroup comprising red, green (including land plants), and glaucophyte algae. Critical to plastid establishment was the transfer of endosymbiont genes to the host nucleus (i.e., endosymbiotic gene transfer [EGT]). It has been postulated that plastid-derived EGT played a significant role in plant nuclear-genome evolution, with 18% (or 4,500) of all nuclear genes in Arabidopsis thaliana having a cyanobacterial origin with about one-half of these recruited for nonplastid functions. Here, we determine whether the level of cyanobacterial gene recruitment proposed for Arabidopsis is of the same magnitude in the algal sisters of plants by analyzing expressed-sequence tag (EST) data from the glaucophyte alga Cyanophora paradoxa. Bioinformatic analysis of 3,576 Cyanophora nuclear genes shows that 10.8% of these with significant database hits are of cyanobacterial origin and one-ninth of these have nonplastid functions. Our data indicate that unlike plants, early-diverging algal groups appear to retain a smaller number of endosymbiont genes in their nucleus, with only a minor proportion of these recruited for nonplastid functions.

Gene discovery in an invasive tephritid model pest species, the Mediterranean fruit fly, Ceratitis capitata.

BACKGROUND: The medfly, Ceratitis capitata, is a highly invasive agricultural pest that has become a model insect for the development of biological control programs. Despite research into the behavior and classical and population genetics of this organism, the quantity of sequence data available is limited. We have utilized an expressed sequence tag (EST) approach to obtain detailed information on transcriptome signatures that relate to a variety of physiological systems in the medfly; this information emphasizes on reproduction, sex determination, and chemosensory perception, since the study was based on normalized cDNA libraries from embryos and adult heads. RESULTS: A total of 21,253 high-quality ESTs were obtained from the embryo and head libraries. Clustering analyses performed separately for each library resulted in 5201 embryo and 6684 head transcripts. Considering an estimated 19% overlap in the transcriptomes of the two libraries, they represent about 9614 unique transcripts involved in a wide range of biological processes and molecular functions. Of particular interest are the sequences that share homology with Drosophila genes involved in sex determination, olfaction, and reproductive behavior. The medfly transformer2 (tra2) homolog was identified among the embryonic sequences, and its genomic organization and expression were characterized. CONCLUSION: The sequences obtained in this study represent the first major dataset of expressed genes in a tephritid species of agricultural importance. This resource provides essential information to support the investigation of numerous questions regarding the biology of the medfly and other related species and also constitutes an invaluable tool for the annotation of complete genome sequences. Our study has revealed intriguing findings regarding the transcript regulation of tra2 and other sex determination genes, as well as insights into the comparative genomics of genes implicated in chemosensory reception and reproduction.

A pig multi-tissue normalised cDNA library: large-scale sequencing, cluster analysis and 9K micro-array resource generation.

BACKGROUND: Domestic animal breeding and product quality improvement require the control of reproduction, nutrition, health and welfare in these animals. It is thus necessary to improve our knowledge of the major physiological functions and their interactions. This would be greatly enhanced by the availability of expressed gene sequences in the databases and by cDNA arrays allowing the transcriptome analysis of any function.The objective within the AGENAE French program was to initiate a high-throughput cDNA sequencing program of a 38-tissue normalised library and generate a diverse microarray for transcriptome analysis in pig species. RESULTS: We constructed a multi-tissue cDNA library, which was normalised and subtracted to reduce the redundancy of the clones. Expressed Sequence Tags were produced and 24449 high-quality sequences were released in EMBL database. The assembly of all the public ESTs (available through SIGENAE website) resulted in 40786 contigs and 54653 singletons. At least one Agenae sequence is present in 11969 contigs (12.5%) and in 9291 of the deeper-than-one-contigs (22.8%). Sequence analysis showed that both normalisation and subtraction processes were successful and that the initial tissue complexity was maintained in the final libraries. A 9K nylon cDNA microarray was produced and is available through CRB-GADIE. It will allow high sensitivity transcriptome analyses in pigs. CONCLUSION: In the present work, a pig multi-tissue cDNA library was constructed and a 9K cDNA microarray designed. It contributes to the Expressed Sequence Tags pig data, and offers a valuable tool for transcriptome analysis.

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.

MicroRNA-10a Regulation of Proinflammatory Mediators: An Important Component of Untreated Juvenile Dermatomyositis.

OBJECTIVE: To identify differentially expressed microRNA (miRNA) in muscle biopsies (MBx) from 15 untreated children with juvenile dermatomyositis (JDM) compared with 5 controls. METHODS: Following MBx miRNA profiling, differentially expressed miRNA and their protein targets were validated by quantitative real-time PCR (qRT-PCR) and immunological assay. The association of miRNA-10a and miRNA-10b with clinical data was evaluated, including Disease Activity Score (DAS), von Willebrand factor antigen (vWF:Ag), nailfold capillary end row loops, duration of untreated disease, and tumor necrosis factor (TNF)-α-308A allele. RESULTS: In JDM, 16/362 miRNA were significantly differentially expressed [false discovery rate (FDR) < 0.05]. Among these, miRNA-10a was the most downregulated miRNA in both FDR and ranking of fold change: miRNA-10a = -2.27-fold, miRNA-10b = -1.80-fold. Decreased miRNA-10a and miRNA-10b expressions were confirmed using q RT-PCR: -4.16 and -2.59 fold, respectively. The qRT-PCR documented that decreased miRNA-10a expression was related to increased vascular cell adhesion molecule 1 in 13 of these JDM cases (correlation -0.67, p = 0.012), unlike miRNA-10b data (not significant). Concurrent JDM plasma contained increased levels of interleukin (IL) 6 (p = 0.0363), IL-8 (p = 0.0005), TNF-α (p = 0.0011), and monocyte chemoattractant proteins 1 (p = 0.0139). Decreased miRNA-10a, but not miRNA-10b, was associated with the TNF-α-308A allele (p = 0.015). In the 15 JDM, a trend of association of miRNA-10a (but not miRNA-10b) with vWF:Ag and DAS was observed. CONCLUSION: MiRNA-10a downregulation is an important element in untreated JDM muscle pathophysiology. We speculate that muscle miRNA expression in adult dermatomyositis differs from muscle miRNA expression in untreated childhood JDM.

Insights into a dinoflagellate genome through expressed sequence tag analysis.

BACKGROUND: Dinoflagellates are important marine primary producers and grazers and cause toxic "red tides". These taxa are characterized by many unique features such as immense genomes, the absence of nucleosomes, and photosynthetic organelles (plastids) that have been gained and lost multiple times. We generated EST sequences from non-normalized and normalized cDNA libraries from a culture of the toxic species Alexandrium tamarense to elucidate dinoflagellate evolution. Previous analyses of these data have clarified plastid origin and here we study the gene content, annotate the ESTs, and analyze the genes that are putatively involved in DNA packaging. RESULTS: Approximately 20% of the 6,723 unique (11,171 total 3'-reads) ESTs data could be annotated using Blast searches against GenBank. Several putative dinoflagellate-specific mRNAs were identified, including one novel plastid protein. Dinoflagellate genes, similar to other eukaryotes, have a high GC-content that is reflected in the amino acid codon usage. Highly represented transcripts include histone-like (HLP) and luciferin binding proteins and several genes occur in families that encode nearly identical proteins. We also identified rare transcripts encoding a predicted protein highly similar to histone H2A.X. We speculate this histone may be retained for its role in DNA double-strand break repair. CONCLUSION: This is the most extensive collection to date of ESTs from a toxic dinoflagellate. These data will be instrumental to future research to understand the unique and complex cell biology of these organisms and for potentially identifying the genes involved in toxin production.

Human rhinovirus infection causes different DNA methylation changes in nasal epithelial cells from healthy and asthmatic subjects.

BACKGROUND: Mechanisms underlying the development of virus-induced asthma exacerbations remain unclear. To investigate if epigenetic mechanisms could be involved in virus-induced asthma exacerbations, we undertook DNA methylation profiling in asthmatic and healthy nasal epithelial cells (NECs) during Human Rhinovirus (HRV) infection in vitro. METHODS: Global and loci-specific methylation profiles were determined via Alu element and Infinium Human Methylation 450 K microarray, respectively. Principal components analysis identified the genomic loci influenced the most by disease-status and infection. Real-time PCR and pyrosequencing were used to confirm gene expression and DNA methylation, respectively. RESULTS: HRV infection significantly increased global DNA methylation in cells from asthmatic subjects only (43.6% to 44.1%, p = 0.04). Microarray analysis revealed 389 differentially methylated loci either based on disease status, or caused by virus infection. There were disease-associated DNA methylation patterns that were not affected by HRV infection as well as HRV-induced DNA methylation changes that were unique to each group. A common methylation locus stood out in response to HRV infection in both groups, where the small nucleolar RNA, H/ACA box 12 (SNORA12) is located. Further analysis indicated that a relationship existed between SNORA12 DNA methylation and gene expression in response to HRV infection. CONCLUSIONS: We describe for the first time that Human rhinovirus infection causes DNA methylation changes in airway epithelial cells that differ between asthmatic and healthy subjects. These epigenetic differences may possibly explain the mechanism by which respiratory viruses cause asthma exacerbations.

ANALYSIS OF ALEXANDRIUM TAMARENSE (DINOPHYCEAE) GENES REVEALS THE COMPLEX EVOLUTIONARY HISTORY OF A MICROBIAL EUKARYOTE().

Microbial eukaryotes may extinguish much of their nuclear phylogenetic history due to endosymbiotic/horizontal gene transfer (E/HGT). We studied E/HGT in 32,110 contigs of expressed sequence tags (ESTs) from the dinoflagellate Alexandrium tamarense (Dinophyceae) using a conservative phylogenomic approach. The vast majority of predicted proteins (86.4%) in this alga are novel or dinoflagellate-specific. We searched for putative homologs of these predicted proteins against a taxonomically broadly sampled protein database that includes all currently available data from algae and protists and reconstructed a phylogeny from each of the putative homologous protein sets. Of the 2,523 resulting phylogenies, 14-17% are potentially impacted by E/HGT involving both prokaryote and eukaryote lineages, with 2-4% showing clear evidence of reticulate evolution. The complex evolutionary histories of the remaining proteins, many of which may also have been affected by E/HGT, cannot be interpreted using our approach with currently available gene data. We present empirical evidence of reticulate genome evolution that combined with inadequate or highly complex phylogenetic signal in many proteins may impede genome-wide approaches to infer the tree of microbial eukaryotes.

Analysis of expressed sequence tags from the four main developmental stages of Trypanosoma congolense.

Trypanosoma congolense is one of the most economically important pathogens of livestock in Africa. Culture-derived parasites of each of the three main insect stages of the T. congolense life cycle, i.e., the procyclic, epimastigote and metacyclic stages, and bloodstream stage parasites isolated from infected mice, were used to construct stage-specific cDNA libraries and expressed sequence tags (ESTs or cDNA clones) in each library were sequenced. Thirteen EST clusters encoding different variant surface glycoproteins (VSGs) were detected in the metacyclic library and 26 VSG EST clusters were found in the bloodstream library, 6 of which are shared by the metacyclic library. Rare VSG ESTs are present in the epimastigote library, and none were detected in the procyclic library. ESTs encoding enzymes that catalyze oxidative phosphorylation and amino acid metabolism are about twice as abundant in the procyclic and epimastigote stages as in the metacyclic and bloodstream stages. In contrast, ESTs encoding enzymes involved in glycolysis, the citric acid cycle and nucleotide metabolism are about the same in all four developmental stages. Cysteine proteases, kinases and phosphatases are the most abundant enzyme groups represented by the ESTs. All four libraries contain T. congolense-specific expressed sequences not present in the Trypanosoma brucei and Trypanosoma cruzi genomes. Normalized cDNA libraries were constructed from the metacyclic and bloodstream stages, and found to be further enriched for T. congolense-specific ESTs. Given that cultured T. congolense offers an experimental advantage over other African trypanosome species, these ESTs provide a basis for further investigation of the molecular properties of these four developmental stages, especially the epimastigote and metacyclic stages for which it is difficult to obtain large quantities of organisms. The T. congolense EST databases are available at: http://www.sanger.ac.uk/Projects/T_congolense/EST_index.shtml. The sequence data have been submitted to EMBL under the following accession numbers: FN263376-FN292969.

Global demethylation of rat chondrosarcoma cells after treatment with 5-aza-2'-deoxycytidine results in increased tumorigenicity.

Abnormal patterns of DNA methylation are observed in several types of human cancer. While localized DNA methylation of CpG islands has been associated with gene silencing, the effect that genome-wide loss of methylation has on tumorigenesis is not completely known. To examine its effect on tumorigenesis, we induced DNA demethylation in a rat model of human chondrosarcoma using 5-aza-2-deoxycytidine. Rat specific pyrosequencing assays were utilized to assess the methylation levels in both LINEs and satellite DNA sequences following 5-aza-2-deoxycytidine treatment. Loss of DNA methylation was accompanied by an increase in invasiveness of the rat chondrosarcoma cells, in vitro, as well as by an increase in tumor growth in vivo. Subsequent microarray analysis provided insight into the gene expression changes that result from 5-aza-2-deoxycytidine induced DNA demethylation. In particular, two genes that may function in tumorigenesis, sox-2 and midkine, were expressed at low levels in control cells but upon 5-aza-2-deoxycytidine treatment these genes became overexpressed. Promoter region DNA analysis revealed that these genes were methylated in control cells but became demethylated following 5-aza-2-deoxycytidine treatment. Following withdrawal of 5-aza-2-deoxycytidine, the rat chondrosarcoma cells reestablished global DNA methylation levels that were comparable to that of control cells. Concurrently, invasiveness of the rat chondrosarcoma cells, in vitro, decreased to a level indistinguishable to that of control cells. Taken together these experiments demonstrate that global DNA hypomethylation induced by 5-aza-2-deoxycytidine may promote specific aspects of tumorigenesis in rat chondrosarcoma cells.

Genome sequence of Aedes aegypti, a major arbovirus vector.

We present a draft sequence of the genome of Aedes aegypti, the primary vector for yellow fever and dengue fever, which at approximately 1376 million base pairs is about 5 times the size of the genome of the malaria vector Anopheles gambiae. Nearly 50% of the Ae. aegypti genome consists of transposable elements. These contribute to a factor of approximately 4 to 6 increase in average gene length and in sizes of intergenic regions relative to An. gambiae and Drosophila melanogaster. Nonetheless, chromosomal synteny is generally maintained among all three insects, although conservation of orthologous gene order is higher (by a factor of approximately 2) between the mosquito species than between either of them and the fruit fly. An increase in genes encoding odorant binding, cytochrome P450, and cuticle domains relative to An. gambiae suggests that members of these protein families underpin some of the biological differences between the two mosquito species.

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.

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

ESTprep: preprocessing cDNA sequence reads.

MOTIVATION: High accuracy of data always governs the large-scale gene discovery projects. The data should not only be trustworthy but should be correctly annotated for various features it contains. Sequence errors are inherent in single-pass sequences such as ESTs obtained from automated sequencing. These errors further complicate the automated identification of EST-related sequencing. A tool is required to prepare the data prior to advanced annotation processing and submission to public databases. RESULTS: This paper describes ESTprep, a program designed to preprocess expressed sequence tag (EST) sequences. It identifies the location of features present in ESTs and allows the sequence to pass only if it meets various quality criteria. Use of ESTprep has resulted in substantial improvement in accurate EST feature identification and fidelity of results submitted to GenBank. AVAILABILITY: The program is freely available for download from http://genome.uiowa.edu/pubsoft/software.html

1274 full-open reading frames of transcripts expressed in the developing mouse nervous system.

As part of the trans-National Institutes of Health (NIH) Mouse Brain Molecular Anatomy Project (BMAP), and in close coordination with the NIH Mammalian Gene Collection Program (MGC), we initiated a large-scale project to clone, identify, and sequence the complete open reading frame (ORF) of transcripts expressed in the developing mouse nervous system. Here we report the analysis of the ORF sequence of 1274 cDNAs, obtained from 47 full-length-enriched cDNA libraries, constructed by using a novel approach, herein described. cDNA libraries were derived from size-fractionated cytoplasmic mRNA isolated from brain and eye tissues obtained at several embryonic stages and postnatal days. Altogether, including the full-ORF MGC sequences derived from these libraries by the MGC sequencing team, NIH_BMAP full-ORF sequences correspond to approximately 20% of all transcripts currently represented in mouse MGC. We show that NIH_BMAP clones comprise 68% of mouse MGC cDNAs > or =5 kb, and 54% of those > or =4 kb, as of March 15, 2004. Importantly, we identified transcripts, among the 1274 full-ORF sequences, that are exclusively or predominantly expressed in brain and eye tissues, many of which encode yet uncharacterized proteins.

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