Deep sequencing of the zebrafish transcriptome response to mycobacterium infection

De novo characterization of Japanese scallop Mizuhopecten yessoensis transcriptome and analysis of its gene expression following cadmium exposure

Background

Japanese scallop has been cultured on a large-scale in China for many years. However, serious marine pollution in recent years has resulted in considerable loss to this industry. Moreover, due to the lack of genomic resources, limited research has been carried out on this species. To facilitate the understanding at molecular level immune and stress response mechanism, an extensive transcriptomic profiling and digital gene expression (DGE) database of Japanese scallop upon cadmium exposure was carried out using the Illumina sequencing platform.

Results

RNA-seq produced about 112 million sequencing reads from the tissues of adult Japanese scallops. These reads were assembled into 194,839 non-redundant sequences with open reading frame (ORF), of which 14,240 putative amino acid sequences were assigned biological function annotation and were annotated with gene ontology and eukaryotic orthologous group terms. In addition, we identified 720 genes involved in response to stimulus and 302 genes involved in immune-response pathways. Furthermore, we investigated the transcriptomic changes in the gill and digestive gland of Japanese scallops following cadmium exposure using a tag-based DGE system. A total of 7,556 and 3,002 differentially expressed genes were detected, respectively, and functionally annotated with KEGG pathway annotations.

Conclusion

This study provides a comprehensive transcripts sequence resource for the Japanese scallop and presents a survey of gene expression in response to heavy metal exposure in a non-model marine invertebrate via the Illumina sequencing platform. These results may contribute to the in-depth elucidation of the molecular mechanisms involved in bivalve responses to marine pollutants.

Deep mRNA sequencing analysis to capture the transcriptome landscape of zebrafish embryos and larvae

Transcriptome analysis is a powerful tool to obtain large amount genome-scale gene expression profiles. Despite its extensive usage to diverse biological problems in the last decade, transcriptomic researches approaching the zebrafish embryonic development have been very limited. Several recent studies have made great progress in this direction, yet the large gap still exists, especially regarding to the transcriptome dynamics of embryonic stages from early gastrulation onwards. Here, we present a comprehensive analysis about the transcriptomes of 9 different stages covering 7 major periods (cleavage, blastula, gastrula, segmentation, pharyngula, hatching and early larval stage) in zebrafish development, by recruiting the RNA-sequencing technology. We detected the expression for at least 24,065 genes in at least one of the 9 stages. We identified 16,130 genes that were significantly differentially expressed between stages and were subsequently classified into six clusters. Each revealed gene cluster had distinct expression patterns and characteristic functional pathways, providing a framework for the understanding of the developmental transcriptome dynamics. Over 4000 genes were identified as preferentially expressed in one of the stages, which could be of high relevance to stage-specific developmental and molecular events. Among the 68 transcription factor families active during development, most had enhanced average expression levels and thus might be crucial for embryogenesis, whereas the inactivation of the other families was likely required by the activation of the zygotic genome. We discussed our RNA-seq data together with previous findings about the Wnt signaling pathway and some other genes with known functions, to show how our data could be used to advance our understanding about these developmental functional elements. Our study provides ample information for further study about the molecular and cellular mechanisms underlying vertebrate development.

Analysis of the transcriptional differences between indigenous and invasive whiteflies reveals possible mechanisms of whitefly invasion

BACKGROUND

The whitefly Bemisa tabaci is a species complex of more than 31 cryptic species which include some of the most destructive invasive pests of crops worldwide. Among them, Middle East-Asia Minor 1 (MEAM1) and Mediterranean have invaded many countries and displaced the native whitefly species. The successful invasion of the two species is largely due to their wide range of host plants, high resistance to insecticides and remarkable tolerance to environmental stresses. However, the molecular differences between invasive and indigenous whiteflies remain largely unknown.

METHODOLOGY/PRINCIPAL FINDINGS

Here the global transcriptional difference between the two invasive whitefly species (MEAM1, MED) and one indigenous whitefly species (Asia II 3) were analyzed using the Illumina sequencing. Our analysis indicated that 2,422 genes between MEAM1 and MED; 3,073 genes between MEAM1 and Asia II 3; and 3,644 genes between MED and Asia II 3 were differentially expressed. Gene Ontology enrichment analysis revealed that the differently expressed genes between the invasive and indigenous whiteflies were significantly enriched in the term of 'oxidoreductase activity'. Pathway enrichment analysis showed that carbohydrate, amino acid and glycerolipid metabolisms were more active in MEAM1 and MED than in Asia II 3, which may contribute to their differences in biological characteristics. Our analysis also illustrated that the majority of genes involved in 'drug metabolic pathway' were expressed at a higher level in MEAM1 and MED than in Asia II 3. Taken together, these results revealed that the genes related to basic metabolism and detoxification were expressed at an elevated level in the invasive whiteflies, which might be responsible for their higher resistance to insecticides and environmental stresses.

CONCLUSIONS/SIGNIFICANCE

The extensive comparison of MEAM1, MED and Asia II 3 gene expression may serve as an invaluable resource for revealing the molecular mechanisms underlying their biological differences and the whitefly invasion.

De novo transcriptome sequencing and comparative analysis of differentially expressed genes in Gossypium aridum under salt stress

Salinity stress is one of the most serious factors that impede the growth and development of various crops. Wild Gossypium species, which are remarkably tolerant to salt water immersion, are valuable resources for understanding salt tolerance mechanisms of Gossypium and improving salinity resistance in upland cotton. To generate a broad survey of genes with altered expression during various stages of salt stress, a mixed RNA sample was prepared from the roots and leaves of Gossypium aridum plants subjected to salt stress. The transcripts were sequenced using the Illumina sequencing platform. After cleaning and quality checks, approximately 41.5 million clean reads were obtained. Finally, these reads were eventually assembled into 98,989 unigenes with a mean size of 452 bp. All unigenes were compared to known cluster of orthologous groups (COG) sequences to predict and classify the possible functions of these genes, which were classified into at least 25 molecular families. Variations in gene expression were then examined after exposing the plants to 200 mM NaCl for 3, 12, 72 or 144 h. Sequencing depths of approximately six million raw tags were achieved for each of the five stages of salt stress. There were 2634 (1513 up-regulated/1121 down-regulated), 2449 (1586 up-regulated/863 down-regulated), 2271 (946 up-regulated/1325 down-regulated) and 3352 (933 up-regulated/2419 down-regulated) genes that were differentially expressed after exposure to NaCl for 3, 12, 72 and 144 h, respectively. Digital gene expression analysis indicated that pathways involved in "transport", "response to hormone stimulus" and "signaling" play important roles during salt stress, while genes involved in "protein kinase activity" and "transporter activity" undergo major changes in expression during early and later stages of salt stress, respectively.

Genome expression profile analysis reveals important transcripts in maize roots responding to the stress of heavy metal Pb

