De novo transcriptome analysis in radish (Raphanus sativus L.) and identification of critical genes involved in bolting and flowering

Development of SNP and InDel markers by genome resequencing and transcriptome sequencing in radish (Raphanus sativus L.)

BACKGROUND

Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) are the most abundant genetic variations and widely distribute across the genomes in plant. Development of SNP and InDel markers is a valuable tool for genetics and genomic research in radish (Raphanus sativus L.).

RESULTS

In this study, a total of 366,679 single nucleotide polymorphisms (SNPs) and 97,973 insertion-deletion (InDel) markers were identified based on genome resequencing between 'YZH' and 'XHT'. In all, 53,343 SNPs and 4,257 InDels were detected in two cultivars by transcriptome sequencing. Among the InDel variations, 85 genomic and 15 transcriptomic InDels were newly developed and validated PCR. The 100 polymorphic InDels markers generated 207 alleles among 200 Chinese radish germplasm, with an average 2.07 of the number of alleles (Na) and with an average 0.33 of the polymorphism information content (PIC). Population structure and phylogenetic relationship revealed that the radish cultivars from northern China were clustered together and the southwest China cultivars were clustered together. RNA-Seq analysis revealed that 11,003 differentially expressed genes (DEGs) were identified between the two cultivars, of which 5,020 were upregulated and 5,983 were downregulated. In total, 145 flowering time-related DGEs were detected, most of which were involved in flowering time integrator, circadian clock/photoperiod autonomous, and vernalization pathways. In flowering time-related DGEs region, 150 transcriptomic SNPs and 9 InDels were obtained.

CONCLUSIONS

The large amount of SNPs and InDels identified in this study will provide a valuable marker resource for radish genetic and genomic studies. The SNPs and InDels within flowering time-related DGEs provide fundamental insight into for dissecting molecular mechanism of bolting and flowering in radish.

Genetic Mechanisms for Hybrid Breeding in Vegetable Crops

To address the complex challenges faced by our planet such as rapidly changing climate patterns, food and nutritional insecurities, and the escalating world population, the development of hybrid vegetable crops is imperative. Vegetable hybrids could effectively mitigate the above-mentioned fundamental challenges in numerous countries. Utilizing genetic mechanisms to create hybrids not only reduces costs but also holds significant practical implications, particularly in streamlining hybrid seed production. These mechanisms encompass self-incompatibility (SI), male sterility, and gynoecism. The present comprehensive review is primarily focused on the elucidation of fundamental processes associated with floral characteristics, the genetic regulation of floral traits, pollen biology, and development. Specific attention is given to the mechanisms for masculinizing and feminizing cucurbits to facilitate hybrid seed production as well as the hybridization approaches used in the biofortification of vegetable crops. Furthermore, this review provides valuable insights into recent biotechnological advancements and their future utilization for developing the genetic systems of major vegetable crops.

Epistatic interactions among multiple copies of FLC genes with naturally occurring insertions correlate with flowering time variation in radish

Brassicaceae crops, which underwent whole-genome triplication during their evolution, have multiple copies of flowering-related genes. Interactions among multiple gene copies may be involved in flowering time regulation; however, this mechanism is poorly understood. In this study, we performed comprehensive, high-throughput RNA sequencing analysis to identify candidate genes involved in the extremely late-bolting (LB) trait in radish. Then, we examined the regulatory roles and interactions of radish FLOWERING LOCUS C (RsFLC) paralogs, the main flowering repressor candidates. Seven flowering integrator genes, five vernalization genes, nine photoperiodic/circadian clock genes and eight genes from other flowering pathways were differentially expressed in the early-bolting (EB) cultivar 'Aokubinagafuto' and LB radish cultivar 'Tokinashi' under different vernalization conditions. In the LB cultivar, RsFLC1 and RsFLC2 expression levels were maintained after 40 days of cold exposure. Bolting time was significantly correlated with the expression rates of RsFLC1 and RsFLC2. Using the EB × LB F2 population, we performed association analyses of genotypes with or without 1910- and 1627-bp insertions in the first introns of RsFLC1 and RsFLC2, respectively. The insertion alleles prevented the repression of their respective FLC genes under cold conditions. Interestingly, genotypes homozygous for RsFLC2 insertion alleles maintained high RsFLC1 and RsFLC3 expression levels under cold conditions, and two-way analysis of variance revealed that RsFLC1 and RsFLC3 expression was influenced by the RsFLC2 genotype. Our results indicate that insertions in the first introns of RsFLC1 and RsFLC2 contribute to the late-flowering trait in radish via different mechanisms. The RsFLC2 insertion allele conferred a strong delay in bolting by inhibiting the repression of all three RsFLC genes, suggesting that radish flowering time is determined by epistatic interactions among multiple FLC gene copies.

