- Research article
- Open Access
Transcriptomic analysis reveals the mechanism of thermosensitive genic male sterility (TGMS) of Brassica napus under the high temperature inducement
BMC Genomics volume 20, Article number: 644 (2019)
The thermo-sensitive genic male sterility (TGMS) of Brassica napus facilitates reproductive researches and hybrid seed production. Considering the complexity and little information about the molecular mechanism involved in B. napus TGMS, comparative transcriptomic analyses were peroformed for the sterile (160S-MS) and fertile (160S-MF) flowers to identify potential crucial genes and pathways associated with TGMS.
In total, RNA-seq analysis showed that 2202 genes (561 up-regulated and 1641 down-regulated) were significantly differentially expressed in the fertile flowers of 160S-MF at 25 °C when compared the sterile flower of 160S-MS at 15 °C. Detailed analysis revealed that expression changes in genes encoding heat shock proteins, antioxidant, skeleton protein, GTPase and calmodulin might be involved in TGMS of B. napus. Moreover, gene expression of some key members in plant hormone signaling pathways, such as auxin, gibberellins, jasmonic acid, abscisic acid, brassinosteroid signalings, were significantly surppressed in the flowers of 160S, suggesting that these genes might be involved in the regulation in B. napus TGMS. Here, we also found that transcription factor MADS, NFY, HSF, MYB/C and WRKY might play a crucial role in male fertility under the high temperature condition.
High temperature can significant affect gene expression in the flowers. The findings in the current study improve our understanding of B. napus TGMS at the molecular level and also provide an effective foundation for male fertility researches in other important economic crops.
Heterosis has been applied to increase the crop yield, alleviate the contradiction between quality and production and stress resistance. Male sterile (MS) is one of the most important method utilized in the heterosis. As a common biological phenomenon, plant MS system is also valuable for studies on anther development, organogenesis, cytoplasmic inheritance, and nucleo- cytoplasmic interactions [1,2,3]. Presently, MS types used in the hybrid breeding mainly include cytoplasmic male sterility (CMS), genic male sterility (GMS), self incompatibility (SI) and chemical gametocide (CG) [4,5,6,7]. CMS is characterized by the production of non-functional pollen and has been used to generate hybrids with significant heterosis in some crops. However, this sterility type is not stable in sterile line (sometimes carrying a trace of fertile pollen). Moreover, parental screening is affected by the restore line and/or the maintenance line, and also negatively affected by the cytoplasm. Although, three-line system has been applied for the hybrid seed production, a pure CMS line is difficult to be obtained. The CG method is influenced by the weather conditions and the contents of chemical gametocide agents and application time are hard to control. Therefore, the development of new male sterile materials and understanding the molecular mechanism are essential for rapeseed cultivation and other economic crops. Since the first report of photo/thermo-sensitive genic male sterility (P/TGMS) of rice 58S was reported , many P/TGMS lines in some plant species have been developed , which is characterized by the fertility that is not only genetically controlled but also affected by the temperature, photoperiod and other environmental factors. Thermo-sensitive genic male sterility (TGMS) is a special sterile system, in which anther development can be depressed under the sterile conditions and non-functional anthers are produced. While under the fertilizing conditions, the depression is totally or partially removed and functional anthers can be produced.
Great progresses have been made in understanding the molecular mechanisms of rice P/TGMS. A mutation in a short-form RNase Z produce a TGMS rice and currently has been used for hybrid rice production . The TGMS trait of Hengnong S-1 is controlled by tms9–1, encoding a PHD-finger protein [11, 12]. Over-expression of UPD-glucose pyrophosphorylase 1(Ugp-1) in rice can produce TGMS trait under the normal temperatures and this trait can be reverted at a low temperature. In addition, many genes encoding transcription factors, members of some signal transduction pathway or proteins involved in metabolisms have been identified from low-temperature treated rice TGMS Y58S and Pei ai 64S . Interestingly, methylation level is significantly higher in the sterile Pei ai 64S than in fertile Pei ai 64S [14, 15]. These results enriched our understanding of the molecular mechanisms of TGMS.
Several ecotypical TGMS rapeseed lines have been identified and applied to produce hybrids, but development and application of novel TGMS lines will enrich the utilization of the heterosis to increase yield and oil seeds quality. Moreover, it is also helpful for understanding the molecular mechanism of the male sterility. So far, numerous differentially expressed genes have been identified in flowers or anthers of Brassica oleracea [16, 17], Brassica rapa , Brassica juncea , and Brassica napus [20,21,22], which outlines an intricate transcriptional regulation that occurs in TGMS system when the temperature is changed. Although numerous efforts have been made in molecular dissections of the fertility of sterile materials, conversion between the fertility and sterility, and anther abortion is not well-known. Due to the complexity of TGMS, previous studies on the regulatory mechanism are far from enough in the depth and breadth. To further explore the mechanisms and promote their wider applications, it is necessary to study it more deeply and systematically.
B. napus is a world widely economic crop and an important source of edible oil. Our previous study showed that the temperature below 20 °C can induce the male fertility (160S-MF) of 160S, whereas the temperature above 25 °C can induce the sterility (160S-MS). The fertility rate, pod-setting ratio and seed number per pod correspondingly decreased with the increased temperature , suggesting that the temperature is a key factor to affect the morphology of flower and the fertility of pollens. However, the molecular regulatory mechanisms are unknown. Therefore, as the first step towards understanding the molecular mechanism of the TGMS in rapeseed, we aim to mine the sterility-related gene by RNA-seq approach, which has been widely used to identify deferentially expressed genes from two or more compared samples. In this study, the flowers of 160S-MF and 160S-MS were used to characterize gene expression patterns. Eight TGMS-related candidate genes were further validated by quantitative realtime-PCR (qPCR). In this way, this study explored the relationship between the male sterility/fertility and gene expression patterns, to obtain more insights in TGMS and offered valuable genetic information for breed improvement and the application of B. napus.
Results and discussion
Two transcriptome libraries were constructed from the flowers of 160S-MF and 160S-MS. High throughput sequencing generated 57,818,154 raw reads for 160S-MF and 68,606,002 raw reads for 160S-MS. After removing the adaptor, low-quality and ambiguous nucleotides, 45,263,606 and 60,008,172 high quality reads were obtained for 160S-MF and 160S-MS, respectively. The average length of the clean reads is 98.66 bp in 160S-MF and 98.7 bp in 160S-MS (Additional file 1: Table S1). All the raw transcriptome data have been deposited at the GenBank (PRJNA. 513,109). Reads assembly produced 101,040 unigenes with 46% GC content.
All assembled unigenes were analyzed by BLAST search against the NR, Swiss-Prot, PFAM, TREMBL and KOG databases (E < 1e-5). The results showed that 65.35% unigenes significantly matched to the genes in TREMBL database, followed by 60.0% in NR database, 44.47% in PFAM, 35.23% in SWISS-PROT database and 10.65% in KOG database (Additional file 2: Figure S1A). Unmatched unigenes could be attributable to the short sequence reads, or might be unique to B. napus, or the relatively short sequences lacked conserved functional domains. The hitted homologue distribution analysis revealed that 36% unigenes was matched to the homologues in Arabidopsis thaliana, 28% matched to Arabidopsis lyrata subsp. lyrata, and 23% matched to Capsella rubella (Additional file 2: Figure S1B).
GO assignments were used to classify the gene functions using the Blast2GO. In all, 31,482 unigenes were mapped and categorized into “biological process”, “cellular component” and “biological process” (Fig. 1). In biological process category, “metabolic process” (GO:008152) presented the largest cluster, accounting for 22.0%, followed by “cellular process” (GO:009987; 18.1%) were the most abundant groups. In cell component category, “cell” (GO:005623; 10.3%) and cell apart (GO:0044464; 10.3%) were predominant. Under the molecular function category, genes related to “binding” (GO:005488; 19.3%) and “catalytic activities” (GO:003824; 18.1%) were the most highly represented. Detailed GO data were shown in Additional file 3: Table S2.
KOG and KEGG pathway enrichment
KOG classification indicated that 31,482 unigenes were clustered into 25 functional categories (Fig. 2). Unigenes in “post-translational modification, protein turnover, chaperones” category were the most abundant (1568 ungienes, 1.55%), followed by “signal transduction mechanisms” (1345 unigenes, 1.33%), “general function prediction only” (1207 unigenes, 1.19%), “intracellular trafficking, secretion, and vesicular transport” (811 unigenes, 0.8%) and “carbohydrate transport and metabolism” (809 unigenes, 0.80%), “nuclear structure” (73 unigenes, 0.07%), “defense mechanisms” (65 unigenes, 0.06%), “extracellular structures” (19 unigenes, 0.02%) and “cell motility” (62 unigenes, 0.009%) was the smallest group. Detailed KOG data were shown in Additional file 4: Table S3. To better understand the gene functional differences, the up- and down-regulated differentially expressed genes (DEGs) were grouped by KEGG (Fig. 3). The results revealed that 27,817 (88.4%) transcripts were assigned to 295 pathways (Additional file 5: Table S4). The top 5 pathways were associated with ribosome (ko03010, 1281 genes), plant hormone signal transduction (ko04075, 1116 genes), protein processing in endoplasmic reticulum (ko04141, 800 genes), spliceosome (ko03040, 731 genes), and starch and sucrose metabolism (ko00500, 720 genes).
Differentially expressed genes between 160S-MF and 160S-MS
The expression levels were measured using reads per kilobase of exon model per million mapped sequence reads (RPKM). RPKM values for each gene in 160S-MS and 160S-MF were compared to identify DEGs and 2202 DEGs (561 up-regulated and 1641 down-regulated) were produced.
To further determine how these DEGs were involved in the biological processes, KEGG analysis was performed and 204 pathways were enriched (Additional file 6: Table S5). For up-regulated genes in the flowers of 160S-MS, the top 5 pathways were associated with plant hormone signal transduction (ko04075, 13 genes), glyoxylate and dicarboxylate metabolism (ko00630, 12 genes), protein processing in endoplasmic reticulum (ko04141, 11 genes), carbon fixation in photosynthetic organisms (ko00710, 9 genes), and starch and spliceosome (ko03040, 8 genes). For the down-regulated genes in 160S-MS flowers, the top 5 pathway mainly associated with starch and sucrose metabolism (ko00500, 31 genes), regulation of actin cytoskeleton (ko04810, 26 genes), pentose and glucuronate interconversions (ko00040, 25 genes), Fc gamma R-mediated phagocytosis (ko04666, 21 genes) and shigellosis (ko05131, 16 genes).
In this study, it was worthy to be mentioned that unigenes in the pathways responsible for thermo-sensitive genic male sterility included oocyte meiosis, plant hormone signal transduction, RNA degradation, ubiquitin mediated proteolysis, protein processing in endoplasmic reticulum, spliceosome, and mRNA surveillance pathway may provide valuable resource for the identification of unique genes involved in TGMS.
Validation of differentially expressed genes
To confirm the expression profiling obtained by RNA-seq analysis, the expression levels of 8 genes (encoding auxin-induced protein, auxin-induced protein6B, HSP70–1, HSP70–2, proflin-1, proflin-2, CML17 and probable pectate lyase 3) putatively involved in MS and showing differential expression patterns between 160S-MS and 160S-MF were further validated by qPCR. Our results revealed that the expression level change of all the selected genes based on qPCR analysis agreed with those detected by RNA-seq analysis (Fig. 4). However, the fold changes in expression level of some genes were different depending on the detection method. The transcript levels of auxin-induce protein and proflin-2 were only detected in 160S-MS. Profilin-1, CML17 and probable pectate lyase 3 were significantly up-regulated in 160S-MS, while auxin-induced protein 6B, HSP70–1 and HSP70–2 were remarkably up-regulated in 160S-MF.
Possible molecular mechanisms of male sterility of B. napus under the higher temperature
Developing pollen and the surrounding tapetal cells show a high sensitivity to heat stress that often results in premature degeneration of tapetal cells and aberrant developmental or programmed cell death. While the cause of this sensitivity remains largely unknown, we suggest a schematic model to illustrate the possible roles of some key players in the developing pollen that might be related to TGMS under heat stress (Fig. 5).
Abiotic stress and sterility
Heat shock proteins
Natural floral organ abortion or degeneration results in unisexual or fully sterile flower, while abiotic factors lead to sterility after initiation of floral reproductive organs. Even a mild or short term of abiotic stress can cause a remarkable decrease in fertility. From floral emergence to 3 days post anthesis, heat stress (36 (day)/31 °C (night)) results in male sterility of wheat due to abnormal pollen grains . Barley grown under the high temperature (30–35 °C day/20–25 °C night) displayed tapetum degradation . Rices that receives a short-term (33.7 °C, 1 h) or a long-term high temperature explosure (38 °C, 41 °C, 5 days) displayed significantly lower fertility, but a better fertility was observed when heat stress was applied before or after anthesis . Although the mechanism of the heat-induced sterility is not clear, it might be related to heat shock proteins (HSP). In Arabidopsis, disruption of HSP40 caused pollen tubes to burst and decrease pollen tubes length .
Similarly, mutaion in a small heat shock protein gene BOB1 results in irregular flowers and sterilize siliques of Arabidopsis . In maize, pollen sterility may be due to the missing protection of some HSPs [29, 30], supporting the fact that some HSP genes did not express at dehiscence . In rice, heat stress induces many HSP genes, including HSP70, HSP26 and HSP17 . Interestingly, six HSP genes (HSP17.6, 4 HSP70s and HSP90) were significantly down-regulated in 160S-MS flowers (Additional file 7: Table S6), implying that low transcription level of HSP genes may contribute male sterility of B. napus. Of course, more work needs to be done to understand the detail role of HSP in male sterility.
Another major mechanism underlining high temperature-induced damages in developing pollen and/or tapetal cells is through the accumulation of reactive oxygen species (ROS), which typically causes membrane degradation, reduced translation and transcription, and eventually programmed cell death (PCD) [33,34,35]. Plants have evolved an effective non-enzymatic system (e.g. glutathione, ascorbic acid) and enzymatic system (e.g. SOD, APX, CAT, GR and GST) to scavenge the accumulated ROS. In many cytoplasmic male sterile crop varieties (e.g. rice, peper, cotton), pollen PCD is caused by an excessive mitochondrial accumulation of ROS together with a decreased scavenging capacity in the developing spores. Under the high temperature condition (25 °C), decreased expressions of eight antioxidant-related genes (encoding SOD, CAT, GST and POD) were observed in 160S-MS (Additional file 7: Table S6), which was consistent with our previous results that antioxidative enzyme activities were higher 160S-MF . It is possible that high temperature leads to higher respiration rate and excessive ROS is accumulated in 160S-MS flowers, thus low activties of antioxidative enzymes scavenge them defeciencely and result in the pollen abortion.
Role of hormones in male sterility
Plant hormones have strong influences on flower and fertility. Some hormones are essential for both male and female organ development, while others are male or female specific . To investigate the relationship between the hormone and the male sterility in B. napus, the enriched hormone signaling pathways were analyzed (Additional file 7: Table S6).
Due to the importance of auxin in floral organs formation, disruption of genes associated with auxin biosynthesis, transport and signaling leads to flowers with various abnormalities. In barley and Arabidopsis anthers, heat stress (5 day at 30/25 °C day/night in barley and 31–33 °C in Arabidopsis) significantly reduced endogenous auxin level in developing pollen mother cells (PMCs) and tapetal cells through suppressing the transcript level of Yuc2 andYuc6 (encoding auxin biosynthesis enzymes) which leads to a premature abortion of microspore development [37, 38]. Auxin response factors (ARFs) repress or activate expression of auxin response genes. In Abrabidopsis, the arf6/arf8 double mutant displays shorten petal and stamen filaments, failure to release pollen, and abnormal ovules [39, 40]. Similarly, double mutation of arf3(attin)/let causes a decreased number of stamens , while mutation in arf5/mp leads to either small or absent lateral flowers . The differential expression of arf7 in 16S-MS and 160S-MF might be related with male sterility.
Gibberellins (GA) is necessary for filament elongation and pollen development in Arabidopsis and Oryza sativa. GA signaling functions through the degradation of DELLA proteins, a class of transcriptional repressors that inhibit GA-dependent variations in the expression of down-stream target genes [43, 44]. DELLA loss-function mutant displays sterility phenotype due to abnormal microsporogenesis and delayed growth of all floral orangs, e.g. stamens with shortened filament that cannot pollinate pistils . As a key member of GA signaling pathway, four DELLA gene was significantly up-regulated in 160S-MS flowers (Additional file 7: Table S6), highlighting that GA is required for male sterility in B. napus. Another analysis showed that DELLA proteins could suppress the synthesis of JA that is required for the activation of transcription factor MYB21/24/57, which is essential for late stamen filament growth in Arabdopsis .
Jasmonic acid (JA) has been known as a growth regulator acting as a signal for plant growth, development and fertility. JA signaling is activated jasmonate-ZIM domain proteins (JAZ) that binds to and suppress the activities of transcription factors, such as MYC family, that regulating the expression of JA-responsive genes [47, 48]. In Arabidopsis and maize tassel, JA is also identified as a factor affecting anther dehiscence, stamen and pollen maturation [49, 50]. The opr mutant is deficient in JA synthesis and produces staments with abnormal filament elongation and dehiscence in Arabidopsis . Silimarly, dobule mutation of opr7/opr8 in maize forms feminized tassel devoid of stamen formation, and extreme elongation of ear shanks . Jasmonate acid carboxyl methyltransferase (JMT) is a key enzyme responsible of convertion of JA to MeJA. The down-expression of JMT and JAZ genes in 160S-MS flowers (Additional file 7: Table S6) would result in a decreased JA level. This deficiency might contribute male sterility in B. napus. Further study looking into the roles of these enzymes on sterility could be valuable.
As a key hormone involved in abiotic stress resposnes, abscisic acid (ABA) plays an important role in male fertility during reproductive stress . In rice, ABA regulates the expression of tapetum cell wall bound invertase and monosaccharide transport genes, leading to a disturbed carbohydrate accumulation in the anther and thus, in pollen sterility . Regulatory component of ABA receptor PYR/PYL senses ABA and initiates the ABA signaling. The expression level of two PYR/PYL genes were significantly decreased in 160S-MS (Additional file 7: Table S6), suggesting that the lack of ABA signal might be involved in B. napus sterility. Because both ABA biosynthesis and catabolism afftect ABA signaling in plant tissues, therefore further studies focusing on the mechanism that control ABA homeostasis during the reproductive stages would provide useful information for male sterility study.
Brassinosteroids (BRs) are perceived by a plasma membrane-localized receptor, brassinosteroid insensitive 1(BRI1), and essential for male fertility in plants. In Arabidopsis, BR-deficient cpd mutant displays male sterility due to the failure of pollen tube elongation. Similarly, mutation of dwf4 results in reduced filament elongation, leading to a failure of pollen delivery to the stigma [55, 56]. However, the underlying mechanism of BRs in regulating anther/pollen development is completely unknown. As the receptor of BR, two BRI1 genes was remarkably down-regulated in 160S-MS flowers (Additional file 7: Table S6), implying that BR controls male sterility at least in part via directly regulating key genes for anther and/or pollen development in B. napus.
Molecular genetic regulation of male sterility
Many transcription factors (TFs) are well known to control plant growth and development and several TFs have been reported to be specially related to anthers and pollen development . In this study, genes encoding TFs, such as MADS-box transcription factor (MADS), nuclear transcription factor Y(NFY), heat shock transcription factor (HSF), MYB/C and WRKY were identified and deferentially expressed in 160S-MS flowers and 160S-MF flowers (Additional file 7: Table S6). The MADS-box genes encode a family of TFs that control diverse developmental processes such as meristem identity, flowering time determination, floral organ identity and silique morphogenesis . In rice and Arabidopsis, MIKC type MADS-box genes was expressed in the gametophyte, while MADS62, MADS63 and MADS68 were all specifically expressed late in pollen development and their reduced expression result in defective anther development and degenerated pollen grains. Knock-down of MADS2 resulted in anther abortion and pollen development defect in maize . In Oryza sativa, MADS3 is required for stamen identity determination during early flower development and also played a critical role in late anther development by modulating ROS homeostasis . Nine MADS genes showed significantly lower expression in sterile flower, which suggested that them may be involved in floral organ identity or other functions during anther development. The nuclear factor Y(NFY) transcription factors are considered as important regulators of many plant developmental processes, including tolerance/resistance to abiotic stresses and male sterility . Two α-subunit genes (NFYA) were down-regulated in sterile flowers, indicating that they might be involved in thermo-response during the anther development.
The MYB superfamily is characterized by a conserved DNA-binding domain, the MYB domain, and is one of the largest TF families in plants. Some MYB transcription factors have been found to play important role in flower development. In Arabidopsis thaliana, disruption of AtMYB26 results in male sterility in due to non-dehiscent anthers . AtMYB103 is especially expressed in anthers and its decreased expression results in early degeneration of the tapetum and abnormal pollen grains . AtMYB24 is induced by jasmonates and is essential for filament elongation . AtMYB80 regulates programmed cell death of the tapetum in the anther, and thus required for pollen development . However, it is not known whether any other MYB transcription factor also plays an essential role in anther, stamen and pollen development, pollen tube growth, especially in agriculture important crops. In this study, the transcript levels of two MYB genes were increased in sterile flowers, but 17 MYB genes were down-regulated. Further analysis will be performed to determine whether some of them participates the regulation of the male sterility of B. napus.
Other miscellaneous molecules putatively involved in male sterility
The actin cytoskeleton is known to play key roles in the morphogenesis and function of highly specialized cell types. Aberrant regulation of actin- and tubulin-related genes can disrupt the organization of actin and microtubules in meiotic process, which leads to defective cytokinesis in developing pollens and male sterility . It has been suggested that low actin levels in anthers are associated with the male sterility in some plants . In the present study, 17 actin-related genes were up-regulated and one β-tubulin-related gene was down-regulated in the flowers of 160S-MS (Additional file 7: Table S6). Actin depolymerizing factor (ADF), an actin binding protein, plays a role in regulating F-actin filament assembly. In Arabidopsis, the overexpression of AtADF1 resulted in the disappearance of thick actin cables in different cell types and caused irregular cellular and tissue morphogenesis, and development arrest. In contrast, the down-regulation of AtADF could promote the formation of actin cables and resulted in a delay in the flowering and organ growth . High expression of TaADF (4.28- fold change) was also found in SQ-1 treated wheat sterile line compared with the fertilize line, which caused the microfilament desambly and the accumulation of callose . In the present study, 7 actin depolymerization factor/cofilin-like domains (ADF domains) genes were up-regulated in 160S-MF flowers at 25 °C, suggesting that the defective ckeleton dynamics in the anther might be associated with the male sterile.
Rho of plant (ROP)
ROP proteins are members of GTPases and function in many fundamental cellular activities, such as cell polarity establishment in pollen tube, cell morphogenesis, regulation of actin cytoskeleton, and hormonal response, and also been implicated in abiotic and biotic stress response [67,68,69] . In Arabidopsis, a mutation in SPIKE1, coding for a ROP protein, leads to the stunted growth and severe defects in polarized cell growth . Additionally, the epidermal pavement cell morphogenesis is regulated by the counter signaling of two Rop-mediated antagonistic pathways, the coordination between plant hormone auxin and specific Rop GTPases together organizes and restructures the cytoskeletal elements for cell morphogenesis and patterning . We found that the transcript levels of 3 ROP genes were significantly increased in 160S-MS flowers (Additional file 7: Table S6), suggesting that they might be associated with the thermo-sensitive male sterility. Till now, only a few ROPs and their interacting proteins have been characterized and more detailed studies will be needed to generate further insights into the related signaling pathways in TGMS.
In plant, calcium-dependent protein kinase (CPK)-mediated signal has been implicated in many biological processes, such as environmental stress response, pollen tube growth and fertilization . In maize, several CPKs were enriched in pollen and might serve as an essential switch in regulating maize pollen tube growth . For example, ZmCPK32 is specifically accumulated in pollen during shedding and putatively associated with pollen tube development . Despite numerous reports on the role of calcium and calmodulin (CaM) in pollen germination and tube growth, the proteins that mediate calcium/CaM action have not been identified in B.napus. In this study, the expression levels of 10 CaM genes (2 CML3, 2 CML17, 2 CML25, CML7, CML42, Polcalcin Bra r1 and 2) were up-regulated, while 2 genes (CML21, CML35) were down-regulated in 160S-MS flowers (Additional file 7: Table S6), which supported the mechanism that calcium signaling pathway is responsible for male fertility.
In this study, 2202 genes were significantly differentially expressed in the flowers at 25 °C compared with these at 15 °C. Functional enrichment analysis revealed that genes encoding heat shock proteins, antioxidant, skeleton protein, and calmodulin might be involved in TGMS of B. napus. Interestingly, the transcription levels of some key members in signaling pathways, intermediated by auxin, gibberellins, jasmonic acid, abscisic acid, brassinosteroid signalings, were significantly surppressed in sterile flowers. Moreover, well-known transcription factor MADS, HSF, MYB/C and WRKY might be involved in male fertility under the high temperature. Our results are helpful for understanding the molecular mechanisms of B. napus TGMS and also provide an effective foundation for male fertility researches in other important economic crops.
Brassica napus 160S were grown at 15 ± 1 °C under long day length before the bud forming, and then moved into the climate chambers at 25 ± 1 °C with a photoperiod of 14 h light and 10 h dark, and 65% relative humidity (RH) until the sterile flower forming. The naturally growing young sterile flowers were collected, frozen in liquid nitrogen, and stored at − 80 °C for later use.
Illumina sequencing and de novo assembly
Total RNA from the sterile or fertile flowers was extracted using the Trizol reagent (Invitrogen, USA) and mRNA was purified using the RNeasy Plant Mini Kit (Qiagen, Valencia, CA). RNA concertrations were measured using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA) and the integrity was determined using a Bioanalyzer 2100 (Aligent, CA). The RNA integrity number (RIN) values were greater than 8.5 for all samples. Sequencing libraries were prepared according to the manufacturer’s instructions (Illumina, CA). mRNA was subjected to two rounds purification using poly(T) oligo-attached magnetic beads, and fragmented using an RNA fragmentation kit. First strand cDNA was generated using reverse transcriptase and random primers. Following the second strand cDNA synthesis and adaptor ligation, 200 bp cDNA fragments were isolated using gel electrophoresis and amplified by 18 cycles of PCR. The products were loaded onto an Illumina HiSeq2000 instrument and subjected to paired-end sequencing. Images analysis, base-calling, raw tags generation, and tags counting were performed using the Illumina data processing pipeline.
Raw reads were trimmed to obtain high-quality reads by removing adaptors and low-quality tags. All clean reads were assembled and annotated against the following database: Nt (NCBI non-redundant necleotide sequences), Nr (NCBI non-redundant protein sequences), Pfam (Protein family), UniProtKB/Swiss-port (the UniProt Knowledgebase), KOG (euKaryotic Ortholog Groups), GO (Gene Ontology); and KO (KEGG Ortholog database). All BLAST searches were performed with the cut-off of E-value ≤1E− 5.
Differential gene expression, GO and KEGG enrichment analysis
To identify the differentially expressed genes (DEGs) between two samples, the number of expressed tags was calculated and then normalized in terms of RPKM (Mortazavi A, 2008). Genes with a absolute value of ∣log2 ratio∣ ≥ 1, False discovery rate (FDR) ≤ 0.05 and the adjusted P-value < 0.05 were assigned as differentially expressed. Differentially expressed transcripts were classified into functional categories using Blast2GO with default parameters , based on three levels of GO terms: biological processes (BP), molecular function (MF) and cell components (CC). GO and KEGG enrichment analysis of DEGs was performed using the GOseq R packages (http://www.r-project.org/) and KOBAS (http://www.kobas. co.uk/), respectively.
Validation of DEGs by qPCR
To verify the accuracy of RNA-seq results, 8 DEGs were selected for relative qPCR analysis. Total RNA was extracted, treated with DNase, and was reverse transcribed into cDNA using a RevertAid First Strand cDNA Synthesis kit (Fermentas, Vilnius, Lithuania). SYBR-based qPCR reactions were performed on a LightCycler 480 system (Roche, Switzerland) with the following reaction conditions: 95 °C for 1 min, followed by 40 cycles of 95 °C for 10 s, and 60 °C for 30 s. All PCR reactions were performed in triplicate, and the relative expression levels were calculated with the 2-ΔΔCT method  using Ubc9 and Ubc21 as the references. Gene-specific primers were listed in Additional file 8: Table S7.
Availability of data and materials
The sequence raw data from this study have been submitted to the NCBI Sequence Read Archive (SRR8398869) (https://www.ncbi.nlm.nih.gov/sra/?term=SRR8398869).
Phytohormone abscisic acid
Actin depolymerizing factor
Auxin response factor
Indole-3-acetic acid-induced protein
Auxin influx carrier
Calcium and calmodulin
Cytoplasmic male sterility
Calcium-dependent protein kinase
Differentially expressed genes
Indole-3-acetic acid-amido synthetase
Genic male sterility
Heat shock transcription factor
Heat shock proteins
Jasmonate acid carboxyl methyltransferase
Kyoto encyclopedia of genes and genomes
Eukaryotic orthologous groups
MADS-box transcription factor
Nuclear transcription factor Y
Programmed cell death
Rho of plant
Thermo-sensitive genic male sterility
Fan ZG, Stefansson BR. Influence of temperature on sterility of two cytoplasmic male-sterility system in rape (Brassica napus L.). Can J Plant Sci. 1986;66(2):221–7.
Yang LY, Liu PW, Yang GS. Development of polima temperature sensitive cytoplasmic male sterile lines of Brassica napus through isolated microspore culture. Plant Breed. 2010;125(4):368–71.
Cao SH, Zhang XQ, Zhang AM. Review of the molecular regulation mechanism and genetics of photoperiod- and/or thermosensitive male sterility. Chin Bull Bot. 2005;22(1):19–26.
Yang G, Fu T, Ma C, Yang X. Studies on the ecotypical male sterile line of Brassica napus L. environmental effect on the ecotypical male sterile line. J Huazhong Agric University. 1997;16(5):330–4.
Tang JH, Hu YM, Fu ZY, He ZY, He J, Chen WC. Development, identification and genetic analysis of a new thermo-sensitive genic male sterile line in maize. Sci Agric Sin. 2007;40(5):889–94.
Li LJ, Ru ZG, Gao QR, Jiang H, Guo FZ, Wu SW, Sun Z. Male sterility and thermo- photosensitivity characteristics of BNS in wheat. Sci Agric Sin. 2009;42(9):3019–27.
Fu QY, Cao YP, Li YY. Advanced on studies and applications of photo-thermo-sensitive male sterility in wheat. J Triticeae Crops. 2010;30(3):576–80.
Shi MS. The discovery and study of the photosensitive recessive male-sterile rice (Oryza sativa L. subsp. japonica). Sci Agric Sin. 1985;2:44–8.
Kim YJ, Zhang D. Molecular control of male fertility for crop hybrid breeding. Trends Plant Sci. 2018;23(1):53–65.
Zhou H, Zhou M, Yang Y, Li J, Zhu L, Jiang D, Dong J, Liu Q, Gu L, Zhou L, Feng M, Qin P, Hu X, Song C, Shi J, Song X, Ni E, Wu X, Deng Q, Liu Z, Chen M, Liu Y, Cao X, Zhuang C. RNase ZS1 processes UbL40 mRNAs and controls thermosensitive genic male sterility in rice. Nat Commun. 2014;5(1):4884.
Li H, Yuan Z, Vizcay-Barrena G, Yang C, Liang W, Zong J, Wilson ZA, Zhang D. Persistent tapetal cell1 encodes a PHD-finger protein that is required for tapetal cell death and pollen development in rice. Plant Physiol. 2011;156(2):615–30.
Qi Y, Liu Q, Zhang L, Mao B, Yan D, Jin Q, He Z. Fine mapping and candidate gene analysis of the novel thermo-sensitive genic male sterility tms9-1 gene in rice. Theor Appl Genet. 2014;127(5):1173–82.
Bai B, Wu J, Sheng WT, Zhou B, Zhou LJ, Zhuang W, Yao DP, Deng QY. Comparative analysis of anther transcriptome profiles of two different rice male sterile lines genotypes under cold stress. Int J Mol Sci. 2015;16(5):11398–416.
Hu J, Chen X, Zhang H, Ding Y. Genome-wide analysis of DNA methylation in photoperiod- and thermo-sensitive male sterile rice Peiai 64S. BMC Genomics. 2015;16(1):102.
Chen X, Hu J, Zhang H, Ding Y. DNA methylation changes in photoperiod-thermo-sensitive male sterile rice PA64S under two different conditions. Gene. 2014;537(1):143–8.
Kim HA, Lim CJ, Kim S, Choe JK, Jo SH, Baek N, Kwon SY. High-throughput sequencing and de novo assembly of Brassica oleracea var. capitata L. for transcriptome analysis. PLoS One. 2014;9(3):e92087.
Ma Y, Kang J, Wu J, Zhu Y, Wang X. Identification of tapetum-specific genes by comparing global gene expression of four different male sterile lines in Brassica oleracea. Plant Mol Biol. 2015;87(6):541–54.
Dong XS, Feng H, Xu M, Lee J, Kim YK, Lim YP, Piao Z, Park YD, Ma H, Hur Y. Comprehensive analysis of genic male sterility-related genes in Brassica rapa using a newly developed Br300K oligomeric chip. PLoS One. 2013;8(9):e72178.
Paritosh K, Gupta V, Yadava SK, Singh P, Pradhan AK, Pental D. RNA-seq based SNPs for mapping in Brassica juncea (AABB): synteny analysis between the two constituent genomes A (from B. rapa) and B (from B. nigra) shows highly divergent gene block arrangement and unique block fragmentation patterns. BMC Genomics. 2014;15(1):396.
Qu C, Fu F, Liu M, Zhao H, Liu C, Li J, Tang Z, Xu X, Qiu X, Wang R, Lu K. Comparative transcriptome analysis of recessive male sterility (RGMS) in sterile and fertile Brassica napus lines. PLoS One. 2015;10(12):e0144118.
An H, Yang Z, Yi B, Wen J, Shen J, Tu J, Ma C, Fu T. Comparative transcript profiling of the fertile and sterile flower buds of pol CMS in B. napus. BMC Genomics. 2014;15(1):258.
Liu XQ, Liu ZQ, Yu CY, Dong JG, Hu SW, Xu AX. TGMS in rapeseed (Brassica napus) resulted in aberrant transcriptional regulation, asynchronous microsporocyte meiosis, defective tapetum, and fused sexine. Fronti Plant Sci. 2017;8:1268.
Zhang T, Shen LY, Wang RX, Zou Y, Zhao JH, Li RC, Liang JL, Gong HM. Fertility alteration and utilization of male-sterile line 160S in Brassica napus. Acta Botan Boreali-Occiden Sin. 2012;32(1):35–41.
Tashiro T, Wardlaw I. The response to high temperature shock and humidity changes prior to and during the early stages of grain development in wheat. Funct Plant Biol. 1990;17(5):551–61.
Omidi M, Siahpoosh MR, Mamghani R, Modarresi M. The influence of terminal heat stress on meiosis abnormalities in pollen mother cells of wheat. Cytologia. 2014;79(1):49–58.
Jagadish S, Craufurd P, Wheeler T. High temperature stress and spikelet fertility in rice (Oryza sativa L.). J Exp Bot. 2007;58(7):1627–35.
Yang KZ, Xia C, Liu XL, Dou XY, Wang W, Chen LQ, Zhang XQ, Xie LF, He L, Ma X, Ye D. A mutation in thermosensitive male sterile 1, encoding a heat shock protein with DnaJ and PDI domains, leads to thermosensitive gametophytic male sterility in Arabidopsis. Plant J. 2009;57(5):870–82.
Perez DE, Hoyer JS, Johnson AI, Moody ZR, Lopez J, Kaplinsky NJ. BOBBER1 is a noncanonical Arabidopsis small heat shock protein required for both development and thermotolerance. Plant Physiol. 2009;151(1):241–52.
Dupuis I, Dumas C. Influence of temperature stress on in vitro fertilization and heat shock protein synthesis in maize (Zea mays L.) reproductive tissues. Plant Physiol. 1990;94(2):665–70.
Hopf N, Plesofsky-Vig N, Brambl R. The heat shock response of pollen and other tissues of maize. Plant Mol Boil. 1992;19(4):623–30.
Young TE, Ling J, Geisler-Lee CJ, Tanguay RL, Caldwell C, Gallie DR. Developmental and thermal regulation of the maize heat shock protein, HSP101. Plant Physiol. 2001;127(3):777–91.
Kim M, Kim H, Lee W, Lee Y, Kwon SW, Lee J. Quantitative shotgun proteomics analysis of rice anther proteins after exposure to high temperature. Int J Genomics. 2015;2015:238704.
Xie HT, Wan ZY, Li S, Zhang Y. Spatiotemporal production of reactive oxygen species by NADPH oxidase is critical for tapetal programmed cell death and pollen development in Arabidopsis. Plant Cell. 2014;26(5):2007–23.
Yu SX, Feng QN, Xie HT, Li S, Zhang Y. Reactive oxygen species mediate tapetal programmed cell death in tobacco and tomato. BMC Plant Biol. 2017;17(1):76.
Yi J, Moon S, Lee YS, Zhu L, Liang W, Zhang D, Jung KH, An G. Defective tapetum cell death 1 (DTC1) regulates ROS levels by binding to metallothionein during tapetum degeneration. Plant Physiol. 2016;170(3):1611–23.
Smith AR, Zhao D. Sterility caused by floral organ degeneration and abiotic stresses in Arabidopsis and cereal grains. Front Plant Sci. 2016;7:1503.
Higashitani A. High temperature injury and auxin biosynthesis in microsporogenesis. Front Plant Sci. 2013;4:47.
Sakata T, Oshino T, Miura S, Tomabechi M, Tsunaga Y, Higashitani N, Miyazawa Y, Takahashi H, Watanabe M, Higashitani A. Auxins reverse plant male sterility caused by high temperatures. Proc Natl Acad Sci U S A. 2010;107(19):8569–74.
Nagpal P, Ellis CM, Weber H, Ploense SE, Barkawi LS, Guilfoyle TJ, Hagen G, Alonso JM, Cohen JD, Farmer EE. Auxin response factors ARF6 and ARF8 promote jasmonic acid production and flower maturation. Development. 2005;132(18):4107–18.
Wu MF, Tian Q, Reed JW. Arabidopsis microRNA167 controls patterns of ARF6 and ARF8 expression, and regulates both female and male reproduction. Development. 2006;133(21):4211–8.
Sessions A, Nemhauser JL, McColl A, Roe JL, Feldmann KA, Zambryski PC. ETTIN patterns the Arabidopsis floral meristem and reproductive organs. Development. 1997;124(22):4481–91.
Liu Z, Miao L, Huo R, Song X, Cameron J, Kong L, Sundaresan V, Yu X. ARF2-ARF4 and ARF5 are essential for female and male gametophyte development in Arabidopsis. Plant Cell Physiol. 2018;59(1):179–89.
Cao D, Cheng H, Wu W, Soo HM, Peng J. Gibberellin mobilizes distinct DELLA-dependent transcriptomes to regulate seed germination and floral development in Arabidopsis. Plant Physiol. 2006;142(2):509–25.
Hou X, Hu WW, Shen L, Lee LY, Tao Z, Han JH, Yu H. Global identification of DELLA target genes during Arabidopsis flower development. Plant Physiol. 2008;147(3):1126–42.
Plackett AR, Ferguson AC, Powers SJ, Wanchoo-Kohli A, Phillips AL, Wilson ZA, Hedden P, Thomas SG. DELLA activity is required for successful pollen development in the Columbia ecotype of Arabidopsis. New Phytol. 2014;201(3):825–36.
Cheng H, Song S, Xiao L, Soo HM, Cheng Z, Xie D, Peng J. Gibberellin acts through jasmonate to control the expression of MYB21, MYB24, and MYB57 to promote stamen filament growth in Arabidopsis. PLoS Genet. 2009;5(3):e1000440.
Cheng Z, Sun L, Qi T, Zhang B, Peng W, Liu Y, Xie D. The bHLH transcription factor MYC3 interacts with the jasmonate ZIM-domain proteins to mediate jasmonate response in Arabidopsis. Mol Plant. 2011;4(2):279–88.
Fernández-Calvo P, Chini A, Fernández-Barbero G, Chico J-M, Gimenez-Ibanez S, Geerinck J, Eeckhout D, Schweizer F, Godoy M, Franco-Zorrilla JM, Pauwels L, Witters E, Puga MI, Paz-Ares J, Goossens A, Reymond P, De Jaeger G, Solano R. The Arabidopsis bHLH transcription factors MYC3 and MYC4 are targets of JAZ repressors and act additively with MYC2 in the activation of jasmonate responses. Plant Cell. 2011;23(2):701–15.
Ishiguro S, Kawai-Oda A, Ueda J, Nishida I, Okada K. The DEFECTIVE IN ANTHER DEHISCENCE1 gene encodes a novel phospholipase A1 catalyzing the initial step of jasmonic acid biosynthesis, which synchronizes pollen maturation, anther DEHISCENCE, and flower opening in Arabidopsis. Plant Cell. 2001;13(10):2191–209.
Caldelari D, Wang G, Farmer EE, Dong X. Arabidopsis lox3 lox4 double mutants are male sterile and defective in global proliferative arrest. Plant Mol Biol. 2011;75(1–2):25–33.
Stintzi A, Browse J. The Arabidopsis male-sterile mutant, opr3, lacks the 12-oxophytodienoic acid reductase required for jasmonate synthesis. Proc Natl Acad Sci U S A. 2000;97(19):10625–30.
Yan Y, Christensen S, Isakeit T, Engelberth J, Meeley R, Hayward A, Emery RJ, Kolomiets MV. Disruption of OPR7 and OPR8 reveals the versatile functions of jasmonic acid in maize development and defense. Plant Cell. 2012;24(4):1420–36.
Ji X, Dong B, Shiran B, Talbot MJ, Edlington JE, Hughes T, White RG, Gubler F, Dolferus R. Control of abscisic acid catabolism and abscisic acid homeostasis is important for reproductive stage stress tolerance in cereals. Plant Physiol. 2011;156(2):647–62.
Oliver SN, Van Dongen JT, Alfred SC, Mamun EA, Zhao X, Saini HS, Fernandes SF, Blanchard CL, Sutton BG, Geigenberger P. Cold-induced repression of the rice anther-specific cell wall invertase gene OSINV4 is correlated with sucrose accumulation and pollen sterility. Plant Cell Environ. 2005;28(12):1534–51.
Ye Q, Zhu W, Li L, Zhang S, Yin Y, Ma H, Wang X. Brassinosterioids control male fertility by regulating the expression of key genes involved in Arabidopsis anther and pollen development. Proc Natl Acad Sci U S A. 2010;107(13):6100–5.
Yin Y, Wang ZY, Mora-Garcia S, Li J, Yoshida S, Asami T, Chory J. BES1 accumulates in the necleus in response to brassinosteroids to regulate gene expression and promote stem elongation. Cell. 2002;109(2):181–91.
Wilson ZA, Zhang DB. From Arabidopsis to rice: pathways in pollen development. J Exp Bot. 2009;60(5):1479–92.
Liu Y, Cui S, Wu F, Yan S, Lin X, Du X, Chong K, Schilling S, Theißen G, Meng Z. Functional conservation of MIKC*-type MADS box genes in Arabidopsis and rice pollen maturation. Plant Cell. 2013;25(4):1288–303.
Becker A, Theissen G. The major clades of MADS-box genes and their role in the development and evolution of flowering plants. Mol Phylogenet Evol. 2003;29(3):464–89.
Hu L, Liang W, Yin C, Cui X, Zong J, Wang X, Hu J, Zhang D. Rice MADS3 regulates ROS homeostasis during late anther development. Plant Cell. 2011;23(2):515–33.
Song S, Qi T, Huang H, Ren Q, Wu D, Chang C, Peng W, Liu Y, Peng J, Xie D. The jasmonate-ZIM domain proteins interact with the R2R3-MYB transcription factors MYB21 and MYB24 to affect jasmonate-regulated stamen development in Arabidopsis. Plant Cell. 2011;23(3):1000–13.
Phan HA, Li SF, Parish RW. MYB80, a regulator of tapetal and pollen development, is functionally conserved in crops. Plant Mol Biol. 2012;78(1–2):171–83.
Xu C, Liu Z, Zhang L, Zhao C, Yuan S, Zhang F. Organization of actin cytoskeleton during meiosis I in a wheat thermo-sensitive genic male sterile line. Protoplasma. 2013;250(1):415–22.
Xie CT, Wei DM, Tian HQ. Advances in cell biological researches on male sterility of higher plants. J Plant Physiol Mol Biol. 2006;32(1):17–23.
Dong CH, Xia GX, Hong Y, Ramachandran S, Kost B, Chua NH. ADF proteins are involved in the control of flowering and regulate F-actin organization, cell expansion, and organ growth in Arabidopsis. Plant Cell. 2001;13(6):1333–46.
Zhang J, Zhu QD, Ju L, Zhang GS, Yu YA, Niu N, Wang JW, Ma SC. Organization changes of actin cytoskeleton and callose with the expression analysis of related genes in physiological male sterile wheat. Sci Agric Sin. 2015;48(14):2687–96.
Zheng ZL, Nafisi M, Tam A, Li H, Crowell DN, Chary SN, Schroeder JI, Shen J, Yang Z. Plasma membrane-associated ROP10 small GTPase is a specific negative regulator of abscisic acid responses in Arabidopsis. Plant Cell. 2002;14(11):2787–97.
Baxter-Burrell A, Yang Z, Springer PS, Bailey-Serres J. Rop GAP4-dependent Rop GTPase rheostat control of Arabidopsis oxygen deprivation tolerance. Science. 2002;296(5575):2026–8.
Schulthesis H, Dechert C, Kogel KH, Hückelhoven R. Functional analysis of barley RAC/ROP G-protein family members in susceptibility to the powdery mildew fungus. Plant J. 2003;36(5):589–601.
Basu D, Le J, Zakharova T, Mallery EL, Szymanski DB. A SPIKE1 signaling complex controls actin-dependent cell morphogenesis through the heteromeric WAVE and ARP2/3 complexes. Proc Natl Acad Sci U S A. 2008;105(10):4044–9.
Yang Z. Cell polarity signaling in Arabidopsis. Annu Rev Cell Dev Biol. 2008;24(1):551–75.
Steinhorst L, Kudla J. Calcium-a central regulator of pollen germination and tube growth. BBA-Mol Cell Res. 2013;1833(7):1573–81.
Chao Q, Gao ZF, Wang YF, Li Z, Huang XH, Wang YC, Mei YC, Zhao BG, Li L, Jiang YB, Wang BC. The proteome and phosphoproteome of maize pollen uncovers fertility candidate proteins. Plant Mol Biol. 2016;91(3):287–304.
Li J, Li Y, Deng Y, Chen P, Feng F, Chen W, Zhou X, Wang Y. A calcium-dependent protein kinase, ZmCPK32, specifically expressed in maize pollen to regulate pollen tube growth. PLoS One. 2018;13(5):e0195787.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
Authors thank the Chongqing Engineering Research Center of Special Crop Resources for scientific advice.
This study was founded by the National Natural Science Foundation of China (31171588) and the project from Chongqing Science and Technology Bureau (cstc2016shmszx80051). The funding bodies had no roles in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Line 160S were screened by the corresponding author and kept in our lab. We confirm that no field experiments were involved in this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Sumary of the RNA-seq data. (XLS 19 kb)
Figure S1. Gene annotation statistics (A) and species classification (B). (PNG 1091 kb)
Table S2. Gene Ontology (GO) terms assigned to the unigenes. (XLSX 11 kb)
Table S3. KOG anotation of the mapped genes. (XLSX 9 kb)
Table S4. Functional annotation of DEGs based on KEGG categorization. (XLSX 20 kb)
Table S5. KEGG pathway enrichment analysis of differentially expressed genes between 160S-MF and 160S-MS. (XLS 50 kb)
Table S6. Differentially expressed genes putatively involved in TGMS of B. napus. (XLS 36 kb)
Table S7. Primers used for qRT-PCR. (XLSX 9 kb)