Lead (Pb) has become one of the most abundant heavy metal pollutants of the environment. With its large biomass, maize could be an important object for studying the phytoremediation of Pb-contaminated soil. In our previous research, we screened 19 inbred lines of maize for Pb concentration, and line 178 was identified to be a hyperaccumulator for Pb in both the roots and aboveground parts. To identify important genes and metabolic pathways related to Pb accumulation and tolerance, line 178 was underwent genome expression profile under Pb stress and a control (CK). A total of approximately 11 million cDNA tags were sequenced and 4 665 539 and 4 936 038 clean tags were obtained from the libraries of the test and CK, respectively. In comparison to CK, 2379 and 1832 genes were identified up- or downregulated, respectively, more than fivefolds under Pb stress. Interestingly, all the genes were related to cellular processes and signaling, information storage and processing or metabolism functions. Particularly, the genes involved in posttranslational modification, protein turnover and chaperones; signal transduction, carbohydrate transport and metabolism; and lipid transport and metabolism significantly changed under the treatment. In addition, seven pathways including ribosome, photosynthesis, and carbon fixation were affected significantly, with 118, 12, 34, 21, 18, 72 and 43 differentially expressed genes involved. The significant upregulation of the ribosome pathway may reveal an important secret for Pb tolerance of line 178. And the sharp increase of laccase transcripts and metal ion transporters were suggested to account in part for Pb hyperaccumulation in the line.

Transcriptome profiling of the goose (Anser cygnoides) ovaries identify laying and broodiness phenotypes

Background

The geese have strong broodiness and poor egg performance. These characteristics are the key issues that hinder the goose industry development. Yet little is known about the mechanisms responsible for follicle development due to lack of genomic resources. Hence, studies based on high-throughput sequencing technologies are needed to produce a comprehensive and integrated genomic resource and to better understand the biological mechanisms of goose follicle development.

Methodology/Principal Findings

In this study, we performed de novo transcriptome assembly and gene expression analysis using short-read sequencing technology (Illumina). We obtained 67,315,996 short reads of 100 bp, which were assembled into 130,514 unique sequences by Trinity strategy (mean size = 753bp). Based on BLAST results with known proteins, these analyses identified 52,642 sequences with a cut-off E-value above 10⁻⁵. Assembled sequences were annotated with gene descriptions, gene ontology and clusters of orthologous group terms. In addition, we investigated the transcription changes during the goose laying/broodiness period using a tag-based digital gene expression (DGE) system. We obtained a sequencing depth of over 4.2 million tags per sample and identified a large number of genes associated with follicle development and reproductive biology including cholesterol side-chain cleavage enzyme gene and dopamine beta-hydroxylas gene. We confirm the altered expression levels of the two genes using quantitative real-time PCR (qRT-PCR).

Conclusions/Significance

The obtained goose transcriptome and DGE profiling data provide comprehensive gene expression information at the transcriptional level that could promote better understanding of the molecular mechanisms underlying follicle development and productivity.

Sequencing of the first ayu (Plecoglossus altivelis) macrophage transcriptome and microarray development for investigation the effect of LECT2 on macrophages

Macrophages play an important role in first-line host defense of innate immune in fishes. However, it is difficult to investigate cellular mechanism of immune response in fish species with little genomic information available. Here we present the first use of RNA-Sequencing to study the macrophage transcriptome of ayu, Plecoglossus altivelis, which is an economically important fish in East Asia. De novo assembly generated 49,808 non-redundant consensus sequences, among which 23,490 transcripts found respective coding sequences. 15,707 transcripts are predicted to be involved in known metabolic or signaling pathways. The sequences were then used to develop a microarray for measurement the effect of recombinant LECT2 on ayu macrophages. LECT2 altered expression of a variety of genes mainly implicated in actin cytoskeleton, pattern recognition receptors and cytokines. Meanwhile, LECT2 enhanced phagocytosis, bacterial killing, and respiratory burst in ayu macrophages, which supported the thought derived from the microarray data that LECT2 activates macrophages. In conclusion, our results contribute to understanding the specific regulation mechanism of LECT2 in macrophage activation, and the combination of transcriptome analysis and microarray assay is a good method for screening a special tissue or cell response to a stimulus or pathogen in non-model fish species.

De novo assembly and characterization of spotted seal Phoca largha transcriptome using Illumina paired-end sequencing

Spotted seal (Phoca largha) is categorized as a critically endangered species in China. The aim of this study was to investigate spotted seal transcriptome by the approach of Illumina paired-end sequencing technology. We obtained a total of 52,146,394 reads for the mixed tissues of liver and spleen from the spotted seal. The de novo assemblies yielded 354,014 contigs and 178,466 unigenes. In the transcriptome, 193 unigenes were assigned to defense mechanisms. Three unigenes encoded MHC class I and 17 unigenes encoded MHC class II. In addition, bioinformatics analysis revealed a total of 4425 simple sequence repeats (SSRs). Fifty SSRs were randomly selected to validate amplification and determine the degree of polymorphism in the genomic DNA pools. Thirty-five primer pairs successfully amplified the expected DNA fragments and detected significant polymorphism among 28 spotted seal individuals. These results would contribute to the understanding of the genetic makeup of spotted seal transcriptome and provide useful information for functional genomic research in this species.

Deep RNA sequencing of the skeletal muscle transcriptome in swimming fish

Deep RNA sequencing (RNA-seq) was performed to provide an in-depth view of the transcriptome of red and white skeletal muscle of exercised and non-exercised rainbow trout (Oncorhynchus mykiss) with the specific objective to identify expressed genes and quantify the transcriptomic effects of swimming-induced exercise. Pubertal autumn-spawning seawater-raised female rainbow trout were rested (n = 10) or swum (n = 10) for 1176 km at 0.75 body-lengths per second in a 6,000-L swim-flume under reproductive conditions for 40 days. Red and white muscle RNA of exercised and non-exercised fish (4 lanes) was sequenced and resulted in 15–17 million reads per lane that, after de novo assembly, yielded 149,159 red and 118,572 white muscle contigs. Most contigs were annotated using an iterative homology search strategy against salmonid ESTs, the zebrafish Danio rerio genome and general Metazoan genes. When selecting for large contigs (>500 nucleotides), a number of novel rainbow trout gene sequences were identified in this study: 1,085 and 1,228 novel gene sequences for red and white muscle, respectively, which included a number of important molecules for skeletal muscle function. Transcriptomic analysis revealed that sustained swimming increased transcriptional activity in skeletal muscle and specifically an up-regulation of genes involved in muscle growth and developmental processes in white muscle. The unique collection of transcripts will contribute to our understanding of red and white muscle physiology, specifically during the long-term reproductive migration of salmonids.

Different transcriptional profiles of RAW264.7 infected with Mycobacterium tuberculosis H37Rv and BCG identified via deep sequencing

BACKGROUND

The Mycobacterium tuberculosis H37Rv and BCG effects on the host cell transcriptional profile consider a main research point. In the present study the transcriptome profiling analysis of RAW264.7 either infected with Mycobacterium tuberculosis H37Rv or BCG have been reported using Solexa/Illumina digital gene expression (DGE).

RESULTS

