Comparative transcript profiling of fertile and sterile flower buds from multiple-allele-inherited male sterility in Chinese cabbage (Brassica campestris L. ssp. pekinensis)

Lipid and sugar metabolism play an essential role in pollen development and male sterility: a case analysis in Brassica napus

AIMS

The genic male sterility (GMS) system is an important strategy for generating heterosis in plants. To better understand the essential role of lipid and sugar metabolism and to identify additional candidates for pollen development and male sterility, transcriptome and metabolome analysis of a GMS line of 1205AB in B. napus was used as a case study.

DATA RESOURCES GENERATED

To characterize the GMS system, the transcriptome and metabolome profiles were generated for 24 samples and 48 samples of 1205AB in B. napus, respectively. Transcriptome analysis yielded a total of 156.52 Gb of clean data and revealed the expression levels of 109,541 genes and 8,501 novel genes. In addition, a total of 1,353 metabolites were detected in the metabolomic analysis, including 784 in positive ion mode and 569 in negative ion mode.

KEY RESULTS

A total of 15,635 differentially expressed genes (DEGs) and 83 differential metabolites (DMs) were identified from different comparison groups, most of which were involved in lipid and sugar metabolism. The combination of transcriptome and metabolome analysis revealed 49 orthologous GMS genes related to lipid metabolism and 46 orthologous GMS genes related to sugar metabolism, as well as 45 novel genes.

UTILITY OF THE RESOURCE

The transcriptome and metabolome profiles and their analysis provide useful reference data for the future discovery of additional GMS genes and the development of more robust male sterility breeding systems for use in the production of plant hybrids.

Mutations in BrABCG26, encoding an ATP-binding cassette transporter, are responsible for male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

KEY MESSAGE

The gene BrABCG26 responsible for male sterility of Chinese cabbage was confirmed by two allelic mutants. Male-sterile lines are an important way of heterosis utilization in Chinese cabbage. In this study, two allelic male-sterile mutants msm3-1 and msm3-2 were obtained from a Chinese cabbage double haploid (DH) line 'FT' by using EMS-mutagenesis. Compared to the wild-type 'FT,' the stamens of mutants were completely degenerated and had no pollen, and other characters had no obvious differences. Cytological observation revealed that the failure of vacuolation of the mononuclear microspore, accompanied by abnormal tapetal degradation, resulted in anther abortion in mutants. Genetic analysis showed that a recessive gene controlled the mutant trait. MutMap combined with kompetitive allele specific PCR genotyping analyses showed that BraA01g038270.3C, encoding a transporter ABCG26 that played a vital role in pollen wall formation, was the candidate gene for msm3-1, named BrABCG26. Compared with wild-type 'FT,' the mutations existed on the second exon (C to T) and the sixth exon (C to T) of BrABCG26 gene in mutants msm3-1 and msm3-2, leading to the loss-of-function truncated protein, which verified the BrABCG26 function in stamen development. Subcellular localization and expression pattern analysis indicated that BrABCG26 was localized in the nucleus and was expressed in all organs, with the highest expression in flower buds. Compared to the wild-type 'FT,' the expressions of BrABCG26 were significantly reduced in flower buds and anthers of mutants. Promoter activity analysis showed that a strong GUS signal was detected in flower buds. These results indicated that BrABCG26 is responsible for the male sterility of msm3 mutants in Chinese cabbage.

Unveiling molecular mechanisms of pepper resistance to Phytophthora capsici through grafting using iTRAQ-based proteomic analysis

Phytophthora blight severely threatens global pepper production. Grafting bolsters plant disease resistance, but the underlying molecular mechanisms remain unclear. In this study, we used P. capsici-resistant strain 'ZCM334' and susceptible strain 'Early Calwonder' for grafting. Compared to self-rooted 'Early Calwonder' plants, 'ZCM334' grafts exhibited delayed disease onset, elevated resistance, and reduced leaf cell damage, showcasing the potential of grafting in enhancing pepper resistance to P. capsici. Proteomic analysis via the iTRAQ technology unveiled 478 and 349 differentially expressed proteins (DEPs) in the leaves and roots, respectively, between the grafts and self-rooted plants. These DEPs were linked to metabolism and cellular processes, stimulus responses, and catalytic activity and were significantly enriched in the biosynthesis of secondary metabolites, carbon fixation in photosynthetic organizations, and pyruvate metabolism pathways. Twelve DEPs exhibiting consistent expression trends in both leaves and roots, including seven related to P. capsici resistance, were screened. qRT-PCR analysis confirmed a significant correlation between the protein and transcript levels of DEPs after P. capsici inoculation. This study highlights the molecular mechanisms whereby grafting enhances pepper resistance to Phytophthora blight. Identification of key genes provides a foundation for studying the regulatory network governing the resistance of pepper to P. capsici.

