Analysis of floral transcription factors from Lycoris longituba

Transcriptomics analyses reveal the key genes involved in stamen petaloid formation in Alcea rosea L

Alcea rosea L. is a traditional flower with a long cultivation history. It is extensively cultivated in China and is widely planted in green belt parks or used as cut flowers and potted ornamental because of its rich colors and flower shapes. Double-petal A. rosea flowers have a higher aesthetic value compared to single-petal flowers, a phenomenon determined by stamen petaloid. However, the underlying molecular mechanism of this phenomenon is still very unclear. In this study, an RNA-based comparative transcriptomic analysis was performed between the normal petal and stamen petaloid petal of A. rosea. A total of 3,212 differential expressed genes (DEGs), including 2,620 up-regulated DEGs and 592 down-regulated DEGs, were identified from 206,188 unigenes. Numerous DEGs associated with stamen petaloid were identified through GO and KEGG enrichment analysis. Notably, there were 63 DEGs involved in the plant hormone synthesis and signal transduction, including auxin, cytokinin, gibberellin, abscisic acid, ethylene, brassinosteroid, jasmonic acid, and salicylic acid signaling pathway and 56 key transcription factors (TFs), such as MADS-box, bHLH, GRAS, and HSF. The identification of these DEGs provides an important clue for studying the regulation pathway and mechanism of stamen petaloid formation in A. rosea and provides valuable information for molecular plant breeding.

Comparative Transcriptome Analysis Identifies Key Regulatory Genes Involved in Anthocyanin Metabolism During Flower Development in Lycoris radiata

Lycoris is used as a garden flower due to the colorful and its special flowers. Floral coloration of Lycoris is a vital trait that is mainly regulated via the anthocyanin biosynthetic pathway. In this study, we performed a comparative transcriptome analysis of Lycoris radiata petals at four different flower development stages. A total of 38,798 differentially expressed genes (DEGs) were identified by RNA sequencing, and the correlation between the expression level of the DEGs and the anthocyanin content was explored. The identified DEGs are significantly categorized into 'flavonoid biosynthesis,' 'phenylpropanoid biosynthesis,' 'Tropane, piperidine and pyridine alkaloid biosynthesis,' 'terpenoid backbone biosynthesis' and 'plant hormone signal transduction' by Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The candidate genes involved in anthocyanin accumulation in L. radiata petals during flower development stages were also identified, which included 56 structural genes (especially LrDFR1 and LrFLS) as well as 27 key transcription factor DEGs (such as C3H, GATA, MYB, and NAC). In addition, a key structural gene namely LrDFR1 of anthocyanin biosynthesis pathway was identified as a hub gene in anthocyanin metabolism network. During flower development stages, the expression level of LrDFR1 was positively correlated with the anthocyanin content. Subcellular localization revealed that LrDFR1 is majorly localized in the nucleus, cytoplasm and cell membrane. Overexpression of LrDFR1 increased the anthocyanin accumulation in tobacco leaves and Lycoris petals, suggesting that LrDFR1 acts as a positively regulator of anthocyanin biosynthesis. Our results provide new insights for elucidating the function of anthocyanins in L. radiata petal coloring during flower development.

Assessing and integrating the transcriptome analysis with plant development, trichomes, and secondary metabolites yield potential in Mentha arvensisL

Mentha arvensisL. (corn mint) is well known for the production of menthol, a used commodity in flavouring industries, and provides natural fragrances. Glandular trichomes are responsible for producing specific secondary metabolites in vascular plants having species chemistry. Ten cultivars/varieties of M. arvensis, namely, Saksham, Kosi, Himalaya, Gomti, Sambhav, Kalka, Damroo, Kushal, and Shivalik, were used to study the developmental regulation of trichomes, essential oil yield, chemical constituents of essential oil and morphological parameters were estimated with gene expression using a randomized block design. Simultaneously, RNA sequence-based transcriptome analysis was done to reveal the transcription factors and differential gene analysis, which are responsible for the biosynthesis of essential oil as well as trichome development. Plant growth showed the maximum transition between 35 and 50 days stage, while essential oil and its metabolite bioconversion was observed in between 70 and 100 days stage. Glandular trichomes were maximally increased between 50, 70 days, and 100 days stage in var. Kosi followed by var. Saryu which has rapid growth in oil content. Menthol reductase activity was found to be a regulatory element during development, as it follows the inverse trend of menthol content and leads to menthol accumulation in subcuticular spaces. Transcriptional factors, cog, and nonredundant novel genes were identified. The composition of mintessential oils is regulated at multiple levels, including transcript abundance, catalytic properties of enzyme catalysts, and cell type-specific epigenetic processes.

