Microsporogenesis is an indispensable period for investigating microspore development and cytoplasmic male sterility (CMS) occurrence. Radish CMS line plays a critical role in elite F1 hybrid seed production and heterosis utilization. However, the molecular mechanisms of microspore development and CMS occurrence have not been thoroughly uncovered in radish. In this study, a comparative analysis of radish floral buds from a CMS line (NAU-WA) and its maintainer (NAU-WB) was conducted using next generation sequencing (NGS) technology. Digital gene expression (DGE) profiling revealed that 3504 genes were significantly differentially expressed between NAU-WA and NAU-WB library, among which 1910 were upregulated and 1594 were downregulated. Gene ontology (GO) analysis showed that these differentially expressed genes (DEGs) were mainly enriched in extracellular region, catalytic activity, and response to stimulus. KEGG enrichment analysis revealed that the DEGs were predominantly associated with flavonoid biosynthesis, glycolysis, and biosynthesis of secondary metabolites. Real-time quantitative PCR analysis showed that the expression profiles of 13 randomly selected DEGs were in high agreement with results from Illumina sequencing. Several candidate genes encoding ATP synthase, auxin response factor (ARF), transcription factors (TFs), chalcone synthase (CHS), and male sterility (MS) were responsible for microsporogenesis. Furthermore, a schematic diagram for functional interaction of DEGs from NAU-WA vs. NAU-WB library in radish plants was proposed. These results could provide new information on the dissection of the molecular mechanisms underlying microspore development and CMS occurrence in radish.