Identification of NPF Family Genes in Brassica rapa Reveal Their Potential Functions in Pollen Development and Response to Low Nitrate Stress

Nitrate Transporter 1/Peptide Transporter Family (NPF) genes encode membrane transporters involved in the transport of diverse substrates. However, little is known about the diversity and functions of NPFs in Brassica rapa. In this study, 85 NPFs were identified in B. rapa (BrNPFs) which comprised eight subfamilies. Gene structure and conserved motif analysis suggested that BrNFPs were conserved throughout the genus. Stress and hormone-responsive cis-acting elements and transcription factor binding sites were identified in BrNPF promoters. Syntenic analysis suggested that tandem duplication contributed to the expansion of BrNPFs in B. rapa. Transcriptomic profiling analysis indicated that BrNPF2.6, BrNPF2.15, BrNPF7.6, and BrNPF8.9 were expressed in fertile floral buds, suggesting important roles in pollen development. Thirty-nine BrNPFs were responsive to low nitrate availability in shoots or roots. BrNPF2.10, BrNPF2.19, BrNPF2.3, BrNPF5.12, BrNPF5.16, BrNPF5.8, and BrNPF6.3 were only up-regulated in roots under low nitrate conditions, indicating that they play positive roles in nitrate absorption. Furthermore, many genes were identified in contrasting genotypes that responded to vernalization and clubroot disease. Our results increase understanding of BrNPFs as candidate genes for genetic improvement studies of B. rapa to promote low nitrate availability tolerance and for generating sterile male lines based on gene editing methods.

BrACOS5 mutations induced male sterility via impeding pollen exine formation in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BrACOS5 mutations led to male sterility of Chinese cabbage verified in three allelic male-sterile mutants. Chinese cabbage (Brassica rapa L. ssp. pekinensis) is one of the major vegetable crops in East Asia, and the utilization of male-sterile line is an important measure for its hybrid seed production. Herein, we isolated three allelic male-sterile mutants, msm1-1, msm1-2 and msm1-3, from an ethyl methane sulfonate (EMS) mutagenized population of Chinese cabbage double-haploid (DH) line 'FT', whose microspores were completely aborted with severely absent exine, and tapetums were abnormally developed. Genetic analyses indicated that the three male-sterile mutants belonged to allelic mutation and were triggered by the same recessive nuclear gene. MutMap-based gene mapping and kompetitive allele-specific PCR (KASP) analysis demonstrated that three different single-nucleotide polymorphisms (SNPs) of BraA09g012710.3C were responsible for the male sterility of msm1-1/2/3, respectively. BraA09g012710.3C is orthologous of Arabidopsis thaliana ACOS5 (AT1G62940), encoding an acyl-CoA synthetase in sporopollenin biosynthesis, and specifically expressed in anther, so we named BraA09g012710.3C as BrACOS5. BrACOS5 localizes to the endoplasmic reticulum (ER). Mutations of BrACOS5 resulted in decreased enzyme activities and altered fatty acid contents in msm1 anthers. As well as the transcript accumulations of putative orthologs involved in sporopollenin biosynthesis were significantly down-regulated excluding BrPKSA. These results provide strong evidence for the integral role of BrACOS5 in conserved sporopollenin biosynthesis pathway and also contribute to uncovering exine development pattern and underlying male sterility mechanism in Chinese cabbage.

Mutation in BrGGL7 gene encoding a GDSL esterase / lipase causes male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