The DGE analysis showed 1,917 different expressed genes between the BCG and H37Rv group. In addition, approximately 5% of the transcripts appeared to be predicted genes that have never been described before. KEGG Orthology (KO) annotations showed more than 71% of these transcripts are possibly involved in approximately 210 known metabolic or signaling pathways. The gene of the 28 pathways about pathogen recognition receptors and Mycobacterium tuberculosis interaction with macrophages were analyzed using the CLUSTER 3.0 available, the Tree View tool and Gene Orthology (GO). Some genes were randomly selected to confirm their altered expression levels by quantitative real-time PCR (qRT-PCR).

CONCLUSION

The present study used DGE from pathogen recognition receptors and Mycobacterium tuberculosis interaction with macrophages to understand the interplay between Mycobacterium tuberculosis and RAW264.7. Meanwhile find some important host protein which was affected by Mycobacterium tuberculosis to provide evidence for the further improvement of the present efficacy of existing Mycobacterium tuberculosis therapy and vaccine.

Transcriptome characterization and gene expression of Epinephelus spp in endoplasmic reticulum stress-related pathway during betanodavirus infection in vitro

Background

Grouper (Epinephelus spp) is an economically important fish species worldwide. However, viral pathogens such as nervous necrosis virus (NNV) have been causing severe infections in the fish, resulting in great loss in the grouper aquaculture industry. Yet, the understanding of the molecular mechanisms underlying the pathogenicity of NNV is still inadequate, mainly due to insufficient genomic information of the host.

Results

De novo assembly of grouper transcriptome in the grouper kidney (GK) cells was conducted by using short read sequencing technology of Solexa/Illumina. A sum of 66,582 unigenes with mean length of 603 bp were obtained, and were annotated according to Gene Ontology (GO) and Clusters of Orthologous Groups (COG). In addition, the tag-based digital gene expression (DGE) system was used to investigate the gene expression and pathways associated with NNV infection in GK cells. The analysis revealed endoplasmic reticulum (ER) stress response was prominently affected in NNV-infected GK cells. A further analysis revealed an interaction between the NNV capsid protein and the ER chaperone immunoglobulin heavy-chain binding protein (BiP). Furthermore, exogenous expression of NNV capsid protein was able to induce XBP-1 mRNA splicing in vivo, suggesting a role of the capsid protein in the NNV-induced ER stress.

Conclusions

Our data presents valuable genetic information for Epinephelus spp., which will benefit future study in this non-model but economically important species. The DGE profile of ER stress response in NNV-infected cells provides information of many important components associated with the protein processing in ER. Specifically, we showed that the viral capsid protein might play an important role in the ER stress response.

De novo assembly, characterization and functional annotation of pineapple fruit transcriptome through massively parallel sequencing

Background

Pineapple (Ananas comosus var. comosus), is an important tropical non-climacteric fruit with high commercial potential. Understanding the mechanism and processes underlying fruit ripening would enable scientists to enhance the improvement of quality traits such as, flavor, texture, appearance and fruit sweetness. Although, the pineapple is an important fruit, there is insufficient transcriptomic or genomic information that is available in public databases. Application of high throughput transcriptome sequencing to profile the pineapple fruit transcripts is therefore needed.

Methodology/Principal Findings

To facilitate this, we have performed transcriptome sequencing of ripe yellow pineapple fruit flesh using Illumina technology. About 4.7 millions Illumina paired-end reads were generated and assembled using the Velvet de novo assembler. The assembly produced 28,728 unique transcripts with a mean length of approximately 200 bp. Sequence similarity search against non-redundant NCBI database identified a total of 16,932 unique transcripts (58.93%) with significant hits. Out of these, 15,507 unique transcripts were assigned to gene ontology terms. Functional annotation against Kyoto Encyclopedia of Genes and Genomes pathway database identified 13,598 unique transcripts (47.33%) which were mapped to 126 pathways. The assembly revealed many transcripts that were previously unknown.

Conclusions

The unique transcripts derived from this work have rapidly increased of the number of the pineapple fruit mRNA transcripts as it is now available in public databases. This information can be further utilized in gene expression, genomics and other functional genomics studies in pineapple.

Transcriptome analysis of Nicotiana tabacum infected by Cucumber mosaic virus during systemic symptom development

Virus infection of plants may induce a variety of disease symptoms. However, little is known about the molecular mechanism of systemic symptom development in infected plants. Here we performed the first next-generation sequencing study to identify gene expression changes associated with disease development in tobacco plants (Nicotiana tabacum cv. Xanthi nc) induced by infection with the M strain of Cucumber mosaic virus (M-CMV). Analysis of the tobacco transcriptome by RNA-Seq identified 95,916 unigenes, 34,408 of which were new transcripts by database searches. Deep sequencing was subsequently used to compare the digital gene expression (DGE) profiles of the healthy plants with the infected plants at six sequential disease development stages, including vein clearing, mosaic, severe chlorosis, partial and complete recovery, and secondary mosaic. Thousands of differentially expressed genes were identified, and KEGG pathway analysis of these genes suggested that many biological processes, such as photosynthesis, pigment metabolism and plant-pathogen interaction, were involved in systemic symptom development. Our systematic analysis provides comprehensive transcriptomic information regarding systemic symptom development in virus-infected plants. This information will help further our understanding of the detailed mechanisms of plant responses to viral infection.

Transcriptome analysis of the roots at early and late seedling stages using Illumina paired-end sequencing and development of EST-SSR markers in radish

The tuberous root of radish is an important vegetable, but insufficient transcriptomic and genomic data are currently available to understand the molecular mechanisms underlying tuberous root formation and development. High-throughput transcriptomic sequencing is essential to generate a large transcript sequence data set for gene discovery and molecular marker development. In this study, a total of 107.3 million clean reads were generated using Illumina paired-end sequencing technology. De novo assembly generated 61,554 unigenes with an average length of 820 bp. Based on a sequence similarity search with known proteins or nucleotides, 85.51 % (52,634), 90.18 % (55,507) and 54 % (33,242) consensus sequences showed homology with sequences in the Nr, Nt and Swiss-Prot databases, respectively. Of these annotated unigenes, 21,109 and 17,343 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. A total of 27,809 unigenes were assigned to 123 pathways in the Kyoto Encyclopedia of Genes and Genomes database. Analysis of transcript differences between libraries from the early and late seedling developmental stages demonstrated that starch and sucrose metabolism and phenylpropanoid biosynthesis may be the dominant metabolic events during tuberous root formation and plant hormones probably play critical roles in regulation of this developmental process. In total, 14,641 potential EST-SSRs were identified among the unigenes, and 12,733 primer pairs for 2,511 SSR were obtained. Summarily, this study gave us a clue to understand the radish tuberous root formation and development, and also provided us with a valuable sequence resource for novel gene discovery and marker-assisted selective breeding in radish.

KEY MESSAGE

De novo assembled and characterized the radish tuberous root transcriptome; explored the mechanism of radish tuberous root formation; development of EST-SSR markers in radish.

RNA-seq liver transcriptome analysis reveals an activated MHC-I pathway and an inhibited MHC-II pathway at the early stage of vaccine immunization in zebrafish

