De novo assembly of expressed transcripts and global analysis of the Phalaenopsis aphrodite transcriptome

An Outstandingly Rare Occurrence of Mycoviruses in Soil Strains of the Plant-Beneficial Fungi from the Genus Trichoderma and a Novel Polymycoviridae Isolate

In fungi, viral infections frequently remain cryptic causing little or no phenotypic changes. It can indicate either a long history of coevolution or a strong immune system of the host. Some fungi are outstandingly ubiquitous and can be recovered from a great diversity of habitats. However, the role of viral infection in the emergence of environmental opportunistic species is not known. The genus of filamentous and mycoparasitic fungi Trichoderma (Hypocreales, Ascomycota) consists of more than 400 species, which mainly occur on dead wood, other fungi, or as endo- and epiphytes. However, some species are environmental opportunists because they are cosmopolitan, can establish in a diversity of habitats, and can also become pests on mushroom farms and infect immunocompromised humans. In this study, we investigated the library of 163 Trichoderma strains isolated from grassland soils in Inner Mongolia, China, and found only four strains with signs of the mycoviral nucleic acids, including a strain of T. barbatum infected with a novel strain of the Polymycoviridae and named and characterized here as Trichoderma barbatum polymycovirus 1 (TbPMV1). Phylogenetic analysis suggested that TbPMV1 was evolutionarily distinct from the Polymycoviridae isolated either from Eurotialean fungi or from the order Magnaportales. Although the Polymycoviridae viruses were also known from Hypocrealean Beauveria bassiana, the phylogeny of TbPMV1 did not reflect the phylogeny of the host. Our analysis lays the groundwork for further in-depth characterization of TbPMV1 and the role of mycoviruses in the emergence of environmental opportunism in Trichoderma. IMPORTANCE Although viruses infect all organisms, our knowledge of some groups of eukaryotes remains limited. For instance, the diversity of viruses infecting fungi-mycoviruses-is largely unknown. However, the knowledge of viruses associated with industrially relevant and plant-beneficial fungi, such as Trichoderma spp. (Hypocreales, Ascomycota), may shed light on the stability of their phenotypes and the expression of beneficial traits. In this study, we screened the library of soilborne Trichoderma strains because these isolates may be developed into bioeffectors for plant protection and sustainable agriculture. Notably, the diversity of endophytic viruses in soil Trichoderma was outstandingly low. Only 2% of 163 strains contained traces of dsRNA viruses, including the new Trichoderma barbatum polymycovirus 1 (TbPMV1) characterized in this study. TbPMV1 is the first mycovirus found in Trichoderma. Our results indicate that the limited data prevent the in-depth study of the evolutionary relationship between soilborne fungi and is worth further investigation.

PbABCG1 and PbABCG2 transporters are required for the emission of floral monoterpenes in Phalaenopsis bellina

Terrestrial plants emit volatiles into the atmosphere to attract both pollinators and the enemies of herbivores, for defense. Phalaenopsis bellina is a scented orchid species in which the main scent components are monoterpenes, including linalool and geraniol, and their derivatives. Here, we investigated whether ABC transporters are involved in floral scent emission. We carried out whole-genome identification of ABC transporter-related genes using four floral transcriptomics libraries of P. bellina. We identified 86 ABC subfamily G genes related to terpenoid transport. After comparing the gene expression patterns of P. bellina with that of Phalaenopsis aphrodite subsp. formosana, a scentless species, followed by gene-to-gene correlation analysis, PbABCG1 and PbABCG2 were selected. The temporal expression of both PbABCG1 and PbABCG2 was highly correlated with that of the key enzyme PbGDPS and the major transcription factor PbbHLH4 in monoterpene biosynthesis, with optimal expression on day 5 post-anthesis. Spatial gene expression analysis showed that PbABCG1 was highly expressed in sepals, whereas PbABCG2 was expressed in the lip. Subcellular localization with a GFP fusion protein revealed that both PbABCG1 and PbABCG2 are cytoplasmic membrane proteins. Co-downregulation of PbABCG1 and PbABCG2 using both double-strand RNA interference and tobacco rattle virus-based gene silencing led to a significant decrease in monoterpene emission, accompanied by an increase in the internal monoterpene pools. Furthermore, ectopic expression of PbABCG1 and PbABCG2 in an ABC16^- mutant yeast strain rescued its tolerance to geraniol. Altogether, our results indicate that PbABCG1 and PbABCG2 play substantial roles in monoterpene transport/emission in P. bellina floral scent.

Transcriptome Analysis Reveals Endogenous Hormone Changes during Spike Development in Phalaenopsis

Phalaenopsis orchids are popular worldwide due to their high ornamental and economic value; the spike and inflorescence formation of their flowers could be efficiently controlled under proper conditions. In this study, transcriptomic profiles and endogenous hormone changes were investigated to better understand the spike formation of Phalaenopsis. Morphological observations revealed four spike initiation statuses (i.e., S0: the status refers to axillary buds remaining dormant in the leaf axils; S1: the status refers to the 0.5 cm-long initial spike; S2: the status refers to the 1 cm-long spike; S3: the status refers to the 3 cm-long spike) during the process of spike development, while anatomical observations revealed four related statuses of inflorescence primordium differentiation. A total of 4080 differentially expressed genes were identified based on pairwise comparisons of the transcriptomic data obtained from the S0 to S3 samples; high levels of differential gene expression were mostly observed in S1 vs. S2, followed by S0 vs. S1. Then, the contents of 12 endogenous hormones (e.g., irindole-3-acetic acid (IAA), salicylic acid (SA), abscisic acid (ABA), gibberellins, and cytokinins) were measured. The results showed that the ABA content was decreased from S0 to S1, while the gibberellic acid 1 (GA1) content exhibited an opposite trend, indicating the reduction in ABA levels combined with the increase in GA1 levels in S0 promoted the axillary bud dormancy breaking, preparing for the following spike initiation. The GA20 oxidase and ABA 8'-hydroxylase genes, which are involved in endogenous hormone metabolism and signaling pathways, displayed similar expression patterns, suggesting they were probably the key genes participating in the GA and ABA regulation. Taken together, the findings of this study indicate that GA and ABA may be the key endogenous hormones breaking the dormancy and promoting the germination of axillary buds in Phalaenopsis.

Identification and characterization of a new geminivirus from soybean plants and determination of V2 as a pathogenicity factor and silencing suppressor

BACKGROUND

Soybean is one of the four major crops in China. The occurrence of viruses in soybean causes significant economic losses.

RESULTS

In this study, the soybean leaves from stay-green plants showing crinkle were collected for metatranscriptomic sequencing. A novel geminivirus, tentatively named soybean geminivirus A (SGVA), was identified in soybean stay-green plants. Sequence analysis of the full-length SGVA genome revealed a genome of 2762 nucleotides that contain six open reading frames. Phylogenetic analyses revealed that SGVA was located adjacent to the clade of begomoviruses in both the full genome-based and C1-based phylogenetic tree, while in the CP-based phylogenetic tree, SGVA was located adjacent to the clade of becurtoviruses. SGVA was proposed as a new recombinant geminivirus. Agroinfectious clone of SGVA was constructed. Typical systemic symptoms of curly leaves were observed at 11 dpi in Nicotiana benthamiana plants and severe dwarfism was observed after 3 weeks post inoculation. Expression of the SGVA encoded V2 and C1 proteins through a potato virus X (PVX) vector caused severe symptoms in N. benthamiana. The V2 protein inhibited local RNA silencing in co-infiltration assays in GFP transgenic 16C N. benthamiana plants. Further study revealed mild symptoms in N. benthamiana plants inoculated with SGVA-ZZ _V2-STOP and SGVA-ZZ _V2-3738AA mutants. Both the relative viral DNA and CP protein accumulation levels significantly decreased when compared with SGVA-inoculated plants.

CONCLUSIONS

This work identified a new geminivirus in soybean stay-green plants and determined V2 as a pathogenicity factor and silencing suppressor.

Fusarivirus accessory helicases present an evolutionary link for viruses infecting plants and fungi

A significant number of mycoviruses have been identified that are related to plant viruses, but their evolutionary relationships are largely unexplored. A fusarivirus, Rhizoctonia solani fusarivirus 4 (RsFV4), was identified in phytopathogenic fungus Rhizoctonia solani (R. solani) strain XY74 co-infected by an alphaendornavirus. RsFV4 had a genome of 10,833 ​nt (excluding the poly-A tail), and consisted of four non-overlapping open reading frames (ORFs). ORF1 encodes an 825 aa protein containing a conserved helicase domain (Hel1). ORF3 encodes 1550 aa protein with two conserved domains, namely an RNA-dependent RNA polymerase (RdRp) and another helicase (Hel2). The ORF2 and ORF4 likely encode two hypothetical proteins (520 and 542 aa) with unknown functions. The phylogenetic analysis based on Hel2 and RdRp suggest that RsFV4 was positioned within the fusarivirus group, but formed an independent branch with three previously reported fusariviruses of R. solani. Notably, the Hel1 and its relatives were phylogenetically closer to helicases of potyviruses and hypoviruses than fusariviruses, suggesting fusarivirus Hel1 formed an evolutionary link between these three virus groups. This finding provides evidence of the occurrence of a horizontal gene transfer or recombination event between mycoviruses and plant viruses or between mycoviruses. Our findings are likely to enhance the understanding of virus evolution and diversity.

•

Rhizoctonia solani strain XY74 hosts two mycoviruses, fusarivirus (RsFV4) and endornavirus (RsAEV1).

•

RsFV4 consists of four ORFs and is evolutionarily associated to fusariviruses.

•

Two ORFs of RsFV4 encode two helicases belonging to superfamly II.

•

The accessory helicase of RsFV4 and its relatives are phylogenetically related to mycoviruses and plant viruses.

