Full-length Expressed Sequence Tags Research Articles

Development and Evaluation of SSR Markers Based on Large Scale Full-Length Transcriptome Sequencing in Sugarcane

Sugarcane generates the vast majority of sugar consumed worldwide. The large-scale development of EST-SSR markers are crucial for the analysis of genetic mechanisms important in agronomic traits of sugarcane. In this study, based on published 107,598 full-length expressed sequence tags (ESTs) from 22 sugarcane cultivars, we identified a total of 46,043 loci of simple repeat sequences. The SSRs with di-nucleotide (8175 loci) and tri-nucleotide (10,709 loci) repeat motifs account for 40.96% (18,884) of the total SSR loci, with (AT/TA)n and (ACA/TGT)n being found to be the major repeat motifs. A total of 1070 pairs of primers were randomly selected to verify using 12 accessions including S. officinarum, S. robustum, S. spontaneum and cultivars. The results showed that 847 (79.1%) were successfully amplified, of which 349 (32.60%) were polymorphic in the analyzed samples. The polymorphic rates of EST-SSR markers with di-nucleotide, tri-nucleotide, tetra-nucleotide, penta-nucleotide and hexa-nucleotide repeat motifs were 36.82% (81), 49.40% (165), 43.04% (34), 30.19% (32) and 34.26% (37), respectively. Among them, the EST-SSRs comprising five types of repeat motifs (CA)n, (CCG)n, (CCGC)n, (CGGCG)n and (GCCGGC)n presented high degrees of polymorphism. Finally, 20 EST-SSR markers with polymorphisms were selected to perform a cluster analysis of 48 accessions from Saccharum and Erianthus genera. These results demonstrated that Chinese sugarcane cultivars from group IV (25/39) had a lower genetic divergence than cultivars from other countries and EST-SSR markers can also be used in determining the genetic diversity of Erianthus genus. This set of EST-SSR markers is an important molecular resource for analyzing genetic diversity, molecular fingerprinting and genetic mapping.

Generation and characterization of expressed sequence tags in young roots of tea (Camellia assamica)

Assam tea (Camellia assamica) is perennial crop susceptible to moisture stress. We used its tender roots to construct a cDNA library for the identification, functional annotation, and analysis of transcripts. A total of 811 full-length expressed sequence tags were generated. After processing and assembly, 207 unigenes comprising 58 contigs and 149 singletons were registered. Finally, 35.75 % of the unigenes could be assigned to functional categories based on the Arabidopsis proteome. There was 43 % of a coding GC content and 1 272 coding DNA sequences found in the unigenes. Codon usage analysis shows leucine as the highest (9.92 %) and tryptophan (2.0 %) as the lowest coded amino acids. Further, a comparative study with drought-induced genes of young roots (reported earlier) reveals that 4.83 % of genes required for normal growth of roots were also induced by a drought stress. Expressions of 10 unigenes under different abiotic stresses, such as drought, cold, and salinity, were further confirmed by RT-qPCR. The sequence tags generated in this study will be valuable resources for functional genomics study of tea and other woody crop plants in future.

Molecular cloning and characterization of Siamese crocodile (Crocodylus siamensis) copper, zinc superoxide dismutase (CSI-Cu,Zn-SOD) gene

