NCBI Non-redundant Protein Research Articles

Genome-wide transcriptomic and phylogenetic analyses reveal distinct aluminum-tolerance mechanisms in the aluminum-accumulating species buckwheat (Fagopyrum tataricum).

BackgroundSimilar to common buckwheat (Fagopyrum esculentum), tartary buckwheat (Fagopyrum tataricum) shows a high level of aluminum (Al) tolerance and accumulation. However, the molecular mechanisms for Al detoxification and accumulation are still poorly understood. To begin to elucidate the molecular basis of Al tolerance and accumulation, we used the Illumina high-throughput mRNA sequencing (RNA-seq) technology to conduct a genome-wide transcriptome analysis on both tip and basal segments of the roots exposed to Al.ResultsBy using the Trinity method for the de novo assembly and cap3 software to reduce the redundancy and chimeras of the transcripts, we constructed 39,815 transcripts with an average length of 1184 bp, among which 20,605 transcripts were annotated by BLAST searches in the NCBI non-redundant protein database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that expression of genes involved in the defense of cell wall toxicity and oxidative stress was preferentially induced by Al stress. Our RNA-seq data also revealed that organic acid metabolism was unlikely to be a rate-limiting step for the Al-induced secretion of organic acids in buckwheat. We identified two citrate transporter genes that were highly induced by Al and potentially involved in the release of citrate into the xylem. In addition, three of four conserved Al-tolerance genes were found to be duplicated in tartary buckwheat and display diverse expression patterns.ConclusionsNearly 40,000 high quality transcript contigs were de novo assembled for tartary buckwheat, providing a reference platform for future research work in this plant species. Our differential expression and phylogenetic analysis revealed novel aspects of Al-tolerant mechanisms in buckwheat.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-014-0395-z) contains supplementary material, which is available to authorized users.

Identification of heat shock proteins via transcriptome profiling of tree peony leaf exposed to high temperature.

The tree peony leaf is an important vegetative organ that is sensitive to abiotic stress and particularly to high temperature. This sensitivity affects plant growth and restricts tree peony distribution. However, the transcriptomic information currently available on the peony leaf in public databases is limited. In this study, we sequenced the transcriptomes of peony leaves subjected to high temperature using the Illumina HiSeq TM 2000 platform. We performed de novo assembly of 93,714 unigenes (average length of 639.7 bp). By searching the public databases, 22,323 unigenes and 13,107 unigenes showed significant similarities with proteins in the NCBI non-redundant protein database and SWISS-PROT database (E-value < 1e-5), respectively. We assigned 17,340 unigenes to Gene Ontology categories, and we assigned 7618 unigenes to clusters of orthologous groups for eukaryotic complete genomes. By searching the Kyoto Encyclopedia of Genes and Genomes Pathway database, 8014 unigenes were assigned to 6 main categories, including 290 KEGG pathways. To advance research on improving thermotolerance, we identified 24 potential heat shock protein genes with complete open reading frames from the transcriptomic sequences. This is the first study to characterize the leaf transcriptome of tree peony leaf using high-throughput sequencing. The information obtained from the tree peony leaf is valuable for gene discovery, and the identified heat shock protein genes can be used to improve plant stress-tolerance.

Transcriptome sequencing and analysis of the coconut leaf beetle, Brontispa longissima.

The coconut leaf beetle, Brontispa longissima, is a destructive pest of palm plants. Although its ecological and biological characteristics are well understood, its genetic information remains largely unknown. To advance our understanding of its molecular ba-sis of biology and ecology, we sequenced and analyzed its whole transcriptome by using high-throughput Illumina paired-end sequencing technology. Approximately 8.08 Gb of clean reads were generated in a single run, which were assembled by using Trinity into 41,652 unigenes with an average length of 932 bp. By sequence similarity searches for known proteins, 23,077 (55.4%) unigenes were annotated by BLASTx searches against the NCBI non-redundant protein database. Of the unigenes assembled, 18,153 and 13,733 were assigned to Gene Ontology and Clusters of Orthologous Groups of proteins, respectively. In addition, 10,415 unigenes were mapped onto 247 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. These transcriptomic resources will facilitate gene identification and elucidate the molecular mechanisms of biological and ecological aspects under-lying this palm pest, in order to design a new control strategy.

