De Novo Assembly And Transcriptome Analysis Of The Rubber Tree (hevea Brasiliensis) And Snp Markers Development For Rubber Biosynthesis Pathways.

Hevea brasiliensis (Willd. Ex Adr. Juss.) Muell.-Arg. is the primary source of natural rubber that is native to the Amazon rainforest. The singular properties of natural rubber make it superior to and competitive with synthetic rubber for use in several applications. Here, we performed RNA sequencing (RNA-seq) of H. brasiliensis bark on the Illumina GAIIx platform, which generated 179,326,804 raw reads on the Illumina GAIIx platform. A total of 50,384 contigs that were over 400 bp in size were obtained and subjected to further analyses. A similarity search against the non-redundant (nr) protein database returned 32,018 (63%) positive BLASTx hits. The transcriptome analysis was annotated using the clusters of orthologous groups (COG), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Pfam databases. A search for putative molecular marker was performed to identify simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs). In total, 17,927 SSRs and 404,114 SNPs were detected. Finally, we selected sequences that were identified as belonging to the mevalonate (MVA) and 2-C-methyl-D-erythritol 4-phosphate (MEP) pathways, which are involved in rubber biosynthesis, to validate the SNP markers. A total of 78 SNPs were validated in 36 genotypes of H. brasiliensis. This new dataset represents a powerful information source for rubber tree bark genes and will be an important tool for the development of microsatellites and SNP markers for use in future genetic analyses such as genetic linkage mapping, quantitative trait loci identification, investigations of linkage disequilibrium and marker-assisted selection.

MicroRNA transcriptome profiles during swine skeletal muscle development

Background: MicroRNA (miR) are a class of small RNAs that regulate gene expression by inhibiting translation of protein encoding transcripts. To evaluate the role of miR in skeletal muscle of swine, global microRNA abundance was measured at specific developmental stages including proliferating satellite cells, three stages of fetal growth, day-old neonate, and the adult. Results: Twelve potential novel miR were detected that did not match previously reported sequences. In addition, a number of miR previously reported to be expressed in mammalian muscle were detected, having a variety of abundance patterns through muscle development. Muscle-specific miR-206 was nearly absent in proliferating satellite cells in culture, but was the highest abundant miR at other time points evaluated. In addition, miR-1 was moderately abundant throughout developmental stages with highest abundance in the adult. In contrast, miR-133 was moderately abundant in adult muscle and either not detectable or lowly abundant throughout fetal and neonate development. Changes in abundance of ubiquitously expressed miR were also observed. MiR-432 abundance was highest at the earliest stage of fetal development tested (60 day-old fetus) and decreased throughout development to the adult. Conversely, miR-24 and miR-27 exhibited greatest abundance in proliferating satellite cells and the adult, while abundance of miR-368, miR-376, and miR-423-5p was greatest in the neonate. Conclusion: These data present a complete set of transcriptome profiles to evaluate miR abundance at specific stages of skeletal muscle growth in swine. Identification of these miR provides an initial group of miR that may play a vital role in muscle development and growth.

Poly (A)(+) Transcriptome Assessment of ERBB2-Induced Alterations in Breast Cell Lines

We report the first quantitative and qualitative analysis of the poly (A)(+) transcriptome of two human mammary cell lines, differentially expressing (human epidermal growth factor receptor) an oncogene over-expressed in approximately 25% of human breast tumors. Full-length cDNA populations from the two cell lines were digested enzymatically, individually tagged according to a customized method for library construction, and simultaneously sequenced by the use of the Titanium 454-Roche-platform. Comprehensive bioinformatics analysis followed by experimental validation confirmed novel genes, splicing variants, single nucleotide polymorphisms, and gene fusions indicated by RNA-seq data from both samples. Moreover, comparative analysis showed enrichment in alternative events, especially in the exon usage category, in ERBB2 over-expressing cells, data indicating regulation of alternative splicing mediated by the oncogene. Alterations in expression levels of genes, such as LOX, ATP5L, GALNT3, and MME revealed by large-scale sequencing were confirmed between cell lines as well as in tumor specimens with different ERBB2 backgrounds. This approach was shown to be suitable for structural, quantitative, and qualitative assessment of complex transcriptomes and revealed new events mediated by ERBB2 overexpression, in addition to potential molecular targets for breast cancer that are driven by this oncogene.

Crotalus durissus collilineatus venom gland transcriptome: Analysis of gene expression profile