Transcriptome analysis of critical genes related to flowering in Mikania micrantha at different altitudes provides insights for a potential control

BACKGROUND

Mikania micrantha is a vine with strong invasion ability, and its strong sexual reproduction ability is not only the main factor of harm, but also a serious obstacle to control. M. micrantha spreads mainly through seed production. Therefore, inhibiting the flowering and seed production of M. micrantha is an effective strategy to prevent from continuing to spread.

RESULT

The flowering number of M. micrantha is different at different altitudes. A total of 67.01 Gb of clean data were obtained from nine cDNA libraries, and more than 83.47% of the clean reads were mapped to the reference genome. In total, 5878 and 7686 significantly differentially expressed genes (DEGs) were found in E2 vs. E9 and E13 vs. E9, respectively. Based on the background annotation and gene expression, some candidate genes related to the flowering pathway were initially screened, and their expression levels in the three different altitudes in flower bud differentiation showed the same trend. That is, at an altitude of 1300 m, the flower integration gene and flower meristem gene were downregulated (such as SOC1 and AP1), and the flowering inhibition gene was upregulated (such as FRI and SVP). Additionally, the results showed that there were many DEGs involved in the hormone signal transduction pathway in the flower bud differentiation of M. micrantha at different altitudes.

CONCLUSIONS

Our results provide abundant sequence resources for clarifying the underlying mechanisms of flower bud differentiation and mining the key factors inhibiting the flowering and seed production of M. micrantha to provide technical support for the discovery of an efficient control method.

Transcriptome Analysis of Pennisetum glaucum (L.) R. Br. Provides Insight Into Heat Stress Responses

Pennisetum glaucum (L.) R. Br., being widely grown in dry and hot weather, frequently encounters heat stress at various stages of growth. The crop, due to its inherent capacity, efficiently overcomes such stress during vegetative stages. However, the same is not always the case with the terminal (flowering through grain filling) stages of growth, where recovery from stress is more challenging. However, certain pearl millet genotypes such as 841-B are known to overcome heat stress even at the terminal growth stages. Therefore, we performed RNA sequencing of two contrasting genotypes of pearl millet (841-B and PPMI-69) subjected to heat stress (42°C for 6 h) at flowering stages. Over 274 million high quality reads with an average length of 150 nt were generated, which were assembled into 47,310 unigenes having an average length of 1,254 nucleotides, N50 length of 1853 nucleotides, and GC content of 53.11%. Blastx resulted in the annotation of 35,628 unigenes, and functional classification showed 15,950 unigenes designated to 51 Gene Ontology terms. A total of 13,786 unigenes were allocated to 23 Clusters of Orthologous Groups, and 4,255 unigenes were distributed to 132 functional Kyoto Encyclopedia of Genes and Genomes database pathways. A total of 12,976 simple sequence repeats and 305,759 SNPs were identified in the transcriptome data. Out of 2,301 differentially expressed genes, 10 potential candidate genes were selected based on log2 fold change and adjusted p value parameters for their differential gene expression by qRT-PCR. We were able to identify differentially expressed genes unique to either of the two genotypes, and also, some DEGs common to both the genotypes were enriched. The differential expression patterns suggested that 841-B 6 h has better ability to maintain homeostasis during heat stress as compared to PPMI-69 6 h. The sequencing data generated in this study, like the SSRs and SNPs, shall serve as an important resource for the development of genetic markers, and the differentially expressed heat responsive genes shall be used for the development of transgenic crops.

New insight into the pigment composition and molecular mechanism of flower coloration in tulip (Tulipa gesneriana L.) cultivars with various petal colors

Pigmentation of various components leads to different colors in tulip flowers. To understand the molecular basis of the petal coloration in tulip, integrative analyses of the pigment components and transcriptome profiles were conducted on four tulip cultivars with different petal colors. A total of four major anthocyanins and 46 carotenoids were identified. The anthocyanin cyanidin 3-O-galactoside showed markedly higher abundances in the B cultivar than in the other varieties, and among the 46 kinds of carotenoids, (E/Z)-phytoene, violaxanthin myristate and violaxanthin palmitate were the major components. The RNA-seq and qRT-PCR results indicated that the pigment accumulation was linked to the expression of genes involved in the anthocyanin and carotenoid biosynthesis pathways. Yeast two-hybrid (Y2H) assays showed the interaction between different regulator factors in tulip MYB-bHLH-WD40 (MBW) complexes. Co-expression analyses of genes were performed, which include anthocyanin and carotenoid biosynthesis genes and transcription factors involved in MYB, bHLH, WRKY, AUX-IAA and MADS-box. The co-expression network and related analysis provide a basis for the discovery of color regulatory factors. Taken together, our study sheds light on the anthocyanin and carotenoid synthesis pathways and candidate regulatory transcription factors underlying flower coloration and shows the potential of flower breeding or pigments engineering in tulips.

