Digital Gene Expression Profiling Research Articles

A peptide hormone gene, GhPSK promotes fibre elongation and contributes to longer and finer cotton fibre

Cotton fibres, the single-celled trichomes derived from the ovule epidermis, provide the most important natural material for the global textile industry. A number of studies have demonstrated that regulating endogenous hormone levels through transgenic approaches can improve cotton fibre qualities. Phytosulfokine-α (PSK-α) is a novel peptide hormone in plants that is involved in regulating cell proliferation and elongation. However, its potential applications in crop genetic improvement have not been evaluated. In this study, we describe how exogenous PSK-α application promotes cotton fibre cell elongation in vitro. Chlorate, an effective inhibitor of peptide sulfation, suppressed fibre elongation in ovule culture. Exogenously applied PSK-α partly restored the chlorate-induced suppression. A putative PSK gene (GhPSK) was cloned from Gossypium hirsutum. Expression pattern analysis revealed that GhPSK is preferentially expressed in rapidly elongating fibre cells (5-20 days postanthesis). Overexpression of GhPSK in cotton increased the endogenous PSK-α level and promoted cotton fibre cell elongation, resulting in longer and finer fibres. Further results from electrophysiological and physiological analyses suggest that GhPSK affects fibre development through regulation of K(+) efflux. Digital gene expression (DGE) profile analysis of GhPSK overexpression lines indicates that PSK signalling may regulate the respiratory electron-transport chain and reactive oxygen species to affect cotton fibre development. These results imply that peptide hormones are involved in cotton fibre growth and suggest a new strategy for the biotechnological improvement of cotton fibre quality.

Analysis options for high-throughput sequencing in miRNA expression profiling

BackgroundRecently high-throughput sequencing (HTS) using next generation sequencing techniques became useful in digital gene expression profiling.Our study introduces analysis options for HTS data based on mapping to miRBase or counting and grouping of identical sequence reads. Those approaches allow a hypothesis free detection of miRNA differential expression.MethodsWe compare our results to microarray and qPCR data from one set of RNA samples. We use Illumina platforms for microarray analysis and miRNA sequencing of 20 samples from benign follicular thyroid adenoma and malignant follicular thyroid carcinoma. Furthermore, we use three strategies for HTS data analysis to evaluate miRNA biomarkers for malignant versus benign follicular thyroid tumors.ResultsHigh correlation of qPCR and HTS data was observed for the proposed analysis methods. However, qPCR is limited in the differential detection of miRNA isoforms. Moreover, we illustrate a much broader dynamic range of HTS compared to microarrays for small RNA studies. Finally, our data confirm hsa-miR-197-3p, hsa-miR-221-3p, hsa-miR-222-3p and both hsa-miR-144-3p and hsa-miR-144-5p as potential follicular thyroid cancer biomarkers.ConclusionsCompared to microarrays HTS provides a global profile of miRNA expression with higher specificity and in more detail. Summarizing of HTS reads as isoform groups (analysis pipeline B) or according to functional criteria (seed analysis pipeline C), which better correlates to results of qPCR are promising new options for HTS analysis. Finally, data opens future miRNA research perspectives for HTS and indicates that qPCR might be limited in validating HTS data in detail.

Transcriptome Analysis of Silkworm, Bombyx mori, during Early Response to Beauveria bassiana Challenges

Host–pathogen interactions are complex processes and it is a central challenge to reveal these interactions. Fungal infection of silkworm, Bombyx mori, may induce a variety of responsive reaction. However, little is known about the molecular mechanism of silkworm immune response against the fungal infection. To obtain an overview of the interaction between silkworm and an entomopathogenic fungus Beauveria bassiana, Digital Gene Expression profiling, a tag based high-throughput transcriptome sequencing method, was employed to screen and identify differentially expressed genes (DEGs, FDR≤0.001, ∣log2ratio∣≥1) of silkworm larvae during early response against B. bassiana infection. Total 1430 DEGs including 960 up-regulated and 470 down-regulated ones were identified, of which 627 DEGs can be classified into GO categories by Gene Ontology (GO) analysis. KEGG pathways analysis of these DEGs suggested that many biological processes, such as defense and response, signal transduction, phagocytosis, regulation of gene expression, RNA splicing, biosynthesis and metabolism, protein transport etc. were involved in the interaction between the silkworm and B. bassiana. A number of differentially expressed fungal genes were also identified by mapping the sequencing tags to B. bassiana genome. These results provided new insights to the molecular mechanism of silkworm immune response to B. bassiana infection.