Background

Zebrafish (Danio rerio) is a prominent vertebrate model of human development and pathogenic disease and has recently been utilized to study teleost immune responses to infectious agents threatening the aquaculture industry. In this work, to clarify the host immune mechanisms underlying the protective effects of a putative vaccine and improve its immunogenicity in the future efforts, high-throughput RNA sequencing technology was used to investigate the immunization-related gene expression patterns of zebrafish immunized with Edwardsiella tarda live attenuated vaccine.

Results

Average reads of 18.13 million and 14.27 million were obtained from livers of zebrafish immunized with phosphate buffered saline (mock) and E. tarda vaccine (WED), respectively. The reads were annotated with the Ensembl zebrafish database before differential expressed genes sequencing (DESeq) comparative analysis, which identified 4565 significantly differentially expressed genes (2186 up-regulated and 2379 down-regulated in WED; p<0.05). Among those, functional classifications were found in the Gene Ontology database for 3891 and in the Kyoto Encyclopedia of Genes and Genomes database for 3467. Several pathways involved in acute phase response, complement activation, immune/defense response, and antigen processing and presentation were remarkably affected at the early stage of WED immunization. Further qPCR analysis confirmed that the genes encoding the factors involved in major histocompatibility complex (MHC)-I processing pathway were up-regulated, while those involved in MHC-II pathway were down-regulated.

Conclusion

These data provided insights into the molecular mechanisms underlying zebrafish immune response to WED immunization and might aid future studies to develop a highly immunogenic vaccine against gram-negative bacteria in teleosts.

Pathogen recognition and activation of the innate immune response in zebrafish

The zebrafish has proven itself as an excellent model to study vertebrate innate immunity. It presents us with possibilities for in vivo imaging of host-pathogen interactions which are unparalleled in mammalian model systems. In addition, its suitability for genetic approaches is providing new insights on the mechanisms underlying the innate immune response. Here, we review the pattern recognition receptors that identify invading microbes, as well as the innate immune effector mechanisms that they activate in zebrafish embryos. We compare the current knowledge about these processes in mammalian models and zebrafish and discuss recent studies using zebrafish infection models that have advanced our general understanding of the innate immune system. Furthermore, we use transcriptome analysis of zebrafish infected with E. tarda, S. typhimurium, and M. marinum to visualize the gene expression profiles resulting from these infections. Our data illustrate that the two acute disease-causing pathogens, E. tarda and S. typhimurium, elicit a highly similar proinflammatory gene induction profile, while the chronic disease-causing pathogen, M. marinum, induces a weaker and delayed innate immune response.

Characterization of the spotted seal Phoca largha transcriptome using Illumina paired-end sequencing and development of SSR markers

Next-generation sequencing provides a powerful new approach for developing functional genomic tools for non-model species. This study aims to investigate the spotted seal (Phoca largha) transcriptome by the approach of Illumina paired-end sequencing technology. We obtained a total of 52,146,394 reads for the mixed tissues of liver and spleen from spotted seal. The de novo assemblies yielded 354,014 contigs and 178,466 unigenes. In addition, bioinformatics analysis revealed a total of 4425 simple sequence repeats (SSRs). Fifty SSRs were randomly selected to validate amplification and determine the degree of polymorphism in the genomic DNA pools. Thirty-five primer pairs successfully amplified expected DNA fragments and detected significant polymorphism among 28 spotted seal individuals. These results contribute to the understanding of the genetic makeup of spotted seal transcriptome and provide useful information for functional genomic research in this species.

De novo assembly and characterization of bark transcriptome using Illumina sequencing and development of EST-SSR markers in rubber tree (Hevea brasiliensis Muell. Arg.)

Background

In rubber tree, bark is one of important agricultural and biological organs. However, the molecular mechanism involved in the bark formation and development in rubber tree remains largely unknown, which is at least partially due to lack of bark transcriptomic and genomic information. Therefore, it is necessary to carried out high-throughput transcriptome sequencing of rubber tree bark to generate enormous transcript sequences for the functional characterization and molecular marker development.

Results

In this study, more than 30 million sequencing reads were generated using Illumina paired-end sequencing technology. In total, 22,756 unigenes with an average length of 485 bp were obtained with de novo assembly. The similarity search indicated that 16,520 and 12,558 unigenes showed significant similarities to known proteins from NCBI non-redundant and Swissprot protein databases, respectively. Among these annotated unigenes, 6,867 and 5,559 unigenes were separately assigned to Gene Ontology (GO) and Clusters of Orthologous Group (COG). When 22,756 unigenes searched against the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) database, 12,097 unigenes were assigned to 5 main categories including 123 KEGG pathways. Among the main KEGG categories, metabolism was the biggest category (9,043, 74.75%), suggesting the active metabolic processes in rubber tree bark. In addition, a total of 39,257 EST-SSRs were identified from 22,756 unigenes, and the characterizations of EST-SSRs were further analyzed in rubber tree. 110 potential marker sites were randomly selected to validate the assembly quality and develop EST-SSR markers. Among 13 Hevea germplasms, PCR success rate and polymorphism rate of 110 markers were separately 96.36% and 55.45% in this study.

Conclusion

By assembling and analyzing de novo transcriptome sequencing data, we reported the comprehensive functional characterization of rubber tree bark. This research generated a substantial fraction of rubber tree transcriptome sequences, which were very useful resources for gene annotation and discovery, molecular markers development, genome assembly and annotation, and microarrays development in rubber tree. The EST-SSR markers identified and developed in this study will facilitate marker-assisted selection breeding in rubber tree. Moreover, this study also supported that transcriptome analysis based on Illumina paired-end sequencing is a powerful tool for transcriptome characterization and molecular marker development in non-model species, especially those with large and complex genomes.

Flowering time and transcriptome variation in Capsella bursa-pastoris (Brassicaceae)

• Flowering is a major developmental transition and its timing in relation to environmental conditions is of crucial importance to plant fitness. Understanding the genetic basis of flowering time variation is important to determining how plants adapt locally. • Here, we investigated flowering time variation of Capsella bursa-pastoris collected from different latitudes in China. We also used a digital gene expression (DGE) system to generate partial gene expression profiles for 12 selected samples. • We found that flowering time was highly variable and most strongly correlated with day length and winter temperature. Significant differences in gene expression between early- and late-flowering samples were detected for 72 candidate genes for flowering time. Genes related to circadian rhythms were significantly overrepresented among the differentially expressed genes. • Our data suggest that circadian rhythms and circadian clock genes play an important role in the evolution of flowering time, and C. bursa-pastoris plants exhibit expression differences for candidate genes likely to affect flowering time across the broad range of environments they face in China.

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.

Genome expression profile analysis of the immature maize embryo during dedifferentiation

