Transcriptome analysis of Cymbidium sinense and its application to the identification of genes associated with floral development

Sex-Related Gene Network Revealed by Transcriptome Differentiation of Bisexual and Unisexual Flowers of Orchid Cymbidium tortisepalum

Despite extensive research on orchid reproductive strategies, the genetic studies of sex differentiation in the orchid family are still lacking. In this study, we compared three sexual phenotypes of Cymbidium tortisepalum bisexual flowers as well as female and male unisexual mutants. Through comparative transcriptomes, we analyzed the sex-biased differentially expressed genes (DEGs) and gene co-expression networks of sex organs (gynostemium and ovary) among them, identified the candidate genes of sex differentiation, and validated their expression by qRT-PCR. The C. tortisepalum unisexual mutants with degenerated phenotypes were compared to the bisexual plants with respect to both the flower organs and plant morphologies. Totally, 12,145, 10,789, and 14,447 genes were uniquely expressed in the female, male, and hermaphrodite sex organs, respectively. A total of 4291 sex-biased DEGs were detected among them, with 871, 2867, and 1937 DEGs in the comparisons of bisexual vs. female, bisexual vs. male, and male vs. female flowers, respectively. Two co-expressed network modules, with 81 and 419 genes were tightly correlated with female sexual traits, while two others with 265 and 135 genes were highly correlated with male sexual traits. Two female-biased hub genes (CtSDR3b and CtSDR3b-like) nested in the female modules, the homologs of maize sex determinant tasselseed2, may control the feminization of C. tortisepalum. At the same time, two male-biased hub genes (CtYAB2 and CtYAB5) nested in the male modules, the homologs of grape sex determinant VviYABBY3, may control the androphany of C. tortisepalum. This study discovered the molecular regulation networks and proposed a model for orchid sex differentiation, therefore providing for the first time the genetic basis for the sex separation in the orchid family.

Transcriptomic and Hormone Analyses Provide Insight into the Regulation of Axillary Bud Outgrowth of Eucommia ulmoides Oliver

An essential indicator of Eucommia ulmoides Oliver (E. ulmoides) is the axillary bud; the growth and developmental capacity of axillary buds could be used to efficiently determine the structural integrity of branches and plant regeneration. We obtained axillary buds in different positions on the stem, including upper buds (CK), tip buds (T1), and bottom buds (T2), which provided optimal materials for the study of complicated regulatory networks that control bud germination. This study used transcriptomes to analyze the levels of gene expression in three different types of buds, and the results showed that 12,131 differentially expressed genes (DEGs) were discovered via the pairwise comparison of transcriptome data gathered from CK to T2, while the majority of DEGs (44.38%) were mainly found between CK and T1. These DEGs were closely related to plant hormone signal transduction and the amino acid biosynthesis pathway. We also determined changes in endogenous hormone contents during the process of bud germination. Interestingly, except for indole-3-acetic acid (IAA) content, which showed a significant upward trend (p < 0.05) in tip buds on day 4 compared with day 0, the other hormones showed no significant change during the process of germination. Then, the expression patterns of genes involved in IAA biosynthesis and signaling were examined through transcriptome analysis. Furthermore, the expression levels of genes related to IAA biosynthesis and signal transduction were upregulated in tip buds. Particularly, the expression of the IAA degradation gene Gretchen Hagen 3 (GH3.1) was downregulated on day 4, which may support the concept that endogenous IAA promotes bud germination. Based on these data, we propose that IAA synthesis and signal transduction lead to morphological changes in tip buds during the germination process. On this basis, suggestions to improve the efficiency of the production and application of E. ulmoides are put forward to provide guidance for future research.

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

BACKGROUND

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

RESULT

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

CONCLUSIONS

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

Apelin-13 Protects Neurons by Attenuating Early-Stage Postspinal Cord Injury Apoptosis In Vitro

Apelin is a 77-amino-acid peptide that is an endogenous ligand for the G protein-coupled receptor APJ (Apelin receptor, APJ). Apelin-13, as the most bioactive affinity fragment of apelin, plays a role in energy metabolism, myocardial ischemia-reperfusion injury, and the regulation of the inflammatory response during oxidative stress, but its role in spinal cord injury is still unclear. This research identified and verified the differential expression of apelin in rat spinal cord injured tissues and normal spinal cord tissues by transcriptome sequencing in vivo and proved that apelin-13 protects neurons by strengthening autophagy and attenuating early-stage postspinal cord injury apoptosis in vitro. After constructing the model concerning a rat spinal cord hemisection damage, transcriptome sequencing was performed on the injured and normal spinal cord tissues of rats, which identified the differentially expressed gene apelin, with qRT-PCR detecting the representative level of apelin. The oxygen-glucose deprivation (OGD) model of PC12 cells was constructed in vitro to simulate spinal cord injury. The OGD injury times were 2 h, 4 h, 6 h, 8 h, and 12 h, and the non-OGD injury group was used as the control. The expression of apelin at each time point was observed by Western blotting. The expression of apelin was the lowest in the 6 h OGD injury group (p < 0.05). Therefore, the OGD injury time of 6 h was used in subsequent experiments. The noncytotoxic drug concentration of apelin-13 was determined with a Cell Counting Kit-8 (CCK-8) assay. An appropriate dose of apelin-13 (1 μM) significantly improved cell survival (p < 0.05). Thus, subsequent experiments selected a concentration of 1 μM apelin-13 as it significantly increased cell viability. Finally, we divided the experimental groups into four groups according to whether they received drugs (1 μM apelin-13, 24 h) or OGD (6 h): (1) control group: without apelin-13 or OGD injury; (2) apelin-13 group: with apelin-13 but no OGD injury; (3) OGD group: with OGD injury but without apelin-13; and (4) OGD + apelin-13 group: with apelin-13 and OGD injury. The TUNEL assay and flow cytometry results showed that compared with the OGD group, apoptosis in the OGD+Apelin-13 group was significantly reduced (p < 0.001). Determination of cell viability under different conditions by CCK-8 assay results displays that Apelin-13 can significantly improve the cell viability percentage under OGD conditions (p < 0.001). Western blotting results showed that apelin-13 decreased the expression ratios of apoptosis-related proteins Bax/Bcl-2 and cleaved-caspase3/caspase3 (p < 0.05), increasing the key to Beclin1-dependent autophagy pathway expression of the protein Beclin1. This finding indicates that apelin-13 protects neurons by strengthening autophagy and attenuating early-stage postspinal cord injury apoptosis in vitro.

De Novo Transcriptome Assembly and Analysis of Longevity Genes Using Subterranean Termite (Reticulitermes chinensis) Castes

The longevity phenomenon is entirely controlled by the insulin signaling pathway (IIS-pathway). Both vertebrates and invertebrates have IIS-pathways that are comparable to one another, though no one has previously described de novo transcriptome assembly of IIS-pathway-associated genes in termites. In this research, we analyzed the transcriptomes of both reproductive (primary kings “PK” and queens “PQ”, secondary worker reproductive kings “SWRK” and queens “SWRQ”) and non-reproductive (male “WM” and female “WF” workers) castes of the subterranean termite Reticulitermes chinensis. The goal was to identify the genes responsible for longevity in the reproductive and non-reproductive castes. Through transcriptome analysis, we annotated 103,589,264 sequence reads and 184,436 (7G) unigenes were assembled, GC performance was measured at 43.02%, and 64,046 sequences were reported as CDs sequences. Of which 35 IIS-pathway-associated genes were identified, among 35 genes, we focused on the phosphoinositide-dependent kinase-1 (Pdk1), protein kinase B2 (akt2-a), tuberous sclerosis-2 (Tsc2), mammalian target of rapamycin (mTOR), eukaryotic translation initiation factor 4E (EIF4E) and ribosomal protein S6 (RPS6) genes. Previously these genes (Pdk1, akt2-a, mTOR, EIF4E, and RPS6) were investigated in various organisms, that regulate physiological effects, growth factors, protein translation, cell survival, proliferation, protein synthesis, cell metabolism and survival, autophagy, fecundity rate, egg size, and follicle number, although the critical reason for longevity is still unclear in the termite castes. However, based on transcriptome profiling, the IIS-pathway-associated genes could prolong the reproductive caste lifespan and health span. Therefore, the transcriptomic shreds of evidence related to IIS-pathway genes provide new insights into the maintenance and relationships between biomolecular homeostasis and remarkable longevity. Finally, we propose a strategy for future research to decrypt the hidden costs associated with termite aging in reproductive and non-reproductive castes.

