Analysis of the differential expression of the genes related to Brassica napus seed development

Transcriptome analysis reveals important candidate genes involved in grain-size formation at the stage of grain enlargement in common wheat cultivar "Bainong 4199"

Grain-size is one of the yield components, and the first 14 days after pollination (DAP) is a crucial stage for wheat grain-size formation. To understand the mechanism of grain-size formation at the whole gene expression level and to identify the candidate genes related to grain pattern formation, cDNA libraries from immature grains of 5 DAP and 14 DAP were constructed. According to transcriptome analysis, a total of 12,555 new genes and 9,358 differentially expressed genes (DEGs) were obtained. In DEGs, 2,876, 3,357 and 3,125 genes were located on A, B and D subgenome respectively. 9,937 (79.15%) new genes and 9,059 (96.80%) DEGs were successfully annotated. For DEGs, 4,453 were up-regulated and 4,905 were down-regulated at 14 DAP. The Gene Ontology (GO) database indicated that most of the grain-size-related genes were in the same cluster. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database analysis showed that 130, 129 and 20 DEGs were respectively involved in starch and sucrose metabolism, plant hormone signal transduction and ubiquitin-mediated proteolysis. Expression levels of 8 randomly selected genes were confirmed by qRT-PCR, which was consistent with the transcriptome data. The present database will help us understand the molecular mechanisms underlying early grain development and provide the foundation for increasing grain-size and yield in wheat breeding programs.

Combined QTL mapping, physiological and transcriptomic analyses to identify candidate genes involved in Brassica napus seed aging

Seed aging is an inevitable problem in the germplasm conservation of oil crops. Thus, clarifying the genetic mechanism of seed aging is important for rapeseed breeding. In this study, Brassica napus seeds were exposed to an artificial aging environment (40 °C and 90% relative humidity). Using a population of 172 recombinant inbred lines, 13 QTLs were detected on 8 chromosomes, which explained ~ 9.05% of the total phenotypic variation. The QTLs q2015AGIA-C08 and q2016AGI-C08-2 identified in the two environments were considered the same QTL. After artificial aging, lower germination index, increased relative electrical conductivity, malondialdehyde and proline content, and reduced soluble sugar, protein content and antioxidant enzyme activities were detected. Furthermore, seeds of extreme lines that were either left untreated (R0 and S0) or subjected to 15 days of artificial aging (R15 and S15) were used for transcriptome sequencing. In total, 2843, 1084, 429 and 1055 differentially expressed genes were identified in R15 vs. R0, S15 vs. S0, R0 vs. S0 and R15 vs. S15, respectively. Through integrated QTL mapping and RNA-sequencing analyses, seven genes, such as BnaA03g37460D, encoding heat shock transcription factor C1, and BnaA03g40360D, encoding phosphofructokinase 4, were screened as candidate genes involved in seed aging. Further researches on these candidate genes could broaden our understanding of the regulatory mechanisms of seed aging.

Identification of genomic ATP binding cassette (ABC) transporter genes and Cd-responsive ABCs in Brassica napus

The plant ATP binding cassette (ABC) transporters are one of the integral membrane proteins responsible for uptake and allocation of a wide range of metabolites and xenobiotics including heavy metals (e.g. zinc, manganese and cadmium). They play multiple roles in plant growth, development and environmental stress responses. Although the ABC transporters have been identified in model plants such as Arabidopsis and rice, they are have not been annotated and identified in rapeseed (Brassica napus) and also, little is known about functionality of these metal transporters. B. napus is an important oil crop ranking the third largest source of vegetable oil worldwide. Importantly, it is long considered as a desirable candidate for phytoremediation owning to its massive dry weight productivity and moderate Cd accumulation. In this study, we identified 314 ABC protein genes from B. napus using bioinformatics and high-throughput sequencing. Eight subfamilies including ABCA-G and ABCI have been categorized. The ABCG proteins constitute the largest subfamily with 116 members, and the ABCB and ABCC subfamilies ranks second and third with 69 and 47 members, respectively. Analyses of ABCs in B. napus genome reveal that their evolutional expansion was through localized allele duplications. Most of the ABC genes (74.2%, 233/314) were validated by RNA-sequencing rapeseed seedlings. Among the 233 profiled BnaABCs, 132 genes were differentially expressed (>1.5 fold change, p < 0.05) and 84 genes were significantly induced under Cd stress. Analyses of specific cis-elements in the upstream of eight representative genes show diverse motifs, which potentially respond to environmental stress, hormone responsiveness and other development signals.