Maize is one of the most important cereal crops worldwide and one of the primary targets of genetic manipulation, which provides an excellent way to promote its production. However, the obvious difference of the dedifferentiation frequency of immature maize embryo among various genotypes indicates that its genetic transformation is dependence on genotype and immature embryo-derived undifferentiated cells. To identify important genes and metabolic pathways involved in forming of embryo-derived embryonic calli, in this study, DGE (differential gene expression) analysis was performed on stages I, II, and III of maize inbred line 18-599R and corresponding control during the process of immature embryo dedifferentiation. A total of ∼21 million cDNA tags were sequenced, and 4,849,453, 5,076,030, 4,931,339, and 5,130,573 clean tags were obtained in the libraries of the samples and the control, respectively. In comparison with the control, 251, 324 and 313 differentially expressed genes (DEGs) were identified in the three stages with more than five folds, respectively. Interestingly, it is revealed that all the DEGs are related to metabolism, cellular process, and signaling and information storage and processing functions. Particularly, the genes involved in amino acid and carbohydrate transport and metabolism, cell wall/membrane/envelope biogenesis and signal transduction mechanism have been significantly changed during the dedifferentiation. To our best knowledge, this study is the first genome-wide effort to investigate the transcriptional changes in dedifferentiation immature maize embryos and the identified DEGs can serve as a basis for further functional characterization.

Identification of vascular and hematopoietic genes downstream of etsrp by deep sequencing in zebrafish

The transcription factor etsrp/Er71/Etv2 is a master control gene for vasculogenesis in all species studied to date. It is also required for hematopoiesis in zebrafish and mice. Several novel genes expressed in vasculature have been identified through transcriptional profiling of zebrafish embryos overexpressing etsrp by microarrays. Here we re-examined this transcriptional profile by Illumina RNA-sequencing technology, revealing a substantially increased number of candidate genes regulated by etsrp. Expression studies of 50 selected candidate genes from this dataset resulted in the identification of 39 new genes that are expressed in vascular cells. Regulation of these genes by etsrp was confirmed by their ectopic induction in etsrp overexpressing and decreased expression in etsrp deficient embryos. Our studies demonstrate the effectiveness of the RNA-sequencing technology to identify biologically relevant genes in zebrfish and produced a comprehensive profile of genes previously unexplored in vascular endothelial cell biology.

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.

Genome sequencing of ovine isolates of Mycobacterium avium subspecies paratuberculosis offers insights into host association

Background

The genome of Mycobacterium avium subspecies paratuberculosis (MAP) is remarkably homogeneous among the genomes of bovine, human and wildlife isolates. However, previous work in our laboratories with the bovine K-10 strain has revealed substantial differences compared to sheep isolates. To systematically characterize all genomic differences that may be associated with the specific hosts, we sequenced the genomes of three U.S. sheep isolates and also obtained an optical map.

Results

Our analysis of one of the isolates, MAP S397, revealed a genome 4.8 Mb in size with 4,700 open reading frames (ORFs). Comparative analysis of the MAP S397 isolate showed it acquired approximately 10 large sequence regions that are shared with the human M. avium subsp. hominissuis strain 104 and lost 2 large regions that are present in the bovine strain. In addition, optical mapping defined the presence of 7 large inversions between the bovine and ovine genomes (~ 2.36 Mb). Whole-genome sequencing of 2 additional sheep strains of MAP (JTC1074 and JTC7565) further confirmed genomic homogeneity of the sheep isolates despite the presence of polymorphisms on the nucleotide level.

Conclusions

Comparative sequence analysis employed here provided a better understanding of the host association, evolution of members of the M. avium complex and could help in deciphering the phenotypic differences observed among sheep and cattle strains of MAP. A similar approach based on whole-genome sequencing combined with optical mapping could be employed to examine closely related pathogens. We propose an evolutionary scenario for M. avium complex strains based on these genome sequences.

Global transcriptome profiling of the pine shoot beetle, Tomicus yunnanensis (Coleoptera: Scolytinae)

Background

The pine shoot beetle Tomicus yunnanensis (Coleoptera: Scolytinae) is an economically important pest of Pinus yunnanensis in southwestern China. Developed resistance to insecticides due to chemical pesticides being used for a long time is a factor involved in its serious damage, which poses a challenge for management. In addition, highly efficient adaptation to divergent environmental ecologies results in this pest posing great potential threat to pine forests. However, the molecular mechanisms remain unknown as only limited nucleotide sequence data for this species is available.

Methodology/Principal Findings

In this study, we applied next generation sequencing (Illumina sequencing) to sequence the adult transcriptome of T. yunnanensis. A total of 51,822,230 reads were obtained. They were assembled into 140,702 scaffolds, and 60,031 unigenes. The unigenes were further functionally annotated with gene descriptions, Gene Ontology (GO), Clusters of Orthologous Groups (COG), and Kyoto Encyclopedia of Genes and Genome (KEGG). In total, 80,932 unigenes were classified into GO, 13,599 unigenes were assigned to COG, and 33,875 unigenes were found in KO categories. A biochemical pathway database containing 219 predicted pathways was also created based on the annotations. In depth analysis of the data revealed a large number of genes related to insecticides resistance and heat shock protein genes associated with environmental stress.

Conclusions/Significance

The results facilitate the investigations of molecular resistance mechanisms to insecticides and environmental stress. This study lays the foundation for future functional genomics studies of important biological questions of this pest.

Zebrafish: model for the study of inflammation and the innate immune response to infectious diseases

The zebrafish (Danio rerio) has been extensively used in biomedical research as a model to study vertebrate development and hematopoiesis and recently, it has been adopted into varied fields including immunology. After fertilization, larvae survive with only the innate immune responses because adaptive immune system is morphologically and functionally mature only after 4-6 weeks postfertilization. This temporal separation provides a suitable system to study the vertebrate innate immune response in vivo, independently from the adaptive immune response. The transparency of early life stages allows a useful real-time visualization. Adult zebrafish which have complete (innate and adaptative) immune systems offer also advantages over other vertebrate infection models: small size, relatively rapid life cycle, ease of breeding, and a growing list of molecular tools for the study of infectious diseases. In this review, we have tried to give some examples of the potential of zebrafish as a valuable model in innate immunity and inflammation studies.

The challenges of implementing pathogen control strategies for fishes used in biomedical research

Over the past several decades, a number of fish species, including the zebrafish, medaka, and platyfish/swordtail, have become important models for human health and disease. Despite the increasing prevalence of these and other fish species in research, methods for health maintenance and the management of diseases in laboratory populations of these animals are underdeveloped. There is a growing realization that this trend must change, especially as the use of these species expands beyond developmental biology and more towards experimental applications where the presence of underlying disease may affect the physiology animals used in experiments and potentially compromise research results. Therefore, there is a critical need to develop, improve, and implement strategies for managing health and disease in aquatic research facilities. The purpose of this review is to report the proceedings of a workshop entitled "Animal Health and Disease Management in Research Animals" that was recently held at the 5th Aquatic Animal Models for Human Disease in September 2010 at Corvallis, Oregon to discuss the challenges involved with moving the field forward on this front.

Differential expression and intrachromosomal evolution of the sghC1q genes in zebrafish (Danio rerio)