Pathogen resistance in Sphagneticola trilobata (Singapore daisy): molecular associations and differentially expressed genes in response to disease from a widespread fungus

Understanding the molecular associations underlying pathogen resistance in invasive plant species is likely to provide useful insights into the effective control of alien plants, thereby facilitating the conservation of native biodiversity. In the current study, we investigated pathogen resistance in an invasive clonal plant, Sphagneticola trilobata, at the molecular level. Sphagneticola trilobata (i.e., Singapore daisy) is a noxious weed that affects both terrestrial and aquatic ecosystems, and is less affected by pathogens in the wild than co-occurring native species. We used Illumina sequencing to investigate the transcriptome of S. trilobata following infection by a globally distributed generalist pathogen (Rhizoctonia solani). RNA was extracted from leaves of inoculated and un-inoculated control plants, and a draft transcriptome of S. trilobata was generated to examine the molecular response of this species following infection. We obtained a total of 49,961,014 (94.3%) clean reads for control (un-inoculated plants) and 54,182,844 (94.5%) for the infected treatment (inoculated with R. solani). Our analyses facilitated the discovery of 117,768 de novo assembled contigs and 78,916 unigenes. Of these, we identified 3506 differentially expressed genes and 60 hormones associated with pathogen resistance. Numerous genes, including candidate genes, were associated with plant-pathogen interactions and stress response in S. trilobata. Many recognitions, signaling, and defense genes were differentially regulated between treatments, which were confirmed by qRT-PCR. Overall, our findings improve our understanding of the genes and molecular associations involved in plant defense of a rapidly spreading invasive clonal weed, and serve as a valuable resource for further work on mechanism of disease resistance and managing invasive plants.

Physiological and transcriptomic analyses characterized high temperature stress response mechanisms in Sorbus pohuashanensis

Sorbus pohuashanensis (Hance) Hedl. is a Chinese native alpine tree species, but the problem of introducing S. pohuashanensis to low altitude areas has not been solved. In this study, we aimed to explore the molecular regulatory network of S. pohuashanensis in response to high-temperature stress using RNA-Sequencing technology and physiological and biochemical determination. Based on transcriptomic data, we obtained 1221 genes (752 up-regulated and 469 down-regulated) that were differentially expressed during 8 h 43℃ treatment and candidate genes were related to calcium signaling pathway, plant hormone signal transduction, heat shock factors, chaperones, ubiquitin mediated proteolysis, cell wall modification, ROS scavenging enzymes, detoxification and energy metabolism. The analysis of high temperature response at the physiological level and biochemical level were performed. The chlorophyll fluorescence parameters of leaf cells decreased, the content of osmotic regulators increased, and the activity of ROS scavenging enzymes decreased. The molecular regulatory network of S. pohuashanensis in response to high-temperature stress was preliminarily revealed in this study, which provides fundamental information improving introducing methods and discovering heat-tolerant genes involved in high-temperature stress in this species and provides a reference for other plants of the genus Sorbus.

Two distant helicases in one mycovirus: evidence of horizontal gene transfer between mycoviruses, coronaviruses and other nidoviruses

Nidovirales, which accommodates viruses with the largest RNA genomes, includes the notorious coronaviruses; however, the evolutionary route for nidoviruses is not well understood. We have characterized a positive-sense (+) single-stranded (ss) RNA mycovirus, Rhizoctonia solani hypovirus 2 (RsHV2), from the phytopathogenic fungus Rhizoctonia solani. RsHV2 has the largest RNA genome size of 22,219 nucleotides, excluding the poly(A) tail, in all known mycoviruses, and contains two open reading frames (ORF1 and ORF2). ORF1 encodes a protein of 2,009 amino acid (aa) that includes a conserved helicase domain belonging to helicase superfamily I (SFI). In contrast, ORF2 encodes a polyprotein of 4459 aa containing the hallmark genes of hypoviruses. The latter includes a helicase belonging to SFII. Following phylogenetic analysis, the ORF1-encoded helicase (Hel1) unexpectedly clustered in an independent evolutionary branch together with nidovirus helicases, including coronaviruses, and bacteria helicases. Thus, Hel1 presence indicates the occurrence of horizontal gene transfer between viruses and bacteria. These findings also suggest that RsHV2 is most likely a recombinant arising between hypoviruses and nidoviruses.

Identification of a Novel NtLRR-RLK and Biological Pathways That Contribute to Tolerance of TMV in Nicotiana tabacum

Tobacco mosaic virus (TMV) infection can causes serious damage to tobacco crops. To explore the approach of preventing TMV infection of plants, two tobacco cultivars with different resistances to TMV were used to analyze transcription profiling before and after TMV infection. The involvement of biological pathways differed between the tolerant variety (Yuyan8) and the susceptible variety (NC89). In particular, the plant-virus interaction pathway was rapidly activated in Yuyan8, and specific resistance genes were enriched. Liquid chromatography tandem mass spectrometry analysis detected large quantities of antiviral substances in the tolerant Yuyan8. A novel Nicotiana tabacum leucine-rich repeat receptor kinase (NtLRR-RLK) gene was identified as being methylated and this was verified using bisulfite sequencing. Transient expression of TMV-green fluorescent protein in pRNAi-NtLRR-RLK transgenic plants confirmed that NtLRR-RLK was important for susceptibility to TMV. The specific protein interaction map generated from our study revealed that levels of BIP1, E3 ubiquitin ligase, and LRR-RLK were significantly elevated, and all were represented at node positions in the protein interaction map. The same expression tendency of these proteins was also found in pRNAi-NtLRR-RLK transgenic plants at 24 h after TMV inoculation. These data suggested that specific genes in the infection process can activate the immune signal cascade through different resistance genes, and the integration of signal pathways could produce resistance to the virus. These results contribute to the overall understanding of the molecular basis of plant resistance to TMV and in the long term could identify new strategies for prevention and control virus infection.

Recent advances in biotechnological studies on wild grapevines as valuable resistance sources for smart viticulture

Cultivated grapevines, Vitis vinifera subsp. sativa, are thought to have been domesticated from wild populations of Vitis vinifera subsp. sylvestris in Central Asia. V. vinifera subsp. sativa is one of the most economically important fruit crops worldwide. Since cultivated grapevines are susceptible to multiple biotic and abiotic soil factors, they also need to be grafted on resistant rootstocks that are mostly developed though hybridization between American wild grapevine species (V. berlandieri, V. riparia, and V. rupestris). Therefore, wild grapevine species are essential genetic materials for viticulture to face biotic and abiotic stresses in both cultivar and rootstock parts. Actually, viticulture faces several environmental constraints that are further intensified by climate change. Recently, several reports on biotic and abiotic stresses-response in wild grapevines revealed accessions tolerant to different constraints. The emergence of advanced techniques such as omics technologies, marker-assisted selection (MAS), and functional analysis tools allowed a more detailed characterization of resistance mechanisms in these wild grapevines and suggest a number of species (V. rotundifolia, V. rupestris, V. riparia, V. berlandieri and V. amurensis) have untapped potential for new resistance traits including disease resistance loci and key tolerance genes. The present review reports on the importance of different biotechnological tools in exploring and examining wild grapevines tolerance mechanisms that can be employed to promote elite cultivated grapevines under climate change conditions.

Genome-wide analysis of small RNAs from Odontoglossum ringspot virus and Cymbidium mosaic virus synergistically infecting Phalaenopsis

Cymbidium mosaic virus (CymMV) and Odontoglossum ringspot virus (ORSV) are the two most prevalent viruses infecting orchids and causing economic losses worldwide. Mixed infection of CymMV and ORSV could induce intensified symptoms as early at 10 days post-inoculation in inoculated Phalaenopsis amabilis, where CymMV pathogenesis was unilaterally enhanced by ORSV. To reveal the antiviral RNA silencing activity in orchids, we characterized the viral small-interfering RNAs (vsiRNAs) from CymMV and ORSV singly or synergistically infecting P. amabilis. We also temporally classified the inoculated leaf-tip tissues and noninoculated adjacent tissues as late and early stages of infection, respectively. Regardless of early or late stage with single or double infection, CymMV and ORSV vsiRNAs were predominant in 21- and 22-nt sizes, with excess positive polarity and under-represented 5'-guanine. While CymMV vsiRNAs mainly derived from RNA-dependent RNA polymerase-coding regions, ORSV vsiRNAs encompassed the coat protein gene and 3'-untranslated region, with a specific hotspot residing in the 3'-terminal pseudoknot. With double infection, CymMV vsiRNAs increased more than 5-fold in number with increasing virus titres. Most vsiRNA features remained unchanged with double inoculation, but additional ORSV vsiRNA hotspot peaks were prominent. The potential vsiRNA-mediated regulation of the novel targets in double-infected tissues thereby provides a different view of CymMV and ORSV synergism. Hence, temporally profiled vsiRNAs from taxonomically distinct CymMV and ORSV illustrate active antiviral RNA silencing in their natural host, Phalaenopsis, during both early and late stages of infection. Our findings provide insights into offence-defence interactions among CymMV, ORSV and orchids.

Molecular Characterization and Expression Profile of PaCOL1, a CONSTANS-like Gene in Phalaenopsis Orchid