Bolting reduces ferulic acid and flavonoid biosynthesis and induces root lignification in Angelica sinensis

Angelica sinensis is a perennial herbaceous species that produces the bioactive metabolites ferulic acid and alkylphthalides widely applied in the treatment of cardio-cerebrovascular diseases. While the effects of bolting on plant biomass and metabolites accumulation have been partly investigated, the mechanism of bolting reducing metabolites biosynthesis is still limited. In this study, the root biomass, accumulations of ferulic acid, flavonoids and lignin, antioxidant capacity, and related genes expression at four different bolting stages were investigated. The results showed that there was a 2.2-, 2.4- and 2.9-fold decrease of the root biomass, ferulic acid and flavonoids contents, while a 2.9-fold increase of lignin content on a per plant basis during the bolting stages. The antioxidant capacity also exhibited significant decrease with growth and development. The differential expression levels of the 20 genes, which are involved in biosynthesis of ferulic acid (e.g. AsPAL1, As4CLs and AsHCT), flavonoids (e.g. AsCHS, AsCHI and AsI3'H) and lignin (e.g. AsCAD1 and AsLACs), were consistent with changes in the above metabolites accumulation. The findings will provide useful references for improving the production of bioactive metabolites in A. sinensis.

Identification of key gene networks controlling vernalization development characteristics of Isatis indigotica by full-length transcriptomes and gene expression profiles

Isatis indigotica Fort., as a common Chinese medicinal raw material, will lose its medicinal value if it blooms early, so it is highly valuable to clarify the induction mechanism of the vernalization of I. indigotica at low temperature. In this study, the concentrations of soluble sugar, proline, glutathione and zeatin in two germplasms of I. indigotica with different degrees of low temperature tolerance (Y1 and Y2) were determined at 10 days, 20 days and 30 days of low-temperature treatment, and the full-length transcriptome of 24 samples was sequenced by Nanopore sequencing with Oxford Nanopore Technologies (ONT). After that, the data of transcripts involved in the vernalization of I. indigotica at low temperature were obtained, and these transcripts were identified using weighted gene co-expression network analysis (WGCNA). The results revealed the massive accumulation of soluble sugar and proline in Y1 and Y2 after low temperature induction. A total of 18,385 new transcripts, 6168 transcription factors and 470 lncRNAs were obtained. Differential expression analysis showed that gibberellin, flavonoids, fatty acids and some processes related to low temperature response were significantly enriched. Eight key transcripts were identified by WGCNA, among which ONT.14640.1, ONT.9119.1, ONT.13080.2 and ONT.16007.1 encodes a flavonoid transporter, 9-cis-epoxycarotenoid dioxygenase 3 (NCED3), growth factor gene and L-aspartate oxidase in plants, respectively. It indicated that secondary metabolites such as hormones and flavonoids play an important role in the vernalization of I. indigotica. qRT-PCR proved the reliability of transcriptome results. These results provide important insights on the low-temperature vernalization of I. indigotica, and provide a research basis for analyzing the vernalization mechanism of I. indigotica.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01110-2.

Global identification of long non-coding RNAs involved in the induction of spinach flowering

Background

Spinach is a beneficial annual vegetable species and sensitive to the bolting or early flowering, which causes a large reduction in quality and productivity. Indeed, bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components. Although some key flowering responsive genes have been identified in spinach, non-coding RNA molecules like long non-coding RNAs (lncRNAs) were not investigated yet. Herein, we used bioinformatic approaches to analyze the transcriptome datasets from two different accessions Viroflay and Kashan at two vegetative and reproductive stages to reveal novel lncRNAs and the construction of the lncRNA-mRNA co-expression network. Additionally, correlations among gene expression modules and phenotypic traits were investigated; day to flowering was chosen as our interesting trait.

Results

In the present study, we identified a total of 1141 lncRNAs, of which 111 were differentially expressed between vegetative and reproductive stages. The GO and KEGG analyses carried out on the cis target gene of lncRNAs showed that the lncRNAs play an important role in the regulation of flowering spinach. Network analysis pinpointed several well-known flowering-related genes such as ELF, COL1, FLT, and FPF1 and also some putative TFs like MYB, WRKY, GATA, and MADS-box that are important regulators of flowering in spinach and could be potential targets for lncRNAs.

Conclusions

This study is the first report on identifying bolting and flowering-related lncRNAs based on transcriptome sequencing in spinach, which provides a useful resource for future functional genomics studies, genes expression researches, evaluating genes regulatory networks and molecular breeding programs in the regulation of the genetic mechanisms related to bolting in spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07989-1.

Association between RsFT, RsFLC and RsCOL5 (A&B) expression and flowering regulation in Japanese wild radish

Flowering is an important step in the life cycle of plants and indicates adaptability to external climatic cues such as temperature and photoperiod. We investigated the expression patterns of core genes related to flowering-time regulation in Japanese wild radish (Raphanus sativus var. raphanistroides) with different vernalization requirements (obligate and facultative) and further identified climatic cues that may act as natural selective forces. Specifically, we analysed flowering-time variation under different cold and photoperiod treatments in Japanese wild radish collected from the Hokkaido (northern lineage) and Okinawa (southern lineage) islands, which experience contrasting climatic cues. The cultivation experiment verified the obligate and facultative vernalization requirements of the northern and southern wild radish accessions, respectively. The expression of major genes involved in flowering time indicated that RsFLC and RsCOL5 (A&B) may interact to regulate flowering time. Notably, floral initiation in the northern lineage was strongly correlated with RsFLC expression, whereas flowering in the southern linage was correlated with induction of RsCOL5-A expression, despite high RsFLC transcript levels. These results suggested that the northern accessions are more sensitive to prolonged cold exposure, whereas the southern accessions are more sensitive to photoperiod. These different mechanisms ultimately confer an optimal flowering time in natural populations in response to locally contrasting climatic cues. This study provides new insights into the variant mechanisms underlying floral pathways in Japanese wild radish from different geographic locations.

Transcriptome architecture reveals genetic networks of bolting regulation in spinach

BACKGROUND

Bolting refers to the early flowering stem production on agricultural and horticultural crops before harvesting. Indeed, bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components, which cause a large reduction in the quality and productivity of vegetable crops like spinach. However, little is known about the signaling pathways and molecular functions involved in bolting mechanisms in spinach. The genetic information regarding the transition from vegetative growth to the reproductive stage in spinach would represent an advantage to regulate bolting time and improvement of resistant cultivars to minimize performance loss.

RESULTS

To investigate the key genes and their genetic networks controlling spinach bolting, we performed RNA-seq analysis on early bolting accession Kashan and late-bolting accession Viroflay at both vegetative and reproductive stages and found a significant number of differentially expressed genes (DEGs) ranging from 195 to 1230 in different comparisons. These genes were mainly associated with the signaling pathways of vernalization, photoperiod/circadian clock, gibberellin, autonomous, and aging pathways. Gene ontology analysis uncovered terms associated with carbohydrate metabolism, and detailed analysis of expression patterns for genes of Fructose-1, 6-bisphosphate aldolase, TREHALOSE-6-PHOSPHATE SYNTHASE 1, FLOWERING PROMOTING FACTOR 1, EARLY FLOWERING, GIGANTEA, and MADS-box proteins revealed their potential roles in the initiating or delaying of bolting.

CONCLUSION

This study is the first report on identifying bolting and flowering-related genes based on transcriptome sequencing in spinach, which provides insight into bolting control and can be useful for molecular breeding programs and further study in the regulation of the genetic mechanisms related to bolting in other vegetable crops.

Temperature-related cline in the root mass fraction in East Asian wild radish along the Japanese archipelago

Wild plants with a wide distribution, including those exposed to a wide variety of environmental conditions, may have variations in key functional traits relevant for agricultural applications. The East Asian wild radish (Raphanus sativus var. raphanistroides) is an appropriate model plant because it is widely distributed and has outstanding sink capacity as well as two cultivars within the species. Multiple common garden trials with 14 populations and three testing sites were conducted across the Japanese archipelago to quantify variations in yield and allocation. Significant inter-population variations and interaction effects with testing sites were detected for the root and shoot mass and the root mass fraction (RMF). While the rank order of the population changed drastically among sites and the variance components of genetic effects were small in yield traits (2.4%-4.7%), RMF displayed a large genetic variance (23.2%) and was consistently higher in the northern populations at all sites. Analyses revealed that the mean temperature of growing season of the seed origin was the most prominent factor explaining variation in RMF, irrespective of the sites. We concluded that the trait of resource allocation had a temperature-related cline and plants in cooler climates could invest more resources into their roots.

Gibberellin Promotes Bolting and Flowering via the Floral Integrators RsFT and RsSOC1-1 under Marginal Vernalization in Radish

Gibberellic acid (GA) is one of the factors that promotes flowering in radish (Raphanus Sativus L.), although the mechanism mediating GA activation of flowering has not been determined. To identify this mechanism in radish, we compared the effects of GA treatment on late-flowering (NH-JS1) and early-flowering (NH-JS2) radish lines. GA treatment promoted flowering in both lines, but not without vernalization. NH-JS2 plants displayed greater bolting and flowering pathway responses to GA treatment than NH-JS1. This variation was not due to differences in GA sensitivity in the two lines. We performed RNA-seq analysis to investigate GA-mediated changes in gene expression profiles in the two radish lines. We identified 313 upregulated, differentially expressed genes (DEGs) and 207 downregulated DEGs in NH-JS2 relative to NH-JS1 in response to GA. Of these, 21 and 8 genes were identified as flowering time and GA-responsive genes, respectively. The results of RNA-seq and quantitative PCR (qPCR) analyses indicated that RsFT and RsSOC1-1 expression levels increased after GA treatment in NH-JS2 plants but not in NH-JS1. These results identified the molecular mechanism underlying differences in the flowering-time genes of NH-JS1 and NH-JS2 after GA treatment under insufficient vernalization conditions.

Analysis of the quinoa genome reveals conservation and divergence of the flowering pathways

Quinoa (Chenopodium quinoa Willd.) is a grain crop grown in the Andes renowned as a highly nutritious plant exhibiting tolerance to abiotic stress such as drought, cold and high salinity. Quinoa grows across a range of latitudes corresponding to differing day lengths, suggesting regional adaptations of flowering regulation. Improved understanding and subsequent modification of the flowering process, including flowering time, ensuring high yields, is one of the key factors behind expansion of cultivation zones and goals of the crop improvement programs worldwide. However, our understanding of the molecular basis of flower initiation and development in quinoa is limited. Here, we use a computational approach to perform genome-wide identification and analysis of 611 orthologues of the Arabidopsis thaliana flowering genes. Conservation of the genes belonging to the photoperiod, gibberellin and autonomous pathways was observed, while orthologues of the key genes found in the vernalisation pathway (FRI, FLC) were absent from the quinoa genome. Our analysis indicated that on average each Arabidopsis flowering gene has two orthologous copies in quinoa. Several genes including orthologues of MIF1, FT and TSF were identified as homologue-rich genes in quinoa. We also identified 459 quinoa-specific genes uniquely expressed in the flower and/or meristem, with no known orthologues in other species. The genes identified provide a resource and framework for further studies of flowering in quinoa and related species. It will serve as valuable resource for plant biologists, crop physiologists and breeders to facilitate further research and establishment of modern breeding programs for quinoa.

Crosstalk in the darkness: bulb vernalization activates meristem transition via circadian rhythm and photoperiodic pathway

BACKGROUND

Geophytes possess specialized storage organs - bulbs, tubers, corms or rhizomes, which allow their survival during unfovarable periods and provide energy support for sprouting and sexual and vegetative reproduction. Bulbing and flowering of the geophyte depend on the combined effects of the internal and external factors, especially temperature and photoperiod. Many geophytes are extensively used in agriculture, but mechanisms of regulation of their flowering and bulbing are still unclear.

RESULTS

Comparative morpho-physiological and transcriptome analyses and quantitative validation of gene expression shed light on the molecular regulation of the responses to vernalization in garlic, a typical bulbous plant. Long dark cold exposure of bulbs is a major cue for flowering and bulbing, and its interactions with the genetic makeup of the individual plant dictate the phenotypic expression during growth stage. Photoperiod signal is not involved in the initial nuclear and metabolic processes, but might play role in the later stages of development, flower stem elongation and bulbing. Vernalization for 12 weeks at 4 °C and planting in November resulted in flower initiation under short photoperiod in December-January, and early blooming and bulbing. In contrast, non-vernalized plants did not undergo meristem transition. Comparisons between vernalized and non-vernalized bulbs revealed ~ 14,000 differentially expressed genes.

CONCLUSIONS

Low temperatures stimulate a large cascades of molecular mechanisms in garlic, and a variety of flowering pathways operate together for the benefit of meristem transition, annual life cycle and viable reproduction results.The circadian clock appears to play a central role in the transition of the meristem from vegetative to reproductive stage in bulbous plant, serving as integrator of the low-temperature signals and the expression of the genes associated with vernalization, photoperiod and meristem transition. The reserved photoperiodic pathway is integrated at an upstream point, possibly by the same receptors. Therefore, in bulb, low temperatures stimulate cascades of developmental mechanisms, and several genetic flowering pathways intermix to achieve successful sexual and vegetative reproduction.

De novo sequencing of tree peony (Paeonia suffruticosa) transcriptome to identify critical genes involved in flowering and floral organ development

BACKGROUND

Tree peony (Paeonia suffruticosa Andrews) is a globally famous ornamental flower, with large and colorful flowers and abundant flower types. However, a relatively short and uniform flowering period hinders the applications and production of ornamental tree peony. Unfortunately, the molecular mechanism of regulating flowering time and floral organ development in tree peony has yet to be elucidated. Because of the absence of genomic information, 454-based transcriptome sequence technology for de novo transcriptomics was used to identify the critical flowering genes using re-blooming, non-re-blooming, and wild species of tree peonies.

RESULTS

A total of 29,275 unigenes were obtained from the bud transcriptome, with an N50 of 776 bp. The average length of unigenes was 677.18 bp, and the longest sequence was 5815 bp. Functional annotation showed that 22,823, 17,321, 13,312, 20,041, and 9940 unigenes were annotated by NCBI-NR, Swiss-Prot, COG, GO, and KEGG, respectively. Within the differentially expressed genes (DEGs) 64 flowering-related genes were identified and some important flowering genes were also characterized by bioinformatics methods, reverse transcript polymerase chain reaction (RT-PCR), and rapid-amplification of cDNA ends (RACE). Then, the putative genetic network of flowering induction pathways and a floral organ development model were put forward, according to the comparisons of DEGs in any two samples and expression levels of the important flowering genes in differentiated buds, buds from different developmental stages, and with GA or vernalization treated. In tree peony, five pathways (long day, vernalization, autonomous, age, and gibberellin) regulated flowering, and the floral organ development followed an ABCE model. Moreover, it was also found that the genes PsAP1, PsCOL1, PsCRY1, PsCRY2, PsFT, PsLFY, PsLHY, PsGI, PsSOC1, and PsVIN3 probably regulated re-blooming of tree peony.

CONCLUSION

This study provides a comprehensive report on the flowering-related genes in tree peony for the first time and investigated the expression levels of the critical flowering related genes in buds of different cultivars, developmental stages, differentiated primordium, and flower parts. These results could provide valuable insights into the molecular mechanisms of flowering time regulation and floral organ development.

Non-vernalization Flowering and Seed Set of Cabbage Induced by Grafting Onto Radish Rootstocks

Cabbage (Brassica oleracea var. capitata) requires a long-term low-temperature exposure for floral induction, causing a delay in the breeding cycle. The objective of this study is to develop a method to induce flowering in cabbage without low-temperature treatment, using a grafting method. We conducted grafting experiments using two flower-induced Chinese kale cultivars (B. oleracea var. alboglabra) and seven radish cultivars/accessions as rootstocks and investigated the flowering response of grafted cabbage scions without low-temperature treatment. "Watanabe-seiko No.1" cabbage, when grafted onto the two Chinese kale cultivars, did not formed flower buds. Flowering was successfully induced in "Watanabe-seiko No.1" by grafting onto three out of the seven tested radish cultivars, and in "Kinkei No.201" and "Red cabbage" by grafting onto one tested radish cultivar. In "Watanabe-seiko No.1," the earliest flower bud appearance was observed at 29 days after grafting. Seeds were also obtained from the three cabbage cultivars that flowered by grafting. Gene expression analysis of "Watanabe-seiko No.1" cabbage scions which formed flower buds by grafting, revealed high expression of the homolog of the floral integrator, SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (BoSOC1), at the time of flower bud appearance. However, in the same leaf samples, we observed low expression of two homologs of florigen, FLOWERING LOCUS T (BoFT.C2 and BoFT.C6). In addition, two homologs of the floral repressor FLOWERING LOCUS C (BoFLC3 and BoFLC4), which are known to be down-regulated before flower bud differentiation in the vernalization pathway, were highly expressed, indicating that grafting onto radish induces cabbage flowering independently of the vernalization pathway. The expression level of the radish FT homolog (RsFT) in "Rat's tail-G2," which had highly induced flowering in the grafted cabbage scion, was higher than in the other radish cultivars. However, although "Rat's tail-CH" effectively induced flowering in the cabbage scion, the expression of RsFT was low in this cultivar. In this study, floral induction of non-vernalized cabbage cannot be explained by the expression levels of RsFT in rootstock plants, alone. The flowering of non-vernalized cabbage would be induced by transmissible agents from rootstocks and not by the expression of cabbage FT, BoFT, from the scion itself.

Transcriptome sequencing assisted discovery and computational analysis of novel SNPs associated with flowering in Raphanus sativus in-bred lines for marker-assisted backcross breeding

The sequencing of radish genome aids in the better understanding and tailoring of traits associated with economic importance. In order to accelerate the genomics assisted breeding and genetic selection, transcriptomes of 33 radish inbred lines with diverse traits were sequenced for the development of single nucleotide polymorphic (SNP) markers. The sequence reads ranged from 2,560,543,741 bp to 20,039,688,139 bp with the GC (%) of 47.80-49.34 and phred quality score (Q30) of 96.47-97.54%. A total of 4951 polymorphic SNPs were identified among the accessions after stringent filtering and 298 SNPs with efficient marker assisted backcross breeding (MAB) markers were generated from the polymorphic SNPs. Further, functional annotations of SNPs revealed the effects and importance of the SNPs identified in the flowering process. The SNPs were predominantly associated with the four major flowering related transcription factors such as MYB, MADS box (AG), AP2/EREB, and bHLH. In addition, SNPs in the vital flowering integrator gene (FT) and floral repressors (EMBRYONIC FLOWER 1, 2, and FRIGIDA) were identified among the radish inbred lines. Further, 50 SNPs were randomly selected from 298 SNPs and validated using Kompetitive Allele Specific PCR genotyping system (KASP) in 102 radish inbred lines. The homozygosity of the inbred lines varied from 56 to 96% and the phylogenetic analysis resulted in the clustering of inbred lines into three subgroups. Taken together, the SNP markers identified in the present study can be utilized for the discrimination, seed purity test, and adjusting parental combinations for breeding in radish.

Comprehensive transcriptome analysis reveals distinct regulatory programs during vernalization and floral bud development of orchardgrass (Dactylis glomerata L.)

BACKGROUND

Vernalization and the transition from vegetative to reproductive growth involve multiple pathways, vital for controlling floral organ formation and flowering time. However, little transcription information is available about the mechanisms behind environmental adaption and growth regulation. Here, we used high-throughput sequencing to analyze the comprehensive transcriptome of Dactylis glomerata L. during six different growth periods.

RESULTS

During vernalization, 4689 differentially expressed genes (DEGs) significantly increased in abundance, while 3841 decreased. Furthermore, 12,967 DEGs were identified during booting stage and flowering stage, including 7750 up-regulated and 5219 down-regulated DEGs. Pathway analysis indicated that transcripts related to circadian rhythm, photoperiod, photosynthesis, flavonoid biosynthesis, starch, and sucrose metabolism changed significantly at different stages. Coexpression and weighted correlation network analysis (WGCNA) analysis linked different stages to transcriptional changes and provided evidence of inner relation modules associated with signal transduction, stress responses, cell division, and hormonal transport.

CONCLUSIONS

We found enrichment in transcription factors (TFs) related to WRKY, NAC, AP2/EREBP, AUX/IAA, MADS-BOX, ABI3/VP1, bHLH, and the CCAAT family during vernalization and floral bud development. TFs expression patterns revealed intricate temporal variations, suggesting relatively separate regulatory programs of TF modules. Further study will unlock insights into the ability of the circadian rhythm and photoperiod to regulate vernalization and flowering time in perennial grass.

Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.)

Background

Alternative splicing (AS) plays important roles in gene expression and proteome diversity. Single nucleotide polymorphism (SNP) and insertion/deletion (InDel) are abundant polymorphisms and co-dominant inheritance markers, which have been widely used in germplasm identification, genetic mapping and marker-assisted selection in plants. So far, however, little information is available on utilization of AS events and development of SNP and InDel markers from transcriptome in radish.

Results

In this study, three radish transcriptome datasets were collected and aligned to the reference radish genome. A total of 56,530 AS events were identified from three radish genotypes with intron retention (IR) being the most frequent AS type, which accounted for 59.4% of the total expressed genes in radish. In all, 22,412 SNPs and 9436 InDels were identified with an average frequency of 1 SNP/17.9 kb and 1 InDel/42.5 kb, respectively. A total of 43,680 potential SSRs were identified in 31,604 assembled unigenes with a density of 1 SSR/2.5 kb. The ratio of SNPs with nonsynonymous/synonymous mutations was 1.05:1. Moreover, 35 SNPs and 200 InDels were randomly selected and validated by Sanger sequencing, 83.9% of the SNPs and 70% of the InDels exhibited polymorphism among these three genotypes. In addition, the 15 SNPs and 125 InDels were found to be unevenly distributed on 9 linkage groups. Furthermore, 40 informative InDel markers were successfully used for the genetic diversity analysis on 32 radish accessions.

