Identification of Genes Differentially Expressed Between Ochratoxin-Producing and Non-Producing Strains of Aspergillus westerdijkiae

Approximately 70 % of Aspergillus westerdijkiae strains are able to produce ochratoxin A (OTA), a nephrotoxic and carcinogenic mycotoxin which have been found in cereal and food commodities. Despite of its importance there is, up to now, no information available about which genes are differentially expressed between A. westerdijkiae ochratoxin-producing and non-producing strains. Using cDNA RDA approach we successfully sequenced 231 raw ESTs expected to be enriched in the ochratoxin-producing strain. BLASTX searches against the public databases showed that of these, 205 ESTs (79 %) exhibited significant similarities with proteins of known functions, 28 ESTs (11 %) had matches to hypothetical proteins, and the remaining 27 ESTs (10 %) had no significant hits. EST alignment resulted in a total of 14 non-redundant consensus sequences. Three putative genes encoding oxidoreductases were validated as up-expressed in the OTA producer strain using RT-qPCR approach. The expression of the putative genes encoding a cytochrome P450 family protein, 3-hydroxyphenylacetate-6-hydroxylase, and endoplasmic reticulum oxidoreductin were higher (32-, 2.8- and 20-fold respectively) in the OTA producer strain compared to the non-producer strain.     

Transcriptome generation and analysis from spleen of Indian catfish, Clarias batrachus (Linnaeus, 1758) through normalized cDNA library

Catfishes are commercially important fish for both the fisheries and aquaculture industry. Clarias batrachus, an Indian catfish species is economically important owing to its high demand. A normalized cDNA library was constructed from spleen of the Indian catfish to identify genes associated with immune function. One thousand nine hundred thirty seven ESTs were submitted to the GenBank with an average read length of approximately 700 bp. Clustering analysis of ESTs yielded 1,698 unique sequences, including 184 contigs and 1,514 singletons. Significant homology to known genes was found by homology searches against data in GenBank in 576 (34 %) ESTs, including similarity to functionally annotated unigenes for 158 ESTs. Additionally, 433 ESTs revealed similarity to unigenes and ESTs in the dbEST but the remaining 658 EST sequences (39 %) did not match any sequence in GenBank. Of a total of 1,698 ESTs generated, 65 ESTs were found to be associated with immune functions. Gene Ontology and KEGG pathway analyses of C. batrachus ESTs collectively revealed a preponderance of immune relevant pathways apart from the presence of pathways involved in protein processing, localization, folding and protein degradation. This study constitutes first EST analysis of lymphoid organ in aquaculturally important Indian catfish species and could pave the way for further research of immune-related genes and functional genomics in this catfish.     

Comparative gene expression in the symbiotic and aposymbiotic Aiptasia pulchella by expressed sequence tag analysis

Intracellular symbiotic relationships are prevalent between cnidarians, such as corals and sea anemones, and the photosynthetic dinoflagellate symbionts. However, there is little understanding about how the genes express when the symbiotic relationship is set up. To characterize genes involved in this association, the endosymbiosis between sea anemone, Aiptasia pulchella, and dinoflagellate zooxanthellae, Symbiodinium spp., was employed as a model. Two complementary DNA (cDNA) libraries were constructed from RNA isolated from symbiotic and aposymbiotic A. pulchella. Using single-pass sequencing of cDNA clones, a total of 870 expressed sequence tags (ESTs) clones were generated from the two libraries: 474 from symbiotic animal and 396 from aposymbiotic animal. The initial ESTs consisted of 143 clusters and 231 singletons. A BLASTX search revealed that 147 unique genes had similarities with protein sequences available from databases; 120 of these clones were categorized according to their putative function. However, many ESTs could not assign functionally. The putative roles of some of the identified genes relative to endosymbiosis were discussed. This is the first report of the use of EST analysis to examine the gene expression in symbiotic and aposymbiotic states of the cnidarians. The systematic analysis of EST from this study provides a useful database for future investigations of the molecular mechanisms involved in algal-cnidarian symbiosis.     

Generation of expressed sequence tags from low-CO2 and high-CO2 adapted cells of Chlamydomonas reinhardtii.

To characterize genes whose expression is induced in carbon-stress conditions, 12,969 and 13,450 5'-end expressed sequence tags (ESTs) were generated from cells grown in low-CO2 and high-CO2 conditions of the unicellular green alga, Chlamydomonas reinhardtii. These ESTs were clustered into 4436 and 3566 non-redundant EST groups, respectively. Comparison of their sequences with those of 3433 non-redundant ESTs previously generated from the cells under the standard growth condition indicated that 2665 and 1879 EST groups occurred only in the low-CO2 and high-CO2 populations, respectively. It was also noted that 96.2% and 96.0% of the cDNA species respectively obtained from the low-CO2 and high-CO2 conditions had no similar EST sequence deposited in the public databases. The EST species identified only in the low-CO2 treated cells included genes previously reported to be expressed specifically in low-CO2 acclimatized cells, suggesting that the ESTs generated in this study will be a useful source for analysis of genes related to carbon-stress acclimatization. The sequence information and search results of each clone will appear at the web site: http://www.kazusa.or.jp/en/plant/chlamy/EST/.     

Development of simple sequence repeat markers and diversity analysis in alfalfa (Medicago sativa L.)

Efficient and robust molecular markers are essential for molecular breeding in plant. Compared to dominant and bi-allelic markers, multiple alleles of simple sequence repeat (SSR) markers are particularly informative and superior in genetic linkage map and QTL mapping in autotetraploid species like alfalfa. The objective of this study was to enrich SSR markers directly from alfalfa expressed sequence tags (ESTs). A total of 12,371 alfalfa ESTs were retrieved from the National Center for Biotechnology Information. Total 774 SSR-containing ESTs were identified from 716 ESTs. On average, one SSR was found per 7.7 kb of EST sequences. Tri-nucleotide repeats (48.8 %) was the most abundant motif type, followed by di—(26.1 %), tetra—(11.5 %), penta—(9.7 %), and hexanucleotide (3.9 %). One hundred EST–SSR primer pairs were successfully designed and 29 exhibited polymorphism among 28 alfalfa accessions. The allele number per marker ranged from two to 21 with an average of 6.8. The PIC values ranged from 0.195 to 0.896 with an average of 0.608, indicating a high level of polymorphism of the EST–SSR markers. Based on the 29 EST–SSR markers, assessment of genetic diversity was conducted and found that Medicago sativa ssp. sativa was clearly different from the other subspecies. The high transferability of those EST–SSR markers was also found for relative species.     

Analysis of Jatropha curcas transcriptome for oil enhancement and genic markers

Oil-rich seeds of Jatropha curcas are being focussed as a source of bio-diesel. However, prior to its industrial use, a lot of crop improvement efforts are required in Jatropha. Availability of a large number of EST sequences of Jatropha in public domain allow identification of candidate genes for several agronomic characters including oil content in seeds. Here, we have analysed 42,477 ESTs of Jatropha spanning 22.9 Mbp for microsatellites and fatty acid metabolism related sequences. Unigene sequences were built using CAP 3 programme resulted in 12,358 contigs and 5,730 singlets. Nearly, 8 % unigenes showed presence of microsatellites, slightly over-represented compared to their occurrence in ESTs. Most of the microsatellites were either di- or tri-nucleotide repeats, while other categories of tetra-, penta- and hexa-nucleotide repeats together constituted ~4 % of total microsatellites. Assessment of functional relevance of unigenes was carried out using Blast2GO using its default settings. The overall sequence similarity level against sequences in ‘nr’ database was >80 %. A total of 931 sequences that participated in any of the pathways related to fatty acid or lipid metabolism were found at GO level 6. Among these, GO terms “Fatty acid metabolic process” and “Fatty acid biosynthetic process” were most over-represented. Overall, our work has due relevance in identifying molecular markers for the candidate genes for oil content in Jatropha seeds, and will prove to be an important reference for further studies for identification of trait specific markers in Jatropha.     

Analysis of expressed sequence tags from oil palm (Elaeis guineensis)

This is the first report of a systematic study of genes expressed by means of expressed sequence tag (EST) analysis in oil palm, a species of the Arecales order, a phylogenetically key clade of monocotyledons that is not widely represented in the sequence databases. Five different cDNA libraries were generated from male and female inflorescences, shoot apices and zygotic embryos and unidirectional systematic sequencing was performed. A total of 2411 valid EST sequences were thus obtained. Cluster analysis enabled the identification of 209 groups of related sequences and 1874 singletons. Putative functions were assigned to 1252 of the set of 2083 non-redundant ESTs obtained. The EST database described here is a first step towards gene discovery and cDNA array-based expression analysis in oil palm.     

ANALYSIS OF EXPRESSED SEQUENCE TAGS (ESTS) FROM THE POLAR DIATOM FRAGILARIOPSIS CYLINDRUS1

Analysis of expressed sequence tags (ESTs) was performed to gain insights into cold adaptation in the polar diatom Fragilariopsis cylindrus Grunow. The EST library was generated from RNA isolated 5 days after F. cylindrus cells were shifted from approximately +5° C to ?1.8°C. A total of 1376 ESTs were sequenced from a non‐normalized cDNA library and assembled into 996 tentative unique sequences. About 27% of the ESTs displayed similarity (tBLASTX, e‐value of ≤10^?4) to predicted proteins in the centric diatom Thalassiosira pseudonana Hasle & Heindal. Eleven additional algae and plant data bases were used for annotation of sequences not covered by Thalassiosira sequences (7%). Most of the ESTs were similar to genes encoding proteins responsible for translation, ribosomal structure, and biogenesis (3%), followed by genes encoding proteins for amino acid transport and metabolism and post‐translational modifications. Interestingly, 66% of all the EST sequences from F. cylindrus displayed no similarity ( e ‐value ≤10^?4) to sequences from the 12 non‐redundant databases. Even 6 of the 10 strong to moderately expressed sequences in this EST library could not be identified. Adaptation of F. cylindrus to freezing temperatures of seawater may require a complex protein metabolism and possibly also genes, which were highly expressed but still unknown. However, it could also mean that due to low temperatures, there might have been a stronger pressure to adapt amino acid sequences, making it more difficult to identify these unknown sequences and/or that there are still few protist sequences available for comparison.     

Analysis of expressed sequence tags of flower buds in Lotus japonicus.

In order to study gene expression in a reproductive organ, we constructed a cDNA library of mature flower buds in Lotus japonicus, and characterized expressed sequence tags (ESTs) of 842 clones randomly selected. The EST sequences were clustered into 718 non-redundant groups. From BLAST and FASTA search analyses of both protein and DNA databases, 58.5% of the EST groups showed significant sequence similarities to known genes. Several genes encoding these EST clones were identified as pollen-specific genes, such as pectin methylesterase, ascorbate oxidase, and polygalacturonase, and as homologous genes involved in pollen-pistil interaction. Comparison of these EST sequences with those derived from the whole plant of L. japonicus, revealed that 64.8% of EST sequences from the flower buds were not found in EST sequences of the whole plant. Taken together, the EST data from flower buds generated in this study is useful in dissecting gene expression in floral organ of L. japonicus.     

Expressed sequence tag analysis of khat (Catha edulis) provides a putative molecular biochemical basis for the biosynthesis of phenylpropylamino alkaloids

Khat (Catha edulis Forsk.) is a flowering perennial shrub cultivated for its neurostimulant properties resulting mainly from the occurrence of (S)-cathinone in young leaves. The biosynthesis of (S)-cathinone and the related phenylpropylamino alkaloids (1S,2S)-cathine and (1R,2S)-norephedrine is not well characterized in plants. We prepared a cDNA library from young khat leaves and sequenced 4,896 random clones, generating an expressed sequence tag (EST) library of 3,293 unigenes. Putative functions were assigned to > 98% of the ESTs, providing a key resource for gene discovery. Candidates potentially involved at various stages of phenylpropylamino alkaloid biosynthesis from L-phenylalanine to (1S,2S)-cathine were identified.     

EST-based in silico identification and in vitro test of antimicrobial peptides in Brassica napus

Background

Brassica napus is the third leading source of vegetable oil in the world after soybean and oil palm. The accumulation of gene sequences, especially expressed sequence tags (ESTs) from plant cDNA libraries, has provided a rich resource for genes discovery including potential antimicrobial peptides (AMPs). In this study, we used ESTs including those generated from B. napus cDNA libraries of seeds, pathogen-challenged leaves and deposited in the public databases, as a model, to perform in silico identification and consequently in vitro confirmation of putative AMP activities through a highly efficient system of recombinant AMP prokaryotic expression.

Results

In total, 35,788 were generated from cDNA libraries of pathogen-challenged leaves and 187,272 ESTs from seeds of B. napus, and the 644,998 ESTs of B. napus were downloaded from the EST database of PlantGDB. They formed 201,200 unigenes. First, all the known AMPs from the AMP databank (APD2 database) were individually queried against all the unigenes using the BLASTX program. A total of 972 unigenes that matched the 27 known AMP sequences in APD2 database were extracted and annotated using Blast2GO program. Among these unigenes, 237 unigenes from B. napus pathogen-challenged leaves had the highest ratio (1.15 %) in this unigene dataset, which is 13 times that of the unigene datasets of B. napus seeds (0.09 %) and 2.3 times that of the public EST dataset. About 87 % of each EST library was lipid-transfer protein (LTP) (32 % of total unigenes), defensin, histone, endochitinase, and gibberellin-regulated proteins. The most abundant unigenes in the leaf library were endochitinase and defensin, and LTP and histone in the pub EST library. After masking of the repeat sequence, 606 peptides that were orthologous matched to different AMP families were found. The phylogeny and conserved structural motifs of seven AMPs families were also analysed. To investigate the antimicrobial activities of the predicted peptides, 31 potential AMP genes belonging to different AMP families were selected to test their antimicrobial activities after bioinformatics identification. The AMP genes were all optimized according to Escherichia coli codon usage and synthetized through one-step polymerase chain reaction method. The results showed that 28 recombinant AMPs displayed expected antimicrobial activities against E. coli and Micrococcus luteus and Sclerotinia sclerotiorum strains.

Conclusion

The study not only significantly expanded the number of known/predicted peptides, but also contributed to long-term plant genetic improvement for increased resistance to diverse pathogens of B.napus. These results proved that the high-throughput method developed that combined an in silico procedure with a recombinant AMP prokaryotic expression system is considerably efficient for identification of new AMPs from genome or EST sequence databases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1849-x) contains supplementary material, which is available to authorized users.     

Comparison of gene representation between diploid and triploid Paragonimus westermani by expressed sequence tag analyses

Expressed sequence tag (EST) analysis of the diploid and triploid Paragonimus westermani genes was done to have a rapid and informative outlook of the gene-expression profiles of the parasites. Totals of 506 and 505 ESTs were generated from the diploid and triploid P. westermani cDNA libraries. Based on the BLASTx search results of the diploid P. westermani ESTs, 308 (60.9%) matched significantly with formerly identified genes and 198 (39.1%) showed no significant homology in the GenBank database. A similar homology pattern was shown from the triploid EST BLASTx search results with 346 (68.5%) sharing homology with previously identified genes and 159 (31.5%) showing no significant homology. The EST data from both libraries were analyzed and grouped into 9 categories. Comparison of the 2 EST pools revealed high similarities among the categories of the significantly matched genes. Single genes matched repeatedly were also observed in the 2 EST data. Some genes were found that are not yet characterized in P. westermani; these genes were matched by both the diploid and triploid ESTs. Further study of these genes may provide us with more understanding on the parasite's biology and their specific functions in the 2 strains.