Next Article in Journal
Neuropeptide Substance-P-Conjugated Chitosan Nanofibers as an Active Modulator of Stem Cell Recruiting
Next Article in Special Issue
Genome-Wide Identification, Characterization and Expression Analysis of the Chalcone Synthase Family in Maize
Previous Article in Journal
Earthworm Lumbricus rubellus MT-2: Metal Binding and Protein Folding of a True Cadmium-MT
Previous Article in Special Issue
Anion Channel Inhibitor NPPB-Inhibited Fluoride Accumulation in Tea Plant (Camellia sinensis) Is Related to the Regulation of Ca2+, CaM and Depolarization of Plasma Membrane Potential
 
 
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 1
10.3390/ijms17010042
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

Comparative Transcriptome Profile of the Cytoplasmic Male Sterile and Fertile Floral Buds of Radish (Raphanus sativus L.)

by
Shiyong Mei
1,*,
Touming Liu
2 and
Zhiwei Wang
3
1
College of Life Sciences, Wuhan University, Wuhan 430072, China
2
Institute of Bast Fiber Crops and Center of Southern Economic Crops, Chinese Academy of Agricultural Sciences, Changsha 410205, China
3
CAS Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(1), 42; https://doi.org/10.3390/ijms17010042
Submission received: 27 October 2015 / Revised: 22 December 2015 / Accepted: 23 December 2015 / Published: 6 January 2016
(This article belongs to the Special Issue Plant Molecular Biology)
Browse Figures
Versions Notes

Abstract

:
Radish cytoplasmic male sterility (CMS) has been widely used for breeding in Raphanus and Brassica genera. However, the detailed regulation network of the male sterility remains to be determined. Our previous work has shown that the abnormalities in a CMS radish appeared shortly after the tetrad stage when microspores were malformed and the tapetal cells grew abnormally large. In this work, histological analysis shows that anthers are at the tetrad stage when the radish buds are about 1.5 mm in length. Furthermore, a high throughput RNA sequencing technology was employed to characterize the transcriptome of radish buds with length about 1.5 mm from two CMS lines possessing the CMS-inducing orf138 gene and corresponding near-isogenic maintainer lines. A total of 67,140 unigenes were functionally annotated. Functional terms for these genes are significantly enriched in 55 Gene Ontology (GO) groups and 323 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The transcriptome detected transcripts for 72 out of a total of 79 protein genes encoded in the chloroplast genome from radish. In contrast, the radish mitochondrial genome contains 34 protein genes, but only 16 protein transcripts were detected from the transcriptome. The transcriptome comparison between CMS and near-isogenic maintainer lines revealed 539 differentially expressed genes (DEGs), indicating that the false positive rate for comparative transcriptome profiling was clearly decreased using two groups of CMS/maintainer lines with different nuclear background. The level of 127 transcripts was increased and 412 transcripts were decreased in the CMS lines. No change in levels of transcripts except CMS-inducing orf138 was identified from the mitochondrial and chloroplast genomes. Some DEGs which would be associated with the CMS, encoding MYB and bHLH transcription factors, pentatricopeptide repeat (PPR) proteins, heat shock transcription factors (HSFs) and heat shock proteins (HSPs), are discussed. The transcriptome dataset and comparative analysis will provide an important resource for further understanding anther development, the CMS mechanism and to improve molecular breeding in radish.

Graphical Abstract

1. Introduction

Plant cytoplasmic male sterility (CMS) is a maternally inherited trait preventing the production of functional pollen but maintaining female fertility. CMS has been observed in many crops and widely applied to produce hybrids with significant heterosis [1]. CMS is usually associated with chimeric open reading frames caused by mitochondrial genome reorganization [2,3,4]. In many cases, pollen fertility of the CMS plants can be restored by corresponding nuclear-encoded fertility restorer genes suggesting that CMS is an ideal model system to study nuclear-cytoplasmic gene interaction [5].
Considerable progress has been made in the molecular dissection of CMS and fertility restoration in plants. Many novel chimeric mitochondrial genes have been identified and are responsible for CMS through a variety of mechanisms (such as urf13 from maize, pcf from petunia, orf79/orfH79 from rice, orf224/orf222 from rapeseed etc.) [6,7,8,9,10,11]. Recently, a new mitochondrial gene, WA352, was found to confer CMS in rice. WA352 interacts with the nuclear-encoded mitochondrial protein OsCOX11 which has a role in hydrogen peroxide degradation. The interaction between WA352 and OsCOX11 triggers an early reactive oxygen species (ROS) burst which promotes cytochrome c release to the cytosol, causing premature tapetal programmed cell death (PCD) [3]. Several fertility restorer (Rf) genes were isolated through map-based cloning strategy. The first cloned Rf gene encodes an aldehyde dehydrogenase [12]. Rf17 and Rf2 in rice encode acyl-carrier protein synthase-like protein and glycine-rich protein, respectively [13,14]. Sugar beet Rf1 encodes the mitochondrial protein bvORF20 that resembles a yeast protease involved in mitochondrial protein quality control and interacts with the sterility-inducing factor preSATP6 [15,16]. All the other Rf genes encode pentatricopeptide repeat (PPR) proteins [8,17,18,19,20,21,22,23].
Radish (Raphanus sativus L.), an annual or biennial herb belonging to the family Brassicaceae, is an economically important root vegetable crop produced throughout the world. The edible taproot of radish is an excellent source of dietary fiber and carbohydrates, essential minerals and organic nutrients for humans [24,25,26]. Three different radish CMS germplasms, Ogura/Kosena, NWB and DCGMS CMSs, have been identified [18,20,27,28]. Orf138/orf125 confers Ogura/Kosena CMS [18,20], and orf463 is likely to be a causative factor for DCGMS CMS [28]. In recent years, next-generation sequencing technology obviously increased the efficiency and speed of gene discovery. By high-throughput cDNA sequencing (RNA-Seq), researchers can obtain almost the entire transcriptional landscape of gene expression in a high throughput and quantitative way. In this study, comparative transcriptome analysis of radish buds from CMS and their maintainer lines was conducted by Illumina sequencing, which might provide assistance for further deep analysis about anther development and CMS in radish.

