WEGO: a web tool for plotting GO annotations

Next-generation sequencing-based transcriptome analysis of Cryptolaemus montrouzieri under insecticide stress reveals resistance-relevant genes in ladybirds

As the most efficient natural enemy of mealybugs, the ladybird Cryptolaemus montrouzieri Mulsant plays an important role in integrated pest management. We report here a profiling analysis of C. montrouzieri under insecticide stress to gain a deeper view of insecticide resistance in ladybirds. For transcriptome sequencing, more than 26 million sequencing reads were produced. These reads were assembled into 38,369 non-redundant transcripts (mean size=453 nt). 23,248 transcripts were annotated with their gene description. Using a tag-based DGE (Digital gene expression) system, over 5.7 million tags were sequenced in both the insecticide stress group and the control group, and mapped to 38,369 transcripts. We obtained 993 genes that were significantly up- or down-regulated under insecticide stress in the ladybird transcriptome. These results can contribute to in-depth research into the molecular mechanisms of resistance and enhance our current understanding of the effects of insecticides on natural enemies.

Transcriptome and comparative gene expression analysis of Sogatella furcifera (Horváth) in response to southern rice black-streaked dwarf virus

Background

The white backed planthopper (WBPH), Sogatella furcifera (Horváth), causes great damage to many crops by direct feeding or transmitting plant viruses. Southern rice black-streaked dwarf virus (SRBSDV), transmitted by WBPH, has become a great threat to rice production in East Asia.

Methodology/Principal Findings

By de novo transcriptome assembling and massive parallel pyrosequencing, we constructed two transcriptomes of WBPH and profiled the alternation of gene expression in response to SRBSDV infection in transcriptional level. Over 25 million reads of high-quality DNA sequences and 81388 different unigenes were generated using Illumina technology from both viruliferous and non-viruliferous WBPH. WBPH has a very similar gene ontological distribution to other two closely related rice planthoppers, Nilaparvata lugens and Laodelphax striatellus. 7291 microsatellite loci were also predicted which could be useful for further evolutionary analysis. Furthermore, comparative analysis of the two transcriptomes generated from viruliferous and non-viruliferous WBPH provided a list of candidate transcripts that potentially were elicited as a response to viral infection. Pathway analyses of a subset of these transcripts indicated that SRBSDV infection may perturb primary metabolism and the ubiquitin-proteasome pathways. In addition, 5.5% (181 out of 3315) of the genes in cell cytoskeleton organization pathway showed obvious changes. Our data also demonstrated that SRBSDV infection activated the immunity regulatory systems of WBPH, such as RNA interference, autophagy and antimicrobial peptide production.

Conclusions/Significance

We employed massively parallel pyrosequencing to collect ESTs from viruliferous and non-viruliferous samples of WBPH. 81388 different unigenes have been obtained. We for the first time described the direct effects of a Reoviridae family plant virus on global gene expression profiles of its insect vector using high-throughput sequencing. Our study will provide a road map for future investigations of the fascinating interactions between Reoviridae viruses and their insect vectors, and provide new strategies for crop protection.

Digital gene expression analysis based on integrated de novo transcriptome assembly of sweet potato [Ipomoea batatas (L.) Lam]

BACKGROUND

Sweet potato (Ipomoea batatas L. [Lam.]) ranks among the top six most important food crops in the world. It is widely grown throughout the world with high and stable yield, strong adaptability, rich nutrient content, and multiple uses. However, little is known about the molecular biology of this important non-model organism due to lack of genomic resources. Hence, studies based on high-throughput sequencing technologies are needed to get a comprehensive and integrated genomic resource and better understanding of gene expression patterns in different tissues and at various developmental stages.

METHODOLOGY/PRINCIPAL FINDINGS

Illumina paired-end (PE) RNA-Sequencing was performed, and generated 48.7 million of 75 bp PE reads. These reads were de novo assembled into 128,052 transcripts (≥ 100 bp), which correspond to 41.1 million base pairs, by using a combined assembly strategy. Transcripts were annotated by Blast2GO and 51,763 transcripts got BLASTX hits, in which 39,677 transcripts have GO terms and 14,117 have ECs that are associated with 147 KEGG pathways. Furthermore, transcriptome differences of seven tissues were analyzed by using Illumina digital gene expression (DGE) tag profiling and numerous differentially and specifically expressed transcripts were identified. Moreover, the expression characteristics of genes involved in viral genomes, starch metabolism and potential stress tolerance and insect resistance were also identified.

