Transcriptome de novo assembly and analysis of differentially expressed genes related to cytoplasmic male sterility in onion

Repeat-mediated recombination results in Complex DNA structure of the mitochondrial genome of Trachelospermum jasminoides

BACKGROUND

Trachelospermum jasminoides has medicinal and ornamental value and is widely distributed in China. Although the chloroplast genome has been documented, the mitochondrial genome has not yet been studied.

RESULTS

The mitochondrial genome of T. jasminoides was assembled and functionally annotated using Illumina and nanopore reads. The mitochondrial genome comprises a master circular molecular structure of 605,764 bp and encodes 65 genes: 39 protein-coding genes, 23 transfer RNA (tRNA) genes and 3 ribosomal RNA genes. In addition to the single circular conformation, we found many alternative conformations of the T. jasminoides mitochondrial genome mediated by 42 repetitive sequences. Six repetitive sequences (DRS01-DRS06) were supported by nanopore long reads, polymerase chain reaction (PCR) amplifications, and Sanger sequencing of the PCR products. Eleven homologous fragments were identified by comparing the mitochondrial and chloroplast genome sequences, including three complete tRNA genes. Moreover, 531 edited RNA sites were identified in the protein-coding sequences based on RNA sequencing data, with nad4 having the highest number of sites (54).

CONCLUSION

To our knowledge, this is the first description of the mitochondrial genome of T. jasminoides. Our results demonstrate the existence of multiple conformations. These findings lay a foundation for understanding the genetics and evolutionary dynamics of Apocynaceae.

AlliumDB: a central portal for comparative and functional genomics in Allium