Gene expression profiling during seed-filling process in peanut with emphasis on oil biosynthesis networks

Pod-filling is an important stage of peanut (Arachis hypogaea) seed development. It is partially controlled by genetic factors, as cultivars considerably vary in pod-filling potential. Here, a study was done to detect changes in mRNA levels that accompany pod-filling processes. Four seed developmental stages were sampled from two peanut genotypes differing in their oil content and pod-filling potential. Transcriptome data were generated by RNA-Seq and explored with respect to genic and subgenomic patterns of expression. Very dynamic transcriptomic changes occurred during seed development in both genotypes. Yet, general higher expression rates of transcripts and an enrichment in processes involved "energy generation" and "primary metabolites" were observed in the genotype with the better pod-filling ("Hanoch"). A dataset of 584 oil-related genes was assembled and analyzed, resulting in several lipid metabolic processes highly expressed in Hanoch, including oil storage and FA synthesis/elongation. Homoeolog-specific gene expression analysis revealed that both subgenomes contribute to the oil genes expression. Yet, biases were observed in particular parts of the pathway with possible biological meaning, presumably explaining the genotypic variation in oil biosynthesis and pod-filling. This study provides baseline information and a resource that may be used to understand development and oil biosynthesis in the peanut seeds.

BOG: R-package for Bacterium and virus analysis of Orthologous Groups

BOG (Bacterium and virus analysis of Orthologous Groups) is a package for identifying groups of differentially regulated genes in the light of gene functions for various virus and bacteria genomes. It is designed to identify Clusters of Orthologous Groups (COGs) that are enriched among genes that have gone through significant changes under different conditions. This would contribute to the detection of pathogens, an important scientific research area of relevance in uncovering bioterrorism, among others. Particular statistical analyses include hypergeometric, Mann–Whitney rank sum, and gene set enrichment. Results from the analyses are organized and presented in tabular and graphical forms for ease of understanding and dissemination of results. BOG is implemented as an R-package, which is available from CRAN or can be downloaded from http://www.stat.osu.edu/~statgen/SOFTWARE/BOG/.

Transcriptomic analysis revealed the mechanism of oil dynamic accumulation during developing Siberian apricot (Prunus sibirica L.) seed kernels for the development of woody biodiesel

BACKGROUND

Siberian apricot (Prunus sibirica L.) has emerged as a novel potential source of biodiesel in China, but the molecular regulatory mechanism of oil accumulation in Siberian apricot seed kernels (SASK) is still unknown at present. To better develop SASK oil as woody biodiesel, it is essential to profile transcriptome and to identify the full repertoire of potential unigenes involved in the formation and accumulation of oil SASK during the different developing stages.

RESULTS

We firstly detected the temporal patterns for oil content and fatty acid (FA) compositions of SASK in 7 different developing stages. The best time for obtaining the high quality and quantity of SASK oil was characterized at 60 days after flowering (DAF), and the representative periods (10, 30, 50, 60, and 70 DAF) were selected for transcriptomic analysis. By Illumina/Solexa sequencings, approximately 65 million short reads (average length = 96 bp) were obtained, and then assembled into 124,070 unigenes by Trinity strategy (mean size = 829.62 bp). A total of 3,000, 2,781, 2,620, and 2,675 differentially expressed unigenes were identified at 30, 50, 60, and 70 DAF (10 DAF as the control) by DESeq method, respectively. The relationship between the unigene transcriptional profiles and the oil dynamic patterns in developing SASK was comparatively analyzed, and the specific unigenes encoding some known enzymes and transcription factors involved in acetyl-coenzyme A (acetyl-CoA) formation and oil accumulation were determined. Additionally, 5 key metabolic genes implicated in SASK oil accumulation were experimentally validated by quantitative real-time PCR (qRT-PCR). Our findings could help to construction of oil accumulated pathway and to elucidate the molecular regulatory mechanism of increased oil production in developing SASK.