Conclusions

These results would not only provide new insights into transcriptome complexity and AS regulation, but also furnish large amount of molecular marker resources for germplasm identification, genetic mapping and further genetic improvement of radish in breeding programs.

Electronic supplementary material

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

De Novo RNA Sequencing and Transcriptome Analysis of Monascus purpureus and Analysis of Key Genes Involved in Monacolin K Biosynthesis

Monascus purpureus is an important medicinal and edible microbial resource. To facilitate biological, biochemical, and molecular research on medicinal components of M. purpureus, we investigated the M. purpureus transcriptome by RNA sequencing (RNA-seq). An RNA-seq library was created using RNA extracted from a mixed sample of M. purpureus expressing different levels of monacolin K output. In total 29,713 unigenes were assembled from more than 60 million high-quality short reads. A BLAST search revealed hits for 21,331 unigenes in at least one of the protein or nucleotide databases used in this study. The 22,365 unigenes were categorized into 48 functional groups based on Gene Ontology classification. Owing to the economic and medicinal importance of M. purpureus, most studies on this organism have focused on the pharmacological activity of chemical components and the molecular function of genes involved in their biogenesis. In this study, we performed quantitative real-time PCR to detect the expression of genes related to monacolin K (mokA-mokI) at different phases (2, 5, 8, and 12 days) of M. purpureus M1 and M1-36. Our study found that mokF modulates monacolin K biogenesis in M. purpureus. Nine genes were suggested to be associated with the monacolin K biosynthesis. Studies on these genes could provide useful information on secondary metabolic processes in M. purpureus. These results indicate a detailed resource through genetic engineering of monacolin K biosynthesis in M. purpureus and related species.

