WEGO: a web tool for plotting GO annotations

Proteomic analysis of Trichinella spiralis proteins in intestinal epithelial cells after culture with their larvae by shotgun LC-MS/MS approach

Although it has been known for many years that Trichinella spiralis initiates infection by invading intestinal epithelium, the mechanisms by which the parasite invades the intestinal epithelium are unknown. The purpose of this study was to screen the invasion-related proteins among the increased proteins of intestinal epithelial cells after culture with T. spiralis and to study their molecular functions. The proteins of HCT-8 cells which cultured with T. spiralis infective larvae were analyzed by SDS-PAGE and Western blot. Results showed that compared with proteins of normal HCT-8 cells, four additional protein bands (115, 61, 35 and 24 kDa) of HCT-8 cells cultured with the infective larvae were recognized by sera of the mice infected with T. spiralis, which may be the invasion-related proteins released by the infective larvae. Three bands (61, 35 and 24 kDa) were studied employing shotgun LC-MS/MS. Total 64 proteins of T. spiralis were identified from T. spiralis protein database by using SEQUEST searches, of which 43 (67.2%) proteins were distributed in a range of 10-70 kDa, and 26 proteins (40.6%) were in the range of pI 5-6. Fifty-four proteins were annotated according to Gene Ontology Annotation in terms of molecular function, biological process, and cellular localization. Out of 54 annotated proteins, 43 proteins (79.6%) had binding activity and 23 proteins (42.6%) had catalytic activity (e.g. hydrolase, transferase, etc.), which might be related to the invasion of intestinal epithelial cells by T. spiralis. The protein profile provides a valuable basis for further studies of the invasion-related proteins of T. spiralis.

Pseudorabies virus infected porcine epithelial cell line generates a diverse set of host microRNAs and a special cluster of viral microRNAs

Pseudorabies virus (PRV) belongs to Alphaherpesvirinae subfamily that causes huge economic loss in pig industry worldwide. It has been recently demonstrated that many herpesviruses encode microRNAs (miRNAs), which play crucial roles in viral life cycle. However, the knowledge about PRV-encoded miRNAs is still limited. Here, we report a comprehensive analysis of both viral and host miRNA expression profiles in PRV-infected porcine epithelial cell line (PK-15). Deep sequencing data showed that the ∼4.6 kb intron of the large latency transcript (LLT) functions as a primary microRNA precursor (pri-miRNA) that encodes a cluster of 11 distinct miRNAs in the PRV genome, and 209 known and 39 novel porcine miRNAs were detected. Viral miRNAs were further confirmed by stem-loop RT-PCR and northern blot analysis. Intriguingly, all of these viral miRNAs exhibited terminal heterogeneity both at the 5′ and 3′ ends. Seven miRNA genes produced mature miRNAs from both arms and two of the viral miRNA genes showed partially overlapped in their precursor regions. Unexpectedly, a terminal loop-derived small RNA with high abundance and one special miRNA offset RNA (moRNA) were processed from a same viral miRNA precursor. The polymorphisms of viral miRNAs shed light on the complexity of host miRNA-processing machinery and viral miRNA-regulatory mechanism. The swine genes and PRV genes were collected for target prediction of the viral miRNAs, revealing a complex network formed by both host and viral genes. GO enrichment analysis of host target genes suggests that PRV miRNAs are involved in complex cellular pathways including cell death, immune system process, metabolic pathway, indicating that these miRNAs play significant roles in virus-cells interaction of PRV and its hosts. Collectively, these data suggest that PRV infected epithelial cell line generates a diverse set of host miRNAs and a special cluster of viral miRNAs, which might facilitate PRV replication in cells.

Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-Seq

Background

Chinese bayberry (Myrica rubra Sieb. and Zucc.) is an important subtropical fruit crop and an ideal species for fruit quality research due to the rapid and substantial changes that occur during development and ripening, including changes in fruit color and taste. However, research at the molecular level is limited by a lack of sequence data. The present study was designed to obtain transcript sequence data and examine gene expression in bayberry developing fruit based on RNA-Seq and bioinformatic analysis, to provide a foundation for understanding the molecular mechanisms controlling fruit quality changes during ripening.

Results

RNA-Seq generated 1.92 G raw data, which was then de novo assembled into 41,239 UniGenes with a mean length of 531 bp. Approximately 80% of the UniGenes (32,805) were annotated against public protein databases, and coding sequences (CDS) of 31,665 UniGenes were determined. Over 3,600 UniGenes were differentially expressed during fruit ripening, with 826 up-regulated and 1,407 down-regulated. GO comparisons between the UniGenes of these two types and interactive pathways (Ipath) analysis found that energy-related metabolism was enhanced, and catalytic activity was increased. All genes involved in anthocyanin biosynthesis were up-regulated during the fruit ripening processes, concurrent with color change. Important changes in carbohydrate and acid metabolism in the ripening fruit are likely associated with expression of sucrose phosphate synthase (SPS) and glutamate decarboxylase (GAD).

Conclusions

Mass sequence data of Chinese bayberry was obtained and the expression profiles were examined during fruit ripening. The UniGenes were annotated, providing a platform for functional genomic research with this species. Using pathway mapping and expression profiles, the molecular mechanisms for changes in fruit color and taste during ripening were examined. This provides a reference for the study of complicated metabolism in non-model perennial species.

Heat acclimation and the role of RpoS in prolonged heat shock of Escherichia coli O157

Escherichia coli, a commensal mesophile that primarily inhabits the gastro-intestinal tract, responds to temperature up-shifts with transient expression of stress-response proteins. The goal of this study was to identify adaptive proteins of E. coli O157 crucial for growth resumption of this human pathogen after heat shock, with specific focus on the role of the RpoS sigma factor. Using the comparative proteomic analysis of hyper-thermally acclimatized wild-type strain B-1 and rpoS-mutant strain SV521, we identified 39 proteins that underwent significantly-different induction upon temperature shock at 45°C or rpoS mutation. All identified proteins of the heat post-acclimation stimulon fell into two large sub-groups: (i) stress proteins, including molecular chaperons, proteases, DNA/RNA stabilizing enzymes, and anti-oxidant proteins, and (ii) housekeeping proteins. It was found that in the heat stress stimulon RpoS has significantly (P=0.012) limited control over the key stress proteins involved in translation, translational elongation, protein folding and refolding. However, RpoS showed a significant (P=0.035) control over the cellular metabolic processes that included NADPH regeneration, pentose-phosphate shunt, nicotinamide nucleotide and NADP metabolic processes, reflecting its specific importance in promoting resource utilization (energy, protein synthesis etc.) during proliferation of hyperthermally-adapted cells. Pathogenic strains, like E. coli O157, have the ability to survive a variety of harsh stress conditions, leading to their entry into the food chain, and subsequent pathogenesis. This research offers insights into the physiological response of this pathogen during the critical period following adaptation to thermal stress and subsequent resumption of growth.

Gene discovery in the threatened elkhorn coral: 454 sequencing of the Acropora palmata transcriptome

Background

Cnidarians, including corals and anemones, offer unique insights into metazoan evolution because they harbor genetic similarities with vertebrates beyond that found in model invertebrates and retain genes known only from non-metazoans. Cataloging genes expressed in Acropora palmata, a foundation-species of reefs in the Caribbean and western Atlantic, will advance our understanding of the genetic basis of ecologically important traits in corals and comes at a time when sequencing efforts in other cnidarians allow for multi-species comparisons.

Results

A cDNA library from a sample enriched for symbiont free larval tissue was sequenced on the 454 GS-FLX platform. Over 960,000 reads were obtained and assembled into 42,630 contigs. Annotation data was acquired for 57% of the assembled sequences. Analysis of the assembled sequences indicated that 83–100% of all A. palmata transcripts were tagged, and provided a rough estimate of the total number genes expressed in our samples (∼18,000–20,000). The coral annotation data contained many of the same molecular components as in the Bilateria, particularly in pathways associated with oxidative stress and DNA damage repair, and provided evidence that homologs of p53, a key player in DNA repair pathways, has experienced selection along the branch separating Cnidaria and Bilateria. Transcriptome wide screens of paralog groups and transition/transversion ratios highlighted genes including: green fluorescent proteins, carbonic anhydrase, and oxidative stress proteins; and functional groups involved in protein and nucleic acid metabolism, and the formation of structural molecules. These results provide a starting point for study of adaptive evolution in corals.

Conclusions

Currently available transcriptome data now make comparative studies of the mechanisms underlying coral's evolutionary success possible. Here we identified candidate genes that enable corals to maintain genomic integrity despite considerable exposure to genotoxic stress over long life spans, and showed conservation of important physiological pathways between corals and bilaterians.

Analysis of expression sequence tags from a full-length-enriched cDNA library of developing sesame seeds (Sesamum indicum)

BACKGROUND

Sesame (Sesamum indicum) is one of the most important oilseed crops with high oil contents and rich nutrient value. However, genetic improvement efforts in sesame could not get benefit from molecular biology technology due to poor DNA and RNA sequence resources. In this study, we carried out a large scale of expressed sequence tags (ESTs) sequencing from developing sesame seeds and further conducted analysis on seed storage products-related genes.

RESULTS

A normalized and full-length enriched cDNA library from 5 ~ 30 days old immature seeds was constructed and randomly sequenced, leading to generation of 41,248 expressed sequence tags (ESTs) which then formed 4,713 contigs and 27,708 singletons with 44.9% uniESTs being putative full-length open reading frames. Approximately 26,091 of all these uniESTs have significant matches to the counterparts in Nr database of GenBank, and 21,628 of them were assigned to one or more Gene ontology (GO) terms. Homologous genes involved in oil biosynthesis were identified including some conservative transcription factors regulating oil biosynthesis such as LEAFY COTYLEDON1 (LEC1), PICKLE (PKL), WRINKLED1 (WRI1) and majority of them were found for the first time in sesame seeds. One hundred and 17 ESTs were identified possibly involved in biosynthesis of sesame lignans, sesamin and sesamolin. In total, 9,347 putative functional genes from developing seeds were identified, which accounts for one third of total genes in the sesame genome. Further analysis of the uniESTs identified 1,949 non-redundant simple sequence repeats (SSRs).

CONCLUSIONS

This study has provided an overview of genes expressed during sesame seed development. This collection of sesame full-length cDNAs covered a wide variety of genes in seeds, in particular, candidate genes involved in biosynthesis of sesame oils and lignans. These EST sequences enriched with full length will contribute to comparative genomic studies on sesame and other oilseed plants and serve as an abundant information platform for functional marker development and functional gene study.

Transcriptomics of a giant freshwater prawn (Macrobrachium rosenbergii): de novo assembly, annotation and marker discovery

BACKGROUND

Giant freshwater prawn (Macrobrachium rosenbergii or GFP), is the most economically important freshwater crustacean species. However, as little is known about its genome, 454 pyrosequencing of cDNA was undertaken to characterise its transcriptome and identify genes important for growth.

METHODOLOGY AND PRINCIPAL FINDINGS

A collection of 787,731 sequence reads (244.37 Mb) obtained from 454 pyrosequencing analysis of cDNA prepared from muscle, ovary and testis tissues taken from 18 adult prawns was assembled into 123,534 expressed sequence tags (ESTs). Of these, 46% of the 8,411 contigs and 19% of 115,123 singletons possessed high similarity to sequences in the GenBank non-redundant database, with most significant (E value < 1e(-5)) contig (80%) and singleton (84%) matches occurring with crustacean and insect sequences. KEGG analysis of the contig open reading frames identified putative members of several biological pathways potentially important for growth. The top InterProScan domains detected included RNA recognition motifs, serine/threonine-protein kinase-like domains, actin-like families, and zinc finger domains. Transcripts derived from genes such as actin, myosin heavy and light chain, tropomyosin and troponin with fundamental roles in muscle development and construction were abundant. Amongst the contigs, 834 single nucleotide polymorphisms, 1198 indels and 658 simple sequence repeats motifs were also identified.

CONCLUSIONS

The M. rosenbergii transcriptome data reported here should provide an invaluable resource for improving our understanding of this species' genome structure and biology. The data will also instruct future functional studies to manipulate or select for genes influencing growth that should find practical applications in aquaculture breeding programs.

Pyrosequencing the midgut transcriptome of the brown planthopper, Nilaparvata lugens

The brown planthopper, Nilaparvata lugens, is a serious pest threatening rice production across the world. To identify the main features of the gene expression and the key components of the midgut of N. lugens responsible for nutrition, xenobiotic metabolism and the immune response, we used pyrosequencing to sample the transcriptome. More than 190,000 clean sequences were generated, which led to about 30,000 unique sequences. Sequence analysis indicated that genes with abundant transcripts in the midgut of N. lugens were mainly sugar hydrolyases and transporters, proteases and detoxification-related proteins. Based on the sequence information, we cloned the candidate sucrase gene; this enzyme is likely to interact with the perimicrovillar membrane through its highly hydrophobic C-terminal region. Many proteases were identified, which supported the hypothesis that N. lugens uses the proteolysis system for digestion. Scores of detoxification genes were newly identified, including cytochrome P450s, glutathione S-transferases, caroxylesterases. A wealth of new transcripts possibly participating in the immune response were described as well. The gene encoding a peptidoglycan recognition protein was cloned. Unlike in Acyrthosiphon pisum, the immunodeficiency pathway may be present in N. lugens. This is the first global analysis of midgut transcriptome from N. lugens.

Comparative proteomic analysis of challenged Zhikong scallop (Chlamys farreri): a new insight into the anti-Vibrio immune response of marine bivalves

The current studies on molecular mechanism of bivalves interacting with bacteria are mainly on mRNA and recombinant protein levels. These works provide little information on natural proteins, which limit further understandings. In this study, we conducted a pioneer work to promote researches on the anti-Vibrio immune response of Zhikong Scallop Chlamys farreri through proteomic techniques. Firstly a reference map was constructed for the hepatopancreas of C. farreri. Totally 65 protein spots were included in the reference map, while 46 of them were identified. Gene ontology analysis revealed high activities of metabolism and immunity in hepatopancreas. Furthermore, hepatopancreas of C. farreri injected with Vibrio harveyi at 24 h post-injection (hpi) were used for comparative proteomic analysis. Totally 27 differentially expressed proteins spots after challenge were screened; and 15 were successfully identified. These proteins include some immune-related proteins, metabolism enzymes, and new molecules which were not paid attentions in previous immunity studies in C. farreri. The results indicated that molecular chaperons and the antioxidant system are key elements in the anti-Vibrio immune response of hepatopancreas of C. farreri. The identification of new molecules provides indications for further studies. The results of this work provide a new insight into the anti-Vibrio immune response of marine bivalves.

From model to crop: functional analysis of a STAY-GREEN gene in the model legume Medicago truncatula and effective use of the gene for alfalfa improvement

Medicago truncatula has been developed into a model legume. Its close relative alfalfa (Medicago sativa) is the most widely grown forage legume crop in the United States. By screening a large population of M. truncatula mutants tagged with the transposable element of tobacco (Nicotiana tabacum) cell type1 (Tnt1), we identified a mutant line (NF2089) that maintained green leaves and showed green anthers, central carpels, mature pods, and seeds during senescence. Genetic and molecular analyses revealed that the mutation was caused by Tnt1 insertion in a STAY-GREEN (MtSGR) gene. Transcript profiling analysis of the mutant showed that loss of the MtSGR function affected the expression of a large number of genes involved in different biological processes. Further analyses revealed that SGR is implicated in nodule development and senescence. MtSGR expression was detected across all nodule developmental zones and was higher in the senescence zone. The number of young nodules on the mutant roots was higher than in the wild type. Expression levels of several nodule senescence markers were reduced in the sgr mutant. Based on the MtSGR sequence, an alfalfa SGR gene (MsSGR) was cloned, and transgenic alfalfa lines were produced by RNA interference. Silencing of MsSGR led to the production of stay-green transgenic alfalfa. This beneficial trait offers the opportunity to produce premium alfalfa hay with a more greenish appearance. In addition, most of the transgenic alfalfa lines retained more than 50% of chlorophylls during senescence and had increased crude protein content. This study illustrates the effective use of knowledge gained from a model system for the genetic improvement of an important commercial crop.

