Next Article in Journal
Ligand and Structure-Based Approaches for the Identification of Peptide Deformylase Inhibitors as Antibacterial Drugs
Previous Article in Journal
Genome-Wide Analysis of Codon Usage Bias in Epichloë festucae
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 17
Issue 7
10.3390/ijms17071139
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Transcriptomic Analysis Reveals the Mechanism of Phelipanche aegyptiaca Seed Germination

by
Zhaoqun Yao
,
Fang Tian
,
Xiaolei Cao
,
Ying Xu
,
Meixiu Chen
,
Benchun Xiang
and
Sifeng Zhao
*
Key Laboratory at Universities of Xinjiang Uygur Autonomous Region for Oasis Agricultural Pest Management and Plant Protection Resource Utilization, Shihezi University, Shihezi 832003, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(7), 1139; https://doi.org/10.3390/ijms17071139
Submission received: 10 May 2016 / Revised: 10 July 2016 / Accepted: 11 July 2016 / Published: 15 July 2016
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Phelipanche aegyptiaca is one of the most destructive root parasitic plants of Orobanchaceae. This plant has significant impacts on crop yields worldwide. Conditioned and host root stimulants, in particular, strigolactones, are needed for unique seed germination. However, no extensive study on this phenomenon has been conducted because of insufficient genomic information. Deep RNA sequencing, including de novo assembly and functional annotation was performed on P. aegyptiaca germinating seeds. The assembled transcriptome was used to analyze transcriptional dynamics during seed germination. Key gene categories involved were identified. A total of 274,964 transcripts were determined, and 53,921 unigenes were annotated according to the NR, GO, COG, KOG, and KEGG databases. Overall, 5324 differentially expressed genes among dormant, conditioned, and GR24-treated seeds were identified. GO and KEGG enrichment analyses demonstrated numerous DEGs related to DNA, RNA, and protein repair and biosynthesis, as well as carbohydrate and energy metabolism. Moreover, ABA and ethylene were found to play important roles in this process. GR24 application resulted in dramatic changes in ABA and ethylene-associated genes. Fluridone, a carotenoid biosynthesis inhibitor, alone could induce P. aegyptiaca seed germination. In addition, conditioning was probably not the indispensable stage for P. aegyptiaca, because the transcript level variation of MAX2 and KAI2 genes (relate to strigolactone signaling) was not up-regulated by conditioning treatment.

Graphical Abstract

1. Introduction