MutMap and KASP analyses revealed that the BrGGL7 gene is responsible for the male-sterile trait of ftms1 in Chinese cabbage, with functional verification in Arabidopsis. The application of a male-sterile line is an ideal approach of hybrid seed production in Chinese cabbage. In this study, we obtained a male-sterile mutant (ftms1) from the double haploid line 'FT' using ethyl methane sulfonate (EMS) mutagenesis. The mutant was completely sterile due to abnormal enlargement and vacuolization of the tapetum cells. A single recessive nuclear gene was found to control male sterility in the mutant, while MutMap and KASP analyses identified BraA05g022470.3C (BrGGL7), which encodes a GDSL esterase / lipase, as the candidate mutant gene. A single nucleotide substitution from C to T occurred within the domain of BrGGL7 in ftms1, resulting in premature translation termination in the fourth exon. Meanwhile, qRT-PCR analysis indicated that BrGGL7 was prominently expressed in the anthers, and expression was greater in the wild-type 'FT' than ftms1. Genetic complementation of the orthologous Arabidopsis ggl7 mutant further confirmed the role of BrGGL7 in pollen development. These findings suggest that BrGGL7 plays a fundamental role in pollen formation, providing important insight into the molecular mechanisms underlying male sterility in Chinese cabbage.

Defect in BrMS1, a PHD-finger transcription factor, induces male sterility in ethyl methane sulfonate-mutagenized Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Male sterility is an ideal character for the female parent in commercial hybrid seed production in Chinese cabbages. We identified three allele male sterile mutants msm2-1/2/3 in progenies of ethyl methane sulfonate mutagenized Chinese cabbage. It was proved that their male sterilities were controlled by a same recessive nuclear gene. Cytological observation showed that the delayed tapetal programmed cell death (PCD) as well as the abnormal pollen exine and intine led to pollen abortion in these mutants. MutMap combined with KASP analyses showed that BraA10g019050.3C, a homologous gene of AtMS1 encoding a PHD-finger transcription factor and regulated pollen development, was the causal gene. A single-nucleotide mutation from G to A occurred at the 2443th base of BrMS1 in msm2-1 which results in premature termination of the PHD-finger protein translation; a single-nucleotide mutation from G to A existed at 1372th base in msm2-2 that makes for frameshift mutation; a single-nucleotide mutation from G to A distributed at 1887th base in msm2-3 which issues in the amino acid changed from Asp to Asn. The three allelic mutations in BrMS1 all led to the male sterile phenotype, which revealed its function in stamen development. Quantitative reverse transcription polymerase chain reaction analysis indicated that BrMS1 specially expressed in the anther at the early stage of pollen development and its expression level was higher in msm2-1/2/3 than that in the wild-type "FT." BrMS1 was located at the nucleus and a length of 12 amino acid residues at the C-terminus had transcriptional activation activity. RNA-seq indicated that the mutation in BrMS1 affected the transcript level of genes related to the tapetum PCD and pollen wall formation, which brought out the pollen abortion. These male sterile mutants we developed provided a novel gene resource for hybrid breeding in Chinese cabbage.

Identification of an anther-specific promoter from a male sterile AB line in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

The promoter of the male sterile gene is important for studying male sterility. In this study, BraA08g014780.3C which differentially expressed between male sterile and fertile plants was identified from a genetic male sterile AB line of Chinese cabbage by RNA-seq. qRT-PCR revealed that BraA08g014780.3C was mainly expressed in the early stage of floral bud development in fertile plants, and preferentially expressed in their anthers. The promoter of BraA08g014780.3C was cloned and analyzed. Cis acting element analysis showed that the promoter of BraA08g014780.3C contains POLLEN1LELAT52 and GTGANTG10, which are both pollen-specific expression elements. The BraA08g014780.3Cp::GUS vector was constructed, then transformed to Arabidopsis thaliana Col-0. PCR analysis and sequencing of the transgenic Arabidopsis revealed that the GUS gene driven by the BraA08g014780.3C promoter was successfully transformed to the Arabidopsis. GUS staining indicated that the promoter of BraA08g014780.3C was an anther-specific promoter. Identifying the anther-specific promoter laid a foundation for revealing BraA08g014780.3C function in male sterility of Chinese cabbage.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03160-z.

MicroRNAs Involved in Regulatory Cytoplasmic Male Sterility by Analysis RNA-seq and Small RNA-seq in Soybean