Why mosaic? Gene expression profiling of African cassava mosaic virus‐infected cassava reveals the effect of chlorophyll degradation on symptom development

Cassava mosaic disease, caused by cassava begomoviruses, is the most serious disease for cassava in Africa. However, the pathogenesis of this disease is poorly understood. We employed high throughput digital gene expression profiling based on the Illumina Solexa sequencing technology to investigate the global transcriptional response of cassava to African cassava mosaic virus infection. We found that 3,210 genes were differentially expressed in virus-infected cassava leaves. Gene ontology term and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that genes implicated in photosynthesis were most affected, consistent with the chlorotic symptoms observed in infected leaves. The upregulation of chlorophyll degradation genes, including the genes encoding chlorophyllase, pheophytinase, and pheophorbide a oxygenase, and downregulation of genes encoding the major apoproteins in light-harvesting complex II were confirmed by qRT-PCR. These findings, together with the reduction of chlorophyll b content and fewer grana stacks in the infected leaf cells, reveal that the degradation of chlorophyll plays an important role in African cassava mosaic virus symptom development. This study will provide a road map for future investigations into viral pathogenesis.

Identification and bioinformatic analysis of signal responsive/calmodulin-binding transcription activators gene models in Vitis vinifera

In this study, 10 grapevine (Vitis vinifera) SR/CAMTA (Signal Responsive/Calmodulin-binding Transcription Activators) gene models were identified from three grapevine genome protein datasets. They belong to four gene groups: VvCAMTA1, VvCAMTA3, VvCAMTA4 and VvCAMTA5, which were located on chromosome 5, 7_random, 1 and 5, respectively. Alternative splicing could explain the multiple gene models in one gene group. Subcellular localization using the WoLF tool showed that most of the VvCAMTAs were located in the nucleus, except for VvCAMTA3.1, VvCAMTA3.2 and VvCAMTA5.2, which were located in the chloroplast, chloroplast and cytosol, respectively. Subcellular localization using TargetP showed that most of the VvCAMTAs were not located in the chloroplast, mitochondrion and secretory pathway in cells. VvCAMTA1.1 and VvCAMTA1.2 were located in the mitochondria. The digital gene expression profile showed that VvCAMTAs play important roles in Ca2+ signal transduction. The gene expression patterns of VvCAMTAs were different; for example, VvCAMTA1 was expressed mainly in the bud, while VvCAMTA3 was expressed mainly in fruit and inflorescence, with low expression in the bud. The results of this study make a substantial contribution to our knowledge concerning genes, genome annotation, and cell signal transduction in grapevine.

Transcriptome and Biochemical Analyses Revealed a Detailed Proanthocyanidin Biosynthesis Pathway in Brown Cotton Fiber

Brown cotton fiber is the major raw material for colored cotton industry. Previous studies have showed that the brown pigments in cotton fiber belong to proanthocyanidins (PAs). To clarify the details of PA biosynthesis pathway in brown cotton fiber, gene expression profiles in developing brown and white fibers were compared via digital gene expression profiling and qRT-PCR. Compared to white cotton fiber, all steps from phenylalanine to PA monomers (flavan-3-ols) were significantly up-regulated in brown fiber. Liquid chromatography mass spectrometry analyses showed that most of free flavan-3-ols in brown fiber were in 2, 3-trans form (gallocatechin and catechin), and the main units of polymeric PAs were trihydroxylated on B ring. Consistent with monomeric composition, the transcript levels of flavonoid 3′, 5′-hydroxylase and leucoanthocyanidin reductase in cotton fiber were much higher than their competing enzymes acting on the same substrates (dihydroflavonol 4-reductase and anthocyanidin synthase, respectively). Taken together, our data revealed a detailed PA biosynthesis pathway wholly activated in brown cotton fiber, and demonstrated that flavonoid 3′, 5′-hydroxylase and leucoanthocyanidin reductase represented the primary flow of PA biosynthesis in cotton fiber.

Isolation and expression studies of the ERD15 gene involved in drought-stressed responses.