The secreted globular head C1q (sghC1q) genes can be characterized as a family of genetic loci encoding signal peptides followed by single complement component 1q globular (gC1q) motifs. Members of this family have been referred to as precerebellin-like (Cblnl), C1q-like or ovary specific C1q-like factors, and are transcribed in response to infection and/or during early development. This study was primarily undertaken to identify the zebrafish sghC1q (or DrsghC1q) genes that increase their transcription in response to infection and to examine their transcriptional patterns during early development. Twenty sghC1q genes were found in the zebrafish (Danio rerio) genome (Zv9). Two of the examined twenty genes showed significant up-regulation within 24h of infection with the fish pathogen Streptococcus iniae, and eleven of the examined twenty were expressed during early development. Due to the clustered nature of these genes on chromosomes two and seven, intrachromosomal duplication events are hypothesized and explored.

Deep-sequencing analysis of the mouse transcriptome response to infection with Brucella melitensis strains of differing virulence

Brucella melitensis is an important zoonotic pathogen that causes brucellosis, a disease that affects sheep, cattle and occasionally humans. B. melitensis strain M5-90, a live attenuated vaccine cultured from B. melitensis strain M28, has been used as an effective tool in the control of brucellosis in goats and sheep in China. However, the molecular changes leading to attenuated virulence and pathogenicity in B. melitensis remain poorly understood. In this study we employed the Illumina Genome Analyzer platform to perform genome-wide digital gene expression (DGE) analysis of mouse peritoneal macrophage responses to B. melitensis infection. Many parallel changes in gene expression profiles were observed in M28- and M5-90-infected macrophages, suggesting that they employ similar survival strategies, notably the induction of anti-inflammatory and antiapoptotic factors. Moreover, 1019 differentially expressed macrophage transcripts were identified 4 h after infection with the different B. melitensis strains, and these differential transcripts notably identified genes involved in the lysosome and mitogen-activated protein kinase (MAPK) pathways. Further analysis employed gene ontology (GO) analysis: high-enrichment GOs identified endocytosis, inflammatory, apoptosis, and transport pathways. Path-Net and Signal-Net analysis highlighted the MAPK pathway as the key regulatory pathway. Moreover, the key differentially expressed genes of the significant pathways were apoptosis-related. These findings demonstrate previously unrecognized changes in gene transcription that are associated with B. melitensis infection of macrophages, and the central signaling pathways identified here merit further investigation. Our data provide new insights into the molecular attenuation mechanism of strain M5-90 and will facilitate the generation of new attenuated vaccine strains with enhanced efficacy.

Deep sequencing whole transcriptome exploration of the σE regulon in Neisseria meningitidis

Bacteria live in an ever-changing environment and must alter protein expression promptly to adapt to these changes and survive. Specific response genes that are regulated by a subset of alternative σ⁷⁰-like transcription factors have evolved in order to respond to this changing environment. Recently, we have described the existence of a σ^E regulon including the anti-σ-factor MseR in the obligate human bacterial pathogen Neisseria meningitidis. To unravel the complete σ^E regulon in N. meningitidis, we sequenced total RNA transcriptional content of wild type meningococci and compared it with that of mseR mutant cells (ΔmseR) in which σ^E is highly expressed. Eleven coding genes and one non-coding gene were found to be differentially expressed between H44/76 wildtype and H44/76ΔmseR cells. Five of the 6 genes of the σ^E operon, msrA/msrB, and the gene encoding a pepSY-associated TM helix family protein showed enhanced transcription, whilst aniA encoding a nitrite reductase and nspA encoding the vaccine candidate Neisserial surface protein A showed decreased transcription. Analysis of differential expression in IGRs showed enhanced transcription of a non-coding RNA molecule, identifying a σ^E dependent small non-coding RNA. Together this constitutes the first complete exploration of an alternative σ-factor regulon in N. meningitidis. The results direct to a relatively small regulon indicative for a strictly defined response consistent with a relatively stable niche, the human throat, where N. meningitidis resides.

Analysis of transcriptome differences between resistant and susceptible strains of the citrus red mite Panonychus citri (Acari: Tetranychidae)

BACKGROUND

The citrus red mite is a worldwide citrus pest and a common sensitizing allergen of asthma and rhinitis. It has developed strong resistance to many registered acaricides, However, the molecular mechanisms of resistance remain unknown. we therefore used next generation sequencing technology to investigate the global transcriptomes between resistant strains and susceptible strains.

RESULTS

We obtained 34,159, 30,466 and 32,217 unigenes by assembling the SS reads, RS reads and SS&RS reads respectively. There are total 17,581 annotated unigenes from SS&RS reads by BLAST searching databases of nr, the Clusters of Orthologous Groups (COGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) with an E-value ≤ 1e-5, in which 7,075 unigenes were annotated in the COG database, 12, 712 unigenes were found in the KEGG database and 3,812 unigenes were assigned to Gene ontology (GO). Moreover, 2,701 unigenes were judged to be the differentially expressed genes (DEGs) based on the uniquely mapped reads. There are 219 pathways in all annotated unigenes and 198 pathways in DEGs that mapped to the KEGG database. We identified 211 metabolism genes and target genes related to general insecticide resistance such as P450 and Cytochrome b, and further compared their differences between RS and SS. Meanwhile, we identified 105 and 194 genes related to growth and reproduction, respectively, based on the mode of action of Hexythiazox. After further analyses, we found variation in sequences but not in gene expression related to mite growth and reproduction between different strains.

CONCLUSION

To our knowledge, this is the first comparative transcriptome study to discover candidate genes involved in phytophagous mite resistance. This study identified differential unigenes related to general pesticide resistance and organism growth and reproduction in P. citri. The assembled, annotated transcriptomes provide a valuable genomic resource for further understanding the molecular basis of resistance mechanisms.

Digital gene expression profiling of the Phytophthora sojae transcriptome

The transcriptome of the oomycete plant pathogen Phytophthora sojae was profiled at ten different developmental and infection stages based on a 3'-tag digital gene-expression protocol. More than 90 million clean sequence tags were generated and compared with the P. sojae genome and its 19,027 predicted genes. A total of 14,969 genes were detected, of which 10,044 were deemed reliable because they mapped to unambiguous tags. A comparison of the whole-library genes' expression patterns suggested four groups: i) mycelia and zoosporangia, ii) zoospores and cysts, iii) germinating cysts, and iv) five infection site libraries (IF1.5 to IF24h). The libraries from the different groups showed major transitional shifts in gene expression. From the ten libraries, 722 gene expression?pattern clusters were obtained and the top 16 clusters, containing more than half of the genes, comprised enriched genes with different functions including protein localization, triphosphate metabolism, signaling process, and noncoding RNA metabolism. An evaluation of the average expression level of 30 pathogenesis-related gene families revealed that most were infection induced but with diverse expression patterns and levels. A web-based server named the Phytophthora Transcriptional Database has been established.

PAMPs, PRRs and the genomics of gram negative bacterial recognition in fish