CONSTANS (CO) and CONSTANS-like (COL) genes play important roles in coalescing signals from photoperiod and temperature pathways. However, the mechanism of CO and COLs involved in regulating the developmental stage transition and photoperiod/temperature senescing remains unclear. In this study, we identified a COL ortholog gene from the Taiwan native orchid Phalaenopsis aphrodite. The Phalaenopsis aphrodite CONSTANS-like 1 (PaCOL1) belongs to the B-box protein family and functions in the nucleus and cytosol. Expression profile analysis of Phalaenopsis aphrodite revealed that PaCOL1 was significantly expressed in leaves, but its accumulation was repressed during environmental temperature shifts. We found a differential profile for PaCOL1 accumulation, with peak accumulation at late afternoon and at the middle of the night. Arabidopsis with PaCOL1 overexpression showed earlier flowering under short-day (SD) conditions (8 h/23 °C light and 16 h/23 °C dark) but similar flowering time under long-day (LD) conditions (16 h/23 °C light and 8 h/23 °C dark). Transcriptome sequencing revealed several genes upregulated in PaCOL1-overexpressing Arabidopsis plants that were previously involved in flowering regulation of the photoperiod pathway. Yeast two-hybrid (Y2H) analysis and bimolecular fluorescence complementation (BiFC) analysis revealed that PaCOL1 could interact with a crucial clock-associated regulator, AtCCA1, and a flowering repressor, AtFLC. Furthermore, expressing PaCOL1 in cca1.lhy partially reversed the mutant flowering time under photoperiod treatment, which confirms the role of PaCOL1 function in the rhythmic associated factors for modulating flowering.

High-quality RNA isolation from pigment-rich Dendrobium flowers

Isolation of high-quality RNA from Dendrobium flowers is challenging because of the high levels of pigment, polysaccharides, and polyphenols. In the present study, an efficient CTAB method for RNA extraction from the pigment-rich flowers of Dendrobium was optimised. The optimised method yielded high quantities of RNA (10.1-12.9 µg/g). Spectrophotometric values of A_260/280 in the range of 2.2 to 2.4 and A_260/230 values of 2.0 suggested that the isolated RNA was free of polyphenols, polysaccharides, and protein contaminants. RNA integrity numbers determined by microfluidics were in the range of 7.9-8.9 indicative of intact RNA. In the improved method, the addition of 3 M NaCl and 3% PVP-10 in the extraction buffer, followed by an incubation period of 45 min at 65 °C, eliminated most of the polysaccharides, polyphenolic compounds, and denatured protein. Extraction with phenol:chloroform:isoamyl alcohol (125:24:1) effectively removed pigments from the aqueous phase, while the precipitation of RNA with lithium chloride minimised the co-precipitation of protein, DNA, and polysaccharide and resulted in the extraction of high quality of RNA. The suitability of the RNA for downstream processing was confirmed via RT-PCR amplification of Chalcone synthase gene from cDNA prepared from RNA isolated from different developmental stages of the flower of a Dendrobium hybrid. The present method will be highly useful for the isolation of RNA from pigment, polyphenol, and polysaccharide-rich plant tissues.

Comparative Transcriptome and iTRAQ Proteome Analyses Reveal the Mechanisms of Diapause in Aphidius gifuensis Ashmead (Hymenoptera: Aphidiidae)

Aphidius gifuensis Ashmead (Hymenoptera: Aphidiidae) is a solitary endoparasitoid used in the biological control of various aphids. Diapause plays an important role in the successful production and deployment of A. gifuensis. Diapause can effectively extend the shelf life of biological control agents and solve several practical production problems like long production cycles, short retention periods, and discontinuities between supply and demand. In recent years, studies have been conducted on the environmental regulation and physiological and biochemical mechanisms of diapause in A. gifuensis. Nevertheless, the molecular mechanism of diapause in this species remains unclear. In this study, we compared the transcriptomes and proteomes of diapause and non-diapause A. gifuensis to identify the genes and proteins associated with this process. A total of 557 transcripts and 568 proteins were differentially expressed between the two groups. Among them, (1) genes involved in trehalose synthesis such as glycogen synthase, glycogen phosphorylase, and trehalose 6-phosphate synthase were upregulated in diapause at mRNA or protein level while glycolysis and gluconeogenesis-related genes were downregulated, suggesting that A. gifuensis stores trehalose as an energy resource and cryoprotectant; (2) the expression of immune-related genes like C-type lectins, hemocyanin, and phenoloxidase was increased, which helps to maintain immunity during diapause; (3) a chitin synthase and several cuticular protein genes were upregulated to harden the cuticle of diapausing A. gifuensis larval. These findings improve our understanding of A. gifuensis. diapause and provide the foundation for further pertinent studies.

Development of polymorphic microsatellite markers by using de novo transcriptome assembly of Calanthe masuca and C. sinica (Orchidaceae)

BACKGROUND

Calanthe masuca and C. sinica are two genetically closely related species in Orchidaceae. C. masuca is widely distributed in Asia, whereas C. sinica is restricted to Yunnan and Guangxi Provinces in southwest China. Both play important roles in horticulture and are under the pressure of population decline. Understanding their genetic background can greatly help us develop effective conservation strategies for these species. Simple sequence repeats (SSRs) are useful for genetic diversity analysis, presumably providing key information for the study and preservation of the wild populations of the two species we are interested in.

RESULTS

In this study, we performed RNA-seq analysis on the leaves of C. masuca and C. sinica, obtaining 40,916 and 71,618 unigenes for each species, respectively. In total, 2,019/3,865 primer pairs were successfully designed from 3,764/7,189 putative SSRs, among which 197 polymorphic SSRs were screened out according to orthologous gene pairs. After mononucleotide exclusion, a subset of 129 SSR primers were analysed, and 13 of them were found to have high polymorphism levels. Further analysis demonstrated that they were feasible and effective against C. masuca and C. sinica as well as transferable to another species in Calanthe. Molecular evolutionary analysis revealed functional pathways commonly enriched in unigenes with similar evolutionary rates in the two species, as well as pathways specific to each species, implicating species-specific adaptation. The divergence time between the two closely related species was tentatively determined to be 3.42 ± 1.86 Mya.

CONCLUSIONS

We completed and analysed the transcriptomes of C. masuca and C. sinica, assembling large numbers of unigenes and generating effective polymorphic SSR markers. This is the first report of the development of expressed sequence tag (EST)-SSR markers for Calanthe. In addition, our study could enable further genetic diversity analysis and functional and comparative genomic studies on Calanthe.

PbbHLH4 regulates floral monoterpene biosynthesis in Phalaenopsis orchids

Floral scent is an important factor in attracting pollinators and repelling florivores. In Phalaenopsis bellina (Orchidaceae), the major floral scent components are monoterpenoids. Previously, we determined that expression of GERANYL DIPHOSPHATE SYNTHASE (PbGDPS) is highly correlated with monoterpene biosynthesis in Phalaenosis orchids. Here, we found that both cis- and trans-regulation were present on the GDPS promoters, with trans-regulation playing a key role. To investigate the regulation of biosynthesis of floral scent, we compared the transcriptomic data of two Phalaenopsis orchids with contrasting scent phenotypes. Eight transcription factors (TFs) that exhibited sequential elevation in abundance through floral development in P. bellina were identified, and their transcript levels were higher in the scented orchid than the scentless one. Five of these TFs transactivated several structural genes involved in monoterpene biosynthesis including PbbHLH4, PbbHLH6, PbbZIP4, PbERF1, and PbNAC1. Ectopic transient expression of each of these TFs in scentless orchids resulted in stimulation of terpenoid biosynthesis. PbbHLH4 most profoundly induced monoterpene biosynthesis, with a 950-fold increase of monoterpenoid production in the scentless orchid. In conclusion, we determined that biosynthesis of orchid floral monoterpenes was sequentially regulated, with PbbHLH4 playing a crucial role for monoterpene biosynthesis.

Comparative transcriptomics provides insight into the molecular basis of species diversification of section Trigonopedia (Cypripedium) on the Qinghai-Tibetan Plateau

Deceptive pollination is key to the species richness of Orchidaceae. However, the genetic basis of species diversification is still under study. Section Trigonopedia is a monophyletic clade of genus Cypripedium distributed in the southwest of China. The species of this section are pollinated by different flies. Pollinator differentiation makes section Trigonopedia an ideal group for studying the genetic basis underlying species diversification. Here, we sequenced the transcriptomes of eight species of the genus Cypripedium, including six co-flowering species of section Trigonopedia and two species outside this section as an outgroup. We reconstructed the phylogeny of the section with the combined 1572 single-copy genes extracted from the eight species and produced a highly resolved tree of the section. Furthermore, we combined substitution rate estimation and differential expression analysis to identify candidate genes, including genes related to floral scent synthesis and environmental adaptation, involved in species differentiation. Field investigations showed that these species have adapted to different habitats. We propose that the species diversification in this section is initiated by floral scent differentiation, followed by habitat differentiation, finally leading to speciation. This study sheds novel light on the diversification of closely related orchid species in the Qinghai-Tibetan region.

Post genomics era for orchid research

Among 300,000 species in angiosperms, Orchidaceae containing 30,000 species is one of the largest families. Almost every habitats on earth have orchid plants successfully colonized, and it indicates that orchids are among the plants with significant ecological and evolutionary importance. So far, four orchid genomes have been sequenced, including Phalaenopsis equestris, Dendrobium catenatum, Dendrobium officinale, and Apostaceae shengen. Here, we review the current progress and the direction of orchid research in the post genomics era. These include the orchid genome evolution, genome mapping (genome-wide association analysis, genetic map, physical map), comparative genomics (especially receptor-like kinase and terpene synthase), secondary metabolomics, and genome editing.

Flower Development and Sex Determination between Male and Female Flowers in Vernicia fordii