The genus Allium belongs to the botanical family Amaryllidaceae and includes economically important crops such as onion, garlic, bunching onion, and leek, used as vegetables, spices, and traditional medicines. The large sizes of Allium genomes hamper the genetic dissection of agronomically important traits and molecular breeding. With the growing accumulation of genomic, resequencing, transcriptome, and phenotypic data, the demand for an integrative Allium database is increasing. Here we present a user-friendly database, AlliumDB (https://allium.qau.edu.cn), as a functional genomics hub integrating public and in-house data. The database contains all currently available nuclear and organelle genomes for Allium species, with genes comprehensively annotated based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, orthology, gene families, protein families (Pfam), and non-coding RNA families (Rfam). Transcriptome and variation profiles are integrated into dynamic visualization tools. We took phenotypic photographs and generated trait records for hundreds of Allium germplasms collected worldwide, which are included in the database. We incorporated JBrowse for the visualization of gene structures, RNA sequencing data, and variation data. Analysis tools such as the basic local alignment search tool (BLAST), sequence fetch, enrichment, and motif analyses are available to explore potential gene functions. This database incorporates comprehensive Allium genotypic and phenotypic datasets. As the community assembles new genomes and generates resequencing data for Allium germplasms, the database will be improved and continuously updated with these multi-omics data and comparative genomic studies. We expect the AlliumDB database to become a key resource for the study of Allium crops.

Cloning and Expression Analysis of Onion (Allium cepa L.) MADS-Box Genes and Regulation Mechanism of Cytoplasmic Male Sterility

Flower organ development is one of the most important processes in plant life. However, onion CMS (cytoplasmic male sterility) shows an abnormal development of floral organs. The regulation of MADS-box transcription factors is important for flower development. To further understand the role of MADS-box transcription factors in the regulation of cytoplasmic male sterility onions. We cloned the full-length cDNA of five MADS-box transcription factors from the flowers of onion using RACE (rapid amplification of cDNA ends) technology. We used bioinformatics methods for sequence analysis and phylogenetic analysis. Real-time quantitative PCR was used to detect the expression patterns of these genes in different onion organs. The relative expression levels of five flower development genes were compared in CMS onions and wild onions. The results showed that the full-length cDNA sequences of the cloned MADS-box genes AcFUL, AcDEF, AcPI, AcAG, and AcSEP3 belonged to A, B, C, and E MADS-box genes, respectively. A phylogenetic tree construction analysis was performed on its sequence. Analysis of MADS-box gene expression in wild onion and CMS onion showed that the formation of CMS onion was caused by down-regulation of AcDEF, AcPI, and AcAG gene expression, up-regulation of AcSEP3 gene expression, and no correlation with AcFUL gene expression. This work laid the foundation for further study of the molecular mechanism of onion flower development and the molecular mechanism of CMS onion male sterility.

Comparative Transcriptome Analysis Reveals a Potential Regulatory Network for Ogura Cytoplasmic Male Sterility in Cabbage (Brassica oleracea L.)

Ogura cytoplasmic male sterility (CMS) lines are widely used breeding materials in cruciferous crops and play important roles in heterosis utilization; however, the sterility mechanism remains unclear. To investigate the microspore development process and gene expression changes after the introduction of orf138 and Rfo, cytological observation and transcriptome analysis were performed using a maintainer line, an Ogura CMS line, and a restorer line. Semithin sections of microspores at different developmental stages showed that the degradation of tapetal cells began at the tetrad stage in the Ogura CMS line, while it occurred at the bicellular microspore stage to the tricellular microspore stage in the maintainer and restorer lines. Therefore, early degradation of tapetal cells may be the cause of pollen abortion. Transcriptome analysis results showed that a total of 1287 DEGs had consistent expression trends in the maintainer line and restorer line, but were significantly up- or down-regulated in the Ogura CMS line, indicating that they may be closely related to pollen abortion. Functional annotation showed that the 1287 core DEGs included a large number of genes related to pollen development, oxidative phosphorylation, carbohydrate, lipid, and protein metabolism. In addition, further verification elucidated that down-regulated expression of genes related to energy metabolism led to decreased ATP content and excessive ROS accumulation in the anthers of Ogura CMS. Based on these results, we propose a transcriptome-mediated induction and regulatory network for cabbage Ogura CMS. Our research provides new insights into the mechanism of pollen abortion and fertility restoration in Ogura CMS.

Conjunctive Analyses of BSA-Seq and BSR-Seq Unveil the Msβ-GAL and MsJMT as Key Candidate Genes for Cytoplasmic Male Sterility in Alfalfa (Medicago sativa L.)

Knowing the molecular mechanism of male sterility in alfalfa is important to utilize the heterosis more effectively. However, the molecular mechanisms of male sterility in alfalfa are still unclear. In this study, the bulked segregant analysis (BSA) and bulked segregant RNA-seq (BSR) were performed with F2 separation progeny to study the molecular mechanism of male sterility in alfalfa. The BSA-seq analysis was located in a candidate region on chromosome 5 containing 626 candidate genes which were associated with male sterility in alfalfa, while the BSR-seq analysis filtered seven candidate DEGs related to male sterility, and these candidate genes including EF-Tu, β-GAL, CESA, PHGDH, and JMT. The conjunctive analyses of BSR and BSA methods revealed that the genes of Msβ-GAL and MsJMT are the common detected candidate genes involved in male sterility in alfalfa. Our research provides a theory basis for further study of the molecular mechanism of male sterility in alfalfa and significant information for the genetic breeding of Medicago sativa.

Proteomic analysis of differential anther development from sterile/fertile lines in Capsicum annuum L

Background

Pepper (Capsicum annuum L.) is a major cash crop throughout the world. Male sterility is an important characteristic in crop species that leads to a failure to produce functional pollen, and it has crucial roles in agricultural breeding and the utilization of heterosis.

Objectives

In this study, we identified many crucial factors and important components in metabolic pathways in anther and pollen development, and elucidated the molecular mechanism related to pollen abortion in pepper.

Methods

Pepper pollen was observed at different stages to detect the characteristics associated with male sterility and fertility. The phytohormone and oxidoreductase activities were detected in spectrophotometric and redox reaction assays, respectively. Proteins were extracted from male sterile and fertile pepper lines, and identified by TMT/iTRAQ (tandem mass tags/isobaric tags for relative and absolute quantitation) and LC-MS/MS (liquid chromatograph-mass spectrometer) analysis. Differentially abundant proteins (DAPs) were analyzed based on Gene Ontology annotations and the Kyoto Encyclopedia of Genes and Genomes database according to |fold change)| > 1.3 and P value < 0.05. DAPs were quantified in the meiosis, tetrad, and binucleate stages by parallel reaction monitoring (PRM).

Results

In this study, we screened and identified one male sterile pepper line with abnormal cytological characteristics in terms of pollen development. The peroxidase and catalase enzyme activities were significantly reduced and increased, respectively, in the male sterile line compared with the male fertile line. Phytohormone analysis demonstrated that the gibberellin, jasmonic acid, and auxin contents changed by different extents in the male sterile pepper line. Proteome analysis screened 1,645 DAPs in six clusters, which were mainly associated with the chloroplast and cytoplasm based on their similar expression levels. According to proteome analysis, 45 DAPs were quantitatively identified in the meiosis, tetrad, and binucleate stages by PRM, which were related to monoterpenoid biosynthesis, and starch and sucrose metabolism pathways.