Crotalus durissus rattlesnakes are responsible for the most lethal cases of snakebites in Brazil. Crotalus durissus collilineatus subspecies is related to a great number of accidents in Southeast and Central West regions, but few studies on its venom composition have been carried out to date. In an attempt to describe the transcriptional profile of the C. durissus collilineatus venom gland, we generated a cDNA library and the sequences obtained could be identified by similarity searches on existing databases. Out of 673 expressed sequence tags (ESTs) 489 produced readable sequences comprising 201 singletons and 47 clusters of two or more ESTs. One hundred and fifty reads (60.5%) produced significant hits to known sequences. The results showed a predominance of toxin-coding ESTs instead of transcripts coding for proteins involved in all cellular functions. The most frequent toxin was crotoxin, comprising 88% of toxin-coding sequences. Crotoxin B, a basic phospholipase A(2) (PLA(2)) subunit of crotoxin, was represented in more variable forms comparing to the non-enzymatic subunit (crotoxin A), and most sequences coding this molecule were identified as CB1 isoform from Crotalus durissus terrificus venom. Four percent of toxin-related sequences in this study were identified as growth factors, comprising five sequences for vascular endothelial growth factor (VEGF) and one for nerve growth factor (NGF) that showed 100% of identity with C. durissus terrificus NGF. We also identified two clusters for metalloprotease from PII class comprising 3% of the toxins, and two for serine proteases, including gyroxin (2.5%). The remaining 2.5% of toxin-coding ESTs represent singletons identified as homologue sequences to cardiotoxin, convulxin, angiotensin-converting enzyme inhibitor and C-type natriuretic peptide, Ohanin, crotamin and PLA(2) inhibitor. These results allowed the identification of the most common classes of toxins in C. durissus collilineatus snake venom, also showing some unknown classes for this subspecies and even for C. durissus species, such as cardiotoxins and VEGF. (C) 2009 Published by Elsevier Masson SAS.

Transcriptome changes induced by docetaxel in human mammary cell lines expressing different levels of ERBB2

The taxane docetaxel is currently the most effective chemotherapeutic drug for the treatment of advanced breast cancer. However, a considerable proportion of breast cancer patients do not respond positively to docetaxel. The mechanisms of docetaxel resistance are poorly understood. Overexpression of ERBB2 occurs in 15-30% of breast tumors and is associated with chemoresistance to a variety of anticancer drugs. In the present study, we sought to identify genes involved in ERBB2-mediated chemoresistance to docetaxel. We generated SAGE libraries from two human mammary cell lines expressing basal (HB4a) and high (C5.2) levels of ERBB2 before and after intensive exposure to docetaxel and identified potential ERBB2 target genes implicated in a variety of cellular processes including cell proliferation, cell adhesion, apoptosis and cytoskeleton organization. Comparison of the transcriptome of the cell lines before and after docetaxel exposure revealed substantially different expression patterns. Twenty-one differentially expressed genes between HB4a and C5.2 cell lines, before and after docetaxel treatment, were further analyzed by qPCR. The alterations in the expression patterns in HB4a and C5.2 cell lines in response to docetaxel treatment observed by SAGE analysis were confirmed by qPCR for the majority of the genes analyzed. Our study provides a comprehensive view of the expression changes induced in two human mammary cells expressing different levels of ERBB2 in response to docetaxel that could contribute to the elucidation of the mechanisms involved in ERBB2-mediated chemoresistance in breast cancer.

Analysis of the Nicotiana tabacum Stigma/Style Transcriptome Reveals Gene Expression Differences between Wet and Dry Stigma Species

The success of plant reproduction depends on pollen-pistil interactions occurring at the stigma/style. These interactions vary depending on the stigma type: wet or dry. Tobacco (Nicotiana tabacum) represents a model of wet stigma, and its stigmas/styles express genes to accomplish the appropriate functions. For a large-scale study of gene expression during tobacco pistil development and preparation for pollination, we generated 11,216 high-quality expressed sequence tags (ESTs) from stigmas/styles and created the TOBEST database. These ESTs were assembled in 6,177 clusters, from which 52.1% are pistil transcripts/genes of unknown function. The 21 clusters with the highest number of ESTs (putative higher expression levels) correspond to genes associated with defense mechanisms or pollen-pistil interactions. The database analysis unraveled tobacco sequences homologous to the Arabidopsis (Arabidopsis thaliana) genes involved in specifying pistil identity or determining normal pistil morphology and function. Additionally, 782 independent clusters were examined by macroarray, revealing 46 stigma/style preferentially expressed genes. Real-time reverse transcription-polymerase chain reaction experiments validated the pistil-preferential expression for nine out of 10 genes tested. A search for these 46 genes in the Arabidopsis pistil data sets demonstrated that only 11 sequences, with putative equivalent molecular functions, are expressed in this dry stigma species. The reverse search for the Arabidopsis pistil genes in the TOBEST exposed a partial overlap between these dry and wet stigma transcriptomes. The TOBEST represents the most extensive survey of gene expression in the stigmas/styles of wet stigma plants, and our results indicate that wet and dry stigmas/styles express common as well as distinct genes in preparation for the pollination process.

Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs

Only a small proportion of the mouse genome is transcribed into mature messenger RNA transcripts. There is an international collaborative effort to identify all full-length mRNA transcripts from the mouse, and to ensure that each is represented in a physical collection of clones. Here we report the manual annotation of 60,770 full-length mouse complementary DNA sequences. These are clustered into 33,409 'transcriptional units', contributing 90.1% of a newly established mouse transcriptome database. Of these transcriptional units, 4,258 are new protein-coding and 11,665 are new non-coding messages, indicating that non-coding RNA is a major component of the transcriptome. 41% of all transcriptional units showed evidence of alternative splicing. In protein-coding transcripts, 79% of splice variations altered the protein product. Whole-transcriptome analyses resulted in the identification of 2,431 sense-antisense pairs. The present work, completely supported by physical clones, provides the most comprehensive survey of a mammalian transcriptome so far, and is a valuable resource for functional genomics.

Analysis of the mouse transcriptome for genes involved in the function of the nervous system

We analyzed the mouse Representative Transcript and Protein Set for molecules involved in brain function. We found full-length cDNAs of many known brain genes and discovered new members of known brain gene families, including Family 3 G-protein coupled receptors, voltage-gated channels, and connexins. We also identified previously unknown candidates for secreted neuroactive molecules. The existence of a large number of unique brain ESTs suggests an additional molecular complexity that remains to be explored. A list of genes containing CAG stretches in the coding region represents a first step in the potential identification of candidates for hereditary neurological disorders.

The choroid plexus transcriptome reveals changes in type I and II interferon responses in a mouse model of Alzheimer's disease

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by a marked decline in cognition and memory function. Increasing evidence highlights the essential role of neuroinflammatory and immune-related molecules, including those produced at the brain barriers, on brain immune surveillance, cellular dysfunction and amyloid beta (Aß) pathology in AD. Therefore, understanding the response at the brain barriers may unravel novel pathways of relevance for the pathophysiology of AD. Herein, we focused on the study of the choroid plexus (CP), which constitutes the blood-cerebrospinal fluid barrier, in aging and in AD. Specifically, we used the PDGFB-APPSwInd (J20) transgenic mouse model of AD, which presents early memory decline and progressive Aß accumulation, and littermate age-matched wild-type (WT) mice, to characterize the CP transcriptome at 3, 5-6 and 11-12months of age. The most striking observation was that the CP of J20 mice displayed an overall overexpression of type I interferon (IFN) response genes at all ages. Moreover, J20 mice presented a high expression of type II IFN genes in the CP at 3months, which became lower than WT at 5-6 and 11-12months. Importantly, along with a marked memory impairment and increased glial activation, J20 mice also presented a similar overexpression of type I IFN genes in the dorsal hippocampus at 3months. Altogether, these findings provide new insights on a possible interplay between type I and II IFN responses in AD and point to IFNs as targets for modulation in cognitive decline.

Global transcriptome sequencing identifies chlamydospore specific markers in Candida albicans and Candida dubliniensis.

Candida albicans and Candida dubliniensis are pathogenic fungi that are highly related but differ in virulence and in some phenotypic traits. During in vitro growth on certain nutrient-poor media, C. albicans and C. dubliniensis are the only yeast species which are able to produce chlamydospores, large thick-walled cells of unknown function. Interestingly, only C. dubliniensis forms pseudohyphae with abundant chlamydospores when grown on Staib medium, while C. albicans grows exclusively as a budding yeast. In order to further our understanding of chlamydospore development and assembly, we compared the global transcriptional profile of both species during growth in liquid Staib medium by RNA sequencing. We also included a C. albicans mutant in our study which lacks the morphogenetic transcriptional repressor Nrg1. This strain, which is characterized by its constitutive pseudohyphal growth, specifically produces masses of chlamydospores in Staib medium, similar to C. dubliniensis. This comparative approach identified a set of putatively chlamydospore-related genes. Two of the homologous C. albicans and C. dubliniensis genes (CSP1 and CSP2) which were most strongly upregulated during chlamydospore development were analysed in more detail. By use of the green fluorescent protein as a reporter, the encoded putative cell wall related proteins were found to exclusively localize to C. albicans and C. dubliniensis chlamydospores. Our findings uncover the first chlamydospore specific markers in Candida species and provide novel insights in the complex morphogenetic development of these important fungal pathogens.