Vernicia fordii is a monoecious and diclinous species with male and female flowers on the same inflorescence. Low female to male flower ratio is one of the main reasons for low yield in this species. However, little is known of its floral development and sex determination. Here, according to the results of scanning electron microscopy and histological analysis, the floral development of V. fordii was divided into 12 stages and the first morphological divergence between the male and female flowers was found to occur at stage 7. The male flowers are always unisexual, but the female flowers present bisexual characteristics, with sterile stamen (staminode) restricted to pre-meiosis of mother sporogenous cells and cell death occurring at later development stages. To further elucidate the molecular mechanism underling sex determination at the divergence stage for male and female flowers, comparative transcriptome analysis was performed. In total, 56,065 unigenes were generated and 608 genes were differentially expressed between male and female flowers, among which 310 and 298 DEGs (differentially expressed genes) showed high expression levels in males and females, respectively. The transcriptome data showed that the sexual dimorphism of female flowers was affected by jasmonic acid, transcription factors, and some genes related to the floral meristem activity. Ten candidate genes showed consistent expression in the qRT-PCR validation and DEGs data. In this study, we provide developmental characterization and transcriptomic information for better understanding of the development of unisexual flowers and the regulatory networks underlying the mechanism of sex determination in V. fordii, which would be helpful in the molecular breeding of V. fordii to improve the yield output.

Tissue-Specific Floral Transcriptome Analysis of the Sexually Deceptive Orchid Chiloglottis trapeziformis Provides Insights into the Biosynthesis and Regulation of Its Unique UV-B Dependent Floral Volatile, Chiloglottone 1

The Australian sexually deceptive orchid, Chiloglottis trapeziformis, employs a unique UV-B-dependent floral volatile, chiloglottone 1, for specific male wasp pollinator attraction. Chiloglottone 1 and related variants (2,5-dialkylcyclohexane-1,3-diones), represent a unique class of specialized metabolites presumed to be the product of cyclization between two fatty acid (FA) precursors. However, the genes involved in the biosynthesis of precursors, intermediates, and transcriptional regulation remains to be discovered. Chiloglottone 1 production occurs in the aggregation of calli (callus) on the labellum under continuous UV-B light. Therefore, deep sequencing, transcriptome assembly, and differential expression (DE) analysis were performed across different tissue types and UV-B treatments. Transcripts expressed in the callus and labellum (∼23,000 transcripts) were highly specialized and enriched for a diversity of known and novel metabolic pathways. DE analysis between chiloglottone-emitting callus versus the remainder of the labellum showed strong coordinated induction of entire FA biosynthesis and β-oxidation pathways including genes encoding Ketoacyl-ACP Synthase, Acyl-CoA Oxidase, and Multifunctional Protein. Phylogenetic analysis revealed potential gene duplicates with tissue-specific differential regulation including two Acyl-ACP Thioesterase B and a Ketoacyl-ACP Synthase genes. UV-B treatment induced the activation of UVR8-mediated signaling and large-scale transcriptome changes in both tissues, however, neither FA biosynthesis/β-oxidation nor other lipid metabolic pathways showed clear indications of concerted DE. Gene co-expression network analysis identified three callus-specific modules enriched with various lipid metabolism categories. These networks also highlight promising candidates involved in the cyclization of chiloglottone 1 intermediates (e.g., Bet v I and dimeric α,β barrel proteins) and orchestrating regulation of precursor pathways (e.g., AP2/ERF) given a strong co-regulation with FA biosynthesis/β-oxidation genes. Possible alternative biosynthetic routes for precursors (e.g., aldehyde dehydrogenases) were also indicated. Our comprehensive study constitutes the first step toward understanding the biosynthetic pathways involved in chiloglottone 1 production in Chiloglottis trapeziformis - supporting the roles of FA metabolism in planta, gene duplication as a potential source of new genes, and co-regulation of novel pathway genes in a tissue-specific manner. This study also provides a new and valuable resource for future discovery and comparative studies in plant specialized metabolism of other orchids and non-model plants.

De novo transcriptomic analysis and development of EST-SSRs for Sorbus pohuashanensis (Hance) Hedl

Sorbus pohuashanensis is a native tree species of northern China that is used for a variety of ecological purposes. The species is often grown as an ornamental landscape tree because of its beautiful form, silver flowers in early summer, attractive pinnate leaves in summer, and red leaves and fruits in autumn. However, development and further utilization of the species are hindered by the lack of comprehensive genetic information, which impedes research into its genetics and molecular biology. Recent advances in de novo transcriptome sequencing (RNA-seq) technology have provided an effective means to obtain genomic information from non-model species. Here, we applied RNA-seq for sequencing S. pohuashanensis leaves and obtained a total of 137,506 clean reads. After assembly, 96,213 unigenes with an average length of 770 bp were obtained. We found that 64.5% of the unigenes could be annotated using bioinformatics tools to analyze gene function and alignment with the NCBI database. Overall, 59,089 unigenes were annotated using the Nr database(non-redundant protein database), 35,225 unigenes were annotated using the GO (Gene Ontology categories) database, and 33,168 unigenes were annotated using COG (Cluster of Orthologous Groups). Analysis of the unigenes using the KEGG (Kyoto Encyclopedia of Genes and Genomes) database indicated that 13,953 unigenes were involved in 322 metabolic pathways. Finally, simple sequence repeat (SSR) site detection identified 6,604 unigenes that included EST-SSRs and a total of 7,473 EST-SSRs in the unigene sequences. Fifteen polymorphic SSRs were screened and found to be of use for future genetic research. These unigene sequences will provide important genetic resources for genetic improvement and investigation of biochemical processes in S. pohuashanensis.

Orchidstra 2.0-A Transcriptomics Resource for the Orchid Family

Orchidaceae, the orchid family, encompasses more than 25,000 species and five subfamilies. Due to their beautiful and exotic flowers, distinct biological and ecological features, orchids have aroused wide interest among both researchers and the general public. We constructed the Orchidstra database, a resource for orchid transcriptome assembly and gene annotations. The Orchistra database has been under active development since 2013. To accommodate the increasing amount of orchid transcriptome data and house more comprehensive information, Orchidstra 2.0 has been built with a new database system to store the annotations of 510,947 protein-coding genes and 161,826 noncoding transcripts, covering 18 orchid species belonging to 12 genera in five subfamilies of Orchidaceae. We have improved the N50 size of protein-coding genes, provided new functional annotations (including protein-coding gene annotations, protein domain/family information, pathways analysis, Gene Ontology term assignments, orthologous genes across orchid species, cross-links to the database of model species, and miRNA information), and improved the user interface with better website performance. We also provide new database functionalities for database searching and sequence retrieval. Moreover, the Orchidstra 2.0 database incorporates detailed RNA-Seq gene expression data from various tissues and developmental stages in different orchid species. The database will be useful for gene prediction and gene family studies, and for exploring gene expression in orchid species. The Orchidstra 2.0 database is freely accessible at http://orchidstra2.abrc.sinica.edu.tw.

Transcriptome Analysis to Identify the Putative Biosynthesis and Transport Genes Associated with the Medicinal Components of Achyranthes bidentata Bl

Achyranthes bidentata is a popular perennial medicine herb used for 1000s of years in China to treat various diseases. Although this herb has multiple pharmaceutical purposes in China, no transcriptomic information has been reported for this species. In addition, the understanding of several key pathways and enzymes involved in the biosynthesis of oleanolic acid and ecdysterone, two pharmacologically active classes of metabolites and major chemical constituents of A. bidentata root extracts, is limited. The aim of the present study was to characterize the transcriptome profile of the roots and leaves of A. bidentata to uncover the biosynthetic and transport mechanisms of the active components. In this study, we identified 100,987 transcripts, with an average length of 1146.8 base pairs. A total of 31,634 (31.33%) unigenes were annotated, and 12,762 unigenes were mapped to 303 pathways according to the Kyoto Encyclopedia of Genes and Genomes pathway database. Moreover, we identified a total of 260 oleanolic acid and ecdysterone genes encoding biosynthetic enzymes. Furthermore, the key enzymes involved in the oleanolic acid and ecdysterone synthesis pathways were analyzed using quantitative real-time polymerase chain reaction, revealing that the roots expressed these enzymes to a greater extent than the leaves. In addition, we identified 85 ATP-binding cassette transporters, some of which might be involved in the translocation of secondary metabolites.

Comparative transcriptional profile of the fish parasite Cryptocaryon irritans

BACKGROUND

Cryptocaryon irritans is an obligate ectoparasitic ciliate pathogen of marine fishes. It can infect most marine teleosts and cause heavy economic losses in aquaculture. There is currently no effective method of controlling this disease, and little information is available regarding the genes involved in its development and virulence. We aimed to investigate the distinct features of the three major life-cycle stages of C. irritans in terms of gene transcription level, and identify candidate vaccines/drug targets. We established a reference transcriptome of C. irritans by RNA-seq.

METHODS

Three cDNA libraries using total poly(A)⁺ mRNA isolated from trophonts, tomonts, and theronts was constructed and sequenced, respectively. Clean reads from the three stages were de novo assembled to generated unigene. Annotation of unigenes and transcriptomic comparison of three stages was performed.

RESULTS

Totals of 73.15, 62.23, and 109.57 million clean reads were generated from trophont, tomont, and theront libraries, respectively. After de novo assembly, 49,104 unigenes were obtained, including 9,253 unigenes with significant similarities to proteins from other ciliates. Transcriptomic comparisons revealed that 2,470 genes were differentially expressed among the three stages, including 2,011, 1,404, and 1,797 genes that were significantly differentially expressed in tomont/theront, tomont/trophont, and theront/trophont pairwise comparisons, respectively. Based on the results of hierarchical clustering, all differentially expressed genes (DEGs) were located in five major clusters. DEGs in clusters 1 and 2 were more highly expressed in tomonts than in other stages, DEGs in cluster 3 were dominant in the tomont and trophont stages, whereas clusters 4 and 5 included genes upregulated in the theront stage. In addition, Immobilization antigens (I-antigens) and proteases have long been considered major targets for vaccine development and potential drug targets in parasites, respectively. In the present study, nine putative I-antigens transcripts and 161 protease transcripts were found in the transcriptome of C. irritans.