Conclusions

We screened 1,645 DAPs by proteomic analysis and 45 DAPs were related to anther and pollen development in a male sterile pepper line. In addition, the activities of peroxidase and catalase as well as the abundances of phytohormones such as gibberellin, jasmonic acid, and auxin were related to male sterility. The results obtained in this study provide insights into the molecular mechanism responsible for male sterility and fertility in pepper.

Using Transcriptome Analysis to Explore Gray Mold Resistance-Related Genes in Onion (Allium cepa L.)

Gray mold disease caused by Botrytis in onions (Allium cepa L.) during growth and storage negatively affects their yield and quality. Exploring the genes related to gray mold resistance in onion and their application to the breeding of resistant onion lines will support effective and ecological control methods of the disease. Here, the genetic relationship of 54 onion lines based on random amplified polymorphic DNA (RAPD) and in vitro-cultured onion lines infected with gray mold were used for screening resistance and susceptibility traits. Two genetically related onion lines were selected, one with a resistant and one with a susceptible phenotype. In vitro gray mold infection was repeated with these two lines, and leaf samples were collected for gene expression studies in time series. Transcript sequences obtained by RNA sequencing were subjected to DEG analysis, variant analysis, and KEGG mapping. Among the KEGG pathways, ‘α-linoleic acid metabolism’ was selected because the comparison of the time series expression pattern of Jasmonate resistant 1 (JAR1), Coronatine-insensitive protein 1 (COI 1), and transcription factor MYC2 (MYC2) genes between the resistant and susceptible lines revealed its significant relationship with gray-mold-resistant phenotypes. Expression pattern and SNP of the selected genes were verified by quantitative real-time PCR and high-resolution melting (HRM) analysis, respectively. The results of this study will be useful for the development of molecular marker and finally breeding of gray-mold-resistant onions.

Did doubly uniparental inheritance (DUI) of mtDNA originate as a cytoplasmic male sterility (CMS) system?

Animal and plant species exhibit an astonishing diversity of sexual systems, including environmental and genetic determinants of sex, with the latter including genetic material in the mitochondrial genome. In several hermaphroditic plants for example, sex is determined by an interaction between mitochondrial cytoplasmic male sterility (CMS) genes and nuclear restorer genes. Specifically, CMS involves aberrant mitochondrial genes that prevent pollen development and specific nuclear genes that restore it, leading to a mixture of female (male-sterile) and hermaphroditic individuals in the population (gynodioecy). Such a mitochondrial-nuclear sex determination system is thought to be rare outside plants. Here, we present one possible case of CMS in animals. We hypothesize that the only exception to the strict maternal mtDNA inheritance in animals, the doubly uniparental inheritance (DUI) system in bivalves, might have originated as a mitochondrial-nuclear sex-determination system. We document and explore similarities that exist between DUI and CMS, and we propose various ways to test our hypothesis.

Transcriptomic analysis reveals candidate genes for male sterility in Prunus sibirica

Background

The phenomenon of male sterility widely occurs in Prunus sibirica and has a serious negative impact on yield. We identified the key stage and cause of male sterility and found differentially expressed genes related to male sterility in Prunus sibirica, and we analyzed the expression pattern of these genes. This work aimed to provide valuable reference and theoretical basis for the study of reproductive development and the mechanisms of male sterility in Prunus sibirica.

Method

The microstructures of male sterile flower buds and male fertile flower buds were observed by paraffin section. Transcriptome sequencing was used to screen genes related to male sterility in Prunus sibirica. Quantitative real-time PCR analysis was performed to verify the transcriptome data.

Results

Anther development was divided into the sporogenous cell stage, tetrad stage, microspore stage, and pollen maturity stage. Compared with male fertile flower buds, in the microspore stage, the pollen sac wall tissue in the male sterile flower buds showed no signs of degeneration. In the pollen maturity stage, the tapetum and middle layer were not fully degraded, and anther development stopped. Therefore, the microspore stage was the key stage for anther abortion , and the pollen maturity stage was the post stage for anther abortion. A total of 4,108 differentially expressed genes were identified by transcriptome analysis. Among them, 1,899 were up-regulated, and 2,209 were down-regulated in the transcriptome of male sterile flower buds. We found that "protein kinase activity", "apoptosis process", "calcium binding", "cell death", "cytochrome c oxidase activity", "aspartate peptidase activity", "cysteine peptidase activity" and other biological pathways such as "starch and sucrose metabolism" and "proteasome" were closely related to male sterility in Prunus sibirica. A total of 331 key genes were preliminarily screened.

