Identification of differentially-expressed genes associated with pistil abortion in Japanese apricot by genome-wide transcriptional analysis

PmAGAMOUS recruits polycomb protein PmLHP1 to regulate single-pistil morphogenesis in Japanese apricot

Japanese apricot (Prunus mume Sieb. et Zucc.) is a traditional fruit tree with a long history. Multiple pistils (MP) lead to the formation of multiple fruits, decreasing fruit quality and yield. In this study, the morphology of flowers was observed at 4 stages of pistil development: undifferentiated stage (S1), predifferentiation stage (S2), differentiation stage (S3), and late differentiation stage (S4). In S2 and S3, the expression of PmWUSCHEL (PmWUS) in the MP cultivar was significantly higher than that in the single-pistil (SP) cultivar, and the gene expression of its inhibitor, PmAGAMOUS (PmAG), also showed the same trend, indicating that other regulators participate in the regulation of PmWUS during this period. Chromatin immunoprecipitation-qPCR (ChIP-qPCR) showed that PmAG could bind to the promoter and the locus of PmWUS, and H3K27me3 repressive marks were also detected at these sites. The SP cultivar exhibited an elevated level of DNA methylation in the promoter region of PmWUS, which partially overlapped with the region of histone methylation. This suggests that the regulation of PmWUS involves both transcription factors and epigenetic modifications. Also, the gene expression of Japanese apricot LIKE HETEROCHROMATIN PROTEIN (PmLHP1), an epigenetic regulator, in MP was significantly lower than that in SP in S2 to 3, contrary to the trend in expression of PmWUS. Our results showed that PmAG recruited sufficient PmLHP1 to maintain the level of H3K27me3 on PmWUS during the S2 of pistil development. This recruitment of PmLHP1 by PmAG inhibits the expression of PmWUS at the precise time, leading to the formation of 1 normal pistil primordium.

NBLAST: a graphical user interface-based two-way BLAST software with a dot plot viewer

BLAST, a basic bioinformatics tool for searching local sequence similarity, has been one of the most widely used bioinformatics programs since its introduction in 1990. Users generally use the web-based NCBI-BLAST program for BLAST analysis. However, users with large sequence data are often faced with a problem of upload size limitation while using the web-based BLAST program. This proves inconvenient as scientists often want to run BLAST on their own data, such as transcriptome or whole genome sequences. To overcome this issue, we developed NBLAST, a graphical user interface-based BLAST program that employs a two-way system, allowing the use of input sequences either as "query" or "target" in the BLAST analysis. NBLAST is also equipped with a dot plot viewer, thus allowing researchers to create custom database for BLAST and run a dot plot similarity analysis within a single program. It is available to access to the NBLAST with http://nbitglobal.com/nblast.

Transcriptome analysis of flower bud identified genes associated with pistil abortions between long branches and spur twigs in apricots (Prunus armeniaca L.)

Pistil abortions of flower buds occur frequently in many apricot cultivars, especially in long branches. However, the molecular mechanism underlying pistil abortion in apricots remains unclear. To better understand the molecular mechanism of pistil abortions between long branches and spur twigs, paraffin sections and high-throughput sequencing technology were employed to analyze the expression patterns of genes associated with pistil abortions during later flower bud development stage in ‘Shajinhong’ apricot. The result of stage III (separation of bud scales) was the critical stage of pistil abortion in apricots. A total of 163 differentially expressed genes were identified as candidate genes related to pistil abortion in long branches. These genes are implicated in programmed cell death, hormone signaling, cell wall degeneration, and the carbohydrate metabolism pathway. The results showed that the up-regulation of gene expression of Xyloglucan endotransglucosylase/hydrolase and β-glucosidase in flower buds might be the direct cause of cell wall breakdown and pistil necrosis in long branches. We hypothesize that there is a molecular relationship between pistil abortion before blooming and cellulose degradation, and then carbohydrate transport in the case of carbon deficiency in long branches. Our work provides new insights into cellulose degradation in abortion pistils and valuable information on flower development in apricots, and also provides a useful reference for cultivation regulation in apricot or other fruit crops.

Integrative Analyses of Transcriptomics and Metabolomics in Sex Differentiation of Mulberry Flowers