CONCLUSION

It was concluded that DEGs enriched in tomonts were involved in cell division, to increase the number of theronts and ensure parasite continuity. DEGs enriched in theronts were associated with response to stimuli, whereas genes enriched in trophonts were related to nutrient accumulation and cell growth. In addition, the I-antigen and protease transcripts in our transcriptome could contribute to the development of vaccines or targeted drugs. Together, the results of the present study provide novel insights into the physiological processes of a marine parasitic ciliate.

De novo transcriptome assembly databases for the butterfly orchid Phalaenopsis equestris

Orchids are renowned for their spectacular flowers and ecological adaptations. After the sequencing of the genome of the tropical epiphytic orchid Phalaenopsis equestris, we combined Illumina HiSeq2000 for RNA-Seq and Trinity for de novo assembly to characterize the transcriptomes for 11 diverse P. equestris tissues representing the root, stem, leaf, flower buds, column, lip, petal, sepal and three developmental stages of seeds. Our aims were to contribute to a better understanding of the molecular mechanisms driving the analysed tissue characteristics and to enrich the available data for P. equestris. Here, we present three databases. The first dataset is the RNA-Seq raw reads, which can be used to execute new experiments with different analysis approaches. The other two datasets allow different types of searches for candidate homologues. The second dataset includes the sets of assembled unigenes and predicted coding sequences and proteins, enabling a sequence-based search. The third dataset consists of the annotation results of the aligned unigenes versus the Nonredundant (Nr) protein database, Kyoto Encyclopaedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG) databases with low e-values, enabling a name-based search.

De novo transcriptome sequencing and comparative analysis to discover genes related to floral development in Cymbidium faberi Rolfe

Cymbidium faberi is a traditional orchid flower in China that is highly appreciated for its fragrant aroma from its zygomorphic flowers. One bottleneck of the commercial production of C. faberi is the long vegetative growth phase of the orchid and the difficulty of the regulation of its flowering time. Moreover, its flower size, shape and color are often targeting traits for orchid breeders. Understanding the molecular mechanisms of floral development in C. faberi will ultimately benefit the genetic improvement of this orchid plant. The goal of this study is to identify potential genes and regulatory networks related to the floral development in C. faberi by using transcriptome sequencing, de novo assembly and computational analyses. The vegetative and flower buds of C. faberi were sampled for such comparisons. The RNA-seq yielded about 189,300 contigs that were assembled into 172,959 unigenes. Furthermore, a total of 13,484 differentially expressed unigenes (DEGs) were identified between the vegetative and flower buds. There were 7683 down-regulated and 5801 up-regulated DEGs in the flower buds compared to those in the vegetative buds, among which 3430 and 6556 DEGs were specifically enriched in the flower or vegetative buds, respectively. A total of 173 DEGs orthologous to known genes associated with the floral organ development, floral symmetry and flowering time were identified, including 12 TCP transcription factors, 34 MADS-box genes and 28 flowering time related genes. Furthermore, expression levels of ten genes potentially involved in floral development and flowering time were verified by quantitative real-time PCR. The identified DEGs will facilitate the functional genetic studies for further understanding the flower development of C. faberi.

Application of a modified drop method for high-resolution pachytene chromosome spreads in two Phalaenopsis species

BACKGROUND

Preparation of good chromosome spreads without cytoplasmic contamination is the crucial step in cytogenetic mapping. To date, cytogenetic research in the Orchidaceae family has been carried out solely on mitotic metaphase chromosomes. Well-spread meiotic pachytene chromosomes can provide higher resolution and fine detail for analysis of chromosomal structure and are also beneficial for chromosomal FISH (fluorescence in situ hybridization) mapping. However, an adequate method for the preparation of meiotic pachytene chromosomes in orchid species has not yet been reported.

RESULTS

Two Taiwanese native Phalaenopsis species were selected to test the modified drop method for preparation of meiotic pachytene chromosomes from pollinia. In this modified method, pollinia were ground and treated with an enzyme mixture to completely remove cell walls. Protoplasts were resuspended in ethanol/glacial acetic acid and dropped onto a wet inclined slide of 30° from a height of 0.5 m. The sample was then flowed down the inclined plane to spread the chromosomes. Hundreds of pachytene chromosomes with little to no cytoplasmic contamination were well spread on each slide. We also showed that the resolution of 45S rDNA-containing chromosomes at the pachytene stage was up to 20 times higher than that at metaphase. Slides prepared following this modified drop method were amenable to FISH mapping of both 45S and 5S rDNA on pachytene chromosomes and, after FISH, the chromosomal structure remained intact for further analysis.

CONCLUSION

This modified drop method is suitable for pachytene spreads from pollinia of Phalaenopsis orchids. The large number and high-resolution pachytene spreads, with little or no cytoplasmic contamination, prepared by the modified drop method could be used for FISH mapping of DNA fragments to accelerate the integration of cytogenetic and molecular research in Phalaenopsis orchids.

Transcriptome profiling of the microalga Chlorella pyrenoidosa in response to different carbon dioxide concentrations

To enrich our knowledge of carbon dioxide (CO₂)-concentrating mechanism (CCM) in eukaryotic algae, we used high-throughput sequencing to investigate the transcriptome profiling of the microalga Chlorella pyrenoidosa (Chlorophyta) response to different CO₂ levels. Altogether, 53.86 million (M) and 62.10M clean short reads of 100 nucleotides (nt) were generated from this microalga cultured at 4-fold air CO₂ (control) and air CO₂ concentrations by Illumina sequencing. A total of 32,662 unigenes were assembled from the two pooled samples. With an E-value cut-off of 1e-5, 9590, 6782, 5954, and 9092 unigenes were annotated in NR, Gene Ontology (GO), Eukaryotic Cluster of Orthologous Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. After screening, 51 differentially expressed unigenes were up-regulated and 8 were down-regulated in the air CO₂ group, relative to the control. The transcript levels of eight differentially expressed unigenes were validated by real-time quantitative PCR, which manifested that thioredoxin-like protein, laminin subunit beta-1, and chlorophyll a/b binding protein might be associated with the utilization of inorganic carbon at low CO₂ levels.

The genome and transcriptome of Phalaenopsis yield insights into floral organ development and flowering regulation

The Phalaenopsis orchid is an important potted flower of high economic value around the world. We report the 3.1 Gb draft genome assembly of an important winter flowering Phalaenopsis ‘KHM190’ cultivar. We generated 89.5 Gb RNA-seq and 113 million sRNA-seq reads to use these data to identify 41,153 protein-coding genes and 188 miRNA families. We also generated a draft genome for Phalaenopsis pulcherrima ‘B8802,’ a summer flowering species, via resequencing. Comparison of genome data between the two Phalaenopsis cultivars allowed the identification of 691,532 single-nucleotide polymorphisms. In this study, we reveal that the key role of PhAGL6b in the regulation of labellum organ development involves alternative splicing in the big lip mutant. Petal or sepal overexpressing PhAGL6b leads to the conversion into a lip-like structure. We also discovered that the gibberellin pathway that regulates the expression of flowering time genes during the reproductive phase change is induced by cool temperature. Our work thus depicted a valuable resource for the flowering control, flower architecture development, and breeding of the Phalaenopsis orchids.

Comparative transcriptome analysis of Gastrodia elata (Orchidaceae) in response to fungus symbiosis to identify gastrodin biosynthesis-related genes

BACKGROUND

Gastrodia elata Blume (Orchidaceae) is an important Chinese medicine with several functional components. In the life cycle of G. elata, the orchid develops a symbiotic relationship with two compatible mycorrhizal fungi Mycena spp. and Armillaria mellea during seed germination to form vegetative propagation corm and vegetative growth to develop tubers, respectively. Gastrodin (p-hydroxymethylphenol-beta-D-glucoside) is the most important functional component in G. elata, and gastrodin significantly increases from vegetative propagation corms to tubers. To address the gene regulation mechanism in gastrodin biosynthesis in G. elata, a comparative analysis of de novo transcriptome sequencing among the vegetative propagation corms and tubers of G. elata and A. mellea was conducted using deep sequencing.

RESULTS

Transcriptome comparison between the vegetative propagation corms and juvenile tubers of G. elata revealed 703 differentially expressed unigenes, of which 298 and 405 unigenes were, respectively up-regulated (fold-change ≥ 2, q-value < 0.05, the trimmed mean of M-values (TMM)-normalized fragments per kilobase of transcript per Million mapped reads (FPKM) > 10) and down-regulated (fold-change ≤ 0.5, q-value <0.05, TMM-normalized FPKM > 10) in juvenile tubers. After Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 112 up-regulated unigenes with KEGG Ortholog identifiers (KOids) or enzyme commission (EC) numbers were assigned to 159 isogroups involved in seventy-eight different pathways, and 132 down-regulated unigenes with KOids or EC numbers were assigned to 168 isogroups, involved in eighty different pathways. The analysis of the isogroup genes from all pathways revealed that the two unigenes TRINITY_DN54282_c0_g1 (putative monooxygenases) and TRINITY_DN50323_c0_g1 (putative glycosyltransferases) might participate in hydroxylation and glucosylation in the gastrodin biosynthetic pathway.

CONCLUSIONS

The gene expression of the two unique unigenes encoding monooxygenase and glycosyltransferase significantly increases from vegetative propagation corms to tubers, and the molecular basis of gastrodin biosynthesis in the tubers of G. elata is proposed.

Transcriptome-wide analysis of the MADS-box gene family in the orchid Erycina pusilla