Conclusion

The occurrence of male sterility in Prunus sibirica involved many biological processes and metabolic pathways. According to the results of microstructure observations, related physiological indexes determination and transcriptome analysis, we reveal that the occurrence of male sterility in Prunus sibirica may be caused by abnormal metabolic processes such as the release of cytochrome c in the male sterile flower buds, the imbalance of the antioxidant system being destroyed, and the inability of macromolecular substances such as starch to be converted into soluble small molecules at the correct stage of reproductive development, resulting in energy loss. As a result, the tapetum cannot be fully degraded, thereby blocking anther development, which eventually led to the formation of male sterility.

A comprehensive integrated transcriptome and metabolome analyses to reveal key genes and essential metabolic pathways involved in CMS in kenaf

Numbers of critical genes and pathways were found from the levels of transcriptome and metabolome, which were useful information for understanding of kenaf CMS mechanism. Cytoplasmic male sterility (CMS) is a maternally inherited trait in higher plants that leads to the inability to produce or release functional pollen. However, there is lack of comprehensive studies to reveal the molecular basis of CMS occurrence in kenaf. Herein, we performed transcriptome and UPLC-MS-based metabolome analyses in the anthers of a CMS (UG93A) and its maintainer (UG93B) to sort out essential genes and metabolites responding to CMS in kenaf. Transcriptome characterized 7769 differentially expressed genes (DEGs) between these two materials, and pathway enrichment analysis indicated that these DEGs were involved mainly in pentose and glucuronate interconversions, starch and sucrose metabolism, taurine and hypotaurine metabolism. In the metabolome assay, a total of 116 significantly different metabolites (SDMs) were identified between the CMS and its maintainer line, and these SDMs were involved in eight KEGG pathways, including flavone and flavonol biosynthesis, glycerophospholipid metabolism, flavonoid biosynthesis, glycosylphosphatidylinositol-anchor biosynthesi. Integrated analyses of transcriptome and metabolome showed that 50 genes had strong correlation coefficient values (R² > 0.9) with ten metabolites enriched in six pathways; notably, most genes and metabolites of flavonoid biosynthesis pathways and flavone and flavonol biosynthesis pathways involved in flavonoids biosynthetic pathways were downregulated in CMS compared to those in maintainer. Taken together, the decreased accumulation of flavonoids resulted from the compromised biosynthesis pathways coupled with energy deficiency in the anthers may contribute largely to CMS in UG93A of kenaf.

Differentially Expressed Genes Shared by Two Distinct Cytoplasmic Male Sterility (CMS) Types of Silene vulgaris Suggest the Importance of Oxidative Stress in Pollen Abortion

Cytoplasmic male sterility (CMS), encoded by the interacting mitochondrial and nuclear genes, causes pollen abortion or non-viability. CMS is widely used in agriculture and extensively studied in crops. Much less is known about CMS in wild species. We performed a comparative transcriptomic analysis of male sterile and fertile individuals of Silene vulgaris, a model plant for the study of gynodioecy, to reveal the genes responsible for pollen abortion in this species. We used RNA-seq datasets previously employed for the analysis of mitochondrial and plastid transcriptomes of female and hermaphrodite flower buds, making it possible to compare the transcriptomes derived from three genomes in the same RNA specimen. We assembled de novo transcriptomes for two haplotypes of S. vulgaris and identified differentially expressed genes between the females and hermaphrodites, associated with stress response or pollen development. The gene for alternative oxidase was downregulated in females. The genetic pathways controlling CMS in S. vulgaris are similar to those in crops. The high number of the differentially expressed nuclear genes contrasts with the uniformity of organellar transcriptomes across genders, which suggests these pathways are evolutionarily conserved and that selective mechanisms may shield organellar transcription against changes in the cytoplasmic transcriptome.

Omics approaches in Allium research: Progress and way ahead

BACKGROUND