Superoxide dismutase (SOD, EC 1.15.1.1) is an antioxidant enzyme found in all living cells. It regulates oxidative stress by breaking down superoxide radicals to oxygen and hydrogen peroxide. A gene coding for Cu,Zn-SOD was cloned and characterized from Siamese crocodile (Crocodylus siamensis; CSI). The full-length expressed sequence tag (EST) of this Cu,Zn-SOD gene (designated as CSI-Cu,Zn-SOD) contained 462bp encoding a protein of 154 amino acids without signal peptides, indicated as intracellular CSI-Cu,Zn-SOD. This agreed with the results from the phylogenetic tree, which indicated that CSI-Cu,Zn-SOD belonged to the intracellular Cu,Zn-SOD. Chromosomal location determined that the CSI-Cu,Zn-SOD was localized to the proximal region of the Siamese crocodile chromosome 1p. Several highly conserved motifs, two conserved signature sequences (GFHVHEFGDNT and GNAGGRLACGVI), and conserved amino acid residues for binding copper and zinc (His47, His49, His64, His72, His81, Asp84, and His120) were also identified in CSI-Cu,Zn-SOD. Real-time PCR analysis showed that CSI-Cu,Zn-SOD mRNA was expressed in all the tissues examined (liver, pancreas, lung, kidney, heart, and whole blood), which suggests a constitutively expressed gene in these tissues. Expression of the gene in Escherichia coli cells followed by purification yielded a recombinant CSI-Cu,Zn-SOD, with Km and Vmax values of 6.075mM xanthine and 1.4×10−3mmolmin−1mg−1, respectively. This Vmax value was 40 times lower than native Cu,Zn-SOD (56×10−3mmolmin−1mg−1), extracted from crocodile erythrocytes. This suggests that cofactors, protein folding properties, or post-translational modifications were lost during the protein purification process, leading to a reduction in the rate of enzyme activity in bacterial expression of CSI-Cu,Zn-SOD.

The Babesia bovis gene and promoter model: an update from full-length EST analysis.

BackgroundBabesia bovis is an apicomplexan parasite that causes babesiosis in infected cattle. Genomes of pathogens contain promising information that can facilitate the development of methods for controlling infections. Although the genome of B. bovis is publically available, annotated gene models are not highly reliable prior to experimental validation. Therefore, we validated a preproposed gene model of B. bovis and extended the associated annotations on the basis of experimentally obtained full-length expressed sequence tags (ESTs).ResultsFrom in vitro cultured merozoites, 12,286 clones harboring full-length cDNAs were sequenced from both ends using the Sanger method, and 6,787 full-length cDNAs were assembled. These were then clustered, and a nonredundant referential data set of 2,115 full-length cDNA sequences was constructed. The comparison of the preproposed gene model with our data set identified 310 identical genes, 342 almost identical genes, 1,054 genes with potential structural inconsistencies, and 409 novel genes. The median length of 5' untranslated regions (UTRs) was 152 nt. Subsequently, we identified 4,086 transcription start sites (TSSs) and 2,023 transcriptionally active regions (TARs) by examining 5' ESTs. We identified ATGGGG and CCCCAT sites as consensus motifs in TARs that were distributed around -50 bp from TSSs. In addition, we found ACACA, TGTGT, and TATAT sites, which were distributed periodically around TSSs in cycles of approximately 150 bp. Moreover, related periodical distributions were not observed in mammalian promoter regions.ConclusionsThe observations in this study indicate the utility of integrated bioinformatics and experimental data for improving genome annotations. In particular, full-length cDNAs with one-base resolution for TSSs enabled the identification of consensus motifs in promoter sequences and demonstrated clear distributions of identified motifs. These observations allowed the illustration of a model promoter composition, which supports the differences in transcriptional regulation frameworks between apicomplexan parasites and mammals.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-678) contains supplementary material, which is available to authorized users.

Cloning and characterising an unusual perforin from chicken (Gallus gallus)

In mammals the 67kDa pore-forming protein perforin is essential to the granule exocytosis pathway used by cytotoxic lymphocytes to eliminate virally infected and malignant cells. There is indirect evidence that this pathway exists in lower vertebrates such as teleost fish and birds, although in genome databases for the chicken and other bird species the perforin gene is incomplete and no full length expressed sequence tag has been reported. We present here the full gene and transcript sequence of chicken perforin. The inferred protein product contains an extended C-terminus that is at least 90 amino acids longer than any mammalian perforin, which is also evident in partial genomic sequences from other birds. To determine whether this extension is present in the translated protein, we raised two polyclonal antisera. The antisera identified a protein of just less than 80kDa in both transfected COS-1 cells and concanavalin A stimulated chicken splenocytes, indicating that the extended C-terminus is present in the mature protein. Our findings confirm that perforin exists in birds, and show that it is considerably longer than perforin of non-avian vertebrates.