The early response to the dehydration 15 (ERD15) gene is widely involved in the processes of signal transduction, programmed cell death, gene transcription, and stress tolerance in plants. In a previous study, the ERD15 gene was shown to be an important regulator of the abscisic acid response and salicylic acid-dependent defense pathway, acting as an important negative regulator of abscisic acid. The complete IbERD15 gene (accession No. KF723428) was isolated by reverse transcription-polymerase chain reaction. The IbERD15 gene contains an open reading frame of 504 bp, encodes a peptide of 167 amino acids, and has a molecular mass of 18.725 kDa. The transcript levels of the IbERD15 gene in a variety of tissues were examined by digital gene expression profiling. The roots of the sweet potato were treated by 3 degrees of polyethylene glycol, and the results indicate that the IbERD15 gene might play an important role in the defense response to drought stress. Moreover, the IbERD15 gene was successfully transformed into yeast cells for analysis of drought tolerance in transgenic yeast.

Transcriptome analysis of blood orange (Citrus sinensis) following fruit bagging treatment by digital gene expression profiling

SummaryBagging has been widely used in fruit crops to improve fruit quality. In this study, the internal qualities of blood orange (Citrus sinensis cv. ‘Tarocco’ ) fruit were significantly affected by bagging treatment. Total flavonoid (TFL) and total anthocyanin (TAN) concentrations increased, while total soluble solids (TSS), total phenolics (TPH), and ascorbic acid (AsA) concentrations and titratable acidity (TA) decreased in response to bagging treatment. To elucidate the molecular mode of action of bagging on fruit, high-throughput tag-sequencing (Tag-seq) analysis was used to analyse gene expression profiles in blood orange fruit in response to a bagging treatment. Over 21 × 106 total clean reads were obtained per library. Approx. 53.7 – 71.7% and 3.1 – 6.4% of the clean reads were mapped onto citrus genes or intergenic regions, respectively. However, approx. 25.2 – 39.9% of the clean reads failed to align with the citrus genome. Overall, bagging treatment resulted in an increase in transcripts involved in a range of metabolic pathways such as hormone synthesis and signalling, transcriptional regulation processes, the tricarboxylic acid (TCA) cycle, sucrose and starch metabolism, ascorbate metabolism, and the phenylpropanoid pathway, many of which are related to fruit quality. Analysis of the differentially expressed genes revealed that decreased TSS concentrations, TA, and AsA concentrations could mainly be attributed to the increased expression of genes involved in catabolic pathways, rather than in anabolic pathways. Competition for limited amounts of substrates for the flavonoid, phenolics, and anthocyanin pathways may have led to an increase in TAN and TFL concentrations, and a decrease in TPH concentration. In addition, our results showed that bagging treatment imposed a variety of stresses on fruit, thereby leading to an array of biochemical, physiological, and molecular changes. One gene (bHLH) that may be involved in anthocyanin biosynthesis in blood orange fruit has been identified and is discussed. These results reveal the extensive biochemical and genetic mechanisms operating in citrus fruit in response to a bagging treatment.

Novel expression profiles of microRNAs suggest that specific miRNAs regulate gene expression for the sexual maturation of female Schistosoma japonicum after pairing