Orchids exhibit a range of unique flower shapes and are a valuable ornamental crop. MADS-box transcription factors are key regulatory components in flower initiation and development. Changing the flower shape and flowering time can increase the value of the orchid in the ornamental horticulture industry. In this study, 28 MADS-box genes were identified from the transcriptome database of the model orchid Erycina pusilla. The full-length genomic sequences of these MADS-box genes were obtained from BAC clones. Of these, 27 were MIKC-type EpMADS (two truncated forms) and one was a type I EpMADS. Eleven EpMADS genes contained introns longer than 10 kb. Phylogenetic analysis classified the 24 MIKC(c) genes into nine subfamilies. Three specific protein motifs, AG, FUL and SVP, were identified and used to classify three subfamilies. The expression profile of each EpMADS gene correlated with its putative function. The phylogenetic analysis was highly correlated with the protein domain identification and gene expression results. Spatial expression of EpMADS6, EpMADS12 and EpMADS15 was strongly detected in the inflorescence meristem, floral bud and seed via in situ hybridization. The subcellular localization of the 28 EpMADS proteins was also investigated. Although EpMADS27 lacks a complete MADS-box domain, EpMADS27-YFP was localized in the nucleus. This characterization of the orchid MADS-box family genes provides useful information for both orchid breeding and studies of flowering and evolution.

De novo transcriptome sequencing and analysis of Coccinella septempunctata L. in non-diapause, diapause and diapause-terminated states to identify diapause-associated genes

BACKGROUND

The most common ladybird beetle, Coccinella septempunctata L., is an excellent predator of crop pests such as aphids and white flies, and it shows a wide range of adaptability, a large appetite and a high reproductive ability. Diapause research plays an important role in the artificial propagation and shelf-life extension of insect products. Although this lady beetle's regulatory, physiological and biochemical characteristics in the diapause period are well understood, the molecular mechanism of diapause remains unknown. Therefore, we collected female adults in three different states, i.e., non-diapause, diapause and diapause termination, for transcriptome sequencing.

RESULTS

After transcriptome sequencing using the Illumina HiSeq 2500 platform with pretreatment, a total of 417.6 million clean reads from nine samples were filtered using the program FASTX (version 0.0). Additionally, 106,262 contigs were assembled into 82,820 unigenes with an average length of 921 bp and an N50 of 1,241 bp. All of the unigenes were annotated through BLASTX alignment against the Nr or UniProt database, and 37,872 unigenes were matched. We performed further analysis of these unigenes using the Clusters of Orthologous Groups of proteins (COG), Gene Ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Through pairwise comparisons of the non-diapause (ND), diapause (D), and diapause-terminated (DT) groups, 3,501 and 1,427 differentially expressed genes (DEGs) were identified between D and ND and between DT and D, respectively. Moreover, 443 of the DEGs were specifically expressed during the diapause period (i.e., DEGs that were expressed at the highest or lowest levels during diapause compared with the other stages). GO function and KEGG pathway enrichment were performed on all DEGs and showed that RNA-directed DNA polymerase activity and fatty acid metabolism were significantly affected. Furthermore, eight specific expressed genes were selected for validation using qRT-PCR. Among these eight genes, seven genes were up-regulated, and one gene was down-regulated; the change trends of the eight genes were the same between the qRT-PCR and RNA-seq analysis results.

CONCLUSIONS

In this study, a new method for collecting and identifying diapause insects was described. We generated a vast quantity of transcriptome data from C. septempunctata L., providing a resource for gene function research. The diapause-associated genes that we identified establish a foundation for future studies on the molecular mechanisms of diapause.

Transcriptome Analysis of Syringa oblata Lindl. Inflorescence Identifies Genes Associated with Pigment Biosynthesis and Scent Metabolism

Syringa oblata Lindl. is a woody ornamental plant with high economic value and characteristics that include early flowering, multiple flower colors, and strong fragrance. Despite a long history of cultivation, the genetics and molecular biology of S. oblata are poorly understood. Transcriptome and expression profiling data are needed to identify genes and to better understand the biological mechanisms of floral pigments and scents in this species. Nine cDNA libraries were obtained from three replicates of three developmental stages: inflorescence with enlarged flower buds not protruded, inflorescence with corolla lobes not displayed, and inflorescence with flowers fully opened and emitting strong fragrance. Using the Illumina RNA-Seq technique, 319,425,972 clean reads were obtained and were assembled into 104,691 final unigenes (average length of 853 bp), 41.75% of which were annotated in the NCBI non-redundant protein database. Among the annotated unigenes, 36,967 were assigned to gene ontology categories and 19,956 were assigned to eukaryoticorthologous groups. Using the Kyoto Encyclopedia of Genes and Genomes pathway database, 12,388 unigenes were sorted into 286 pathways. Based on these transcriptomic data, we obtained a large number of candidate genes that were differentially expressed at different flower stages and that were related to floral pigment biosynthesis and fragrance metabolism. This comprehensive transcriptomic analysis provides fundamental information on the genes and pathways involved in flower secondary metabolism and development in S. oblata, providing a useful database for further research on S. oblata and other plants of genus Syringa.

Functional Characterization of PhapLEAFY, a FLORICAULA/LEAFY Ortholog in Phalaenopsis aphrodite

The plant-specific transcription factor LEAFY (LFY) is considered to be a master regulator of flower development in the model plant, Arabidopsis. This protein plays a dual role in plant growth, integrating signals from the floral inductive pathways and acting as a floral meristem identity gene by activating genes for floral organ development. Although LFY occupies an important position in flower development, the functional divergence of LFY homologs has been demonstrated in several plants including monocots and gymnosperms. In particular, the functional roles of LFY genes from orchid species such as Phalaenopsis that contain unique floral morphologies with distinct expression patterns of floral organ identity genes remain elusive. Here, PhapLFY, an ortholog of Arabidopsis LFY from Phalaenopsis aphrodite subsp. formosana, a Taiwanese native monopodial orchid, was isolated and characterized through analyses of expression and protein activity. PhapLFY transcripts accumulated in the floral primordia of developing inflorescences, and the PhapLFY protein had transcriptional autoactivation activity forming as a homodimer. Furthermore, PhapLFY rescues the aberrant floral phenotypes of Arabidopsis lfy mutants. Overexpression of PhapLFY alone or together with PhapFT1, a P. aphrodite subsp. formosana homolog of Arabidopsis FLOWERING LOCUS T (FT) in rice, caused precocious heading. Consistently, a higher Chl content in the sepals and morphological changes in epidermal cells were observed in the floral organs of PhapLFY knock-down orchids generated by virus-induced gene silencing. Taken together, these results suggest that PhapLFY is functionally distinct from RICE FLORICAULA/LEAFY (RFL) but similar to Arabidopsis LFY based on phenotypes of our transgenic Arabidopsis and rice plants.

Functional Characterization of Phalaenopsis aphrodite Flowering Genes PaFT1 and PaFD

We show that the key flowering regulators encoded by Phalaenopsis aphrodite FLOWERING LOCUS T1 (PaFT1) and PaFD share high sequence homologies to these from long-day flowering Arabidopsis and short-day flowering rice. Interestingly, PaFT1 is specifically up-regulated during flowering inductive cooling treatment but is not subjected to control by photoperiod in P. aphrodite. Phloem or shoot apex-specific expression of PaFT1 restores the late flowering of Arabidopsis ft mutants. Moreover, PaFT1 can suppress the delayed flowering caused by SHORT VEGATATIVE PHASE (SVP) overexpression as well as an active FRIGIDA (FRI) allele, indicating the functional conservation of flowering regulatory circuit in different plant species. PaFT1 promoter:GUS in Arabidopsis showed similar staining pattern to that of Arabidopsis FT in the leaves and guard cells but different in the shoot apex. A genomic clone or heat shock-inducible expression of PaFT1 is sufficient to the partial complementation of the ft mutants. Remarkably, ectopic PaFT1 expression also triggers precocious heading in rice. To further demonstrate the functional conservation of the flowering regulators, we show that PaFD, a bZIP transcription factor involved in flowering promotion, interacts with PaFT1, and PaFD partially complemented Arabidopsis fd mutants. Transgenic rice expressing PaFD also flowered early with increased expression of rice homologues of APETALA1 (AP1). Consistently, PaFT1 knock-down Phalaenopsis plants generated by virus-induced gene silencing exhibit delayed spiking. These studies suggest functional conservation of FT and FD genes, which may have evolved and integrated into distinct regulatory circuits in monopodial orchids, Arabidopsis and rice that promote flowering under their own inductive conditions.

dsRNA silencing of an R2R3-MYB transcription factor affects flower cell shape in a Dendrobium hybrid

BACKGROUND

The R2R3-MYB genes regulate pigmentation and morphogenesis of flowers, including flower and cell shape, and therefore have importance in the development of new varieties of orchids. However, new variety development is limited by the long breeding time required in orchids. In this study, we identified a cDNA, DhMYB1, that is expressed during flower development in a hybrid orchid, Dendrobium hybrida (Dendrobium bobby messina X Dendrobium chao phraya) then used the direct application of dsRNA to observe the effect of gene silencing on flower phenotype and floral epidermal cell shape.

RESULTS

Flower bud development in the Dendrobium hybrid was characterised into seven stages and the time of meiosis was determined as between stages 3 to 5 when the bud is approximately half of the mature size. Scanning electron microscopy characterisation of adaxial epidermal cells of the flower perianth, showed that the petals and sepals each are divided into two distinct domains based on cell shape and size, while the labellum comprises seven domains. Thirty-two partial cDNA fragments representing R2R3-MYB gene sequences were isolated from D. hybrida. Phylogenetic analysis revealed that nine of the translated sequences were clustered with MYB sequences that are known to be involved in cell shape development and from these, DhMYB1 was selected for full length cDNA cloning and functional study. Direct application of a 430 bp dsRNA from the 3' region of DhMYB1 to emerging orchid flower buds reduced expression of DhMYB1 RNA compared with untreated control. Scanning electron microscopy of adaxial epidermal cells within domain one of the labellum of flowers treated with DhMYB1 dsRNA showed flattened epidermal cells whilst those of control flowers were conical.