First Report of stem rot in Cymbidium sinense Caused by Fusarium oxysporum (FOSC) in China

Cymbidium sinense (Jackson ex Andr.) Willd is a perennial terrestrial plant in the orchid family mainly distributed in China, Japan, India and Southeast Asia that occupies a strong position in the flower market due to its bright green leaves and fragrant flowers (Zhang et al. 2013). Cymbidium sinense is not only valued by people for its ornamental and economic value, but its roots have antiasthmatic medicinal properties (Ke et al. 2004). In August 2020, about 15% stem rot on two-year old C. sinense with varying severity was observed in five nursery gardens located in Enshi city (N 30° 16', E 109° 29'), Hubei province, China. Typical symptoms of C. sinense included roots and inner part of the pseudobulbs changing from white to brown and rotting. Leaves became brown and withered from bottom to top, and there was an obvious blight yellow halo at the junction of diseased and healthy tissue, which eventually caused the whole plant to wilt and die (Fig. 1d). To isolate the pathogen, a total of 15 leaf tissues from the disease-health junction (3 × 3 mm) from 5 individual plants (3 leaves/plant) with symptoms were surface sterilized with 75% ethanol for 30 s and 2% sodium hypochlorite (NaOCl) for 3 min. The sterilized tissue was rinsed three times with sterilized water, and then placed on potato dextrose agar (PDA) for incubation at 28°C in the dark for 5 days. Isolated colonies were subcultured by a hyphal tip protocol. Thirteen fungal isolates were obtained. Through preliminary pathogenicity tests, we found that ten isolates induced leaf blight. These ten isolates with pathogenicity showed similar morphological characteristics, with initial white-flocculent aerial mycelium that secreted a lavender pigment and produced colonies with an irregular edge after 3 days on PDA. The ten strains were cultured on PDA plates at 28℃ for 5 and 15 days to observe colony and conidial characteristics. The ten strains were identified as Fusarium based on morphological characteristics (Leslie and Summerell 2006). Strain ML0303 was selected for further identification. Macroconidia were falciform, hyaline, slightly pointed at both ends with two to four septa, 24.0 ± 5.6 µm × 4.7 ± 0.8 µm (n = 50). Microconidia were hyaline, oval, globose, with zero to one septum, 5.5 ± 1.3 µm × 2.2 ± 0.5 µm (n = 50) (Fig. 1c). Total genomic DNA of strain ML0303 was extracted with a CTAB protocol (Stenglein and Balatti 2006). The translation elongation factor (EF-1α), RNA polymerase II second largest subunit (RPB2) and β-tubulin (Tub2) genes were amplified respectively using primer pairs EF1/EF2, RPB2-5F2/RPB2-7cR and T1/T22 respectively (O'Donnell. et al. 2010, O'Donnell. et al. 1997). The EF-1α, RPB2 and Tub2 (accession numbers-MW719874, OL614838, OL689398, respectively) gene sequences were submitted to GenBank. EF-1α, RPB2 and Tub2 sequences of ML0303 showed 99.5% - 100% identity respectively with Fusarium oxysporum in the Genbank and FUSARIUM-ID databases. The multilocus sequence data was used to infer a phylogenetic tree via a Neighbor-joining (NJ), Maximum-likelihood (ML) and Maximum-Parsimony(MP) together with reference sequences from GenBank. The topology of the three trees was similar; only the NJ tree is presented here. Strain ML0303 and F. oxysporum formed a clade supported with high values (NJ/ML/MP: 96,95,97). The results indicated that the fungus was F. oxysporum based on the phylogenetic analysis and BLASTn queries. For pathogenicity tests, conidia of strain ML0303 were collected by rinsing PDA plates. Two-year-old C. sinense grown in plastic pots filled with sterilized autoclaved sandy loam soil were used for the tests. Three pots (two plants/pot) were included in each treatment. Spore suspensions (106spores/ml) of strain ML0303 were used to irrigate the stem-zone of the plants, and sterile water was used as control. The two treatments were placed in a greenhouse and incubated at 28±2℃ with a 14-hour light/10-hour dark cycle. The experiment was repeated twice. After three weeks, stem rot symptoms were observed on C. sinense inoculated with ML0303, that were the as same as observed in the nursery (Fig. 1e-h). No symptoms were observed on the negative control. Fusarium oxysporum was re-isolated from the infected plants to fulfill Koch's postulates. Partial EF-1α and RPB2 gene sequences were used for molecular identification. Members of the FOSC are notorious for causing many diseases, which includes stem rot of Sulcorebutia heliosa and root rot of Torreya grandis (Garibaldi et al. 2020; Zhang et al. 2016). To our knowledge, this is the first report of stem rot by F. oxysporum on C. sinense in China. The finding of this pathogen provides a clear target for stem rot control.

The Antennal Sensilla and Expression Patterns of Olfactory Genes in the Lower Termite Reticulitermes aculabialis (Isoptera: Rhinotermitidae)

The insect olfactory system plays pivotal roles in insect survival and reproduction through odor detection. Morphological and physiological adaptations are caste-specific and evolved independently in workers, soldiers, and reproductives in termites. However, it is unclear whether the olfactory system is involved in the division of labor in termite colonies. In the present study, the antennal sensilla of alates, workers, soldiers, nymphs, and larvae of the termite Reticulitermes aculabialis Tsai et Hwang ( Isoptera: Rhinotermitidae) were investigated. Transcriptomes were used to detect olfactory genes, and differential expression levels of olfactory genes were confirmed in various castes by qRT-PCR analysis. Nine types of sensilla were identified on the antennae of R. aculabialis, and soldiers possessed all 9 types. In 89,475 assembled unigenes, we found 16 olfactory genes, including 6 chemosensory protein (CSP) and 10 odorant-binding protein (OBP) genes. These OBP genes included 8 general odorant-binding protein genes (GOBPs) and 2 pheromone-binding protein-related protein (PBP) genes. Five CSP genes were more highly expressed in alates than in workers, soldiers, larvae, and nymphs, and the expression levels of CSP6 were significantly higher in nymphs. Seven GOBP and two PBP genes exhibited significantly higher expression levels in alates, and there were no significant differences in the expression levels of GOBP2 among workers, soldiers, alates, and larvae. These results suggest that alates, as primary reproductives, have unique expression patterns of olfactory genes, which play key roles in nuptial flight, mate seeking, and new colony foundation.

Orchid Phylotranscriptomics: The Prospects of Repurposing Multi-Tissue Transcriptomes for Phylogenetic Analysis and Beyond

The Orchidaceae is rivaled only by the Asteraceae as the largest plant family, with the estimated number of species exceeding 25,000 and encompassing more than 700 genera. To gain insights into the mechanisms driving species diversity across both global and local scales, well-supported phylogenies targeting different taxonomic groups and/or geographical regions will be crucial. High-throughput sequencing technologies have revolutionized the field of molecular phylogenetics by simplifying the process of obtaining genome-scale sequence data. Consequently, there has been an explosive growth of such data in public repositories. Here we took advantage of this unprecedented access to transcriptome data from predominantly non-phylogenetic studies to assess if it can be repurposed to gain rapid and accurate phylogenetic insights across the orchids. Exhaustive searches revealed transcriptomic data for more than 100 orchid species spanning 5 subfamilies, 13 tribes, 21 subtribes, and 50 genera that were amendable for exploratory phylotranscriptomic analysis. Next, we performed re-assembly of the transcriptomes before strategic selection of the final samples based on a gene completeness evaluation. Drawing on these data, we report phylogenetic analyses at both deep and shallow evolutionary scales via maximum likelihood and shortcut coalescent species tree methods. In this perspective, we discuss some key outcomes of this study and conclude by highlighting other complementary, albeit rarely explored, insights beyond phylogenetic analysis that repurposed multi-tissue transcriptome can offer.