Sex determination and sex differentiation of plants are important physiological processes of plant development. Mulberry (Morus indica L.) is an important economic tree being cultivated in sericulture countries, and mulberry leaf is commonly used for sericulture. The transcriptomic and metabolomic differences between the staminate flowers (SFs) and pistillate flowers (PFs) of mulberry were investigated by RNA sequencing and ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Overall, we uncovered 4,230 genes and 209 metabolites are significantly differentially expressed between the SFs and PFs of mulberry. The combined transcriptomic and metabolomic analysis revealed these differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) are involved in flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism, and these detected DEGs and DEMs may be associated with sex differentiation of mulberry through the regulation of the enrichment pathways, such as the MAPK pathway, flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism. This study will provide a rich source for the analysis of the molecular mechanism of mulberry sex differentiation processes.

Dosage-Dependent Gynoecium Development and Gene Expression in Brassica napus-Orychophragmus violaceus Addition Lines

Distant hybridization usually leads to female sterility of the hybrid but the mechanism behind this is poorly understood. Complete pistil abortion but normal male fertility was shown by one Brassica napus-Orychophragmus violaceus monosomic alien addition line (MA, AACC + 1 IO, 2n = 39) produced previously. To study the effect of a single O. violaceus chromosome addition on pistil development in different genetic backgrounds, hybrids between the MA and B. carinata (BBCC), B. juncea (AABB), and two synthetic hexaploids (AABBCC) were firstly produced in this study which show complete female sterility. A microspore culture was further performed to produce the haploid monosomic alien addition line (HMA, AC + 1 IO, 2n = 20) and disomic addition line (DA, AACC + 2 IO, 2n = 40) together with haploid (H, AC, 2n = 19) and double haploid (DH, AACC, 2n = 38) plants of B. napus from MA to investigate the dosage effect of the alien O. violaceus chromosome on pistil development and gene expression. Compared to MA, the development of the pistils of DA and HMA was completely or partially recovered, in which the pistils could swell and elongate to a normal shape after open pollination, although no seeds were produced. Comparative RNA-seq analyses revealed that the numbers of the differentially expressed genes (DEGs) were significantly different, dosage-dependent, and consistent with the phenotypic difference in pairwise comparisons of HMA vs. H, DA vs. DH, MA vs. DH, MA vs. DA, and MA vs. HMA. The gene ontology (GO) enrichment analysis of DEGs showed that a number of genes involved in the development of the gynoecium, embryo sac, ovule, and integuments. Particularly, several common DEGs for pistil development shared in HMA vs. H and DA vs. DH showed functions in genotoxic stress response, auxin transport, and signaling and adaxial/abaxial axis specification. The results provided updated information for the molecular mechanisms behind the gynoecium development of B. napus responding to the dosage of alien O. violaceus chromosomes.

Genome-wide association analysis for yield-related traits at the R6 stage in a Chinese soybean mini core collection

BACKGROUND

Soybean (Glycine max (L.) Merr.) is an economically important crop for vegetable oil and protein production, and yield is a critical trait for grain/vegetable uses of soybean. However, our knowledge of the genes controlling the vegetable soybean yield remains limited.

OBJECTIVE

To better understand the genetic basis of the vegetable soybean yield.

METHODS

The 100-pod fresh weight (PFW), 100-seed fresh weight (SFW), kernel percent (KP) and moisture content of fresh seeds (MCFS) at the R6 stage are four yield-related traits for vegetable soybean. We investigated a soybean mini core collection composed of 224 germplasm accessions for four yield-related traits in two consecutive years. Based on 1514 single nucleotide polymorphisms (SNPs), genome-wide association studies (GWAS) were conducted using a mixed linear model (MLM).

RESULTS

Extensive phenotypic variation existed in the soybean mini core collection and significant positive correlations were shown among most of traits. A total of 16 SNP markers for PFW, SFW, KP and MCFS were detected in all environments via GWAS. Nine SNP markers were repeatedly identified in two environments. Among these markers, eight were located in or near regions where yield-related QTLs have been reported in previous studies, and one was a novel genetic locus identified in this study. In addition, we conducted candidate gene analysis to the large-effect SNP markers, a total of twelve genes were proposed as potential candidate genes of soybean yield at the R6 stage.

CONCLUSION

These results will be beneficial for understanding the genetic basis of soybean yield at the R6 stage and facilitating the pyramiding of favourable alleles for future high-yield breeding by marker-assisted selection in vegetable soybean.