BackgroundSchistosoma japonicum is one of the major causative agents of schistosomiasis. The pairing of males and females leads to female sexual maturation and maintains this mature state. However, the mechanisms by which pairing facilitates sexual maturation are yet to be investigated.MethodsParasites isolated from single- and double-sex cercariae-infected mice were analyzed by Solexa to uncover pair-regulated miRNA profiles. To reveal the biological functions of differentially expressed miRNAs among the samples, we predicted the target genes of these differentially expressed miRNAs and compared the gene expression between 23-d-old female schistosomula from double-sex infections (23DSI) and 23-d-old female schistosomula from single-sex infections (23SSI) by analyzing digital gene expression profiling (DGE). KEGG pathway analysis was used to investigate the relevant biological processes of these target genes to understand the significance of differentially expressed miRNAs after pairing.ResultsThe differentially expressed miRNA profiles of female 18- and 23-d post-single- and double-sex infections were analysed by Solexa. Similar miRNA profiles were observed in 18SSI and 18DSI, with the presence of identically expressed high-abundance miRNA, such as miRNA-1, miRNA-71b-5p and let-7. By contrast, in 23DSI and 23SSI, most of these high-abundance miRNAs were down-regulated. Furthermore, among all samples, bantam was distinctly up-regulated in 23 DSI, and miR-1, miR-71, miR-7-5p, and miR-7 were distinctly up-regulated in 23SSI. The transcriptomes of 23DSI and 23SSI revealed that the predicted target genes of miRNA-1, miRNA-71, miRNA-7, and miR-7-5p were associated with the ribonucleoprotein complex assembly and microtubule-based process. Conversely, the predicted target genes of bantam were related to the embryo development, development of primary sexual characteristics and regulation of transcription. KEGG pathway analysis revealed that in unpaired females, the highly-expressed miRNA-1, miRNA-71, miRNA-7, and miR-7-5p only inhibited the limited pathways, such as proteasome and ribosome assembly. Meanwhile, in paired mature females, highly-expressed bantam inhibited more biological pathways, such as the citrate cycle, glycolysis, fatty acid biosynthesis and RNA degradation.ConclusionsThe differentially expressed miRNAs between 23SSI and 23DSI and their different functions indicated that more genes or metabolic pathways in paired mature females were inhibited than those in unpaired ones. The results suggested that after pairing, specific miRNAs regulated gene expression to lead to female sexual maturation.

Analysis of banana transcriptome and global gene expression profiles in banana roots in response to infection by race 1 and tropical race 4 of Fusarium oxysporum f. sp. cubense.

BackgroundCavendish, the most widely grown banana cultivar, is relatively resistant to Race 1 of Fusarium oxysporum f. sp. cubense (Foc1) which caused widespread Panama disease during the first half of the 20th century but is susceptible to Tropical Race 4 of Foc (Foc TR4) which is threatening world banana production. The genome of the diploid species Musa acuminata which is the ancestor of a majority of triploid banana cultivars has recently been sequenced. Availability of banana transcriptomes will be highly useful for improving banana genome annotation and for biological research. The knowledge of global gene expression patterns influenced by infection of different Foc races will help to understand the host responses to the infection.ResultsRNA samples from different organs of the Cavendish cultivar were pooled for deep sequencing using the Illumina technology. Analysis of the banana transcriptome led to identification of over 842 genes that were not annotated by the Musa genome project. A large number of simple nucleotide polymorphisms (SNPs) and short insertions and deletion (indels) were identified from the transcriptome data. GFP-expressing Foc1 and Foc TR4 were used to monitor the infection process. Both Foc1 and Foc TR4 were found to be able to invade banana roots and spread to root vascular tissues in the first two days following inoculation. Digital gene expression (DGE) profiling analysis reveal that the infection by Foc1 and Foc TR4 caused very similar changes in the global gene expression profiles in the banana roots during the first two days of infection. The Foc infection led to induction of many well-known defense-related genes. Two genes encoding the ethylene biosynthetic enzyme ACC oxidase and several ethylene-responsive transcription factors (ERF) were among the strongly induced genes by both Foc1 and Foc TR4.ConclusionsBoth Foc1 and Foc TR4 are able to spread into the vascular system of banana roots during the early infection process and their infection led to similar gene expression profiles in banana roots. The transcriptome profiling analysis indicates that the ethylene synthetic and signalling pathways were activated in response to the Foc infection.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-14-851) contains supplementary material, which is available to authorized users.

Phylogenetic and transcriptional analysis of an expanded bZIP transcription factor family in Phytophthora sojae