Genome-Wide Characterization of the MADS-Box Gene Family in Radish (Raphanus sativus L.) and Assessment of Its Roles in Flowering and Floral Organogenesis

The MADS-box gene family is an important transcription factor (TF) family that is involved in various aspects of plant growth and development, especially flowering time and floral organogenesis. Although it has been reported in many plant species, the systematic identification and characterization of MADS-box TF family is still limited in radish (Raphanus sativus L.). In the present study, a comprehensive analysis of MADS-box genes was performed, and a total of 144 MADS-box family members were identified from the whole radish genome. Meanwhile, a detailed list of MADS-box genes from other 28 plant species was also investigated. Through the phylogenetic analysis between radish and Arabidopsis thaliana, all the RsMADS genes were classified into two groups including 68 type I (31 Mα, 12 Mβ and 25Mγ) and 76 type II (70 MIKC^C and 6 MIKC^∗). Among them, 41 (28.47%) RsMADS genes were located in nine linkage groups of radish from R1 to R9. Moreover, the homologous MADS-box gene pairs were identified among radish, A. thaliana, Chinese cabbage and rice. Additionally, the expression profiles of RsMADS genes were systematically investigated in different tissues and growth stages. Furthermore, quantitative real-time PCR analysis was employed to validate expression patterns of some crucial RsMADS genes. These results could provide a valuable resource to explore the potential functions of RsMADS genes in radish, and facilitate dissecting MADS-box gene-mediated molecular mechanisms underlying flowering and floral organogenesis in root vegetable crops.