Physiology and transcriptome of the polycyclic aromatic hydrocarbon-degrading Sphingomonas sp. LH128 after long-term starvation.

The survival, physiology and gene expression profile of the phenanthrene-degrading Sphingomonas sp. LH128 was examined after an extended period of complete nutrient starvation and compared with a non-starved population that had been harvested in exponential phase. After 6 months of starvation in an isotonic solution, only 5 % of the initial population formed culturable cells. Microscopic observation of GFP fluorescent cells, however, suggested that a larger fraction of cells (up to 80 %) were still alive and apparently had entered a viable but non-culturable (VBNC) state. The strain displayed several cellular and genetic adaptive strategies to survive long-term starvation. Flow cytometry, microscopic observation and fatty acid methyl ester (FAME) analysis showed a reduction in cell size, a change in cell shape and an increase in the degree of membrane fatty acid saturation. Transcriptome analysis showed decreased expression of genes involved in ribosomal protein biosynthesis, chromosomal replication, cell division and aromatic catabolism, increased expression of genes involved in regulation of gene expression and efflux systems, genetic translocations, and degradation of rRNA and fatty acids. Those phenotypic and transcriptomic changes were not observed after 4 h of starvation. Despite the starvation situation, the polycyclic aromatic hydrocarbon (PAH) catabolic activity was immediate upon exposure to phenanthrene. We conclude that a large fraction of cells maintain viability after an extended period of starvation apparently due to tuning the expression of a wide variety of cellular processes. Due to these survival attributes, bacteria of the genus Sphingomonas, like strain LH128, could be considered as suitable targets for use in remediation of nutrient-poor PAH-contaminated environments.

Insights into the transcriptome of oenocytes from Aedes aegypti pupae

Oenocytes are ectodermic cells present in the fat body of several insect species and these cells are considered to be analogous to the mammalian liver, based on their role in lipid storage, metabolism and secretion. Although oenocytes were identified over a century ago, little is known about their messenger RNA expression profiles. In this study, we investigated the transcriptome of Aedes aegypti oenocytes. We constructed a cDNA library from Ae. aegypti MOYO-R strain oenocytes collected from pupae and randomly sequenced 687 clones. After sequences editing and assembly, 326 high-quality contigs were generated. The most abundant transcripts identified corresponded to the cytochrome P450 superfamily, whose members have roles primarily related to detoxification and lipid metabolism. In addition, we identified 18 other transcripts with putative functions associated with lipid metabolism. One such transcript, a fatty acid synthase, is highly represented in the cDNA library of oenocytes. Moreover, oenocytes expressed several immunity-related genes and the majority of these genes were lysozymes. The transcriptional profile suggests that oenocytes play diverse roles, such as detoxification and lipid metabolism, and increase our understanding of the importance of oenocytes in Ae. aegypti homeostasis and immune competence.

Transcriptome analysis of the mobile genome ICEclc in Pseudomonas knackmussii B13.

Background: Integrative and conjugative elements (ICE) form a diverse group of DNA elements that are integrated in the chromosome of the bacterial host, but can occasionally excise and horizontally transfer to a new host cell. ICE come in different families, typically with a conserved core for functions controlling the element's behavior and a variable region providing auxiliary functions to the host. The ICEclc element of Pseudomonas knackmussii strain B13 is representative for a large family of chromosomal islands detected by genome sequencing approaches. It provides the host with the capacity to degrade chloroaromatics and 2-aminophenol. Results: Here we study the transcriptional organization of the ICEclc core region. By northern hybridizations, reverse-transcriptase polymerase chain reaction (RT-PCR) and Rapid Amplification of cDNA Ends (5'-RACE) fifteen transcripts were mapped in the core region. The occurrence and location of those transcripts were further confirmed by hybridizing labeled cDNA to a semi-tiling micro-array probing both strands of the ICEclc core region. Dot blot and semi-tiling array hybridizations demonstrated most of the core transcripts to be upregulated during stationary phase on 3-chlorobenzoate, but not on succinate or glucose. Conclusions: The transcription analysis of the ICEclc core region provides detailed insights in the mode of regulatory organization and will help to further understand the complex mode of behavior of this class of mobile elements. We conclude that ICEclc core transcription is concerted at a global level, more reminiscent of a phage program than of plasmid conjugation.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill.

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.