BackgroundBasic leucine zipper (bZIP) transcription factors are present exclusively in eukaryotes and constitute one of the largest and most diverse transcription factor families. The proteins are responsible for central developmental and physiological processes in plants, animals, and fungi, including the pathogenicity of fungal plant pathogens. However, there is limited understanding of bZIPs in oomycetes, which are fungus-like organisms in the kingdom Stramenopila. Oomycetes include many destructive plant pathogens, including the well-studied species Phytophthora sojae, which causes soybean stem and root rot.ResultsCandidate bZIPs encoded in the genomes of P. sojae and four other oomycetes, two diatoms, and two fungal species were predicted using bioinformatic methods. Comparative analysis revealed expanded numbers of bZIP candidates in oomycetes, especially the Phytophthora species, due to the expansion of several novel bZIP classes whose highly conserved asparagines in basic DNA-binding regions were substituted by other residues such as cysteine. The majority of these novel bZIP classes were mostly restricted to oomycetes. The large number of novel bZIPs appears to be the result of widespread gene duplications during oomycete evolution. The majority of P. sojae bZIP candidates, including both conventional and novel bZIP classes, were predicted to contain canonical protein secondary structures. Detection of gene transcripts using digital gene expression profiling and qRT-PCR suggested that most of the candidates were not pseudogenes. The major transcriptional shifts of bZIPs occurred during the zoosporangia/zoospore/cyst and host infection stages. Several infection-associated bZIP genes were identified that were positively regulated by H2O2 exposure.ConclusionsThe identification of large classes of bZIP proteins in oomycetes with novel bZIP motif variants, that are conserved and developmentally regulated and thus presumably functional, extends our knowledge of this important family of eukaryotic transcription factors. It also lays the foundation for detailed studies of the roles of these proteins in development and infection in P. sojae and other oomycetes.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-14-839) contains supplementary material, which is available to authorized users.

De novo characterization of the Anthuriumtranscriptome and analysis of its digital gene expression under cold stress

BackgroundAnthurium andraeanum is one of the most popular tropical flowers. In temperate and cold zones, a much greater risk of cold stress occurs in the supply of Anthurium plants. Unlike the freeze-tolerant model plants, Anthurium plants are particularly sensitive to low temperatures. Improvement of chilling tolerance in Anthurium may significantly increase its production and extend its shelf-life. To date, no previous genomic information has been reported in Anthurium plants.ResultsUsing Illumina sequencing technology, we generated over two billion base of high-quality sequence in Anthurium, and demonstrated de novo assembly and annotation of genes without prior genome information. These reads were assembled into 44,382 unigenes (mean length = 560 bp). Based on similarity search with known protein in the non-redundant (nr) protein database, 27396 unigenes (62%) were functionally annotated with a cut-off E-value of 10-5. Further, DGE tags were mapped to the assembled transcriptome for gene expression analysis under cold stress. In total, 4363 differentially expressed genes were identified. Among these genes, 292, 805 and 708 genes were up-regulated after 1-h, 5-h and 24-h cold treatment, respectively. Then we mapped these cold-induced genes to the KEGG database. Specific enrichment was observed in photosynthesis pathway, metabolic pathways and oxidative phosphorylation pathway in 1-h cold-treated plants. After a 5-h cold treatment, the metabolic pathways and oxidative phosphorylation pathway were significantly identified as the top two pathways. After 24-h cold treatment, mRNA surveillance pathway, RNA transport pathway and plant-pathogen interaction pathway were significantly enriched. Together, a total of 39 cold-inducible transcription factors were identified, including subsets of AP2/ERF, Zinc figure, NAC, MYB and bZIP family members.ConclusionOur study is the first to provide the transcriptome sequence resource for Anthurium plants, and demonstrate its digital gene expression profiling under cold conditions using the assembled transcriptome data for reference. These data provides a valuable resource for genetic and genomic studies under abiotic conditions for Anthurium plants.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-14-827) contains supplementary material, which is available to authorized users.

Transcriptome profiling of fruit development and maturation in Chinese white pear (Pyrus bretschneideri Rehd)