Cytoplasmic male sterility (CMS) is an important plant characteristic for exploiting heterosis to enhance crop traits during breeding. However, the CMS regulatory network remains unclear in plants, even though researchers have attempted to isolate genes associated with CMS. In this study, we performed high-throughput sequencing and degradome analyses to identify microRNAs (miRNAs) and their targets in a soybean CMS line (JLCMS9A) and its maintainer line (JLCMS9B). Additionally, the differentially expressed genes during reproductive development were identified using RNA-seq data. A total of 280 miRNAs matched soybean miRNA sequences in miRBase, including mature miRNAs and pre-miRNAs. Of the 280 miRNAs, 30, 23, and 21 belonged to the miR166, miR156, and miR171 families, respectively. Moreover, 410 novel low-abundant miRNAs were identified in the JLCMS9A and JLCMS9B flower buds. Furthermore, 303 and 462 target genes unique to JLCMS9A and JLCMS9B, respectively, as well as 782 common targets were predicted based on the degradome analysis. Target genes differentially expressed between the CMS line and the maintainer line were revealed by an RNA-seq analysis. Moreover, all target genes were annotated with diverse functions related to biological processes, cellular components, and molecular functions, including transcriptional regulation, the nucleus, meristem maintenance, meristem initiation, cell differentiation, auxin-activated signaling, plant ovule development, and anther development. Finally, a network was built based on the interactions. Analyses of the miRNA, degradome, and transcriptome datasets generated in this study provided a comprehensive overview of the reproductive development of a CMS soybean line. The data presented herein represent useful information for soybean hybrid breeding. Furthermore, the study results indicate that miRNAs might contribute to the soybean CMS regulatory network by modulating the expression of CMS-related genes. These findings lay the foundation for future studies on the molecular mechanisms underlying soybean CMS.

Genome-wide analysis of proline-rich extension-like receptor protein kinase (PERK) in Brassica rapa and its association with the pollen development

BACKGROUND

Proline-rich extension-like receptor protein kinases (PERKs) are an important class of receptor kinases located in the plasma membrane, most of which play a vital role in pollen development.

RESULTS

Our study identified 25 putative PERK genes from the whole Brassica rapa genome (AA). Phylogenetic analysis of PERK protein sequences from 16 Brassicaceae species divided them into four subfamilies. The biophysical properties of the BrPERKs were investigated. Gene duplication and synteny analyses and the calculation of Ka/Ks values suggested that all 80 orthologous/paralogous gene pairs between B. rapa and A. thaliana, B. nigra and B. oleracea have experienced strong purifying selection. RNA-Seq data and qRT-PCR analyses showed that several BrPERK genes were expressed in different tissues, while some BrPERKs exhibited high expression levels only in buds. Furthermore, comparative transcriptome analyses from six male-sterile lines of B. rapa indicated that 7 BrPERK genes were downregulated in all six male-sterile lines. Meanwhile, the interaction networks of the BrPERK genes were constructed and 13 PERK coexpressed genes were identified, most of which were downregulated in the male sterile buds.

CONCLUSION

Combined with interaction networks, coexpression and qRT-PCR analyses, these results demonstrated that two BrPERK genes, Bra001723.1 and Bra037558.1 (the orthologs of AtPERK6 (AT3G18810)), were downregulated beginning in the meiosis II period of male sterile lines and involved in anther development. Overall, this comprehensive analysis of some BrPERK genes elucidated their roles in male sterility.

Single-molecule real-time sequencing facilitates the analysis of transcripts and splice isoforms of anthers in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

Anther development has been extensively studied at the transcriptional level, but a systematic analysis of full-length transcripts on a genome-wide scale has not yet been published. Here, the Pacific Biosciences (PacBio) Sequel platform and next-generation sequencing (NGS) technology were combined to generate full-length sequences and completed structures of transcripts in anthers of Chinese cabbage.

RESULTS

Using single-molecule real-time sequencing (SMRT), a total of 1,098,119 circular consensus sequences (CCSs) were generated with a mean length of 2664 bp. More than 75% of the CCSs were considered full-length non-chimeric (FLNC) reads. After error correction, 725,731 high-quality FLNC reads were estimated to carry 51,501 isoforms from 19,503 loci, consisting of 38,992 novel isoforms from known genes and 3691 novel isoforms from novel genes. Of the novel isoforms, we identified 407 long non-coding RNAs (lncRNAs) and 37,549 open reading frames (ORFs). Furthermore, a total of 453,270 alternative splicing (AS) events were identified and the majority of AS models in anther were determined to be approximate exon skipping (XSKIP) events. Of the key genes regulated during anther development, AS events were mainly identified in the genes SERK1, CALS5, NEF1, and CESA1/3. Additionally, we identified 104 fusion transcripts and 5806 genes that had alternative polyadenylation (APA).