The genus Allium (Family: Amaryllidaceae) is an economically important group of crops cultivated worldwide for their use as a vegetable and spices. Alliums are also well known for their nutraceutical properties. Among alliums, onion, garlic, leek, and chives cultivated worldwide. Despite their substantial economic and medicinal importance, the genome sequence of any of the Allium is not available, probably due to their large genome sizes. Recently evolved omics technologies are highly efficient and robust in elucidating molecular mechanisms of several complex life processes in plants. Omics technologies, such as genomics, transcriptomics, proteomics, metabolomics, metagenomics, etc. have the potential to open new avenues in research and improvement of allium crops where genome sequence information is limited. A significant amount of data has been generated using these technologies for various Allium species; it will help in understanding the key traits in Allium crops such as flowering, bulb development, flavonoid biosynthesis, male sterility and stress tolerance at molecular and metabolite level. This information will ultimately assist us in speeding up the breeding in Allium crops.

METHOD

In the present review, major omics approaches, and their progress, as well as potential applications in Allium crops, could be discussed in detail.

RESULTS

Here, we have discussed the recent progress made in Allium research using omics technologies such as genomics, transcriptomics, micro RNAs, proteomics, metabolomics, and metagenomics. These omics interventions have been used in alliums for marker discovery, the study of the biotic and abiotic stress response, male sterility, organ development, flavonoid and bulb color, micro RNA discovery, and microbiome associated with Allium crops. Further, we also emphasized the integrated use of these omics platforms for a better understanding of the complex molecular mechanisms to speed up the breeding programs for better cultivars.

CONCLUSION

All the information and literature provided in the present review throws light on the progress and potential of omics platforms in the research of Allium crops. We also mentioned a few research areas in Allium crops that need to be explored using omics technologies to get more insight. Overall, alliums are an under-studied group of plants, and thus, there is tremendous scope and need for research in Allium species.

PCA via joint graph Laplacian and sparse constraint: Identification of differentially expressed genes and sample clustering on gene expression data

Background

In recent years, identification of differentially expressed genes and sample clustering have become hot topics in bioinformatics. Principal Component Analysis (PCA) is a widely used method in gene expression data. However, it has two limitations: first, the geometric structure hidden in data, e.g., pair-wise distance between data points, have not been explored. This information can facilitate sample clustering; second, the Principal Components (PCs) determined by PCA are dense, leading to hard interpretation. However, only a few of genes are related to the cancer. It is of great significance for the early diagnosis and treatment of cancer to identify a handful of the differentially expressed genes and find new cancer biomarkers.

Results

In this study, a new method gLSPCA is proposed to integrate both graph Laplacian and sparse constraint into PCA. gLSPCA on the one hand improves the clustering accuracy by exploring the internal geometric structure of the data, on the other hand identifies differentially expressed genes by imposing a sparsity constraint on the PCs.

Conclusions

Experiments of gLSPCA and its comparison with existing methods, including Z-SPCA, GPower, PathSPCA, SPCArt, gLPCA, are performed on real datasets of both pancreatic cancer (PAAD) and head & neck squamous carcinoma (HNSC). The results demonstrate that gLSPCA is effective in identifying differentially expressed genes and sample clustering. In addition, the applications of gLSPCA on these datasets provide several new clues for the exploration of causative factors of PAAD and HNSC.

Comprehensive analysis of Ogura cytoplasmic male sterility-related genes in turnip (Brassica rapa ssp. rapifera) using RNA sequencing analysis and bioinformatics

Ogura-type cytoplasmic male sterility (Ogura-CMS) has been widely used in the hybrid breeding industry for cruciferous vegetables. Turnip (Brassica rapa ssp. rapifera) is one of the most important local cruciferous vegetables in China, cultivated for its fleshy root as a flat disc. Here, morphological characteristics of an Ogura-CMS line ‘BY10-2A’ and its maintainer fertile (MF) line ‘BY10-2B’ of turnip were investigated. Ogura-CMS turnip showed a reduction in the size of the fleshy root, and had distinct defects in microspore development and tapetum degeneration during the transition from microspore mother cells to tetrads. Defective microspore production and premature tapetum degeneration during microgametogenesis resulted in short filaments and withered white anthers, leading to complete male sterility of the Ogura-CMS line. Additionally, the mechanism regulating Ogura-CMS in turnip was investigated using inflorescence transcriptome analyses of the Ogura-CMS and MF lines. The de novo assembly resulted in a total of 84,132 unigenes. Among them, 5,117 differentially expressed genes (DEGs) were identified, including 1,339 up- and 3,778 down-regulated genes in the Ogura-CMS line compared to the MF line. A number of functionally known members involved in anther development and microspore formation were addressed in our DEG pool, particularly genes regulating tapetum programmed cell death (PCD), and associated with pollen wall formation. Additionally, 185 novel genes were proposed to function in male organ development based on GO analyses, of which 26 DEGs were genotype-specifically expressed. Our research provides a comprehensive foundation for understanding anther development and the CMS mechanism in turnip.