CONCLUSIONS/SIGNIFICANCE

The combined de novo transcriptome assembly strategy can be applied to other organisms whose reference genomes are not available. The data provided here represent the most comprehensive and integrated genomic resources for cloning and identifying genes of interest in sweet potato. Characterization of sweet potato transcriptome provides an effective tool for better understanding the molecular mechanisms of cellular processes including development of leaves and storage roots, tissue-specific gene expression, potential biotic and abiotic stress response in sweet potato.

Comparative proteomic analysis reveals that caspase-1 and serine protease may be involved in silkworm resistance to Bombyx mori nuclear polyhedrosis virus

The silkworm Bombyx mori is of great economic value. The B. mori nuclear polyhedrosis virus (BmNPV) is one of the most common and severe pathogens for silkworm. Although certain immune mechanisms exist in silkworms, most silkworms are still susceptible to BmNPV infection. Interestingly, BmNPV infection resistance in some silkworm strains is varied and naturally existing. We have previously established a silkworm strain NB by genetic cross, which is highly resistant to BmNPV invasion. To investigate the molecular mechanism of silkworm resistance to BmNPV infection, we employed proteomic approach and genetic cross to globally identify proteins differentially expressed in parental silkworms NB and 306, a BmNPV-susceptible strain, and their F(1) hybrids. In all, 53 different proteins were found in direct cross group (NB♀, 306♂, F(1) hybrid) and 21 in reciprocal cross group (306♀, NB♂, F(1) hybrid). Gene ontology and KEGG pathway analyses showed that most of these different proteins are located in cytoplasm and are involved in many important metabolisms. Caspase-1 and serine protease expressed only in BmNPV-resistant silkworms, but not in BmNPV-susceptible silkworms, which was further confirmed by Western blot. Taken together, our data suggests that both caspase-1 and serine protease play a critical role in silkworm resistance against BmNPV infection.

Transcriptome analysis of stem development in the tumourous stem mustard Brassica juncea var. tumida Tsen et Lee by RNA sequencing

BACKGROUND

Tumourous stem mustard (Brassica juncea var. tumida Tsen et Lee) is an economically and nutritionally important vegetable crop of the Cruciferae family that also provides the raw material for Fuling mustard. The genetics breeding, physiology, biochemistry and classification of mustards have been extensively studied, but little information is available on tumourous stem mustard at the molecular level. To gain greater insight into the molecular mechanisms underlying stem swelling in this vegetable and to provide additional information for molecular research and breeding, we sequenced the transcriptome of tumourous stem mustard at various stem developmental stages and compared it with that of a mutant variety lacking swollen stems.

RESULTS

Using Illumina short-read technology with a tag-based digital gene expression (DGE) system, we performed de novo transcriptome assembly and gene expression analysis. In our analysis, we assembled genetic information for tumourous stem mustard at various stem developmental stages. In addition, we constructed five DGE libraries, which covered the strains Yong'an and Dayejie at various development stages. Illumina sequencing identified 146,265 unigenes, including 11,245 clusters and 135,020 singletons. The unigenes were subjected to a BLAST search and annotated using the GO and KO databases. We also compared the gene expression profiles of three swollen stem samples with those of two non-swollen stem samples. A total of 1,042 genes with significantly different expression levels occurring simultaneously in the six comparison groups were screened out. Finally, the altered expression levels of a number of randomly selected genes were confirmed by quantitative real-time PCR.

CONCLUSIONS

Our data provide comprehensive gene expression information at the transcriptional level and the first insight into the understanding of the molecular mechanisms and regulatory pathways of stem swelling and development in this plant, and will help define new mechanisms of stem development in non-model plant organisms.

Annotation of the transcriptome from Taenia pisiformis and its comparative analysis with three Taeniidae species

Background

Taenia pisiformis is one of the most common intestinal tapeworms and can cause infections in canines. Adult T. pisiformis (canines as definitive hosts) and Cysticercus pisiformis (rabbits as intermediate hosts) cause significant health problems to the host and considerable socio-economic losses as a consequence. No complete genomic data regarding T. pisiformis are currently available in public databases. RNA-seq provides an effective approach to analyze the eukaryotic transcriptome to generate large functional gene datasets that can be used for further studies.

Methodology/Principal Findings

In this study, 2.67 million sequencing clean reads and 72,957 unigenes were generated using the RNA-seq technique. Based on a sequence similarity search with known proteins, a total of 26,012 unigenes (no redundancy) were identified after quality control procedures via the alignment of four databases. Overall, 15,920 unigenes were mapped to 203 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Through analyzing the glycolysis/gluconeogenesis and axonal guidance pathways, we achieved an in-depth understanding of the biochemistry of T. pisiformis. Here, we selected four unigenes at random and obtained their full-length cDNA clones using RACE PCR. Functional distribution characteristics were gained through comparing four cestode species (72,957 unigenes of T. pisiformis, 30,700 ESTs of T. solium, 1,058 ESTs of Eg+Em [conserved ESTs between Echinococcus granulosus and Echinococcus multilocularis]), with the cluster of orthologous groups (COG) and gene ontology (GO) functional classification systems. Furthermore, the conserved common genes in these four cestode species were obtained and aligned by the KEGG database.

Conclusion

This study provides an extensive transcriptome dataset obtained from the deep sequencing of T. pisiformis in a non-model whole genome. The identification of conserved genes may provide novel approaches for potential drug targets and vaccinations against cestode infections. Research can now accelerate into the functional genomics, immunity and gene expression profiles of cestode species.

Transcriptome profiling of testis during sexual maturation stages in Eriocheir sinensis using Illumina sequencing

The testis is a highly specialized tissue that plays dual roles in ensuring fertility by producing spermatozoa and hormones. Spermatogenesis is a complex process, resulting in the production of mature sperm from primordial germ cells. Significant structural and biochemical changes take place in the seminiferous epithelium of the adult testis during spermatogenesis. The gene expression pattern of testis in Chinese mitten crab (Eriocheir sinensis) has not been extensively studied, and limited genetic research has been performed on this species. The advent of high-throughput sequencing technologies enables the generation of genomic resources within a short period of time and at minimal cost. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset for testis of E. sinensis. In two runs, we produced 25,698,778 sequencing reads corresponding with 2.31 Gb total nucleotides. These reads were assembled into 342,753 contigs or 141,861 scaffold sequences, which identified 96,311 unigenes. Based on similarity searches with known proteins, 39,995 unigenes were annotated based on having a Blast hit in the non-redundant database or ESTscan results with a cut-off E-value above 10⁻⁵. This is the first report of a mitten crab transcriptome using high-throughput sequencing technology, and all these testes transcripts can help us understand the molecular mechanisms involved in spermatogenesis and testis maturation.

Transcriptome analysis reveals the time of the fourth round of genome duplication in common carp (Cyprinus carpio)

Background

Common carp (Cyprinus carpio) is thought to have undergone one extra round of genome duplication compared to zebrafish. Transcriptome analysis has been used to study the existence and timing of genome duplication in species for which genome sequences are incomplete. Large-scale transcriptome data for the common carp genome should help reveal the timing of the additional duplication event.

Results

We have sequenced the transcriptome of common carp using 454 pyrosequencing. After assembling the 454 contigs and the published common carp sequences together, we obtained 49,669 contigs and identified genes using homology searches and an ab initio method. We identified 4,651 orthologous pairs between common carp and zebrafish and found 129,984 paralogous pairs within the common carp. An estimation of the synonymous substitution rate in the orthologous pairs indicated that common carp and zebrafish diverged 120 million years ago (MYA). We identified one round of genome duplication in common carp and estimated that it had occurred 5.6 to 11.3 MYA. In zebrafish, no genome duplication event after speciation was observed, suggesting that, compared to zebrafish, common carp had undergone an additional genome duplication event. We annotated the common carp contigs with Gene Ontology terms and KEGG pathways. Compared with zebrafish gene annotations, we found that a set of biological processes and pathways were enriched in common carp.

Conclusions

The assembled contigs helped us to estimate the time of the fourth-round of genome duplication in common carp. The resource that we have built as part of this study will help advance functional genomics and genome annotation studies in the future.

De novo assembly and characterisation of the transcriptome during seed development, and generation of genic-SSR markers in peanut (Arachis hypogaea L.)

Background