CONCLUSIONS

Our work demonstrated the transcriptome diversity and complexity of anther development in Chinese cabbage. The findings provide a basis for further genome annotation and transcriptome research in Chinese cabbage.

Investigation of the genes associated with a male sterility mutant (msm) in Chinese cabbage (Brassica campestris ssp. pekinensis) using RNA-Seq

In Chinese cabbage, hybrid seed production is performed using male sterility lines, an important approach to heterosis utilization. In this study, a stably inherited male sterile mutant msm was obtained from the 'FT'-doubled haploid line of Chinese cabbage using isolated microspore culture combined with ⁶⁰Co γ-ray mutagenesis. The genetic backgrounds of 'FT' and msm were highly consistent; however, compared with wild-type 'FT', msm exhibited completely degenerated stamens and no pollen phenotype. Other characters showed no significant differences. Cytological observations revealed that stamen abortion in msm begins during the tetrad period and that tapetum cells were abnormally expanded and highly vacuolated, leading to microspore abortion. Genetic analysis indicated that the msm mutant phenotype is controlled by a single recessive nuclear gene. Comparative transcriptome analysis of 'FT' and msm flower buds using RNA-Seq technology revealed 1653 differentially expressed genes, among which, a large number associated with male sterility were detected, including 64 pollen development- and pollen tube growth-related genes, 94 pollen wall development-related genes, 11 phytohormone-related genes, and 16 transcription factor-related genes. An overwhelming majority of these genes were down-regulated in msm compared with 'FT'. Furthermore, KEGG pathway analysis indicated that a variety of carbohydrate metabolic and lipid metabolic pathways were significantly enriched, which may be related to pollen abortion. The expression patterns of 24 male sterility-related genes were analyzed using qRT-PCR. In addition, 24,476 single-nucleotide polymorphisms and 413,073 insertion-deletion events were specifically detected in msm. These results will facilitate elucidation of the regulatory mechanisms underlying male sterility in Chinese cabbage.

Transcriptome-Wide Identification and Characterization of Circular RNAs in Leaves of Chinese Cabbage (Brassica rapa L. ssp. pekinensis) in Response to Calcium Deficiency-Induced Tip-burn

Circular RNA (circRNA) is a newly discovered non-coding RNA, which play significant roles in the function and transcriptional regulation of microRNA. To date, in Chinese cabbage, the functional characteristic of circRNAs in response to calcium deficiency-induced tip-burn have not been reported. In this study, 730 circRNAs were isolated from Chinese cabbage leaves, of which 23 and 22 were differentially expressed in different calcium deficiency stages compared with the control. Forty-six host genes of the differentially expressed circRNAs were identified, and one circRNA was found to act as miRNAs sponges. Based on the functional analysis of host genes and target mRNAs of the corresponding miRNAs, the identified circRNAs might participated in response to stimulus, electron carrier activity, ATPase activity, cell wall metabolism, transcription factors and plant hormone signal transduction. ABF2, a positive regulator of the abiotic stress response in the abscisic acid (ABA) pathway, may play a role in calcium deficiency tolerance through a circRNA regulatory pathway. Correspondingly, the concentration of ABA is also increased during the Ca²⁺ deficiency stress. Our results suggest that circRNAs participate in a broad range of biological processes and physiological functions in the response to calcium deficiency-induced tip-burn and provide a basis for further studies of the biological roles that circRNAs play in the plant stress response.

Transcriptome-based analysis of carotenoid accumulation-related gene expression in petals of Chinese cabbage (Brassica rapa L.)