Expression and Putative Function of Innate Immunity Genes under in situ Conditions in the Symbiotic Hydrothermal Vent Tubeworm Ridgeia piscesae

The relationships between hydrothermal vent tubeworms and sulfide-oxidizing bacteria have served as model associations for understanding chemoautotrophy and endosymbiosis. Numerous studies have focused on the physiological and biochemical adaptations that enable these symbioses to sustain some of the highest recorded carbon fixation rates ever measured. However, far fewer studies have explored the molecular mechanisms underlying the regulation of host and symbiont interactions, specifically those mediated by the innate immune system of the host. To that end, we conducted a series of studies where we maintained the tubeworm, Ridgeia piscesae, in high-pressure aquaria and examined global and quantitative changes in gene expression via high-throughput transcriptomics and quantitative real-time PCR (qPCR). We analyzed over 32,000 full-length expressed sequence tags as well as 26 Mb of transcript sequences from the trophosome (the organ that houses the endosymbiotic bacteria) and the plume (the gas exchange organ in contact with the free-living microbial community). R. piscesae maintained under conditions that promote chemoautotrophy expressed a number of putative cell signaling and innate immunity genes, including pattern recognition receptors (PRRs), often associated with recognizing microbe-associated molecular patterns (MAMPs). Eighteen genes involved with innate immunity, cell signaling, cell stress and metabolite exchange were further analyzed using qPCR. PRRs, including five peptidoglycan recognition proteins and a Toll-like receptor, were expressed significantly higher in the trophosome compared to the plume. Although PRRs are often associated with mediating host responses to infection by pathogens, the differences in expression between the plume and trophosome also implicate similar mechanisms of microbial recognition in interactions between the host and symbiont. We posit that regulation of this association involves a molecular “dialogue” between the partners that includes interactions between the host’s innate immune system and the symbiont.

An unusual plant triterpene synthase with predominant α‐amyrin‐producing activity identified by characterizing oxidosqualene cyclases from Malus × domestica

The pentacyclic triterpenes, in particular ursolic acid and oleanolic acid and their derivatives, exist abundantly in the plant kingdom, where they are well known for their anti-inflammatory, antitumour and antimicrobial properties. α-Amyrin and β-amyrin are the precursors of ursolic and oleanolic acids, respectively, formed by concerted cyclization of squalene epoxide by a complex synthase reaction. We identified three full-length expressed sequence tag sequences in cDNA libraries constructed from apple (Malus × domestica 'Royal Gala') that were likely to encode triterpene synthases. Two of these expressed sequence tag sequences were essentially identical (> 99% amino acid similarity; MdOSC1 and MdOSC3). MdOSC1 and MdOSC2 were expressed by transient expression in Nicotiana benthamiana leaves and by expression in the yeast Pichia methanolica. The resulting products were analysed by GC and GC-MS. MdOSC1 was shown to be a mixed amyrin synthase (a 5 : 1 ratio of α-amyrin to β-amyrin). MdOSC1 is the only triterpene synthase so far identified in which the level of α-amyrin produced is > 80% of the total product and is, therefore, primarily an α-amyrin synthase. No product was evident for MdOSC2 when expressed either transiently or in yeast, suggesting that this putative triterpene synthase is either encoded by a pseudogene or does not express well in these systems. Transcript expression analysis in Royal Gala indicated that the genes are mostly expressed in apple peel, and that the MdOSC2 expression level was much lower than that of MdOSC1 and MdOSC3 in all the tissues tested. Amyrin content analysis was undertaken by LC-MS, and demonstrated that levels and ratios differ between tissues, but that the true consequence of synthase activity is reflected in the ursolic/oleanolic acid content and in further triterpenoids derived from them. Phylogenetic analysis placed the three triterpene synthase sequences with other triterpene synthases that encoded either α-amyrin and/or β-amyrin synthase. MdOSC1 and MdOSC3 clustered with the multifunctional triterpene synthases, whereas MdOSC2 was most similar to the β-amyrin synthases.

Improved genome assembly and evidence-based global gene model set for the chordate Ciona intestinalis: new insight into intron and operon populations

