A new RNASeq-based reference transcriptome for sugar beet and its application in transcriptome-scale analysis of vernalization and gibberellin responses

Leaf and shoot apical meristem transcriptomes of quinoa (Chenopodium quinoa Willd.) in response to photoperiod and plant development

Understanding the regulation of flowering time is crucial for adaptation of crops to new environment. In this study, we examined the timing of floral transition and analysed transcriptomes in leaf and shoot apical meristems of photoperiod-sensitive and -insensitive quinoa accessions. Histological analysis showed that floral transition in quinoa initiates 2-3 weeks after sowing. We found four groups of differentially expressed genes in quinoa genome that responded to plant development and floral transition: (i) 222 genes responsive to photoperiod in leaves, (ii) 1812 genes differentially expressed between accessions under long-day conditions in leaves, (iii) 57 genes responding to developmental changes under short-day conditions in leaves and (iv) 911 genes responding to floral transition within the shoot apical meristem. Interestingly, among numerous candidate genes, two putative FT orthologs together with other genes (e.g. SOC1, COL, AP1) were previously reported as key regulators of flowering time in other species. Additionally, we used coexpression networks to associate novel transcripts to a putative biological process based on the annotated genes within the same coexpression cluster. The candidate genes in this study would benefit quinoa breeding by identifying and integrating their beneficial haplotypes in crossing programs to develop adapted cultivars to diverse environmental conditions.

Identification of the Genes Encoding B3 Domain-Containing Proteins Related to Vernalization of Beta vulgaris

Vernalization is the process of exposure to low temperatures, which is crucial for the transition from vegetative to reproductive growth of plants. In this study, the global landscape vernalization-related mRNAs and long noncoding RNAs (lncRNAs) were identified in Beta vulgaris. A total of 22,159 differentially expressed mRNAs and 4418 differentially expressed lncRNAs were uncovered between the vernalized and nonvernalized samples. Various regulatory proteins, such as zinc finger CCCH domain-containing proteins, F-box proteins, flowering-time-related proteins FY and FPA, PHD finger protein EHD3 and B3 domain proteins were identified. Intriguingly, a novel vernalization-related lncRNA-mRNA target-gene co-expression regulatory network and the candidate vernalization genes, VRN1, VRN1-like, VAL1 and VAL2, encoding B3 domain-containing proteins were also unveiled. The results of this study pave the way for further illumination of the molecular mechanisms underlying the vernalization of B. vulgaris.

Genomic and transcriptomic-based analysis of agronomic traits in sugar beet (Beta vulgaris L.) pure line IMA1

Sugar beet (Beta vulgaris L.) is an important sugar-producing and energy crop worldwide. The sugar beet pure line IMA1 independently bred by Chinese scientists is a standard diploid parent material that is widely used in hybrid-breeding programs. In this study, a high-quality, chromosome-level genome assembly for IMA1was conducted, and 99.1% of genome sequences were assigned to nine chromosomes. A total of 35,003 protein-coding genes were annotated, with 91.56% functionally annotated by public databases. Compared with previously released sugar beet assemblies, the new genome was larger with at least 1.6 times larger N50 size, thereby substantially improving the completeness and continuity of the sugar beet genome. A Genome-Wide Association Studies analysis identified 10 disease-resistance genes associated with three important beet diseases and five genes associated with sugar yield per hectare, which could be key targets to improve sugar productivity. Nine highly expressed genes associated with pollen fertility of sugar beet were also identified. The results of this study provide valuable information to identify and dissect functional genes affecting sugar beet agronomic traits, which can increase sugar beet production and help screen for excellent sugar beet breeding materials. In addition, information is provided that can precisely incorporate biotechnology tools into breeding efforts.

Comparative proteomic analysis on chloroplast proteins provides new insights into the effects of low temperature in sugar beet

BACKGROUND

Low temperature, which is one of the main environmental factors that limits geographical distribution and sucrose yield, is a common abiotic stress during the growth and development of sugar beet. As a regulatory hub of plant response to abiotic stress, activity in the chloroplasts is related to many molecular and physiological processes, particularly in response to low temperature stress.

RESULTS

The contents of chlorophyll (Chl) and malondialdehyde (MDA), relative electrical conductivity (REL), and superoxide dismutase (SOD) activity were measured. The results showed that sugar beet could manage low temperature stress by regulating the levels of Chl, REL and MDA, and the activity of SOD. The physiological responses indicated that sugar beets respond positively to low temperature treatments and are not significantly damaged. Moreover, to determine the precise time to response low temperature in sugar beet, well-known abiotic stresses-responsive transcript factor family, namely DEHYDRATION RESPONSIVE ELEMENT BINDING PROTEIN (DREB), was selected as the marker gene. The results of phylogenetic analyses showed that BvDREBA1 and BvDREBA4 were in the same branch as the cold- and drought-responsive AtDREB gene. In addition, the expression of BvDREBs reached its maximum level at 24 h after low temperature by RNA-Seq and qRT-PCR analysis. Furthermore, the changes in chloroplast proteome after low temperature at 24 h were detected using a label-free technique. A total of 416 differentially expressed proteins were identified. GO enrichment analysis showed that 16 GO terms were significantly enriched, particularly chloroplast stroma, chloroplast envelope, and chloroplast thylakoid membrane. It is notable that the transport of photosynthetic proteins (BvLTD and BvTOC100), the formation of starch granules (BvPU1, BvISA3, and BvGWD3) and the scavenging of reactive oxygen species (BvCu/Zn-SOD, BvCAT, BvPrx, and BvTrx) were the pathways used by sugar beets to respond to low temperatures at an early stage.

CONCLUSIONS

These results provide a preliminarily analysis of how chloroplasts of sugar beet respond to low temperature stress at the translational level and provide a theoretical basis for breeding low temperature resistant varieties of sugar beet.

Transcriptional reprogramming of xylem cell wall biosynthesis in tension wood

Tension wood (TW) is a specialized xylem tissue developed under mechanical/tension stress in angiosperm trees. TW development involves transregulation of secondary cell wall genes, which leads to altered wood properties for stress adaptation. We induced TW in the stems of black cottonwood (Populus trichocarpa, Nisqually-1) and identified two significantly repressed transcription factor (TF) genes: class B3 heat-shock TF (HSFB3-1) and MYB092. Transcriptomic analysis and chromatin immunoprecipitation (ChIP) were used to identify direct TF-DNA interactions in P. trichocarpa xylem protoplasts overexpressing the TFs. This analysis established a transcriptional regulatory network in which PtrHSFB3-1 and PtrMYB092 directly activate 8 and 11 monolignol genes, respectively. The TF-DNA interactions were verified for their specificity and transactivator roles in 35 independent CRISPR-based biallelic mutants and overexpression transgenic lines of PtrHSFB3-1 and PtrMYB092 in P. trichocarpa. The gene-edited trees (mimicking the repressed PtrHSFB3-1 and PtrMYB092 under tension stress) have stem wood composition resembling that of TW during normal growth and under tension stress (i.e., low lignin and high cellulose), whereas the overexpressors showed an opposite effect (high lignin and low cellulose). Individual overexpression of the TFs impeded lignin reduction under tension stress and restored high levels of lignin biosynthesis in the TW. This study offers biological insights to further uncover how metabolism, growth, and stress adaptation are coordinately regulated in trees.