To identify genes associated with carotenoid accumulation in petals of Chinese cabbage, the composition and content of carotenoids were analyzed, and comparative transcriptome sequencing was performed between the yellow flower line, 92S105, and the orange flower line, 94C9. High-performance liquid chromatography (HPLC) revealed that petals of 92S105 were high in violaxanthin as well as lutein, whereas petals of 94C9 showed considerable levels of lutein and β-carotene. Transcriptome analysis showed that 3534 and 3833 genes were up- and down-regulated in 94C9, respectively. Among these differentially expressed genes (DEGs), many related to carotenoid accumulation were identified, including 12 carotenoid biosynthesis pathway genes, 4 transcription factor genes, and 1028 specifically expressed genes. β-carotene hydroxylase 1 (BrBCH1), BrBCH2, zeaxanthin epoxidase (BrZEP), and MYB transcription factor gene (BrGAMYB) were down-regulated in petals of 94C9 when compared with petals of 92S105, which caused β-carotene accumulation and may lead to orange petal color in 94C9. Expression levels of 20 DEGs were verified by qPCR and the results were highly consistent with those of transcriptome sequencing. Moreover, Gene Ontology (GO) enrichment analysis revealed that membrane, binding, and metabolic processes were the most significantly enriched GO terms in cellular component, molecular function, and biological process ontologies, respectively. In conclusion, our study analyzed the differences in composition and content of carotenoids between 92S105 and 94C9 and identified potential candidate genes related to carotenoid accumulation in petals, thereby creating a solid foundation for future studies on the mechanism regulating carotenoid accumulation in petals of Chinese cabbage.

Analysis of metabolic pathways related to fertility restoration and identification of fertility candidate genes associated with Aegilops kotschyi cytoplasm in wheat (Triticum aestivum L.)

BACKGROUND

Thermo-sensitive male-sterility based on Aegilops kotschyi cytoplasm (K-TCMS) plays an important role in hybrid wheat breeding. This has important possible applications in two-line hybrid wheat breeding but the genetic basis and molecular regulation mechanism related to fertility restoration are poorly understood. In this study, comparative transcriptome profiling based on RNA sequencing was conducted for two near-isogenic lines comprising KTM3315R and its sterile counterpart KTM3315A, a total of six samples (3 repetitions per group), in order to identify fertility restoration genes and their metabolic pathways.

RESULTS

In total, 2642 significant differentially expressed genes (DEGs) were detected, among which 1238 were down-regulated and 1404 were up-regulated in fertile anthers. Functional annotation enrichment analysis identified important pathways related to fertility restoration, such as carbohydrate metabolism, phenylpropanoid metabolism and biosynthesis, as well as candidate genes encoding pectin methylesterase and flavanone 3-hydroxylase. Moreover, transcription factor analysis showed that a large number of DEGs were mainly involved with the WRKY, bHLH, and MYB transcription factor families. Determination of total soluble sugar and flavonoid contents demonstrated that important metabolic pathways and candidate genes are associated with fertility restoration. Twelve DEGs were selected and detected by quantitative reverse-transcribed PCR, and the results indicated that the transcriptome sequencing results were reliable.

CONCLUSIONS

Our results indicate that identified DEGs were related to the fertility restoration and they proved to be crucial in Aegilops kotschyi cytoplasm. These findings also provide a basis for exploring the molecular regulation mechanism associated with wheat fertility restoration as well as screening and cloning related genes.

Profiling of indole metabolic pathway in thermo-sensitive Bainong male sterile line in wheat (Triticum aestivum L.)

Bainong male sterile (BNS) wheat (Triticum aestivum L.) is a thermo-sensitive genic male sterile line with excellent sterility and self-restoration. We focused on transcriptional profiles of differentially expressed probes between BNS sterile and fertile anthers. Anthers, rachis and spikes from sterile line and fertile line were collected. Extracted RNA was assayed using wheat expression microarray and Gene Ontology was analyzed using Cytoscape with ClueGO. An indole (indole-3-acetic acid: IAA) metabolism pathway sub-network was almost formed in all differentially expressed profiles between sterile and fertile samples. IAA sub-network contained four nodes of indole and alkaloid metabolism connecting main network via indole compounds. This sub-network was absent in rachis and intact in transformed fertile anthers, which was the main differently expressed metabolism pathway in F₁ anthers with restorer genes. Alkaloid metabolism was absent in sterile anthers. Abnormal metabolism of IAA may be involved in BNS sterility. BNS transformation may be regulated by the production of IAA and alkaloid metabolism pathway, which favor the safe utilization of the sterile line in hybrid wheat production.

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.