A total of 4100 plant species in approximately 19 families of flowering plants have evolved to parasitize other plants. These plants are widely distributed in Africa, Europe, Asia, America, and Australia. Haustoria are used by these plants to obtain water and nutrients from other living plants or hosts [1,2]. Parasitic weeds cause serious agricultural problems, resulting in remarkable crop losses in many regions worldwide. About $1.3 to 2.6 billion of annual food crop losses are caused by broomrape infestation [3]. Phelipanche spp. are holoparasites that lack chlorophyll. They parasitize more-temperate climate crops, such as sunflower, tomato, potato, tobacco, carrot, clovers, cucumber, rapeseed and legumes [4]. Orobanche cumana has damaged 20,000 ha of farmland in Greece and China, with estimated yield losses of 60% in Greece and 20%–50% in China. In 1994, Phelipanche aegyptiaca had extensive infestations of muskmelon and watermelon, leading to 20%–70% yield losses in Xinjiang Province, China [5]. The life cycle of Phelipanche spp. has a number of mechanisms that coordinate the life cycles of parasites to that of their host. The main steps in the life cycle are conditioning of seeds, germination under stimulants secreted by hosts, adhesion and formation of appressorium, penetration through host tissues, formation of haustorium to connect the host vascular tissues, development of a tubercle and apex, stem growth and emergence, and flowering and seed production [6,7,8]. The seeds of Orobanche spp. contain only little reserves. These seeds can survive for a few days only after germination unless they reach a host root to establish a xylem connection.
The Orobanche spp. parasitic strategy generally succeeds by coordinating early developmental stages with chemical signals from hosts. An important step in the life cycle of Orobanche spp. is their germination at the right place and time, enabling them to establish the connection they require to survive. Orobanche spp. usually use so-called germination stimulants secreted by roots of their hosts. To date, three different types of compounds, namely, dihydroquinones (dihydrosorgoleone), sesquiterpene lactones, and strigolactones (SLs), have been identified as chemical signals or germination stimulants for Striga spp. and Orobanche spp. Among these germination stimulants, SLs are the most active in inducing germination at 10−7 to 10−15 mol/L [9,10]. SLs are new plant hormones that control shoot branching, root architecture, cambial growth, and senescence [11,12]. SLs are synthesized from carlactone, which is derived from all-trans β-carotene via the action of an isomerase (D27) and two carotenoid cleavage dioxygenases (CCD7 and CCD8). Then, further ring closures and functionalizations involves in members of the CYP711 family (MAX1). Once synthesized, SLs may be transported by PhPDR1, a member of the ABC family within the plant and in the rhizosphere. Finally, MAX2 interacts with D14/KAI2 in an SLs-dependent manner, and this leads to SL ubiquitylation dependent degradation of D53 by the SCFMAX2 complex [10,11]. Further, these hormones also serve as extra organismal signals in soil that recruit symbioses with arbuscular mycorrhizal fungi and trigger Orobanchaceae plant seed germination [13].
Orobanche spp. exerts the greatest damage prior to their emergence, and the majority of field loss may occur before diagnosis of infection. Numerous physical, cultural, chemical, and biological approaches have been explored against root parasites. However, none of these techniques are effective and provide economical results [14,15,16,17,18]. SLs are regarded as potential new strategies to control Orobanche spp. [19]. The tomato SL-deficient mutant (SL-ORT1) displays resistance to P. aegyptiaca because of the inability of SL-ORT1 roots to produce and secrete natural germination stimulants (SLs) to the rhizosphere. Silencing of the tomato MAX3/CCD7 gene, which is the critical gene for SL production, reduces the number of P. ramose infection [20,21]. AM symbiosis in tomato also reduces SL production and P. ramosa infection [22]. In 2016, it was successful to reduce O. ramosa in tobacco fields by using SL analogues via suicidal germination approach [23]. However, genomic and molecular resources for Phelipanche are limited [24], and the mechanism of SLs inducing Phelipanche seed germination remains unclear.
Transcriptome and proteome technology can facilitate the understanding of the molecular basis of complex developmental processes [25]. De novo assembly and characterization of the transcriptome of Cuscuta pentagona and three parasites of Orobanchaceae have uncovered genes associated with plant parasitism [26,27,28]. Transcriptome sequencing successfully provided new insight into seed germination processes in Arabidopsis, wheat, and cotton [29,30,31]. In the present study, dormant, conditioned, and GR24 (a synthetic analog of SL)-treated P. aegyptiaca seeds were used to assemble and annotate a reference transcriptome. The transcriptome data were then used to analyze different expressions of genes during different seed germination stages. Finally, the role of plant hormones in P. aegyptiaca seed germination was investigated by physiological tests using differentially expressed genes (DEGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis.

2. Results

2.1. De Novo Assembly of P. aegyptiaca Transcriptome

RNA-sequencing (RNA-seq) library was prepared from dormant, conditioned, and GR24-treated seeds and sequenced using the Hiseq 2500 platform. A total of 78,540,698 reads were obtained. The percentage of Q30 base in all samples was over 88.36%. Illumina adapters were trimmed, and low-quality bases were filtered. Trinity software package [32] was used to assemble clean data. A total of 274,964 transcripts (longer than 200 bp) were found. The average length of N50 of transcripts was 1066 bp, and the evaluated number of unigenes was 94,419 (Table 1 and Figure 1).

2.2. Transcriptome Functional Annotation

Various approaches were used to acquire comprehensive information on the assembled transcriptomes. First, unigenes were used for homology searches against the non-redundant (NR), Gene Ontology (GO), Cluster of Orthologous Groups of proteins (COG), KOG, and KEGG databases using the BLASTX algorithm, with an E-value ≤ 10−5. Subsequently, amino acid sequences from predicted unigenes were aligned against the Pfam database with an E-value ≤ 10−6 using HMMER software. Thus, 53,921 unigenes were annotated. However, numerous unigenes (40,498) remained unannotated. Annotation statistics of BLASTX hits and domain hits are summarized in Table S1.
BLASTX against the NR database provided insight into the taxonomic distribution of the transcripts. High mapping rates of contigs to the protein databases suggested that most of the contigs could be translated into proteins. BLASTX hits and top hits in terms of the total numbers of hits of all transcripts were mostly observed in Sesamum indicum (11,865 hits), Mimulus guttatus (3410 hits), and Vitis vinifera (2261 hits) (Figure 2).
GO is an international classification system for standardized vocabulary, which is used in the comprehensive description of functions of uncharacterized genes. The three main independent GO categories are biological processes, molecular functions, and cellular components. Our annotated unigenes were assigned to GO terms for functional classification. A total of 21,851 (27.3%), 21,961 (27.5%), and 23,287 (29.2%) unigenes were annotated in the biological process, molecular function, and cellular component categories, respectively (Figure 3). The majority of unigenes under cellular component were involved in cells (23.42%), cell parts (23.68%), and organelles (17.05%). Approximately 43.14% and 41.33% of unigenes within the molecular function category were clustered in catalytic activity and binding, respectively. Sub-categories metabolic (28.19%), cellular (23.76%), and single-organism (17.93%) processes were the major sub-categories of biological process.

2.3. qRT-PCR Validation

We also used RNA samples isolated for RNA sequencing to perform qRT-PCR analysis (Figures S1 and S2). The 12 randomly chosen genes and 8 genes relating to plant hormone (GA, ABA, ethylene and BR) biosynthesis were used to qRT-PCR validation. From the supplemental Figure S2, the 12 gene dynamic changes of qPCR and the RNA-seq results were consistent. Furthermore, expression levels of the 8 plant hormone genes are also supported by the changes of each hormone content change (Figure 4 and Figure 5). These results suggested the applicability of RNA-Seq to Phelipanche transcriptome analysis is an accurate and reliable way to find DEGs during P. aegyptiaca seed germination.

2.4. Comparative Analysis of Differential Expression during Germination

We used RPKM (reads aligned to gene per kilobase of exon per million mapped reads) to estimate gene expression values to better understand the dynamic performance among the three different transcriptomes during seed germination. False discovery rate (FDR) <0.01 and fold change ≥2 were used as cutoff values to identify differentially regulated transcripts. A total of 5324 DEGs among dormant, conditioned, and GR24-treatedseeds were identified. A total of 633 and 61 unigenes were up-regulated and down-regulated in dormant seeds, respectively, compared with conditioned seeds. A total of 3764 unigenes were up-regulated and 866 unigenes were down-regulated between conditioned and GR24-treated seeds (Figure S3).
Based on the DEGs among the three samples, clustering analysis was performed via MeV edition 4.90 using the algorithm of “hierarchical clustering” (Figure 6). Groups of unigenes with similar expression patterns under different stages were measured in nine clusters. Gradual changes in the gene expression patterns of clusters 1, 2, and 3 varied, and nearly half of DEGs belonged to these three clusters. DEGs of cluster 4 were not active until GR24 application. These genes were probably regulated by GR24, and play crucial roles in the last germination stage. Cluster 5 contained 1840 DEGs that were up-regulated during seed germination. These genes were evidently not only involved in plant hormones but also related to many metabolic pathways. However, the roles of these genes during P. aegyptiaca germination have not been elucidated. DEGs in cluster 6 were examined in dormant seeds, so these genes inhibited P. aegyptiaca seed germination, contrary to the expression patterns of genes of clusters 7 and 8. Thus, cluster 7 genes probably inhibited germination, in contrast to the promotion of germination by cluster 8 genes. The inhibition or promotion of the genes in clusters 7 and 8 was released or induced by conditioning, respectively. DEGs in the last group, cluster 9, were down-regulated gradually in chronological order. These genes played more important roles in dry seed than during seed germination.
Furthermore, to determine the biological significance of the DEGs, GO terms enrichment analysis of the total DEGs was also carried out (p ≤ 0.05). As shown in Figure S4, the biological processes of translation (GO:0006412), response to salt stress (GO:0009651), oxidation-reduction process (GO:0055114) were involved in 221, 44 and 43 DEGs, respectively, indicating mRNA translation, stress resistance and energy metabolism played a vital role in conditioned phase. Between the conditioned and GR24 treatment stage, 151, 119, 47 and 43 DEGs were enriched in the oxidation-reduction process (GO:0055114), translation (GO:0006412), metabolic process (GO:0008152), and protein phosphorylation (GO:0006468), respectively. This suggested energy metabolism and mRNA translation still play dominate roles, and protein and other metabolic process were carried out.

2.5. Gene Change during P. aegyptiaca Seed Germination

Removal of the seed coat permeability barrier, leaching of germination inhibitors, facilitating the access of stimulants to their putative cellular target, biosynthesis of plant hormones, and formation or activation of receptors for germination stimulants have been inferred to occur during conditioning [33,34]. DEGs between dormant and conditioned seeds showed the two biggest group of genes involved from DNA to protein repair and biosynthesis, as well as carbohydrate and energy metabolism (Figure 5). Tested P. aegyptiaca seeds and endogenous germination inhibitors have not been proven to inhibit the germination of this plant. We found that the gibberellin 20 oxidase 1 (GA20OX1) transcript was up-regulated and GA3 levels increased during conditioning. In addition, the 3-epi-6-deoxocathasterone 23-monooxygenase (CYP90C1) transcript was up-regulated. This cytochrome P450 enzyme plays an important role in early C-22 oxidation of brassinosteroid (BR) biosynthesis. The expression levels of the two genes were supported by quantitative real-time polymerase chain reaction (qRT-PCR) (Figure 6).
The second phase started after stimulation, during which sharp transition in Phelipanche seeds occurred [35]. Dominant DEGs were associated with DNA, RNA, protein, carbohydrate, and energy metabolism. More genes and higher expression levels were found during this phase. Additionally, genes involved in plant hormone biosynthesis and metabolism were dramatically up-regulated. Gibberellin 2-beta-dioxygenase 1 (GA2OX1) gene and gibberellin 3-beta-dioxygenase 1 (GA3OX1) transcript were up-regulated significantly, and GA3 levels increased. Cytochrome P450 90B1 (CYP90B1) gene, which encodes a BR biosynthesis enzyme, was activated dramatically, and the amount of BRs increased. The tendencies of the two phytohormones levels were apparently similar to those during the conditioning stage. Numerous studies have proven that ethylene can promote seed germination. Our data showed that in methionine adenosyltransferase (MAT), aminocyclopropanecarboxylate synthase (ACS), and aminocyclopropanecarboxylate oxidase (ACO), which are three key enzymes of ethylene biosynthesis, corresponding genes were up-regulated significantly after GR24 treatment (Figure 6). The 9-cis-epoxycarotenoid dioxygenase transcript, which plays a crucial role in abscisic acid (ABA) biosynthesis, was dramatically down-regulated. In addition, in ABA 8′-hydroxylase 1 (CYP707A1), which catalyzes hydroxylation at the C-8′ position of ABA, the transcript was up-regulated. These results coincided with the qRT-PCR results (Figure 6) and ABA level reduction.
As Kusumoto reported, receptors for germination stimulants are perhaps formed or activated during conditioning [34]. Transcriptome results showed that the F-box protein MAX2 expression level remained high in the whole germination process, which is associate with perception of SLs [36]. Furhter, two α/β-hydrolase KAI2-like genes were identified in our results, the expression of one was high, and the other expression level even down-regulated during the germination process (Figure S5). These results demonstrated that the SL receptor complex was not only active during conditioning or after GR24 application, but it was more likely active all the time.

2.6. Endogenous Hormone Levels and Germination Test

Changes in the endogenous plant hormone levels were tested to elucidate the roles of hormones. GA3, BRs, and ethylene levels continuously increased, similar to those in previous reports on seed germination. Interestingly, GA3 and BR levels gradually increased, contrary to the increase in ethylene level after GR24 application. An initial slight decrease in the ABA level was followed by a significant reduction after GR24 treatment (Figure 7).
A germination experiment was designed to explore the key changes in hormone levels for P. aegyptiaca seed germination. The results showed that individual application of the four hormones failed to induce P. aegyptiaca seed germination. However, fluridone, the inhibitor of carotenoid biosynthesis, promoted seed germination without GR24 application. GA3, BRs, and ethephon (it can be converted into ethylene at pH < 7) significantly help fluridone increase germination rates. GA3, fluridone, or GA3, as well as ethephon and fluridone, induced seed germination without GR24 application (Table 2 and Figure S6). These results were probably due to the enhanced GA/ABA ratio caused by exogenous fluridone, as well as the GA3 and ethylene-induced changes in the GA/ABA ratio value by adjusting the ABA content. Therefore, the change in ABA level may be the key factor in P. aegyptiaca seed germination, which is controlled by GR24.

3. Discussion

Transcriptome sequencing was used to identify the role of GR24 in P. aegyptiaca seed germination. Three developmental stages from dormant to GR24-treated seeds were selected. Several key genes and pathways were determined. Protein biosynthesis and metabolism, as well as energy metabolism, were dominant in P. aegyptiaca seed germination. Additionally, many genes were involved in plant hormone biosynthesis and transduction.

3.1. Transcriptome Assembly and Annotation

RNA-seq reads of dormant, conditioned, and GR24-treated seeds were used to assemble a reliable transcriptome. However, only 61% of the filtered transcripts were annotated. This result was similar to many other de novo-assembled plant transcriptomes, such as Capsicum frutescens, Nicotiana benthamiana, and Cuscuta pentagona [27,37,38]. We found that 35% hits belonged to S. indicum, the genome of which was sequenced in 2014 [39]. Although the taxonomic position of Phelipanche and Striga is closer compared with Sesamum and Mimulus, belonging to Lamiales, the two top hits results are S. indicum and M. guttatus. This resulted from a lack of Striga genome data.

3.2. DEGs Associated with Broomrape Seed Germination

3.2.1. DNA, RNA, and Protein

Numerous stored mRNAs, such as for seed reserve synthesis, late embryogenesis abundant proteins, protein synthesis, and degradation in mature dry seeds, are necessary for germination [40,41,42]. Degradation of stored mRNAs is a prerequisite to germination [43,44]. In our study, Trf4/5, CONT3, and PABP1 transcripts related to RNA degradation were up-regulated in the conditioned phase. By contrast, CNOT3 and DDX6 transcripts related to RNA degradation were down-regulated after GR24 application. Stored proteins were used to restart cellular activity, followed by mobilization of the stored mRNA pool [42].
Hunt and Macovei reported the DNA repair and replication in Arabidopsis and Medicagotruncatula seed germination [45,46]. We observed that xerodermapigmentosum complementation group C (XPC) and the HR23B genes of the DNA repair XPC complex were up-regulated. Replication protein A, which is a DNA replication protein, and ligase Lig I genes were up-regulated after GR24 treatment. Poly (ADP-ribose) polymerases (PARPs) were implicated in DNA repair in Arabidopsis [45]. We found that PARP transcripts were very high in dormant and conditioned seeds, but their expression was dramatically down-regulated after GR24 treatment. Thus, DNA repair may have been completed during the conditioning stage.
De novo RNA synthesis started after DNA repair and replication. α-amanitin, which is a highly potent and specific inhibitor of DNA-dependent RNA polymerase II, strongly affected the speed and uniformity of Arabidopsis seed germination [47]. Our data also showed that the DNA-dependent RNA polymeraseII transcripts were up-regulated during the entire germination period. Pentatricopeptide repeat proteins genes were up-regulated, which also participate in RNA splicing, cleavage, editing, stability, and translation [48].
De novo protein synthesis during early seed germination was indispensable. As Rajjou demonstrated, the germination of Arabidopsis seeds is entirely blocked in the presence of cycloheximide, a translation inhibitor [47]. However, free amino acids in dry seeds are not sufficient for protein synthesis during germination [49]. Therefore, storage protein degradation is an optimal way to solve this problem. The protein L-isoaspartyl-O-methyltransferase (PIMT) repair pathway facilitates the active elimination of deleterious protein products during germination [50]. Although we did not observe PIMT transcript activity, the two main protein degradation KEGG pathways, namely, proteasome- and ubiquitin-mediated proteolysis, were activated during the conditioning stage. Furthermore, majority of transcripts in both pathways were down-regulated in the presence of GR24.

3.2.2. Energy Metabolism

Energy is supplied from cellular respiration, including glycolysis of carbohydrates and β-oxidation of fatty acid, even amino acids, during seed germination. All biochemical reaction products feed into the tricarboxylic acid (TCA) cycle and provide power for mitochondrial electron transport chain to produce ATP [51].
The two main pathways for carbohydrate metabolism are the oxidative pentose phosphate pathway (OPPP) and glycolysis (Figure S7). Glucose-6-phosphate dehydrogenase (G6PDH) plays a key role in OPPP, and its expression level was stable during germination. By contrast, two key enzymes in the glycolysis pathway, namely, 6-phosphofructokinase (PFK) and pyruvate kinase (PK) transcripts, were up-regulated during the conditioning stage. After treatment with GR24, α-amylase and PFK gene expression levels was down-regulated. This phenomenon suggested that carbohydrates were depleted after conditioned treatment.
Thus, fatty acid or protein acted as an energy supplement in the next germination phase (Figure S7). Many lipoxygenases and fatty acyl-CoA synthetase, acyl-CoA dehydrogenase, enoyl-CoA hydratase, β-hydroxyacyl-CoA dehydrogenase, and acetyl-CoA acetyltransferase transcripts involved in fatty acid β-oxidation were active. Isocitratelyase (ICL), which is a key enzyme in the glyoxylate cycle, and its transcript was active. TAG provides energy to the embryo before the initiation of photosynthesis [52]. Although our results demonstrated that two crucial enzymes, namely, glycerol-3-phosphate dehydrogenase (G3PDHc) and triacylglycerol (TAG) lipase genes, were not activated, glycerol kinase expression increased after GR24 application. Therefore, fatty acid was the energy donor processed through both the glyoxylate cycle and glycerol shunt pathway (Figure S7).
Aspartate and alanine (AspAT) and aminotransferases (AlaAT) genes were activated during imbibition and inferred to participate in respiratory pathways [51]. Our data illustrated that the corresponding genes of both enzymes were not activated during the conditioning phase. However, only the AspAT transcript was activated after GR24 application. The activity and transcript levels of glutamate decarboxylase (GAD), which is the entry enzyme of the aminobutyrate (GABA) shunt, increase in germinating seeds [43,53]. GAD transcript expression level in P. aegyptiaca was high at the conditioning stage and upon GR24 application. Oxaloacetate is another product of aspartate. This product can feed into the TCA cycle either by transforming to PEP under PEPCK or by conversion to malate by cytosolic malate dehydrogenase (MDHc) (Figure S7). Our result showed that the PEPCK gene was activated during the whole developmental process, but MDHc transcripts were down-regulated after GR24 treatment. Thus, further studies are necessary to determine the manner by which the two pathways of amino acids supplying energy cooperate with each other during P. aegyptiaca seed germination at different stages.
Pyruvate, as the energy metabolism hub, feeds acetyl-CoA to the TCA cycle through pyruvate dehydrogenase complex (PDC) and alcohol dehydrogenase (ADH) or lactate dehydrogenase (LDH) under aerobic the condition, contributing to the transformation of pyruvate into ethanol or lactate, respectively. PDC-related genes were activated after imbibition and GR24 application. The genes of the entire coded TCA cycle enzymes, including malate dehydrogenase in mitochondria, citrate synthase, aconitatehydratase, isocitrate dehydrogenase, oxoglutarate dehydrogenase complex, succinate-CoA ligase, succinate dehydrogenase (SDH), and fumaratehydratase transcripts were up-regulated (Figure S7). The activity of these genes also highly increased during the early germination stage in many plants [51]. However, GR24 treatment resulted in the down-regulation of SDH and fumaratehydratase transcripts. This phenomenon needs further investigation. Additionally, the fermentation-related transcript, pyruvate dehydrogenase transcript, was up-regulated, contrary to LDH.

3.2.3. Phytohormone Activity

ABA inhibits both embryo growth and endosperm rupture to establish seed dormancy and inhibit germination. Nevertheless, GAs counteracts ABA responses and promotes seed germination [54]. Several studies have shown that ABA levels decrease during P. aegyptiaca seed conditioning [55]. ABA biosynthesis inhibitors, namely, fluridone and norflurazon, can promote dormancy release and germination [56,57]. Application of exogenous ABA inhibits P. ramosa seed germination [57]. In this study, changes in NCED and CYP707A1 transcripts coincided with the change in ABA level (Figure 7) and results from previous studies [57,58,59]. A recent study on P. ramosa seeds demonstrated that the seed response to SL is controlled by ABA-independent DNA methylation [60]. Thus, exogenous ABA biosynthesis inhibitor fluridone was examined and found to induce P. aegyptiaca seed germination without GR24. Although a large number of previous studies showed SLs act in a positive role with respect to seed germination, as ABA and SLs share a biosynthesis pathway, so this phenomenon resulting from the effect of ABA or SLs needs further discussion. In many studies on Striga and Orobanche, fluridone alone was not able to induce their seed germination [34,55,61,62,63], but our results showed that application of fluridone without GR24 could induce seed germination. We speculate the different results occurred because the Phelipanche seed used in studies was different, and this exciting phenomenon needs more investigation.
Several studies provided evidence that GA synthesis occurs during conditioning [64]. In our experiment, the GA20OX1 transcript was up-regulated after conditioned treatment and GR24 induced GA2OX1 and GA3OX1 transcript up-regulation. Some studies showed that application of GA biosynthesis inhibitors (e.g., uniconazole) during seed conditioning inhibits the subsequent germination of Striga hermonthica, P. ramosa, P. aegyptiaca, and O. minor in response to GR24 [57,62,65]. BR application after conditioning increased the rate of seed germination induced by stimulants [56]. CYP90C1and CYP90B1 transcripts were up-regulated during different germination stages. Although BR level increased during the whole process, BRs were inferred to alter the balance of other endogenous phytohormones, such as indole-3-acetic acid (IAA), ABA, GAs, and ethylene rather than acting directly [66]. However, we found that exogenous GA3 or BR promoted P. aegyptiaca seed response to GR24, but individual application of these two hormones failed to induce seed germination.
Ethylene, which is another plant hormone, also plays an important role in seed germination by promoting the weakening of the micropylar endosperm [67,68]. This hormone can neutralize many of the negative functions of ABA during germination [69,70,71]. Ethylene, ethephon, or aminocyclopropanecarboxylate (ACC) stimulates seed germination in numerous species, including several parasitic plants, such as P. ramosa [72] and some Striga species [73,74]. The results showed that methionine adenosyltransferase (MAT), ACS, and ACO expression levels were induced by GR24. This result was supported by Zehhar’s study, in which ethylene synthesis was required for the induction of P. ramosa seed germination by GR24 [57]. However, we failed to quantify ethylene accurately; low endogenous ethylene was probably released into the medium. However, ethylene down-regulated ABA accumulation by both inhibiting its synthesis and promoting its inactivation, as well as negatively regulated ABA signaling [75]. Thus, indirect degradation by controlling ethylene requires further examination.
The ABA biosynthesis pathway was inhibited, and its degradation was promoted by GR24. This is consistent with the research on CYP707A in the seed germination of P. ramosa [58,60]. By contrast, the ethylene biosynthesis pathway was activated. Therefore, the bottleneck of P. aegyptiaca seed germination was endogenous ethylene and the ABA level, which were probably controlled by SLs.
SLs carry out their physiological functions via a specific reception system similar to other plant hormones. SLs bind to the D14/MAX2/SCF complex, and then lead to SMAX1, a D53 homolog, degradation, finally acting on seed germination, and this model is widely accepted [11]. However, recently, it was reported that the α/β-hydrolase KAI2 gene (D14-like, also known as HTL) plays a central role in the perception of SLs [36,76]; this gene was thought to be the key gene in reception of the strigolactone-analogous compound karrikin [11], originally found in forest-fire smoke. The transcriptome and qRT-PCR analysis of MAX2, D14, and two KAI2-like genes illustrated their expression levels are stable during seed germination. We only found two KAI2-like genes in our data; however, there are 11 KAI2 genes in S. hermonthica [76]. According to the results, we conservatively speculate the transcripts of partial components of the SL perception complex existed prior to conditioning.

4. Materials and Methods

4.1. Sample Preparation

Mature seeds of P. aegyptiaca were collected from Cucumis melo L. in the Hami area of the Xinjiang Uygur Autonomous Region, China, in 2012. Dormant, conditioned, and GR24-treated seeds were used in this study. All seeds were surface-sterilized by 75% ethanol for 1 min and 1% sodium hypochlorite for 20 min followed by three rinses in sterile distilled water [77]. Subsequently, 10 g of seeds were frozen in liquid nitrogen as the dormant seed sample and then stored at −80 °C until RNA was extracted. The conditioned seed sample was treated as follows: 10 g of sterilized seeds were spread on sterilized glass fiber filter paper moistened with sterile distilled water. The seeds were placed into 10 cm Petri dishes and incubated in the dark at 25 °C/20 °C for 7 days [60]. The seeds were then frozen in liquid nitrogen and stored at −80 °C until RNA was extracted. The GR24-treated seed sample consisted of 10 g of conditioned seeds induced by applying GR24 at 10−6 mol/L (0.1% acetone, v/v) for 1 day, frozen in liquid nitrogen, and stored at −80 °C until RNA was extracted [78].

4.2. RNA-seq Library Preparation and Sequencing

Total RNA of each sample was isolated using Trizol Reagent (Invitrogen, Gaithersburg, MD, USA) according to the manufacturer’s instructions. RNA quality was characterized by agarose gel electrophoresis and spectrophotometry. High-quality RNA, with 28S:18S of more than 1.5 and absorbance 260/280 ratio between 1.8 and 2.2, was used for library construction and sequencing. Illumina HiSeq2500 library was constructed according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). Magnetic beads with oligo (dT) (Thermo Scientific, Waltham, MA, USA) were used to enrich the mRNA from total RNA. mRNA was randomly cleaved by adding fragmentation buffer. The fragments were used to synthesize first-strand cDNA with random hexamers. After first-strand cDNA was synthesized, the buffer, dNTPs, RNase H, and DNA polymerase I were added to synthesize double-stranded cDNA (dscDNA), which was purified using AMPure XP beads. The purified dscDNA was end repaired, added to an A-tail, and linked with sequencing adapters [79]. AMPureXP beads (Beckman Coulter, CA, USA) were used to choose suitable fragments. The sequencing library was enriched with PCR amplification. After verification with Qubit 2.0 and quantification with Agilent 2100, the effective concentration was accurately quantified using qRT-PCR. The library was then sequenced with an Illumina HiSeq2500 platform (Biomarker Technology Co., Ltd., Beijing, China).