Transcriptome analysis of sugar beet (Beta vulgaris L.) in response to alkaline stress

RNA-seq was used to analyze the transcriptional changes in sugar beet (Beta vulgaris L.) triggered by alkaline solution to elucidate the molecular mechanism underlying alkaline tolerance in sugar beet. Several differentially expressed genes related to stress tolerance were identified. Our results provide a valuable resource for the breeding of new germplasms with high alkaline tolerance. Alkalinity is a highly stressful environmental factor that limits plant growth and production. Sugar beet own the ability to acclimate to various abiotic stresses, especially salt and alkaline stress. Although substantial previous studies on response of sugar beet to saline stress has been conducted, the expressions of alkali-responsive genes in sugar beet have not been comprehensively investigated. In this study, we conducted transcriptome analysis of leaves in sugar beet seedlings treated with alkaline solutions for 0 day (control, C), 3 days (short-term alkaline treatment, ST) and 7 days (long-term alkaline treatment, LT). The clean reads were obtained and assembled into 25,507 unigenes. Among them, 975 and 383 differentially expressed genes (DEGs) were identified in the comparison groups ST_vs_C and LT_vs_C, respectively. Gene ontology (GO) analysis revealed that oxidation-reduction process and lipid metabolic process were the most enriched GO term among the DEGs in ST_vs_C and LT_vs_C, respectively. According to Kyoto Encyclopedia of Genes and Genomes pathway, carbon fixation in photosynthetic organisms pathway were significantly enriched under alkaline stress. Besides, expression level of genes encoding D-3-phosphoglycerate dehydrogenase 1, glutamyl-tRNA reductase 1, fatty acid hydroperoxide lyase, ethylene-insensitive protein 2, metal tolerance protein 11 and magnesium-chelatase subunit ChlI, etc., were significantly altered under alkaline stress. Additionally, among the DEGs, 136 were non-annotated genes and 24 occurred with differential alternative splicing. Our results provide a valuable resource on alkali-responsive genes and should benefit the improvement of alkaline stress tolerance in sugar beet.

Mechanisms of Sugar Beet Response to Biotic and Abiotic Stresses

Sugar beet is used not only in the sugar production, but also in a wide range of industries including the production of bioethanol as a source of renewable energy, extraction of pectin and production of molasses. The red beetroot has attracted much attention as health-promoting and disease-preventing functional food. The negative effects of environmental stresses, including abiotic and biotic ones, significantly decrease the cash crop sugar beet productivity. In this paper, we outline the mechanisms of sugar beet response to biotic and abiotic stresses at the levels of physiological change, the genes' functions, transcription and translation. Regarding the physiological changes, most research has been carried out on salt and drought stress. The functions of genes from sugar beet in response to salt, cold and heavy metal stresses were mainly investigated by transgenic technologies. At the transcriptional level, the transcriptome analysis of sugar beet in response to salt, cold and biotic stresses were conducted by RNA-Seq or SSH methods. At the translational level, more than 800 differentially expressed proteins in response to salt, K⁺/Na⁺ ratio, iron deficiency and resupply and heavy metal (zinc) stress were identified by quantitative proteomics techniques. Understanding how sugar beet respond and tolerate biotic and abiotic stresses is important for boosting sugar beet productivity under these challenging conditions. In order to minimize the negative impact of these stresses, studying how the sugar beet has evolved stress coping mechanisms will provide new insights and lead to novel strategies for improving the breeding of stress-resistant sugar beet and other crops.

Comprehensive transcriptome analysis of faba bean in response to vernalization

This study unravels the transcriptional response of a highly productive faba bean cultivar under vernalization treatment. Faba bean (Vicia faba L.) is a member of the Leguminosae family and an important food crop worldwide providing valuable nutrients for humans. However, genome-wide studies and comprehensive sequencing resources of faba bean remain limited. Vernalization is crucial for enhanced yields in a number of winter-sown crops. However, the effects of vernalization on faba bean remain unknown. In this study, we generated a high-quality transcriptome assembly and functional annotation source for vernalized faba bean (Vicia faba L.) cv. Tongxian-2, a domesticated cultivar from southern China. A total of 369.9 million clean Illumina paired-end RNA-Seq reads were generated, and the transcriptome was assembled into 68,683 unigene sequences, with an average length of 1018 bp and an N50 of 1652 bp. Comprehensive functional annotation provided putative functional descriptions for more than 70% of the faba bean transcripts. We annotated a total of 1560 faba bean transcripts encoding transcription factors (TFs) belonging to 55 distinct TF families. The bHLH (168 transcripts), ERF (123 transcripts) and WRKY (105 transcripts) contained the largest number of TFs in response to vernalization. Genome-wide transcript changes comparing vernalized and unvernalized seedlings were investigated using bioinformatics approaches, which revealed a strong repression of photosynthesis and carbon metabolism, while genes participating in 'response to stress' were significantly induced. We also specifically identified vernalization-induced twenty-two 'pollen-pistil interaction' genes. A detailed functional annotation and expression profile analyses unveiled a number of protein kinases, which were specifically induced in vernalized seedlings. We also identified a total of 6852 simple sequence repeats (SSRs) in 6552 transcripts, representing a valuable genomic molecular marker resource for faba bean. In summary, this study provides new insights into the vernalization process in this economically valuable crop. The transcriptome data obtained provides us with a valuable candidate gene resource for future functional and molecular breeding studies. These data will contribute to the genome annotation for ensuing genome projects.

Towards sex identification of Asian Palmyra palm (Borassus flabellifer L.) by DNA fingerprinting, suppression subtractive hybridization and de novo transcriptome sequencing

BACKGROUND

Asian Palmyra palm, the source of palm-sugar, is dioecious with a long juvenile period requiring at least 12 years to reach its maturity. To date, there is no reliable molecular marker for identifying sexes before the first bloom, limiting crop designs and utilization. We aimed to identify sex-linked markers for this palm using PCR-based DNA fingerprinting, suppression subtractive hybridization (SSH) and transcriptome sequencing.

METHODS

DNA fingerprints were generated between males and females based on RAPD, AFLP, SCoT, modified SCoT, ILP, and SSR techniques. Large-scale cloning and screening of SSH libraries and de novo transcriptome sequencing of male and female cDNA from inflorescences were performed to identify sex-specific genes for developing sex-linked markers.

RESULTS

Through extensive screening and re-testing of the DNA fingerprints (up to 1,204 primer pairs) and transcripts from SSH (>10,000 clones) and transcriptome data, however, no sex-linked marker was identified. Although de novo transcriptome sequencing of male and female inflorescences provided ∼32 million reads and 187,083 assembled transcripts, PCR analysis of selected sex-highly represented transcripts did not yield any sex-linked marker. This result may suggest the complexity and small sex-determining region of the Asian Palmyra palm. To this end, we provide the first global transcripts of male and female inflorescences of Asian Palmyra palm. Interestingly, sequence annotation revealed a large proportion of transcripts related to sucrose metabolism, which corresponds to the sucrose-rich sap produced in the inflorescences, and these transcripts will be useful for further understanding of sucrose production in sugar crop plants. Provided lists of sex-specific and differential-expressed transcripts would be beneficial to the further study of sexual development and sex-linked markers in palms and related species.