Genome-wide analysis of PmTCP4 transcription factor binding sites by ChIP-Seq during pistil abortion in Japanese apricot

The TCP4 transcription factor plays an important role in plant growth and development, especially in flower development. PmTCP4 is involved in the process of pistil abortion in Japanese apricot, but its molecular mechanism, particularly the DNA binding sites and co-regulatory genes, are quite unknown. Therefore, to identify the genome-wide binding sites of PmTCP4 transcription factors and their co-regulatory genes, chromatin immunoprecipitation sequencing (ChIP-Seq) was carried out. ChIP-Seq data produced the maximum enriched peaks in two Japanese apricot cultivars 'Daqiandi' (DQD) and 'Longyan' (LY), which showed that the majority of DNA-protein interactions are relevant and have a significant function in binding sites. Moreover, 720 and 251 peak-associated genes regulated by PmTCP4 were identified in DQD and LY, respectively, and most of them were involved in the flower and pistil development process. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that photosynthesis and oxidative phosphorylation were the most enriched pathways in both cultivars and all identified genes related to these pathways were down-regulated. This study will provide a reference for a better understanding of the PmTCP4 regulatory mechanism during pistil abortion in Japanese apricot.

Photosynthesis and yield response to elevated CO2, C4 plant foxtail millet behaves similarly to C3 species

Foxtail millet (Setaria italica) is a nutrient-rich food source traditionally grown in arid and semi-arid areas, as it is well adapted to drought climate. Yet there is limited information as how the crop responses to the changing climate. In order to investigate the response of foxtail millet to elevated [CO₂] and the underlying mechanism, the crop was grown at ambient [CO₂] (400 μmol mol^-1) and elevated [CO₂] (600 μmol mol^-1) in an open-top chamber (OTC) experimental facility in North China. The changes in leaf photosynthesis, chlorophyll fluorescence, biomass, yield and global gene expression in response to elevated [CO₂] were determined. Despite foxtail millet being a C₄ photosynthetic crop, photosynthetic rates (P_N) and intrinsic water-use efficiency (WUEi), were increased under elevated [CO₂]. Similarly, grain yield and above-ground biomass also significantly increased (P <  0.05) for the two years of experimentation under elevated [CO₂]. Increases in seeds and tiller number, spike and stem weight were the main contributors to the increased grain yield and biomass. Using transcriptomic analyses, this study further identified some genes which play a role in cell wall reinforcement, shoot initiation, stomatal conductance, carbon fixation, glycolysis / gluconeogenesis responsive to elevated [CO₂]. Changes in these genes reduced plant height, increased stem diameters, and promote CO₂ fixation. Higher photosynthetic rates at elevated [CO₂] demonstrated that foxtail millet was not photosynthetically saturated at elevated [CO₂] and its photosynthesis response to elevated [CO₂] were analogous to C₃ plants.

Transcriptional regulatory networks controlling taste and aroma quality of apricot (Prunus armeniaca L.) fruit during ripening

BACKGROUND

Taste and aroma, which are important organoleptic qualities of apricot (Prunus armeniaca L.) fruit, undergo rapid and substantial changes during ripening. However, the associated molecular mechanisms remain unclear. The goal of this study was to identify candidate genes for flavor compound metabolism and to construct a regulatory transcriptional network.

RESULTS

We characterized the transcriptome of the 'Jianali' apricot cultivar, which exhibits substantial changes in flavor during ripening, at 50 (turning), 73 (commercial maturation) and 91 (full ripe) days post anthesis (DPA) using RNA sequencing (RNA-Seq). A weighted gene co-expression network analysis (WGCNA) revealed that four of 19 modules correlated highly with flavor compound metabolism (P < 0.001). From them, we identified 1237 differentially expressed genes, with 16 intramodular hubs. A proposed pathway model for flavor compound biosynthesis is presented based on these genes. Two SUS1 genes, as well as SPS2 and INV1 were correlated with sugar biosynthesis, while NADP-ME4, two PK-like and mitochondrial energy metabolism exerted a noticeable effect on organic acid metabolism. CCD1 and FAD2 were identified as being involved in apocarotenoid aroma volatiles and lactone biosynthesis, respectively. Five sugar transporters (Sweet10, STP13, EDR6, STP5.1, STP5.2), one aluminum-activated malate transporter (ALMT9) and one ABCG transporter (ABCG11) were associated with the transport of sugars, organic acids and volatiles, respectively. Sixteen transcription factors were also highlighted that may also play regulatory roles in flavor quality development.