CONCLUSIONS

This is the first study of oil temporal patterns, transcriptome sequencings, and differential profiles in developing SASK. All our results will serve as the important foundation to further deeply explore the regulatory mechanism of SASK high-quality oil accumulation, and may also provide some reference for researching the woody biodiesel plants.

Identification of desiccation tolerance transcripts potentially involved in rape (Brassica napus L.) seeds development and germination

To investigate regulatory processes and protective mechanisms leading to desiccation tolerance (DT) in seeds, cDNA amplified fragment length polymorphism (cDNA-AFLP) in conjunction with 128 primer combinations was used to detect differential gene expression in rape seeds in response to DT during seed development and germination. We obtained approximately 8000 transcript-derived fragments (TDFs), of which 394 TDFs with differential expression patterns ("sustained expression", "up-regulated", "couple with seed DT", and "down-regulated") were excised from gels and re-amplified by polymerase chain reaction (PCR). After sequencing and comparison with the National Center for Biotechnology Information database, 176 TDFs presented significant similarity with known genes that could be classified into the following categories: metabolism and energy, stress resistance and defense, storage, signal transduction, and other functional categories. Using semiquantitative reverse-transcription PCR and real-time PCR approaches, the significance of the differences was further confirmed in fresh seeds and dehydrated seeds. The genes that encode superoxide dismutase, peroxiredoxin, caleosin, oleosin S3, steroleosin, late embryogenesis abundant protein, glutathione reductase, β-glucosidase, S23 transcriptional repressor, and some heat-shock proteins could be associated with DT. The results of this study will aid in the identification of candidate genes for future experiments that seek to understand seed DT.

Choreography of Transcriptomes and Lipidomes of Nannochloropsis Reveals the Mechanisms of Oil Synthesis in Microalgae

To reveal the molecular mechanisms of oleaginousness in microalgae, transcriptomic and lipidomic dynamics of the oleaginous microalga Nannochloropsis oceanica IMET1 under nitrogen-replete (N+) and N-depleted (N-) conditions were simultaneously tracked. At the transcript level, enhanced triacylglycerol (TAG) synthesis under N- conditions primarily involved upregulation of seven putative diacylglycerol acyltransferase (DGAT) genes and downregulation of six other DGAT genes, with a simultaneous elevation of the other Kennedy pathway genes. Under N- conditions, despite downregulation of most de novo fatty acid synthesis genes, the pathways that shunt carbon precursors from protein and carbohydrate metabolic pathways into glycerolipid synthesis were stimulated at the transcript level. In particular, the genes involved in supplying carbon precursors and energy for de novo fatty acid synthesis, including those encoding components of the pyruvate dehydrogenase complex (PDHC), glycolysis, and PDHC bypass, and suites of specific transporters, were substantially upregulated under N- conditions, resulting in increased overall TAG production. Moreover, genes involved in the citric acid cycle and β-oxidation in mitochondria were greatly enhanced to utilize the carbon skeletons derived from membrane lipids and proteins to produce additional TAG or its precursors. This temporal and spatial regulation model of oil accumulation in microalgae provides a basis for improving our understanding of TAG synthesis in microalgae and will also enable more rational genetic engineering of TAG production.

Identification of key drought stress-related genes in the hyacinth bean