BackgroundPear (Pyrus spp) is an important fruit species worldwide; however, its genetics and genomic information is limited. Combining the Solexa/Illumina RNA-seq high-throughput sequencing approach (RNA-seq) with Digital Gene Expression (DGE) analysis would be a powerful tool for transcriptomic study. This paper reports the transcriptome profiling analysis of Chinese white pear (P. bretschneideri) using RNA-seq and DGE to better understand the molecular mechanisms in fruit development and maturation of Chinese white pear.ResultsDe novo transcriptome assembly and gene expression analysis of Chinese white pear were performed in an unprecedented depth (5.47 gigabase pairs) using high-throughput Illumina RNA-seq combined with a tag-based Digital Gene Expression (DGE) system. Approximately, 60.77 million reads were sequenced, trimmed, and assembled into 90,227 unigenes. These unigenes comprised 17,619 contigs and 72,608 singletons with an average length of 508 bp and had an N50 of 635 bp. Sequence similarity analyses against six public databases (Uniprot, NR, and COGs at NCBI, Pfam, InterPro, and KEGG) found that 61,636 unigenes can be annotated with gene descriptions, conserved protein domains, or gene ontology terms. By BLASTing all 61,636 unigenes in KEGG, a total of 31,215 unigenes were annotated into 121 known metabolic or signaling pathways in which a few primary, intermediate, and secondary metabolic pathways are directly related to pear fruit quality. DGE libraries were constructed for each of the five fruit developmental stages. Variations in gene expression among all developmental stages of pear fruit were significantly different in a large amount of unigenes.ConclusionExtensive transcriptome and DGE profiling data at five fruit developmental stages of Chinese white pear have been obtained from a deep sequencing, which provides comprehensive gene expression information at the transcriptional level. This could facilitate understanding of the molecular mechanisms in fruit development and maturation. Such a database can also be used as a public information platform for research on molecular biology and functional genomics in pear and other related species.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-14-823) contains supplementary material, which is available to authorized users.

Genome-wide survey and expression profiles of the AP2/ERF family in castor bean (Ricinus communis L.)

BackgroundThe AP2/ERF transcription factor, one of the largest gene families in plants, plays a crucial role in the regulation of growth and development, metabolism, and responses to biotic and abiotic stresses. Castor bean (Ricinus communis L., Euphobiaceae) is one of most important non-edible oilseed crops and its seed oil is broadly used for industrial applications. The available genome provides a great chance to identify and characterize the global information on AP2/ERF transcription factors in castor bean, which might provide insights in understanding the molecular basis of the AP2/ERF family in castor bean.ResultsA total of 114 AP2/ERF transcription factors were identified based on the genome in castor bean. According to the number of the AP2/ERF domain, the conserved amino acid residues within AP2/ERF domain, the conserved motifs and gene organization in structure, and phylogenetical analysis, the identified 114 AP2/ERF transcription factors were characterized. Global expression profiles among different tissues using high-throughput sequencing of digital gene expression profiles (DGEs) displayed diverse expression patterns that may provide basic information in understanding the function of the AP2/ERF gene family in castor bean.ConclusionsThe current study is the first report on identification and characterization of the AP2/ERF transcription factors based on the genome of castor bean in the family Euphobiaceae. Results obtained from this study provide valuable information in understanding the molecular basis of the AP2/ERF family in castor bean.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-14-785) contains supplementary material, which is available to authorized users.

Genomic surveys and expression analysis of bZIP gene family in castor bean (Ricinus communis L.)

The basic leucine zipper (bZIP) transcription factors comprise a family of transcriptional regulators present extensively in plants, involved in regulating diverse biological processes such as flower and vascular development, seed maturation, stress signaling and pathogen defense. Castor bean (Ricinus communis L. Euphorbiaceae) is one of the most important non-edible oilseed crops and its seed oil is broadly used for industrial applications. We performed a comprehensive genome-wide identification and analysis of the bZIP transcription factors that exist in the castor bean genome in this study. In total, 49 RcbZIP transcription factors were identified, characterized and categorized into 11 groups (I-XI) based on their gene structure, DNA-binding sites, conserved motifs, and phylogenetic relationships. The dimerization properties of 49 RcbZIP proteins were predicted on the basis of the characteristic features in the leucine zipper. Global expression profiles of 49 RcbZIP genes among different tissues were examined using high-throughput sequencing of digital gene expression profiles, and resulted in diverse expression patterns that may provide basic information to further reveal the function of the 49 RcbZIP genes in castor bean. The results obtained from this study would provide valuable information in understanding the molecular basis of the RcbZIP transcription factor family and their potential function in regulating the growth and development, particularly in seed filling of castor bean.

A Comparison of Digital Gene Expression Profiling and Methyl DNA Immunoprecipitation as Methods for Gene Discovery in Honeybee (Apis mellifera) Behavioural Genomic Analyses