Transcriptome characterization and SSR discovery in large-scale loach Paramisgurnus dabryanus (Cobitidae, Cypriniformes)

The large-scale loach (Paramisgurnus dabryanus, Cypriniformes) is a bottom-dwelling freshwater species of fish found mainly in eastern Asia. The natural germplasm resources of this important aquaculture species has been recently threatened due to overfishing and artificial propagation. The objective of this study is to obtain the first functional genomic resource and candidate molecular markers for future conservation and breeding research. Illumina paired-end sequencing generated over one hundred million reads that resulted in 71,887 assembled transcripts, with an average length of 1465bp. 42,093 (58.56%) protein-coding sequences were predicted; and 43,837 transcripts had significant matches to NCBI nonredundant protein (Nr) database. 29,389 and 14,419 transcripts were assigned into gene ontology (GO) categories and Eukaryotic Orthologous Groups (KOG), respectively. 22,102 (31.14%) transcripts were mapped to 302 KEGG pathways. In addition, 15,106 candidate SSR markers were identified, with 11,037 pairs of PCR primers designed. 400 primers pairs of SSR selected randomly were validated, of which 364 (91%) pairs of primers were able to produce PCR products. Further test with 41 loci and 20 large-scale loach specimens collected from the four largest lakes in China showed that 36 (87.8%) loci were polymorphic. The transcriptomic profile and SSR repertoire obtained in this study will facilitate population genetic studies and selective breeding of large-scale loach in the future.

Identification of regulatory genes implicated in continuous flowering of longan (Dimocarpus longan L.).

Longan (Dimocarpus longan L.) is a tropical/subtropical fruit tree of significant economic importance in Southeast Asia. However, a lack of transcriptomic and genomic information hinders research on longan traits, such as the control of flowering. In this study, high-throughput RNA sequencing (RNA-Seq) was used to investigate differentially expressed genes between a unique longan cultivar ‘Sijimi’(S) which flowers throughout the year and a more typical cultivar ‘Lidongben’(L) which flowers only once in the season, with the aim of identifying candidate genes associated with continuous flowering. 36,527 and 40,982 unigenes were obtained by de novo assembly of the clean reads from cDNA libraries of L and S cultivars. Additionally 40,513 unigenes were assembled from combined reads of these libraries. A total of 32,475 unigenes were annotated by BLAST search to NCBI non-redundant protein (NR), Swiss-Prot, Clusters of Orthologous Groups (COGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Of these, almost fifteen thousand unigenes were identified as significantly differentially expressed genes (DEGs) by using Reads Per kb per Million reads (RPKM) method. A total of 6,415 DEGs were mapped to 128 KEGG pathways, and 8,743 DEGs were assigned to 54 Gene Ontology categories. After blasting the DEGs to public sequence databases, 539 potential flowering-related DEGs were identified. In addition, 107 flowering-time genes were identified in longan, their expression levels between two longan samples were compared by RPKM method, of which the expression levels of 15 were confirmed by real-time quantitative PCR. Our results suggest longan homologues of SHORT VEGETATIVE PHASE (SVP), GIGANTEA (GI), F-BOX 1 (FKF1) and EARLY FLOWERING 4 (ELF4) may be involved this flowering trait and ELF4 may be a key gene. The identification of candidate genes related to continuous flowering will provide new insight into the molecular process of regulating flowering time in woody plants.

Genes regulating lipid and protein metabolism are highly expressed in mammary gland of lactating dairy goats.

Dairy goats serve as an important source of milk and also fulfill agricultural and economic roles in developing countries. Understanding the genetic background of goat mammary gland is important for research on the regulatory mechanisms controlling tissue function and the synthesis of milk components. We collected tissue at four different stages of goat mammary gland development and generated approximately 25 GB of data from Illumina de novo RNA sequencing. The combined reads were assembled into 51,361 unigenes, and approximately 60.07 % of the unigenes had homology to other proteins in the NCBI non-redundant protein database (NR). Functional classification through eukaryotic Ortholog Groups of Protein (KOG), gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that the unigenes from goat mammary glands are involved in a wide range of biological processes and metabolic pathways, including lipid metabolism and lactose metabolism. The results of qPCR revealed that genes encoding FABP3, FASN, SCD, PLIN2, whey proteins (LALBA and BLG), and caseins (CSN1S1, CSN1S2, CSN2 and CSN3) at 100 and 310 days postpartum increased significantly compared with the non-lactating period. In addition to their role in lipid and protein synthesis, the higher expression at 310 days postpartum could contribute to mammary cell turnover during pregnancy. In conclusion, this is the first study to characterize the complete transcriptome of goat mammary glands and constitutes a comprehensive genomic resource available for further studies of ruminant lactation.

De novo assembly and characterization of the transcriptome of seagrass Zostera marina using Illumina paired-end sequencing.

BackgroundThe seagrass Zostera marina is a monocotyledonous angiosperm belonging to a polyphyletic group of plants that can live submerged in marine habitats. Zostera marina L. is one of the most common seagrasses and is considered a cornerstone of marine plant molecular ecology research and comparative studies. However, the mechanisms underlying its adaptation to the marine environment still remain poorly understood due to limited transcriptomic and genomic data.Principal FindingsHere we explored the transcriptome of Z. marina leaves under different environmental conditions using Illumina paired-end sequencing. Approximately 55 million sequencing reads were obtained, representing 58,457 transcripts that correspond to 24,216 unigenes. A total of 14,389 (59.41%) unigenes were annotated by blast searches against the NCBI non-redundant protein database. 45.18% and 46.91% of the unigenes had significant similarity with proteins in the Swiss-Prot database and Pfam database, respectively. Among these, 13,897 unigenes were assigned to 57 Gene Ontology (GO) terms and 4,745 unigenes were identified and mapped to 233 pathways via functional annotation against the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). We compared the orthologous gene family of the Z. marina transcriptome to Oryza sativa and Pyropia yezoensis and 11,667 orthologous gene families are specific to Z. marina. Furthermore, we identified the photoreceptors sensing red/far-red light and blue light. Also, we identified a large number of genes that are involved in ion transporters and channels including Na+ efflux, K+ uptake, Cl− channels, and H+ pumping.ConclusionsOur study contains an extensive sequencing and gene-annotation analysis of Z. marina. This information represents a genetic resource for the discovery of genes related to light sensing and salt tolerance in this species. Our transcriptome can be further utilized in future studies on molecular adaptation to abiotic stress in Z. marina.

A comprehensive transcriptome provides candidate genes for sex determination/differentiation and SSR/SNP markers in yellow catfish.

Sex dimorphic growth pattern has significant theory and application implications in fish. Recently, a Y- and X-specific allele marker-assisted sex control technique has been developed for mass production of all-male population in yellow catfish (Pelteobagrus fulvidraco), but the genetic information for sex determination and sex control breeding has remained unclear. Here, we attempted to provide the first insight into a comprehensive transcriptome covering multiple tissues from XX females, XY males, and YY super-males of yellow catfish by using 454 GS-FLX platform, for a better assembly and gene coverage. A total of 1,202,933 high quality reads (about 540Mbp) were obtained and assembled into 28,297 contigs and 141,951 singletons. BLASTX searches against the NCBI non-redundant protein database (nr) led a total of 52,564 unique sequences including 18,748 contigs and 33,816 singletons to match 25,669 known or predicted unique proteins. All of them with annotated function were categorized by gene ontology (GO) analysis, and 712 were assigned to reproduction and reproductive process. Some potential genes relevant to reproductive system including steroid hormone biosynthesis and GnRH (gonadotropin-releasing hormone) signaling pathway were further identified by Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis; and at least 21 sex determination and differentiation-related genes, such as Dmrt1, Sox9a/b, Cyp19b, WT1, and AMH were identified and characterized. Additionally, a total of 82,794 simple sequence repeats (SSRs), 26,450 single nucleotide polymorphisms (SNPs), and 4,145 insertions and deletions (INDELs) were revealed from the transcriptome data. Therefore, the current transcriptome resources highlight further studies on sex-control breeding in yellow catfish and will benefit future studies on reproduction and sex determination in teleost fish.

De novo characterization of the Lycium chinense Mill. leaf transcriptome and analysis of candidate genes involved in carotenoid biosynthesis

Lycium chinense Mill. (Chinese wolfberry), enriching in carotenoids, is an important Chinese herbal medicine. However, studies on the functional genomics research, especially the carotenoid biosynthesis and accumulation, are limited because of insufficiently available datasets. RNA-Seq was performed by the Illumina sequencing platform. Approximately 26 million clean reads were generated after filtering. Clean reads were assembled by SOAPdenovo and subsequently annotated. Among all 61,595 unigenes, 37,816 (61.39%), 25,266 (41.02%), and 17,598 (28.57%) unigenes were annotated in NCBI non-redundant protein, Swiss-Prot, and Kyoto Encyclopedia of Genes and Genomes (KEGG) database, respectively. A total of 16,073 and 11,394 unigenes were assigned to Gene Ontology and Cluster of Orthologous Group, respectively. Furthermore, the majority of genes encoding the enzymes in the carotenoid biosynthesis pathway were identified in the unigene datasets. We first found several genes related to L. chinense carotenoid biosynthesis. The expression levels and the biological functions of these genes involved in carotenoid biosynthesis in the leaf and the green ripening fruit were further confirmed by qPCR and high performance liquid chromatography (HPLC). In the present study, we first characterized the transcriptome of L. chinense leaf, which may provide useful data for functional genomics investigations in L. chinense in the future. And essential genes involved in the carotenoid biosynthesis pathway may contribute to elucidate the expression patterns in different stages of development and fruit ripening and the specific mechanisms of carotenoid biosynthesis/accumulation in L. chinense.

De novo assembly of Eugenia uniflora L. transcriptome and identification of genes from the terpenoid biosynthesis pathway

Pitanga (Eugenia uniflora L.) is a member of the Myrtaceae family and is of particular interest due to its medicinal properties that are attributed to specialized metabolites with known biological activities. Among these molecules, terpenoids are the most abundant in essential oils that are found in the leaves and represent compounds with potential pharmacological benefits. The terpene diversity observed in Myrtaceae is determined by the activity of different members of the terpene synthase and oxidosqualene cyclase families. Therefore, the aim of this study was to perform a de novo assembly of transcripts from E. uniflora leaves and to annotation to identify the genes potentially involved in the terpenoid biosynthesis pathway and terpene diversity. In total, 72,742 unigenes with a mean length of 1048bp were identified. Of these, 43,631 and 36,289 were annotated with the NCBI non-redundant protein and Swiss-Prot databases, respectively. The gene ontology categorized the sequences into 53 functional groups. A metabolic pathway analysis with KEGG revealed 8,625 unigenes assigned to 141 metabolic pathways and 40 unigenes predicted to be associated with the biosynthesis of terpenoids. Furthermore, we identified four putative full-length terpene synthase genes involved in sesquiterpenes and monoterpenes biosynthesis, and three putative full-length oxidosqualene cyclase genes involved in the triterpenes biosynthesis. The expression of these genes was validated in different E. uniflora tissues.

Characterization of the global transcriptome using Illumina sequencing and novel microsatellite marker information in seashore paspalum

Seashore paspalum (Paspalum vaginatum O. Swartz) is a halophytic warm-season perennial grass, which plays an important role in protection of environment. However, very limited transcriptomic and genomic information are currently available for seashore paspalum. In this study, the seashore paspalum transcriptome were sequenced using Illumina paired-end sequencing technology. In total, 32,603 unigenes with an average length of 970 bp were obtained by de novo assembly. Of the unigenes, 25,411 unigenes (77.94 %) had significant similarity with known proteins in the NCBI non-redundant protein and Swiss-Prot databases. Of these annotated unigenes, 20,962 and 10,620 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. A total of 4,699 unigenes were mapped into 118 Kyoto Encyclopedia of Genes and Genomes pathways. In addition, 3,010 potential simple sequence repeats (SSRs) were predicted for microsatellite analyse. Tri-nucleotide was the dominant repeat (1,583, 52.58 %), followed by di-nucleotide (975, 32.39 %) and tetra-nucleotide (217, 7.21 %). Fifty SSR sites were randomly selected for validation and development of EST–SSR markers. This study provided the global sequence data for seashore paspalum and demonstrated that the Illumina paired-end sequencing is a fast and cost-effective approach to gene discovery and molecular marker development.

Transcriptome sequencing and identification of cold tolerance genes in hardy Corylus species (C. heterophylla Fisch) floral buds.

BackgroundThe genus Corylus is an important woody species in Northeast China. Its products, hazelnuts, constitute one of the most important raw materials for the pastry and chocolate industry. However, limited genetic research has focused on Corylus because of the lack of genomic resources. The advent of high-throughput sequencing technologies provides a turning point for Corylus research. In the present study, we performed de novo transcriptome sequencing for the first time to produce a comprehensive database for the Corylus heterophylla Fisch floral buds.ResultsThe C. heterophylla Fisch floral buds transcriptome was sequenced using the Illumina paired-end sequencing technology. We produced 28,930,890 raw reads and assembled them into 82,684 contigs. A total of 40,941 unigenes were identified, among which 30,549 were annotated in the NCBI Non-redundant (Nr) protein database and 18,581 were annotated in the Swiss-Prot database. Of these annotated unigenes, 25,311 and 10,514 unigenes were assigned to gene ontology (GO) categories and clusters of orthologous groups (COG), respectively. We could map 17,207 unigenes onto 128 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) database. Additionally, based on the transcriptome, we constructed a candidate cold tolerance gene set of C. heterophylla Fisch floral buds. The expression patterns of selected genes during four stages of cold acclimation suggested that these genes might be involved in different cold responsive stages in C. heterophylla Fisch floral buds.ConclusionThe transcriptome of C. heterophylla Fisch floral buds was deep sequenced, de novo assembled, and annotated, providing abundant data to better understand the C. heterophylla Fisch floral buds transcriptome. Candidate genes potentially involved in cold tolerance were identified, providing a material basis for future molecular mechanism analysis of C. heterophylla Fisch floral buds tolerant to cold stress.

Identification of novel cupredoxin homologs using overlapped conserved residues based approach.

Cupredoxin-like proteins are mainly copper-binding proteins that conserve a typical rigid Greek-key arrangement consisting of an eight-stranded β-sandwich, even though they share as little as 10-15% sequence similarity. The electron transport function of the Cupredoxins is critical for respiration and photosynthesis, and the proteins have therapeutic potential. Despite their crucial biological functions, the identification of the distant Cupredoxin homologs has been a difficult task due to their low sequence identity. In this study, the overlapped conserved residue (OCR) fingerprint for the Cupredoxin superfamily, which consists of conserved residues in three aspects (i.e., the sequence, structure, and intramolecular interaction), was used to detect the novel Cupredoxin homologs in the NCBI non-redundant protein sequence database. The OCR fingerprint could identify 54 potential Cupredoxin sequences, which were validated by scanning them against the conserved Cupredoxin motif near the Cu-binding site. This study also attempted to model the 3D structures and to predict the functions of the identified potential Cupredoxins. This study suggests that the OCR-based approach can be used efficiently to detect novel homologous proteins with low sequence identity, such as Cupredoxins.

ESTs analysis of putative genes engaged in Polyporus umbellatus sclerotial development.

Polyporus umbellatus is one of the most widely used and precious medicinal fungi and the underground sclerotia are known to be with great medicinal value. However, the molecular mechanisms involved in sclerotial development are poorly understood. In the present study, we constructed a forward suppression subtractive hybridization (SSH) cDNA library of Polyporus umbellatus to identify genes expressing differently between mycelium and sclerotia. In this library, a total of 1202 clones were sequenced, assembled into 222 contigs and 524 singletons which were further searched against the NCBI nonredundant (NR) protein database (E-value cutoff, 10−5). Based on sequence similarity with known proteins, 378 sequences between mycelium and sclerotial were identified and classified into different functional categories through Gene Ontology (GO), Clusters of orthologous Groups of proteins (COGs). We have finally identified a majority of differentially expressed genes (constituting 5.6% of the present library) between the two different periods. An expression level of 32 selected expressed sequence tags (ESTs) generated from the above SSH cDNA library was studied through RT-PCR. This study provides the first global overview of genes putatively involved in Polyporus umbellatus sclerotial development and provides a preliminary basis for further functional research in terms of regulated gene expression in sclerotial production.

Microbial community diversity and physical-chemical features of the Southwestern Atlantic Ocean.

Microbial oceanography studies have demonstrated the central role of microbes in functioning and nutrient cycling of the global ocean. Most of these former studies including at Southwestern Atlantic Ocean (SAO) focused on surface seawater and benthic organisms (e.g., coral reefs and sponges). This is the first metagenomic study of the SAO. The SAO harbors a great microbial diversity and marine life (e.g., coral reefs and rhodolith beds). The aim of this study was to characterize the microbial community diversity of the SAO along the depth continuum and different water masses by means of metagenomic, physical-chemical and biological analyses. The microbial community abundance and diversity appear to be strongly influenced by the temperature, dissolved organic carbon, and depth, and three groups were defined [1. surface waters; 2. sub-superficial chlorophyll maximum (SCM) (48-82 m) and 3. deep waters (236-1,200 m)] according to the microbial composition. The microbial communities of deep water masses [South Atlantic Central water, Antarctic Intermediate water and Upper Circumpolar Deep water] are highly similar. Of the 421,418 predicted genes for SAO metagenomes, 36.7 % had no homologous hits against 17,451,486 sequences from the North Atlantic, South Atlantic, North Pacific, South Pacific and Indian Oceans. From these unique genes from the SAO, only 6.64 % had hits against the NCBI non-redundant protein database. SAO microbial communities share genes with the global ocean in at least 70 cellular functions; however, more than a third of predicted SAO genes represent a unique gene pool in global ocean. This study was the first attempt to characterize the taxonomic and functional community diversity of different water masses at SAO and compare it with the microbial community diversity of the global ocean, and SAO had a significant portion of endemic gene diversity. Microbial communities of deep water masses (236-1,200 m) are highly similar, suggesting that these water masses have very similar microbiological attributes, despite the common knowledge that water masses determine prokaryotic community and are barriers to microbial dispersal. The present study also shows that SCM is a clearly differentiated layer within Tropical waters with higher abundance of phototrophic microbes and microbial diversity.

A comprehensive reference transcriptome resource for the common house spider Parasteatoda tepidariorum.

Parasteatoda tepidariorum is an increasingly popular model for the study of spider development and the evolution of development more broadly. However, fully understanding the regulation and evolution of P. tepidariorum development in comparison to other animals requires a genomic perspective. Although research on P. tepidariorum has provided major new insights, gene analysis to date has been limited to candidate gene approaches. Furthermore, the few available EST collections are based on embryonic transcripts, which have not been systematically annotated and are unlikely to contain transcripts specific to post-embryonic stages of development. We therefore generated cDNA from pooled embryos representing all described embryonic stages, as well as post-embryonic stages including nymphs, larvae and adults, and using Illumina HiSeq technology obtained a total of 625,076,514 100-bp paired end reads. We combined these data with 24,360 ESTs available in GenBank, and 1,040,006 reads newly generated from 454 pyrosequencing of a mixed-stage embryo cDNA library. The combined sequence data were assembled using a custom de novo assembly strategy designed to optimize assembly product length, number of predicted transcripts, and proportion of raw reads incorporated into the assembly. The de novo assembly generated 446,427 contigs with an N50 of 1,875 bp. These sequences obtained 62,799 unique BLAST hits against the NCBI non-redundant protein data base, including putative orthologs to 8,917 Drosophila melanogaster genes based on best reciprocal BLAST hit identity compared with the D. melanogaster proteome. Finally, we explored the utility of the transcriptome for RNA-Seq studies, and showed that this resource can be used as a mapping scaffold to detect differential gene expression in different cDNA libraries. This resource will therefore provide a platform for future genomic, gene expression and functional approaches using P. tepidariorum.

De novo transcriptome assembly from inflorescence of Orchis italica: analysis of coding and non-coding transcripts.

The floral transcriptome of Orchis italica, a wild orchid species, was obtained using Illumina RNA-seq technology and specific de novo assembly and analysis tools. More than 100 million raw reads were processed resulting in 132,565 assembled transcripts and 86,079 unigenes with an average length of 606 bp and N50 of 956 bp. Functional annotation assigned 38,984 of the unigenes to records present in the NCBI non-redundant protein database, 32,161 of them to Gene Ontology terms, 15,775 of them to Eukaryotic Orthologous Groups (KOG) and 7,143 of them to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The in silico expression analysis based on the Fragments Per Kilobase of transcript per Million mapped reads (FPKM) was confirmed by real-time RT-PCR experiments on 10 selected unigenes, which showed high and statistically significant positive correlation with the RNA-seq based expression data. The prediction of putative long non-coding RNAs was assessed using two different software packages, CPC and Portrait, resulting in 7,779 unannotated unigenes that matched the threshold values for both of the analyses. Among the predicted long non-coding RNAs, one is the homologue of TAS3, a long non-coding RNA precursor of trans-acting small interfering RNAs (ta-siRNAs). The differential expression pattern observed for the selected putative long non-coding RNAs suggests their possible functional role in different floral tissues.

Transcriptome analysis of the differences in gene expression between testis and ovary in green mud crab (Scylla paramamosain).

BackgroundThe green mud crab (Scylla paramamosain) is the most prevalent crustacean on the southeast coast of China. The molecular regulatory mechanism of sex determination and gonadal differentiation in this species has received considerable attention in recent years because of the huge differences—both biological and economic—between male and female crabs. In this study, next-generation sequencing technology was used to develop deep-coverage transcriptomic sequencing data for the testis and ovary of S. paramamosain.ResultsA total of 365,116 reads (testis 171,962, ovary 193,154) with an average sequence length of 285 bp were produced from testis and ovary cDNA libraries. After filtering out contaminating reads, the clean reads were assembled, producing a total of 21,791 isotigs and leaving 22,814 reads as singlets. Using the BLASTX program, 3,471 unique sequences (2,275 isotigs and 1,196 singletons) were annotated with known protein sequences from the NCBI non-redundant (Nr) protein sequence database. The Gene Ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses allowed the 224 unique sequences that were annotated with enzyme code (EC) numbers to be mapped into 174 KEGG pathways. After comparing the ovary and testis libraries, 4,021 gonad-differentially, 10,522 ovary-specifically, and 19,013 testis-specifically expressed genes were identified. Moreover, 33 ovary-specific, 14 testis-specific, and 34 gonad-differential transcripts were confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, 8,610 putative simple sequence repeats (SSRs) and 23,879 potential single nucleotide polymorphisms (SNPs) were identified.ConclusionThis is the first large-scale RNA sequencing of S. paramamosain to be reported. We have identified many important functional genes and made a preliminary attempt to construct the regulatory network involved in the gonadal development of crustaceans. The annotated transcriptome data will provide fundamental support for future research into the reproduction biology of S. paramamosain. A large number of candidate SSRs and SNPs were detected, which could be used as genetic markers for population genetics and functional genomics in this species.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-585) contains supplementary material, which is available to authorized users.