An Integrated Analysis of Transcriptome and miRNA Sequencing Provides Insights into the Dynamic Regulations during Flower Morphogenesis in Petunia

Published genome sequences can facilitate multiple genome sequencing studies of flower development, which can serve as the basis for later analysis of variation in flower phenotypes. To identify potential regulators related to flower morphology, we captured dynamic expression patterns under five different developmental stages of petunia flowers, a popular bedding plant, using transcriptome and miRNA sequencing. The significant transcription factor (TF) families, including MYB, MADS, and bHLH, were elucidated. MADS-box genes exhibited co-expression patterns with BBR-BPC, GATA, and Dof genes in different modules according to a weighted gene co-expression network analysis. Through miRNA sequencing, a total of 45 conserved and 26 novel miRNAs were identified. According to GO and KEGG enrichment analysis, the carbohydrate metabolic process, photosynthesis, and phenylalanine metabolism were significant at the transcriptomic level, while the response to hormone pathways was significantly enriched by DEmiR-targeted genes. Finally, an miRNA–RNA network was constructed, which suggested the possibility of novel miRNA-mediated regulation pathways being activated during flower development. Overall, the expression data in the present study provide novel insights into the developmental gene regulatory network facilitated by TFs, miRNA, and their target genes.

The Genetic and Hormonal Inducers of Continuous Flowering in Orchids: An Emerging View

Orchids are the flowers of magnetic beauty. Vivid and attractive flowers with magnificent shapes make them the king of the floriculture industry. However, the long-awaited flowering is a drawback to their market success, and therefore, flowering time regulation is the key to studies about orchid flower development. Although there are some rare orchids with a continuous flowering pattern, the molecular regulatory mechanisms are yet to be elucidated to find applicable solutions to other orchid species. Multiple regulatory pathways, such as photoperiod, vernalization, circadian clock, temperature and hormonal pathways are thought to signalize flower timing using a group of floral integrators. This mini review, thus, organizes the current knowledge of floral time regulators to suggest future perspectives on the continuous flowering mechanism that may help to plan functional studies to induce flowering revolution in precious orchid species.

Genetic insights into the regulatory pathways for continuous flowering in a unique orchid Arundina graminifolia

BACKGROUND

Manipulation of flowering time and frequency of blooming is key to enhancing the ornamental value of orchids. Arundina graminifolia is a unique orchid that flowers year round, although the molecular basis of this flowering pattern remains poorly understood.

RESULTS

We compared the A. graminifolia transcriptome across tissue types and floral developmental stages to elucidate important genetic regulators of flowering and hormones. Clustering analyses identified modules specific to floral transition and floral morphogenesis, providing a set of candidate regulators for the floral initiation and timing. Among candidate floral homeotic genes, the expression of two FT genes was positively correlated with flower development. Assessment of the endogenous hormone levels and qRT-PCR analysis of 32 pathway-responsive genes supported a role for the regulatory networks in floral bud control in A. graminifolia. Moreover, WGCNA showed that flowering control can be delineated by modules of coexpressed genes; especially, MEgreen presented group of genes specific to flowering.

CONCLUSIONS

Candidate gene selection coupled with hormonal regulators brings a robust source to understand the intricate molecular regulation of flowering in precious orchids.

The genome of Cymbidium sinense revealed the evolution of orchid traits

The Orchidaceae is of economic and ecological importance and constitutes ˜10% of all seed plant species. Here, we report a genome physical map for Cymbidium sinense, a well-known species belonging to genus Cymbidium that has thousands of natural variation varieties of flower organs, flower and leaf colours and also referred as the King of Fragrance, which make it arose into a unique cultural symbol in China. The high-quality chromosome-scale genome assembly was 3.52 Gb in size, 29 638 protein-coding genes were predicted, and evidence for whole-genome duplication shared with other orchids was provided. Marked amplification of cytochrome- and photosystem-related genes was observed, which was consistent with the shade tolerance and dark green leaves of C. sinense. Extensive duplication of MADS-box genes, and the resulting subfunctional and expressional differentiation, was associated with regulation of species-specific flower traits, including wild-type and mutant-type floral patterning, seasonal flowering and ecological adaption. CsSEP4 was originally found to positively regulate gynostemium development. The CsSVP genes and their interaction proteins CsAP1 and CsSOC1 were significantly expanded and involved in the regulation of low-temperature-dependent flowering. Important genetic clues to the colourful leaf traits, purple-black flowers and volatile trait in C. sinense were also found. The results provide new insights into the molecular mechanisms of important phenotypic traits of Cymbidium and its evolution and serve as a powerful platform for future evolutionary studies and molecular breeding of orchids.

Transcriptome Analysis of Flower Development and Mining of Genes Related to Flowering Time in Tomato (Solanum lycopersicum)

Flowering is a morphogenetic process in which angiosperms shift from vegetative growth to reproductive growth. Flowering time has a strong influence on fruit growth, which is closely related to productivity. Therefore, research on crop flowering time is particularly important. To better understand the flowering period of the tomato, we performed transcriptome sequencing of early flower buds and flowers during the extension period in the later-flowering “Moneymaker” material and the earlier-flowering “20965” homozygous inbred line, and we analyzed the obtained data. At least 43.92 million clean reads were obtained from 12 datasets, and the similarity with the tomato internal reference genome was 92.86–94.57%. Based on gene expression and background annotations, 49 candidate genes related to flowering time and flower development were initially screened, among which the greatest number belong to the photoperiod pathway. According to the expression pattern of candidate genes, the cause of early flowering of “20965” is predicted. The modes of action of the differentially expressed genes were classified, and the results show that they are closely related to hormone regulation and participated in a variety of life activities in crops. The candidate genes we screened and the analysis of their expression patterns provide a basis for future functional verification, helping to explore the molecular mechanism of tomato flowering time more comprehensively.

Comparative Transcriptomic Analysis of Differentially Expressed Transcripts Associated with Flowering Time of Loquat (Eriobotya japonica Lindl.)

Flowering is an important phenophase of plant species, however, knowledge about the regulatory mechanism controlling flowering cues in loquat is limited. To identify candidate genes regulating flowering time in loquat, we used RNA-Seq technology to conduct a comparative transcriptome analysis of differentiating apical buds collected from the early-flowering variety ‘Baiyu’ and the late-flowering variety ‘Huoju’. A total of 28,842 differentially expressed transcripts (DETs) were identified. Of these, 42 DETs controlled flowering time while 17 other DETs were associated with the ABA signaling pathway. Compared with those in ‘Huoju’, EjFT, EjFY, EjFLK, and EjCAL1-like were significantly upregulated in ‘Baiyu’. Moreover, transcripts of the ABA 8′-hydroxylases (EjABH2, EjABH4, and EjABH4-like2), the ABA receptors (EjPYL4/8), and the bZIP transcription factor EjABI5-like were upregulated in ‘Baiyu’ compared with ‘Huoju’. Hence, they might regulate loquat flowering time. There was no significant difference between ‘Baiyu’ and ‘Huoju’ in terms of IAA content. However, the ABA content was about ten-fold higher in the apical buds of ‘Baiyu’ than in those of ‘Huoju’. The ABA:IAA ratio sharply rose and attained a peak during bud differentiation. Thus, ABA is vital in regulating floral bud formation in loquat. The results of the present study help clarify gene transcription during loquat flowering.

Transcriptional Cascade in the Regulation of Flowering in the Bamboo Orchid Arundina graminifolia