Transcriptomic and Proteomic Analysis of Shaan2A Cytoplasmic Male Sterility and Its Maintainer Line in Brassica napus

Cytoplasmic male sterility (CMS) lines are widely used for hybrid production in Brassica napus. The Shaan2A CMS system is one of the most important in China and has been used for decades; however, the male sterility mechanism underlying Shaan2A CMS remains unknown. Here, we performed transcriptomic and proteomic analysis, combined with additional morphological observation, in the Shaan2A CMS. Sporogenous cells, endothecium, middle layer, and tapetum could not be clearly distinguished in Shaan2A anthers. Furthermore, Shaan2A anther chloroplasts contained fewer starch grains than those in Shaan2B (a near-isogenic line of Shaan2A), and the lamella structure of chloroplasts in Shaan2A anther wall cells was obviously aberrant. Transcriptomic analysis revealed differentially expressed genes (DEGs) mainly related to carbon metabolism, lipid and flavonoid metabolism, and the mitochondrial electron transport/ATP synthesis pathway. Proteomic results showed that differentially expressed proteins were mainly associated with carbohydrate metabolism, energy metabolism, and genetic information processing pathways. Importantly, nine gene ontology categories associated with anther and pollen development were enriched among down-regulated DEGs at the young bud (YB) stage, including microsporogenesis, sporopollenin biosynthetic process, and tapetal layer development. Additionally, 464 down-regulated transcription factor (TF) genes were identified at the YB stage, including some related to early anther differentiation such as SPOROCYTELESS (SPL, also named NOZZLE, NZZ), DYSFUNCTIONAL TAPETUM 1 (DYT1), MYB80 (formerly named MYB103), and ABORTED MICROSPORES (AMS). These results suggested that the sterility gene in the Shaan2A mitochondrion might suppress expression of these TF genes in the nucleus, affecting early anther development. Finally, we constructed an interaction network of candidate proteins based on integrative analysis. The present study provides new insights into the molecular mechanism of Shaan2A CMS in B. napus.

Comparative Transcriptome Analysis between Fertile and CMS Flower Buds in Wucai (Brassica campestris L.)

Background

Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) is a variant of nonheading Chinese cabbage (Brassica campestris L.), which is one of the major vegetables in China. Cytoplasmic male sterility (CMS) has been used for Wucai breeding in recent years. However, the underlying molecular mechanism of Wucai CMS remains unclear. In this study, the phenotypic and cytological features of Wucai CMS were observed by anatomical analysis, and a comparative transcriptome analysis was carried out to identify genes related to male sterility using Illumina RNA sequencing technology (RNA-Seq).

Results

Microscopic observation demonstrated that tapetum development was abnormal in the CMS line, which failed to produce fertile pollen. Bioinformatics analysis detected 4430 differentially expressed genes (DEGs) between the fertile and sterile flower buds. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to better understand the functions of these DEGs. Among the DEGs, 35 genes (53 DEGS) were implicated in anther and pollen development, and 11 genes were involved in pollen cell wall formation and modification; most of these showed downregulated expression in sterile buds. In addition, several genes related to tapetum development (A6, AMS, MS1, MYB39, and TSM1) and a few genes annotated to flowering (CO, AP3, VIN3, FLC, FT, and AGL) were detected and confirmed by qRT-PCR as being expressed at the meiosis, tetrad, and uninucleate microspore stages, thus implying possible roles in specifying or determining the fate and development of the tapetum, male gametophyte and stamen. Moreover, the top four largest transcription factor families (MYB, bHLH, NAC and WRKY) were analyzed, and most showed reduced expression in sterile buds. These differentially expressed transcription factors might result in abortion of pollen development in Wucai.

Conclusion

The present comparative transcriptome analysis suggested that many key genes and transcription factors involved in anther development show reduced gene expression patterns in the CMS line, which might contribute to male sterility in Wucai. This study provides valuable information for a better understanding of CMS molecular mechanisms and functional genome studies in Wucai.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5331-4) contains supplementary material, which is available to authorized users.