CONCLUSIONS

DhMYB1 is expressed throughout flower bud development and is involved in the development of the conical cell shape of the epidermal cells of the Dendrobium hybrida flower labellum. The direct application of dsRNA changed the phenotype of floral cells, thus, this technique may have application in floriculture biotechnology.

Application of an Integrated Omics Approach for Identifying Host Proteins That Interact With Odontoglossum ringspot virus Capsid Protein

The glutamic acid at position 100 (E(100)) in the capsid protein (CP) of Odontoglossum ringspot virus (ORSV) plays an important role in long-distance viral movement in Nicotiana benthamiana. The ORSV(E100A) mutant, which has a glutamic acid to alanine substitution, shows a loss of systemic infectivity in N. benthamiana. Transmission electron microscopy and size-exclusion chromatography assays showed that E(100) is essential for CP-CP interaction and viral particle assembly. To identify the ORSV triggering or response genes and CP-interacting proteins (CP-IP), an integrated omics approach based on next-generation sequencing and proteomics profiling was used in this study. The whole-transcriptomes of healthy and ORSV-infected leaves of N. benthamiana were analyzed, and the gene information was used to create a N. benthamiana protein database that was used for protein identification following mass spectrometry analysis. The integrated omics approach identified several putative host proteins that interact with ORSV CP(WT) and were categorized as photosystem subunits, defense-associated proteins, and cell division components. The expression pattern and CP interaction of these CP-IP were examined by semiquantitative reverse transcription polymerase chain reaction and an in vitro binding assay, respectively, to verify the in silico data. Among these proteins, a proteinase inhibitor of N. benthamiana (NbPI2) was highly associated with CP(E100A) as compared with CP(WT), and NbPI1 and NbPI2 were highly induced in ORSV-infected plants. NbPI1- and NbPI2-silenced plants (via a Tobacco rattle virus-induced gene-silencing system) did not exhibit a difference in ORSV infection. Thus, whether NbPI1 and NbPI2 play a role in plant immunity requires further investigation. In summary, the integrated omics approach provides massive and valuable information to identify the ORSV CP-IP and these CP-IP will help us to understand the movement of this virus and plant-virus interaction.

A de novo floral transcriptome reveals clues into Phalaenopsis orchid flower development

Phalaenopsis has a zygomorphic floral structure, including three outer tepals, two lateral inner tepals and a highly modified inner median tepal called labellum or lip; however, the regulation of its organ development remains unelucidated. We generated RNA-seq reads with the Illumina platform for floral organs of the Phalaenopsis wild-type and peloric mutant with a lip-like petal. A total of 43,552 contigs were obtained after de novo assembly. We used differentially expressed gene profiling to compare the transcriptional changes in floral organs for both the wild-type and peloric mutant. Pair-wise comparison of sepals, petals and labellum between peloric mutant and its wild-type revealed 1,838, 758 and 1,147 contigs, respectively, with significant differential expression. PhAGL6a (CUFF.17763), PhAGL6b (CUFF.17763.1), PhMADS1 (CUFF.36625.1), PhMADS4 (CUFF.25909) and PhMADS5 (CUFF.39479.1) were significantly upregulated in the lip-like petal of the peloric mutant. We used real-time PCR analysis of lip-like petals, lip-like sepals and the big lip of peloric mutants to confirm the five genes' expression patterns. PhAGL6a, PhAGL6b and PhMADS4 were strongly expressed in the labellum and significantly upregulated in lip-like petals and lip-like sepals of peloric-mutant flowers. In addition, PhAGL6b was significantly downregulated in the labellum of the big lip mutant, with no change in expression of PhAGL6a. We provide a comprehensive transcript profile and functional analysis of Phalaenopsis floral organs. PhAGL6a PhAGL6b, and PhMADS4 might play crucial roles in the development of the labellum in Phalaenopsis. Our study provides new insights into how the orchid labellum differs and why the petal or sepal converts to a labellum in Phalaenopsis floral mutants.

Genome-wide transcriptome and expression profile analysis of Phalaenopsis during explant browning

Background

Explant browning presents a major problem for in vitro culture, and can lead to the death of the explant and failure of regeneration. Considerable work has examined the physiological mechanisms underlying Phalaenopsis leaf explant browning, but the molecular mechanisms of browning remain elusive. In this study, we used whole genome RNA sequencing to examine Phalaenopsis leaf explant browning at genome-wide level.

Methodology/Principal Findings

We first used Illumina high-throughput technology to sequence the transcriptome of Phalaenopsis and then performed de novo transcriptome assembly. We assembled 79,434,350 clean reads into 31,708 isogenes and generated 26,565 annotated unigenes. We assigned Gene Ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations, and potential Pfam domains to each transcript. Using the transcriptome data as a reference, we next analyzed the differential gene expression of explants cultured for 0, 3, and 6 d, respectively. We then identified differentially expressed genes (DEGs) before and after Phalaenopsis explant browning. We also performed GO, KEGG functional enrichment and Pfam analysis of all DEGs. Finally, we selected 11 genes for quantitative real-time PCR (qPCR) analysis to confirm the expression profile analysis.

Conclusions/Significance

Here, we report the first comprehensive analysis of transcriptome and expression profiles during Phalaenopsis explant browning. Our results suggest that Phalaenopsis explant browning may be due in part to gene expression changes that affect the secondary metabolism, such as: phenylpropanoid pathway and flavonoid biosynthesis. Genes involved in photosynthesis and ATPase activity have been found to be changed at transcription level; these changes may perturb energy metabolism and thus lead to the decay of plant cells and tissues. This study provides comprehensive gene expression data for Phalaenopsis browning. Our data constitute an important resource for further functional studies to prevent explant browning.

De novo transcriptome assembly of pummelo and molecular marker development

Pummelo (Citrus grandis) is an important fruit crop worldwide because of its nutritional value. To accelerate the pummelo breeding program, it is essential to obtain extensive genetic information and develop relative molecular markers. Here, we obtained a 12-Gb transcriptome dataset of pummelo through a mixture of RNA from seven tissues using Illumina pair-end sequencing, assembled into 57,212 unigenes with an average length of 1010 bp. The annotation and classification results showed that a total of 39,584 unigenes had similar hits to the known proteins of four public databases, and 31,501 were classified into 55 Gene Ontology (GO) functional sub-categories. The search for putative molecular markers among 57,212 unigenes identified 10,276 simple sequence repeats (SSRs) and 64,720 single nucleotide polymorphisms (SNPs). High-quality primers of 1174 SSR loci were designed, of which 88.16% were localized to nine chromosomes of sweet orange. Of 100 SSR primers that were randomly selected for testing, 87 successfully amplified clear banding patterns. Of these primers, 29 with a mean PIC (polymorphic information content) value of 0.52 were effectively applied for phylogenetic analysis. Of the 20 SNP primers, 14 primers, including 54 potential SNPs, yielded target amplifications, and 46 loci were verified via Sanger sequencing. This new dataset will be a valuable resource for molecular biology studies of pummelo and provides reliable information regarding SNP and SSR marker development, thus expediting the breeding program of pummelo.

RNA-Seq-based analysis of changes in Borrelia burgdorferi gene expression linked to pathogenicity

Background

Lyme disease is a global public health problem caused by the spirochaete Borrelia burgdorferi. Our previous studies found differences in disease severity between B. burgdorferi B31- and B. garinii SZ-infected mice. We hypothesized that genes that are differentially expressed between Borrelia isolates encode bacterial factors that contribute to disease diversity.

Methods

The present study used high-throughput sequencing technology to characterize and compare the transcriptional profiles of B. burgdorferi B31 and B. garinii SZ cultured in vitro. Real-time quantitative RT-PCR was used to validate selected data from RNA-seq experiments.

Results

A total of 731 genes were differentially expressed between B. burgdorferi B31 and B. garinii SZ isolates, including those encoding lipoproteins and purine transport proteins. The fold difference in expression for B. garinii SZ versus B. burgdorferi B31 ranged from 22.07 to 1.01. Expression of the OspA, OspB and DbpB genes were significantly lower in B. garinii SZ compared to B. burgdorferi B31.

Conclusions

The results support the hypothesis that global changes in gene expression underlie differences in Borrelia pathogenicity. The findings also provide an empirical basis for studying the mechanism of action of specific genes as well as their potential usefulness for the diagnosis and management of Lyme disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-014-0623-2) contains supplementary material, which is available to authorized users.

Integrated transcriptome catalogue and organ-specific profiling of gene expression in fertile garlic (Allium sativum L.)

Background

Garlic is cultivated and consumed worldwide as a popular condiment and green vegetable with medicinal and neutraceutical properties. Garlic cultivars do not produce seeds, and therefore, this plant has not been the subject of either classical breeding or genetic studies. However, recent achievements in fertility restoration in a number of genotypes have led to flowering and seed production, thus enabling genetic studies and breeding in garlic.

Results

A transcriptome catalogue of fertile garlic was produced from multiplexed gene libraries, using RNA collected from various plant organs, including inflorescences and flowers. Over 32 million 250-bp paired-end reads were assembled into an extensive transcriptome of 240,000 contigs. An abundant transcriptome assembled separately from 102,000 highly expressed contigs was annotated and analyzed for gene ontology and metabolic pathways. Organ-specific analysis showed significant variation of gene expression between plant organs, with the highest number of specific reads in inflorescences and flowers. Analysis of the enriched biological processes and molecular functions revealed characteristic patterns for stress response, flower development and photosynthetic activity. Orthologues of key flowering genes were differentially expressed, not only in reproductive tissues, but also in leaves and bulbs, suggesting their role in flower-signal transduction and the bulbing process. More than 100 variants and isoforms of enzymes involved in organosulfur metabolism were differentially expressed and had organ-specific patterns. In addition to plant genes, viral RNA of at least four garlic viruses was detected, mostly in the roots and cloves, whereas only 1–4% of the reads were found in the foliage leaves.