CONCLUSIONS

Apricot RNA-Seq data were obtained and used to generate an annotated set of predicted expressed genes, providing a platform for functional genomic research. Using network analysis and pathway mapping, putative molecular mechanisms for changes in apricot fruit taste and aroma during ripening were elucidated.

RNA-seq reveals transcriptome changes in goats following myostatin gene knockout

Myostatin (MSTN) is a powerful negative regulator of skeletal muscle mass in mammalian species that is primarily expressed in skeletal muscles, and mutations of its encoding gene can result in the double-muscling trait. In this study, the CRISPR/Cas9 technique was used to edit MSTN in Shaanbei Cashmere goats and generate knockout animals. RNA sequencing was used to determine and compare the transcriptome profiles of the muscles from three wild-type (WT) goats, three fibroblast growth factor 5 (FGF5) knockout goats (FGF5+/- group) and three goats with disrupted expression of both the FGF5 and MSTN genes (FM+/- group). The sequence reads were obtained using the Illumina HiSeq 2000 system and mapped to the Capra hircus reference genome using TopHat (v2.0.9). In total, 68.93, 62.04 and 66.26 million clean sequencing reads were obtained from the WT, FM+/- and FGF5+/- groups, respectively. There were 201 differentially expressed genes (DEGs) between the WT and FGF5+/- groups, with 86 down- and 115 up-regulated genes in the FGF5+/- group. Between the WT and FM+/- groups, 121 DEGs were identified, including 81 down- and 40 up-regulated genes in the FM+/- group. A total of 198 DEGs were detected between the FGF5+/- group and FM+/- group, with 128 down- and 70 up-regulated genes in the FM+/- group. At the transcriptome level, we found substantial changes in genes involved in fatty acid metabolism and the biosynthesis of unsaturated fatty acids, such as stearoyl-CoA dehydrogenase, 3-hydroxyacyl-CoA dehydratase 2, ELOVL fatty acid elongase 6 and fatty acid synthase, suggesting that the expression levels of these genes may be directly regulated by MSTN and that these genes are likely downstream targets of MSTN with potential roles in lipid metabolism in goats. Moreover, five randomly selected DEGs were further validated with qRT-PCR, and the results were consistent with the transcriptome analysis. The present study provides insight into the unique transcriptome profile of the MSTN knockout goat, which is a valuable resource for studying goat genomics.

De novo RNA sequencing transcriptome of Rhododendron obtusum identified the early heat response genes involved in the transcriptional regulation of photosynthesis

Rhododendron spp. is an important ornamental species that is widely cultivated for landscape worldwide. Heat stress is a major obstacle for its cultivation in south China. Previous studies on rhododendron principally focused on its physiological and biochemical processes, which are involved in a series of stress tolerance. However, molecular or genetic properties of rhododendron’s response to heat stress are still poorly understood. The phenotype and chlorophyll fluorescence kinetics parameters of four rhododendron cultivars were compared under normal or heat stress conditions, and a cultivar with highest heat tolerance, “Yanzhimi” (R. obtusum) was selected for transcriptome sequencing. A total of 325,429,240 high quality reads were obtained and assembled into 395,561 transcripts and 92,463 unigenes. Functional annotation showed that 38,724 unigenes had sequence similarity to known genes in at least one of the proteins or nucleotide databases used in this study. These 38,724 unigenes were categorized into 51 functional groups based on Gene Ontology classification and were blasted to 24 known cluster of orthologous groups. A total of 973 identified unigenes belonged to 57 transcription factor families, including the stress-related HSF, DREB, ZNF, and NAC genes. Photosynthesis was significantly enriched in the Kyoto Encyclopedia of Genes and Genomes pathway, and the changed expression pattern was illustrated. The key pathways and signaling components that contribute to heat tolerance in rhododendron were revealed. These results provide a potentially valuable resource that can be used for heat-tolerance breeding.

De novo sequencing and comparative transcriptome analysis of the male and hermaphroditic flowers provide insights into the regulation of flower formation in andromonoecious taihangia rupestris