Maize transcriptomic repertoires respond to gibberellin stimulation

Phytohormone gibberellin (GA) serves as hub modulator of diverse biological events. Understanding the transcriptomic features of GA-mediated processes has scientific significance. The transcriptomic landscapes of cereal crops upon GA stimulation remains largely unknown. Herein, to reveal the transcriptomic changes in cereal crop maize under GA treatment, we first selected normal height and GA-sensitive maize dwarf plants from advanced backcross population for GA treatment. RNA-seq analysis discovered multiple protein-coding transcripts that were differentially expressed in GA-treated samples compared to distilled water-treated ones. Some differentially expressed transcripts, namely GA-responsive transcripts in this study, encoded the components of GA pathway, including CPS, KS, and KO enzymes for GA biosynthesis, GA2ox enzymes for GA degradation, DELLA repressors and GID1 receptor for GA signaling. A total of 214 shared GA-responsive transcripts were identified both in GA₃-treated normal height and GA-sensitive dwarf samples. Shared GA-responsive transcripts were involved in GA signaling, auxin biosynthesis, ethylene response, the composition and structure of cell wall, chlorophyll biogenesis, and sugar homeostasis. In addition, the convergence and divergence in expression of shared GA-responsive transcripts were observed in GA₃-treated normal height and GA-sensitive dwarf plants. Interaction network modeling indicated that some shared GA-responsive transcripts tended to be co-regulated, which increases the complexity of GA-triggered regulation at transcriptomic layer. Results presented here will extend our knowledge of GA-mediated regulatory cascade, and enhance our ability to apply hormone GA knowledge in agricultural practice.

Beta vulgaris mitovirus 1 in diverse cultivars of beet and chard

Recent results indicate that mitoviruses, which replicate persistently in host mitochondria, are not restricted to fungi, but instead are found also in plants. Beta vulgaris mitovirus 1 (BevuMV1) is an example first discovered in sugar beet cultivars. For the current study, complete coding sequences of 42 BevuMV1 strains were newly determined, derived from not only sugar beet but also fodder beet, table beet, and Swiss chard cultivars of Beta vulgaris, as well as wild sea beet. BevuMV1 is thus a common phytobiome component of this valuable crop species. Most of the new BevuMV1 sequences originated from RNA extracted from B. vulgaris seed clusters, consistent with vertical transmission of this virus. Results suggest that BevuMV1 entered the B. vulgaris lineage prior to human cultivation and also provides a marker for tracing the maternal ancestry of B. vulgaris cultivars. Especially notable is the monophyletic relationship and limited sequence divergence among BevuMV1 strains from cultivars that are thought or shown to share the "Owen" trait for cytoplasmic male sterility, which is transmitted by maternal mitochondria and has been broadly established in commercial breeding lines of B. vulgaris since the mid-20^th century.

RNA-Seq analysis of differential gene expression in Betula luminifera xylem during the early stages of tension wood formation

Background

Betula luminifera H. Winkler, which is widely distributed in southern China, is an economically important broadleaf tree species. However, little genomic information of B. luminifera is available, and little is known about the molecular mechanisms of wood formation in this species. Meanwhile, few efforts have focused on investigating the early transcriptional changes during tension wood formation in woody plants.

Results

A reference transcriptome dataset was first generated containing 45,700 Unigenes, and 35,135 (76.9%) Unigenes were annotated by a BLAST similarity search against four public databases. Then, based on an anatomical investigation, the global gene expression changes during the early stages of tension wood formation were analyzed. Gene expression profiling showed that a total of 13,273 Unigenes were differentially regulated during the early stages of tension wood formation. Most genes involved in cellulose and lignin biosynthesis were highlighted to reveal their biological importance in tension wood formation. In addition, the transcription levels of many genes involved in the auxin response pathway were significantly changed during the early stages of tension wood formation. Furthermore, 18 TFs co-expressed with key enzymes of cellulose synthesis were identified.

Conclusions

Our results revealed the transcriptional changes associated with TW formation and identified potential key genes in the regulation of this process. These results will help to dissect the molecular mechanism of wood formation and provide key candidate genes for marker-assisted selection in B. luminifera.

Evidence for contemporary plant mitoviruses

Mitoviruses have small RNA(+) genomes, replicate in mitochondria, and have been shown to infect only fungi to date. For this report, sequences that appear to represent nearly complete plant mitovirus genomes were recovered from publicly available transcriptome data. Twenty of the refined sequences, 2684-2898 nt long and derived from 10 different species of land plants, appear to encompass the complete coding regions of contemporary plant mitoviruses, which furthermore constitute a monophyletic cluster within genus Mitovirus. Complete coding sequences of several of these viruses were recovered from multiple transcriptome (but not genome) studies of the same plant species and also from multiple plant tissues. Crop plants among implicated hosts include beet and hemp. Other new results suggest that such genuine plant mitoviruses were immediate ancestors to endogenized mitovirus elements now widespread in land plant genomes. Whether these mitoviruses are wholly cryptic with regard to plant health remains to be investigated.

Investigating the molecular basis for heterophylly in the aquatic plant Potamogeton octandrus (Potamogetonaceae) with comparative transcriptomics

Many plant species exhibit different leaf morphologies within a single plant, or heterophylly. The molecular mechanisms regulating this phenomenon, however, have remained elusive. In this study, the transcriptomes of submerged and floating leaves of an aquatic heterophyllous plant, Potamogeton octandrus Poir, at different stages of development, were sequenced using high-throughput sequencing (RNA-Seq), in order to aid gene discovery and functional studies of genes involved in heterophylly. A total of 81,103 unigenes were identified in submerged and floating leaves and 6,822 differentially expressed genes (DEGs) were identified by comparing samples at differing time points of development. KEGG pathway enrichment analysis categorized these unigenes into 128 pathways. A total of 24,025 differentially expressed genes were involved in carbon metabolic pathways, biosynthesis of amino acids, ribosomal processes, and plant-pathogen interactions. In particular, KEGG pathway enrichment analysis categorized a total of 70 DEGs into plant hormone signal transduction pathways. The high-throughput transcriptomic results presented here highlight the potential for understanding the molecular mechanisms underlying heterophylly, which is still poorly understood. Further, these data provide a framework to better understand heterophyllous leaf development in P. octandrus via targeted studies utilizing gene cloning and functional analyses.

Difference of proteomics vernalization-induced in bolting and flowering transitions of Beta vulgaris

Sugar beet (Beta vulgaris) is a biennial crop that accounts for 30% sugar production of the world. Vernalization is an essential factor for sugar beet reproductative growth under long days. Although genes association with bolting and flowering were well explored, the difference of proteomics in the two growth stages were still poorly understood. To address the molecular mechanism at the level of proteins, an isobaric tags for relative and absolute quantification (iTRAQ)-based quantitative proteomics approach was employed to the three different growth stages (germination, bolting, flowering) of vernalized samples and the corresponding stage germination (17W weeks), 19W and 20W of nonvernalized samples. A total of 1110 peptides, 842 unique peptides and 570 proteins were identified. Most of them were assigned to phenylpropanoid biosynthesis, hormone metabolism and protein processing pathway. IAA and Gibberellins (GA3) promoted growth and development in a threshold manner at growth stage germination after vernalization. A novel discovery was that IAA biosynthetic pathway of sugar beet was the Trp-dependent. In addition, two predominant pathways of protein processing association with vernalization were also identified in sugar beet at growth stage flowering. This study provided an in-depth understanding of the molecular mechanism of vernalization at the level of proteomics.