Transcriptomic Analysis Identifies Differentially Expressed Genes (DEGs) Associated with Bolting and Flowering in Radish (Raphanus sativus L.)

The transition of vegetative growth to bolting and flowering is an important process in the life cycle of plants, which is determined by numerous genes forming an intricate network of bolting and flowering. However, no comprehensive identification and profiling of bolting and flowering-related genes have been carried out in radish. In this study, RNA-Seq technology was applied to analyze the differential gene expressions during the transition from vegetative stage to reproductive stage in radish. A total of 5922 differentially expressed genes (DEGs) including 779 up-regulated and 5143 down-regulated genes were isolated. Functional enrichment analysis suggested that some DEGs were involved in hormone signaling pathways and the transcriptional regulation of bolting and flowering. KEGG-based analysis identified 37 DEGs being involved in phytohormone signaling pathways. Moreover, 95 DEGs related to bolting and flowering were identified and integrated into various flowering pathways. Several critical genes including FT, CO, SOC1, FLC, and LFY were characterized and profiled by RT-qPCR analysis. Correlation analysis indicated that 24 miRNA-DEG pairs were involved in radish bolting and flowering. Finally, a miRNA-DEG-based schematic model of bolting and flowering regulatory network was proposed in radish. These outcomes provided significant insights into genetic control of radish bolting and flowering, and would facilitate unraveling molecular regulatory mechanism underlying bolting and flowering in root vegetable crops.

MADS-Box Genes and Gibberellins Regulate Bolting in Lettuce (Lactuca sativa L.)

Bolting in lettuce is promoted by high temperature and bolting resistance is of great economic importance for lettuce production. But how bolting is regulated at the molecular level remains elusive. Here, a bolting resistant line S24 and a bolting sensitive line S39 were selected for morphological, physiological, transcriptomic and proteomic comparisons. A total of 12204 genes were differentially expressed in S39 vs. S24. Line S39 was featured with larger leaves, higher levels of chlorophyll, soluble sugar, anthocyanin and auxin, consistent with its up-regulation of genes implicated in photosynthesis, oxidation-reduction and auxin actions. Proteomic analysis identified 30 differentially accumulated proteins in lines S39 and S24 upon heat treatment, and 19 out of the 30 genes showed differential expression in the RNA-Seq data. Exogenous gibberellins (GA) treatment promoted bolting in both S39 and S24, while 12 flowering promoting MADS-box genes were specifically induced in line S39, suggesting that although GA regulates bolting in lettuce, it may be the MADS-box genes, not GA, that plays a major role in differing the bolting resistance between these two lettuce lines.