BACKGROUND

Taihangia rupestris, an andromonoecious plant species, bears both male and hermaphroditic flowers within the same individual. However, the establishment and development of male and hermaphroditic flowers in andromonoecious Taihangia remain poorly understood, due to the limited genetic and sequence information. To investigate the potential molecular mechanism in the regulation of Taihangia flower formation, we used de novo RNA sequencing to compare the transcriptome profiles of male and hermaphroditic flowers at early and late developmental stages.

RESULTS

Four cDNA libraries, including male floral bud, hermaphroditic floral bud, male flower, and hermaphroditic flower, were constructed and sequenced by using the Illumina RNA-Seq method. Totally, 84,596,426 qualified Illumina reads were obtained and then assembled into 59,064 unigenes, of which 24,753 unigenes were annotated in the NCBI non-redundant protein database. In addition, 12,214, 7,153, and 8,115 unigenes were assigned into 53 Gene Ontology (GO) functional groups, 25 Clusters of Orthologous Group (COG) categories, and 126 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. By pairwise comparison of unigene abundance between the samples, we identified 1,668 differential expressed genes (DEGs), including 176 transcription factors (TFs) between the male and hermaphroditic flowers. At the early developmental stage, we found 263 up-regulated genes and 436 down-regulated genes expressed in hermaphroditic floral buds, while 844 up-regulated genes and 314 down-regulated genes were detected in hermaphroditic flowers at the late developmental stage. GO and KEGG enrichment analyses showed that a large number of DEGs were associated with a wide range of functions, including cell cycle, epigenetic processes, flower development, and biosynthesis of unsaturated fatty acid pathway. Finally, real-time quantitative PCR was conducted to validate the DEGs identified in the present study.

CONCLUSION

In this study, transcriptome data of this rare andromonoecious Taihangia were reported for the first time. Comparative transcriptome analysis revealed the significant differences in gene expression profiles between male and hermaphroditic flowers at early and late developmental stages. The transcriptome data of Taihangia would be helpful to improve the understanding of the underlying molecular mechanisms in regulation of flower formation and unisexual flower establishment in andromonoecious plants.

Comparative Transcriptome Analysis Reveals Effects of Exogenous Hematin on Anthocyanin Biosynthesis during Strawberry Fruit Ripening

Anthocyanin in strawberries has a positive effect on fruit coloration. In this study, the role of exogenous hematin on anthocyanin biosynthesis was investigated. Our result showed that the white stage of strawberries treated with exogenous hematin had higher anthocyanin content, compared to the control group. Among all treatments, 5 μM of hematin was the optimal condition to promote color development. In order to explore the molecular mechanism of fruit coloring regulated by hematin, transcriptomes in the hematin- and non-hematin-treated fruit were analyzed. A large number of differentially expressed genes (DEGs) were identified in regulating anthocyanin synthesis, including the DEGs involved in anthocyanin biosynthesis, hormone signaling transduction, phytochrome signaling, starch and sucrose degradation, and transcriptional pathways. These regulatory networks may play an important role in regulating the color process of strawberries treated with hematin. In summary, exogenous hematin could promote fruit coloring by increasing anthocyanin content in the white stage of strawberries. Furthermore, transcriptome analysis suggests that hematin-promoted fruit coloring occurs through multiple related metabolic pathways, which provides valuable information for regulating fruit color via anthocyanin biosynthesis in strawberries.

Identification and Functional Annotation of Genes Differentially Expressed in the Reproductive Tissues of the Olive Tree (Olea europaea L.) through the Generation of Subtractive Libraries