Characterisation of the complete chloroplast genome of Lycoris longituba (Amaryllidaceae)

Lycoris longituba is a species in the Amaryllidaceae with high ornamental and medicinal value. It is also an endangered species in East China due to the narrow distribution. Here, we assembled and characterised the complete chloroplast (cp) genome of this species using high throughput sequencing and bioinformatics analysis. As a result, the whole cp genome size is 158,484 bp, including two inverted repeat (IR) regions of 26,765 bp, large single-copy region (LSC) of 86,458 bp and small single-copy region (SSC) of 18,496 bp. A total of 137 genes were identified, including 87 protein-coding genes, 42 tRNA genes, and eight rRNA genes. Phylogenomic analysis was carried out using complete cp genome of 18 species in five families, supporting the closer relationship between L. longituba and L. squamigera than L. radiata and L. sprengeri in Amaryllidaceae.

Genome-wide analysis of spatiotemporal gene expression patterns during floral organ development in Brassica rapa

Flowering is a key agronomic trait that directly influences crop yield and quality and serves as a model system for elucidating the molecular basis that controls successful reproduction, adaptation, and diversification of flowering plants. Adequate knowledge of continuous series of expression data from the floral transition to maturation is lacking in Brassica rapa. To unravel the genome expression associated with the development of early small floral buds (< 2 mm; FB2), early large floral buds (2-4 mm; FB4), stamens (STs) and carpels (CPs), transcriptome profiling was carried out with a Br300K oligo microarray. The results showed that at least 6848 known nonredundant genes (30% of the genes of the Br300K) were differentially expressed during the floral transition from vegetative tissues to maturation. Functional annotation of the differentially expressed genes (DEGs) (fold change ≥ 5) by comparison with a close relative, Arabidopsis thaliana, revealed 6552 unigenes (4579 upregulated; 1973 downregulated), including 131 Brassica-specific and 116 functionally known floral Arabidopsis homologs. Additionally, 1723, 236 and 232 DEGs were preferentially expressed in the tissues of STs, FB2, and CPs. These DEGs also included 43 transcription factors, mainly AP2/ERF-ERF, NAC, MADS-MIKC, C2H2, bHLH, and WRKY members. The differential gene expression during flower development induced dramatic changes in activities related to metabolic processes (23.7%), cellular (22.7%) processes, responses to the stimuli (7.5%) and reproduction (1%). A relatively large number of DEGs were observed in STs and were overrepresented by photosynthesis-related activities. Subsequent analysis via semiquantitative RT-PCR, histological analysis performed with in situ hybridization of BrLTP1 and transgenic reporter lines (BrLTP promoter::GUS) of B. rapa ssp. pekinensis supported the spatiotemporal expression patterns. Together, these results suggest that a temporally and spatially regulated process of the selective expression of distinct fractions of the same genome leads to the development of floral organs. Interestingly, most of the differentially expressed floral transcripts were located on chromosomes 3 and 9. This study generated a genome expression atlas of the early floral transition to maturation that represented the flowering regulatory elements of Brassica rapa.

Integrating Transcriptomic and GC-MS Metabolomic Analysis to Characterize Color and Aroma Formation during Tepal Development in Lycoris longituba

Lycoris longituba, belonging to the Amaryllidaceae family, is a perennial bulb bearing flowers with diverse colors and fragrance. Selection of cultivars with excellent colored and scented flowers has always been the breeding aim for ornamental plants. However, the molecular mechanisms underlying color fading and aroma production during flower expansion in L. longituba remain unclear. Therefore, to systematically investigate these important biological phenomena, the tepals of L. longituba from different developmental stages were used to screen and analyze the metabolic components and relevant genes. Utilizing the Illumina platform, a total of 144,922 unigenes were obtained from the RNA-Seq libraries. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that the phenylpropanoid biosynthesis and flavonoid biosynthesis pathways might play important roles during color and aroma changes. Metabolomic analysis identified 29 volatile organic components (VOCs) from different developmental stages of L. longituba tepals, and orthogonal partial least-squares discriminate analysis (OPLS-DA) revealed that trans-β-ocimene—a terpene—was the most important aroma compound. Meanwhile, we found the content of anthocyanin was significantly reduced during the tepal color fading process. Then, we identified two dihydroflavonol-4-reductase (DFR) and three terpene synthase (TPS) genes, for which expression changes coincided with the production patterns of anthocyanins and trans-β-ocimene, respectively. Furthermore, a number of MYB and bHLH transcription factors (TFs) which might be involved in color- and aroma-formation were also identified in L. longituba tepal transcriptomes. Taken together, this is the first comprehensive report of the color and fragrance in tepals of L. longituba and these results could be helpful in understanding these characteristics and their regulation networks.

Revealing the complex genetic structure of cultivated amaryllis (Hippeastrum hybridum) using transcriptome-derived microsatellite markers

Although amaryllis (Hippeastrum hybridum) plants are commonly used in physiological and ecological research, the extent of their genomic and genetic resources remains limited. The development of molecular markers is therefore of great importance to accelerate genetic improvements in Hippeastrum species. In this study, a total of 269 unique genes were defined that might regulate the flower spathe development of amaryllis. In addition, 2000 simple sequence repeats (SSRs) were detected based on 171,462 de novo assembled unigenes from transcriptome data, and 66,4091 single nucleotide polymorphisms (SNPs) were also detected as putative molecular markers. Twenty-one SSR markers were screened to evaluate the genetic diversity and population structure of 104 amaryllis accessions. A total of 98 SSR loci were amplified for all accessions. The results reveal that Nei’s gene diversity (H) values of these markers ranged between 0.055 and 0.394, whereas the average values of Shannon’s Information index (I) ranged between 0.172 and 0.567. Genetic tree analysis further demonstrates that all accessions can be grouped into three main clusters, which can be further divided into two subgroups. STRUCTURE-based analysis revealed that the highest ΔK values were observed when K = 5, K = 6, K = 7 and K = 8. The results of this study enable large-scale transcriptomics and classification of Hippeastrum genetic polymorphisms and will be useful in the future for resource conservation and production.

De novo sequence assembly and characterization of Lycoris aurea transcriptome using GS FLX titanium platform of 454 pyrosequencing

Background

Lycoris aurea, also called Golden Magic Lily, is an ornamentally and medicinally important species of the Amaryllidaceae family. To date, the sequencing of its whole genome is unavailable as a non-model organism. Transcriptomic information is also scarce for this species. In this study, we performed de novo transcriptome sequencing to produce the first comprehensive expressed sequence tag (EST) dataset for L. aurea using high-throughput sequencing technology.

Methodology and Principal Findings

Total RNA was isolated from leaves with sodium nitroprusside (SNP), salicylic acid (SA), or methyl jasmonate (MeJA) treatment, stems, and flowers at the bud, blooming, and wilting stages. Equal quantities of RNA from each tissue and stage were pooled to construct a cDNA library. Using 454 pyrosequencing technology, a total of 937,990 high quality reads (308.63 Mb) with an average read length of 329 bp were generated. Clustering and assembly of these reads produced a non-redundant set of 141,111 unique sequences, comprising 24,604 contigs and 116,507 singletons. All of the unique sequences were involved in the biological process, cellular component and molecular function categories by GO analysis. Potential genes and their functions were predicted by KEGG pathway mapping and COG analysis. Based on our sequence analysis and published literatures, many putative genes involved in Amaryllidaceae alkaloids synthesis, including PAL, TYDC OMT, NMT, P450, and other potentially important candidate genes, were identified for the first time in this Lycoris. Furthermore, 6,386 SSRs and 18,107 high-confidence SNPs were identified in this EST dataset.

Conclusions

The transcriptome provides an invaluable new data for a functional genomics resource and future biological research in L. aurea. The molecular markers identified in this study will provide a material basis for future genetic linkage and quantitative trait loci analyses, and will provide useful information for functional genomic research in future.

Comparison of the photosynthetic characteristics of four Lycoris species with leaf appearing in autumn under field conditions

The diurnal trends of gas exchange and chlorophyll fluorescence parameters in four Lycoris species (L. houdyshelii, L. aurea, L. radiata var. pumila and L. albiflora) were determined and compared with a portable photosynthesis analysis system. Our study revealed that L. houdyshelii had the lowest light compensation point (LCP), while the other three species had higher LCP (12.37-14.99 μmol m^-2 s^-1); L. aurea had the highest light saturation point (LSP) (1,189 μmol m^-2 s^-1), and L. houdyshelii and L. albiflora had lower LSP with the values being 322 and 345 μmol m^-2 s^-1, respectively, and L. radiata var. pumila showed the intermediate LSP. Both the species L. houdyshelii and L. albiflora exhibited a typical and obvious decline in net photosynthetic rate (P _N) during midday, which was not observed in L. aurea. This indicated a possible photoinhibition in L. houdyshelii and L. albiflora as the ratio of variable to maximum fluorescence (Fv/Fm) values were higher in these two species. The minimal fluorescence (F₀) values were lower in L. aurea and L. radiata var. pumila. The diurnal changes of transpiration rate (E) in all four species presented only one peak, appearing between 11:00 h or 13:00 h. By using simple correlation analyses, it was observed that the environmental factors affecting P _N were different among four species and the main factors were photosynthetic photon flux density (PPFD) and relative humidity especially for L. aurea and L. radiata. The results of studying indicated that the four species could be divided into two groups. The species L. radiata var. pumila and L. aurea were more adapted to a relatively high irradiance, and L. houdyshelii and L. albiflora could be grown in moderate-shade environment in order to scale up their growth and productivity.

Tissue-specific transcriptome profiling of Plutella xylostella third instar larval midgut

The larval midgut of diamondback moth, Plutella xylostella, is a dynamic tissue that interfaces with a diverse array of physiological and toxicological processes, including nutrient digestion and allocation, xenobiotic detoxification, innate and adaptive immune response, and pathogen defense. Despite its enormous agricultural importance, the genomic resources for P. xylostella are surprisingly scarce. In this study, a Bt resistant P. xylostella strain was subjected to the in-depth transcriptome analysis to identify genes and gene networks putatively involved in various physiological and toxicological processes in the P. xylostella larval midgut. Using Illumina deep sequencing, we obtained roughly 40 million reads containing approximately 3.6 gigabases of sequence data. De novo assembly generated 63,312 ESTs with an average read length of 416 bp, and approximately half of the P. xylostella sequences (45.4%, 28,768) showed similarity to the non-redundant database in GenBank with a cut-off E-value below 10(-5). Among them, 11,092 unigenes were assigned to one or multiple GO terms and 16,732 unigenes were assigned to 226 specific pathways. In-depth analysis identified genes putatively involved in insecticide resistance, nutrient digestion, and innate immune defense. Besides conventional detoxification enzymes and insecticide targets, novel genes, including 28 chymotrypsins and 53 ABC transporters, have been uncovered in the P. xylostella larval midgut transcriptome; which are potentially linked to the Bt toxicity and resistance. Furthermore, an unexpectedly high number of ESTs, including 46 serpins and 7 lysozymes, were predicted to be involved in the immune defense.As the first tissue-specific transcriptome analysis of P. xylostella, this study sheds light on the molecular understanding of insecticide resistance, especially Bt resistance in an agriculturally important insect pest, and lays the foundation for future functional genomics research. In addition, current sequencing effort greatly enriched the existing P. xylostella EST database, and makes RNAseq a viable option in the future genomic analysis.

Agave tequilana MADS genes show novel expression patterns in meristems, developing bulbils and floral organs

Agave tequilana is a monocarpic perennial species that flowers after 5-8 years of vegetative growth signaling the end of the plant's life cycle. When fertilization is unsuccessful, vegetative bulbils are induced on the umbels of the inflorescence near the bracteoles from newly formed meristems. Although the regulation of inflorescence and flower development has been described in detail for monocarpic annuals and polycarpic species, little is known at the molecular level for these processes in monocarpic perennials, and few studies have been carried out on bulbils. Histological samples revealed the early induction of umbel meristems soon after the initiation of the vegetative to inflorescence transition in A. tequilana. To identify candidate genes involved in the regulation of floral induction, a search for MADS-box transcription factor ESTs was conducted using an A. tequilana transcriptome database. Seven different MIKC MADS genes classified into 6 different types were identified based on previously characterized A. thaliana and O. sativa MADS genes and sequences from non-grass monocotyledons. Quantitative real-time PCR analysis of the seven candidate MADS genes in vegetative, inflorescence, bulbil and floral tissues uncovered novel patterns of expression for some of the genes in comparison with orthologous genes characterized in other species. In situ hybridization studies using two different genes showed expression in specific tissues of vegetative meristems and floral buds. Distinct MADS gene regulatory patterns in A. tequilana may be related to the specific reproductive strategies employed by this species.