2. Results and Discussion

2.1. The Floral Bud Size Is Associated with the Microspore Developmental Stage

No morphological difference between the CMS line and the corresponding maintainer line was observed except for the anther. CMS lines have empty yellow anthers whereas maintainer lines have full and dehiscent anthers (Figure 1). No histological differences were observed between the CMS line 9802A1 and the corresponding maintainer line 9802B1 in the anther at the tetrad stage when tetrads were enclosed by a thick callose wall in our previous study [29]. After the tetrad stage, microspores released from tetrads underwent cell expansion and cell wall synthesis in male fertile anthers. In contrast, most microspores from CMS anthers grew slowly and were malformed, and the tapetal cell grew abnormally large with increased vacuoles [29]. The cytoplasm of 9802A1 was introgressed into HYBP-B and YH-B at more than 16 times of backcrossing. The resulting CMS lines, HYBP-A and YH-A, display different morphologies, indicating that they differ in their nuclear genomes. In this work, the stage of pollen abortion in HYBP-A and YH-A is the same with that of 9802A1 suggesting that different nuclear backgrounds could not change the pollen abortion phase from the same CMS cytoplasm. In the four lines, HYBP-A/B and YH-A/B, when the buds reached about 0.8 mm, anthers were at the microspore mother cell stage (Figure 2A). Microspore mother cells were undergoing meiosis when the buds were about 1.1 mm in length (Figure 2B). When the bud length was about 1.5 mm, individual tetrads could be seen (Figure 2C). In HYBP-B and YH-B, the phase of the vacuolate microspore appeared in buds about 2.0 mm in length (Figure 2D), with a single large vacuole displacing the single nucleus to one side of the cell.
Ijms 17 00042 g001 1024
Figure 1. Partially dissected radish flowers. Male fertile plant (A); and male sterile plant (B). Bar = 1 mm.
Figure 1. Partially dissected radish flowers. Male fertile plant (A); and male sterile plant (B). Bar = 1 mm.
Ijms 17 00042 g001
Ijms 17 00042 g002 1024
Figure 2. Development of anthers in the maintainer line HYBP-B (A–D) and the CMS line HYBP-A (E–H). Anther from floral bud at about 0.8 mm long (A,E); Anther from floral bud at about 1.1 mm long (B,F); Anther from floral bud at about 1.5 mm long (C,G); Anther from floral bud at about 2.0 mm long (D,H). Bar = 20 μm.
Figure 2. Development of anthers in the maintainer line HYBP-B (A–D) and the CMS line HYBP-A (E–H). Anther from floral bud at about 0.8 mm long (A,E); Anther from floral bud at about 1.1 mm long (B,F); Anther from floral bud at about 1.5 mm long (C,G); Anther from floral bud at about 2.0 mm long (D,H). Bar = 20 μm.
Ijms 17 00042 g002

2.2. Illumina Sequencing

Since distinct anther defects were histologically observed shortly after the tetrad stage in the CMS lines, we reckoned that the gene expression profiles at the tetrad stage could already change between the CMS lines and maintainer lines. The tetrad stage corresponded to floral buds about 1.5 mm in length from HYBP-A, HYBP-B, YH-A and YH-B. Thereafter, RNA of floral buds about 1.5 mm long was isolated from these four lines, and used for RNA-Seq. After trimming the raw data (93,800,000 reads), a total of 83,932,141 clean reads above the minimum quality threshold of Q20 per basewere selected with an average GC content of 47.7%.

2.3. De Novo Assembly and Functional Classification

High-quality clean reads were de novo assembled into 130,240 contigs with an average length of 570 bp and a N50 length of 821 bp. Sequences that did not extend on either end were defined as unigenes. In total, 67,140 unigenes showed similarity to known genes in at least one of the five public databases including NR, Swiss-Prot, GO, KOG (eukaryotic ortholog groups) and KEGG.
GO annotation including three ontologies (molecular function, cellular component and biological process) is an international gene classification system that provides a dynamically-updated standardized vocabulary for assigning functions of genes and their products [30]. BLAST2GO program was employed to find GO terms for assembled unigenes and 29,031 unigenes were assigned at least one GO term in 55 functional groups. Of these, most were assigned to the biological process ontology (86,785; 41.2%), followed by cellular component (78,749; 37.4%) and molecular function (45,090; 21.4%) ontologies (Figure 3).
The KOGs, based on orthologous relationships between genes, include proteins encoded by seven sequenced eukaryotic genomes [31]. In total, 26,527 unigenes were assigned to the KOG classification. Because some unigenes were annotated with multiple KOG functions, a total of 29,353 functional annotations were obtained. The four largest groups annotated were “posttranslational modification, protein turnover, chaperones” (3171 unigenes), “general function prediction only” (3110 unigenes), “signal transduction mechanisms” (2893 unigenes) and “Translation, ribosomal structure and biogenesis” (2237 unigenes) (Figure 4).
Ijms 17 00042 g003 1024
Figure 3. Gene ontology classifications of unigenes from the radish floral bud transcriptome. The x-axis indicates the sub-categories; the y-axis indicates the number of unigenes in a sub-category.
Figure 3. Gene ontology classifications of unigenes from the radish floral bud transcriptome. The x-axis indicates the sub-categories; the y-axis indicates the number of unigenes in a sub-category.
Ijms 17 00042 g003
Ijms 17 00042 g004 1024
Figure 4. KOG function classification of unigenes from the radish floral bud transcriptome. The x-axis indicates the number of unigenes in a function class; the y-axis indicates the function classes.
Figure 4. KOG function classification of unigenes from the radish floral bud transcriptome. The x-axis indicates the number of unigenes in a function class; the y-axis indicates the function classes.
Ijms 17 00042 g004
The KEGG pathway records molecular interaction networks such as pathways and complexes. Sequence similarity search against KEGG GENES identified 5706 unigenes assigned to 323 KEGG pathways. The result showed that the three largest pathway categories were “carbohydrate metabolism” (803 members), “translation” (595 members) and “signal transduction” (563 members) suggesting that these pathways could be crucial for radish anther development (Figure 5).
Ijms 17 00042 g005 1024
Figure 5. Classification based on categories of KEGG pathways. Cellular Processes (A); Environmental Information Processing (B); Genetic Information Processing (C); and Metabolism (D). The x-axis indicates the number of unigenes in a function class; the y-axis indicates the function classes.
Figure 5. Classification based on categories of KEGG pathways. Cellular Processes (A); Environmental Information Processing (B); Genetic Information Processing (C); and Metabolism (D). The x-axis indicates the number of unigenes in a function class; the y-axis indicates the function classes.
Ijms 17 00042 g005

2.4. Comparison of Gene Expression Levels between CMS and Maintainer Lines

RNA-Seq gene expression levels were measured using reads per kilobase of exon model per million mapped sequence reads (RPKM) [32]. RPKM values for each gene were compared between the CMS and corresponding maintainer lines to screen for DEGs. We identified 3843 DEGs (2487 up-regulated and 1356 down-regulated genes) in HYBP-A floral buds compared with HYBP-B and 2041 DEGs (676 up-regulated and 1365 down-regulated genes) in YH-A floral buds compared with YH-B. Of these, 539 DEGs exhibited the same direction of expression change in the two CMS lines relative to their corresponding maintainer lines. The transcription level of a total of 127 transcripts was increased and 412 transcripts was decreased in both HYBP-A and YH-A compared with their corresponding maintainer lines (Table S2). These data suggest that combination analysis of two groups of CMS/maintainer lines with different nuclear background was useful to decrease false positives. GO annotation for these DEGs showed that 268 genes were categorized into 50 functional groups with 1629 functional terms indicating that some genes were annotated with multiple terms (Table S3). The DEGs were assigned into nine subcategories in the “molecular function” category, and the two most abundant terms were ‘binding” (149 members) and “catalytic activity” (122 members) (Figure 6). Among “cellular component”, “cell” (199 members), “cell part” (199 members) and “organelle” (108 members) were the three most strongly represented terms (Figure 6). For “biological process”, the three predominant terms for these DEGs were “metabolic process” (141 members), “cellular process” (135 members), “response to stimulus” (83 members) followed by “biological regulation” (50 members) (Figure 6).
There are a total of 79 protein-encoding genes from the radish chloroplast (cp) genome [33], and we retrieved transcripts for 72 cp protein-encoding genes from radish floral buds at this work (transcripts for cp protein genes, psaJ, petG, ndhK, psaB, rps14, psbC and psbD, were not detected). In contrast, mitochondrial (mt) genome from radish contains 34 protein genes [34], but only 16 protein transcripts (atp1, atp4, atp6, atp9, ccmFC, cob, cox2, cox3, nad4, nad4L, nad5, nad9, orf138, rpL16, rps3 and rps4) were detected from the transcriptome. No difference in levels of transcripts between CMS and maintainer lines was found from these annotated mt and cp transcripts in the present work, except CMS-inducing mt orf138, which was only expressed in the CMS lines.
Ijms 17 00042 g006 1024
Figure 6. Gene Ontology classification of down- and up-regulated unigenes in floral buds in CMS lines compared with maintainer lines. The x-axis indicates the sub-categories; the y-axis indicates the number of unigenes in a sub-category.
Figure 6. Gene Ontology classification of down- and up-regulated unigenes in floral buds in CMS lines compared with maintainer lines. The x-axis indicates the sub-categories; the y-axis indicates the number of unigenes in a sub-category.
Ijms 17 00042 g006
Transcription factors (TFs) play important roles in complex biological processes under a wide range of environmental signals by regulating gene transcription through binding to specific DNA sequences in the promoters of multiple target genes [35]. At least 608 transcription factors (TFs) from Arabidopsis divided into 34 families were identified in the male gametophyte [36]. Here, 25 TF genes from 11 gene families were down regulated in the radish CMS lines. The tapetum directly surrounding the developing male gametophyte plays a crucial role in the microspore development and maturation by supplying necessary nutrients and structural components. In both rice and Arabidopsis at least two bHLH factors (UDT1 and TDR, and DYT1 and AMS, respectively) were shown to be important in tapetum development [37,38,39,40]. The timing of tapetal PCD is crucial for microspore development. T-DNA knockout of Arabidopsis MYB80 (encoding an MYB transcription factor) triggers premature tapetal PCD resulting in male sterility, indicating that MYB80 is important for tapetal PCD at the proper developmental stage [41]. In the present study, three bHLH and three MYB genes were found to be down-regulated in CMS lines. Since functions of their corresponding Arabidopsis genes (AT1G26610, AT3G19580, AT5G04340, AT3G50060, AT5G52660, and AT4G37260) are still unkown, it is worth studying whether these genes are involved in radish CMS.
To date, two C2H2 zinc-finger transcription factors (DAZ1 and DAZ2) and two WRKY transcription factors (WRKY34 and WRKY2) in Arabidopsis are responsible for male gametogenesis. DAZ1 and DAZ2 promote G2- to M-phase transition in pollen [42]. WRKY34 and WRKY2 are required for pollen development at early stage [43]. Since the abortive phenotype for our radish CMS is cytologically observed shortly after the tetrad stage, it remains to be determined whether three DEGs encoding a C2H2 zinc-finger transcription factor and two WRKYs are involved in early events of the radish CMS.
Heat shock transcription factors (HSFs) and heat shock proteins (HSPs) play key roles in signal transduction in the heat stress response. Several studies suggest that some HSFs and HSPs are related to gametophyte development in the absence of heat stress [44,45,46]. In current work, five HSF and 19 HSP genes were down-regulated in radish CMS lines. Especially, AtHSP23.6 is targeted to mitochondria [47], and the maize homolog hsp81 is strongly expressed at the pre-meiotic and meiotic stages of microspore development without heat shock [45].
Many fertility restorer genes for the plant CMSs gene were identified as belonging to the PPR family [8,17,18,19,20,21,22]. In this work, three PPR genes were downregulated in the CMS lines. Their orthologs (GenBank accession numbers AT2G20540, AT1G33350 and AT1G71060) from Arabidopsis were predicted to target to the mitochondrion. Interestingly, one ortholog MEF21 from AT2G20540 is involved in the C-to-U transcript editing of the mitochondrial gene cox3 [48]. These PPR DEGs would be very valuable for further study of their roles in the radish CMS.

2.5. Validation of the Expression Difference of DEGs

In order to validate the expression profiling by Illumina sequencing, the expression levels of 10 genes, which showed differential expression between the CMS line HYBP-A and maintainer line HYBP-B, were further analyzed by qRT-PCR. The results showed that the direction of expression change of all the selected genes based on qRT-PCR agreed with those detected by Illumina sequencing (Figure S1). However, the fold changes in expression levels of some genes were different depending on the detection method.

3. Experimental Section

3.1. Plant Materials and RNA Extraction

The male sterile line 9802A1 possesses the CMS-inducing mitochondrial gene orf138 [49]. CMS of 9802A1 was introduced in two maintainer lines (HYBP-B and YH-B) by a repeated backcross program (16 generations) using the maintainer lines as recurrent parents. Thus, the resulting two lines, HYBP-A and YH-A, have the same nuclear genomes as HYBP-B and YH-B, respectively. HYBP-B has round red-skinned root and YH-B is an elongated white variety, indicating a different genetic background between the two lines.

3.2. Histological Analysis

Fresh floral buds were stored overnight in a formalin-acetic acid-alcohol (FAA) fixative solution. Fixed buds were dehydrated through a series of differently graded ethanol ranging from 50% to 100%. After the length was rapidly measured, dehydrated buds were embedded into Technovit 7100 (Heraeus Kulzer, Wehrheim, Germany) and sectioned with glass knives in a Leica ultracut R ultramicrotome (Leica Microsystems, Wetzlar, Germany). Serial sections (1.5 μm in thickness) of anther tissues were heat fixed to glass slides and stained with 2% toluidine blue. Sectioned anthers were observed with an Olympus BX61 microscope equipped with a colour charge-coupled device (CCD) camera.

3.3. RNA Extraction and RNA-Seq

Floral buds about 1.5 mm in length from HYBP-A, HYBP-B, YH-A and YH-B, were sampled, immediately frozen in liquid nitrogen and kept at −80 °C for RNA extraction. Total RNA was isolated using an E.Z.N.A. Plant RNA Kit (OMEGA Bio-Tek, Norcross, GA, USA) according to the manufacturer’s protocol. RNA samples were treated with DNase I at 37 °C for 30 min to remove DNA. A fragmentation buffer was added to break the mRNA-enriched RNAs into short pieces. With these short fragments as templates, first strands of cDNAs were synthesized using random hexamer primers. The second-strand cDNA was generated using the Second Strand Synthesis Reaction Buffer and the Second Strand Synthesis Enzyme Mix from NEBNext® Ultra™ RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). Short fragments were purified with AMPure XP Beads (Agencourt, Beverly, MA, USA) for end repair and tailing A. The resulting A-tailed fragments were mixed with NEBNext Adaptor and Blunt/TA Ligase Master Mix for adaptor ligation. Finally, the sequencing libraries were constructed following PCR amplification with these adaptor-ligated fragments. Paired-end sequencing was then performed using the Illumina sequencing platform (HiSeq™ 2500) according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). Illumina sequencing was performed at Shanghai Hanyu Bio-Tech Co., Ltd., in Shanghai (China).

3.4. De Novo Assembly and Database Search

Reads from each library were assembled together. After trimming adapter sequences and filtering low-quality reads, the clean reads were assembled using the software package Trinity consisting of three components: Inchworm, Chrysalis and Butterfly [50]. For radish transcriptome unigene annotation, the assembled unigenes were searched against five public databases using NCBI BLASTX: the Swiss-Prot protein, Gene Ontology (GO), NCBI non-redundant protein sequences (NR), the Kyoto Encyclopedia of Genes and Genomes (KEGG) and eukaryotic ortholog groups (KOG), with an E-value threshold of 10−5.

3.5. Differential Gene Expression

Gene expression levels were measured during the RNA-Seq analysis in RPKM [32]. The DEseq software was used to identify differentially expressed genes in pair-wise comparisons and the results of all statistical tests were corrected for multiple testing using the Benjamini-Hochberg false discovery rate (FDR < 0.001) [51]. Sequences were deemed to be significantly differentially expressed if the adjusted p-value obtained by this method (i.e., Q-value) was <0.001 and there was at least a two-fold change (>1 or <−1 in log2-ratio value which was calculated dividing the RPKM value of the CMS line by the RPKM value of the maintainer line). The data was deposited at the Gene Expression Omnibus (GEO) database under accession number GSE65142.
GO function enrichment analyses were used to identify markedly enriched function categories, in differentially expressed genes compared with the whole genome background. The formula for the calculation is as follows:
p = 1 − ∑ i = 0 m − 1 ( M i ) ( N − M n − i ) ( N n )
where, N is the number of all genes with GO annotation, n is the number of DEGs in N, M is the number of all genes annotated to specific GO functional category, and m is the number of DEGs in M.

3.6. Real-Time Quantitative PCR (qRT-PCR) Analysis

qRT-PCR was performed using an optical 96-well plate with an iQ5 multicolor real-time PCR system (Bio-RAD, Hercules, CA, USA). Each reaction contained 1.0 μL of cDNA template from the reverse-transcribed reaction mentioned above, 10 nM of gene-specific primers, and 10 μL of iTaq™ Universal SYBR Green supermix (Bio-RAD) in a final volume of 20 μL. The primer sequence of DEGs and the control gene were listed in Table S1. The thermal cycle used was as follows: 95 °C for 30 s, followed by 40 cycles of 95 °C for 15 s, 60 °C for 30 s. qRT-PCR was performed in triplicate for each sample. The comparative threshold (Ct) cycle method was used for determination of relative transcript levels normalized to the actin gene.

4. Conclusions

A global view of differential expression of genes from floral buds between radish CMS and corresponding maintainer lines was analyzed in this work. These results provide support for the study of the molecular mechanism of radish CMS.

Supplementary Materials

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

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (31470412 and 31100231). We kindly thank members from Shanghai Hanyu Bio-Tech Co., Ltd. for their assistance with original data processing and related bioinformatic analysis. We thank technician Yibin Fu for his excellent fieldwork.

Author Contributions

Shiyong Mei and Zhiwei Wang designed the research; Shiyong Mei, Touming Liu performed the experiments; Shiyong Mei, Touming Liu and Zhiwei Wang analyzed data; Shiyong Mei, Touming Liu and Zhiwei Wang wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, L.; Liu, Y.G. Male sterility and fertility restoration in crops. Annu. Rev. Plant Biol. 2014, 65, 579–606. [Google Scholar] [CrossRef] [PubMed]
  2. Ivanov, M.; Dymshits, G. Cytoplasmic male sterility and restoration of pollen fertility in higher plants. Russ. J. Genet. 2007, 43, 354–368. [Google Scholar] [CrossRef]
  3. Luo, D.; Xu, H.; Liu, Z.; Guo, J.; Li, H.; Chen, L.; Fang, C.; Zhang, Q.; Bai, M.; Yao, N. A detrimental mitochondrial-nuclear interaction causes cytoplasmic male sterility in rice. Nat. Genet. 2013, 45, 573–577. [Google Scholar] [CrossRef] [PubMed]
  4. Schnable, P.S.; Wise, R.P. The molecular basis of cytoplasmic male sterility and fertility restoration. Trends Plant Sci. 1998, 3, 175–180. [Google Scholar] [CrossRef]
  5. Hanson, M.R.; Bentolila, S. Interactions of mitochondrial and nuclear genes that affect male gametophyte development. Plant Cell 2004, 16, S154–S169. [Google Scholar] [CrossRef] [PubMed]
  6. L’Homme, Y.; Stahl, R.J.; Li, X.Q.; Hameed, A.; Brown, G.G. Brassica nap cytoplasmic male sterility is associated with expression of a mtDNA region containing a chimeric gene similar to the pol CMS-associated orf224 gene. Curr. Genet. 1997, 31, 325–335. [Google Scholar] [CrossRef] [PubMed]
  7. Nivison, H.T.; Sutton, C.A.; Wilson, R.K.; Hanson, M.R. Sequencing, processing, and localization of the petunia CMS-associated mitochondrial protein. Plant J. 1994, 5, 613–623. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, Z.; Zou, Y.; Li, X.; Zhang, Q.; Chen, L.; Wu, H.; Su, D.; Chen, Y.; Guo, J.; Luo, D.; et al. Cytoplasmic male sterility of rice with Boro II cytoplasm is caused by a cytotoxic peptide and is restored by two related PPR motif genes via distinct modes of mRNA silencing. Plant Cell 2006, 18, 676–687. [Google Scholar] [CrossRef] [PubMed]
  9. Wise, R.; Fliss, A.; Pring, D.; Gengenbach, B. urf13-T of T cytoplasm maize mitochondria encodes a 13 kD polypeptide. Plant Mol. Biol. 1987, 9, 121–126. [Google Scholar] [CrossRef] [PubMed]
  10. Wolf-Litman, O.; Soferman, O.; Tabib, Y.; Izhar, S. Interaction of the mitochondrial S-Pcf locus for cytoplasmic male sterility in Petunia with multiple fertility-restoration genes in somatic hybrid plants. Theor. Appl. Genet. 1992, 84, 829–834. [Google Scholar] [CrossRef] [PubMed]
  11. Wang, K.; Gao, F.; Ji, Y.; Liu, Y.; Dan, Z.; Yang, P.; Zhu, Y.; Li, S. ORFH79 impairs mitochondrial function via interaction with a subunit of electron transport chain complex III in Honglian cytoplasmic male sterile rice. New Phytol. 2013, 198, 408–418. [Google Scholar] [CrossRef] [PubMed]
  12. Cui, X.; Wise, R.P.; Schnable, P.S. The rf2 nuclear restorer gene of male-sterile T-cytoplasm maize. Science 1996, 272, 1334–1335. [Google Scholar] [CrossRef] [PubMed]
  13. Fujii, S.; Toriyama, K. Suppressed expression of RETROGRADE-REGULATED MALE STERILITY restores pollen fertility in cytoplasmic male sterile rice plants. Proc. Natl. Acad. Sci. USA 2009, 106, 9513–9518. [Google Scholar] [CrossRef] [PubMed]
  14. Itabashi, E.; Iwata, N.; Fujii, S.; Kazama, T.; Toriyama, K. The fertility restorer gene, Rf2, for Lead Rice-type cytoplasmic male sterility of rice encodes a mitochondrial glycine-rich protein. Plant J. 2011, 65, 359–367. [Google Scholar] [CrossRef] [PubMed]
  15. Kitazaki, K.; Arakawa, T.; Matsunaga, M.; Yui-Kurino, R.; Matsuhira, H.; Mikami, T.; Kubo, T. Post-translational mechanisms are associated with fertility restoration of cytoplasmic male sterility in sugar beet. Plant J. 2015, 83, 290–299. [Google Scholar] [CrossRef] [PubMed]
  16. Matsuhira, H.; Kagami, H.; Kurata, M.; Kitazaki, K.; Matsunaga, M.; Hamaguchi, Y.; Hagihara, E.; Ueda, M.; Harada, M.; Muramatsu, A. Unusual and typical features of a novel restorer-of-fertility gene of sugar beet (Beta vulgaris L.). Genetics 2012, 192, 1347–1358. [Google Scholar] [CrossRef] [PubMed]
  17. Bentolila, S.; Alfonso, A.A.; Hanson, M.R. A pentatricopeptide repeat-containing gene restores fertility to cytoplasmic male-sterile plants. Proc. Natl. Acad. Sci. USA 2002, 99, 10887–10892. [Google Scholar] [CrossRef] [PubMed]
  18. Brown, G.G.; Formanová, N.; Jin, H.; Wargachuk, R.; Dendy, C.; Patil, P.; Laforest, M.; Zhang, J.; Cheung, W.Y.; Landry, B.S. The radish Rfo restorer gene of Ogura cytoplasmic male sterility encodes a protein with multiple pentatricopeptide repeats. Plant J. 2003, 35, 262–272. [Google Scholar] [CrossRef] [PubMed]
  19. Desloire, S.; Gherbi, H.; Laloui, W.; Marhadour, S.; Clouet, V.; Cattolico, L.; Falentin, C.; Giancola, S.; Renard, M.; Budar, F. Identification of the fertility restoration locus, Rfo, in radish, as a member of the pentatricopeptide-repeat protein family. EMBO Rep. 2003, 4, 588–594. [Google Scholar] [CrossRef] [PubMed]
  20. Koizuka, N.; Imai, R.; Fujimoto, H.; Hayakawa, T.; Kimura, Y.; Kohno-Murase, J.; Sakai, T.; Kawasaki, S.; Imamura, J. Genetic characterization of a pentatricopeptide repeat protein gene, orf687, that restores fertility in the cytoplasmic male-sterile Kosena radish. Plant J. 2003, 34, 407–415. [Google Scholar] [CrossRef] [PubMed]
  21. Hu, J.; Wang, K.; Huang, W.; Liu, G.; Gao, Y.; Wang, J.; Huang, Q.; Ji, Y.; Qin, X.; Wan, L. The rice pentatricopeptide repeat protein RF5 restores fertility in Hong-Lian cytoplasmic male-sterile lines via a complex with the glycine-rich protein GRP162. Plant Cell 2012, 24, 109–122. [Google Scholar] [CrossRef] [PubMed]
  22. Tang, H.; Luo, D.; Zhou, D.; Zhang, Q.; Tian, D.; Zheng, X.; Chen, L.; Liu, Y.G. The rice restorer Rf4 for wild-abortive cytoplasmic male sterility encodes a mitochondrial-localized PPR protein that functions in reduction of WA352 transcripts. Mol. Plant Pathol. 2014, 7, 1497–1500. [Google Scholar] [CrossRef] [PubMed]
  23. Huang, W.; Yu, C.; Hu, J.; Wang, L.; Dan, Z.; Zhou, W.; He, C.; Zeng, Y.; Yao, G.; Qi, J. Pentatricopeptide-repeat family protein RF6 functions with hexokinase 6 to rescue rice cytoplasmic male sterility. Proc. Natl. Acad. Sci. USA 2015, 112, 14984–14989. [Google Scholar] [CrossRef] [PubMed]
  24. Curtis, I.S. Genetic engineering of radish: Current achievements and future goals. Plant Cell Rep. 2011, 30, 733–744. [Google Scholar] [CrossRef] [PubMed]
  25. Khanum, F.; Swamy, M.S.; Krishna, K.S.; Santhanam, K.; Viswanathan, K. Dietary fiber content of commonly fresh and cooked vegetables consumed in India. Plant Foods Hum. Nutr. 2000, 55, 207–218. [Google Scholar] [CrossRef] [PubMed]
  26. Vardhini, B.V.; Sujatha, E.; Rao, S.; Ram, S. Studies on the effect of brassinosteroids on the qualitative changes in the storage roots of radish. Bulg. J. Agric. Sci. 2012, 18, 63–69. [Google Scholar]
  27. Nahm, S.H.; Lee, H.J.; Lee, S.W.; Joo, G.Y.; Harn, C.H.; Yang, S.G.; Min, B.W. Development of a molecular marker specific to a novel CMS line in radish (Raphanus sativus L.). Theor. Appl. Genet. 2005, 111, 1191–1200. [Google Scholar] [CrossRef] [PubMed]
  28. Park, J.Y.; Lee, Y.P.; Lee, J.; Choi, B.S.; Kim, S.; Yang, T.J. Complete mitochondrial genome sequence and identification of a candidate gene responsible for cytoplasmic male sterility in radish (Raphanus sativus L.) containing DCGMS cytoplasm. Theor. Appl. Genet. 2013, 126, 1763–1774. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, Z.W.; Wang, C.D.; Gao, L.; Mei, S.Y.; Zhou, Y.; Xiang, C.P.; Wang, T. Heterozygous alleles restore male fertility to cytoplasmic male-sterile radish (Raphanus sativus L.): A case of overdominance. J. Exp. Bot. 2013, 64, 2041–2048. [Google Scholar] [CrossRef] [PubMed]
  30. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T. Gene Ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed]
  31. Tatusov, R.L.; Fedorova, N.D.; Jackson, J.D.; Jacobs, A.R.; Kiryutin, B.; Koonin, E.V.; Krylov, D.M.; Mazumder, R.; Mekhedov, S.L.; Nikolskaya, A.N.; et al. The COG database: An updated version includes eukaryotes. BMC Bioinform. 2003, 4, 41. [Google Scholar] [CrossRef] [PubMed]
  32. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [Google Scholar] [CrossRef] [PubMed]
  33. Jeong, Y.M.; Chung, W.H.; Mun, J.H.; Kim, N.; Yu, H.J. De novo assembly and characterization of the complete chloroplast genome of radish (Raphanus sativus L.). Gene 2014, 551, 39–48. [Google Scholar] [CrossRef] [PubMed]
  34. Tanaka, Y.; Tsuda, M.; Yasumoto, K.; Yamagishi, H.; Terachi, T. A complete mitochondrial genome sequence of Ogura-type male-sterile cytoplasm and its comparative analysis with that of normal cytoplasm in radish (Raphanus sativus L.). BMC Genom. 2012, 13, 352. [Google Scholar] [CrossRef] [PubMed]
  35. Franco-Zorrilla, J.M.; López-Vidriero, I.; Carrasco, J.L.; Godoy, M.; Vera, P.; Solano, R. DNA-binding specificities of plant transcription factors and their potential to define target genes. Proc. Natl. Acad. Sci. USA 2014, 111, 2367–2372. [Google Scholar] [CrossRef] [PubMed]
  36. Honys, D.; Twell, D. Transcriptome analysis of haploid male gametophyte development in Arabidopsis. Genome Biol. 2004, 5, R85. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Jung, K.H.; Han, M.J.; Lee, Y.S.; Kim, Y.W.; Hwang, I.; Kim, M.J.; Kim, Y.K.; Nahm, B.H.; An, G. Rice Undeveloped Tapetum1 is a major regulator of early tapetum development. Plant Cell 2005, 17, 2705–2722. [Google Scholar] [CrossRef] [PubMed]
  38. Sorensen, A.M.; Kröber, S.; Unte, U.S.; Huijser, P.; Dekker, K.; Saedler, H. The Arabidopsis ABORTED MICROSPORES (AMS) gene encodes a MYC class transcription factor. Plant J. 2003, 33, 413–423. [Google Scholar] [CrossRef] [PubMed]
  39. Zhang, W.; Sun, Y.; Timofejeva, L.; Chen, C.; Grossniklaus, U.; Ma, H. Regulation of Arabidopsis tapetum development and function by DYSFUNCTIONAL TAPETUM1 (DYT1) encoding a putative bHLH transcription factor. Development 2006, 133, 3085–3095. [Google Scholar] [CrossRef] [PubMed]
  40. Li, N.; Zhang, D.S.; Liu, H.S.; Yin, C.S.; Li, X.X.; Liang, W.Q.; Yuan, Z.; Xu, B.; Chu, H.W.; Wang, J.; et al. The rice tapetum degeneration retardation gene is required for tapetum degradation and anther development. Plant Cell 2006, 18, 2999–3014. [Google Scholar] [CrossRef] [PubMed]
  41. Phan, H.A.; Iacuone, S.; Li, S.F.; Parish, R.W. The MYB80 transcription factor is required for pollen development and the regulation of tapetal programmed cell death in Arabidopsis thaliana. Plant Cell 2011, 23, 2209–2224. [Google Scholar] [CrossRef] [PubMed]
  42. Borg, M.; Rutley, N.; Kagale, S.; Hamamura, Y.; Gherghinoiu, M.; Kumar, S.; Sari, U.; Esparza-Franco, M.A.; Sakamoto, W.; Rozwadowski, K. An EAR-dependent regulatory module promotes male germ cell division and sperm fertility in Arabidopsis. Plant Cell 2014, 26, 2098–2113. [Google Scholar] [CrossRef] [PubMed]
  43. Guan, Y.F.; Meng, X.Z.; Khanna, R.; LaMontagne, E.; Liu, Y.D.; Zhang, S.Q. Phosphorylation of a WRKY transcription factor by MAPKs is required for pollen development and function in Arabidopsis. PLoS Genet. 2014, 10, e1004384. [Google Scholar] [CrossRef] [PubMed]
  44. Wunderlich, M.; Groß-Hardt, R.; Schöffl, F. Heat shock factor HSFB2a involved in gametophyte development of Arabidopsis thaliana and its expression is controlled by a heat-inducible long non-coding antisense RNA. Plant Mol. Biol. 2014, 85, 541–550. [Google Scholar] [CrossRef] [PubMed]
  45. Marrs, K.A.; Casey, E.S.; Capitant, S.A.; Bouchard, R.A.; Dietrich, P.S.; Mettler, I.J.; Sinibaldi, R.M. Characterization of two maize HSP90 heat shock protein genes: Expression during heat shock, embryogenesis, and pollen development. Dev. Genet. 1993, 14, 27–41. [Google Scholar] [CrossRef] [PubMed]
  46. Giorno, F.; Wolters-Arts, M.; Grillo, S.; Scharf, K.D.; Vriezen, W.H.; Mariani, C. Developmental and heat stress-regulated expression of HsfA2 and small heat shock proteins in tomato anthers. J. Exp. Bot. 2010, 61, 453–462. [Google Scholar] [CrossRef] [PubMed]
  47. Bonen, L.; Waters, E.R.; Nguyen, S.L.; Eskandar, R.; Behan, J.; Sanders-Reed, Z. The recent evolution of a pseudogene: Diversity and divergence of a mitochondria-localized small heat shock protein in Arabidopsis thaliana. Genome 2008, 51, 177–186. [Google Scholar] [CrossRef] [PubMed]
  48. Takenaka, M.; Verbitskiy, D.; Zehrmann, A.; Brennicke, A. Reverse genetic screening identifies five E-class PPR proteins involved in RNA editing in mitochondria of Arabidopsis thaliana. J. Biol. Chem. 2010, 285, 27122–27129. [Google Scholar] [CrossRef] [PubMed]
  49. Wang, Z.W.; Zhang, Y.J.; Xiang, C.P.; Mei, S.Y.; Zhou, Y.; Chen, G.P.; Wang, T. A new fertility restorer locus linked closely to the Rfo locus for cytoplasmic male sterility in radish. Theor. Appl. Genet. 2008, 117, 313–320. [Google Scholar] [CrossRef] [PubMed]
  50. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 2013, 8, 1494–1512. [Google Scholar] [CrossRef] [PubMed]
  51. Wang, L.; Feng, Z.; Wang, X.; Wang, X.; Zhang, X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 2010, 26, 136–138. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Mei, S.; Liu, T.; Wang, Z. Comparative Transcriptome Profile of the Cytoplasmic Male Sterile and Fertile Floral Buds of Radish (Raphanus sativus L.). Int. J. Mol. Sci. 2016, 17, 42. https://doi.org/10.3390/ijms17010042

AMA Style

Mei S, Liu T, Wang Z. Comparative Transcriptome Profile of the Cytoplasmic Male Sterile and Fertile Floral Buds of Radish (Raphanus sativus L.). International Journal of Molecular Sciences. 2016; 17(1):42. https://doi.org/10.3390/ijms17010042

Chicago/Turabian Style

Mei, Shiyong, Touming Liu, and Zhiwei Wang. 2016. "Comparative Transcriptome Profile of the Cytoplasmic Male Sterile and Fertile Floral Buds of Radish (Raphanus sativus L.)" International Journal of Molecular Sciences 17, no. 1: 42. https://doi.org/10.3390/ijms17010042

APA Style

Mei, S., Liu, T., & Wang, Z. (2016). Comparative Transcriptome Profile of the Cytoplasmic Male Sterile and Fertile Floral Buds of Radish (Raphanus sativus L.). International Journal of Molecular Sciences, 17(1), 42. https://doi.org/10.3390/ijms17010042

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