The olive tree is a crop of high socio-economical importance in the Mediterranean area. Sexual reproduction in this plant is an essential process, which determines the yield. Successful fertilization is mainly favored and sometimes needed of the presence of pollen grains from a different cultivar as the olive seizes a self-incompatibility system allegedly determined of the sporophytic type. The purpose of the present study was to identify key gene products involved in the function of olive pollen and pistil, in order to help elucidate the events and signaling processes, which happen during the courtship, pollen grain germination, and fertilization in olive. The use of subtractive SSH libraries constructed using, on the one hand one specific stage of the pistil development with germinating pollen grains, and on the other hand mature pollen grains may help to reveal the specific transcripts involved in the cited events. Such libraries have also been created by subtracting vegetative mRNAs (from leaves), in order to identify reproductive sequences only. A variety of transcripts have been identified in the mature pollen grains and in the pistil at the receptive stage. Among them, those related to defense, transport and oxidative metabolism are highlighted mainly in the pistil libraries where transcripts related to stress, and response to biotic and abiotic stimulus have a prominent position. Extensive lists containing information as regard to the specific transcripts determined for each stage and tissue are provided, as well as functional classifications of these gene products. Such lists were faced up to two recent datasets obtained in olive after transcriptomic and genomic approaches. The sequences and the differential expression level of the SSH-transcripts identified here, highly matched the transcriptomic information. Moreover, the unique presence of a representative number of these transcripts has been validated by means of qPCR approaches. The construction of SSH libraries using pistil and pollen, considering the high interaction between male-female counterparts, allowed the identification of transcripts with important roles in stigma physiology. The functions of many of the transcripts obtained are intimately related, and most of them are of pivotal importance in defense, pollen-stigma interaction and signaling.

Transcript Analysis and Regulative Events during Flower Development in Olive (Olea europaea L.)

The identification and characterization of transcripts involved in flower organ development, plant reproduction and metabolism represent key steps in plant phenotypic and physiological pathways, and may generate high-quality transcript variants useful for the development of functional markers. This study was aimed at obtaining an extensive characterization of the olive flower transcripts, by providing sound information on the candidate MADS-box genes related to the ABC model of flower development and on the putative genetic and molecular determinants of ovary abortion and pollen-pistil interaction. The overall sequence data, obtained by pyrosequencing of four cDNA libraries from flowers at different developmental stages of three olive varieties with distinct reproductive features (Leccino, Frantoio and Dolce Agogia), included approximately 465,000 ESTs, which gave rise to more than 14,600 contigs and approximately 92,000 singletons. As many as 56,700 unigenes were successfully annotated and provided gene ontology insights into the structural organization and putative molecular function of sequenced transcripts and deduced proteins in the context of their corresponding biological processes. Differentially expressed genes with potential regulatory roles in biosynthetic pathways and metabolic networks during flower development were identified. The gene expression studies allowed us to select the candidate genes that play well-known molecular functions in a number of biosynthetic pathways and specific biological processes that affect olive reproduction. A sound understanding of gene functions and regulatory networks that characterize the olive flower is provided.

The Quest for Molecular Regulation Underlying Unisexual Flower Development

The understanding of the molecular mechanisms responsible for the making of a unisexual flower has been a long-standing quest in plant biology. Plants with male and female flowers can be divided mainly into two categories: dioecious and monoecious, and both sexual systems co-exist in nature in ca of 10% of the angiosperms. The establishment of male and female traits has been extensively described in a hermaphroditic flower and requires the interplay of networks, directly and indirectly related to the floral organ identity genes including hormonal regulators, transcription factors, microRNAs, and chromatin-modifying proteins. Recent transcriptomic studies have been uncovering the molecular processes underlying the establishment of unisexual flowers and there are many parallelisms between monoecious, dioecious, and hermaphroditic individuals. Here, we review the paper entitled "Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber" published in 2014 in the Frontiers of Plant Science (volume 5 |Article 599) and discussed it in the context of recent studies with other dioecious and monoecious plants that utilized high-throughput platforms to obtain transcriptomic profiles of male and female unisexual flowers. In some unisexual flowers, the developmental programs that control organ initiation fail and male or female organs do not form, whereas in other species, organ initiation and development occur but they abort or arrest during different species-specific stages of differentiation. Therefore, a direct comparison of the pathways responsible for the establishment of unisexual flowers in different species are likely to reveal conserved modules of gene regulatory hubs involved in stamen or carpel development, as well as differences that reflect the different stages of development in which male and/or female organ arrest or loss-of-function occurs.

De novo transcriptome sequencing and gene expression profiling of spinach (Spinacia oleracea L.) leaves under heat stress