An improved assembly of the Ciona intestinalis genome reveals that it contains non-canonical introns and that about 20% of Ciona genes reside in operons.

Identification and cloning of a novel tetraspanin (TSP) homologue fromBrugia malayi

This is the first report of a tetraspanin (TSP)-like molecule in the lymphatic filarial parasites. Expressed sequence tag (EST) database search for TSP like molecules in the filarial genome resulted in three significant EST hits (two partial ESTs from Brugia malayi and one full length EST from Wuchereria bancrofti). The full length gene cloned from B. malayi showed significant similarity to Caenorhabditis elegans TSP and human TSP and hence the gene was named B. malayi TSP (BmTSP). Subsequent Genbank analysis with the predicted ORF of BmTSP showed additional homologous genes reported from Schistosoma mansoni and Taenia solium parasites. Structural analyses showed that BmTSP has four transmembrane domains and other conserved domains such as CCG and two other critical cysteine residues present within the large extracellular loop similar to other reported TSPs. In addition, putative post-translational modifications such as N-glycosylation, protein kinase c phosphorylation, casein kinase II phosphorylation and N-myristoylation sites have been found in BmTSP sequence. Further, PCR analyses showed that BmTSP is differentially transcribed, with highest level of expression being present in the adult stages followed by L3 and mf stages. This study thus describes a novel TSP cloned from B. malayi, its putative functions in cuticle biogenesis and role in protective immunity.

Gene expression profiling of Duchenne muscular dystrophy skeletal muscle.

The primary cause of Duchenne muscular dystrophy (DMD) is a mutation in the dystrophin gene, leading to absence of the corresponding protein, disruption of the dystrophin-associated protein complex, and substantial changes in skeletal muscle pathology. Although the primary defect is known and the histological pathology well documented, the underlying molecular pathways remain in question. To clarify these pathways, we used expression microarrays to compare individual gene expression profiles for skeletal muscle biopsies from DMD patients and unaffected controls. We have previously published expression data for the 12,500 known genes and full-length expressed sequence tags (ESTs) on the Affymetrix HG-U95Av2 chips. Here we present comparative expression analysis of the 50,000 EST clusters represented on the remainder of the Affymetrix HG-U95 set. Individual expression profiles were generated for biopsies from 10 DMD patients and 10 unaffected control patients. Two methods of statistical analysis were used to interpret the resulting data (t-test analysis to determine the statistical significance of differential expression and geometric fold change analysis to determine the extent of differential expression). These analyses identified 183 probe sets (59 of which represent known genes) that differ significantly in expression level between unaffected and disease muscle. This study adds to our knowledge of the molecular pathways that are altered in the dystrophic state. In particular, it suggests that signaling pathways might be substantially involved in the disease process. It also highlights a large number of unknown genes whose expression is altered and whose identity therefore becomes important in understanding the pathogenesis of muscular dystrophy.

A New Member of Plant CS-lyases: A CYSTINE LYASE FROM ARABIDOPSIS THALIANA

Cystine lyases catalyze the breakdown of l-cystine to thiocysteine, pyruvate, and ammonia. Until now there are no reports of the identification of a plant cystine lyase at a molecular level, and it is not clear what biological role this class of enzymes have in plants. A cystine lyase was isolated from Brassica oleracea (L.), and partial amino acid sequencing allowed the corresponding full-length cDNA (BOCL3) to be cloned. The deduced amino acid sequence of BOCL3 showed highest homology to the deduced amino acid sequences of several Arabidopsis thaliana genes annotated as tyrosine aminotransferase-like, including a coronatine, jasmonic acid, and salt stress-inducible gene, CORI3 (78.8% identity), and the unidentified rooty/superroot1 gene (44.8% identity). A full-length expressed sequence tag clone of CORI3 was obtained and recombinant CORI3 was synthesized in Escherichia coli. Isolated recombinant CORI3 catalyzed a cystine lyase reaction, but no aminotransferase reactions. The present study identifies, for the first time, a cystine lyase from plants at a molecular level and redefines the functional assignment of the only functionally identified member of a group of A. thaliana genes annotated as tyrosine aminotransferase-like.