An integrated analysis of miRNA and mRNA expressions in non-small cell lung cancers

Using DNA microarrays, we generated both mRNA and miRNA expression data from 6 non-small cell lung cancer (NSCLC) tissues and their matching normal control from adjacent tissues to identify potential miRNA markers for diagnostics. We demonstrated that hsa-miR-96 is significantly and consistently up-regulated in all 6 NSCLCs. We validated this result in an independent set of 35 paired tumors and their adjacent normal tissues, as well as their sera that are collected before surgical resection or chemotherapy, and the results suggested that hsa-miR-96 may play an important role in NSCLC development and has great potential to be used as a noninvasive marker for diagnosing NSCLC. We predicted potential miRNA target mRNAs based on different methods (TargetScan and miRanda). Further classification of miRNA regulated genes based on their relationship with miRNAs revealed that hsa-miR-96 and certain other miRNAs tend to down-regulate their target mRNAs in NSCLC development, which have expression levels permissive to direct interaction between miRNAs and their target mRNAs. In addition, we identified a significant correlation of miRNA regulation with genes coincide with high density of CpG islands, which suggests that miRNA may represent a primary regulatory mechanism governing basic cellular functions and cell differentiations, and such mechanism may be complementary to DNA methylation in repressing or activating gene expression.

Quantitative proteomics reveals cellular targets of celastrol

Celastrol, a natural substance isolated from plant extracts used in traditional Chinese medicine, has been extensively investigated as a possible drug for treatment of cancer, autoimmune diseases, and protein misfolding disorders. Although studies focusing on celastrol's effects in specific cellular pathways have revealed a considerable number of targets in a diverse array of in vitro models there is an essential need for investigations that can provide a global view of its effects. To assess cellular effects of celastrol and to identify target proteins as biomarkers for monitoring treatment regimes, we performed large-scale quantitative proteomics in cultured human lymphoblastoid cells, a cell type that can be readily prepared from human blood samples. Celastrol substantially modified the proteome composition and 158 of the close to 1800 proteins with robust quantitation showed at least a 1.5 fold change in protein levels. Up-regulated proteins play key roles in cytoprotection with a prominent group involved in quality control and processing of proteins traversing the endoplasmic reticulum. Increased levels of proteins essential for the cellular protection against oxidative stress including heme oxygenase 1, several peroxiredoxins and thioredoxins as well as proteins involved in the control of iron homeostasis were also observed. Specific analysis of the mitochondrial proteome strongly indicated that the mitochondrial association of certain antioxidant defense and apoptosis-regulating proteins increased in cells exposed to celastrol. Analysis of selected mRNA transcripts showed that celastrol activated several different stress response pathways and dose response studies furthermore showed that continuous exposure to sub-micromolar concentrations of celastrol is associated with reduced cellular viability and proliferation. The extensive catalog of regulated proteins presented here identifies numerous cellular effects of celastrol and constitutes a valuable biomarker tool for the development and monitoration of disease treatment strategies.

De Novo assembly of the Manila clam Ruditapes philippinarum transcriptome provides new insights into expression bias, mitochondrial doubly uniparental inheritance and sex determination

Males and females share the same genome, thus, phenotypic divergence requires differential gene expression and sex-specific regulation. Accordingly, the analysis of expression patterns is pivotal to the understanding of sex determination mechanisms. Many bivalves are stable gonochoric species, but the mechanism of gonad sexualization and the genes involved are still unknown. Moreover, during the period of sexual rest, a gonad is not present and sex cannot be determined. A mechanism associated with germ line differentiation in some bivalves, including the Manila clam Ruditapes philippinarum, is the doubly uniparental inheritance (DUI) of mitochondria, a variation of strict maternal inheritance. Two mitochondrial lineages are present, one transmitted through eggs and the other through sperm, as well as a mother-dependent sex bias of the progeny. We produced a de novo annotation of 17,186 transcripts from R. philippinarum and compared the transcriptomes of males and females and identified 1,575 genes with strong sex-specific expression and 166 sex-specific single nucleotide polymorphisms, obtaining preliminary information about genes that could be involved in sex determination. Then we compared the transcriptomes between a family producing predominantly females and a family producing predominantly males to identify candidate genes involved in regulation of sex-specific aspects of DUI system, finding a relationship between sex bias and differential expression of several ubiquitination genes. In mammalian embryos, sperm mitochondria are degraded by ubiquitination. A modification of this mechanism is hypothesized to be responsible for the retention of sperm mitochondria in male embryos of DUI species. Ubiquitination can additionally regulate gene expression, playing a role in sex determination of several animals. These data enable us to develop a model that incorporates both the DUI literature and our new findings.

RNA-Seq analysis and de novo transcriptome assembly of Hevea brasiliensis

Hevea brasiliensis, being the only source of commercial natural rubber, is an extremely economically important crop. In an effort to facilitate biological, biochemical and molecular research in rubber biosynthesis, here we report the use of next-generation massively parallel sequencing technologies and de novo transcriptome assembly to gain a comprehensive overview of the H. brasiliensis transcriptome. The sequencing output generated more than 12 million reads with an average length of 90 nt. In total 48,768 unigenes (mean size = 436 bp, median size = 328 bp) were assembled through de novo transcriptome assembly. Out of 13,807 H. brasiliensis cDNA sequences deposited in Genbank of the National Center for Biotechnology Information (NCBI) (as of Feb 2011), 11,746 sequences (84.5%) could be matched with the assembled unigenes through nucleotide BLAST. The assembled sequences were annotated with gene descriptions, Gene Ontology (GO) and Clusters of Orthologous Group (COG) terms. In all, 37,432 unigenes were successfully annotated, of which 24,545 (65.5%) aligned to Ricinus communis proteins. Furthermore, the annotated uingenes were functionally classified according to the GO, COG and Kyoto Encyclopedia of Genes and Genomes databases. Our data provides the most comprehensive sequence resource available for the study of rubber trees as well as demonstrates effective use of Illumina sequencing and de novo transcriptome assembly in a species lacking genomic information.

Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers

BACKGROUND

Sesame is an important oil crop, but limited transcriptomic and genomic data are currently available. This information is essential to clarify the fatty acid and lignan biosynthesis molecular mechanism. In addition, a shortage of sesame molecular markers limits the efficiency and accuracy of genetic breeding. High-throughput transcriptomic sequencing is essential to generate a large transcriptome sequence dataset for gene discovery and molecular marker development.

RESULTS

Sesame transcriptomes from five tissues were sequenced using Illumina paired-end sequencing technology. The cleaned raw reads were assembled into a total of 86,222 unigenes with an average length of 629 bp. Of the unigenes, 46,584 (54.03%) had significant similarity with proteins in the NCBI nonredundant protein database and Swiss-Prot database (E-value < 10-5). Of these annotated unigenes, 10,805 and 27,588 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. In total, 22,003 (25.52%) unigenes were mapped onto 119 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Furthermore, 44,750 unigenes showed homology to 15,460 Arabidopsis genes based on BLASTx analysis against The Arabidopsis Information Resource (TAIR, Version 10) and revealed relatively high gene coverage. In total, 7,702 unigenes were converted into SSR markers (EST-SSR). Dinucleotide SSRs were the dominant repeat motif (67.07%, 5,166), followed by trinucleotide (24.89%, 1,917), tetranucleotide (4.31%, 332), hexanucleotide (2.62%, 202), and pentanucleotide (1.10%, 85) SSRs. AG/CT (46.29%) was the dominant repeat motif, followed by AC/GT (16.07%), AT/AT (10.53%), AAG/CTT (6.23%), and AGG/CCT (3.39%). Fifty EST-SSRs were randomly selected to validate amplification and to determine the degree of polymorphism in the genomic DNA pools. Forty primer pairs successfully amplified DNA fragments and detected significant amounts of polymorphism among 24 sesame accessions.

CONCLUSIONS

This study demonstrates that Illumina paired-end sequencing is a fast and cost-effective approach to gene discovery and molecular marker development in non-model organisms. Our results provide a comprehensive sequence resource for sesame research.

Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers

Background

Sesame is an important oil crop, but limited transcriptomic and genomic data are currently available. This information is essential to clarify the fatty acid and lignan biosynthesis molecular mechanism. In addition, a shortage of sesame molecular markers limits the efficiency and accuracy of genetic breeding. High-throughput transcriptomic sequencing is essential to generate a large transcriptome sequence dataset for gene discovery and molecular marker development.

Results

Sesame transcriptomes from five tissues were sequenced using Illumina paired-end sequencing technology. The cleaned raw reads were assembled into a total of 86,222 unigenes with an average length of 629 bp. Of the unigenes, 46,584 (54.03%) had significant similarity with proteins in the NCBI nonredundant protein database and Swiss-Prot database (E-value < 10^-5). Of these annotated unigenes, 10,805 and 27,588 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. In total, 22,003 (25.52%) unigenes were mapped onto 119 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Furthermore, 44,750 unigenes showed homology to 15,460 Arabidopsis genes based on BLASTx analysis against The Arabidopsis Information Resource (TAIR, Version 10) and revealed relatively high gene coverage. In total, 7,702 unigenes were converted into SSR markers (EST-SSR). Dinucleotide SSRs were the dominant repeat motif (67.07%, 5,166), followed by trinucleotide (24.89%, 1,917), tetranucleotide (4.31%, 332), hexanucleotide (2.62%, 202), and pentanucleotide (1.10%, 85) SSRs. AG/CT (46.29%) was the dominant repeat motif, followed by AC/GT (16.07%), AT/AT (10.53%), AAG/CTT (6.23%), and AGG/CCT (3.39%). Fifty EST-SSRs were randomly selected to validate amplification and to determine the degree of polymorphism in the genomic DNA pools. Forty primer pairs successfully amplified DNA fragments and detected significant amounts of polymorphism among 24 sesame accessions.

Conclusions

This study demonstrates that Illumina paired-end sequencing is a fast and cost-effective approach to gene discovery and molecular marker development in non-model organisms. Our results provide a comprehensive sequence resource for sesame research.

The transcriptome of Trichuris suis--first molecular insights into a parasite with curative properties for key immune diseases of humans

BACKGROUND

Iatrogenic infection of humans with Trichuris suis (a parasitic nematode of swine) is being evaluated or promoted as a biological, curative treatment of immune diseases, such as inflammatory bowel disease (IBD) and ulcerative colitis, in humans. Although it is understood that short-term T. suis infection in people with such diseases usually induces a modified Th2-immune response, nothing is known about the molecules in the parasite that induce this response.

METHODOLOGY/PRINCIPAL FINDINGS

As a first step toward filling the gaps in our knowledge of the molecular biology of T. suis, we characterised the transcriptome of the adult stage of this nematode employing next-generation sequencing and bioinformatic techniques. A total of ∼65,000,000 reads were generated and assembled into ∼20,000 contiguous sequences ( = contigs); ∼17,000 peptides were predicted and classified based on homology searches, protein motifs and gene ontology and biological pathway mapping.

CONCLUSIONS

These analyses provided interesting insights into a number of molecular groups, particularly predicted excreted/secreted molecules (n = 1,288), likely to be involved in the parasite-host interactions, and also various molecules (n = 120) linked to chemokine, T-cell receptor and TGF-β signalling as well as leukocyte transendothelial migration and natural killer cell-mediated cytotoxicity, which are likely to be immuno-regulatory or -modulatory in the infected host. This information provides a conceptual framework within which to test the immunobiological basis for the curative effect of T. suis infection in humans against some immune diseases. Importantly, the T. suis transcriptome characterised herein provides a curated resource for detailed studies of the immuno-molecular biology of this parasite, and will underpin future genomic and proteomic explorations.

Conserved and divergent rhythms of crassulacean acid metabolism-related and core clock gene expression in the cactus Opuntia ficus-indica

The cactus Opuntia ficus-indica is a constitutive Crassulacean acid metabolism (CAM) species. Current knowledge of CAM metabolism suggests that the enzyme phosphoenolpyruvate carboxylase kinase (PPCK) is circadian regulated at the transcriptional level, whereas phosphoenolpyruvate carboxylase (PEPC), malate dehydrogenase (MDH), NADP-malic enzyme (NADP-ME), and pyruvate phosphate dikinase (PPDK) are posttranslationally controlled. As little transcriptomic data are available from obligate CAM plants, we created an expressed sequence tag database derived from different organs and developmental stages. Sequences were assembled, compared with sequences in the National Center for Biotechnology Information nonredundant database for identification of putative orthologs, and mapped using Kyoto Encyclopedia of Genes and Genomes Orthology and Gene Ontology. We identified genes involved in circadian regulation and CAM metabolism for transcriptomic analysis in plants grown in long days. We identified stable reference genes for quantitative polymerase chain reaction and found that OfiSAND, like its counterpart in Arabidopsis (Arabidopsis thaliana), and OfiTUB are generally appropriate standards for use in the quantification of gene expression in O. ficus-indica. Three kinds of expression profiles were found: transcripts of OfiPPCK oscillated with a 24-h periodicity; transcripts of the light-active OfiNADP-ME and OfiPPDK genes adapted to 12-h cycles, while transcript accumulation patterns of OfiPEPC and OfiMDH were arrhythmic. Expression of the circadian clock gene OfiTOC1, similar to Arabidopsis, oscillated with a 24-h periodicity, peaking at night. Expression of OfiCCA1 and OfiPRR9, unlike in Arabidopsis, adapted best to a 12-h rhythm, suggesting that circadian clock gene interactions differ from those of Arabidopsis. Our results indicate that the evolution of CAM metabolism could be the result of modified circadian regulation at both the transcriptional and posttranscriptional levels.

Sequencing and validation of the genome of a Campylobacter concisus reveals intra-species diversity

Campylobacter concisus is an emerging pathogen of the human gastrointestinal tract. Its role in different diseases remains a subject of debate; this may be due to strain to strain genetic variation. Here, we sequence and analyze the genome of a C. concisus from a biopsy of a child with Crohn's disease (UNSWCD); the second such genome for this species. A 1.8 Mb genome was assembled with paired-end reads from a next-generation sequencer. This genome is smaller than the 2.1 Mb C. concisus reference BAA-1457. While 1593 genes were conserved across UNSWCD and BAA-1457, 138 genes from UNSWCD and 281 from BAA-1457 were unique when compared against the other. To further validate the genome assembly and annotation, comprehensive shotgun proteomics was performed. This confirmed 78% of open reading frames in UNSWCD and, importantly, provided evidence of expression for 217 proteins previously defined as 'hypothetical' in Campylobacter. Substantial functional differences were observed between the UNSWCD and the reference strain. Enrichment analysis revealed differences in membrane proteins, response to stimulus, molecular transport and electron carriers. Synteny maps for the 281 genes not present in UNSWCD identified seven functionally associated gene clusters. These included one associated with the CRISPR family and another which encoded multiple restriction endonucleases; these genes are all involved in resistance to phage attack. Many of the observed differences are consistent with UNSWCD having adapted to greater surface interaction with host cells, as opposed to BAA-1457 which may prefer a free-living environment.

Proteomic analysis of peritrophic membrane (PM) from the midgut of fifth-instar larvae, Bombyx mori

The insect peritrophic membrane (PM), separating midgut epithelium and intestinal contents, is protective lining for the epithelium and plays the important role in absorption of nutrients, and also is the first barrier to the pathogens ingested through oral feeding. In order to understand the biological function of silkworm larval PM, shotgun liquid chromatography tandem mass spectrometry (LC-MS/MS) approach was applied to investigate its protein composition. Total 47 proteins were identified, of which 51.1% of the proteins had the isoelectric point (pI) within the range of 5-7, and 53.2% had molecular weights within the range 15-45 kDa. Most of them were found to be closely related to larval nutrients metabolism and innate immunity. Furthermore, these identified proteins were annotated according to Gene Ontology Annotation in terms of molecular function, biological process and cell localization. Most of the proteins had catalytic activity, binding activity and transport function. The knowledge obtained from this study will favour us to well understand the role of larval PM in larval physiological activities, and also help us to find the potential target and design better biopesticides to control pest, particularly the Lepidoptera insect.

Genome-wide binding site analysis of FAR-RED ELONGATED HYPOCOTYL3 reveals its novel function in Arabidopsis development

FAR-RED ELONGATED HYPOCOTYL3 (FHY3) and its homolog FAR-RED IMPAIRED RESPONSE1 (FAR1), two transposase-derived transcription factors, are key components in phytochrome A signaling and the circadian clock. Here, we use chromatin immunoprecipitation-based sequencing (ChIP-seq) to identify 1559 and 1009 FHY3 direct target genes in darkness (D) and far-red (FR) light conditions, respectively, in the Arabidopsis thaliana genome. FHY3 preferentially binds to promoters through the FHY3/FAR1 binding motif (CACGCGC). Interestingly, FHY3 also binds to two motifs in the 178-bp Arabidopsis centromeric repeats. Comparison between the ChIP-seq and microarray data indicates that FHY3 quickly regulates the expression of 197 and 86 genes in D and FR, respectively. FHY3 also coregulates a number of common target genes with PHYTOCHROME INTERACTING FACTOR 3-LIKE5 and ELONGATED HYPOCOTYL5. Moreover, we uncover a role for FHY3 in controlling chloroplast development by directly activating the expression of ACCUMULATION AND REPLICATION OF CHLOROPLASTS5, whose product is a structural component of the latter stages of chloroplast division in Arabidopsis. Taken together, our data suggest that FHY3 regulates multiple facets of plant development, thus providing insights into its functions beyond light and circadian pathways.

Identification of genes directly involved in shell formation and their functions in pearl oyster, Pinctada fucata

Mollusk shell formation is a fascinating aspect of biomineralization research. Shell matrix proteins play crucial roles in the control of calcium carbonate crystallization during shell formation in the pearl oyster, Pinctada fucata. Characterization of biomineralization-related genes during larval development could enhance our understanding of shell formation. Genes involved in shell biomineralization were isolated by constructing three suppression subtractive hybridization (SSH) libraries that represented genes expressed at key points during larval shell formation. A total of 2,923 ESTs from these libraries were sequenced and gave 990 unigenes. Unigenes coding for secreted proteins and proteins with tandem-arranged repeat units were screened in the three SSH libraries. A set of sequences coding for genes involved in shell formation was obtained. RT-PCR and in situ hybridization assays were carried out on five genes to investigate their spatial expression in several tissues, especially the mantle tissue. They all showed a different expression pattern from known biomineralization-related genes. Inhibition of the five genes by RNA interference resulted in different defects of the nacreous layer, indicating that they all were involved in aragonite crystallization. Intriguingly, one gene (UD_Cluster94.seq.Singlet1) was restricted to the ‘aragonitic line’. The current data has yielded for the first time, to our knowledge, a suite of biomineralization-related genes active during the developmental stages of P.fucata, five of which were responsible for nacreous layer formation. This provides a useful starting point for isolating new genes involved in shell formation. The effects of genes on the formation of the ‘aragonitic line’, and other areas of the nacreous layer, suggests a different control mechanism for aragonite crystallization initiation from that of mature aragonite growth.

The first insight into the tissue specific taxus transcriptome via Illumina second generation sequencing

Background

Illumina second generation sequencing is now an efficient route for generating enormous sequence collections that represent expressed genes and quantitate expression level. Taxus is a world-wide endangered gymnosperm genus and forms an important anti-cancer medicinal resource, but the large and complex genomes of Taxus have hindered the development of genomic resources. The research of its tissue-specific transcriptome is absent. There is also no study concerning the association between the plant transcriptome and metabolome with respect to the plant tissue type.

Methodology/Principal Findings

We performed the de novo assembly of Taxus mairei transcriptome using Illumina paired-end sequencing technology. In a single run, we produced 13,737,528 sequencing reads corresponding to 2.03 Gb total nucleotides. These reads were assembled into 36,493 unique sequences. Based on similarity search with known proteins, 23,515 Unigenes were identified to have the Blast hit with a cut-off E-value above 10⁻⁵. Furthermore, we investigated the transcriptome difference of three Taxus tissues using a tag-based digital gene expression system. We obtained a sequencing depth of over 3.15 million tags per sample and identified a large number of genes associated with tissue specific functions and taxane biosynthetic pathway. The expression of the taxane biosynthetic genes is significantly higher in the root than in the leaf and the stem, while high activity of taxane-producing pathway in the root was also revealed via metabolomic analyses. Moreover, many antisense transcripts and novel transcripts were found; clusters with similar differential expression patterns, enriched GO terms and enriched metabolic pathways with regard to the differentially expressed genes were revealed for the first time.

Conclusions/Significance

Our data provides the most comprehensive sequence resource available for Taxus study and will help define mechanisms of tissue specific functions and secondary metabolism in non-model plant organisms.

Differential proteomic response of rice (Oryza sativa) leaves exposed to high- and low-temperature stress

Global mean surface temperature has been predicted to increase by 1.8-4°C within this century, accompanied by an increase in the magnitude and frequency of extreme temperature events. Developing rice cultivars better adapted to non-optimal temperatures is essential to increase rice yield in the future and, hence, understanding the molecular response of rice to temperature stress is necessary. In this study, we investigated the proteomic responses of leaves of 24-day-old rice seedlings to sudden temperature changes. Rice seedlings grown at 28/20°C (day/night) were subjected to 3-day exposure to 12/5°C or 20/12°C (day/night) for low-temperature stress, and 36/28°C or 44/36°C (day/night) for high-temperature stress, followed by quantitative label-free shotgun proteomic analysis on biological triplicates of each treatment. Out of over 1100 proteins identified in one or more temperature treatments, more than 400 were found to be responsive to temperature stress. Of these, 43, 126 and 47 proteins were exclusively found at 12/5, 20/12 and 44/36°C (day/night), respectively. Our results showed that a greater change occurs in the rice leaf proteome at 20/12°C (day/night) in comparison to other non-optimal temperature regimes. In addition, our study identified more than 20 novel stress-response proteins.

Protein profile of Bacillus subtilis spore

Natural wild-type strains of Bacillus subtilis spore is regarded as a non-pathogenic for both human and animal, and has been classified as a novel food which is currently being used as probiotics added in the consumption. To identify B. subtilis spore proteins, we have accomplished a preliminary proteomic analysis of B. subtilis spore, with a combination of two-dimensional electrophoretic separations and matrix-assisted laser desorption ionization tandem time of flight mass spectrometry (MALDI-TOF-MS). In this article, we presented a reference map of 158 B. subtilis spore proteins with an isoelectric point (pI) between 4 and 7. Followed by mass spectrometry (MS) analysis, we identified 71 B. subtilis spore proteins with high level of confidence. Database searches, combined with hydropathy analysis and GO analysis revealed that most of the B. subtilis spore proteins were hydrophilic proteins related to catalytic function. These results should accelerate efforts to understand the resistance of spore to harsh conditions.

Predicting housekeeping genes based on Fourier analysis

Housekeeping genes (HKGs) generally have fundamental functions in basic biochemical processes in organisms, and usually have relatively steady expression levels across various tissues. They play an important role in the normalization of microarray technology. Using Fourier analysis we transformed gene expression time-series from a Hela cell cycle gene expression dataset into Fourier spectra, and designed an effective computational method for discriminating between HKGs and non-HKGs using the support vector machine (SVM) supervised learning algorithm which can extract significant features of the spectra, providing a basis for identifying specific gene expression patterns. Using our method we identified 510 human HKGs, and then validated them by comparison with two independent sets of tissue expression profiles. Results showed that our predicted HKG set is more reliable than three previously identified sets of HKGs.

Benefits of random-priming: exhaustive survey of a cDNA library from lung tissue of a SARS patient

The severe acute respiratory syndrome (SARS) leads to severe injury in the lungs with multiple factors, though the pathogenesis is still largely unclear. This paper describes the particular analyses of the transcriptome of human lung tissue that was infected by SARS‐associated coronavirus (SARS‐CoV). Random primers were used to produce ESTs from total RNA samples of the lung tissue. The result showed a high diversity of the transcripts, covering much of the human genome, including loci which do not contain protein coding sequences. 10,801 ESTs were generated and assembled into 267 contigs plus 7,659 singletons. Sequences matching to SARS‐CoV RNAs and other pneumonia‐related microbes were found. The transcripts were well classified by functional annotation. Among the 7,872 assembled sequences that were identified as from human genome, 578 non‐coding genes were revealed by BLAST search. The transcripts were mapped to the human genome with the restriction of identity = 100%, which found a candidate pool of 448 novel transcriptional loci where EST transcriptional signal was never found before. Among these, 13 loci were never reported to be transcriptional by other detection methods such as gene chips, tiling arrays, and paired‐end ditags (PETs). The result showed that random‐priming cDNA library is valid for the investigation of transcript diversity in the virus‐infected tissue. The EST data could be a useful supplemental source for SARS pathology researches. J. Med. Virol. 83:574–586, 2011. © 2011 Wiley‐Liss, Inc.

Elementary research of the formation mechanism of sex-related fluorescent cocoon of silkworm, Bombyx mori

To understand mechanisms for the difference of uptaking and transporting the pigments between the male and female in the silkworm, Bombyx mori strain of sex-related fluorescent cocoon, the fluorescent pigments in the midgut lumen, midgut, blood, silk glands and cocoon were analyzed with thin-layer chromatography, and showed that fluorescent colors of cocoons consisted with that of blood and silk glands. The different fluorescent colors of cocoons between the male and female may be mainly caused by the difference of accumulation and transportation for fluorescent pigments in the midgut and in the silk glands. Furthermore the midgut proteins were separated with Native-PAGE, and the proteins respectively recovered from three fluorescent regions presenting on a Native-PAGE gel for the female silkworms were determined using shotgun proteomics and mass spectrometry sequencing, of which 60, 40 and 18 proteins respectively from the region 1, 2 and 3 were identified. It was found that the several kinds of low molecular mass 30 kDa lipoproteins and the actins could be detected in all three regions, troponin, 30 kDa lipoprotein and 27 kDa glycoprotein precursor could be detected in the region 2 and 3, suggesting these proteins may be fluorescent pigments binding candidates proteins. Analysis of gene ontology indicated that the identified proteins in the three regions linked to the cellular component, molecular function, and biological process categories. These results provide a new clew to understand the formation mechanism of sex-related fluorescent cocoon of silkworm.

Pseudogene-derived small interference RNAs regulate gene expression in African Trypanosoma brucei

Pseudogenes have been shown to acquire unique regulatory roles from more and more organisms. We report the observation of a cluster of siRNAs derived from pseudogenes of African Trypanosoma brucei using high through-put analysis. We show that these pseudogene-derived siRNAs suppress gene expression through RNA interference. The discovery that siRNAs may originate from pseudogenes and regulate gene expression in a unicellular eukaryote provides insights into the functional roles of pseudogenes and into the origin of noncoding small RNAs.

TMPad: an integrated structural database for helix-packing folds in transmembrane proteins

α-Helical transmembrane (TM) proteins play an important role in many critical and diverse biological processes, and specific associations between TM helices are important determinants for membrane protein folding, dynamics and function. In order to gain insights into the above phenomena, it is necessary to investigate different types of helix-packing modes and interactions. However, such information is difficult to obtain because of the experimental impediment and a lack of a well-annotated source of helix-packing folds in TM proteins. We have developed the TMPad (TransMembrane Protein Helix-Packing Database) which addresses the above issues by integrating experimentally observed helix–helix interactions and related structural information of membrane proteins. Specifically, the TMPad offers pre-calculated geometric descriptors at the helix-packing interface including residue backbone/side-chain contacts, interhelical distances and crossing angles, helical translational shifts and rotational angles. The TMPad also includes the corresponding sequence, topology, lipid accessibility, ligand-binding information and supports structural classification, schematic diagrams and visualization of the above structural features of TM helix-packing. Through detailed annotations and visualizations of helix-packing, this online resource can serve as an information gateway for deciphering the relationship between helix–helix interactions and higher levels of organization in TM protein structure and function. The website of the TMPad is freely accessible to the public at http://bio-cluster.iis.sinica.edu.tw/TMPad.

Global expression analysis revealed novel gender-specific gene expression features in the blood fluke parasite Schistosoma japonicum

Background

Schistosoma japonicum is one of the remarkable Platyhelminths that are endemic in China and Southeast Asian countries. The parasite is dioecious and can reside inside the host for many years. Rapid reproduction by producing large number of eggs and count-react host anti-parasite responses are the strategies that benefit long term survival of the parasite. Praziquantel is currently the only drug that is effective against the worms. Development of novel antiparasite reagents and immune-prevention measures rely on the deciphering of parasite biology. The decoding of the genomic sequence of the parasite has made it possible to dissect the functions of genes that govern the development of the parasite. In this study, the polyadenylated transcripts from male and female S. japonicum were isolated for deep sequencing and the sequences were systematically analysed.

Results

First, the number of genes actively expressed in the two sexes of S. japonicum was similar, but around 50% of genes were biased to either male or female in expression. Secondly, it was, at the first time, found that more than 50% of the coding region of the genome was transcribed from both strands. Among them, 65% of the genes had sense and their cognate antisense transcripts co-expressed, whereas 35% had inverse relationship between sense and antisense transcript abundance. Further, based on gene ontological analysis, more than 2,000 genes were functionally categorized and biological pathways that are differentially functional in male or female parasites were elucidated.

Conclusions

Male and female schistosomal parasites differ in gene expression patterns, many metabolic and biological pathways have been identified in this study and genes differentially expressed in gender specific manner were presented. Importantly, more than 50% of the coding regions of the S. japonicum genome transcribed from both strands, antisense RNA-mediated gene regulation might play a critical role in the parasite biology.

De novo transcriptome sequencing in Salvia miltiorrhiza to identify genes involved in the biosynthesis of active ingredients

Medicinal Salvia miltiorrhiza is a Chinese herb commonly used for treating cardiovascular diseases and neuroasthenic insomnia. However, little is known at the genetics level about how its compounds are synthesized in that plant. Here, we obtained 56,774 unigenes (average length = 467 bases) in its transcriptome by performing Solexa deep sequencing over the entire growing cycle. Unigenes (34,340; 60.49%) were annotated and 2545 unigenes were assigned to specific pathways. Unigenes (1539) were identified as part of five major, secondary-metabolite pathways, covering almost all nodes in the phenylpropanoid and terpenoid pathways. Using Blast search against AGRIS, 1341 unigenes were found homologous to 686 Arabidopsis transcription factor genes. Real-time PCR was also used to verify the spatio-temporal expression patterns of several novel transcripts related to biosynthesis of active ingredients in that species. These results not only enrich the gene resource but also benefit research into its molecular genetics and functional genomics.

Quantitative proteomic analysis of cold-responsive proteins in rice

Rice is susceptible to cold stress and with a future of climatic instability we will be unable to produce enough rice to satisfy increasing demand. A thorough understanding of the molecular responses to thermal stress is imperative for engineering cultivars, which have greater resistance to low temperature stress. In this study we investigated the proteomic response of rice seedlings to 48, 72 and 96 h of cold stress at 12-14°C. The use of both label-free and iTRAQ approaches in the analysis of global protein expression enabled us to assess the complementarity of the two techniques for use in plant proteomics. The approaches yielded a similar biological response to cold stress despite a disparity in proteins identified. The label-free approach identified 236 cold-responsive proteins compared to 85 in iTRAQ results, with only 24 proteins in common. Functional analysis revealed differential expression of proteins involved in transport, photosynthesis, generation of precursor metabolites and energy; and, more specifically, histones and vitamin B biosynthetic proteins were observed to be affected by cold stress.

Interaction of the tobacco lectin with histone proteins

The tobacco (Nicotiana tabacum) agglutinin or Nictaba is a member of a novel class of plant lectins residing in the nucleus and the cytoplasm of tobacco cells. Since tobacco lectin expression is only observed after the plant has been subjected to stress situations such as jasmonate treatment or insect attack, Nictaba is believed to act as a signaling protein involved in the stress physiology of the plant. In this paper, a nuclear proteomics approach was followed to identify the binding partners for Nictaba in the nucleus and the cytoplasm of tobacco cv Xanthi cells. Using lectin affinity chromatography and pull-down assays, it was shown that Nictaba interacts primarily with histone proteins. Binding of Nictaba with histone H2B was confirmed in vitro using affinity chromatography of purified calf thymus histone proteins on a Nictaba column. Elution of Nictaba-interacting histone proteins was achieved with 1 m N-acetylglucosamine (GlcNAc). Moreover, mass spectrometry analyses indicated that the Nictaba-interacting histone proteins are modified by O-GlcNAc. Since the lectin-histone interaction was shown to be carbohydrate dependent, it is proposed that Nictaba might fulfill a signaling role in response to stress by interacting with O-GlcNAcylated proteins in the plant cell nucleus.

Diversity in protein glycosylation among insect species

BACKGROUND

A very common protein modification in multicellular organisms is protein glycosylation or the addition of carbohydrate structures to the peptide backbone. Although the Class of the Insecta is the largest animal taxon on Earth, almost all information concerning glycosylation in insects is derived from studies with only one species, namely the fruit fly Drosophila melanogaster.

METHODOLOGY/PRINCIPAL FINDINGS

In this report, the differences in glycoproteomes between insects belonging to several economically important insect orders were studied. Using GNA (Galanthus nivalis agglutinin) affinity chromatography, different sets of glycoproteins with mannosyl-containing glycan structures were purified from the flour beetle (Tribolium castaneum), the silkworm (Bombyx mori), the honeybee (Apis mellifera), the fruit fly (D. melanogaster) and the pea aphid (Acyrthosiphon pisum). To identify and characterize the purified glycoproteins, LC-MS/MS analysis was performed. For all insect species, it was demonstrated that glycoproteins were related to a broad range of biological processes and molecular functions. Moreover, the majority of glycoproteins retained on the GNA column were unique to one particular insect species and only a few glycoproteins were present in the five different glycoprotein sets. Furthermore, these data support the hypothesis that insect glycoproteins can be decorated with mannosylated O-glycans.

CONCLUSIONS/SIGNIFICANCE

The results presented here demonstrate that oligomannose N-glycosylation events are highly specific depending on the insect species. In addition, we also demonstrated that protein O-mannosylation in insect species may occur more frequently than currently believed.

Distinct gene expression profiles in egg and synergid cells of rice as revealed by cell type-specific microarrays

Double fertilization in flowering plants refers to a process in which two sperm cells, carried by the pollen tube, fertilize both the egg and the central cell after their release into a synergid cell of the female gametophyte. The molecular processes by which the female gametophytic cells express their unique functions during fertilization are not well understood. Genes expressed in egg and synergid cells might be important for multiple stages of the plant reproductive process. Here, we profiled genome-wide gene expression in egg and synergid cells in rice (Oryza sativa), a model monocot, using a nonenzymatic cell isolation technique. We found that the expression profiles of the egg and synergid cells were already specified at the micropylar end of the female gametophyte during the short developmental period that comprises the three consecutive mitotic nuclear divisions after megaspore generation. In addition, we identified a large number of genes expressed in the rice egg and synergid cells and characterized these genes using Gene Ontology analysis. The analysis suggested that epigenetic and posttranscriptional regulatory mechanisms are involved in the specification and/or maintenance of these cells. Comparisons between the rice profiles and reported Arabidopsis (Arabidopsis thaliana) profiles revealed that genes enriched in the egg/synergid cell of rice were distinct from those in Arabidopsis.

Less label, more free: approaches in label-free quantitative mass spectrometry

In this review we examine techniques, software, and statistical analyses used in label-free quantitative proteomics studies for area under the curve and spectral counting approaches. Recent advances in the field are discussed in an order that reflects a logical workflow design. Examples of studies that follow this design are presented to highlight the requirement for statistical assessment and further experiments to validate results from label-free quantitation. Limitations of label-free approaches are considered, label-free approaches are compared with labelling techniques, and forward-looking applications for label-free quantitative data are presented. We conclude that label-free quantitative proteomics is a reliable, versatile, and cost-effective alternative to labelled quantitation.

Solexa sequencing analysis of chicken pre-adipocyte microRNAs

MicroRNAs (miRNAs) are small non-coding RNAs that play important roles in a variety of biological processes. Studies of miRNAs in mammals suggest that many are involved in lipid metabolism and adipocyte differentiation, but little is known about miRNA expression profiles during chicken adipogenesis. In this study, the Solexa sequencing approach was used to sequence a small RNA library prepared from Arbor Acres broiler pre-adipocytes, and more than 10⁶ short sequence reads were obtained. From these, 159 known chicken miRNAs and 63 novel miRNAs were identified using a bioinformatics approach. Fifty-nine of these miRNA genes were further organized into 27 compact miRNA genomic clusters, and 34 new chicken mirtrons were also discovered, among which there were 27 mirtron candidates. These findings should serve as a foundation for future research on the functional roles of miRNAs in chicken adipocyte differentiation.

De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweet potato (Ipomoea batatas)

Background

The tuberous root of sweetpotato is an important agricultural and biological organ. There are not sufficient transcriptomic and genomic data in public databases for understanding of the molecular mechanism underlying the tuberous root formation and development. Thus, high throughput transcriptome sequencing is needed to generate enormous transcript sequences from sweetpotato root for gene discovery and molecular marker development.

Results

In this study, more than 59 million sequencing reads were generated using Illumina paired-end sequencing technology. De novo assembly yielded 56,516 unigenes with an average length of 581 bp. Based on sequence similarity search with known proteins, a total of 35,051 (62.02%) genes were identified. Out of these annotated unigenes, 5,046 and 11,983 unigenes were assigned to gene ontology and clusters of orthologous group, respectively. Searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG) indicated that 17,598 (31.14%) unigenes were mapped to 124 KEGG pathways, and 11,056 were assigned to metabolic pathways, which were well represented by carbohydrate metabolism and biosynthesis of secondary metabolite. In addition, 4,114 cDNA SSRs (cSSRs) were identified as potential molecular markers in our unigenes. One hundred pairs of PCR primers were designed and used for validation of the amplification and assessment of the polymorphism in genomic DNA pools. The result revealed that 92 primer pairs were successfully amplified in initial screening tests.

Conclusion

This study generated a substantial fraction of sweetpotato transcript sequences, which can be used to discover novel genes associated with tuberous root formation and development and will also make it possible to construct high density microarrays for further characterization of gene expression profiles during these processes. Thousands of cSSR markers identified in the present study can enrich molecular markers and will facilitate marker-assisted selection in sweetpotato breeding. Overall, these sequences and markers will provide valuable resources for the sweetpotato community. Additionally, these results also suggested that transcriptome analysis based on Illumina paired-end sequencing is a powerful tool for gene discovery and molecular marker development for non-model species, especially those with large and complex genome.

Differential metabolic response of cultured rice (Oryza sativa) cells exposed to high- and low-temperature stress

Global mean temperatures are expected to rise by 2-4.5 degrees C by 2100, accompanied by an increase in frequency and amplitude of extreme temperature events. Greater climatic extremes and an expanded range of cultivation will expose rice to increasing stress in the future. Understanding gene expression in disparate thermal regimes is important for the engineering of cultivars with tolerance to nonoptimal temperatures. Our study investigated the proteomic responses of rice cell suspension cultures to sudden temperature changes. Cell cultures grown at 28 degrees C were subjected to 3-day exposure to 12 or 20 degrees C for low-temperature stress, and 36 or 44 degrees C for high-temperature stress. Quantitative label-free shotgun proteomic analysis was performed on biological triplicates of each treatment. Over 1900 proteins were expressed in one or more temperature treatments, and, of these, more than 850 were found to be responsive to either of the temperature extremes. These temperature-responsive proteins included more than 300 proteins which were uniquely expressed at either 12 or 44 degrees C. Our study also identified 40 novel stress-response proteins and observed that switching between the classical and the alternative pathways of sucrose metabolism occurs in response to extremes of temperature.