The honey bee has a well-organized system of division of labour among workers. Workers typically progress through a series of discrete behavioural castes as they age, and this has become an important case study for exploring how dynamic changes in gene expression can influence behaviour. Here we applied both digital gene expression analysis and methyl DNA immunoprecipitation analysis to nurse, forager and reverted nurse bees (nurses that have returned to the nursing state after a period spent foraging) from the same colony in order to compare the outcomes of these different forms of genomic analysis. A total of 874 and 710 significantly differentially expressed genes were identified in forager/nurse and reverted nurse/forager comparisons respectively. Of these, 229 genes exhibited reversed directions of gene expression differences between the forager/nurse and reverted nurse/forager comparisons. Using methyl-DNA immunoprecipitation combined with high-throughput sequencing (MeDIP-seq) we identified 366 and 442 significantly differentially methylated genes in forager/nurse and reverted nurse/forager comparisons respectively. Of these, 165 genes were identified as differentially methylated in both comparisons. However, very few genes were identified as both differentially expressed and differentially methylated in our comparisons of nurses and foragers. These findings confirm that changes in both gene expression and DNA methylation are involved in the nurse and forager behavioural castes, but the different analytical methods reveal quite distinct sets of candidate genes.

Insights into food preference in hybrid F1 of Siniperca chuatsi (♀) × Siniperca scherzeri (♂) mandarin fish through transcriptome analysis

BackgroundAs economically relevant traits, feeding behavior and food preference domestication determine production cost and profitability. Although there are intensive research efforts on feeding behavior and food intake, little is known about food preference. Mandarin fish accept only live prey fish and refuse dead prey fish or artificial diets. Very little is currently known about the genes regulating this unique food preference.ResultsUsing transcriptome sequencing and digital gene expression profiling, we identified 1,986 and 4,526 differentially expressed genes in feeders and nonfeeders of dead prey fish, respectively. Up-regulation of Crbp, Rgr and Rdh8, and down-regulation of Gc expression, consistent with greater visual ability in feeders, could promote positive phototaxis. Altered expressions of period, casein kinase and Rev-erbα might reset circadian phase. Down-regulation of orexigenic and up-regulation of anorexigenic genes in feeders were associated with lower appetite. The mRNA levels of Creb, c-fos, C/EBP, zif268, Bdnf and Syt were dramatically decreased in feeders, which might result in significant deficiency in memory retention of its natural food preference (live prey fish). There were roughly 100 times more potential SNPs in feeders than in nonfeeders.ConclusionsIn summary, differential expression in the genes identified shed new light on why mandarin fish only feed on live prey fish, with pathways regulating retinal photosensitivity, circadian rhythm, appetite control, learning and memory involved. We also found dramatic difference in SNP abundance in feeders vs nonfeeders. These differences together might account for the different food preferences. Elucidating the genes regulating the unique food preference (live prey fish) in mandarin fish could lead to a better understanding of mechanisms controlling food preference in animals, including mammals.

De novo characterization of Larix gmelinii (Rupr.) Rupr. transcriptome and analysis of its gene expression induced by jasmonates

BackgroundLarix gmelinii is a dominant tree species in China’s boreal forests and plays an important role in the coniferous ecosystem. It is also one of the most economically important tree species in the Chinese timber industry due to excellent water resistance and anti-corrosion of its wood products. Unfortunately, in Northeast China, L. gmelinii often suffers from serious attacks by diseases and insects. The application of exogenous volatile semiochemicals may induce and enhance its resistance against insect or disease attacks; however, little is known regarding the genes and molecular mechanisms related to induced resistance.ResultsWe performed de novo sequencing and assembly of the L. gmelinii transcriptome using a short read sequencing technology (Illumina). Chemical defenses of L. gmelinii seedlings were induced with jasmonic acid (JA) or methyl jasmonate (MeJA) for 6 hours. Transcriptomes were compared between seedlings induced by JA, MeJA and untreated controls using a tag-based digital gene expression profiling system. In a single run, 25,977,782 short reads were produced and 51,157 unigenes were obtained with a mean length of 517 nt. We sequenced 3 digital gene expression libraries and generated between 3.5 and 5.9 million raw tags, and obtained 52,040 reliable reference genes after removing redundancy. The expression of disease/insect-resistance genes (e.g., phenylalanine ammonialyase, coumarate 3-hydroxylase, lipoxygenase, allene oxide synthase and allene oxide cyclase) was up-regulated. The expression profiles of some abundant genes under different elicitor treatment were studied by using real-time qRT-PCR.The results showed that the expression levels of disease/insect-resistance genes in the seedling samples induced by JA and MeJA were higher than those in the control group. The seedlings induced with MeJA elicited the strongest increases in disease/insect-resistance genes.ConclusionsBoth JA and MeJA induced seedlings of L. gmelinii showed significantly increased expression of disease/insect-resistance genes. MeJA seemed to have a stronger induction effect than JA on expression of disease/insect-resistance related genes. This study provides sequence resources for L. gmelinii research and will help us to better understand the functions of disease/insect-resistance genes and the molecular mechanisms of secondary metabolisms in L. gmelinii.

Induced Expression of AtLEC1 and AtLEC2 Differentially Promotes Somatic Embryogenesis in Transgenic Tobacco Plants

Arabidopsis LEAFY COTYLEDON (LEC) genes, AtLEC1 and AtLEC2, are important embryonic regulators that play key roles in morphogenesis and maturation phases during embryo development. Ectopic expression of AtLEC1 and AtLEC2 in tobacco caused abnormality in transgenic seedling. When transgenic seeds germinated on medium containing 30 µM DEX, LEC1 transgenic seedlings were ivory and fleshy, with unexpanded cotyledons, stubby hypocotyls, short roots and no obvious callus formation at the shoot meristem position. While LEC2 transgenic seedlings formed embryonic callus on the shoot apical meristem and somatic embryo-like structures emerged from the surface of the callus. When callus were transferred to hormone free MS0 medium more shoots were regenerated from each callus. However, shoot formation was not observed in LEC1 overexpressors. To investigate the mechanisms of LEC2 in somatic embryogenesis, we studied global gene expression by digital gene expression profiling analysis. The results indicated that ectopic expression of LEC2 genes induced accumulation of embryo-specific proteins such as seed storage proteins, late embryogenesis abundant (LEA) proteins, fatty acid biosynthetic enzymes, products of steroid biosynthesis related genes and key regulatory genes of the embryo development. Genes of plant-specific transcription factors such as NAC domain protein, AP2 and GRAS family, resistance-related as well as salicylic acid signaling related genes were up-regulated in LEC2 transgenic seedlings. Ectopi c expression of LEC2 induced large number of somatic embryo formation and shoot regeneration but 20 d DEX induction of LEC1 is not sufficient to induce somatic embryogenesis and shoot formation. Our data provide new information to understand the mechanisms on LEC2 gene’s induction of somatic embryogenesis.

Global transcriptome analysis of orange wheat blossom midge, Sitodiplosis mosellana (Gehin) (Diptera: Cecidomyiidae) to identify candidate transcripts regulating diapause.

BackgroundMany insects enter a developmental arrest (diapause) that allows them to survive harsh seasonal conditions. Despite the well-established ecological significance of diapause, the molecular basis of this crucial adaptation remains largely unresolved. Sitodiplosis mosellana (Gehin), the orange wheat blossom midge (OWBM), causes serious damage to wheat throughout the northern hemisphere, and sporadic outbreaks occur in the world. Traits related to diapause appear to be important factors contributing to their rapid spread and outbreak. To better understand the diapause mechanisms of OWBM, we sequenced the transcriptome and determined the gene expression profile of this species.Methodology/Principal FindingsIn this study, we performed de novo transcriptome analysis using short-read sequencing technology (Illumina) and gene expression analysis with a tag-based digital gene expression (DGE) system. The sequencing results generated 89,117 contigs, and 45,713 unigenes. These unigenes were annotated by Blastx alignment against the NCBI non-redundant (nr), Clusters of orthologous groups (COG), gene orthology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. 20,802 unigenes (45.5% of the total) matched with protein in the NCBI nr database. Two digital gene expression (DGE) libraries were constructed to determine differences in gene expression profiles during diapause and non-diapause developmental stages. Genes related to diapause were analyzed in detail and in addition, nine diapause-related genes were analyzed by real time PCR.Conclusions/SignificanceThe OWBM transcriptome greatly improves our genetic understanding and provides a platform for functional genomics research of this species. The DGE profiling data provides comprehensive information at the transcriptional level that facilitates our understanding of the molecular mechanisms of various physiological aspects including development and diapause stages in OWBM. From this study it is evident that various genes coding metabolic enzymes are crucial for diapause and metamorphosis.