Understanding the mechanisms that underpin pathogen recognition and subsequent orchestration of the immune response in fish is an area of significant importance for both basic research and management of health in aquaculture. In recent years much attention has been given to the identification of pattern recognition receptors (PRRs) in fish, however, characterisation of interactions with specific pathogen-associated molecular patterns (PAMPs) is still incomplete. Microarray studies have significantly contributed to functional studies and early descriptions of PAMP-PRR driven activation of specific response cassettes in the genome have been obtained although much is left to be done. In this review we will address gram negative (G-negative) bacterial recognition in fish addressing contributing factors such as structure-function relationships between G-negative PAMPs, current knowledge of fish PRRs and the input achieved by microarray-based studies ranging from in vivo infection studies to directed in vitro PAMP-cell studies. Finally we revisit the endotoxic recognition paradigm in fish and suggest a series of future perspectives that could contribute toward the further elucidation of G-negative bacterial recognition across the highly diverse group of vertebrates that encompass the fishes.

A systematic analysis on DNA methylation and the expression of both mRNA and microRNA in bladder cancer

BACKGROUND

DNA methylation aberration and microRNA (miRNA) deregulation have been observed in many types of cancers. A systematic study of methylome and transcriptome in bladder urothelial carcinoma has never been reported.

METHODOLOGY/PRINCIPAL FINDINGS

The DNA methylation was profiled by modified methylation-specific digital karyotyping (MMSDK) and the expression of mRNAs and miRNAs was analyzed by digital gene expression (DGE) sequencing in tumors and matched normal adjacent tissues obtained from 9 bladder urothelial carcinoma patients. We found that a set of significantly enriched pathways disrupted in bladder urothelial carcinoma primarily related to "neurogenesis" and "cell differentiation" by integrated analysis of -omics data. Furthermore, we identified an intriguing collection of cancer-related genes that were deregulated at the levels of DNA methylation and mRNA expression, and we validated several of these genes (HIC1, SLIT2, RASAL1, and KRT17) by Bisulfite Sequencing PCR and Reverse Transcription qPCR in a panel of 33 bladder cancer samples.

CONCLUSIONS/SIGNIFICANCE

We characterized the profiles between methylome and transcriptome in bladder urothelial carcinoma, identified a set of significantly enriched key pathways, and screened four aberrantly methylated and expressed genes. Conclusively, our findings shed light on a new avenue for basic bladder cancer research.

Deep sequencing of the innate immune transcriptomic response of zebrafish embryos to Salmonella infection

Salmonella enterica serovar Typhimurium (S. typhimurium) bacteria cause an inflammatory and lethal infection in zebrafish embryos. To characterize the embryonic innate host response at the transcriptome level, we have extended and validated previous microarray data by Illumina next-generation sequencing analysis. We obtained 10 million sequence reads from control and Salmonella-infected zebrafish embryos using a tag-based sequencing method (DGE or Tag-Seq) and 15 million reads using whole transcript sequencing (RNA-Seq), which respectively mapped to circa 65% and 85% of 28,716 known Ensembl transcripts. Both sequencing methods showed a strong correlation of sequence read counts per transcript and an overlap of 241 transcripts differentially expressed in response to infection. A lower overlap of 165 transcripts was observed with previous microarray data. Based on the combined sequencing-based and microarray-based transcriptome data we compiled an annotated reference set of infection-responsive genes in zebrafish embryos, encoding transcription factors, signal transduction proteins, cytokines and chemokines, complement factors, proteins involved in apoptosis and proteolysis, proteins with anti-microbial activities, as well as many known or novel proteins not previously linked to the immune response. Furthermore, by comparison of the deep sequencing data of S. typhimurium infection in zebrafish embryos with previous deep sequencing data of Mycobacterium marinum infection in adult zebrafish we derived a common set of infection-responsive genes. This gene set consists of known and putative innate host defense genes that are expressed both in the absence and presence of a fully developed adaptive immune system and that provide a valuable reference for future studies of host-pathogen interactions using zebrafish infection models.

Sequencing-based expression profiling in zebrafish

Gene expression profiling is a powerful technique for studying biological processes, especially tissue/organ-specific ones, at the molecular level. With the rapid development of the next-generation sequencing techniques, high throughput sequencing-based expression profiling techniques have been more and more widely adopted in molecular biology studies. In this chapter, we described a protocol for applying one of the sequencing-based expression profiling techniques, Digital Gene Expression (DGE), for zebrafish research. The protocol provides guidelines for wet-bench experimental procedures as well as for bioinformatics data analyses. We also discuss potential issues/challenges with the use of DGE.

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.

Deep sequencing-based transcriptome analysis of Plutella xylostella larvae parasitized by Diadegma semiclausum

Background

Parasitoid insects manipulate their hosts' physiology by injecting various factors into their host upon parasitization. Transcriptomic approaches provide a powerful approach to study insect host-parasitoid interactions at the molecular level. In order to investigate the effects of parasitization by an ichneumonid wasp (Diadegma semiclausum) on the host (Plutella xylostella), the larval transcriptome profile was analyzed using a short-read deep sequencing method (Illumina). Symbiotic polydnaviruses (PDVs) associated with ichneumonid parasitoids, known as ichnoviruses, play significant roles in host immune suppression and developmental regulation. In the current study, D. semiclausum ichnovirus (DsIV) genes expressed in P. xylostella were identified and their sequences compared with other reported PDVs. Five of these genes encode proteins of unknown identity, that have not previously been reported.

Results

De novo assembly of cDNA sequence data generated 172,660 contigs between 100 and 10000 bp in length; with 35% of > 200 bp in length. Parasitization had significant impacts on expression levels of 928 identified insect host transcripts. Gene ontology data illustrated that the majority of the differentially expressed genes are involved in binding, catalytic activity, and metabolic and cellular processes. In addition, the results show that transcription levels of antimicrobial peptides, such as gloverin, cecropin E and lysozyme, were up-regulated after parasitism. Expression of ichnovirus genes were detected in parasitized larvae with 19 unique sequences identified from five PDV gene families including vankyrin, viral innexin, repeat elements, a cysteine-rich motif, and polar residue rich protein. Vankyrin 1 and repeat element 1 genes showed the highest transcription levels among the DsIV genes.

Conclusion

This study provides detailed information on differential expression of P. xylostella larval genes following parasitization, DsIV genes expressed in the host and also improves our current understanding of this host-parasitoid interaction.

Host-pathogen interactions made transparent with the zebrafish model

The zebrafish holds much promise as a high-throughput drug screening model for immune-related diseases, including inflammatory and infectious diseases and cancer. This is due to the excellent possibilities for in vivo imaging in combination with advanced tools for genomic and large scale mutant analysis. The context of the embryo’s developing immune system makes it possible to study the contribution of different immune cell types to disease progression. Furthermore, due to the temporal separation of innate immunity from adaptive responses, zebrafish embryos and larvae are particularly useful for dissecting the innate host factors involved in pathology. Recent studies have underscored the remarkable similarity of the zebrafish and human immune systems, which is important for biomedical applications. This review is focused on the use of zebrafish as a model for infectious diseases, with emphasis on bacterial pathogens. Following a brief overview of the zebrafish immune system and the tools and methods used to study host-pathogen interactions in zebrafish, we discuss the current knowledge on receptors and downstream signaling components that are involved in the zebrafish embryo’s innate immune response. We summarize recent insights gained from the use of bacterial infection models, particularly the Mycobacterium marinum model, that illustrate the potential of the zebrafish model for high-throughput antimicrobial drug screening.