Hyacinth bean (Lablab purpureus [Linn.] Sweet) possesses excellent characteristics for field production, but the response of this plant to drought stress has not been described at the molecular level. Suppression subtraction hybridization (SSH) is an effective way to exploit key factors for plant responses to drought stress that are involved in transcriptional and metabolic activities. In this study, forward and reverse SSH libraries were generated from root tissues of the drought-tolerant hyacinth bean genotype MEIDOU 2012 under water-stress conditions. A total of 1,287 unigenes (94 contigs and 1,193 singletons) were derived from sequence alignment and cluster assembly of 1400 ESTs, and 80.6% of those hit against NCBI non-redundant (nr) database with E value <1E-06. BLASTX analysis revealed that the majority top matches were proteins form Glycine max (L.) Merrill. (61.5%). According to a gene ontology (GO) functional classification, 816 functionally annotated unigenes were assigned to the biological process category (74.1%), and 83.9% of them classified into molecular function and 69.2% involved in cellular component. A total of 168 sequences were further annotated with 207 Enzyme Commission (EC) codes and mapped to 83 different KEGG pathways. Seventeen functionally relevant genes were found to be overrepresented under drought stress using enrichment analysis. Differential expression of unigenes were confirmed by quantitative real-time PCR assays, and their transcript profiles generally divided into three patterns, depending on the expression peaked levels after 6, 8 or 10 days dehydration, which indicated that these genes are functionally associated in the drought-stress response.

Differential transcript profiling through cDNA-AFLP showed complexity of rutin biosynthesis and accumulation in seeds of a nutraceutical food crop (Fagopyrum spp.)

Background

Buckwheat, consisting of two cultivated species Fagopyrum tataricum and F. esculentum, is the richest source of flavonoid rutin. Vegetative tissues of both the Fagopyrum species contain almost similar amount of rutin; however, rutin content in seed of F. tataricum are ~50 folds of that in seed of F. esculentum. In order to understand the molecular basis of high rutin content in F. tataricum, differential transcript profiling through cDNA-AFLP has been utilized to decipher what genetic factors in addition to flavonoid structural genes contribute to high rutin content of F. tataricum compared to F. esculentum.

Results

Differential transcript profiling through cDNA-AFLP in seed maturing stages (inflorescence to seed maturation) with 32 primer combinations generated total of 509 transcript fragments (TDFs). 167 TDFs were then eluted, cloned and sequenced from F. tataricum and F. esculentum. Categorization of TDFs on the basis of their presence/absence (qualitative variation) or differences in the amount of expression (quantitative variation) between both the Fagopyrum species showed that majority of variants are quantitative (64%). The TDFs represented genes controlling different biological processes such as basic and secondary metabolism (33%), regulation (18%), signal transduction (14%), transportation (13%), cellular organization (10%), and photosynthesis & energy (4%). Most of the TDFs except belonging to cellular metabolism showed relatively higher transcript abundance in F. tataricum over F. esculentum. Quantitative RT-PCR analysis of nine TDFs representing genes involved in regulation, metabolism, signaling and transport of secondary metabolites showed that all the tested nine TDFs (Ubiquitin protein ligase, ABC transporter, sugar transporter) except MYB 118 showed significantly higher expression in early seed formation stage (S7) of F. tataricum compared to F. esculentum. qRT-PCR results were found to be consistent with the cDNA-AFLP results.

Conclusions

The present study concludes that in addition to structural genes, other classes of genes such as regulators, modifiers and transporters are also important in biosynthesis and accumulation of flavonoid content in plants. cDNA-AFLP technology was successfully utilized to capture genes that are contributing to differences in rutin content in seed maturing stages of Fagopyrum species. Increased transcript abundance of TDFs during transition from flowers to seed maturation suggests their involvement not only in the higher rutin content of F. tataricum over F. esculentum but also in nutritional superiority of the former.

Large-scale sequencing of normalized full-length cDNA library of soybean seed at different developmental stages and analysis of the gene expression profiles based on ESTs

Although GenBank has now covered over 1,400,000 expressed sequence tags (ESTs) from soybean, most ESTs available to the public have been derived from tissues or environmental conditions rather than developing seeds. It is absolutely necessary for annotating the molecular mechanisms of soybean seed development to analyze completely the gene expression profiles of its immature seed at various stages. Here we have constructed a full-length-enriched cDNA library comprised of a total of 45,408 cDNA clones which cover various stages of soybean seed development. Furthermore, we have sequenced from 5' ends of these clones, 36,656 ESTs were obtained in the present study. These EST sequences could be categorized into 27,982 unigenes, including 22,867 contigs and 5,115 singletons, among which 27,931 could be mapped onto soybean 20 chromosome sequences. Comparative genomic analysis with other plants has revealed that these unigenes include lots of candidate genes specific to dicot, legume and soybean. Approximately 1,789 of these unigenes currently show no homology to known soybean sequences, suggesting that many represent mRNAs specifically expressed in seeds. Novel abundant genes involved in the oil synthesis have been found in this study, may serve as a valuable resource for soybean seed improvement.

Comparative deep transcriptional profiling of four developing oilseeds

Transcriptome analysis based on deep expressed sequence tag (EST) sequencing allows quantitative comparisons of gene expression across multiple species. Using pyrosequencing, we generated over 7 million ESTs from four stages of developing seeds of Ricinus communis, Brassica napus, Euonymus alatus and Tropaeolum majus, which differ in their storage tissue for oil, their ability to photosynthesize and in the structure and content of their triacylglycerols (TAG). The larger number of ESTs in these 16 datasets provided reliable estimates of the expression of acyltransferases and other enzymes expressed at low levels. Analysis of EST levels from these oilseeds revealed both conserved and distinct species-specific expression patterns for genes involved in the synthesis of glycerolipids and their precursors. Independent of the species and tissue type, ESTs for core fatty acid synthesis enzymes maintained a conserved stoichiometry and a strong correlation in temporal profiles throughout seed development. However, ESTs associated with non-plastid enzymes of oil biosynthesis displayed dissimilar temporal patterns indicative of different regulation. The EST levels for several genes potentially involved in accumulation of unusual TAG structures were distinct. Comparison of expression of members from multi-gene families allowed the identification of specific isoforms with conserved function in oil biosynthesis. In all four oilseeds, ESTs for Rubisco were present, suggesting its possible role in carbon metabolism, irrespective of light availability. Together, these data provide a resource for use in comparative and functional genomics of diverse oilseeds. Expression data for more than 350 genes encoding enzymes and proteins involved in lipid metabolism are available at the 'ARALIP' website (http://aralip.plantbiology.msu.edu/).

Differential transcript accumulation in chickpea during early phases of compatible interaction with a necrotrophic fungus Ascochyta rabiei

The initial phases of the disease establishment are very crucial for the compatible interactions. Pathogens must overcome the responses generated by the host for the onset of disease invasion. The compatible interaction is inadequately represented in plant-pathogen interaction studies. To gain broader insight into the early responses elicited by chickpea blight fungus Ascochyta rabiei during compatible interaction; we isolated early responsive genes of chickpea using PCR based suppression subtractive hybridization (SSH) strategy. We obtained ~250 unique genes after homology search and redundancy elimination. Based on their potential cellular functions, these genes were broadly classified into eleven different categories viz. stress, signaling, gene regulation, cellular metabolism and genes of unknown functions. Present study revealed few unexpected genes which have a possible role in induced immunity and disease progression. We employed macroarray, northern blot, real-time PCR and cluster analysis to develop transcript profiles. Most of the genes analyzed were early induced and were transcriptionally upregulated upon 24 h post inoculation. Our approach has rendered the isolation of early responsive genes involved in signaling and regulation of metabolic changes upon fungal infection. The information obtained will help to dissect the molecular mechanisms during compatible chickpea-Ascochyta interactions.