Transcriptome profiling reveals the immune response of goose T cells under selenium stimuli

The goose is an economically important poultry species and a principal natural host of avian viruses. This study aimed to determine the effects of selenium on the immune response of geese. Under selenium stimulation, gene expression profiling was investigated using transcriptome sequencing. The selenoproteins were promoted by selenium stimulation, while the heat shock proteins, interleukin and interferons were mainly down-regulated. After comparison, 2228 differentially expressed genes were primarily involved in immune and environmental response, and infectious disease and genetic information processing related pathways were identified. Specifically, the enzymes of the lysosomes which acted as a safeguard in preventing pathogens were mostly up-regulated and six randomly selected differentially expressed genes were validated by quantitative polymerase chain reaction. In addition, the most proportional increased transcription factor family basic helix-loop-helix (bHLH) located in the 5' flank of selenoprotein P-like protein for selenium metabolism was identified by response to the selenium stimulation in this study. These analyses show that selenium can promote immune function by activating selenoproteins, transcript factors and lysosome pathway related genes, while weakening cytokine content genes in geese.

Vernalisation mediated LncRNA-like gene expression in Beta vulgaris

Sugar beet (Beta vulgaris L.) cannot form reproductive shoots during the first year of their life cycle. Flowering only occurs if plants are vernalised and are subsequently exposed to long days. However, the vernalisation mechanism remains poorly understood in sugar beet. Three putative lncRNAs associated with vernalisation (AGL15X1, AGL15X2 and CAULIFLOWER A) were investigated and the hypothesis that their expression occurred in response to vernalisation was experimentally tested. The regulation mechanisms of BvRAV1-like, lncRNA-like genes, BvFT1 and BvFT2 were also examined. The BvRAV1-like gene associated with vernalisation in sugar beet was validated for the first time. Our data confirmed the hypothesis that AGLX2 was the first candidate lncRNA of sugar beet and the BvRAV1-like gene was expressed in response to vernalisation. BvRAV1-like and AGLX2 genes might be coordinated with BvFT2 to promote reproductive growth by repressing BvFT1 during cold exposure followed by long day conditions. A new complementary flowering model of sugar beet was proposed. Our findings opened up new possibility for future studies and further illuminated the molecular mechanism of vernalisation in sugar beet.

Characterising genes associated with flowering time in carrot (Daucus carota L.) using transcriptome analysis

Carrot is generally regarded as a biennial plant with an obligatory vernalization requirement. Early spring cultivation makes plants vulnerable to premature bolting, which results in a loss of commercial value. However, our knowledge of flowering time genes and flowering mechanisms in carrot remain limited. Bolting behavior of D. carota ssp. carota 'Songzi', a wild species sensitive to flower induction by vernalization and photoperiod, and orange cultivar 'Amsterdam forcing', and their offspring were investigated in different growing conditions. We performed RNA-seq to identify the flowering time genes, and digital gene expression (DGE) analysis to examine their expression levels. The circadian patterns of related genes were identified by qPCR. The results showed bolting behavior of carrot was influenced by low temperature, illumination intensity and photoperiod. A total of 45 flowering time-related unigenes were identified, which were classified into five categories including photoperiod, vernalization, autonomous and gibberellin pathway, and floral integrators. Homologs of LATE ELONGATED HYPOCOTYL (LHY) and CONSTANS-LIKE 2 (COL2) were more highly expressed under short day condition than under long day condition. Homologs of COL2, CONSTANS-LIKE 5 (COL5), SUPPRESSION OF OVEREXPRESSION OF CONSTANS 1 (SOC1), FLOWERING LOCUS C (FLC) and GIBBERELLIC ACID INSENSITIVE (GAI) were differentially expressed between 'Songzi' and 'Amsterdam forcing'. The homolog of COL2 (Dct43207) was repressed by light, but that of COL5 (Dct20940) was induced. A preliminary model of genetic network controlling flowering time was constructed by associating the results of DGE analysis with correlation coefficients between genes. This study provides useful information for further investigating the genetic mechanism of flowering in carrot.

Winter Memory throughout the Plant Kingdom: Different Paths to Flowering

Molecular mechanisms contribute to the memory of winter in different plant groups.

Association genetics and transcriptome analysis reveal a gibberellin-responsive pathway involved in regulating photosynthesis

Gibberellins (GAs) regulate a wide range of important processes in plant growth and development, including photosynthesis. However, the mechanism by which GAs regulate photosynthesis remains to be understood. Here, we used multi-gene association to investigate the effect of genes in the GA-responsive pathway, as constructed by RNA sequencing, on photosynthesis, growth, and wood property traits, in a population of 435 Populus tomentosa By analyzing changes in the transcriptome following GA treatment, we identified many key photosynthetic genes, in agreement with the observed increase in measurements of photosynthesis. Regulatory motif enrichment analysis revealed that 37 differentially expressed genes related to photosynthesis shared two essential GA-related cis-regulatory elements, the GA response element and the pyrimidine box. Thus, we constructed a GA-responsive pathway consisting of 47 genes involved in regulating photosynthesis, including GID1, RGA, GID2, MYBGa, and 37 photosynthetic differentially expressed genes. Single nucleotide polymorphism (SNP)-based association analysis showed that 142 SNPs, representing 40 candidate genes in this pathway, were significantly associated with photosynthesis, growth, and wood property traits. Epistasis analysis uncovered interactions between 310 SNP-SNP pairs from 37 genes in this pathway, revealing possible genetic interactions. Moreover, a structural gene-gene matrix based on a time-course of transcript abundances provided a better understanding of the multi-gene pathway affecting photosynthesis. The results imply a functional role for these genes in mediating photosynthesis, growth, and wood properties, demonstrating the potential of combining transcriptome-based regulatory pathway construction and genetic association approaches to detect the complex genetic networks underlying quantitative traits.

Transcriptomic profiling of the salt stress response in excised leaves of the halophyte Beta vulgaris ssp. maritima

Beta vulgaris ssp. maritima is a halophytic relative of cultivated beets. In the present work a transcriptome response to acute salt stress imposed to excised leaves of sea beet was investigated. Salt treatments consisted of adding NaCl directly to the transpiration stream by immersing the petioles of excised leaves into the salt solutions. Sequencing libraries were generated from leaves subjected to either moderate or strong salt stress. Control libraries were constructed from untreated leaves. Sequencing was performed using the Illumina MiSeq platform. We obtained 32970 unigenes by assembling the pooled reads from all the libraries with Trinity software. Screening the nr database returned 18,362 sequences with functional annotation. Using the reference transcriptome we identified 1,246 genes that were differentially expressed after 48 h of NaCl stress. Genes related to several cellular functions such as membrane transport, osmoprotection, molecular chaperoning, redox metabolism or protein synthesis were differentially expressed in response to salt stress. The response of sea beet leaves to salt treatments was marked out by transcriptomic up-regulation of genes related to photosynthetic carbon fixation, ribosome biogenesis, cell wall-building and cell wall expansion. Furthermore, several novel and undescribed transcripts were responsive to salinity in leaves of sea beet.