Spinach (Spinacia oleracea) has cold tolerant but heat sensitive characteristics. The spinach variety ‘Island,’ is suitable for summer periods. There is lack molecular information available for spinach in response to heat stress. In this study, high throughput de novo transcriptome sequencing and gene expression analyses were carried out at different spinach variety ‘Island’ leaves (grown at 24 °C (control), exposed to 35 °C for 30 min (S1), and 5 h (S2)). A total of 133,200,898 clean reads were assembled into 59,413 unigenes (average size 1259.55 bp). 33,573 unigenes could match to public databases. The DEG of controls vs S1 was 986, the DEG of control vs S2 was 1741 and the DEG of S1 vs S2 was 1587. Gene Ontology (GO) and pathway enrichment analysis indicated that a great deal of heat-responsive genes and other stress-responsive genes were identified in these DEGs, suggesting that the heat stress may have induced an extensive abiotic stress effect. Comparative transcriptome analysis found 896 unique genes in spinach heat response transcript. The expression patterns of 13 selected genes were verified by RT-qPCR (quantitative real-time PCR). Our study found a series of candidate genes and pathways that may be related to heat resistance in spinach.

Comparative transcriptome analysis reveals different molecular mechanisms of Bacillus coagulans 2-6 response to sodium lactate and calcium lactate during lactic acid production

Lactate production is enhanced by adding calcium carbonate or sodium hydroxide during fermentation. However, Bacillus coagulans 2-6 can produce more than 180 g/L L-lactic acid when calcium lactate is accumulated, but less than 120 g/L L-lactic acid when sodium lactate is formed. The molecular mechanisms by which B. coagulans responds to calcium lactate and sodium lactate remain unclear. In this study, comparative transcriptomic methods based on high-throughput RNA sequencing were applied to study gene expression changes in B. coagulans 2-6 cultured in non-stress, sodium lactate stress and calcium lactate stress conditions. Gene expression profiling identified 712 and 1213 significantly regulated genes in response to calcium lactate stress and sodium lactate stress, respectively. Gene ontology assignments of the differentially expressed genes were performed. KEGG pathway enrichment analysis revealed that ‘ATP-binding cassette transporters’ were significantly affected by calcium lactate stress, and ‘amino sugar and nucleotide sugar metabolism’ was significantly affected by sodium lactate stress. It was also found that lactate fermentation was less affected by calcium lactate stress than by sodium lactate stress. Sodium lactate stress had negative effect on the expression of ‘glycolysis/gluconeogenesis’ genes but positive effect on the expression of ‘citrate cycle (TCA cycle)’ genes. However, calcium lactate stress had positive influence on the expression of ‘glycolysis/gluconeogenesis’ genes and had minor influence on ‘citrate cycle (TCA cycle)’ genes. Thus, our findings offer new insights into the responses of B. coagulans to different lactate stresses. Notably, our RNA-seq dataset constitute a robust database for investigating the functions of genes induced by lactate stress in the future and identify potential targets for genetic engineering to further improve L-lactic acid production by B. coagulans.

Transcriptome profiling identifies differentially expressed genes in Huoyan goose ovaries between the laying period and ceased period

The Huoyan goose is famous for its high egg-laying performance and is listed as a nationally protected domestic animal by the Chinese government. To elucidate the key regulatory genes involved in Huoyan goose egg laying, RNA from ovarian tissue during the ceased and laying periods was sequenced using the Illumina HiSeq 2000 sequencing platform. More than 12 million reads were produced in ceased and laying libraries that included 11,896,423 and 12,534,799 clean reads, respectively. More than 20% of the reads were matched to the reference genome, and 23% of the reads were matched to reference genes. Genes with a false discovery rate (FDR) ≤0.001 and log₂ratio ≧1 or ≤−1 were characterized as differentially expressed, and 344 up-regulated and 344 down-regulated genes were classified into functional categories. Twelve genes that are mainly involved in pathways for reproduction regulation, such as steroid hormone biosynthesis, GnRH signaling pathways, oocyte meiosis, progesterone-mediated oocyte maturation, steroid biosynthesis, calcium signaling pathways, and G-protein coupled receptor signaling pathway were selected for validation by a quantitative real-time polymerase chain reaction (qRT-PCR) analysis, the qRT-PCR results are consistent with the general expression patterns of those genes from the Illumina sequencing. These data provide comprehensive gene expression information at the transcriptional level that might increase our understanding of the Huoyan goose's reproductive biology.

Analysis of the Citrullus colocynthis transcriptome during water deficit stress

Citrullus colocynthis is a very drought tolerant species, closely related to watermelon (C. lanatus var. lanatus), an economically important cucurbit crop. Drought is a threat to plant growth and development, and the discovery of drought inducible genes with various functions is of great importance. We used high throughput mRNA Illumina sequencing technology and bioinformatic strategies to analyze the C. colocynthis leaf transcriptome under drought treatment. Leaf samples at four different time points (0, 24, 36, or 48 hours of withholding water) were used for RNA extraction and Illumina sequencing. qRT-PCR of several drought responsive genes was performed to confirm the accuracy of RNA sequencing. Leaf transcriptome analysis provided the first glimpse of the drought responsive transcriptome of this unique cucurbit species. A total of 5038 full-length cDNAs were detected, with 2545 genes showing significant changes during drought stress. Principle component analysis indicated that drought was the major contributing factor regulating transcriptome changes. Up regulation of many transcription factors, stress signaling factors, detoxification genes, and genes involved in phytohormone signaling and citrulline metabolism occurred under the water deficit conditions. The C. colocynthis transcriptome data highlight the activation of a large set of drought related genes in this species, thus providing a valuable resource for future functional analysis of candidate genes in defense of drought stress.

A genome-wide perspective of miRNAome in response to high temperature, salinity and drought stresses in Brassica juncea (Czern) L

Micro RNAs (miRNAs) are involved in diverse biological processes including adaptive response towards abiotic stresses. To unravel small RNAs and more specifically miRNAs that can potentially regulate determinants of abiotic stress tolerance, next generation sequencing of B. juncea seedlings subjected to high temperature, high salt and drought conditions was carried out. With the help of UEA sRNA workbench software package, 51 conserved miRNAs belonging to 30 miRNA families were identified. As there was limited genomic information available for B. juncea, we generated and assembled its genome sequence at a low coverage. Using the generated sequence and other publically available Brassica genomic/transcriptomic resources as mapping reference, 126 novel (not reported in any plant species) were discovered for the first time in B. juncea. Further analysis also revealed existence of 32 and 37 star sequences for conserved and novel miRNAs, respectively. The expression of selected conserved and novel miRNAs under conditions of different abiotic stresses was revalidated through universal TaqMan based real time PCR. Putative targets of identified conserved and novel miRNAs were predicted in B. rapa to gain insights into functional roles manifested by B. juncea miRNAs. Furthermore, SPL2-like, ARF17-like and a NAC domain containing protein were experimentally validated as targets of miR156, miR160 and miR164 respectively. Investigation of gene ontologies linked with targets of known and novel miRNAs forecasted their involvement in various biological functions.

Transcriptome profiling of radish (Raphanus sativus L.) root and identification of genes involved in response to Lead (Pb) stress with next generation sequencing

Lead (Pb), one of the most toxic heavy metals, can be absorbed and accumulated by plant roots and then enter the food chain resulting in potential health risks for human beings. The radish (Raphanus sativus L.) is an important root vegetable crop with fleshy taproots as the edible parts. Little is known about the mechanism by which radishes respond to Pb stress at the molecular level. In this study, Next Generation Sequencing (NGS)-based RNA-seq technology was employed to characterize the de novo transcriptome of radish roots and identify differentially expressed genes (DEGs) during Pb stress. A total of 68,940 assembled unique transcripts including 33,337 unigenes were obtained from radish root cDNA samples. Based on the assembled de novo transcriptome, 4,614 DEGs were detected between the two libraries of untreated (CK) and Pb-treated (Pb1000) roots. Gene Ontology (GO) and pathway enrichment analysis revealed that upregulated DEGs under Pb stress are predominately involved in defense responses in cell walls and glutathione metabolism-related processes, while downregulated DEGs were mainly involved in carbohydrate metabolism-related pathways. The expression patterns of 22 selected genes were validated by quantitative real-time PCR, and the results were highly accordant with the Solexa analysis. Furthermore, many candidate genes, which were involved in defense and detoxification mechanisms including signaling protein kinases, transcription factors, metal transporters and chelate compound biosynthesis related enzymes, were successfully identified in response to heavy metal Pb. Identification of potential DEGs involved in responses to Pb stress significantly reflected alterations in major biological processes and metabolic pathways. The molecular basis of the response to Pb stress in radishes was comprehensively characterized. Useful information and new insights were provided for investigating the molecular regulation mechanism of heavy metal Pb accumulation and tolerance in root vegetable crops.