Skip to main content
Log in

Comparative transcriptome analysis provides insight into the important pathways and key genes related to the pollen abortion in the thermo-sensitive genic male sterile line 373S in Brassica napus L.

  • Original Article
  • Published:
Functional & Integrative Genomics Aims and scope Submit manuscript

Abstract

The thermo-sensitive genic male sterility (TGMS) system plays a key role in the production of two-line hybrids in rapeseed (Brassica napus). To uncover key cellular events and genetic regulation associated with TGMS, a combined study using cytological methods and RNA-sequencing analysis was conducted for the rapeseed TGMS line 373S. Cytological studies showed that microspore cytoplasm of 373S plants was condensed, the microspore nucleus was degraded at an early stage, the exine was irregular, and the tapetum developed abnormally, eventually leading to male sterility. RNA-sequencing analysis identified 430 differentially expressed genes (298 upregulated and 132 downregulated) between the fertile and sterile samples. Gene ontology analysis demonstrated that the most highly represented biological processes included sporopollenin biosynthetic process, pollen exine formation, and extracellular matrix assembly. Kyoto encyclopedia of genes and genomes analysis indicated that the enriched pathways included amino acid metabolism, carbohydrate metabolism, and lipid metabolism. Moreover, 26 transcript factors were identified, which may be associated with abnormal tapetum degeneration and exine formation. Subsequently, 19 key genes were selected, which are considered to regulate pollen development and even participate in pollen exine formation. Our results will provide important insight into the molecular mechanisms underlying TGMS in rapeseed.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5

Similar content being viewed by others

Data availability

The RNA-seq raw read data have been submitted in the Sequence Read Archive of the NCBI (accession number: PRJNA612938)

Code availability

Not applicable.

Abbreviations

CHA:

chemical hybridization agent

CMS:

cytoplasmic male sterility

DAPI:

4′,6-diamidino-2-phenylindole

DEGs:

differentially expressed genes

EMS:

ecological male sterility

FPKM:

fragments per kilobase of transcript per million mapped reads

GMS:

genic male sterility

GO:

Gene Ontology

KEGG:

Kyoto encyclopedia of genes and genomes

PMC:

pollen mother cell

PTGMS:

photoperiod and temperature/thermo- sensitive genic male sterility

P/TGMS:

photoperiod and/or temperature sensitive genic male sterility

RGMS:

recessive genetic male sterility

RIN:

RNA integrity number

RNA-seq:

RNA-sequencing

SI:

self-incompatibility

TFs:

transcript factors

TGMS:

thermo-sensitive genic male sterility

ZS9:

Zhongshuang No. 9

References

Download references

Acknowledgements

The authors thank Prof. Huixian Zhao, College of Life Sciences, Northwest A&F University, for her suggestions in the course of preparing the manuscript. We are also grateful to our anonymous reviewers for their thoughtful and detailed reviews.

Funding

This research was funded by the National Natural Science Foundation of China (3217150949), the Key Research and Development Project in Shaanxi Province of China (2018ZDXM-NY-008; 2022NY-153), Modern Crop Seed Industry Project of Shaanxi Province (20171010000004), and the Key Research and Development Projects of Yangling Seed Industry Innovation Center (YLzy-yc2021-01).

Author information

Authors and Affiliations

Authors

Contributions

SH conceived and designed research. YS and DZ performed the experiments. YS and SH analyzed the data. YS and SH wrote the manuscript. HD, ZW, JW, HL and YG discussed and edited the manuscript. All the authors read and approved the final manuscript.

Corresponding author

Correspondence to Shengwu Hu.

Ethics declarations

Ethics approval

Not applicable.

Consent to participate

Not applicable.

Consent for publication

Not applicable.

Conflict of interest

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

The list of differentially expressed genes in pairwise comparisons between the fertile and sterile parent, as well as between fertile and sterile individuals in BC1 in two stages (bud and anther) (XLSX 1074 kb)

Additional file 2:

The 430 differentially expressed genes and their expression pattern in different stages (XLSX 123 kb)

Additional file 3:

Annotation of the 430 differentially expressed genes (XLSX 42 kb)

Additional file 4:

Gene ontology (GO) terms of each cluster in the heatmap (XLSX 261 kb)

Additional file 5:

Results of gene ontology term enrichment of the 430 differentially expressed genes (XLSX 18 kb)

Additional file 6:

The known and predicted interactions among the differentially expressed genes by String (XLSX 10 kb)

Additional file 7:

Cis-regulatory element analysis of the 430 differentially expressed genes (XLSX 1092 kb)

Additional file 8:

RNA-seq of DEGs in 373S under restrictive condition (high temperature) and permissive condition (low temperature) (XLSX 50 kb)

Additional file 9:

Fig. S1. Analysis of the relationship among all 15 samples according to RNA-seq based gene expression values. (A) The heatmap of sample correlation. The bar in right represents the scale of the relationship among samples, and the value in each pane represents the correlation coefficient between two samples. (B) Principal components analysis. All samples were clustered into two separate groups corresponding to developmental stages, buds and anthers shown by cyan and hermosa pink colors, respectively. The numbers in parentheses after PC (principal component) represent the proportion of variance explained by that PC. (C) The variance explained by the first four PCs. (PNG 298 kb)

Additional file 10:

Fig. S2. The pair-wise comparisons between fertile and sterile samples. The criteria for screening DEGs were probability≥0.8 and log2(fold change)≥1. 5A219F-A, anthers from fertile plants in BC1 (5A214// 5A214/5C95); 5A219F-B, buds from fertile plants in BC1; 5A219S-A, anthers from sterile plants in BC1; 5A219S-B, buds from sterile plants in BC1; 5C95-A, anthers from the fertile parent; 5C95-B, buds from the fertile parent; 5A214-A, anthers from the sterile parent; 5A214-B, buds from the sterile parent. (PNG 91 kb)

High Resolution Image (TIF 482 kb)

Additional file 11:

Fig. S3. Venn diagrams showing distribution of 430 differentially expressed genes between male fertile and sterile samples (PNG 50 kb)

High Resolution Image (TIF 2670 kb)

Additional file 12:

Fig. S4. Classification of gene ontology annotations. The x-axis indicates categories and sub-categories, and the y-axis the number of genes. (PNG 395 kb)

High Resolution Image (TIF 10705 kb)

Additional file 13:

Fig. S5. A network of differentially regulated gene ontology (GO) terms in biological processes. The scheme was drawn by using the BiNGO plugin of Cytoscape software. GO terms are presented as nodes and associations as edges, node size represents gene number, the bigger of the node size representing the more gene. The significant level of nodes is shown by different yellow colors as indicated in the scale bar. The lines with arrows mean the regulating relation. (PNG 491 kb)

High Resolution Image (TIF 9129 kb)

Additional file 14:

Fig. S6. Gene ontology functional enrichment of differential expressed genes in biological process. Coloring indicates p-value (high: yellow, low: red). The lower p-value indicates the more significant enrichment. The lines with arrows mean the regulating relation. (PNG 276 kb)

High Resolution Image (TIF 367 kb)

Additional file 15:

Fig. S7. A network of protein-protein interactions based on analysis of the differentially expressed genes. The scheme was drawn using the STRING software. Genes are shown as nodes and interactions as edges (PNG 2745 kb)

High Resolution Image (TIF 7101 kb)

Additional file 16:

Fig. S8. The most enriched pathways of the differentially expressed genes (DEGs) between sterile and fertile samples. (A) The most enriched metabolic activity between the sterile and fertile samples. (B) Pathway enrichment of the DEGs. The y-axis shows significantly enriched kyoto encyclopedia of genes and genomes pathways, and the x-axis shows the rich factor. Rich factor stands for the ratio of the number of DEGs belonging to a pathway to the number of all the annotated genes located in the pathway. The higher rich factor represents the higher level of enrichment. The size of the dot indicates the number of DEGs in the pathway, and the color of the dot reflects the different p-value range. (PNG 412 kb)

High Resolution Image (TIF 31418 kb)

Additional file 17:

Fig. S9. Cluster analysis of 19 differentially expressed genes associated with pollen development. The color key represents log10 (FPKM+1). Color intensities (from green to red shading) increase with elevated expression level, as indicated at upper left. (PNG 207 kb)

High Resolution Image (TIF 14630 kb)

Additional file 18:

Table S1. Description of the samples in the present study (PDF 39 kb)

Additional file 19:

Table S2. Target genes and their primer sequences for qRT-PCR (PDF 22 kb)

Additional file 20:

Table S3. Tetrad size and callose thickness in fertile control and 373S male sterile plants (PDF 104 kb)

Additional file 21:

Table S4. The statistics of RNA-sequencing data in individual samples (PDF 48 kb)

Additional file 22

Table S5 Verification of 430 DEGs identified between sterile samples and fertile samples by the RNA-seq experiment of 373S line under restrictive condition (high temperature) and permissive condition (low temperature) (PDF 106 kb)

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sun, Y., Zhang, D., Dong, H. et al. Comparative transcriptome analysis provides insight into the important pathways and key genes related to the pollen abortion in the thermo-sensitive genic male sterile line 373S in Brassica napus L.. Funct Integr Genomics 23, 26 (2023). https://doi.org/10.1007/s10142-022-00943-8

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1007/s10142-022-00943-8

Keywords

Navigation