The peanut (Arachis hypogaea L.) is an important oilseed crop in tropical and subtropical regions of the world. However, little about the molecular biology of the peanut is currently known. Recently, next-generation sequencing technology, termed RNA-seq, has provided a powerful approach for analysing the transcriptome, and for shedding light on the molecular biology of peanut.

Results

In this study, we employed RNA-seq to analyse the transcriptomes of the immature seeds of three different peanut varieties with different oil contents. A total of 26.1-27.2 million paired-end reads with lengths of 100 bp were generated from the three varieties and 59,077 unigenes were assembled with N50 of 823 bp. Based on sequence similarity search with known proteins, a total of 40,100 genes were identified. Among these unigenes, only 8,252 unigenes were annotated with 42 gene ontology (GO) functional categories. And 18,028 unigenes mapped to 125 pathways by searching against the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). In addition, 3,919 microsatellite markers were developed in the unigene library, and 160 PCR primers of SSR loci were used for validation of the amplification and the polymorphism.

Conclusion

We completed a successful global analysis of the peanut transcriptome using RNA-seq, a large number of unigenes were assembled, and almost four thousand SSR primers were developed. These data will facilitate gene discovery and functional genomic studies of the peanut plant. In addition, this study provides insight into the complex transcriptome of the peanut and established a biotechnological platform for future research.

Ten years of pathway analysis: current approaches and outstanding challenges

Pathway analysis has become the first choice for gaining insight into the underlying biology of differentially expressed genes and proteins, as it reduces complexity and has increased explanatory power. We discuss the evolution of knowledge base–driven pathway analysis over its first decade, distinctly divided into three generations. We also discuss the limitations that are specific to each generation, and how they are addressed by successive generations of methods. We identify a number of annotation challenges that must be addressed to enable development of the next generation of pathway analysis methods. Furthermore, we identify a number of methodological challenges that the next generation of methods must tackle to take advantage of the technological advances in genomics and proteomics in order to improve specificity, sensitivity, and relevance of pathway analysis.

De novo transcriptome of safflower and the identification of putative genes for oleosin and the biosynthesis of flavonoids

Safflower (Carthamus tinctorius L.) is one of the most extensively used oil crops in the world. However, little is known about how its compounds are synthesized at the genetic level. In this study, Solexa-based deep sequencing on seed, leaf and petal of safflower produced a de novo transcriptome consisting of 153,769 unigenes. We annotated 82,916 of the unigenes with gene annotation and assigned functional terms and specific pathways to a subset of them. Metabolic pathway analysis revealed that 23 unigenes were predicted to be responsible for the biosynthesis of flavonoids and 8 were characterized as seed-specific oleosins. In addition, a large number of differentially expressed unigenes, for example, those annotated as participating in anthocyanin and chalcone synthesis, were predicted to be involved in flavonoid biosynthesis pathways. In conclusion, the de novo transcriptome investigation of the unique transcripts provided candidate gene resources for studying oleosin-coding genes and for investigating genes related to flavonoid biosynthesis and metabolism in safflower.

Transcriptome characterization and sequencing-based identification of salt-responsive genes in Millettia pinnata, a semi-mangrove plant

Semi-mangroves form a group of transitional species between glycophytes and halophytes, and hold unique potential for learning molecular mechanisms underlying plant salt tolerance. Millettia pinnata is a semi-mangrove plant that can survive a wide range of saline conditions in the absence of specialized morphological and physiological traits. By employing the Illumina sequencing platform, we generated ~192 million short reads from four cDNA libraries of M. pinnata and processed them into 108,598 unisequences with a high depth of coverage. The mean length and total length of these unisequences were 606 bp and 65.8 Mb, respectively. A total of 54,596 (50.3%) unisequences were assigned Nr annotations. Functional classification revealed the involvement of unisequences in various biological processes related to metabolism and environmental adaptation. We identified 23,815 candidate salt-responsive genes with significantly differential expression under seawater and freshwater treatments. Based on the reverse transcription-polymerase chain reaction (RT-PCR) and real-time PCR analyses, we verified the changes in expression levels for a number of candidate genes. The functional enrichment analyses for the candidate genes showed tissue-specific patterns of transcriptome remodelling upon salt stress in the roots and the leaves. The transcriptome of M. pinnata will provide valuable gene resources for future application in crop improvement. In addition, this study sets a good example for large-scale identification of salt-responsive genes in non-model organisms using the sequencing-based approach.

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.