4.3. Preprocessing of Illumina Reads

The cDNA library of each sample was sequenced by Illumina HiSeq2500 based on sequencing by synthesis technology. This library produced amounts of high-quality reads (most bases reached or had excess Q30). Preprocessing of reads involved Q20-quality trimming (removal of low-quality reads with an average Phred quality score of 20, trimming of low-quality bases from the 3′ ends of the reads) and removal of adapter/primer contamination and duplicated reads using custom Perl scripts. Preprocessed reads were sorted into individual samples based on barcodes using fastx_barcode_splitter. Barcodes were trimmed using fastx_trimmer from Fastx_toolkit.
Prior to assembly, the raw reads were cleaned by removing adapter sequences using standard Illumina pipeline, namely, the CASSAVA program (http://support.illumina.com/sequencing/sequencing_software/casava.ilmn). Adaptor contamination was removed, and reads were screened from the 3′ to 5′ ends to trim the bases with a quality score of Q < 20 using 5 bp windows. Reads with final length < 25 bp were removed.

4.4. De Novo Transcriptome Assembly

The Trinity software package (http://trinityrnaseq.github.io/) was used for de novo assembly of P. aegyptiaca transcriptome from preprocessed RNA-seq reads [32]. A K-mer library was built using the sequencing reads, and K-mer, which probably contained errors, was deleted. The highest frequency of K-mer was used as seed to extend greedily to both ends. This process was continuously cycled until the K-mer ran out. The contigs were clustered to obtain components. A simplified De Bruijn diagram of each contig of each component was constructed. Finally, the graph was unlocked with true reads to obtain transcript sequences. Trinity combines reads of a certain length of overlap to form longer fragments and then processes them for sequencing. The resultant sequences were defined as unigenes for downstream analysis.

4.5. Functional Annotation and Analysis

The assembled contig sets were compared against various databases (NCBINR, Swissprot, GO, COG, KOG, and KEGG databases) using BLASTX (http://blast.ncbi.nlm.nih.gov/Blast.cgi) with an E-value threshold of 1 × 10−5 [80]. Predicted amino acid sequences of unigenes were searched against the Pfam [81] database using HMMER software (http://hmmer.org/) with an E-value ≤ 1 × 10−10 to gain the unigene annotation information.
GO is a controlled vocabulary of terms split into three related ontologies consisting of molecular function, biological processes, and cellular components. It has a hierarchical structure that forms a directed acyclic graph in which each term has defined relationships to one or more other terms in the same domain. This characteristic can be described as a parent–child relationship. Each GO term is represented by a node in this graph, and the nodes are annotated with a set of genes. GO annotation was performed from gene expression data using TopGO [82].
The KEGG pathways were assigned to the assembled unigenes using the KAAS server automatic annotator to summarize the pathway information involved in P. aegyptiaca seed germination [83]. The transcripts were aligned to sequences in the COG database (http://www.ncbi.nlm.nih.gov/COG) and categorized under different functions accordingly.

4.6. Differential Expression Analysis

EBseq is a widely accepted and accurate analysis of RNA-seq data [84]. This process was performed to identify DEGs. RPKM values were used to measure gene expression levels. FDR was used to determine the threshold of the p-value in multiple tests [85]. FDR < 0.01 and fold change ≥2 were considered the cutoff threshold to determine the significance of expression. Cluster of significantly regulated genes was assessed with Cluster [86], MeV [87] and TreeView [88].

4.7. qRT-PCR Validation

Total RNA for RNA-seq samples was also used for qRT-PCR analysis. First-strand cDNA was synthesized using M-MLV reverse transcriptase (Roche, Switzerland) and oligo (dT) 18. We chose 8 plant hormone-related genes, 3 SL perception-related genes and randomly selected 12 genes for further study and validation of RNA-seq using qRT-PCR. Each transcript primer (Table S1) was designed using Primer Premier 5.0 (Premier Biosoft International, Palo Alto, CA, USA). OaTublin1 was used as an internal control [89]. qRT-PCR was conducted using SYBR GreenER™ qPCRSuperMix Universal (Invitrogen, Gaithersburg, MD, USA) according to the manufacturer’s instructions. Thermal cycle conditions for PCR were as follows: 94 °C for 3 min, 40 cycles including 94 °C for 15 s and 60 °C for 30 s. Relative expression levels were calculated using the 2−ΔΔCt method [90] and one-way ANOVA (IBM SPSS Statistics 19.0, Armonk, NY, USA).

4.8. Endogenous Hormone Level Analysis

Approximately 1 g (weighed after 5 min of air-drying) of seeds for each sample (for RNA-seq) was collected. Fresh material was frozen in liquid nitrogen and ground to a fine powder. ABA, IAA, and GAs were extracted according to the method of Durgbanshi with slight modification [91]. BRs were extracted using the method of Kim [92]. The four plant hormones were quantified using appropriate ELISA kits (Chengling, Beijing, China) according to the manufacturer’s instructions [93]. The released ethylene level was quantified according to the method of Cristescu [94]. All quantifications were performed with triplicate samples of each treatment, and data were analyzed by one-way ANOVA (IBM SPSS Statistics 19.0, Armonk, NY, USA).

4.9. Seed Hormone Treatments

P. aegyptiaca seeds were prepared as previously described. IAA (0.1% ethanol, v/v), GA3 (0.1% ethanol, v/v), epi-brassinolide (0.1% ethanol, v/v), ethephon, ABA (0.1% ethanol, v/v), fluridone (0.1% acetone, v/v) (inhibitors of carotenoid biosynthesis), or distilled water (as a control) was added to each Petri dish.
Each hormone, inhibitor, or water was added before conditioning in the first set of treatments. Hormones, inhibitors, or water were applied after 7 days of conditioning as the second group of treatments. Germination percentage was examined under a binocular microscope (Zeiss, Oberkochen, Germany) 3 days after GR24 treatment. Each treatment was performed with three replicate Petri dishes. Data were analyzed by one-way ANOVA (IBM SPSS Statistics 19.0, Armonk, NY, USA).

4.10. Accession Numbers

The quality filtered and trimmed short read data set, which was used in this study, was deposited to the NCBI Short Read Archive (SRA) under accession numbers: SRR2477710, SRR2477711, and SRR2477712.

5. Conclusions

This study generated a well-annotated transcriptome for the root parasitic weed, P. aegyptiaca. The taxonomy of the transcripts of this plant was highly similar to that of S. indicum. Thus, as expected, the parasite belongs to the same family of S. indicum. We clarified DNA repair and replication, RNA biosynthesis, and energy metabolism, as well as the transcription, translation, and biosynthesis of proteins, in this plant. We also found several key genes involved in seed germination, and these genes were associated with changes in the plant hormones. These genes may be developed as new targets to control P. aegyptiaca. These genes were validated by q-PCR and hormones levels were verified in seed germination experiments. Fluridone alone induced P. aegyptiaca seed germination, whereas ABA and ethylene were inferred as the key phytohormones during P. aegyptiaca seed germination. Our results may provide a perspective for suicidal germination to control parasitic weeds. In addition, we found that the expression levels of MAX2 and a KAI2 gene, which were involved in SL signal transduction, were always stable during the whole germination process. Thus, we deduced that transcripts of partial components of the SL perception complex existed prior to conditioning.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/17/7/1139/s1.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 31460467).

Author Contributions

Zhaoqun Yao and Sifeng Zhao conceived and designed the experiments; Fang Tian, Xiaolei Cao and Meixiu Chen performed the experiments; Zhaoqun Yao and Ying Xu analyzed the data; Benchun Xiang and Sifeng Zhao contributed reagents/materials/analysis tools; Zhaoqun Yao wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Parker, C. Parasitic weeds: A world challenge. Weed Sci. 2012, 60, 269–276. [Google Scholar] [CrossRef]
  2. Yoder, J.I.; Scholes, J.D. Host plant resistance to parasitic weeds; recent progress and bottlenecks. Curr. Opin. Plant Biol. 2010, 13, 478–484. [Google Scholar] [CrossRef] [PubMed]
  3. Ghannam, I.; Al-Masri, M.; Barakat, R. The effect of herbicides on the egyptian broomrape (Orobanche aegyptiaca) in tomato fields. Am. J. Plant Sci. 2012, 3, 346. [Google Scholar] [CrossRef]
  4. Westwood, J.H.; Yoder, J.I.; Timko, M.P. The evolution of parasitism in plants. Trends Plant Sci. 2010, 15, 227–235. [Google Scholar] [CrossRef] [PubMed]
  5. Parker, C. Observations on the current status of Orobanche and Striga problems worldwide. Pest. Manag. Sci. 2009, 65, 453–459. [Google Scholar] [CrossRef] [PubMed]
  6. Joel, D.; Hershenhorn, J.; Eizenberg, H.; Aly, R.; Ejeta, G.; Rich, P.; Ransom, J.; Sauerborn, J.; Rubiales, D. Biology and management of weedy root parasites. Horticult. Rev. Westport N. Y. 2007, 33, 267. [Google Scholar]
  7. Parker, C.; Riches, C.R. Parasitic Weeds of the World: Biology and Control; CAB International: Wallingford, UK, 1993. [Google Scholar]
  8. Press, M.C.; Graves, J.D. Parasitic Plants; Chapman & Hall Ltd: London, UK, 1995. [Google Scholar]
  9. Joel, D.M. The long-term approach to parasitic weeds control: Manipulation of specific developmental mechanisms of the parasite. Crop Prot. 2000, 19, 753–758. [Google Scholar] [CrossRef]
  10. Xie, X.; Yoneyama, K.; Yoneyama, K. The strigolactone story. Annu. Rev. Phytopathol. 2010, 48, 93–117. [Google Scholar] [CrossRef] [PubMed]
  11. Koltai, H. Receptors, repressors, pins: A playground for strigolactone signaling. Trends Plant Sci. 2014, 19, 727–733. [Google Scholar] [CrossRef] [PubMed]
  12. Waldie, T.; Hayward, A.; Beveridge, C.A. Axillary bud outgrowth in herbaceous shoots: How do strigolactones fit into the picture? Plant Mol. Biol. 2010, 73, 27–36. [Google Scholar] [CrossRef] [PubMed]
  13. Ruyter-Spira, C.; Al-Babili, S.; van der Krol, S.; Bouwmeester, H. The biology of strigolactones. Trends Plant Sci. 2013, 18, 72–83. [Google Scholar] [CrossRef] [PubMed]
  14. Babaei, S.; Alizadeh, H.; Jahansouz, M.R.; Mashhadi, H.R.; Moeini, M.M. Management of Phelipanche aegyptiaca pomel: Using trap crops in rotation with tomato (Solanum lycopersicom L.). Aust. J. Crop Sci. 2010, 4, 437–442. [Google Scholar]
  15. Grenz, J.; Iştoc, V.; Manschadi, A.; Sauerborn, J. Interactions of sunflower (Helianthus annuus) and sunflower broomrape (Orobanche cumana) as affected by sowing date, resource supply and infestation level. Field Crop. Res. 2008, 107, 170–179. [Google Scholar] [CrossRef]
  16. Kannan, C.; Aditi, P.; Zwanenburg, B. Quenching the action of germination stimulants using borax and thiourea, a new method for controlling parasitic weeds: A proof of concept. Crop Prot. 2015, 70, 92–98. [Google Scholar] [CrossRef]
  17. MÜLLER-STÖVER, D.; Kohlschmid, E.; Sauerborn, J. A novel strain of Fusarium oxysporum from germany and its potential for biocontrol of Orobanche ramosa. Weed Res. 2009, 49, 175–182. [Google Scholar] [CrossRef]
  18. Sillero, J.C.; Rojas-Molina, M.M.; Ávila, C.M.; Rubiales, D. Induction of systemic acquired resistance against rust, ascochyta blight and broomrape in faba bean by exogenous application of salicylic acid and benzothiadiazole. Crop Prot. 2012, 34, 65–69. [Google Scholar] [CrossRef]
  19. López-Ráez, J.A.; Matusova, R.; Cardoso, C.; Jamil, M.; Charnikhova, T.; Kohlen, W.; Ruyter-Spira, C.; Verstappen, F.; Bouwmeester, H. Strigolactones: Ecological significance and use as a target for parasitic plant control. Pest Manag. Sci. 2009, 65, 471–477. [Google Scholar] [CrossRef] [PubMed]
  20. Dor, E.; Yoneyama, K.; Wininger, S.; Kapulnik, Y.; Yoneyama, K.; Koltai, H.; Xie, X.; Hershenhorn, J. Strigolactone deficiency confers resistance in tomato line SL-ORT1 to the parasitic weeds Phelipanche and Orobanche spp. Phytopathology 2011, 101, 213–222. [Google Scholar] [CrossRef] [PubMed]
  21. Vogel, J.T.; Walter, M.H.; Giavalisco, P.; Lytovchenko, A.; Kohlen, W.; Charnikhova, T.; Simkin, A.J.; Goulet, C.; Strack, D.; Bouwmeester, H.J. SlCCD7 controls strigolactone biosynthesis, shoot branching and mycorrhiza-induced apocarotenoid formation in tomato. Plant J. 2010, 61, 300–311. [Google Scholar] [CrossRef] [PubMed]
  22. López-Ráez, J.A.; Charnikhova, T.; Fernández, I.; Bouwmeester, H.; Pozo, M.J. Arbuscular mycorrhizal symbiosis decreases strigolactone production in tomato. J. Plant Physiol. 2011, 168, 294–297. [Google Scholar] [CrossRef] [PubMed]
  23. Zwanenburg, B.; Mwakaboko, A.S.; Kannan, C. Perspective of suicidal germination for parasitic weeds control. Pest Manag. Sci. 2016. [Google Scholar] [CrossRef] [PubMed]
  24. Westwood, J.H.; Depamphilis, C.W.; Das, M.; Fernándezaparicio, M.; Honaas, L.A.; Timko, M.P.; Wafula, E.K.; Wickett, N.J.; Yoder, J.I. The parasitic plant genome project: New tools for understanding the biology of Orobanche and Striga. Weed Sci. 2012, 60, 295–306. [Google Scholar] [CrossRef]
  25. Ruiu, F.; Picarella, M.; Imanishi, S.; Mazzucato, A. A transcriptomic approach to identify regulatory genes involved in fruit set of wild-type and parthenocarpic tomato genotypes. Plant Mol. Biol. 2015, 89, 263–278. [Google Scholar] [CrossRef] [PubMed]
  26. Honaas, L.A.; Wafula, E.K.; Yang, Z.; Der, J.P.; Wickett, N.J.; Altman, N.S.; Taylor, C.G.; Yoder, J.I.; Timko, M.P.; Westwood, J.H. Functional genomics of a generalist parasitic plant: Laser microdissection of host-parasite interface reveals host-specific patterns of parasite gene expression. BMC Plant Biol. 2013, 13, 9. [Google Scholar] [CrossRef] [PubMed]
  27. Ranjan, A.; Ichihashi, Y.; Farhi, M.; Zumstein, K.; Townsley, B.; David-Schwartz, R.; Sinha, N.R. De novo assembly and characterization of the transcriptome of the parasitic weed dodder identifies genes associated with plant parasitism. Plant Physiol. 2014, 166, 1186–1199. [Google Scholar] [CrossRef] [PubMed]
  28. Wickett, N.J.; Honaas, L.A.; Wafula, E.K.; Das, M.; Huang, K.; Wu, B.; Landherr, L.; Timko, M.P.; Yoder, J.; Westwood, J.H. Transcriptomes of the parasitic plant family orobanchaceae reveal surprising conservation of chlorophyll synthesis. Curr. Biol. 2011, 21, 2098–2104. [Google Scholar] [CrossRef] [PubMed]
  29. Nambara, E. A comparative transcriptome analysis on growth regulation between seeds and buds of arabidopsis. Available online: https://pag.confex.com/pag/xxiii/webprogram/Paper15834.html (accessed on 13 July 2016).
  30. Tao, T.; Zhao, L.; Lv, Y.; Chen, J.; Hu, Y.; Zhang, T.; Zhou, B. Transcriptome sequencing and differential gene expression analysis of delayed gland morphogenesis in Gossypium australe during seed germination. PLoS ONE 2013, 8, e75323. [Google Scholar] [CrossRef] [PubMed]
  31. Yu, Y.; Guo, G.; Lv, D.; Hu, Y.; Li, J.; Li, X.; Yan, Y. Transcriptome analysis during seed germination of elite chinese bread wheat cultivar Jimai 20. BMC Plant Biol. 2014, 14, 1–19. [Google Scholar] [CrossRef] [PubMed]
  32. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q. Full-Length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  33. Gonzalez-Verdejo, C.I.; Barandiaran, X.; Moreno, M.T.; Cubero, J.I.; di Pietro, A. An improved axenic system for studying pre-infection development of the parasitic plant Orobanche ramosa. Ann. Bot.-Lond. 2005, 96, 1121–1127. [Google Scholar] [CrossRef] [PubMed]
  34. Kusumoto, D.; Chae, S.H.; Mukaida, K.; Yoneyama, K.; Yoneyama, K.; Joel, D.M.; Takeuchi, Y. Effects of fluridone and norflurazon on conditioning and germination of Striga asiatica seeds. Plant Growth Regul. 2006, 48, 73–78. [Google Scholar] [CrossRef]
  35. Joel, D.; Steffens, J.; Matthews, D. Germination of weedy root parasites. Seed Dev. Germination 1995, 1, 567–597. [Google Scholar]
  36. Tsuchiya, Y.; Yoshimura, M.; Sato, Y.; Kuwata, K.; Toh, S.; Holbrook-Smith, D.; Zhang, H.; McCourt, P.; Itami, K.; Kinoshita, T. Probing strigolactone receptors in Striga hermonthica with fluorescence. Science 2015, 349, 864–868. [Google Scholar] [CrossRef] [PubMed]
  37. Liu, S.; Li, W.; Wu, Y.; Chen, C.; Lei, J. De novo transcriptome assembly in chili pepper (Capsicum frutescens) to identify genes involved in the biosynthesis of capsaicinoids. PLoS ONE 2013, 8, e48156. [Google Scholar] [CrossRef] [PubMed]
  38. Nakasugi, K.; Crowhurst, R.N.; Bally, J.; Wood, C.C.; Hellens, R.P.; Waterhouse, P.M. De novo transcriptome sequence assembly and analysis of rna silencing genes of Nicotiana benthamiana. PLoS ONE 2013, 8, e59534. [Google Scholar] [CrossRef] [PubMed]
  39. Wang, L.; Yu, S.; Tong, C.; Zhao, Y.; Liu, Y.; Song, C.; Zhang, Y.; Zhang, X.; Wang, Y.; Hua, W. Genome sequencing of the high oil crop sesame provides insight into oil biosynthesis. Genome Biol. 2014, 15, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Dal Degan, F.; Rocher, A.; Cameron-Mills, V.; Von Wettstein, D. The expression of serine carboxypeptidases during maturation and germination of the barley grain. Proc. Natl. Acad. Sci. USA 1994, 91, 8209–8213. [Google Scholar] [CrossRef] [PubMed]
  41. Dietrich, R.A.; Maslyar, D.J.; Heupel, R.C.; Harada, J.J. Spatial patterns of gene expression in Brassica napus seedlings: Identification of a cortex-specific gene and localization of mrnas encoding isocitrate lyase and a polypeptide homologous to proteinases. Plant Cell 1989, 1, 73–80. [Google Scholar] [CrossRef] [PubMed]
  42. Rajjou, L.; Duval, M.; Gallardo, K.; Catusse, J.; Bally, J.; Job, C.; Job, D. Seed germination and vigor. Annu. Rev. Plant Biol. 2012, 63, 507–533. [Google Scholar] [CrossRef] [PubMed]
  43. Howell, K.A.; Narsai, R.; Carroll, A.; Ivanova, A.; Lohse, M.; Usadel, B.; Millar, A.H.; Whelan, J. Mapping metabolic and transcript temporal switches during germination in rice highlights specific transcription factors and the role of rna instability in the germination process. Plant Physiol. 2009, 149, 961–980. [Google Scholar] [CrossRef] [PubMed]
  44. Xu, J.; Chua, N.-H. Arabidopsis decapping 5 is required for mrna decapping, p-body formation, and translational repression during postembryonic development. Plant Cell 2009, 21, 3270–3279. [Google Scholar] [CrossRef] [PubMed]
  45. Hunt, L.; Holdsworth, M.J.; Gray, J.E. Nicotinamidase activity is important for germination. Plant J. 2007, 51, 341–351. [Google Scholar] [CrossRef] [PubMed]
  46. Macovei, A.; Balestrazzi, A.; Confalonieri, M.; Faé, M.; Carbonera, D. New insights on the barrel medic MtOOG1 and MtFPQ functions in relation to oxidative stress response in planta and during seed imbibition. Plant Physiol. Biochem. 2011, 49, 1040–1050. [Google Scholar] [CrossRef] [PubMed]
  47. Rajjou, L.; Gallardo, K.; Debeaujon, I.; Vandekerckhove, J.; Job, C.; Job, D. The effect of α-amanitin on the Arabidopsis seed proteome highlights the distinct roles of stored and neosynthesized mrnas during germination. Plant Physiol. 2004, 134, 1598–1613. [Google Scholar] [CrossRef] [PubMed]
  48. Schmitz-Linneweber, C.; Small, I. Pentatricopeptide repeat proteins: A socket set for organelle gene expression. Trends Plant Sci. 2008, 13, 663–670. [Google Scholar] [CrossRef] [PubMed]
  49. Bewley, J.D.; Bradford, K.J.; Hilhorst, H.W.; Nonogaki, H. Development and maturation. In Seeds; Springer: Berlin, Germany, 2013; pp. 27–83. [Google Scholar]
  50. Ogé, L.; Bourdais, G.; Bove, J.; Collet, B.; Godin, B.; Granier, F.; Boutin, J.-P.; Job, D.; Jullien, M.; Grappin, P. Protein repair l-isoaspartyl methyltransferase1 is involved in both seed longevity and germination vigor in Arabidopsis. Plant Cell 2008, 20, 3022–3037. [Google Scholar] [CrossRef] [PubMed]
  51. Weitbrecht, K.; Müller, K.; Leubner-Metzger, G. First off the mark: Early seed germination. J. Exp. Bot. 2011, 62, 3289–3309. [Google Scholar] [CrossRef] [PubMed]
  52. Wind, J.J.; Peviani, A.; Snel, B.; Hanson, J.; Smeekens, S.C. Abi4: Versatile activator and repressor. Trends Plant Sci. 2013, 18, 125–132. [Google Scholar] [CrossRef] [PubMed]
  53. Yang, R.; Guo, Q.; Gu, Z. Gaba shunt and polyamine degradation pathway on γ-aminobutyric acid accumulation in germinating fava bean (Vicia faba L.) under hypoxia. Food Chem. 2013, 136, 152–159. [Google Scholar] [CrossRef] [PubMed]
  54. Linkies, A.; Leubner-Metzger, G. Beyond gibberellins and abscisic acid: How ethylene and jasmonates control seed germination. Plant Cell Rep. 2012, 31, 253–270. [Google Scholar] [CrossRef] [PubMed]
  55. Chae, S.H.; Yoneyama, K.; Takeuchi, Y.; Joel, D.M. Fluridone and norflurazon, carotenoid-biosynthesis inhibitors, promote seed conditioning and germination of the holoparasite Orobanche minor. Physiol. Plant. 2004, 120, 328–337. [Google Scholar] [CrossRef] [PubMed]
  56. Takeuchi, Y.; Omigawa, Y.; Ogasawara, M.; Yoneyama, K.; Konnai, M.; Worsham, A.D. Effects of brassinosteroids on conditioning and germination of clover broomrape (Orobanche minor) seeds. Plant Growth Regul. 1995, 16, 153–160. [Google Scholar] [CrossRef]
  57. Zehhar, N.; Ingouff, M.; Bouya, D.; Fer, A. Possible involvement of gibberellins and ethylene in Orobanche ramosa germination. Weed Res. 2002, 42, 464–469. [Google Scholar] [CrossRef]
  58. Lechat, M.-M.; Pouvreau, J.-B.; Péron, T.; Gauthier, M.; Montiel, G.; Véronési, C.; Todoroki, Y.; Le Bizec, B.; Monteau, F.; Macherel, D. PrCYP707A1, an aba catabolic gene, is a key component of Phelipanche ramosa seed germination in response to the strigolactone analogue GR24. J. Exp. Bot. 2012, 63, 5311–5322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Lee, H.; Lee, K.; Seo, P. The Arabidopsis MYB96 transcription factor plays a role in seed dormancy. Plant Mol. Biol. 2015, 87, 371–381. [Google Scholar] [CrossRef] [PubMed]
  60. Lechat, M.-M.; Brun, G.; Montiel, G.; Véronési, C.; Simier, P.; Thoiron, S.; Pouvreau, J.-B.; Delavault, P. Seed response to strigolactone is controlled by abscisic acid-independent DNA methylation in the obligate root parasitic plant, Phelipanche ramosa L. Pomel. J. Exp. Bot. 2015, 66, 1329–1340. [Google Scholar] [CrossRef] [PubMed]
  61. Matusova, R.; Rani, K.; Verstappen, F.W.; Franssen, M.C.; Beale, M.H.; Bouwmeester, H.J. The strigolactone germination stimulants of the plant-parasitic Striga and Orobanche spp. Are derived from the carotenoid pathway. Plant Physiol. 2005, 139, 920–934. [Google Scholar] [CrossRef] [PubMed]
  62. Song, W.; Zhou, W.; Jin, Z.; Cao, D.; Joel, D.; Takeuchi, Y.; Yoneyama, K. Germination response of Orobanche seeds subjected to conditioning temperature, water potential and growth regulator treatments. Weed Res. 2005, 45, 467–476. [Google Scholar] [CrossRef]
  63. Song, W.; Zhou, W.; Jin, Z.; Zhang, D.; Yoneyama, K.; Takeuchi, Y.; Joel, D. Growth regulators restore germination of Orobanche seeds that are conditioned under water stress and suboptimal temperature. Crop Pasture Sci. 2006, 57, 1195–1201. [Google Scholar] [CrossRef]
  64. Joel, D.; Back, A.; Kleifeld, Y.; Gepstein, S. Seed conditioning and its role in orobanche seed germination: Inhibition by paclobutrazol. In Progress in Orobanche Research, Proceedings of the International Workshop on Orobanche Research, Obermarchtal, Germany, 19–22 August 1989.
  65. Uematsu, K.; Nakajima, M.; Yamaguchi, I.; Yoneyama, K.; Fukui, Y. Role of camp in gibberellin promotion of seed germination in Orobanche minor Smith. J. Plant Growth Regul. 2007, 26, 245–254. [Google Scholar] [CrossRef]
  66. Krizek, D.T.; Mandava, N.B. Influence of spectral quality on the growth response of intact bean plants to brassinosteroid, a growth-promoting steroidal lactone. I. Stem elongation and morphogenesis. Physiol. Plant. 1983, 57, 317–323. [Google Scholar] [CrossRef]
  67. Graeber, K.; Linkies, A.; Müller, K.; Wunchova, A.; Rott, A.; Leubner-Metzger, G. Cross-Species approaches to seed dormancy and germination: Conservation and biodiversity of aba-regulated mechanisms and the brassicaceae DOG1 genes. Plant Mol. Biol. 2010, 73, 67–87. [Google Scholar] [CrossRef] [PubMed]
  68. Morris, K.; Linkies, A.; Müller, K.; Oracz, K.; Wang, X.; Lynn, J.R.; Leubner-Metzger, G.; Finch-Savage, W.E. Regulation of seed germination in the close arabidopsis relative Lepidium sativum: A global tissue-specific transcript analysis. Plant Physiol. 2011, 155, 1851–1870. [Google Scholar] [CrossRef] [PubMed]
  69. Kucera, B.; Cohn, M.A.; Leubner-Metzger, G. Plant hormone interactions during seed dormancy release and germination. Seed Sci. Res. 2005, 15, 281–307. [Google Scholar] [CrossRef]
  70. Linkies, A.; Müller, K.; Morris, K.; Turečková, V.; Wenk, M.; Cadman, C.S.; Corbineau, F.; Strnad, M.; Lynn, J.R.; Finch-Savage, W.E. Ethylene interacts with abscisic acid to regulate endosperm rupture during germination: A comparative approach using Lepidium sativum and Arabidopsis thaliana. Plant Cell 2009, 21, 3803–3822. [Google Scholar] [CrossRef] [PubMed]
  71. Matilla, A.; Matilla-Vazquez, M. Involvement of ethylene in seed physiology. Plant Sci. 2008, 175, 87–97. [Google Scholar] [CrossRef]
  72. Chun, D.; Wilhelm, S.; Sagen, J.; Musselman, L.; Worsham, A.; Eplee, R. Components of record germination in vitro of branched broomrape, orobanche ramosa L. In Supplement to the Proceedings of the Second International Symposium on Parasitic Weeds, Raleigh, NC, USA, 16–19 July 1979.
  73. Bebawi, F.F.; Eplee, R.E. Efficacy of ethylene as a germination stimulant of Striga hermonthica seed. Weed Sci. 1986, 34, 694–698. [Google Scholar]
  74. Egley, G.; Dale, J. Ethylene, 2-chloroethylphosphonic acid, and witchweed germination. Weed Sci. 1970, 18, 586–589. [Google Scholar]
  75. Arc, E.; Sechet, J.; Corbineau, F.; Rajjou, L.; Marion-Poll, A. ABA crosstalk with ethylene and nitric oxide in seed dormancy and germination. Front. Plant Sci. 2013, 4, 63. [Google Scholar] [CrossRef] [PubMed]
  76. Toh, S.; Holbrook-Smith, D.; Stogios, P.J.; Onopriyenko, O.; Lumba, S.; Tsuchiya, Y.; Savchenko, A.; McCourt, P. Structure-function analysis identifies highly sensitive strigolactone receptors in Striga. Science 2015, 350, 203–207. [Google Scholar] [CrossRef] [PubMed]
  77. Westwood, J.H. Characterization of the Orobanche–Arabidopsis system for studying parasite–host interactions. Weed Sci. 2009, 48, 742–748. [Google Scholar]
  78. Fernández-Aparicio, M.; Rubiales, D.; Bandaranayake, P.C.; Yoder, J.I.; Westwood, J.H. Transformation and regeneration of the holoparasitic plant Phelipanche aegyptiaca. Plant Methods 2011, 7, 683–696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Margulies, M.; Egholm, M.; Altman, W.E.; Attiya, S.; Bader, J.S.; Bemben, L.A.; Berka, J.; Braverman, M.S.; Chen, Y.-J.; Chen, Z. Genome sequencing in microfabricated high-density picolitre reactors. Nature 2005, 437, 376–380. [Google Scholar] [CrossRef] [PubMed]
  80. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped blast and psi-blast: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [PubMed]
  81. Finn, R.D.; Bateman, A.; Clements, J.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Heger, A.; Hetherington, K.; Holm, L.; Mistry, J. Pfam: The protein families database. Nucleic Acids Res. 2013, 42, D222–D230. [Google Scholar] [CrossRef] [PubMed]
  82. Alexa, A.; Rahnenfuhrer, J. topGo: Enrichment analysis for gene ontology. Available online: http://bioconductor.uib.no/2.7/bioc/html/topGO.html (accessed on 13 July 2016).
  83. Moriya, Y.; Itoh, M.; Okuda, S.; Yoshizawa, A.C.; Kanehisa, M. Kaas: An automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35, W182–W185. [Google Scholar] [CrossRef] [PubMed]
  84. Leng, N.; Dawson, J.; Kendziorski, C. Ebseq: An R package for gene and isoform differential expression analysis of rna-seq data, R package version 1.5.3. Available online: https://www.biostat.wisc.edu/~kendzior/EBSEQ/ (accessed on 13 July 2016).
  85. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar]
  86. Eisen, M.B.; Spellman, P.T.; Brown, P.O.; Botstein, D. Cluster analysis and display of genome-wide expression patterns. Proc. Natl. Acad. Sci. USA 1998, 95, 14863–14868. [Google Scholar] [CrossRef] [PubMed]
  87. Howe, E.A.; Sinha, R.; Schlauch, D.; Quackenbush, J. Rna-seq analysis in mev. Bioinformatics 2011, 27, 3209–3210. [Google Scholar] [CrossRef] [PubMed]
  88. Saldanha, A.J. Java treeview—extensible visualization of microarray data. Bioinformatics 2004, 20, 3246–3248. [Google Scholar] [CrossRef] [PubMed]
  89. Gonzalez-Verdejo, C.; Die, J.; Nadal, S.; Jimenez-Marin, A.; Moreno, M.; Roman, B. Selection of housekeeping genes for normalization by real-time RT-PCR: Analysis of Or-MYB1 gene expression in Orobanche ramosa development. Anal. Biochem. 2008, 379, 176–181. [Google Scholar] [CrossRef] [PubMed]
  90. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−∆∆Ct method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  91. Durgbanshi, A.; Arbona, V.; Pozo, O.; Miersch, O.; Sancho, J.V.; Gómez-Cadenas, A. Simultaneous determination of multiple phytohormones in plant extracts by liquid chromatography-electrospray tandem mass spectrometry. J. Agric. Food Chem. 2005, 53, 8437–8442. [Google Scholar] [CrossRef] [PubMed]
  92. Kim, T.-W.; Hwang, J.-Y.; Kim, Y.-S.; Joo, S.-H.; Chang, S.C.; Lee, J.S.; Takatsuto, S.; Kim, S.-K. Arabidopsis CYP85A2, a cytochrome P450, mediates the baeyer-villiger oxidation of castasterone to brassinolide in brassinosteroid biosynthesis. Plant Cell 2005, 17, 2397–2412. [Google Scholar] [CrossRef] [PubMed]
  93. Zhou, M.; Xu, M.; Wu, L.; Shen, C.; Ma, H.; Lin, J. CbCBF from Capsella bursa-pastoris enhances cold tolerance and restrains growth in Nicotiana tabacum by antagonizing with gibberellin and affecting cell cycle signaling. Plant Mol. Biol. 2014, 85, 259–275. [Google Scholar] [CrossRef] [PubMed]
  94. Cristescu, S.M.; Woltering, E.; Hermans, C.; Harren, F.J.; te Lintel Hekkert, S. Research tools: Ethylene detection. In Ethylene in Plants; Springer: Berlin, Germany, 2015; pp. 263–286. [Google Scholar]
Ijms 17 01139 g001 1024
Figure 1. Length distribution of Trinity assembly for unigenes of Phelipanche aegyptiaca.
Figure 1. Length distribution of Trinity assembly for unigenes of Phelipanche aegyptiaca.
Ijms 17 01139 g001
Ijms 17 01139 g002 1024
Figure 2. Hit species distribution of BLASTX matches of P. aegyptiaca transcriptome. Number and proportion of P. aegyptiaca transcripts with similarity to sequences from the NR database.
Figure 2. Hit species distribution of BLASTX matches of P. aegyptiaca transcriptome. Number and proportion of P. aegyptiaca transcripts with similarity to sequences from the NR database.
Ijms 17 01139 g002
Ijms 17 01139 g003 1024
Figure 3. GO category distribution of P. aegyptiaca transcripts among level 1 GO categories: Biological process, molecular function, and cellular component.
Figure 3. GO category distribution of P. aegyptiaca transcripts among level 1 GO categories: Biological process, molecular function, and cellular component.
Ijms 17 01139 g003
Ijms 17 01139 g004 1024
Figure 4. Top 20 KEGG pathways in P. aegyptiaca germinating seeds according to the numbers of DEGs. This Figure shows the KEGG metabolic pathways of plants that were well represented by unique sequences of P. aegyptiaca.
Figure 4. Top 20 KEGG pathways in P. aegyptiaca germinating seeds according to the numbers of DEGs. This Figure shows the KEGG metabolic pathways of plants that were well represented by unique sequences of P. aegyptiaca.
Ijms 17 01139 g004
Ijms 17 01139 g005 1024
Figure 5. Validation of hormone-related genes at the GR24 treatment stage by qRT-PCR. Fold changes were originally estimated in the RNA-seq experiment (brown bar) and validated by qRT-PCR (white bar). The results are the average of three biological replicates. Abbreviations of genes: PaGA2OX1, gibberellin 2-beta-dioxygenase 1; PaGA3OX1, gibberellin 3-beta-dioxygenase 1; PaNCED, 9-cis-epoxycarotenoid dioxygenase; PaCYP707A1, ABA 8′-hydroxylase 1; PaMAT, methionine adenosyltransferase; PaACS, amino cyclopropanecarboxylate synthase; PaACO, aminocyclopropanecarboxylate oxidase; PaCYP90B1, cytochrome P450 90B1.
Figure 5. Validation of hormone-related genes at the GR24 treatment stage by qRT-PCR. Fold changes were originally estimated in the RNA-seq experiment (brown bar) and validated by qRT-PCR (white bar). The results are the average of three biological replicates. Abbreviations of genes: PaGA2OX1, gibberellin 2-beta-dioxygenase 1; PaGA3OX1, gibberellin 3-beta-dioxygenase 1; PaNCED, 9-cis-epoxycarotenoid dioxygenase; PaCYP707A1, ABA 8′-hydroxylase 1; PaMAT, methionine adenosyltransferase; PaACS, amino cyclopropanecarboxylate synthase; PaACO, aminocyclopropanecarboxylate oxidase; PaCYP90B1, cytochrome P450 90B1.
Ijms 17 01139 g005
Ijms 17 01139 g006 1024
Figure 6. Heat map from hierarchical clustering of differentially expressed genes during P. aegyptiaca seed germination. Expression changes and cluster analysis of 5324 genes that were differentially expressed between any two of the three samples. Each row represents a differentially expressed gene (DEG), whereas each column represents a sample. Changes in expression levels are shown in color scales with saturation at >2.0-fold changes. Green and red color gradients indicate a decrease and increase in transcript abundance, respectively.
Figure 6. Heat map from hierarchical clustering of differentially expressed genes during P. aegyptiaca seed germination. Expression changes and cluster analysis of 5324 genes that were differentially expressed between any two of the three samples. Each row represents a differentially expressed gene (DEG), whereas each column represents a sample. Changes in expression levels are shown in color scales with saturation at >2.0-fold changes. Green and red color gradients indicate a decrease and increase in transcript abundance, respectively.
Ijms 17 01139 g006
Ijms 17 01139 g007 1024
Figure 7. Time course of ABA, GA3, and brassinosteroid (BR) contents in P. aegyptiaca seed germination. Seeds were conditioned in water at 25 °C in the dark for 7 days and then conditioned in GR24 (10−6 mol/L) at 25 °C for 1 day.
Figure 7. Time course of ABA, GA3, and brassinosteroid (BR) contents in P. aegyptiaca seed germination. Seeds were conditioned in water at 25 °C in the dark for 7 days and then conditioned in GR24 (10−6 mol/L) at 25 °C for 1 day.
Ijms 17 01139 g007
Table
Table 1. Assembly statistics for Phelipanche aegyptiaca transcriptome.
Table 1. Assembly statistics for Phelipanche aegyptiaca transcriptome.
Length RangeContigTranscriptUnigene
Total number11,469,131274,96494,419
Total length (bp)670,369,345293,096,15666,808,988
N50 (bp)531558836
Average (bp)58.451065.94707.58
Table
Table 2. P. aegyptiaca seed germination rate under different conditions.
Table 2. P. aegyptiaca seed germination rate under different conditions.
TreatmentGermination (%)
No ConditionedConditionedConditioned and GR24 Treatment
BR0 a0 a83.3 ± 3.8 fg
GA30 a2.2 ± 2.2 a94.4 ± 1.1 h
ABA0 a0 a0 a
Ethephon0 a0 a90.0 ± 1.9 gh
Fluridone26.5 ± 1.8 b7.8 ± 2.9 a86.7 ± 5.7 gh
GA3 + Fluridone35.6 ± 1.1 c10.0 ± 3.3 a91.1 ± 2.2 gh
Ethephon + Fluridone24.4 ± 2.2 b7.8 ± 1.1 a80.0 ± 5.0 f
GA3 + Ethephon + Fluridone47.8 ± 2.9 e20.0 ± 1.9 b91.1 ± 4.0 gh
Water0 a0 a75.6 ± 9.1 f
“No conditioned” means that the dormant seeds were treated for 9 days directly without conditioning. “Conditioned” represents seeds conditioned in water at 25 °C in the dark for 7 days and then treated with plant hormones or inhibitor for 2 days. “Conditioned and GR24 treatment” involved seeds conditioned in water at 25 °C in the dark for 7 days and then treated with GR24 and other plant hormones or inhibitor for 2 days. The concentrations of hormones or inhibitor were as follows: BR, 1 mg/L; GA3, 10 mg/L; ABA, 1 mg/L; fluridone, 100 mg/L; ethephon, 500 mg/L at pH 5; and GR24, 10−6 mol/L, all the hormones were dissolved in 0.1% ethanol (v/v), except GR24 and fluridone were dissolved in 0.1% acetone (v/v). Means within a column followed by the same letter do not differ significantly according to Fisher’s protected LSD test (p < 0.05).

Share and Cite

MDPI and ACS Style

Yao, Z.; Tian, F.; Cao, X.; Xu, Y.; Chen, M.; Xiang, B.; Zhao, S. Global Transcriptomic Analysis Reveals the Mechanism of Phelipanche aegyptiaca Seed Germination. Int. J. Mol. Sci. 2016, 17, 1139. https://doi.org/10.3390/ijms17071139

AMA Style

Yao Z, Tian F, Cao X, Xu Y, Chen M, Xiang B, Zhao S. Global Transcriptomic Analysis Reveals the Mechanism of Phelipanche aegyptiaca Seed Germination. International Journal of Molecular Sciences. 2016; 17(7):1139. https://doi.org/10.3390/ijms17071139

Chicago/Turabian Style

Yao, Zhaoqun, Fang Tian, Xiaolei Cao, Ying Xu, Meixiu Chen, Benchun Xiang, and Sifeng Zhao. 2016. "Global Transcriptomic Analysis Reveals the Mechanism of Phelipanche aegyptiaca Seed Germination" International Journal of Molecular Sciences 17, no. 7: 1139. https://doi.org/10.3390/ijms17071139

APA Style

Yao, Z., Tian, F., Cao, X., Xu, Y., Chen, M., Xiang, B., & Zhao, S. (2016). Global Transcriptomic Analysis Reveals the Mechanism of Phelipanche aegyptiaca Seed Germination. International Journal of Molecular Sciences, 17(7), 1139. https://doi.org/10.3390/ijms17071139

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop