Transcriptome profile analysis of young floral buds of fertile and sterile plants from the self-pollinated offspring of the hybrid between novel restorer line NR1 and Nsa CMS line in Brassica napus
- Xiaohong Yan†1,
- Caihua Dong†1,
- Jingyin Yu†1,
- Wanghui Liu1,
- Chenghong Jiang1,
- Jia Liu1,
- Qiong Hu1,
- Xiaoping Fang1 and
- Wenhui Wei1Email author
© Yan et al.; licensee BioMed Central Ltd. 2013
Received: 18 September 2012
Accepted: 2 January 2013
Published: 16 January 2013
The fertile and sterile plants were derived from the self-pollinated offspring of the F1 hybrid between the novel restorer line NR1 and the Nsa CMS line in Brassica napus. To elucidate gene expression and regulation caused by the A and C subgenomes of B. napus, as well as the alien chromosome and cytoplasm from Sinapis arvensis during the development of young floral buds, we performed a genome-wide high-throughput transcriptomic sequencing for young floral buds of sterile and fertile plants.
In this study, equal amounts of total RNAs taken from young floral buds of sterile and fertile plants were sequenced using the Illumina/Solexa platform. After filtered out low quality data, a total of 2,760,574 and 2,714,441 clean tags were remained in the two libraries, from which 242,163 (Ste) and 253,507 (Fer) distinct tags were obtained. All distinct sequencing tags were annotated using all possible CATG+17-nt sequences of the genome and transcriptome of Brassica rapa and those of Brassica oleracea as the reference sequences, respectively. In total, 3231 genes of B. rapa and 3371 genes of B. oleracea were detected with significant differential expression levels. GO and pathway-based analyses were performed to determine and further to understand the biological functions of those differentially expressed genes (DEGs). In addition, there were 1089 specially expressed unknown tags in Fer, which were neither mapped to B. oleracea nor to B. rapa, and these unique tags were presumed to arise basically from the added alien chromosome of S. arvensis. Fifteen genes were randomly selected and their expression levels were confirmed by quantitative RT-PCR, and fourteen of them showed consistent expression patterns with the digital gene expression (DGE) data.
A number of genes were differentially expressed between the young floral buds of sterile and fertile plants. Some of these genes may be candidates for future research on CMS in Nsa line, fertility restoration and improved agronomic traits in NR1 line. Further study of the unknown tags which were specifically expressed in Fer will help to explore desirable agronomic traits from wild species.
The novel allo-cytoplasmic male sterility (CMS) system, Nsa CMS line, and the corresponding restorer system, NR1 line, have been successfully developed from somatic hybrids between Brassica napus (oilseed rape) and its wild relative Sinapis arvensis (Yeyou 18, Xinjiang wild mustard from northwestern China) by fusing mesophyll protoplasts [1, 2]. Yeyou 18, a Chinese wild population cataloged into S. arvensis based on genetic analyses , possesses valuable agricultural traits such as enhanced resistance to Sclerotinia sclerotiorum, Leptosphaeria maculans and insects, greater tolerance to low temperatures and drought, low contents of erucic acid and glucosinolates , as well as a low incidence of pod shattering [5, 6]. The Nsa CMS line contains the S. arvensis cytoplasm, which is essentially different from other rapeseed CMS systems such as ogu, nap, pol, tour and hau, based on their origins and molecular characterization . And the cytoplasmic sterile gene in the Nsa CMS line was likely derived from the S. arvensis parent . Furthermore, the Nsa CMS line is more stable to temperature changes compared to pol and nap. NR1, as a B. napus S. arvensis disomic alien addition line, carries one pair of homologous chromosomes from S. arvensis and 19 chromosome pairs from B. napus, and displays important agricultural characters which arise from the alien chromosomes, such as fertility restoration ability to Nsa CMS line, low erucic acid and low glucosinolate contents, S. sclerotiorum resistance and pod shattering resistance .
Fertile and sterile plants were derived from the self-pollinated offspring of the F1 hybrid between the novel restorer line NR1 and the Nsa CMS line. Because NR1 contains one S. arvensis homologous chromosome pair, on which the restorer genes reside, F1 hybrids from NR1 crossed to Nsa CMS line are monosomic . The fertility segregation was observed in self-pollinated plants of F1 hybrid because of the loss of added chromosome, producing fertile and sterile plants, which possess the identical cytoplasmic genetic background arising from Nsa CMS line and similar nuclear genetic background arising from B. napus, except one or two members of the added S. arvensis alien chromosome pair in fertile plants. Floral morphology of fertile plants are normal, whereas sterile plants have stamens reduced in size, abnormal anthers and no pollen produced.
To elucidate gene expression and regulation caused by the A and C subgenomes, alien chromosome and cytoplasm from S. arvensis during the development of young floral bud, especially stamens, we performed a genome-wide high-throughput transcriptomic sequencing for young floral buds of sterile and fertile plants. The transcriptome is the complete set and quantity of transcripts in a cell at a specific developmental stage or under a physiological condition, providing information on gene expression, gene regulation, and amino acid content of proteins . Therefore, transcriptome analysis is essential to interpret the functional elements of the genome and reveals the molecular constituents of cells and tissues. Because of the deep coverage and single base-pair resolution provided by the next-generation sequencing instrument, digital gene-expression (DGE), driven by Solexa/Illumina technology, is an efficient method to analyze transcriptome data. Base on genome-wide expression profiles by sequencing, DGE is able to identify, quantify, and annotate expressed genes on the whole genome level without prior sequence knowledge, opening doors to higher confidence target discovery and pathway studies. This technique has also been widely used in plant research. DGE analysis using Solexa sequencing was performed to identify candidate genes encoding enzymes responsible for the Siraitia grosvenorii triterpene biosynthesis . High-throughput tag-sequencing analysis based on the Solexa Genome Analyzer platform was applied to analyze the gene expression profiling of cucumber plant and revealed the comprehensive mechanisms of waterlogging-responsive transcription . Using the Solexa sequencing system, the transcriptomes were compared between seedlings of two soybean varieties to find genes associated with nitrogen use efficiency . Early developing cotton fiber was analyzed by deep-sequencing, and differential expression patterns of genes in a fuzzless/lintless mutant were revealed . DGE signatures were also used to study maize development, and the results from that study provided a basis for the analysis of short-read expression data and resolved specific expression signatures that will help define mechanisms of action of the maize RA3 gene . In addition, Solexa / Illumina technology was used to analyze gene expression during female flower development  and gene expression of Sinapis alba leaves under drought stress and rewatering growth conditions . Overall, the DGE approach has provided more valuable tools for qualitative and quantitative gene expression analysis than the previous micro array-based assays.
Currently, no progress has been made in identifying the genes that specify CMS and fertility restoration in Nsa CMS and NR1 system, and little is understood about how expression of these genes influences the physiological, developmental or biochemical processes such that pollen formation is disrupted. In this study, the transcriptomes of young floral buds were compared between fertile and sterile plants using the Solexa sequencing system. By investigating the expression of genes related to young floral bud development, a number of candidate genes that are important in this process were identified. Here, to our knowledge, this study is the first to characterize genome-wide comparative analysis of gene expression in young floral buds between fertile and sterile plants. And these sequencing datasets will serve as a valuable resource for novel gene discovery from Yeyou 18, the Chinese wild population of S. arvensis, laying the foundation for elucidating the mechanisms of cytoplasmic male sterility, fertility restoration and improved resistance.
Selection of fertile and sterile plants
CMS-related abnormalities may involve different stages during the formation of male organs and pollen. In many CMS types, the anther development or pollen maturation are impaired. The defect can appear before, during or after meiosis of pollen mother cells resulting in gametophytic or sporophytic CMS types . Although the morphology of normal anther and flower development had been characterized in this study, no detailed morphological comparison of the development of fertile versus sterile anthers was available. To address this, serial sections of fertile and CMS floral buds, representing different development stages of anther, would be compared in further work.
Solexa sequencing tags and the distribution of distinct clean tags in Ste and Fer libraries
total raw data
tags containing N
tag copy number =1
total clean tag
total distinct tag
Annotation analysis of the distinct tag
Annotation of distinct Solexa tags against reference sequences of B. rapa and B. oleracea, respectively
Clean distinct tag
1 tag->1 gene
1 tag->n gene
1 bp Mismatch(Sense)
1 tag->1 gene
1 tag->n gene
1 tag->1 gene
1 tag->n gene
1 bp Mismatch(Antisense)
1 tag->1 gene
1 tag->n gene
Mapping to genome
1 tag->1 position
1 tag->n position
1 bp Mismatch
1 tag->1 position
1 tag->n position
Unambiguously matched genes
Unknown distinct tag
To estimate whether or not the sequencing depth was sufficient for the transcriptome coverage, the sequencing saturation was analyzed in two libraries. The genes that were mapped by all clean tags and unambiguous clean tags increased with the total number of tags. However, when the sequencing counts reached 2 million tags or higher, the number of detected genes was saturated (Additional file 1: Figure S1).
Differential gene expression between the Ste and Fer libraries
Gene Ontology functional analysis of DEGs
To reveal significantly enriched GO terms in DEGs comparing to the genome background, GO enrichment analysis of functional significance was performed. The GO term with P ≤ 0.05 was defined as significantly DEGs enriched GO term. This analysis allowed us to determine the major biological functions, biological processes and cellular components with which DEGs were associated. For molecular function, significantly enriched GO terms of DEGs in B. rapa (Additional file 5) and B. oleracea (Additional file 6) included transferase activity (transferring acyl groups), aldehyde-lyase activity and oxidoreductase activity, and zinc ion transmembrane transporter activity was specific to B. rapa. For enriched biological processes, there were 19 GO terms of DEGs from B. oleracea (Additional file 7) and 28 GO terms of DEGs from B. rapa (Additional file 8), which participated in morphogenesis (pollen wall assembly, pollen exine formation, cellular component assembly involved in morphogenesis, anatomical structure formation involved in morphogenesis), development (pollen development, gametophyte development), metabolic process (small molecule metabolism, fatty acid biosynthetic metabolism, alcohol metabolic metabolism), ion transport (cadmium ion transport, di-, tri-valent inorganic cation transport), lipid localization and response to stimulus (response to metal ion, response to inorganic substance).
Pathway analysis for DEGs
Expressed genes specific to sterile and fertile plants
DEGs encoding transcription factors
Transcription factors are essential for the regulation of gene expression. Changes in gene transcription are associated with changes in expression of transcription factors. Our DGE results showed that among the significantly DEGs, there were forty-five transcription factors for B. rapa, including fifteen up-regulated and thirty down-regulated genes (Additional file 2). Among up-regulated B. rapa transcriptional factors, six ones were specific to Fer plants, including calmodulin-binding protein, mutant TFIIF-alpha, TBP1 (TATA binding protein 1), TBP1 (Telomere-binding protein 1, DNA binding factor (AT-HSFB2B, ATHSFA2) (Additional file 2). For B. oleracea, there were fifty-seven transcription factors whose expression were dramatically altered, including thirty-nine down-regulated and eighteen up-regulated genes among which six ones were only expressed in Fer. And these transcription factors were ATNAC2, ATHSFA2, transcription factor IID-1 and three unknown proteins (Additional file 2). Further study would be needed to reveal the molecular basis of gene expressional regulation.
Genes involved in carbon metabolism
Many genes involved in carbon metabolism were differentially expressed between Fer and Ste. Altered expressions of numerous genes involved in pentose phosphate pathway, starch and sucrose metabolism, fructose and mannose metabolism and carbon fixation in photosynthetic organisms were observed. For example, in B. rapa, twenty genes involved in pentose phosphate pathway showed changed expression. Those genes with increased transcript abundance in Fer encoded glucose-6-phosphate dehydrogenase (Bra028474), three proteins transferring phosphorus-containing groups (Bra025751, Bra025751 and Bra025751), transaldolase-like pro-tein (Bra006197), ribose-phosphate pyrophosphokinase 1/phosphoribosyl diphosphate synthetase 1 (PRSI) (Bra005344), members performing aldehyde-lyase activity (Bra028543, Bra017233, Bra009517, Bra019639 and Bra006960), and the TPMs for these transcripts were up-regulated by 1.0 to 2.1-fold. The eighty-four genes involved in starch and sucrose metabolism were differentially expressed including two genes encoding UDP-glucosyltransferase proteins (Bra009830, Bra015364) and six genes encoding carbohydrate phosphatase proteins (Bra019370, Bra031879, Bra015497, Bra022549, Bra038548, Bra028254). Their expressions were increased by 1.4 to 10.2-fold in Fer. In fructose and mannose metabolism, there were five down-regulated genes, including 6-phosphofructokinase protein family (Bra030394, Bra038519, Bra011089, Bra012940) and phosphoglucomutase-like protein (Bra016184 ), and five up-regulated gene including pfkB-type carbohydrate kinase family protein (Bra014540, Bra007452, Bra029158, Bra028279) and fructose-1,6-bisphosphatase (Bra014841).
DEGs involved in citrate cycle (TCA cycle) and oxidative phosphorylation
Genes involved in oxidoreductase activity
In enriched GO terms, oxidoreductase activity was the most enriched term. In B. rapa, 311 DEGs were involved in oxidoreductase activity, including 26 members acting on the aldehyde or oxo group of donors, NAD or NADP as acceptor. For instance, fatty acyl-CoA reductase (alcohol-forming)/ oxidoreductase (Bra034793, Bra038691), also named MS2 (male sterility 2), acting on the CH-CH group of donors, NAD or NADP as acceptor, was specific to Fer. Others were oxidoreductase, acting on the CH-CH group of donors (Bra006625), cinnamoyl-CoA reductase family (Bra021863, Bra035150), 3-chloroallyl aldehyde dehydrogenase/ aldehyde dehydrogenase (Bra010553, Bra017753, Bra011674, Bra018090, Bra024619 and Bra016330), NAD dependent epimerase/dehydratase family (Bra024073), aldehyde oxidase (Bra002347), NAD or NADH binding/catalytic/ glyceraldehyde-3-phosphate dehydrogenase (Bra026068, Bra026904, Bra040146 and Bra019797), cytosolic factor family protein/phosphoglyceride transfer family protein (Bra026185), cobalamin-independent methionine synthase (Bra023645), and unnamed oxidoreductase members (Bra032939, Bra011470, Bra011869, Bra005012, Bra008743, Bra040146 and Bra023487).
Differential expressed PPR proteins
Proteins which encode pentatricopeptide repeat (PPR) motif are proposed to function as site-specific, RNA-binding adaptor proteins that mediate interactions between RNA substrates and the enzymes that act on them [34–36]. PPR proteins have a central role in plant mitochondrial biogenesis . Most cloned restorer genes encode mitochondria-targeted PPR proteins [38–44]. In this study, when used B. rapa as reference sequences, there were nine PPR proteins showed differential expression, including two up-regulated PPR proteins (Bra011721 and Bra017411) and seven down-regulated PPR proteins (Bra033892, Bra008587, Bra025761, Bra036629, Bra031484, Bra033866 and Bra017881) in Fer. These PPR proteins were good candidates for studying the link between nuclear and mitochondrial gene expression.
Confirmation of tag-mapped genes by qRT-PCR
The corresponding primers of qRT-PCR
Actin of B. napus
Global gene transcription changes in Fer and Ste plants
A global analysis of transcriptome could facilitate the identification of systemic gene expression and regulatory mechanisms. In the present study, a transcriptome profiling of young buds was performed to identify genes that are differentially expressed in Fer and Ste plants. A sequencing depth of 4.5 million tags per library was reached (Table 1), and more than 71% unique tags were matched with the genome and transcriptome of B. rapa and B. oleracea, respectively, suggesting that the database selected was relatively complete. Because fertile plants carried one chromosome from S. arvensis besides 19 B. napus homologous chromosome pairs, exhibiting normal flower morphology and male fertility, more distinct tags were detected in the Fer library than in the Ste library. In addition, for distinct tags, tags with five folds or greater differences in accumulation accounted for 2.31%, however, among which 2.07% tags increased by at least five folds in Fer. This suggested that more genes were involved in gene expression and gene regulation in Fer. Through annotation of distinct tags, there were 1089 specially expressed unknown tags in Fer, which were neither mapped to B. rapa nor mapped to B. oleracea (Additional file 4). We presumed that these tags arose basically from the added alien chromosome of S. arvensis. Further study of these tags would help to dig genes with desirable agronomic traits from wild species. To identify corresponding genes that were associated with the development of stamens (including anthers) and to elucidate the molecular basis of cytoplasmic male sterility, fertility restoration and improved resistance, we analyzed the most differentially regulated tags with a log2 ratio > 1 or < − 1 using a greater statistically significant value (P < 0.001) as well as false discovery rates (FDR < 0.001). And based on the rigorous algorithm, there were 3231 genes of B. rapa (Additional file 2) and 3371 genes of B. oleracea (Additional file 3) which were detected with significant differential expression levels. These DEGs were candidates for further research.
DEGs specific to sterile and fertile plants
Because fertile plants had normal flower morphology and male fertility versus sterile plants, it was proposed that DEGs specific to fertile plants were associated with stamen development and male fertility, whereas DEGs specific to sterile plants may be involved in CMS. Using B. rapa as reference, a total of 760 DEGs were only expressed in fertile plants, and 44 DEGs were specific to sterile plants. The specific DEGs could be categorized into 56 categories through GO functional analysis, which were involved in biochemistry, metabolism, growth, development, apoptosis and so on (Figure 6). For fertile specific DEGs, 11 DEGs were involved in pollen wall assembly (GO:0010208), including four chalcone synthase genes (Bra000559, Bra034658, Bra011566, Bra017681), three beta-1,3-glucanase genes (Bra028343, Bra037057, Bra038969), two male sterility 2 (MS2) members (Bra034793, Bra038691) and two gycosyl hydrolases family proteins (Bra014979, Bra001918). Chalcone synthase was the key enzyme in the flavonoid biosynthesis pathway, which played an important role in many biological functions in plants, including male fertility . Beta-1,3-glucanase genes were involved in male gametophyte development and pollination. Silencing the expression of beta-1,3-glucanase gene by RNA interference in transgenic rice resulted in male sterility due to disruption of degradation of callose in the anther locules at the early microspore stage . For sterile specific DEGs, 40% of the DEGs were assigned to binding, followed by catalytic activity (33.3%), transporter activity (4.4%), antioxidant activity and molecular transducer activity (2.2% for both) with respect to molecule function category.
DEGs involved in carbon metabolism and energy metabolism
Mitochondrial gene/genome rearrangements may alter the expression of common mitochondrial genes coding for proteins involved in respiration/ATP synthesis, affecting ATP production and/or other processes within mitochondria. A lowered ATP production  or a reduced carbohydrate accumulation  was observed in ‘late stage’ CMS flowers. In this study, 11 of 12 DEGs involved in TCA cycle were down-regulated in Ste (Additional file 10). DEGs were also involved in respiration/ATP synthesis and oxidoreductase activity. A total of 20 genes were involved in the oxidative phosphorylation pathway, having effects on NADH dehydrogenase, succinate dehydrogenase/fumarate redutase and ATP synthase (Additional file 10). Oxidoreductase activity was the most enriched GO term (Additional file 5), whereas oxidative stress during microsporogenesis was thought to induce a premature abortion of tapetal cells due to programmed cell death (PCD) in CMS sunflower  and rice . In addition, altered expressions of numerous genes involved in pentose phosphate pathway, starch and sucrose metabolism, fructose and mannose metabolism and carbon fixation in photosynthetic organisms were also observed (Additional file 10). The nature of the mitochondrial genes that influence the expression of nuclear genes were unclear. And thus further investigations on these candidate genes help to illuminate the primary targets and downstream components of CMS-associated mitochondrial signaling pathways in Nsa CMS line.
Differentially expressed PPR proteins
Nuclear genomes have crucial roles in mitochondrial biogenesis and function [31, 51]. An estimated 10% of eukaryotic nuclear genes encode proteins that are targeted to mitochondria following synthesis on cytosolic ribosomes [32, 52]. It has been proposed that the majority of PPR proteins are targeted to plastid or mitochondria , and many genetic and biochemical studies conclude that PPRs directly bind to a specific RNA sequence and promote anterograde regulation such as post-transcriptional splicing, processing, editing or regulating mRNA stability [53–60]. From about 450 PPR proteins in Arabidopsis, more than 60% are predicted to be targeted to mitochondria . PPR motifs have been suggested to possess binding properties to proteins  as well as to RNA . In humans, a PPR related protein with RNA-binding activity is associated with both nuclear and mitochondrial mRNAs and discussed as a potential candidate for a factor that may link nuclear and mitochondrial gene expression . Mutant versions of this protein were shown to be associated with cytochrome c oxidase deficiency in humans [64, 65]. Hence, PPR proteins are good candidates for nuclear-encoded factors controlling distinct steps of transcript maturation in mitochondria. Many of them may turn out to be essential players in CMS systems. In addition, PPR proteins have been identified as targets of different miRNAs [66, 67]. In this study, when the genome and transcriptome of B. rapa were used as reference sequences, there were nine PPR proteins showed differential expression, including two up-regulated PPR proteins and seven down-regulated PPR proteins in Fer. These PPR proteins were good candidates for studying the link between nuclear and mitochondrial gene expression.
In addition to the genes described above, several other transcript profiles were altered between Ste and Fer. Transcriptional factors play important roles in plant growth and development as well as in hormone and stress responses. A single transcription factor can regulate expression of multiple genes. Therefore, alteration in the expression of transcription factor genes normally results in dramatic changes to a plant [68–70]. As a practical consequence, engineering of transcription factor genes provides a valuable means for manipulation of plants . Our DGE results showed that among the significantly DEGs, there were forty-five transcription factors for B. rapa, these transcription factors might therefore be potent tools to engineer enhanced stress tolerance in oil crops [71, 72].
Two male sterility 2 (MS2) members (Bra034793,Bra038691) were found in Fer, using B. rapa as referece (Additional file 2). MS2 gene determines a nuclear male sterile mutant (monogenic recessive) phenotype in Arabidopsis and was isolated and characterized using the En/Spm-I/dSpm transposon-tagging system . The possible function of the MS2 protein was proposed as a fatty acyl reductase in the formation of pollen wall substances . The homologue of MS2 in B. napus, BnapMS2, was isolated by cold plague screening from B. napus anther specific cDNA library . Li et al.  isolated a fragment homologous to BnapMS2 gene from digenic recessive GMS line S45AB of B. napus using RT-PCR technique, and found that there existed an amino acid divergence between fertile S45B and sterile S45A, which may be the putative male sterility site in S45A. Male sterility gene homolog, designated BcMS2, expressed only in stage III flower buds of male fertile Chinese cabbage-pak-choi “ZUBajh97-01B” and there were no detection in all organs of Polima cytoplasmic male sterility (CMS) line ‘Bpol97-05A’ and Ogura CMS line ‘Bogu97-06A’ . In this study, MS2 members (Bra034793, Bra038691) were only expressed in Fer, whose TPM was 1766.48 and 480.76, respectively, and then it was suggested that the two genes were involved in the development of Fer pollens.
Three up-regulated aldehyde dehydrogenase members (Bra024619, Bra018090, Bra016330) were found in Fer, using B. rapa as referece (Additional file 2). Aldehyde dehydrogenases (ALDHs) have been considered as general detoxifying enzymes which eliminate biogenic and xenobiotic aldehydes in an NAD(P)-dependent manner . Maize nuclear fertility restorer gene Rf2 as “biochemical restorer,” encoding an aldehyde dehydrogenase, plays an important role in oxidizing both aliphatic and aromatic aldehydes, reducing the amount of a toxic aldehyde and rescuing the developing pollen through a metabolic effect [79, 80].
This study has demonstrated the usefulness of the digital gene expression (DGE) approach to identify the differentially expressed genes (DEGs) between the young floral buds of sterile and fertile plants. In total, there were 3231 genes of B. rapa and 3371 genes of B. oleracea which were detected with significant differential expression levels. GO and pathway-based analyses indicate that these DEGs were related to many kinds of molecular functions. In addition, there were 1089 specially expressed unknown tags in Fer, which were neither mapped to B. oleracea nor mapped to B. rapa, and these unique tags were presumed to arise basically from the added alien chromosome of S. arvensis. The results provide a strong basis for future research on CMS in Nsa line, fertility restoration and improved agronomic traits in NR1 line. Further work should focus on characterizing these candidate target genes.
The sterile and fertile individual plants were derived from the self-pollinated offspring of the F1 hybrid between restorer line NR1 and Nsa CMS line, which were grown at the experimental station of the Oil Crops Research Institute, Chinese Academy of Agricultural Sciences, Wuhan, China. For Solexa sequencing, young floral buds (less than 2 mm in diameter) were pooled. A total of 0.5 g materials was snap-frozen in liquid nitrogen and stored at −80°C.
Stamens were fixed with FAA (10% formalin / 5% acetic acid/45% ethanol/0.01% Triton X-100) for 45 min, and rehydrated through an ethanol series. Stamens were then embedded in OCT compound and frozen, sectioned at 10-μm thickness with frozen section machine Leica CM1850 (Germany). Observations were made with an Olympus microscope and photographed.
RNA extraction and Solexa/Illumina sequencing
Total RNA from five sterile and fertile individual plants, respectively was extracted using the TRIzol reagent (Invitrogen). After precipitation, RNA was purified with Qiagen’s RNeasy Kit with on-column DNase digestion according to the manufacturer’s instructions. Purified RNA samples were dissolved in diethylpyrocarbonate-treated H2O, and the concentration determined spectroscopycally. The quality of the RNA was assessed on 1.0% denaturing agarose gel in combination with the Bioanalyzer 2100 (Agilent). For Solexa sequencing, at least 6 μg of total RNA (≥400 ng/μL) was sent to BGI-Shenzhen, China, following the protocol offered by Illumina for sequencing. The main reagents and supplies were the Illumina Gene Expression Sample Prep Kit and Illumina Sequencing Chip (flow cell), and the main instruments were the Illumina Cluster Station and the Illumina HiSeq™ 2000 System. In details, initially, poly-(A) mRNA was isolated from total RNA sample by poly-(T) oligo-attached magnetic beads. Then, the first- and second-strand cDNAs were synthesized and digested with restriction enzyme Nla III, which recognized the CATG sites. Following be washed, the Illumina adaptor 1 was ligated to the 5’ ends of cDNA fragments through the sticky CATG site. Subsequently, cDNA fragments were digested by Mme I, as a type of endonuclease with separated recognition sites and digestion sites, which cut at 17 bp downstream of the CATG site, producing tags with adaptor 1. After removing 3'- fragments with magnetic beads precipitation, Illumina adaptor 2 was ligated to the 3'- ends of tags, acquiring tags with different adaptors of both ends to form a tag library. After PCR amplification, purification and denaturation, the single-chain molecules were fixed onto the Illumina Sequencing Chip (Flowcell). Then, four types of nucleotides labeled by four colors were added in, and sequencing was performed via the sequencing by synthesis (SBS) method.
Processing of sequencing tags and gene expression annotation
Sequencing-received raw image data was transformed by base calling into sequence data, (raw data or raw reads), and was stored in FASTQ format. This type of files stored information about read sequences and quality. Each read was described in four lines in FASTQ files. Raw sequences had 3’- adaptor fragments as well as a few low-quality sequences and several types of impurities. Raw sequences were transformed into clean tags after certain data-processing steps. All clean sequencing tags were annotated using the reference sequences which covered all possible CATG+17-nt tag sequences of the genome and transcriptome of B. rapa and those of B. oleracea, respectively using blastn, allowing only a 1-bp mismatch. Clean tags mapped to reference sequences from multiple genes were filtered, and the remaining clean tags were designated as unambiguous clean tags. The number of unambiguous clean tags for each gene was calculated and then normalized to TPM (number of transcripts per million clean tags) [81, 82].
Analysis and screening of DEGs
A rigorous algorithm  was used to identify the differentially expressed genes (DEGs) between the two samples. The P-value corresponds to the differential gene expression test. The FDR (False Discovery Rate) is used to determine the threshold of P-value in multiple tests and analyses by manipulating the FDR value. Assume that R differentially expressed genes have been selected, among which S genes truly show differential expression and V genes are false positives. If we decide that the error ratio “Q = V/R” must stay below a cutoff (e.g. 1%), we should preset the FDR to a number no larger than 0.01. FDR ≤ 0.001 and the absolute value of | log2Ratio |≥ 1 were used as thresholds to judge the significance of differences in transcript abundance . More stringent criteria with smaller FDR and greater fold-change value can be used to identify DEGs.
Real-time quantitative RT-PCR (qRT-PCR) analysis
Real-time quantitative RT-PCR (qRT-PCR) analysis was used to verify the DGE results. The RNA samples used for the qRT-PCR assays were the same as for the DGE experiments. Gene-specific primers were designed according to the reference unigene sequences using the Primer Premier 5.0 (Table S4). qRT-PCR was performed according to the TaKaRa manufacturer specifications (TaKaRa SYBR® PrimeScript™ RT-PCR Kit, Dalian, China). SYBR Green PCR cycling was denatured using a program of 95°C for 10 s, and 40 cycles of 95°C for 5 s and 55°C for 30 s and performed on an iQ™ 5 Multicolor Real-time PCR Detection System (Bio- RAD, USA). B. rapa actin gene was used as a normalizer, and the relative expression levels of genes were presented by 2-ΔΔCT.
GO and pathway enrichment analysis
For GO enrichment analysis, 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 that are annotated to the certain GO terms; and m is the number of DEGs in M. The p value is corrected by Bonferroni, and we chose a corrected-p value ≤ 0.05 as the threshold value. The GO term (P ≤ 0.05) is defined as significantly DEGs enriched GO term. For pathway enrichment analysis, N is the number of genes with a KEGG annotation, n is the number of DEGs in N, M is the number of genes annotated to specific pathways, and m is the number of DEGs in M. The pathways with a Q value of ≤ 0.05 are defined as those with significantly differentially expressed (enriched) genes.
This work was financially supported by the National Natural Science Foundation of China (grant no. 30671312, 30671119, 31000725), the National High Technology Research and Development Program of China (grant no. 2011AA10A104), the Natural Science Foundation of Hubei Province (grant no. 2008CDA083 and 2009CDB191) and the Chenguang Program of Wuhan City (grant no. 201050231022). Thanks to Assistant Prof. Shunbin Ning (Division of Hematology/Oncology, Department of Medicine, Miller School of Medicine, University of Miami, Miami, FL 33136, USA) for his valuable comments on this manuscript.
- Mei DS, Li YC, Hu Q: Investigation of male sterile lines derived from intergeneric somatic hybrids of Brassica napus (+) Orychophragmus violaceus and B. napus (+) Sinapis arvensis. Chin J Oil Crop Sci. 2003, 25: 72-75.Google Scholar
- Hu Q, Li YC, Mei DS, Fang XP, Hansen LN, Andorsen SB: Establishment and Identification of cytoplasmic male sterility in Brassica napus by intergeneric somatic hybridization. Sci Agric Sin. 2004, 37: 333-338.Google Scholar
- Guan CY: Comparative studies of genetic characters in Xinjiang wild rape and Sinapis arvensis L. Acta Agron Sin. 1996, 22: 214-219.Google Scholar
- Hu Q, Qian XZ, Li YC: Committee of Oil Crops (ed) Advance of Oil Crops Research in China. Analysis of composition and correlation of fatty acids in Xinjiang wild rape. 1996, Beijing: Chinese Agricultural Sciences and Technology Press, 93-96.Google Scholar
- Qian XZ, Guan CY: Review on investigation and collection of Xinjiang and Yunnan wild rape. Oil crops in China. Supplement on the investigation and collection of Xinjiang and Yunnan wild rape. Edited by: Qian XZ. 1988, 1-3.Google Scholar
- Hu Q, Andersen SB, Dixelius C, Hansen LN: Production of fertile intergeneric somatic hybrids between Brassica napus and Sinapis arvensis for the enrichment of the rapeseed gene pool. Plant Cell Rep. 2002, 21: 147-152. 10.1007/s00299-002-0491-7.View ArticleGoogle Scholar
- Ogura H: Studies on the new male-sterility in Japanese radish with special reference to the utilization of this sterility towards the practical raising of hybrid seeds. Mem Fac Agric Kagoshima Univ. 1968, 6: 39-78.Google Scholar
- Shiga T, Baba S: Cytoplasmic male sterility in rape plants (Brassica napus L.). Jpn J Breed. 1971, 21: 16-17.Google Scholar
- Fu TD: Production and research of rapeseed in the People’s Republic of China. Eucarpia Crucif Newsl. 1981, 6: 6-7.Google Scholar
- Mathias R: Transfer of cytoplasmic male sterility from brown mustard (Brassica juncea (L.) Coss.) into rapeseed (Brassica napus L.). Z PXanzenzüchtg. 1985, 95: 371-374.Google Scholar
- Wan ZJ, Jing B, Tu JX, Ma CZ, Shen JX, Yi B, Wen J, Huang T, Wang XJ, Fu TD: Genetic characterization of a new cytoplasmic male sterility system (hau) in Brassica juncea and its transfer to B. napus. Theor Appl Genet. 2008, 116: 355-362. 10.1007/s00122-007-0673-3.View ArticlePubMedGoogle Scholar
- Cheng JH, Li YC, Hu Q, Mei DS, Li YD, Xu YS, Wang WM: Molecular identification and distinctness of Nsa male sterile cytoplasm in Brassica napus. Acta Agron Sin. 2008, 34: 1946-1952.View ArticleGoogle Scholar
- Wang WM: Master theses of the Chinese Academy of Agricultural Sciences. Molecular identification of Nsa cytoplasm in Brassica napus L. 2008Google Scholar
- Wei WH, Li YC, Wang LJ, Liu SY, Yan XH, Mei DS, Li YD, Xu YS, Peng PF, Hu Q: Development of a novel Sinapis arvensis disomic addition line in Brassica napus containing the restorer gene for Nsa CMS and improved resistance to Sclerotinia sclerotiorum and pod shattering. Theor Appl Genet. 2010, 120: 1089-1097. 10.1007/s00122-009-1236-6.View ArticlePubMedGoogle Scholar
- Wei WL, Qi XQ, Wang LH, Zhang YX, Hua W, Li DH, Lv HX, Zhang X: Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers. BMC Genomics. 2011, 12: 451-10.1186/1471-2164-12-451.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang Q, Ma XJ, Mo CM, Wilson IW, Song C, Zhao H, Yang YF, Fu W, Qiu DY: An efficient approach to finding Siraitia grosvenorii triterpene biosynthetic genes by RNA-seq and digital gene expression analysis. BMC Genomics. 2011, 12: 343-10.1186/1471-2164-12-343.PubMed CentralView ArticlePubMedGoogle Scholar
- Qi XH, Xu XW, Lin XJ, Zhang WJ, Chen XH: Identification of differentially expressed genes in cucumber (Cucumis sativus L.) root under waterlogging stress by digital gene expression profile. Genomics. 2012, 99: 160-168. 10.1016/j.ygeno.2011.12.008.View ArticlePubMedGoogle Scholar
- Hao QN, Zhou XA, Sha AH, Wang C, Zhou R, Chen SL: Identification of genes associated with nitrogen use efficiency by genome-wide transcriptional analysis of two soybean genotypes. BMC Genomics. 2011, 12: 525-10.1186/1471-2164-12-525.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang QQ, Liu F, Chen XS, Ma XJ, Zeng HQ, Yang ZM: Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/lintless mutant. Genomics. 2010, 6 (96): 369-376.View ArticleGoogle Scholar
- Eveland AL, Satoh-Nagasawa N, Goldshmidt A, Jackson D: Digital gene expression signatures for maize development. Plant Physiol. 2010, 154: 1024-1039. 10.1104/pp.110.159673.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu T, Qin Z, Zhou X, Feng Z, Du Y: Transcriptome profile analysis of floral sex determination in cucumber. J Plant Physiol. 2010, 167 (11): 905-913. 10.1016/j.jplph.2010.02.004.View ArticlePubMedGoogle Scholar
- Dong CH, Li C, Yan XH, Huang SM, Huang JY, Wang LJ, Guo RX, Lu GY, Zhang XK, Fang XP, Wei WH: Gene expression profiling of Sinapis alba leaves under drought stress and rewatering growth conditions with Illumina deep sequencing. Mol Biol Rep. 2012, 39 (5): 5851-5857. 10.1007/s11033-011-1395-9.View ArticlePubMedGoogle Scholar
- Kempken F, Pring DR: Male sterility in higher plants: fundamentals and applications. Prog Bot. 1999, 60: 139-166. 10.1007/978-3-642-59940-8_6.View ArticleGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-10. 10.1093/nar/30.1.207.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagaharu U: Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. Jap J Bot. 1935, 7: 389-452.Google Scholar
- Kapila R, Negi MS, This P, Delseny M, Srivastava PS, Lakshmikumaran M: A new family of dispersed repeats from Brassica nigra: characterization and localization a new family of dispersed repeats from Brassica nigra: characterization and localization. Theor Appl Genet. 1996, 93: 1123-1129. 10.1007/BF00230135.View ArticlePubMedGoogle Scholar
- Wang XW, Wang HZ, Wang J, Sun RF, Wu J, Liu SY, Bai YQ, Mun JH: The genome of the mesopolyploid crop species Brassica rapa. Nature Genet. 2011, 43: 1035-1039. 10.1038/ng.919.View ArticlePubMedGoogle Scholar
- Audic S, Claverie J: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29: 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36: D480-484.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrientos A: Yeast models of human mitochondrial diseases. IUBMB Life. 2003, 55: 83-95.View ArticlePubMedGoogle Scholar
- Reichert AS, Neupert W: Mitochondriomics or what makes us breathe. Trends Genet. 2004, 20: 555-562. 10.1016/j.tig.2004.08.012.View ArticlePubMedGoogle Scholar
- Logan DC: The mitochondrial compartment. J Exp Bot. 2006, 57: 1225-1243. 10.1093/jxb/erj151.View ArticlePubMedGoogle Scholar
- Lurin C, Andrés C, Aubourg S, Bellaoui M, Bitton F, Bruyère C, Caboche M, Debast C, Gualberto J, Hoffmann B, Lecharny A, Le Ret M, Martin-Magniette ML, Mireau H, Peeters N, Renou JP, Szurek B, Taconnat L, Small I: Genome-wide analysis of Arabidopsis pentatricopeptide repeat proteins reveals their essential role in organelle biogenesis. Plant Cell. 2004, 16: 2089-2103. 10.1105/tpc.104.022236.PubMed CentralView ArticlePubMedGoogle Scholar
- Stern DB, Hanson MR, Barkan A: Genetics and genomics of chloroplast biogenesis: maize as a model system. Trends Plant Sci. 2004, 9: 293-301. 10.1016/j.tplants.2004.04.001.View ArticlePubMedGoogle Scholar
- Shikanai T: RNA editing in plant organelles: machinery, physiological function and evolution. Cell Mol Life Sci. 2006, 63: 698-708. 10.1007/s00018-005-5449-9.View ArticlePubMedGoogle Scholar
- Wen L, Ruesch KL, Ortega VM, Kamps TL, Gabay-Laughnan S, Chase CD: A nuclear restorer of fertility mutation disrupts accumulation of mitochondrial ATP synthase subunit A in developing pollen of S male-sterile maize. Genetics. 2003, 165: 771-779.PubMed CentralPubMedGoogle Scholar
- Wang Z, Zou Y, Li X, Zhang Q, Chen L, Wu H, Su D, Chen Y, Guo J, Luo D, Long Y, Zhong Y, Liu YG: 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. 10.1105/tpc.105.038240.PubMed CentralView ArticlePubMedGoogle Scholar
- Akagi H, Nakamura A, Yokozeki-Misono Y, Inagaki A, Takahashi H, Mori K, Fujimura T: Positional cloning of the rice Rf-1 gene, a restorer of BT-type cytoplasmic male sterility that encodes a mitochondriatargeting PPR protein. Theor Appl Genet. 2004, 108: 1449-1457. 10.1007/s00122-004-1591-2.View ArticlePubMedGoogle Scholar
- Kazama T, Toriyama K: A pentatricopeptide repeatcontaining gene that promotes the processing of aberrant atp6 RNA of cytoplasmic male-sterile rice. FEBS Lett. 2003, 544: 99-102. 10.1016/S0014-5793(03)00480-0.View ArticlePubMedGoogle Scholar
- Brown GG, Formanová N, Jin H, Wargachuk R, Dendy C, Patil P, Laforest M, Zhang J, Cheung WY, Landry BS: The radish Rfo restorer gene of Ogura cytoplasmic male sterility encodes a protein with multiple pentatricopeptide repeats. Plant J. 2003, 35: 262-272. 10.1046/j.1365-313X.2003.01799.x.View ArticlePubMedGoogle Scholar
- Desloire S, Gherbi H, Laloui W, Marhadour S, Clouet V, Cattolico L, Falentin C, Giancola S, Renard M, Budar F, Small I, Caboche M, Delourme R, Bendahmane A: Identification of the fertility restoration locus, Rfo, in radish, as a member of the pentatricopeptide-repeat protein family. EMBO Rep. 2003, 4: 588-594. 10.1038/sj.embor.embor848.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1365-313X.2003.01735.x.View ArticlePubMedGoogle Scholar
- Klein RR, Klein PE, Mullet JE, Minx P, Rooney WL, Schertz KF: Fertility restorer locus Rf1 of sorghum (Sorghum bicolor L.) encodes a pentatricopeptide repeat protein not present in the colinear region of rice chromosome 12. Theor Appl Genet. 2005, 111: 994-1012. 10.1007/s00122-005-2011-y.View ArticlePubMedGoogle Scholar
- van der Meer IM, Stam ME, van Tunen AJ, Moi JNM, Stuitje AR: Antisense inhibition of flavonoid biosynthesis in petunia anthers results in male sterility. Plant Cell. 1992, 4: 253-262.PubMed CentralView ArticlePubMedGoogle Scholar
- Wan L, Zha W, Cheng X, Liu C, Lv L, Liu C, Wang Z, Du B, Chen R, Zhu L, He G: A rice β-1,3-glucanase gene Osg1 is required for callose degradation in pollen development. Planta. 2010, 233: 309-323.View ArticlePubMedGoogle Scholar
- Bergman P, Edqvist J, Farbos I, Glimelius K: Malesterile tobacco displays abnormal mitochondrial atp1 transcript accumulation and reduced floral ATP/ADP ratio. Plant Mol Biol. 2000, 42: 531-544. 10.1023/A:1006388814458.View ArticlePubMedGoogle Scholar
- Datta R, Chamusco KC, Chourey PS: Starch biosynthesis during pollen maturation is associated with altered patterns of gene expression in maize. Plant Physiol. 2002, 130: 1645-1656. 10.1104/pp.006908.PubMed CentralView ArticlePubMedGoogle Scholar
- Balk J, Leaver CJ: The PET1-CMS mitochondrial mutation in sunflower is associated with premature programmed cell death and cytochrome c release. Plant Cell. 2001, 13: 1803-1818.PubMed CentralView ArticlePubMedGoogle Scholar
- Li SQ, Wan CX, Kong J, Zhang ZJ, Li YS, Zhu YG: Programmed cell death during microgenesis in a Honglian CMS line of rice is correlated with oxidative stress in mitochondria. Funct Plant Biol. 2004, 31: 369-376. 10.1071/FP03224.View ArticleGoogle Scholar
- Chinnery PF: Searching for nuclear–mitochondrial genes. Trends Genet. 2003, 19: 60-62. 10.1016/S0168-9525(02)00030-6.View ArticlePubMedGoogle Scholar
- Andersson SG, Karlberg O, Canbäck B, Kurland CG: On the origin of mitochondria: a genomics perspective. Philos Trans R Soc London Ser B. 2002, 358: 165-179.View ArticleGoogle Scholar
- Fisk DG, Walker MB, Barkan A: Molecular cloning of the maize gene crp1 reveals similarity between regulators of mitochondrial and chloroplast gene expression. EMBO J. 2000, 18: 2621-2630.View ArticleGoogle Scholar
- Nakamura T, Meirhoff K, Westhoff P, Schuster G: RNA binding properties of HCF152, an Arabidopsis PPR protein involved in the processing of chloroplast RNA. Eur J Biochem. 2003, 270: 4070-4081. 10.1046/j.1432-1033.2003.03796.x.View ArticlePubMedGoogle Scholar
- Kotera E, Tasaka M, Shikanai T: A pentatricopeptide repeat protein is essential for RNA editing in chloroplasts. Nature. 2004, 433: 326-330.View ArticleGoogle Scholar
- Schmitz-Linneweber C, Williams-Carrier R, Barkan A: RNA immunoprecipitation and microarray analysis show a chloroplast pentatricopeptide repeat protein to be associated with the 5' region of mRNAs whose translation it activates. Plant Cell. 2005, 17: 2791-2804. 10.1105/tpc.105.034454.PubMed CentralView ArticlePubMedGoogle Scholar
- Okuda K, Nakamura T, Sugita M, Shimizu T, Shikanai T: A pentatricopeptide repeat protein is a site recognition factor in chloroplast RNA editing. J Bio Chem. 2006, 281: 37661-37667. 10.1074/jbc.M608184200.View ArticleGoogle Scholar
- Okuda K, Myouga F, Motohashi R, Shinozaki K, Shikanai T: Conserved domain structure of pentatricopeptide repeat proteins involved in chloroplast RNA editing. Proc Natl Acad Sci USA. 2007, 104: 8178-8183. 10.1073/pnas.0700865104.PubMed CentralView ArticlePubMedGoogle Scholar
- Hattori M, Miyake H, Sugita M: A pentatricopeptide repeat protein is required for RNA processing of clpP pre-mRNA in moss chloroplasts. J Biol Chem. 2007, 282: 10773-10782. 10.1074/jbc.M608034200.View ArticlePubMedGoogle Scholar
- de Longevialle AF, Meyer EH, Andrés C, Taylor NL, Lurin C, Millar AH, Small ID: The pentatricopeptide repeat gene OTP43 is required for trans-splicing of the mitochondrial nad1 intron 1 in Arabidopsis thaliana. Plant Cell. 2007, 19: 3256-3265. 10.1105/tpc.107.054841.PubMed CentralView ArticlePubMedGoogle Scholar
- Williams PM, Barkan A: A chloroplast-localized PPR protein required for plastid ribosome accumulation. Plant J. 2003, 36: 675-686. 10.1046/j.1365-313X.2003.01915.x.View ArticlePubMedGoogle Scholar
- Small ID, Peeters N: The PPR motif: a TPR-related motif prevalent in plant organellar proteins. Trends Biochem Sci. 2000, 25: 46-47.View ArticlePubMedGoogle Scholar
- Mili S, Pinol-Roma S: LRP130, a pentatricopeptide motif protein with a noncanonical RNA-binding domain, is bound in vivo to mitochondrial and nuclear RNAs. Mol Cell Biol. 2003, 23: 4972-4982. 10.1128/MCB.23.14.4972-4982.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Mootha VK, Lepage P, Miller K, Bunkenborg J, Reich M, Hjerrild M, Delmonte T, Villeneuve A, Sladek R, Xu F, Mitchell GA, Morin C, Mann M, Hudson TJ, Robinson B, Rioux JD, Lander ES: Identification of a gene causing human cytochrome c oxidase deficiency by integrative genomics. Proc Natl Acad Sci USA. 2003, 100: 605-610. 10.1073/pnas.242716699.PubMed CentralView ArticlePubMedGoogle Scholar
- Shadel GS: Coupling the mitochondrial transcription machinery to human disease. Trends Genet. 2004, 20: 513-519. 10.1016/j.tig.2004.08.005.View ArticlePubMedGoogle Scholar
- Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant microRNA targets. Cell. 2002, 110: 513-520. 10.1016/S0092-8674(02)00863-2.View ArticlePubMedGoogle Scholar
- Sunkar R, Zhu JK: Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004, 16: 2001-2019. 10.1105/tpc.104.022830.PubMed CentralView ArticlePubMedGoogle Scholar
- Kater MM, Colombo L, Franken J, Busscher M, Masiero S, Van Lookeren Campagne MM, Angenent GC: Multiple AGAMOUS homologs from cucumber and petunia differ in their ability to induce reproductive organ fate. Plant Cell. 1998, 10: 171-182.PubMed CentralView ArticlePubMedGoogle Scholar
- Grotewold E, Chamberlin M, Snook M, Siame B, Butler L, Swenson J, Maddock S, Clair GS, Bowen B: Engineering secondary metabolism in maize cells by ectopic expression of transcription factors. Plant Cell. 1998, 10: 721-740.PubMed CentralPubMedGoogle Scholar
- Chen JJ, Janssen BJ, Williams A, Sinha N: A gene fusion at a homeobox locus: alterations in leaf shape and implications for morphological evolution. Plant Cell. 1997, 9: 1289-1304.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen W, Provart NJ, Glazebrook J, Zhu T: Expression profile matrix of Arabidopsis transcription factor genes suggests their putative functions in response to environmental stresses. Plant Cell. 2002, 14: 559-574. 10.1105/tpc.010410.PubMed CentralView ArticlePubMedGoogle Scholar
- Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF: Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998, 280: 104-106. 10.1126/science.280.5360.104.View ArticlePubMedGoogle Scholar
- Aarts MGM, Dirkse WG, Stiekema WJ, Pereira A: Transposon tagging of a male sterility gene in Arabidopsis. Nature. 1993, 363 (24): 715-717.View ArticlePubMedGoogle Scholar
- Aarts MGM, Hodge R, Kalantidis K, Florack D, Wilson ZA, Mulligan BJ, Stiekema WJ, Scott R, Pereira A: The Arabidopsis MALE STERILITY 2 protein shares similarity with reductase in elongation/condensation complexes. Plant J. 1997, 12 (3): 615-623. 10.1046/j.1365-313X.1997.d01-8.x.View ArticlePubMedGoogle Scholar
- Hodge R, Paul W, Draper J, Scott R: Cold-plaque screening: a simple technique for the isolation of low abundance, diVerentially expressed transcripts from conventional cDNA libraries. Plant J. 1992, 2: 257-260. 10.1111/j.1365-313X.1992.00257.x.View ArticleGoogle Scholar
- Li DM, Hou L, Luo XY, Pei Y, Yang GW: Cloning and sequence analysis of fragment homologous to MS2Bnap gene in rapeseed recessive genic male sterile line S45AB (Brassica napus L). Acta Agron Sin. 2002, 28 (1): 1-5.Google Scholar
- Jiang M, Cao JS: Isolation and characterization of a male sterility gene homolog BcMS2 from Chinese cabbagge-pak-choi that expressing in an anther-specific manner. Mol Biol Rep. 2008, 35: 299-305. 10.1007/s11033-007-9086-2.View ArticlePubMedGoogle Scholar
- Yoshida A, Rzhetsky A, Hsu LC, Chang C: Human aldehyde dehydrogenase gene family. Eur J Biochem. 1998, 251: 549-557. 10.1046/j.1432-1327.1998.2510549.x.View ArticlePubMedGoogle Scholar
- Liu F, Cui X, Horner HT, Weiner H, Schnable PS: Mitochondrial aldehyde dehydrogenase activity is required for male fertility in maize. Plant Cell. 2001, 13: 1063-1078.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu F, Schnable PS: Functional specialization of maize mitochondrial aldehyde dehydrogenases. Plant Physiol. 2002, 130: 1657-1674. 10.1104/pp.012336.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoen PA T, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GJ, den Dunnen JT: Deep sequencing based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36: e141-10.1093/nar/gkn705.PubMed CentralView ArticlePubMedGoogle Scholar
- Morrissy AS, Morin RD, Delaney A, Zeng T, McDonald H, Jones S, Zhao Y, Hirst M, Marra MA: Next-generation tag sequencing for cancer gene expression profiling. Genome Res. 2009, 19: 1825-1835. 10.1101/gr.094482.109.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.