Epigenomics and bolting tolerance in sugar beet genotypes

In sugar beet (Beta vulgaris altissima), bolting tolerance is an essential agronomic trait reflecting the bolting response of genotypes after vernalization. Genes involved in induction of sugar beet bolting have now been identified, and evidence suggests that epigenetic factors are involved in their control. Indeed, the time course and amplitude of DNA methylation variations in the shoot apical meristem have been shown to be critical in inducing sugar beet bolting, and a few functional targets of DNA methylation during vernalization have been identified. However, molecular mechanisms controlling bolting tolerance levels among genotypes are still poorly understood. Here, gene expression and DNA methylation profiles were compared in shoot apical meristems of three bolting-resistant and three bolting-sensitive genotypes after vernalization. Using Cot fractionation followed by 454 sequencing of the isolated low-copy DNA, 6231 contigs were obtained that were used along with public sugar beet DNA sequences to design custom Agilent microarrays for expression (56k) and methylation (244k) analyses. A total of 169 differentially expressed genes and 111 differentially methylated regions were identified between resistant and sensitive vernalized genotypes. Fourteen sequences were both differentially expressed and differentially methylated, with a negative correlation between their methylation and expression levels. Genes involved in cold perception, phytohormone signalling, and flowering induction were over-represented and collectively represent an integrative gene network from environmental perception to bolting induction. Altogether, the data suggest that the genotype-dependent control of DNA methylation and expression of an integrative gene network participate in bolting tolerance in sugar beet, opening up perspectives for crop improvement.

The Pokeweed Leaf mRNA Transcriptome and Its Regulation by Jasmonic Acid

The American pokeweed plant, Phytolacca americana, is recognized for synthesizing pokeweed antiviral protein (PAP), a ribosome inactivating protein (RIP) that inhibits the replication of several plant and animal viruses. The plant is also a heavy metal accumulator with applications in soil remediation. However, little is known about pokeweed stress responses, as large-scale sequencing projects have not been performed for this species. Here, we sequenced the mRNA transcriptome of pokeweed in the presence and absence of jasmonic acid (JA), a hormone mediating plant defense. Trinity-based de novo assembly of mRNA from leaf tissue and BLASTx homology searches against public sequence databases resulted in the annotation of 59 096 transcripts. Differential expression analysis identified JA-responsive genes that may be involved in defense against pathogen infection and herbivory. We confirmed the existence of several PAP isoforms and cloned a potentially novel isoform of PAP. Expression analysis indicated that PAP isoforms are differentially responsive to JA, perhaps indicating specialized roles within the plant. Finally, we identified 52 305 natural antisense transcript pairs, four of which comprised PAP isoforms, suggesting a novel form of RIP gene regulation. This transcriptome-wide study of a Phytolaccaceae family member provides a source of new genes that may be involved in stress tolerance in this plant. The sequences generated in our study have been deposited in the SRA database under project # SRP069141.

Comprehensive Analysis of the Triterpenoid Saponins Biosynthetic Pathway in Anemone flaccida by Transcriptome and Proteome Profiling

BACKGROUND

Anemone flaccida Fr. Shmidt (Ranunculaceae), commonly known as 'Di Wu' in China, is a perennial herb with limited distribution. The rhizome of A. flaccida has long been used to treat arthritis as a tradition in China. Studies disclosed that the plant contains a rich source of triterpenoid saponins. However, little is known about triterpenoid saponins biosynthesis in A. flaccida.

RESULTS

In this study, we conducted the tandem transcriptome and proteome profiling of a non-model medicinal plant, A. flaccida. Using Illumina HiSeq 2000 sequencing and iTRAQ technique, a total of 46,962 high-quality unigenes were obtained with an average sequence length of 1,310 bp, along with 1473 unique proteins from A. flaccida. Among the A. flaccida transcripts, 36,617 (77.97%) showed significant similarity (E-value < 1e (-5)) to the known proteins in the public database. Of the total 46,962 unigenes, 36,617 open reading frame (ORFs) were predicted. By the fragments per kilobases per million reads (FPKM) statistics, 14,004 isoforms/unigenes were found to be upregulated, and 14,090 isoforms/unigenes were down-regulated in the rhizomes as compared to those in the leaves. Based on the bioinformatics analysis, all possible enzymes involved in the triterpenoid saponins biosynthetic pathway of A. flaccida were identified, including cytosolic mevalonate pathway (MVA) and the plastidial methylerythritol pathway (MEP). Additionally, a total of 126 putative cytochrome P450 (CYP450) and 32 putative UDP glycosyltransferases were selected as the candidates of triterpenoid saponins modifiers. Among them, four of them were annotated as the gene of CYP716A subfamily, the key enzyme in the oleanane-type triterpenoid saponins biosynthetic pathway. Furthermore, based on RNA-Seq and proteome analysis, as well as quantitative RT-PCR verification, the expression level of gene and protein committed to triterpenoids biosynthesis in the leaf versus the rhizome was compared.

CONCLUSION

A combination of the de novo transcriptome and proteome profiling based on the Illumina HiSeq 2000 sequencing platform and iTRAQ technique was shown to be a powerful method for the discovery of candidate genes, which encoded enzymes that were responsible for the biosynthesis of novel secondary metabolites in a non-model plant. The transcriptome data of our study provides a very important resource for the understanding of the triterpenoid saponins biosynthesis of A. flaccida.

AtRAV1 and AtRAV2 overexpression in cotton increases fiber length differentially under drought stress and delays flowering

There is a longstanding problem of an inverse relationship between cotton fiber qualities versus high yields. To better understand drought stress signaling and adaptation in cotton (Gossypium hirsutum) fiber development, we expressed the Arabidopsis transcription factors RELATED_TO_ABA-INSENSITIVE3/VIVIPAROUS1/(RAV1) and AtRAV2, which encode APETALA2-Basic3 domain proteins shown to repress transcription of FLOWERING_LOCUS_T (FT) and to promote stomatal opening cell-autonomously. In three years of field trials, we show that AtRAV1 and AtRAV2-overexpressing cotton had ∼5% significantly longer fibers with only marginal decreases in yields under well-watered or drought stress conditions that resulted in 40-60% yield penalties and 3-7% fiber length penalties in control plants. The longer transgenic fibers from drought-stressed transgenics could be spun into yarn which was measurably stronger and more uniform than that from well-watered control fibers. The transgenic AtRAV1 and AtRAV2 lines flowered later and retained bolls at higher nodes, which correlated with repression of endogenous GhFT-Like (FTL) transcript accumulation. Elevated expression early in development of ovules was observed for GhRAV2L, GhMYB25-Like (MYB25L) involved in fiber initiation, and GhMYB2 and GhMYB25 involved in fiber elongation. Altered expression of RAVs controlling critical nodes in developmental and environmental signaling hierarchies has the potential for phenotypic modification of crops.

Root hydraulic conductivity and adjustments in stomatal conductance: hydraulic strategy in response to salt stress in a halotolerant species

Recent advances at the molecular level are introducing a new scenario that needs to be integrated into the analysis of plant hydraulic properties. Although it is not yet clear to what extent this scenario alters the current proposal for the hydraulic circuit models, it introduces new insights when studying plants that are able to easily overcome water restrictions. In this context, our aim was to explore water adjustments in a halotolerant model (Beta vulgaris) by studying the coordination between the root in terms of root hydraulic conductivity (Lpr) and the shoot as reflected in the stomatal conductance (gs). The root water pathways were also analysed in terms of root suberization (apoplastic barrier) and aquaporin transcript levels (cell-to-cell pathway). Beta vulgaris showed the ability to rapidly lose (4 h) and gain (24 h) turgor when submitted to salt stress (200 mM). The reduction profile observed in Lpr and gs was consistent with a coupled process. The tuning of the root water flow involved small variations in the studied aquaporin's transcripts before anatomical modifications occurred. Exploring Lpr enhancement after halting the stress contributed to show not only a different profile in restoring Lpr but also the capacity to uncouple Lpr from gs. Beta vulgaris root plays a key role and can anticipate water loss before the aerial water status is affected.

Exploiting single-molecule transcript sequencing for eukaryotic gene prediction

We develop a method to predict and validate gene models using PacBio single-molecule, real-time (SMRT) cDNA reads. Ninety-eight percent of full-insert SMRT reads span complete open reading frames. Gene model validation using SMRT reads is developed as automated process. Optimized training and prediction settings and mRNA-seq noise reduction of assisting Illumina reads results in increased gene prediction sensitivity and precision. Additionally, we present an improved gene set for sugar beet (Beta vulgaris) and the first genome-wide gene set for spinach (Spinacia oleracea). The workflow and guidelines are a valuable resource to obtain comprehensive gene sets for newly sequenced genomes of non-model eukaryotes.

Whole transcriptome profiling of the vernalization process in Lilium longiflorum (cultivar White Heaven) bulbs

Background

Vernalization is an obligatory requirement of extended exposure to low temperatures to induce flowering in certain plants. It is the most important factor affecting flowering time and quality in Easter lily (Lilium longiflorum). Exposing the bulbs to 4 °C gradually decreases flowering time up to 50 % compared to non-vernalized plants. We aim to understand the molecular regulation of vernalization in Easter lily, for which we characterized the global expression in lily bulb meristems after 0, 2, 5, 7 and 9 weeks of incubation at 4 °C.

Results

We assembled de-novo a transcriptome which, after filtering, yielded 121,572 transcripts and 42,430 genes which hold 15,414 annotated genes, with up to 3,657 GO terms. This extensive annotation was mapped to the more general GO slim plant with a total of 94 terms. The response to cold exposure was summarized in 6 expression clusters, providing useful patterns for dissecting the dynamics of vernalization in lily. The functional annotation (GO and GO slim plant) was used to group transcripts in gene sets. Analysis of these gene sets and profiles revealed that most of the enriched functions among genes up-regulated by cold exposure were related to epigenetic processes and chromatin remodeling. Candidate vernalization genes in lily were selected based on their sequence similarity to known regulators of flowering in other species.

Conclusions

We present a detailed analysis of gene expression dynamics during vernalization in Lilium, covering several time points and accounting for biological variation by the use of replicates. The resulting collection of transcripts and novel isoforms provides a useful resource for studying the changes occurring during vernalization at a fine level. The selected potential candidate genes can shed light on the regulation of this process.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1675-1) contains supplementary material, which is available to authorized users.

Early transcriptional changes in Beta vulgaris in response to low temperature

Major metabolic pathways and genes affected by low-temperature treatment were identified and a thorough picture of the early transcriptional changes in sugar beet plantlets upon cold stress was given. Sugar beet (Beta vulgaris L.) is an important source of sugar and bioethanol production in temperate areas worldwide. In these areas, plantlet survival and sucrose yield of mature plants can be seriously limited by low temperatures, especially when plantlets are exposed to freezing temperatures (below 0 °C) at the early developmental stages. This frequently occurs when the crop is sown in early spring or even in autumn (autumn sowing) to escape drought at maturity and pathogen outbreaks. The knowledge of molecular responses induced in plantlets early upon exposure to low temperature is necessary to understand mechanisms that allow the plant to survive and to identify reactions that can influence other late-appearing traits. In this work, a wide study of sugar beet transcriptome modulation after a short exposure to a cold stress, mimicking what is experienced in vivo by young plantlets when temperature drops in the early spring nights, was carried out by high-throughput sequencing of leaves and root RNAs (RNA-Seq). A significant picture of the earliest events of temperature sensing was achieved for the first time for sugar beet: the retrieval of a great amount of transcription factors and the intensity of modulation of a large number of genes involved in several metabolic pathways suggest a fast and deep rearrangement of sugar beet plantlets metabolism as early response to cold stress, with both similarities and specificities between the two organs.

Reliable in silico identification of sequence polymorphisms and their application for extending the genetic map of sugar beet (Beta vulgaris)

Molecular markers are a highly valuable tool for creating genetic maps. Like in many other crops, sugar beet (Beta vulgaris L.) breeding is increasingly supported by the application of such genetic markers. Single nucleotide polymorphism (SNP) based markers have a high potential for automated analysis and high-throughput genotyping. We developed a bioinformatics workflow that uses Sanger and 2nd-generation sequence data for detection, evaluation and verification of new transcript-associated SNPs from sugar beet. RNAseq data from one parent of an established mapping population were produced by 454-FLX sequencing and compared to Sanger ESTs derived from the other parent. The workflow established for SNP detection considers the quality values of both types of reads, provides polymorphic alignments as well as selection criteria for reliable SNP detection and allows painless generation of new genetic markers within genes. We obtained a total of 14,323 genic SNPs and InDels. According to empirically optimised settings for the quality parameters, we classified these SNPs into four usability categories. Validation of a subset of the in silico detected SNPs by genotyping the mapping population indicated a high success rate of the SNP detection. Finally, a total of 307 new markers were integrated with existing data into a new genetic map of sugar beet which offers improved resolution and the integration of terminal markers.

Deep sequencing-based comparative transcriptional profiles of Cymbidium hybridum roots in response to mycorrhizal and non-mycorrhizal beneficial fungi

Background

The Orchidaceae is one of the largest families in the plant kingdom and orchid mycorrhizae (OM) are indispensable in the life cycle of all orchids under natural conditions. In spite of this, little is known concerning the mechanisms underlying orchid- mycorrhizal fungi interactions. Our previous work demonstrated that the non-mycorrhizal fungus Umbelopsis nana ZH3A-3 could improve the symbiotic effects of orchid mycorrhizal fungus Epulorhiza repens ML01 by co-cultivation with Cymbidium hybridum plantlets. Thus, we investigated the C. hybridum transcript profile associated with different beneficial fungi.

Results

More than 54,993,972 clean reads were obtained from un-normalized cDNA library prepared from fungal- and mock- treated Cymbidium roots at four time points using RNA-seq technology. These reads were assembled into 16,798 unique transcripts, with a mean length of 1127 bp. A total of 10,971 (65.31%) sequences were annotated based on BLASTX results and over ninety percent of which were assigned to plant origin. The digital gene expression profiles in Cymbidium root at 15 days post inoculation revealed that 1674, 845 and 1743 genes were sigificantly regulated in response to ML01, ZH3A-3 and ML01+ ZH3A-3 treatments, respectively. Twenty-six genes in different regulation patterns were validated using quantitative RT-PCR. Our analysis showed that general defense responses were co- induced by three treatments, including cell wall modification, reactive oxygen species detoxification, secondary biosynthesis and hormone balance. Genes involved in phosphate transport and root morphogenesis were also detected to be up-regulated collectively. Among the OM specifically induced transcripts, genes related to signaling, protein metabolism and processing, defense, transport and auxin response were identifed. Aside from these orchid transcripts, some putative fungal genes were also identified in symbiotic roots related to plant cell wall degradation, remodeling the fungal cell wall and nutrient transport.

Conclusion

The orchid root transcriptome will facilitate our understanding of orchid - associated biological mechanism. The comparative expression profiling revealed that the transcriptional reprogramming by OM symbiosis generally overlapped that of arbuscular mycorrhizas and ectomycorrhizas. The molecular basis of OM formation and function will improve our knowledge of plant- mycorrhzial fungi interactions, and their effects on plant and fungal growth, development and differentiation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-747) contains supplementary material, which is available to authorized users.

A scalable and accurate targeted gene assembly tool (SAT-Assembler) for next-generation sequencing data

Gene assembly, which recovers gene segments from short reads, is an important step in functional analysis of next-generation sequencing data. Lacking quality reference genomes, de novo assembly is commonly used for RNA-Seq data of non-model organisms and metagenomic data. However, heterogeneous sequence coverage caused by heterogeneous expression or species abundance, similarity between isoforms or homologous genes, and large data size all pose challenges to de novo assembly. As a result, existing assembly tools tend to output fragmented contigs or chimeric contigs, or have high memory footprint. In this work, we introduce a targeted gene assembly program SAT-Assembler, which aims to recover gene families of particular interest to biologists. It addresses the above challenges by conducting family-specific homology search, homology-guided overlap graph construction, and careful graph traversal. It can be applied to both RNA-Seq and metagenomic data. Our experimental results on an Arabidopsis RNA-Seq data set and two metagenomic data sets show that SAT-Assembler has smaller memory usage, comparable or better gene coverage, and lower chimera rate for assembling a set of genes from one or multiple pathways compared with other assembly tools. Moreover, the family-specific design and rapid homology search allow SAT-Assembler to be naturally compatible with parallel computing platforms. The source code of SAT-Assembler is available at https://sourceforge.net/projects/sat-assembler/. The data sets and experimental settings can be found in supplementary material.

Next-generation sequencing (NGS) provides an efficient and affordable way to sequence the genomes or transcriptomes of a large amount of organisms. With fast accumulation of the sequencing data from various NGS projects, the bottleneck is to efficiently mine useful knowledge from the data. As NGS platforms usually generate short and fragmented sequences (reads), one key step to annotate NGS data is to assemble short reads into longer contigs, which are then used to recover functional elements such as protein-coding genes. Short read assembly remains one of the most difficult computational problems in genomics. In particular, the performance of existing assembly tools is not satisfactory on complicated NGS data sets. They cannot reliably separate genes of high similarity, recover under-represented genes, and incur high computational time and memory usage. Hence, we propose a targeted gene assembly tool, SAT-Assembler, to assemble genes of interest directly from NGS data with low memory usage and high accuracy. Our experimental results on a transcriptomic data set and two microbial community data sets showed that SAT-Assembler used less memory and recovered more target genes with better accuracy than existing tools.

Genetic and molecular bases of yield-associated traits: a translational biology approach between rice and wheat

Transferring the knowledge bases between related species may assist in enlarging the yield potential of crop plants. Being cereals, rice and wheat share a high level of gene conservation; however, they differ at metabolic levels as a part of the environmental adaptation resulting in different yield capacities. This review focuses on the current understanding of genetic and molecular regulation of yield-associated traits in both crop species, highlights the similarities and differences and presents the putative knowledge gaps. We focus on the traits associated with phenology, photosynthesis, and assimilate partitioning and lodging resistance; the most important drivers of yield potential. Currently, there are large knowledge gaps in the genetic and molecular control of such major biological processes that can be filled in a translational biology approach in transferring genomics and genetics informations between rice and wheat.

Related to ABA-Insensitive3(ABI3)/Viviparous1 and AtABI5 transcription factor coexpression in cotton enhances drought stress adaptation

Drought tolerance is an important trait being pursued by the agbiotech industry. Abscisic acid (ABA) is a stress hormone that mediates a multitude of processes in growth and development, water use efficiency (WUE) and gene expression during seed development and in response to environmental stresses. Arabidopsis B3-domain transcription factor Related to ABA-Insensitive3 (ABI3)/Viviparous1 (namely AtRAV2) and basic leucine zipper (bZIPs) AtABI5 or AtABF3 transactivated ABA-inducible promoter:GUS reporter expression in a maize mesophyll protoplast transient assay and showed synergies in reporter transactivation when coexpressed. Transgenic cotton (Gossypium hirsutum) expressing AtRAV1/2 and/or AtABI5 showed resistance to imposed drought stress under field and greenhouse conditions and exhibited improved photosynthesis and WUEs associated with absorption through larger root system and greater leaf area. We observed synergy for root biomass accumulation in the greenhouse, intrinsic WUE in the field and drought tolerance in stacked AtRAV and AtABI5 double-transgenic cotton. We assessed AtABI5 and AtRAV1/2 involvement in drought stress adaptations through reactive oxygen species scavenging and osmotic adjustment by marker gene expression in cotton. Deficit irrigation-grown AtRAV1/2 and AtABI5 transgenics had 'less-stressed' molecular and physiological phenotypes under drought, likely due to improved photoassimilation and root and shoot sink strengths and enhanced expression of endogenous GhRAV and genes for antioxidant and osmolyte biosynthesis. Overexpression of bZIP and RAV TFs could impact sustainable cotton agriculture and potentially other crops under limited irrigation conditions.

RNA sequencing read depth requirement for optimal transcriptome coverage in Hevea brasiliensis

BACKGROUND

One of the concerns of assembling de novo transcriptomes is determining the amount of read sequences required to ensure a comprehensive coverage of genes expressed in a particular sample. In this report, we describe the use of Illumina paired-end RNA-Seq (PE RNA-Seq) reads from Hevea brasiliensis (rubber tree) bark to devise a transcript mapping approach for the estimation of the read amount needed for deep transcriptome coverage.

FINDINGS

We optimized the assembly of a Hevea bark transcriptome based on 16 Gb Illumina PE RNA-Seq reads using the Oases assembler across a range of k-mer sizes. We then assessed assembly quality based on transcript N50 length and transcript mapping statistics in relation to (a) known Hevea cDNAs with complete open reading frames, (b) a set of core eukaryotic genes and (c) Hevea genome scaffolds. This was followed by a systematic transcript mapping process where sub-assemblies from a series of incremental amounts of bark transcripts were aligned to transcripts from the entire bark transcriptome assembly. The exercise served to relate read amounts to the degree of transcript mapping level, the latter being an indicator of the coverage of gene transcripts expressed in the sample. As read amounts or datasize increased toward 16 Gb, the number of transcripts mapped to the entire bark assembly approached saturation. A colour matrix was subsequently generated to illustrate sequencing depth requirement in relation to the degree of coverage of total sample transcripts.

CONCLUSIONS

We devised a procedure, the "transcript mapping saturation test", to estimate the amount of RNA-Seq reads needed for deep coverage of transcriptomes. For Hevea de novo assembly, we propose generating between 5-8 Gb reads, whereby around 90% transcript coverage could be achieved with optimized k-mers and transcript N50 length. The principle behind this methodology may also be applied to other non-model plants, or with reads from other second generation sequencing platforms.

Transcriptome-wide analysis of vernalization reveals conserved and species-specific mechanisms in Brachypodium

Several temperate cereals need vernalization to promote flowering. Little, however, is known about the vernalization-memory-related genes, and almost no comparative analysis has been performed. Here, RNA-Seq was used for transcriptome analysis in non-vernalized, vernalized and post-vernalized Brachypodium distachyon (L.) Beauv. seedlings. In total, the expression of 1,665 genes showed significant changes (fold change ≥4) in response to vernalization. Among them, 674 putative vernalization-memory-related genes with a constant response to vernalization were significantly enriched in transcriptional regulation and monooxygenase-mediated biological processes. Comparative analysis of vernalization-memory-related genes with barley demonstrated that the oxidative-stress response was the most conserved pathway between these two plant species. Moreover, Brachypodium preferred to regulate transcription and protein phosphorylation processes, while vernalization-memory-related genes, whose products are cytoplasmic membrane-bound-vesicle-located proteins, were preferred to be regulated in barley. Correlation analysis of the vernalization-related genes with barley revealed that the vernalization mechanism was conserved between these two plant species. In summary, vernalization, including its memory mechanism, is conserved between Brachypodium and barley, although several species-specific features also exist. The data reported here will provide primary resources for subsequent functional research in vernalization.

Tracing the transcriptomic changes in synthetic Trigenomic allohexaploids of Brassica using an RNA-Seq approach

Polyploidization has played an important role in plant evolution and speciation, and newly formed allopolyploids have experienced rapid transcriptomic changes. Here, we compared the transcriptomic differences between a synthetic Brassica allohexaploid and its parents using a high-throughput RNA-Seq method. A total of 35,644,409 sequence reads were generated, and 32,642 genes were aligned from the data. Totals of 29,260, 29,060, and 29,697 genes were identified in Brassicarapa, Brassicacarinata, and Brassica allohexaploid, respectively. We compared 7,397 differentially expressed genes (DEGs) between Brassica hexaploid and its parents, as well as 2,545 nonadditive genes of Brassica hexaploid. We hypothesized that the higher ploidy level as well as secondary polyploidy might have influenced these changes. The majority of the 3,184 DEGs between Brassica hexaploid and its paternal parent, B. rapa, were involved in the biosynthesis of secondary metabolites, plant–pathogen interactions, photosynthesis, and circadian rhythm. Among the 2,233 DEGs between Brassica hexaploid and its maternal parent, B. carinata, several played roles in plant–pathogen interactions, plant hormone signal transduction, ribosomes, limonene and pinene degradation, photosynthesis, and biosynthesis of secondary metabolites. There were more significant differences in gene expression between the allohexaploid and its paternal parent than between it and its maternal parent, possibly partly because of cytoplasmic and maternal effects. Specific functional categories were enriched among the 2,545 nonadditive genes of Brassica hexaploid compared with the additive genes; the categories included response to stimulus, immune system process, cellular process, metabolic process, rhythmic process, and pigmentation. Many transcription factor genes, methyltransferases, and methylation genes showed differential expression between Brassica hexaploid and its parents. Our results demonstrate that the Brassica allohexaploid can generate extensive transcriptomic diversity compared with its parents. These changes may contribute to the normal growth and reproduction of allohexaploids.

Analysis and annotation of the hexaploid oat seed transcriptome

BACKGROUND

Next generation sequencing provides new opportunities to explore transcriptomes. However, challenges remain for accurate differentiation of homoeoalleles and paralogs, particularly in polyploid organisms with no supporting genome sequence. In this study, RNA-Seq was employed to generate and characterize the first gene expression atlas for hexaploid oat.

RESULTS

The software packages Trinity and Oases were used to produce a transcript assembly from nearly 134 million 100-bp paired-end reads from developing oat seeds. Based on the quality-parameters employed, Oases assemblies were superior. The Oases 67-kmer assembly, denoted dnOST (de novo Oat Seed Transcriptome), is over 55 million nucleotides in length and the average transcript length is 1,043 nucleotides. The 74.8× sequencing depth was adequate to differentiate a large proportion of putative homoeoalleles and paralogs. To assess the robustness of dnOST, we successfully identified gene transcripts associated with the biosynthetic pathways of three compounds with health-promoting properties (avenanthramides, tocols, β-glucans), and quantified their expression.

CONCLUSIONS

To our knowledge, this study provides the first direct performance comparison between two major assemblers in a polyploid organism. The workflow we developed provides a useful guide for comparable analyses in other organisms. The transcript assembly developed here is a major advance. It expands the number of oat ESTs 3-fold, and constitutes the first comprehensive transcriptome study in oat. This resource will be a useful new tool both for analysis of genes relevant to nutritional enhancement of oat, and for improvement of this crop in general.

Adaptation to seasonality and the winter freeze

Flowering plants initially diversified during the Mesozoic era at least 140 million years ago in regions of the world where temperate seasonal environments were not encountered. Since then several cooling events resulted in the contraction of warm and wet environments and the establishment of novel temperate zones in both hemispheres. In response, less than half of modern angiosperm families have members that evolved specific adaptations to cold seasonal climates, including cold acclimation, freezing tolerance, endodormancy, and vernalization responsiveness. Despite compelling evidence for multiple independent origins, the level of genetic constraint on the evolution of adaptations to seasonal cold is not well understood. However, the recent increase in molecular genetic studies examining the response of model and crop species to seasonal cold offers new insight into the evolutionary lability of these traits. This insight has major implications for our understanding of complex trait evolution, and the potential role of local adaptation in response to past and future climate change. In this review, we discuss the biochemical, morphological, and developmental basis of adaptations to seasonal cold, and synthesize recent literature on the genetic basis of these traits in a phylogenomic context. We find evidence for multiple genetic links between distinct physiological responses to cold, possibly reinforcing the coordinated expression of these traits. Furthermore, repeated recruitment of the same or similar ancestral pathways suggests that land plants might be somewhat pre-adapted to dealing with temperature stress, perhaps making inducible cold traits relatively easy to evolve.

RNA-Seq Assembly - Are We There Yet?

Transcriptomic sequence resources represent invaluable assets for research, in particular for non-model species without a sequenced genome. To date, the Next Generation Sequencing technologies 454/Roche and Illumina have been used to generate transcriptome sequence databases by mRNA-Seq for more than fifty different plant species. While some of the databases were successfully used for downstream applications, such as proteomics, the assembly parameters indicate that the assemblies do not yet accurately reflect the actual plant transcriptomes. Two different assembly strategies have been used, overlap consensus based assemblers for long reads and Eulerian path/de Bruijn graph assembler for short reads. In this review, we discuss the challenges and solutions to the transcriptome assembly problem. A list of quality control parameters and the necessary scripts to produce them are provided.