Massively parallel sequencing and analysis of expressed sequence tags in a successful invasive plant

BACKGROUND

Invasive species pose a significant threat to global economies, agriculture and biodiversity. Despite progress towards understanding the ecological factors associated with plant invasions, limited genomic resources have made it difficult to elucidate the evolutionary and genetic factors responsible for invasiveness. This study presents the first expressed sequence tag (EST) collection for Senecio madagascariensis, a globally invasive plant species.

METHODS

We used pyrosequencing of one normalized and two subtractive libraries, derived from one native and one invasive population, to generate an EST collection. ESTs were assembled into contigs, annotated by BLAST comparison with the NCBI non-redundant protein database and assigned gene ontology (GO) terms from the Plant GO Slim ontologies.

KEY RESULTS

Assembly of the 221,746 sequence reads resulted in 12,442 contigs. Over 50 % (6183) of 12,442 contigs showed significant homology to proteins in the NCBI database, representing approx. 4800 independent transcripts. The molecular transducer GO term was significantly over-represented in the native (South African) subtractive library compared with the invasive (Australian) library. Based on NCBI BLAST hits and literature searches, 40 % of the molecular transducer genes identified in the South African subtractive library are likely to be involved in response to biotic stimuli, such as fungal, bacterial and viral pathogens.

CONCLUSIONS

This EST collection is the first representation of the S. madagascariensis transcriptome and provides an important resource for the discovery of candidate genes associated with plant invasiveness. The over-representation of molecular transducer genes associated with defence responses in the native subtractive library provides preliminary support for aspects of the enemy release and evolution of increased competitive ability hypotheses in this successful invasive. This study highlights the contribution of next-generation sequencing to better understanding the molecular mechanisms underlying ecological hypotheses that are important in successful plant invasions.

Comparative transcriptional profiling and preliminary study on heterosis mechanism of super-hybrid rice

Heterosis is a biological phenomenon whereby the offspring from two parents show improved and superior performance than either inbred parental lines. Hybrid rice is one of the most successful apotheoses in crops utilizing heterosis. Transcriptional profiling of F(1) super-hybrid rice Liangyou-2186 and its parents by serial analysis of gene expression (SAGE) revealed 1183 differentially expressed genes (DGs), among which DGs were found significantly enriched in pathways such as photosynthesis and carbon-fixation, and most of the key genes involved in the carbon-fixation pathway exhibited up-regulated expression in F(1) hybrid rice. Moreover, increased catabolic activity of corresponding enzymes and photosynthetic efficiency were also detected, which combined to indicate that carbon fixation is enhanced in F(1) hybrid, and might probably be associated with the yield vigor and heterosis in super-hybrid rice. By correlating DGs with yield-related quantitative trait loci (QTL), a potential relationship between differential gene expression and phenotypic changes was also found. In addition, a regulatory network involving circadian-rhythms and light signaling pathways was also found, as previously reported in Arabidopsis, which suggest that such a network might also be related with heterosis in hybrid rice. Altogether, the present study provides another view for understanding the molecular mechanism underlying heterosis in rice.

Shotgun proteomic analysis of the fat body during metamorphosis of domesticated silkworm (Bombyx mori)

Protein expression profiles in the fat bodies of larval, pupal, and moth stages of silkworm were determined using shotgun proteomics and MS sequencing. We identified 138, 217, and 86 proteins from the larval, pupal and moth stages, respectively, of which 12 were shared by the 3 stages. There were 92, 150, and 45 specific proteins identified in the larval, pupal and moth stages, respectively, of which 17, 68, and 9 had functional annotations. Among the specific proteins identified in moth fat body, sex-specific storage-protein 1 precursor and chorion protein B8 were unique to the moth stage, indicating that the moth stage fat body is more important for adult sexual characteristics. Many ribosomal proteins (L23, L4, L5, P2, S10, S11, S15A and S3) were found in pupal fat bodies, whereas only three (L14, S20, and S7) and none were identified in larval and moth fat bodies, respectively. Twenty-three metabolic enzymes were identified in the pupal stage, while only four and two were identified in the larval and moth stages, respectively. In addition, an important protein, gloverin2, was only identified in larval fat bodies. Gene ontology (GO) analysis of the proteins specific to the three stages linked them to the cellular component, molecular function, and biological process categories. The most diverse GO functional classes were involved by the relatively less specific proteins identified in larva. GO analysis of the proteins shared among the three stages showed that the pupa and moth stages shared the most similar protein functions in the fat body.

Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune-relevant genes in marine fish

Background

Systematic research on fish immunogenetics is indispensable in understanding the origin and evolution of immune systems. This has long been a challenging task because of the limited number of deep sequencing technologies and genome backgrounds of non-model fish available. The newly developed Solexa/Illumina RNA-seq and Digital gene expression (DGE) are high-throughput sequencing approaches and are powerful tools for genomic studies at the transcriptome level. This study reports the transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus using RNA-seq and DGE in an attempt to gain insights into the immunogenetics of marine fish.

Results

RNA-seq analysis generated 169,950 non-redundant consensus sequences, among which 48,987 functional transcripts with complete or various length encoding regions were identified. More than 52% of these transcripts are possibly involved in approximately 219 known metabolic or signalling pathways, while 2,673 transcripts were associated with immune-relevant genes. In addition, approximately 8% of the transcripts appeared to be fish-specific genes that have never been described before. DGE analysis revealed that the host transcriptome profile of Vibrio harveyi-challenged L. japonicus is considerably altered, as indicated by the significant up- or down-regulation of 1,224 strong infection-responsive transcripts. Results indicated an overall conservation of the components and transcriptome alterations underlying innate and adaptive immunity in fish and other vertebrate models. Analysis suggested the acquisition of numerous fish-specific immune system components during early vertebrate evolution.

Conclusion

This study provided a global survey of host defence gene activities against bacterial challenge in a non-model marine fish. Results can contribute to the in-depth study of candidate genes in marine fish immunity, and help improve current understanding of host-pathogen interactions and evolutionary history of immunogenetics from fish to mammals.

From raw materials to validated system: the construction of a genomic library and microarray to interpret systemic perturbations in Northern bobwhite

The limited availability of genomic tools and data for nonmodel species impedes computational and systems biology approaches in nonmodel organisms. Here we describe the development, functional annotation, and utilization of genomic tools for the avian wildlife species Northern bobwhite (Colinus virginianus) to determine the molecular impacts of exposure to 2,6-dinitrotoluene (2,6-DNT), a field contaminant of military concern. Massively parallel pyrosequencing of a normalized multitissue library of Northern bobwhite cDNAs yielded 71,384 unique transcripts that were annotated with gene ontology (GO), pathway information, and protein domain analysis. Comparative genome analyses with model organisms revealed functional homologies in 8,825 unique Northern bobwhite genes that are orthologous to 48% of Gallus gallus protein-coding genes. Pathway analysis and GO enrichment of genes differentially expressed in livers of birds exposed for 60 days (d) to 10 and 60 mg/kg/d 2,6-DNT revealed several impacts validated by RT-qPCR including: prostaglandin pathway-mediated inflammation, increased expression of a heme synthesis pathway in response to anemia, and a shift in energy metabolism toward protein catabolism via inhibition of control points for glucose and lipid metabolic pathways, PCK1 and PPARGC1, respectively. This research effort provides the first comprehensive annotated gene library for Northern bobwhite. Transcript expression analysis provided insights into the metabolic perturbations underlying several observed toxicological phenotypes in a 2,6-DNT exposure case study. Furthermore, the systemic impact of dinitrotoluenes on liver function appears conserved across species as PPAR signaling is similarly affected in fathead minnow liver tissue after exposure to 2,4-DNT.