Evidence for Proteomic and Metabolic Adaptations Associated with Alterations of Seed Yield and Quality in Sulfur-limited Brassica napus L

In Brassica napus, seed yield and quality are related to sulfate availability, but the seed metabolic changes in response to sulfate limitation remain largely unknown. To address this question, proteomics and biochemical studies were carried out on mature seeds obtained from plants grown under low sulfate applied at the bolting (LS32), early flowering (LS53), or start of pod filling (LS70) stage. The protein quality of all low-sulfate seeds was reduced and associated with a reduction of S-rich seed storage protein accumulation (as Cruciferin Cru4) and an increase of S-poor seed storage protein (as Cruciferin BnC1). This compensation allowed the protein content to be maintained in LS70 and LS53 seeds but was not sufficient to maintain the protein content in LS32 seeds. The lipid content and quality of LS53 and LS32 seeds were also affected, and these effects were primarily associated with a reduction of C18-derivative accumulation. Proteomics changes related to lipid storage, carbohydrate metabolism, and energy (reduction of caleosins, phosphoglycerate kinase, malate synthase, ATP-synthase β-subunit, and thiazole biosynthetic enzyme THI1 and accumulation of β-glucosidase and citrate synthase) provide insights into processes that may contribute to decreased oil content and altered lipid composition (in favor of long-chain fatty acids in LS53 and LS32 seeds). These data indicate that metabolic changes associated with S limitation responses affect seed storage protein composition and lipid quality. Proteins involved in plant stress response, such as dehydroascorbate reductase and Cu/Zn-superoxide dismutase, were also accumulated in LS53 and LS32 seeds, and this might be a consequence of reduced glutathione content under low S availability. LS32 treatment also resulted in (i) reduced germination vigor, as evidenced by lower germination indexes, (ii) reduced seed germination capacity, related to a lower seed viability, and (iii) a strong decrease of glyoxysomal malate synthase, which is essential for the use of fatty acids during seedling establishment.

Towards Decrypting Cryptobiosis—Analyzing Anhydrobiosis in the Tardigrade Milnesium tardigradum Using Transcriptome Sequencing

BackgroundMany tardigrade species are capable of anhydrobiosis; however, mechanisms underlying their extreme desiccation resistance remain elusive. This study attempts to quantify the anhydrobiotic transcriptome of the limno-terrestrial tardigrade Milnesium tardigradum.ResultsA prerequisite for differential gene expression analysis was the generation of a reference hybrid transcriptome atlas by assembly of Sanger, 454 and Illumina sequence data. The final assembly yielded 79,064 contigs (>100 bp) after removal of ribosomal RNAs. Around 50% of them could be annotated by SwissProt and NCBI non-redundant protein sequences. Analysis using CEGMA predicted 232 (93.5%) out of the 248 highly conserved eukaryotic genes in the assembly. We used this reference transcriptome for mapping and quantifying the expression of transcripts regulated under anhdydrobiosis in a time-series during dehydration and rehydration. 834 of the transcripts were found to be differentially expressed in a single stage (dehydration/inactive tun/rehydration) and 184 were overlapping in two stages while 74 were differentially expressed in all three stages. We have found interesting patterns of differentially expressed transcripts that are in concordance with a common hypothesis of metabolic shutdown during anhydrobiosis. This included down-regulation of several proteins of the DNA replication and translational machinery and protein degradation. Among others, heat shock proteins Hsp27 and Hsp30c were up-regulated in response to dehydration and rehydration. In addition, we observed up-regulation of ployubiquitin-B upon rehydration together with a higher expression level of several DNA repair proteins during rehydration than in the dehydration stage.ConclusionsMost of the transcripts identified to be differentially expressed had distinct cellular function. Our data suggest a concerted molecular adaptation in M. tardigradum that permits extreme forms of ametabolic states such as anhydrobiosis. It is temping to surmise that the desiccation tolerance of tradigrades can be achieved by a constitutive cellular protection system, probably in conjunction with other mechanisms such as rehydration-induced cellular repair.