Flowering in orchids is the most important horticultural trait regulated by multiple mechanisms. Arundina graminifolia flowers throughout the year unlike other orchids with a narrow flowering span. However, little is known of the genetic regulation of this peculiar flowering pattern. This study identifies a number of transcription factor (TF) families in five stages of flower development and four tissue types through RNA-seq transcriptome. About 700 DEGs were annotated to the transcription factor category and classified into 35 TF families, which were involved in multiple signaling pathways. The most abundant TF family was bHLH, followed by MYB and WRKY. Some important members of the bHLH, WRKY, MYB, TCP, and MADS-box families were found to regulate the flowering genes at transcriptional levels. Particularly, the TFs WRKY34 and ERF12 possibly respond to vernalization and photoperiod signaling, MYB108, RR9, VP1, and bHLH49 regulate hormonal balance, and CCA1 may control the circadian pathway. MADS-box TFs including MADS6, 14, 16, AGL5, and SEP may be important regulators of flowering in A. graminifolia. Therefore, this study provides a theoretical basis for understanding the molecular mechanism of flowering in A. graminifolia.

Comparative Transcriptomic and Metabolic Analyses Reveal the Molecular Mechanism of Ovule Development in the Orchid, Cymbidium sinense

Ovule development is pivotal to plant reproduction and seed development. Cymbidium sinense (Orchidaceae) has high ornamental value due to its pleasant aroma and elegant floral morphology. The regulatory mechanism underlying ovule development in orchids, especially C. sinense, is largely unknown and information on the C. sinense genome is very scarce. In this study, a combined analysis was performed on the transcriptome and non-targeted metabolomes of 18 C. sinense 'Qi Jian Hei Mo' ovule samples. Transcriptome analysis assembled gene-related information related to six growth stages of C. sinense ovules (S1-S6, equivalent to 30, 35, 42, 46, 53, and 60 days after pollination). Illumina sequencing technology was used to obtain the complete set of transcriptome sequences of the 18 samples. A total of 81,585 unigene sequences were obtained after assembly, 24,860 (30.47%) of which were functionally annotated. Using transcriptome sequencing technology, a total of 9845 differentially expressed unigenes (DEUs) were identified in C. sinense ovules that were assigned to specific metabolic pathways according to the Kyoto Encyclopedia of Genes and Genomes (KEGG). DEUs associated with transcription factors (TFs) and phytohormones were identified and analyzed. The TFs homeobox and MADS-box were associated with C. sinense ovule development. In particular, the phytohormones associated with DEUs such as indole-3-acetic acid (IAA), cytokinin (CK), gibberellin (GA), abscisic acid (ABA), brassinosteroid (BR), and jasmonate (JA), may have important regulatory effects on C. sinense ovule development. Metabolomic analysis showed an inconsistent number of KEGG annotations of differential metabolites across comparisons (S2_vs_S4, S2_vs_S5, and S4_vs_S5 contained 23, 26, and 3 annotations, respectively) in C. sinense ovules. This study provides a valuable foundation for further understanding the regulation of orchid ovule development and formation, and establishes a theoretical background for future practical applications during orchid cultivation.

Comparative transcriptome analysis reveals the regulatory effects of acetylcholine on salt tolerance of Nicotiana benthamiana

Salinity is a major cause of crop losses worldwide. Acetylcholine (ACh) can ameliorate the adverse effects of abiotic stresses on plant growth, including salinity stress; however, the underlying molecular mechanisms of this process are unclear. Here, seedlings of Nicotiana benthamiana grown under normal conditions or exposed to 150 mmol L-1 NaCl salinity stress were then treated with a root application of 10 μM ACh. Exogenous ACh application resulted in the downregulation of the activity of the antioxidant enzymes, ascorbate peroxidase, and catalase. ACh-treated plants had lower levels of reactive oxygen species, including the superoxide anion radical and hydrogen peroxide. Transcriptome analysis indicated that ACh treatment under salt stress promoted the differential expression of 658 genes in leaves of N. benthamiana (527 were upregulated and 131 were downregulated). Gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed that exogenous ACh application was associated with a substantial increase in the transcripts of genes related to cell wall peroxidases, xyloglucan endotransglucosylases or hydrolases, and expansins, indicating that ACh activates cell wall biosynthesis in salt-stressed plants. ACh also enhanced the expression of genes associated with the auxin, gibberellin, brassinosteroid, and salicylic acid signalling pathways, indicating that ACh induces the activation of these pathways under salt stress. Collectively, these findings indicate that ACh-induced salt tolerance in N. benthamiana seedlings is mediated by the inhibition of antioxidant enzymes, activation of cell wall biosynthesis, and hormone signalling pathways. Stress-induced genes involved in osmotic regulation and oxidation resistance were induced by ACh under salt stress. The genes whose transcript levels were elevated by ACh treatment in salt-stressed N. benthamiana could be used as molecular markers of the physiological status of plants under salt stress.

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

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

De novo assembly and Transcriptome characterization of an endemic species of Vietnam, Panax vietnamensis Ha et Grushv., including the development of EST-SSR markers for population genetics

BACKGROUND

Understanding the genetic diversity in endangered species that occur inforest remnants is necessary to establish efficient strategies for the species conservation, restoration and management. Panax vietnamensis Ha et Grushv. is medicinally important, endemic and endangered species of Vietnam. However, genetic diversity and structure of population are unknown due to lack of efficient molecular markers.

RESULTS

In this study, we employed Illumina HiSeq™ 4000 sequencing to analyze the transcriptomes of P. vietnamensis (roots, leaves and stems). Raw reads total of 23,741,783 was obtained and then assembled, from which the generated unigenes were 89,271 (average length = 598.3191 nt). The 31,686 unigenes were annotated in different databases i.e. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Nucleotide Collection (NR/NT) and Swiss-Prot for functional annotation. Further, 11,343 EST-SSRs were detected. From 7774 primer pairs, 101 were selected for polymorphism validation, in which; 20 primer pairs were successfully amplified to DNA fragments and significant amounts of polymorphism was observed within population. The nine polymorphic microsatellite loci were used for population structure and diversity analyses. The obtained results revealed high levels of genetic diversity in populations, the average observed and expected heterozygosity were HO = 0.422 and HE = 0.479, respectively. During the Bottleneck analysis using TPM and SMM models (p < 0.01) shows that targeted population is significantly heterozygote deficient. This suggests sign of the bottleneck in all populations. Genetic differentiation between populations was moderate (FST = 0.133) and indicating slightly high level of gene flow (Nm = 1.63). Analysis of molecular variance (AMOVA) showed 63.17% of variation within individuals and 12.45% among populations. Our results shows two genetic clusters related to geographical distances.

CONCLUSION

Our study will assist conservators in future conservation management, breeding, production and habitats restoration of the species.

Involvement of CsERF2 in leaf variegation of Cymbidium sinense 'Dharma'

MAIN CONCLUSION

CsERF2, an ethylene response factor, plays a role in leaf variegation. Leaf variegation is a main ornamental characteristic in Cymbidium sinense (C. sinense). However, the mechanisms of leaf color variegation remain largely unclear. In the present study, we analyzed the cytological and physiological features, as well as molecular analyses of leaves from wild-type (WT) and leaf variegation mutants of Cymbidium sinense 'Dharma'. Chloroplasts with typical and functional structures were discovered in WT and green sectors of the mutants leaves (MG), but not in yellow sectors of the mutant leaves (MY). The activities of key enzymes involved in chlorophyll (Chl) degradation and their substrate contents were significantly increased in MY. Genes related to Chl degradation also showed a significant up-regulation in MY. Transcriptomic analysis showed that the expression of all identified ethylene response factors (ERFs) was significantly up-regulated, and the 1-aminocyclopropane-1-carboxylic acid (ACC) content in MY was significantly higher compared with MG. QRT-PCR analysis validated that the expression levels of genes related to Chl degradation could be positively affected by ethylene (ETH) treatment. Stable overexpression of CsERF2 in Nicotiana tabacum (N. tabacum) led to a decrease in Chl content and abnormal chloroplast. Transcriptomic analysis and qRT-PCR results showed that the KEGG pathway related to chloroplast development and function showed significant change in transgenic N. tabacum. Therefore, the leaf color formation of C. sinense was greatly affected by chloroplast development and Chl metabolism. CsERF2 played an important role in leaf variegation of C. sinense.

Transcriptomic evidence that insulin signalling pathway regulates the ageing of subterranean termite castes

Insulin is a protein hormone that controls the metabolism of sugar, fat and protein via signal transduction in cells, influencing growth and developmental processes such as reproduction and ageing. From nematodes to fruit flies, rodents and other animals, glucose signalling mechanisms are highly conserved. Reproductive termites (queens and kings) exhibit an extraordinarily long lifespan relative to non-reproductive individuals such as workers, despite being generated from the same genome, thus providing a unique model for the investigation of longevity. The key reason for this molecular mechanism, however, remains unclear. To clarify the molecular mechanism underlying this phenomenon, we sequenced the transcriptomes of the primary kings (PKs), primary queens (PQs), male (WMs) and female (WFs) workers of the lower subterranean termite Reticulitermes chinensis. We performed RNA sequencing and identified 33 insulin signalling pathway-related genes in R. chinensis. RT-qPCR analyses revealed that EIF4E and RPS6 genes were highly expressed in WMs and WFs workers, while mTOR expression was lower in PKs and PQs than in WMs and WFs. PQs and PKs exhibited lower expression of akt2-a than female workers. As the highly conserved insulin signalling pathway can significantly prolong the healthspan and lifespan, so we infer that the insulin signalling pathway regulates ageing in the subterranean termite R. chinensis. Further studies are recommended to reveal the biological function of insulin signalling pathway-related genes in the survival of termites to provide new insights into biomolecular homeostasis maintenance and its relationship to remarkable longevity.

Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast- and slow-growing fish

Grass carp is one of the most important freshwater aquaculture species in China. However, the mechanisms underlying the growth of muscle tissue in the fish are unclear. High-throughput RNA-Seq was used to analyze the transcriptome of grass carp muscle tissue between fast- and slow-growing fish family groups. Twenty-four individuals each from 4 fast-growing families and 4 slow-growing families were used to reduce background noise. 71 up-regulated and 35 down-regulated genes were identified in the differentially expressed genes (DEGs). GO and KEGG enrichment analyses revealed the DEGs were involved in the GH/IGF axis, calcium metabolism, protein and glycogen synthesis, oxygen transport, cytoskeletal and myofibrillar components. IGFBP1 was up-regulated in big fish while GHR2 was down-regulated. Glutamic pyruvate transaminase 2, an indicator of liver tissue damage, was down-regulated in big grass carp, which indicates that the fish was better adapted to an artificially formulated diet. GAPDH, the rate-limiting enzyme in glycolytic flux was highly expressed in fast-growing grass carp, reflecting enhanced carbohydrate metabolism. Higher expression of ALAS2 and myoglobin 1 in big grass carp, related to oxygen transport might promote aerobic exercise along with food intake and muscle growth. Genes for cytoskeletal and myofibrillar components such as tropomyosin, meromyosin, and troponin I were also up-regulated in big grass carp. These results provide valuable information about the key genes for use as biomarkers of growth in selective breeding programs for grass carp and contribute to our understanding of the molecular mechanisms and regulative pathways regulating growth in fish.

Transcriptome analysis to identify candidate genes associated with the yellow-leaf phenotype of a Cymbidium mutant generated by γ-irradiation

Leaf color is an important agronomic trait in flowering plants, including orchids. However, factors underlying leaf phenotypes in plants remain largely unclear. A mutant displaying yellow leaves was obtained by the γ-ray-based mutagenesis of a Cymbidium orchid and characterized using RNA sequencing. A total of 144,918 unigenes obtained from over 25 million reads were assigned to 22 metabolic pathways in the Kyoto Encyclopedia of Genes and Genomes database. In addition, gene ontology was used to classify the predicted functions of transcripts into 73 functional groups. The RNA sequencing analysis identified 2,267 differentially expressed genes between wild-type and mutant Cymbidium sp. Genes involved in the chlorophyll biosynthesis and degradation, as well as ion transport, were identified and assayed for their expression levels in wild-type and mutant plants using quantitative real-time profiling. No critical expression changes were detected in genes involved in chlorophyll biosynthesis. In contrast, seven genes involved in ion transport, including two metal ion transporters, were down-regulated, and chlorophyllase 2, associated with chlorophyll degradation, was up-regulated. Together, these results suggest that alterations in chlorophyll metabolism and/or ion transport might contribute to leaf color in Cymbidium orchids.

Functional annotation of mulberry (Morus spp.) transcriptome, differential expression of genes related to growth and identification of putative genic SSRs, SNPs and InDels

Growth is a complex trait associated with mulberry leaf yield and controlled by several genes. In this study, we have explored the molecular basis underlying growth using Transcriptome profiling of contrasting genotypes. A total of 66.6 Mbp of primary transcriptomes from high growth (HGG)-Jalalgarah-3 and M. laevigata (H) and, low growth genotypes (LGG)-Harmutty and Vadagaraparai-2; resulting in 24210, 27998, 28085 and 28764 final transcripts respectively. Out of the 34096 pooled transcripts, 20249 transcripts matched with at least one sequence of the non-redundant database. Functional annotation resulted in the categorization of 18970 transcripts into 3 gene ontology (GO) terms and 7440 were assigned to 23 Kyoto encyclopaedia of genes and genomes (KEGG) pathway. Based on the differentially expressed genes and gene enrichment analysis, over expression of photosynthetic related transcripts in HGG and defence related transcripts in LGG were noted. Simple sequence repeats were mined from unique transcripts and the most abundant motifs were tri- (1883) followed by di- (1710), tetra- (192), penta- (68) and hexa- (40) repeats. Further, a total of 390897 high quality SNPs and 8081 InDels were identified by mapping onto Morus notabilis reference genome. The study provides an insight into the expression of genes involved in growth and further research on utilization in gentic improvement of the crop.

RNA sequencing analysis of Cymbidium goeringii identifies floral scent biosynthesis related genes

BACKGROUND

Cymbidium goeringii belongs to the Orchidaceae, which is one of the most abundant angiosperm families. Cymbidium goeringii consist with high economic value and characteristics include fragrance and multiple flower colors. Floral scent is one of the important strategies for ensuring fertilization. However, limited genetic data is available in this non-model plant, and little known about the molecular mechanism responsible for floral scent in this orchid. Transcriptome and expression profiling data are needed to identify genes and better understand the biological mechanisms of floral scents in this species. Present transcriptomic data provides basic information on the genes and enzymes related to and pathways involved in flower secondary metabolism in this plant.

RESULTS

In this study, RNA sequencing analyses were performed to identify changes in gene expression and biological pathways related scent metabolism. Three cDNA libraries were obtained from three developmental floral stages: closed bud, half flowering stage and full flowering stage. Using Illumina technique 159,616,374 clean reads were obtained and were assembled into 85,868 final unigenes (average length 1194 nt), 33.85% of which were annotated in the NCBI non redundant protein database. Among this unigenes 36,082 were assigned to gene ontology and 23,164 were combined with COG groups. Total 33,417 unigenes were assigned in 127 pathways according to the Kyoto Encyclopedia of Genes and Genomes pathway database. According these transcriptomic data we identified number of candidates genes which differentially expressed in different developmental stages of flower related to fragrance biosynthesis. In q-RT-PCR most of the fragrance related genes highly expressed in half flowering stage.

CONCLUSIONS

RNA-seq and DEG data provided comprehensive gene expression information at the transcriptional level that could be facilitate the molecular mechanisms of floral biosynthesis pathways in three developmental phase's flowers in Cymbidium goeringii, moreover providing useful information for further analysis on C. goeringii, and other plants of genus Cymbidium.

De novo assembly of the Platycladus orientalis (L.) Franco transcriptome provides insight into the development and pollination mechanism of female cone based on RNA-Seq data

For seed-bearing plants, the basis of seed and fruit formation is pollination. The normal progression of pollination is through advances in continuous signal exchange and material transfer, which occur mainly in female reproductive organs; thus, the molecular mechanism of development in female reproductive organs is vital for understanding the principle of pollination. However, molecular biology studies on the development of female cones related to pollination are rare and unclear in gymnosperms, especially in Cupressaceae. In this study, Platycladus orientalis, a monotypic genus within Cupressaceae, was chosen to examine female cone transcriptomes at pre-pollination and pollination stages by Illumina paired-end sequencing technology to de novo sequence six libraries with 3 biological replicates. These libraries were used to construct a P. orientalis transcriptome database containing 71,669 unigenes (4,963 upregulated unigenes and 11,747 downregulated unigenes at the pollination stage) for subsequent analysis. Based on the annotations and expression levels, the functions of differentially expressed unigenes and enriched pathways between the developmental processes of female cones were analysed to detail the preliminary development and pollination mechanism of the female cone. Targeted investigations were specifically performed to determine the elementary mechanism of secretion and functioning of the pollination drop, a vital ovule secretion at the pollination stage. Ultimately, the expression of 15 unigenes selected between two stages were further assessed and confirmed using qRT-PCR, which demonstrated reliable data and significant differences in the expression profiles of key genes. As one of the largest available transcriptomic resources of this species, the database is constructed to prospectively adapt to the physiological and genomic data of woody plants. This work provided the first transcriptome profile of P. orientalis female cones at different developmental stages, and will promote the illumination of the pollination mechanism of P. orientalis, and will serve as the basis for in-depth genomic study in the Cupressaceae family. This initiative will arouse the interest and attention of scholars and pave the way for future studies.

Global transcriptome analysis uncovers the gene co-expression regulation network and key genes involved in grain development of wheat (Triticum aestivum L.)

Wheat grain development is a robust biological process that largely determines grain quality and yield. In this study, we investigated the grain transcriptome of winter wheat cv. Xiaoyan-6 at four developmental stages (5, 10, 15, and 20 days post-anthesis), using high-throughput RNA sequencing (RNA-Seq). We identified 427 grain-specific transcription factors (TFs) and 1653 differentially expressed TFs during grain development as well as a grain co-expression regulation network (GrainNet) of the TFs and their predicted co-expressed genes. Our study identified ten putative key TFs and the predicted regulatory genes of these TFs in wheat grain development of Xiaoyan-6. The analysis was given a firm basis through the study of additional wheat tissues, including root, stem, leaf, flag leaf, grain, spikes (from wheat plants at booting or heading stages) to provide a dataset of 92,478 high-confidence protein-coding genes that were mostly evenly distributed among subgenomes, but unevenly distributed across each of the chromosomes or each of the seven homeologous groups. Within this larger framework of the transcriptomes, we identified 4659 grain-specific genes (SEGs) and 26,500 differentially expressed genes (DEGs) throughout grain development stages tested. The SEGs identified mainly associate with regulation and signaling-related biological processes, while the DEGs mainly involve in cellular component organization or biogenesis and nutrient reservoir activity during grain development of Xiaoyan-6. This study establishes new targets for modifying genes related to grain development and yield, to fine-tune expression in different varieties and environments.

Transcriptional Study Revealed That Boron Supplementation May Alter the Immune-Related Genes Through MAPK Signaling in Ostrich Chick Thymus

The objective of this study is to construct a digital gene expression tag profile to identify genes potentially related to immune response in the ostrich. Exposure to boron leads to an immune response in the ostrich, although the underlying mechanism remains obscure. Thus, a dire need of biological resource in the form of transcriptomic data for ostriches arises to key out genes and to gain insights into the function of boron on the immune response of thymus. For this purpose, RNA-Seq analysis was performed using the Illumina technique to investigate differentially expressed genes in ostrich thymuses treated with different boric acid concentrations (0, 80, and 640 mg/L). Compared with the control group, we identified 309 upregulated and 593 downregulated genes in the 80 mg/L treated sample and 228 upregulated and 1816 downregulated genes in 640 mg/L treated sample, respectively. Trend analysis of these differentially expressed genes uncovers three statistically significant trends. Functional annotation analysis of the differentially expressed genes verifies multiple functions associated with immune response. When ostrich thymuses were treated with boron, expression changes were observed in genes predominantly associated with MAPK and calcium signaling pathways. The results of this study provide all-inclusive information on gene expression at the transcriptional level that further enhances our apprehension for the molecular mechanisms of boron on the ostrich immune system. The calcium and MAPK signaling pathways might play a pivotal role in regulating the immune response of boron-treated ostriches.

Development of polymorphic EST-SSR markers and characterization of the autotetraploid genome of sainfoin (Onobrychis viciifolia)

Background

Sainfoin (Onobrychis viciifolia) is a highly nutritious, tannin-containing, and tetraploid forage legume. Due to the lack of detailed transcriptomic and genomic information on this species, genetic and breeding projects for sainfoin improvement have been significantly hindered.

Methods

In this study, a total of 24,630,711 clean reads were generated from 14 different sainfoin tissues using Illumina paired-end sequencing technology and deposited in the NCBI SRA database (SRX3763386). From these clean reads, 77,764 unigene sequences were obtained and 6,752 EST-SSRs were identified using de novo assembly. A total of 2,469 primer pairs were designed, and 200 primer pairs were randomly selected to analyze the polymorphism in five sainfoin wild accessions.

Results

Further analysis of 40 sainfoin individuals from the five wild populations using 61 EST-SSR loci showed that the number of alleles per locus ranged from 4 to 15, and the expected heterozygosity varied from 0.55 to 0.91. Additionally, by counting the EST-SSR band number and sequencing the three or four bands in one sainfoin individual, sainfoin was confirmed to be autotetraploid. This finding provides a high level of information about this plant.

Discussion

Through this study, 61 EST-SSR markers were successfully developed and shown to be useful for genetic studies and investigations of population genetic structures and variabilities among different sainfoin accessions.

Comparative de novo flower transcriptome analysis of polygamodioecious tree Garcinia indica

To extend our understanding of molecular mechanism of sex determination in agro-economically important, polygamodioecious tree Garcinia indica (Kokum), high-throughput, next-generation flower transcriptome sequencing (NGS), and comparative analyses were performed to investigate differentially expressed gene in bisexual, female, and male flowers. A total of 49414 unigenes in BS, 45944 unigenes in FL, and 49028 unigenes in ML flowers were annotated. KO annotations revealed that 25 functional categories were large number of genes which were annotated to 'signal transduction'. We identified 33 genes for 'auxin response factor' and 50 for 'ethylene-responsive factor' whose expression changed significantly in all the three paired library combinations. Furthermore, key regulators of floral development such as FLC, SVP, AP1, AP2, AP3, AG, AGL2, AGL4, AGL9, and PI were identified. A total of 327 differentially expressed MADS-box genes were identified in G. indica transcriptome. Analysis of MADS-box genes identified five genes such as MADS AGL11, CRS2-associated factor chloroplastic, conserved hypothetical protein, uncharacterized protein LOC104422218, and MADS-box JOINTLESS-like isoform X3 significantly expressed in only FL flower. In addition, number of DEGs like dynamin 2A, auxin response factor, and spermidine synthase involved in sex expression and reproduction were discovered. The expression patterns of selected genes matched well with the expression levels of unigenes by transcriptome sequencing. Our large-scale comparative analyses may provide valuable hints for the next insights into the molecular mechanism of sex determination in G. indica.

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.

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.

Transcriptome and proteome analyses of the molecular mechanisms associated with coix seed nutritional quality in the process of breeding

In present study, the content of main nutrients in wild coix seed were significantly higher than cultivated coix seed. Transcriptome and proteome were combined to provide new insight of the molecular mechanisms linked to nutritional quality of wild coix seed and cultivated coix seed by RNA sequencing and isobaric tags for relative and absolute quantification techniques. A total of 20,039 genes and 471 proteins exhibited differential expression level in cultivated coix seed when compared with wild coix seed. These genes and proteins revealed that the pathway of flavonoids biosynthesis, starch and sucrose metabolism, lipid metabolism and amino acid metabolism were linked to nutritional quality of coix seed. Our results provided transcriptomics and proteomics information with respect to the molecular mechanisms of nutritional changes of coix seed, identified key genes and proteins that associated with the metabolism and accumulation of nutrients, and helped clarify the mechanisms of nutrient differences.