Comparative transcriptome analyses indicate molecular homology of zebrafish swimbladder and mammalian lung

The fish swimbladder is a unique organ in vertebrate evolution and it functions for regulating buoyancy in most teleost species. It has long been postulated as a homolog of the tetrapod lung, but the molecular evidence is scarce. In order to understand the molecular function of swimbladder as well as its relationship with lungs in tetrapods, transcriptomic analyses of zebrafish swimbladder were carried out by RNA-seq. Gene ontology classification showed that genes in cytoskeleton and endoplasmic reticulum were enriched in the swimbladder. Further analyses depicted gene sets and pathways closely related to cytoskeleton constitution and regulation, cell adhesion, and extracellular matrix. Several prominent transcription factor genes in the swimbladder including hoxc4a, hoxc6a, hoxc8a and foxf1 were identified and their expressions in developing swimbladder during embryogenesis were confirmed. By comparison of enriched transcripts in the swimbladder with those in human and mouse lungs, we established the resemblance of transcriptome of the zebrafish swimbladder and mammalian lungs. Based on the transcriptomic data of zebrafish swimbladder, the predominant functions of swimbladder are in its epithelial and muscular tissues. Our comparative analyses also provide molecular evidence of the relatedness of the fish swimbladder and mammalian lung.

Zebrafish immunoglobulin IgD: unusual exon usage and quantitative expression profiles with IgM and IgZ/T heavy chain isotypes

The zebrafish is an emerging model for comparative immunology and biomedical research. In contrast to the five heavy chain isotype system of mice and human (IgD, IgM, IgA, IgG, IgE), zebrafish harbor gene segments for IgD, IgM, and novel heavy chain isotype called IgZ/T which appears restricted to bony fishes. The purpose of this study was to design and validate a suite of quantitative real time RT-PCR protocols to measure IgH expression in a vertebrate model which has considerable promise for modeling both pathogenic infection and chronic conditions leading to immune dysfunction. Specific primers were designed and following verification of their specificty, relative expression levels of IgD, IgM, and IgZ/T were measured in triplicate for zebrafish raised under standard laboratory conditions. During embryonic stages, low levels of each heavy chain isotype (IgH) were detected with each increasing steadily between 2 and 17 weeks post fertilization. Overall IgM>IgZ>IgD throughout zebrafish development with the copy number of IgM being several fold higher than that of IgD or IgZ/T. IgD exon usage was also characterized, as its extremely long size and presence of a stop codon in the second IgD exon in zebrafish, raised questions as to how this antibody might be expressed. Zebrafish IgD was found to be a chimeric immunoglobulin, with the third IgD exon spliced to the first IgM constant exon thereby circumventing the first and second IgD exons. Collectively, the qRT-PCR results represent the first comparative profile of IgD, IgM, IgZ/T expression over the lifespan of any fish species and the primers and assay parameters reported should prove useful in enabling researchers to rapidly quantify changes in IgH expression in zebrafish models of disease where altered IgH expression is manifested.

Rapid screening of innate immune gene expression in zebrafish using reverse transcription - multiplex ligation-dependent probe amplification

Background

With the zebrafish increasingly being used in immunology and infectious disease research, there is a need for efficient molecular tools to evaluate immune gene expression in this model species. RT-MLPA (reverse transcription - multiplex ligation-dependent probe amplification) provides a sensitive and reproducible method, in which fluorescently labelled amplification products of unique lengths are produced for a defined set of target transcripts. The method employs oligonucleotide probes that anneal to adjacent sites on a target sequence and are then joined by a heat-stable ligase. Subsequently, multiplex PCR with universal primers gives rise to amplicons that can be analyzed with standard sequencing equipment and relative quantification software. Allowing the simultaneous quantification of around 40 selected markers in a one-tube assay, RT-MLPA is highly useful for high-throughput screening applications.

Findings

We employed a dual-colour RT-MLPA probe design for chemical synthesis of probe pairs for 34 genes involved in Toll-like receptor signalling, transcriptional activation of the immune response, cytokine and chemokine production, and antimicrobial defence. In addition, six probe pairs were included for reference genes unaffected by infections in zebrafish. First, we established assay conditions for adult zebrafish infected with different strains of Mycobacterium marinum causing acute and chronic disease. Addition of competitor oligonucleotides was required to achieve peak heights in a similar range for genes with different expression levels. For subsequent analysis of embryonic samples it was necessary to adjust the amounts of competitor oligonucleotides, as the expression levels of several genes differed to a large extent between adult and embryonic tissues. Assay conditions established for one-day-old Salmonella typhimurium-infected embryos could be transferred without further adjustment to five-day-old M. marinum-infected larvae. RT-MLPA results were compared with results of previous transcriptome analyses and with real-time PCR data, demonstrating a good correlation between all expression analysis methods.

Conclusions

The RT-MLPA assay developed in this study provides a rapid, cheap, and robust analysis tool for simultaneous quantification of a set of 34 innate immune response genes. With adjustment of conditions, the assay is suitable for infection studies in both adult and embryonic zebrafish. Application of RT-MLPA will facilitate high-throughput screening of immune responses in the zebrafish model.

Transcriptome characterization of the South African abalone Haliotis midae using sequencing-by-synthesis

BACKGROUND

Worldwide, the genus Haliotis is represented by 56 extant species and several of these are commercially cultured. Among the six abalone species found in South Africa, Haliotis midae is the only aquacultured species. Despite its economic importance, genomic sequence resources for H. midae, and for abalone in general, are still scarce. Next generation sequencing technologies provide a fast and efficient tool to generate large sequence collections that can be used to characterize the transcriptome and identify expressed genes associated with economically important traits like growth and disease resistance.

RESULTS

More than 25 million short reads generated by the Illumina Genome Analyzer were de novo assembled in 22,761 contigs with an average size of 260 bp. With a stringent E-value threshold of 10-10, 3,841 contigs (16.8%) had a BLAST homologous match against the Genbank non-redundant (NR) protein database. Most of these sequences were annotated using the gene ontology (GO) and eukaryotic orthologous groups of proteins (KOG) databases and assigned to various functional categories. According to annotation results, many gene families involved in immune response were identified. Thousands of simple sequence repeats (SSR) and single nucleotide polymorphisms (SNP) were detected. Setting stringent parameters to ensure a high probability of amplification, 420 primer pairs in 181 contigs containing SSR loci were designed.

CONCLUSION

This data represents the most comprehensive genomic resource for the South African abalone H. midae to date. The amount of assembled sequences demonstrated the utility of the Illumina sequencing technology in the transcriptome characterization of a non-model species. It allowed the development of several markers and the identification of promising candidate genes for future studies on population and functional genomics in H. midae and in other abalone species.