Conclusions

The de novo transcriptome of fertile garlic represents a new resource for research and breeding of this important crop, as well as for the development of effective molecular markers for useful traits, including fertility and seed production, resistance to pests and neutraceutical characteristics.

Electronic supplementary material

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

Phenotype and transcriptome analysis reveals chloroplast development and pigment biosynthesis together influenced the leaf color formation in mutants of Anthurium andraeanum 'Sonate'

Leaf color is one of the well-sought traits in breeding program for Anthurium andraeanum Lind. Knowledge of mechanisms in anthuriums to produce leaves with different shades of green would help to effectively select desirable traits. In this study, the micro- and ultra-structural and physiological features of leaves on wild type and leaf color mutants (dark green, rubescent, etiolated, albino) in A. andraeanum 'Sonate' were analyzed. Results show that chloroplasts of leaf color mutants exhibited abnormal morphology and distribution. Using next generation sequencing technology followed by de novo assembly, leaf transcriptomes comprising of 41,017 unigenes with an average sequence length of 768 bp were produced from wild type and rubescent mutant. From the 27,539 (67.1%) unigenes with annotated functions, 858 significantly differently expressed genes (DEGs) were identified, consisting of 446 up-regulated genes and 412 down-regulated genes. Genes that affect chloroplasts development and division, and chlorophyll biosynthesis were included in the down-regulated DEGs. Quantitative real-time PCR (qRT-PCR) analysis validated that the expression level of those genes was significantly lower in the rubescent, etiolated, and albino mutant compared to wild type plants, which concurs with the differences in micro- and ultra-structures and physiological features between these two types of plants. Conclusively, the leaf color formation is greatly affected by the activity of chloroplast development and pigment biosynthesis. And the possible formation pathway of leaf color mutant of A. andraeanum 'Sonate' is deduced based on our results.

A deep transcriptomic analysis of pod development in the vanilla orchid (Vanilla planifolia)

BACKGROUND

Pods of the vanilla orchid (Vanilla planifolia) accumulate large amounts of the flavor compound vanillin (3-methoxy, 4-hydroxy-benzaldehyde) as a glucoside during the later stages of their development. At earlier stages, the developing seeds within the pod synthesize a novel lignin polymer, catechyl (C) lignin, in their coats. Genomic resources for determining the biosynthetic routes to these compounds and other flavor components in V. planifolia are currently limited.

RESULTS

Using next-generation sequencing technologies, we have generated very large gene sequence datasets from vanilla pods at different times of development, and representing different tissue types, including the seeds, hairs, placental and mesocarp tissues. This developmental series was chosen as being the most informative for interrogation of pathways of vanillin and C-lignin biosynthesis in the pod and seed, respectively. The combined 454/Illumina RNA-seq platforms provide both deep sequence coverage and high quality de novo transcriptome assembly for this non-model crop species.

CONCLUSIONS

The annotated sequence data provide a foundation for understanding multiple aspects of the biochemistry and development of the vanilla bean, as exemplified by the identification of candidate genes involved in lignin biosynthesis. Our transcriptome data indicate that C-lignin formation in the seed coat involves coordinate expression of monolignol biosynthetic genes with the exception of those encoding the caffeoyl coenzyme A 3-O-methyltransferase for conversion of caffeoyl to feruloyl moieties. This database provides a general resource for further studies on this important flavor species.

De novo transcriptome characterization of Lilium 'Sorbonne' and key enzymes related to the flavonoid biosynthesis

Lily is an important cut-flower and bulb crop in the commercial market. Here, transcriptome profiling of Lilium 'Sorbonne' was conducted through de novo sequencing based on Illumina platform. This research aims at revealing basic information and data that can be used for applied purposes especially the molecular regulatory information on flower color formation in lily. In total, 36,920,680 short reads which corresponded to 3.32 GB of total nucleotides, were produced through transcriptome sequencing. These reads were assembled into 39,636 Unigenes, of which 30,986 were annotated in Nr, Nt, Swiss-Prot, KEGG, COG, GO databases. Based on the three public protein databases, a total of 32,601 coding sequences were obtained. Meanwhile, 19,242 Unigenes were assigned to 128 KEGG pathways. Those with the greatest representation by unique sequences were for ''metabolic pathways'' (5,406 counts, 28.09 %). Our transcriptome revealed 156 Unigenes that encode key enzymes in the flavonoid biosynthesis pathway including CHS, CHI, F3H, FLS, DFR, etc. MISA software identified 2,762 simple sequence repeats, from which 1,975 primers pairs were designed. Over 2,762 motifs were identified, of which the most frequent was AG/CT (659, 23.86 %), followed by A/T (615, 22.27 %) and CCG/CGG (416, 15.06 %). Based on the results, we believe that the color formation of the Lilium 'Sorbonne' flower was mainly controlled by the flavonoid biosynthesis pathway. Additionally, this research provides initial genetic resources that will be valuable to the lily community for other molecular biology research, and the SSRs will facilitate marker-assisted selection in lily breeding.

A comparison of next generation sequencing technologies for transcriptome assembly and utility for RNA-Seq in a non-model bird

De novo assembled transcriptomes, in combination with RNA-Seq, are powerful tools to explore gene sequence and expression level in organisms without reference genomes. Investigators must first choose which high throughput sequencing platforms will provide data most suitable for their experimental goals. In this study, we explore the utility of 454 and Illumina sequences for de novo transcriptome assembly and downstream RNA-Seq applications in a reproductive gland from a non-model bird species, the Japanese quail (Coturnix japonica). Four transcriptomes composed of either pure 454 or Illumina reads or mixtures of read types were assembled and evaluated for the same cost. Illumina assemblies performed best for de novo transcriptome characterization in terms of contig length, transcriptome coverage, and complete assembly of gene transcripts. Improvements over the Hybrid assembly were marginal, with the exception that the addition of 454 data significantly increased the number of genes annotated. The Illumina assembly provided the best reference to align an independent set of RNA-Seq data as ∼84% of reads mapped to single genes in the transcriptome. Contigs constructed solely from 454 data may impose problems for RNA-Seq as our 454 transcriptome revealed a high number of indels and many ambiguously mapped reads. Correcting the 454 transcriptome with Illumina reads was an effective strategy to deal with indel and frameshift errors inherent to the 454 transcriptome, but at the cost of transcriptome coverage. In the absence of a reference genome, we find that Illumina reads alone produced a high quality transcriptome appropriate for RNA-Seq gene expression analyses.

De novo transcriptome assembly from inflorescence of Orchis italica: analysis of coding and non-coding transcripts

The floral transcriptome of Orchis italica, a wild orchid species, was obtained using Illumina RNA-seq technology and specific de novo assembly and analysis tools. More than 100 million raw reads were processed resulting in 132,565 assembled transcripts and 86,079 unigenes with an average length of 606 bp and N50 of 956 bp. Functional annotation assigned 38,984 of the unigenes to records present in the NCBI non-redundant protein database, 32,161 of them to Gene Ontology terms, 15,775 of them to Eukaryotic Orthologous Groups (KOG) and 7,143 of them to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The in silico expression analysis based on the Fragments Per Kilobase of transcript per Million mapped reads (FPKM) was confirmed by real-time RT-PCR experiments on 10 selected unigenes, which showed high and statistically significant positive correlation with the RNA-seq based expression data. The prediction of putative long non-coding RNAs was assessed using two different software packages, CPC and Portrait, resulting in 7,779 unannotated unigenes that matched the threshold values for both of the analyses. Among the predicted long non-coding RNAs, one is the homologue of TAS3, a long non-coding RNA precursor of trans-acting small interfering RNAs (ta-siRNAs). The differential expression pattern observed for the selected putative long non-coding RNAs suggests their possible functional role in different floral tissues.

BeMADS1 is a key to delivery MADSs into nucleus in reproductive tissues-De novo characterization of Bambusa edulis transcriptome and study of MADS genes in bamboo floral development

BACKGROUND

The bamboo Bambusa edulis has a long juvenile phase in situ, but can be induced to flower during in vitro tissue culture, providing a readily available source of material for studies on reproductive biology and flowering. In this report, in vitro-derived reproductive and vegetative materials of B. edulis were harvested and used to generate transcriptome databases by use of two sequencing platforms: Illumina and 454. Combination of the two datasets resulted in high transcriptome quality and increased length of the sequence reads. In plants, many MADS genes control flower development, and the ABCDE model has been developed to explain how the genes function together to create the different whorls within a flower.

RESULTS

As a case study, published floral development-related OsMADS proteins from rice were used to search the B. edulis transcriptome datasets, identifying 16 B. edulis MADS (BeMADS). The BeMADS gene expression levels were determined qRT-PCR and in situ hybridization. Most BeMADS genes were highly expressed in flowers, with the exception of BeMADS34. The expression patterns of these genes were most similar to the rice homologs, except BeMADS18 and BeMADS34, and were highly similar to the floral development ABCDE model in rice. Transient expression of MADS-GFP proteins showed that only BeMADS1 entered leaf nucleus. BeMADS18, BeMADS4, and BeMADS1 were located in the lemma nucleus. When co-transformed with BeMADS1, BeMADS15, 16, 13, 21, 6, and 7 translocated to nucleus in lemmas, indicating that BeMADS1 is a key factor for subcellular localization of other BeMADS.

CONCLUSION

Our study provides abundant B. edulis transcriptome data and offers comprehensive sequence resources. The results, molecular materials and overall strategy reported here can be used for future gene identification and for further reproductive studies in the economically important crop of bamboo.