Comparative transcriptome analysis provides global insight into gene expression differences between two orchid cultivars

The orchids GL and YL are two cultivars of Cymbidium longibracteatum. YL displays an obviously yellowing rhizome and yellow leaves, while GL ('Longchangsu') shows dark green leaves and greenish rhizome. But the molecular mechanism for the differences between the two cultivars is poorly understood. In the present study, we showed that the structure of chloroplasts was significantly damaged in YL. Biochemical analysis uncovered the contents of chlorophyll a, chlorophyll b, total chlorophyll and carotenoid were notably decreased in YL. Using RNA-Seq technology, more than 38 million clean reads were generated in each pool, and 116,422 unigenes were assembled de novo. 6,660 unigenes with differential expression patterns (FDR≤0.01 and |log2 ratio|≥1) were totally identified between the two cultivars. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of differentially expressed unigenes (DEGs) suggested 33 KEGG pathways were notably enriched, including biological processes such as “phenylpropanoid biosynthesis”, “phagosome”, “starch and sucrose metabolism”, “drug metabolism—cytochrome P450”, “fatty acid elongation”, and “flavone and flavonol biosynthesis”. Further analysis revealed that chlorophyll degeneration related unigene (c48794_g1) and flavonoid biosynthesis related unigenes (c16388_g1, c48963_g1, c63571_g1, c4492_g1, c52282_g1, c78740_g1, c4645_g1) were up-regulated while carotenoid biosynthesis related unigene (c7212_g1) were down-regulated in YL. Additionally, six of NAC, R2R3-MYB, bHLH transcription factors (c42861_g1, c105949_g1, c61265_g1, c42659_g1, c82171_g1, c19158_g1) might be involved in regulation of pigment biosynthesis. The chlorophyll degeneration and the flavonoid biosynthesis related unigenes up-regulation together with the carotenoid biosynthesis related unigenes down-regulation may contribute to the yellowing phenotype of YL.

Development of genic-SSR markers and genetic diversity of Indian lettuce (Lactuca indica L.) in South Korea

Indian lettuce (Lactuca indica L.) is an undomesticated wild plant with high economic potential. We know little about the plant's genome, such as its DNA markers, making genetic research using this plant difficult. In this study, 100 genic simple sequence repeat (SSR) primers with a 99-250 bp target amplicon were synthesized from L. indica transcriptomic sequences. These primers were examined in 8 diverse L. indica accessions, and 90 polymorphic SSRs were obtained. Twenty-three of the 90 polymorphic SSRs were used to investigate transferability to another two Lactuca species, Lactuca serriola and Lactuca sativa. Genetic diversity was investigated in 77 Lactuca accessions, including 73 L. indica collected from across South Korea, 2 L. serriola, and 2 L. sativa. Our genic-SSR markers were highly polymorphic with a mean polymorphic information content of 0.61 and, on average, 10.83 alleles per locus. The average expected heterozygosity (0.76) was higher than the observed heterozygosity. An analysis of molecular variance revealed that most of the total variance in our population is attributable to genetic variation among accessions, rather than among provinces. STRUCTURE, unweighted neighbor-joining phylogenetic trees, and principal coordinate analyses resulted in three clusters, where northern and central-southern L. indica accessions were grouped into two clusters with some admixture. The L. serriola and L. sativa accessions did not produce a separate cluster due to a small sample size. These results show our SSR markers will be useful in germplasm assessment and genetic studies of L. indica and other Lactuca species.

Digital gene expression profiling analysis and its application in the identification of genes associated with improved response to neoadjuvant chemotherapy in breast cancer

BACKGROUND

This study aimed to screen sensitive biomarkers for the efficacy evaluation of neoadjuvant chemotherapy in breast cancer.

METHODS

In this study, Illumina digital gene expression sequencing technology was applied and differentially expressed genes (DEGs) between patients presenting pathological complete response (pCR) and non-pathological complete response (NpCR) were identified. Further, gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were then performed. The genes in significant enriched pathways were finally quantified by quantitative real-time PCR (qRT-PCR) to confirm that they were differentially expressed. Additionally, GSE23988 from Gene Expression Omnibus database was used as the validation dataset to confirm the DEGs.

RESULTS

After removing the low-quality reads, 715 DEGs were finally detected. After mapping to KEGG pathways, 10 DEGs belonging to the ubiquitin proteasome pathway (HECTD3, PSMB10, UBD, UBE2C, and UBE2S) and cytokine-cytokine receptor interactions (CCL2, CCR1, CXCL10, CXCL11, and IL2RG) were selected for further analysis. These 10 genes were finally quantified by qRT-PCR to confirm that they were differentially expressed (the log2 fold changes of selected genes were - 5.34, 7.81, 6.88, 5.74, 3.11, 19.58, 8.73, 8.88, 7.42, and 34.61 for HECTD3, PSMB10, UBD, UBE2C, UBE2S, CCL2, CCR1, CXCL10, CXCL11, and IL2RG, respectively). Moreover, 53 common genes were confirmed by the validation dataset, including downregulated UBE2C and UBE2S.

CONCLUSION

Our results suggested that these 10 genes belonging to these two pathways might be useful as sensitive biomarkers for the efficacy evaluation of neoadjuvant chemotherapy in breast cancer.

Transcriptome de novo assembly and analysis of differentially expressed genes related to cytoplasmic male sterility in onion

Onion (Allium cepa L.) is one of the major vegetables in China and accounts for a large proportion of China's vegetable exports. Onion cytoplasmic male sterility, which is often used in onion breeding, is caused by the interaction between the nuclear genes and the cytoplasm. However, the underlying molecular mechanism of onion cytoplasmic male sterility remains unclear. In this study, we analysed the anther microstructure of the onion cytoplasmic male sterile line SA2 and the onion maintainer line SB2. We found that the pollen abortion in SA2 occurred at the tetrad stage during the microspore development, which was very different from that in SB2. We used the Illumina HiSeq platform to sequence RNA from anthers at the tetrad stage collected from the SA2 and SB2 lines. The RNA sequencing and transcriptome assembly produced 146,413 All-Unigenes. Based on an analysis of the differentially expressed genes, we identified two cytoplasmic control genes, atp9 and cox1, and three nuclear-related genes, SERK1, AG and AMS. These transcriptomic results were also verified by fluorescence quantitative PCR. Our study provides important information about genes related to onion cytoplasmic male sterility, and it will help improve the understanding of the molecular mechanism of onion cytoplasmic male sterility.

RNA sequencing on Amomum villosum Lour. induced by MeJA identifies the genes of WRKY and terpene synthases involved in terpene biosynthesis

Amomum villosum Lour. is an important Chinese medicinal plant that has diverse medicinal functions, and mainly contains volatile terpenes. This study aims to explore the WRKY transcription factors (TFs) and terpene synthase (TPS) unigenes that might be involved in terpene biosynthesis in A. villosum, and thus providing some new information on the regulation of terpenes in plants. RNA sequencing of A. villosum induced by methyl jasmonate (MeJA) revealed that the WRKY family was the second largest TF family in the transcriptome. Thirty-six complete WRKY domain sequences were expressed in response to MeJA. Further, six WRKY unigenes were highly correlated with eight deduced TPS unigenes. Ultimately, we combined the terpene abundance with the expression of candidate WRKY TFs and TPS unigenes to presume a possible model wherein AvWRKY61, AvWRKY28, and AvWRKY40 might coordinately trans-activate the AvNeoD promoter. We propose an approach to further investigate TF unigenes that might be involved in terpenoid biosynthesis, and identified four unigenes for further analyses.

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.

Comprehensive transcriptome profiling of Salvia miltiorrhiza for discovery of genes associated with the biosynthesis of tanshinones and phenolic acids

Tanshinones and phenolic acids are crucial bioactive compounds biosynthesized in Salvia miltiorrhiza. Methyl jasmonate (MeJA) is an effective elicitor to enhance the production of phenolic acids and tanshinones simultaneously, while yeast extract (YE) is used as a biotic elicitor that only induce tanshinones accumulation. However, little was known about the different molecular mechanism. To identify the downstream and regulatory genes involved in tanshinone and phenolic acid biosynthesis, we conducted comprehensive transcriptome profiling of S. miltiorrhiza hairy roots treated with either MeJA or YE. Total 55588 unigenes were assembled from about 1.72 billion clean reads, of which 42458 unigenes (76.4%) were successfully annotated. The expression patterns of 19 selected genes in the significantly upregulated unigenes were verified by quantitative real-time PCR. The candidate downstream genes and other cytochrome P450s involved in the late steps of tanshinone and phenolic acid biosynthesis pathways were screened from the RNA-seq dataset based on co-expression pattern analysis with specific biosynthetic genes. Additionally, 375 transcription factors were identified to exhibit a significant up-regulated expression pattern in response to induction. This study can provide us a valuable gene resource for elucidating the molecular mechanism of tanshinones and phenolic acids biosynthesis in hairy roots of S. miltiorrhiza.

Succinyl-proteome profiling of Dendrobium officinale, an important traditional Chinese orchid herb, revealed involvement of succinylation in the glycolysis pathway

BACKGROUND

Lysine succinylation is a ubiquitous and important protein post-translational modification in various eukaryotic and prokaryotic cells. However, its functions in Dendrobium officinale, an important traditional Chinese orchid herb with high polysaccharide contents, are largely unknown.

RESULTS

In our study, LC-MS/MS was used to identify the peptides that were enriched by immune-purification with a high-efficiency succinyl-lysine antibody. In total, 314 lysine succinylation sites in 207 proteins were identified. A gene ontology analysis showed that these proteins are associated with a wide range of cellular functions, from metabolic processes to stimuli responses. Moreover, two types of conserved succinylation motifs, '***Ksuc******K**' and '****EKsuc***', were identified. Our data showed that lysine succinylation occurred on five key enzymes in the glycolysis pathway. The numbers of average succinylation sites on these five enzymes in plants were lower than those in bacteria and mammals. Interestingly, two active site amino acids residues, K103 and K225, could be succinylated in fructose-bisphosphate aldolase, indicating a potential function of lysine succinylation in the regulation of glycolytic enzyme activities. Furthermore, the protein-protein interaction network for the succinylated proteins showed that several functional terms, such as glycolysis, TCA cycle, oxidative phosphorylation and ribosome, are consisted.

CONCLUSIONS

Our results provide the first comprehensive view of the succinylome of D. officinale and may accelerate future biological investigations of succinylation in the synthesis of polysaccharides, which are major active ingredients.

Genetic mapping of the LOBED LEAF 1 (ClLL1) gene to a 127.6-kb region in watermelon (Citrullus lanatus L.)

The lobed leaf character is a unique morphologic trait in crops, featuring many potential advantages for agricultural productivity. Although the majority of watermelon varieties feature lobed leaves, the genetic factors responsible for lobed leaf formation remain elusive. The F2:3 leaf shape segregating population offers the opportunity to study the underlying mechanism of lobed leaf formation in watermelon. Genetic analysis revealed that a single dominant allele (designated ClLL1) controlled the lobed leaf trait. A large-sized F3:4 population derived from F2:3 individuals was used to map ClLL1. A total of 5,966 reliable SNPs and indels were identified genome-wide via a combination of BSA and RNA-seq. Using the validated SNP and indel markers, the location of ClLL1 was narrowed down to a 127.6-kb region between markers W08314 and W07061, containing 23 putative ORFs. Expression analysis via qRT-PCR revealed differential expression patterns (fold-changes above 2-fold or below 0.5-fold) of three ORFs (ORF3, ORF11, and ORF18) between lobed and non-lobed leaf plants. Based on gene annotation and expression analysis, ORF18 (encoding an uncharacterized protein) and ORF22 (encoding a homeobox-leucine zipper-like protein) were considered as most likely candidate genes. Furthermore, sequence analysis revealed no polymorphisms in cDNA sequences of ORF18; however, two notable deletions were identified in ORF22. This study is the first report to map a leaf shape gene in watermelon and will facilitate cloning and functional characterization of ClLL1 in future studies.

Norovirus contamination and the glycosphingolipid biosynthesis pathway in Pacific oyster: A transcriptomics study

Noroviruses are the primary pathogens associated with shellfish-borne gastroenteritis outbreaks. These viruses remain stable in oysters, suggesting an active mechanism for virus concentration. In this study, a deep RNA sequencing technique was used to analyze the transcriptome profiles of Pacific oysters at different time points after inoculation with norovirus (GII.4). We obtained a maximum of 65, 294, 698 clean sample reads. When aligned to the reference genome, the average mapping ratio of clean data was approximately 65%. In the samples harvested at 12, 24, and 48 h after contamination, 2,223, 2,990, and 2020 genes, respectively, were differentially expressed in contaminated and non-contaminated oyster digestive tissues, including 500, 1748, and 1039 up-regulated and 1723, 1242, and 981 down-regulated genes, respectively. In particular, FUT2 and B3GNT4, genes encoding the signaling components of glycosphingolipid biosynthesis, were significantly up-regulated in contaminated samples. In addition, we found up-regulation of some immune- and disease-related genes in the MHC I pathway (PA28, HSP 70, HSP90, CANX, BRp57, and CALR) and MHC II pathway (GILT, CTSBLS, RFX, and NFY), although NoVs did not cause diseases in the oysters. We detected two types of HBGA-like molecules with positive-to-negative ratios similar to type A and H1 HBGA-like molecules in digestive tissues that were significantly higher in norovirus-contaminated than in non-contaminated oysters. Thus, our transcriptome data analysis indicated that a human pathogen (GII.4 Norovirus) was likely concentrated in the digestive tissues of oysters via HBGA-like molecules that were synthesized by the glycosphingolipid biosynthesis pathway. The identified differentially expressed genes also provide potential candidates for functional analysis to identify genes involved in the accumulation of noroviruses in oysters.

Transcriptome Profiling of Tomato Uncovers an Involvement of Cytochrome P450s and Peroxidases in Stigma Color Formation

Stigma is a crucial structure of female reproductive organ in plants. Stigma color is usually regarded as an important trait in variety identification in some species, but the molecular mechanism of stigma color formation remains elusive. Here, we characterized a tomato mutant, yellow stigma (ys), that shows yellow rather than typical green color in the stigma. Analysis of pigment contents revealed that the level of flavonoid naringenin chalcone was increased in the ys stigma, possibly as a result of higher accumulation of p-coumaric acid, suggesting that naringenin chalcone might play a vital role in yellow color control in tomato stigma. To understand the genes and gene networks that regulate tomato stigma color, RNA-sequencing (RNA-Seq) analyses were performed to compare the transcriptomes of stigmas between ys mutant and wild-type (WT). We obtained 507 differentially expressed genes, in which, 84 and 423 genes were significantly up-regulated and down-regulated in the ys mutant, respectively. Two cytochrome P450 genes, SlC3H1 and SlC3H2 which encode p-coumarate 3-hydroxylases, and six peroxidase genes were identified to be dramatically inhibited in the yellow stigma. Further bioinformatic and biochemical analyses implied that the repression of the two SlC3Hs and six PODs may indirectly lead to higher naringenin chalcone level through inhibiting lignin biosynthesis, thereby contributing to yellow coloration in tomato stigma. Thus, our data suggest that two SlC3Hs and six PODs are involved in yellow stigma formation. This study provides valuable information for dissecting the molecular mechanism of stigma color control in tomato. Statement: This study reveals that two cytochrome P450s (SlC3H1 and SlC3H2) and six peroxidases potentially regulate the yellow stigma formation by indirectly enhancing biosynthesis of yellow-colored naringenin chalcone in the stigma of tomato.