Brems1 mutation induced tapetum deficiency leading to male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Brems1 mutation induced tapetum deficiency leading to male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis) | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Brems1 mutation induced tapetum deficiency leading to male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis) Chuanhong Liu, Lin Wang, Chong Tan, Di Zhao, Zhiyong Liu This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4489236/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 24 Feb, 2025 Read the published version in Theoretical and Applied Genetics → Version 1 posted 5 You are reading this latest preprint version Abstract Male sterile lines are ideal for hybrid seed production in Chinese cabbage. Herein, the complete male sterile mutants M5026 and M5073 were obtained through ethyl methanesulfonate (EMS) mutagenesis in the cabbage double haploid line ‘FT’. Observation of paraffin sections showed that M5026 lacked the tapetum and had excessive microsporocytes. Transmission electron microscopy (TEM) revealed abnormal exine formation in M5026 microspores. Genetic analysis revealed a single recessive nuclear gene caused the male sterility phenotype of M5026. Using Mutmap sequencing and Kompetitive allele-specific PCR (KASP) identification and gene cloning, BraA10g029920.3.5C, encoding EMS1 (Excess microsporocytes 1), a leucine-rich repeat receptor-like protein kinase (LRR-RLK), was identified as the candidate gene of M5026 and named Brems1. A nonsynonymous G-to-A mutation in an exon of the Brems1 gene in M5026 resulted in the substitution of glycine with arginine. Employing Mutmap and cloning approaches, a C-to-T SNP was identified within Brems1 of M5073, resulting in premature translation termination. Both BrEMS1 and Brems1 were subcellularly localized at the cell membrane. qRT-PCR analysis indicated Brems1 exhibited the highest expression level in flower buds, while no expression was detected in roots. Transcriptomic analysis revealed that mutation in Brems1 reduced the expression levels of genes associated with the tapetum, pollen tube, and LRR-RLK family. These results suggested that Brems1 plays a critical role in pollen development and contributes to elucidating the molecular mechanisms underlying tapetum development and male sterility in Chinese cabbage. Chinese cabbage male sterility Brems1 transcriptome Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Key Message The mutation in Brems1 resulting in male sterility in Chinese cabbage were validated through two allelic mutations. Introduction Male sterility in plants refers to the phenomenon in which the pistil develops normally while the stamen degenerates and fails to produce viable pollen. This is an effective method for utilizing hybrid vigor and hybrid breeding. Based on the inheritance pattern of sterility genes, male sterility can be classified into cytoplasmic male sterility and nuclear male sterility (Chen and Liu 2014 ). Cytoplasmic male sterility was typically caused by an interaction between mitochondrial or chloroplast genes and nuclear genes, and fertility could be restored by nuclear-encoded restorer genes (Yamagishi and Bhat 2014 ). Nuclear male sterility was controlled via the nuclear genome, where the molecular regulatory mechanisms of plant male sterility involved genes primarily participating in microspore development, tapetum formation, exine deposition, intine formation, and pollen wall rupture (Wilson et al. 2011 ). The tapetum, located on the innermost layer of the anther wall, directly adjacent to the microspore mother cells, played a crucial role in the nutrient supply and regulating microsporogenesis. After cellular differentiation, the tapetum transformed into polarized secretory cells, synthesizing and secreting callose enzymes via programmed cell death to release uninucleate microspores, while also promoting pollen wall formation and serving as a significant source of nutrition for microspore development (Wu et al. 2000). Plants with mutations that lacked the tapetum or experienced genetic ablation of the tapetum were unable to produce pollen grains (Mariani et al. 1990 ). The Arabidopsis thaliana male-sterile mutants EMS1 (Zhao et al. 2002 ) and TPD1 (Yang et al. 2003 ), which lacked the tapetum and resulted in male infertility. Premature or delayed degeneration tapetum could also lead to male sterility. Dysfunction of the tapetum in the Arabidopsis thaliana male-sterile mutant DYT1 resulted in abnormal pollen development (Zhang et al. 2006 ). Several mutants with abnormal tapetum development have been identified in Chinese cabbage, including BrMS1 (Dong et al. 2022 ), BrACOS5 (Zou et al. 2023 ), BrGGL7 (Zhao et al. 2022 ), and BrABCG26 (Xu et al. 2024 ). The leucine-rich repeat receptor-like kinase (LRR-RLK) consisting of an intracellular domain with kinase activity, a transmembrane domain and extracellular domains responsible for ligand recognition (Shiu et al. 2004 ). Some of the LRR-RLK members family participated in plant stamen development. In Arabidopsis , the receptor-like kinase ERECTA mutant exhibited abnormal flower development in the form of absent anthers (Shpak et al. 2004 ); the triple mutant er-105 , erl1-2 , erl2-1 , encoding an LRR-RLK, resulted in undifferentiated anthers, leading to male sterility (Hord et al. 2008 ); EMS1 gene encoded an LRR-RLK, and its mutant exhibited male sterility due to the absence of middle layer and tapetum cells (Zhao et al. 2002 ). In rice, a mutation in the LRR-RLK gene MSP1 resulted in absent of the middle and tapetum cell layers, leading to male sterility (Yang et al. 2016 ); In addition, a mutation in the ERECTA homologous gene, OsERL , resulted in significant defects in anther development (Liu et al. 2021b ); and knocking down the expression level of the LRR-RLK family gene OsPEX1 using RNAi resulted in male sterility (Hang et al. 2021 ). In cotton, plants with RNAi of GhSERK1 exhibited reduced pollen viability (Shi et al. 2014 ), and editing the GhEMS1 gene using CRISPR/Cas9 produced, male-sterile mutants (Zhang et al. 2022 ). Currently, no reports detailing the involvement of members of the LRR-RLK family in the molecular mechanisms of male fertility in Chinese cabbage have been made. Here, two male sterile mutants, M5026 and M5073 , were obtained using ethyl methanesulfonate (EMS)-induced mutagenesis of ‘FT’. Cytological observations revealed that M5026 lacked tapetum, exhibited excessive microsporocytes production, and showed abnormal microspore exine morphology. Through MutMap, Kompetitive allele-specific PCR (KASP), and cloning sequencing, the predicted candidate gene was identified as BraA10g029920.3.5C ( Brems1 ), encoding an LRR-RLK. Employing Mutmap and cloning approaches, a C-to-T SNP was identified within Brems1 of M5073 , resulting in premature translation termination. BrEMS1 and Brems1 were subcellularly localized at the cell membrane. Brems1 exhibited high expression levels in buds. Transcriptomic analysis indicated that the mutation in Brems1 reduced the expression of genes associated with trichomes, pollen tubes and the LRR-RLK family. This study cloned the Brems1 gene in Chinese cabbage for the first time, contributing to elucidating the molecular mechanisms underlying tapetal development and male sterility in Chinese cabbage. Materials and methods Plant materials The DH line ‘FT’ of Chinese cabbage obtained via microspore culturing was used as the wild-type (Huang et al. 2016 ). The male-sterile mutants M5026 and M5073 were generated using EMS-induced mutagenesis of ‘FT’. Two parents were used to generate the F 1 and F 2 generations for genetic analysis. The F 2 generation was further used for gene mapping and KASP analysis. Floral organ observation Floral organ morphology (pistils, long stamens, short stamens, petals, sepals, and flower buds) of ‘FT’ and M5026 were observed at the reproductive growth stage using a stereomicroscope (Nikon SMZ800, Japan). Pollen viability analysis The pollen viabilities of ‘FT’ and M5026 were assessed using the 2,3,5-triphenyltetrazolium chloride (TTC) staining method. Anthers were collected from flower buds and placed on a glass slide. 50 µL of TTC (0.2%) staining solution was added to the slide, and a coverslip was placed on top. The slides were incubated at 37°C for 20–30 min. Pollen viability was observed using a fluorescence inverted microscope (OLYMPUS DP80, Japan). Paraffin section observation The buds at different developmental stages of ‘FT’ and M5026 were individually fixed in Formalin-Aceto-Alcohol (FAA) fixative for 12 hours, followed by dehydration using an ethanol solution gradient of 50%-100%. Subsequently, the samples were subjected to substitution with a mixed solution of ethanol and xylene (volume ratios = 3:1, 1:1, 1:3), and different mass ratios of xylene and paraffin were then impregnated into the samples. The samples were embedded in paraffin, and sections were cut using a microtome (Leica RM2016, Germany). After staining the anther sections with safranin and fast green, these sections were observed using an optical microscope (ECLIPSE80i, Nikon). TEM observation Buds of ‘FT’ and M5026 were fixed via immersion in a 1:1 mixture of glutaraldehyde and paraformaldehyde. The fixed samples were then stored in the dark at 4°C. The detailed steps for transmission electron microscopy (TEM) were conducted according to the previously described protocols (Ji et al. 2017 ). Genetic analysis Constructing F 1 and F 2 generations were constructed using parental strains M5026 and ‘FT’ for genetic analysis. The number of male fertile and sterile plants in the F 1 and F 2 generations was observed and statistically analyzed. The chi-square (χ 2 ) test was employed to analyze the segregation ratio in the F 2 population. Identification of the candidate gene Candidate genes in M5026 and M5073 were identified using an improved Mutmap method. (Abe et al. 2012 ). A pool of male sterile mutants was constructed by collecting genomic DNA from 58 F 2 plants exhibiting male sterility. Re-sequencing was performed using NovaSeq 6000 platform (Illumina, San Diego, California, USA) following DNA extraction from ‘FT’ and male sterile mutant pool using a DNAsecure Plant Kit (Tiangen Biotech Ltd, Beijing, China). The sliding window method implemented in fastp (v0.20.0) was utilized for filtering raw data to obtain high-quality data (Liu et al. 2021a ). Subsequently, these filtered high-quality data were aligned to the reference genome using bwa (0.7.12-r1039) (Li and Durbin 2010 ). Single-nucleotide polymorphisms (SNPs) and insertion-deletions (INDELs) were detected using the Genome Analysis Toolkit (GATK) (McKenna et al. 2010 ). ANNOVAR was employed for functional annotation (Wang et al. 2010 ). Circos was used to visualize the mapping of mutation information onto the genome (Krzywinski et al. 2009 ). SNP genotype by KASP A total of 70 F 2 male sterile individuals and two ‘FT’ plants were used as materials to conduct a KASP assay to detect SNP co-segregation and identify the Brems1 candidate gene. The thermal cycling conditions for KASP were set as described previously (Xi et al. 2018 ). The primers used in the KASP assay were listed in Table S1 . Cloning and sequencing of Brems1 The full-length DNA and coding sequence (CDS) of the Brems1 were amplified from ‘FT’, M5026 and M5073 . The specific operational steps were conducted according to a previous study (Liu et al. 2021a ). Sanger sequencing was performed by Sangon Biotech (Shanghai, China). Bioinformatic characterization and phylogenetic analysis of the BrEMS1 protein The full-length amino acid sequence of BrEMS1 was submitted to the SMART ( https://smart . embl- heidelberg. de/) website for protein domain analysis. The three-dimensional protein structures of BrEMS1 and Brems1 were analyzed using SwissModel. A total of 1000 bootstrap replicates were established in MEGA6.0 software, and the neighbor-joining method was employed to construct the phylogenetic tree. DNAMAN v6.0 was used to perform multiple sequence alignment of the protein sequences. Subcellular localization of BrEMS1 and Brems1 The CDS sequences of BrEMS1 and Brems1 were cloned from ‘FT’ and M5026 , respectively, and inserted into the ProCAMV35S vector and transformed into Agrobacterium tumefaciens GV3101, and impregnated into tobacco leaf mesophyll cells at an optical density (OD) of 0.6–0.8. After 24 hours of dark treatment followed by 24 hours of light treatment, fluorescence signals were observed using a confocal laser scanning microscope (Leica Microsystems, Wetzlar, Germany). Green fluorescence protein (GFP) fluorescence signals were detected at a wavelength of 496 ~ 540 nm. Quantitative Real-Time PCR (qRT-PCR) The total RNA isolated from flowers, leaves, roots, buds, and stems of ‘FT’ and M5026 plants was subsequently reverse-transcribed into cDNA. qRT-PCR was performed using Ultra SYBR Green Mix (Kangwei Century, Beijing, China) on a QuantStudio 6 PCR system. Actin was selected as the internal control. The relative expression levels were calculated using the 2 −ΔΔCt method. Transcriptome analysis Buds from ‘FT’ and M5026 were individually sampled in triplicate, followed by RNA extraction. Subsequently, the RNA samples were subjected to rigorous quality control using an Agilent 2100 bioanalyzer, and libraries were constructed accordingly. The concentration of the libraries was first preliminarily quantified using a Qubit 2.0 Fluorometer, and then the insert size of the libraries was detected using an Agilent 2100 bioanalyzer. The effective concentrations of the libraries were accurately quantified via qRT-PCR to ensure the quality of the libraries. Once the libraries passed quality control, different libraries were pooled based on their effective concentrations and the desired amount of data for sequencing on the Illumina platform. These raw sequencing data were then filtered as follows: reads with adapters were removed; reads containing N (representing unknown bases) were removed; and low-quality reads (more than 50% of the bases with Qphred scores less than or equal to 5) were removed (Yan et al. 2013 ). These sequencing data were aligned to the reference genome Brara_Chiifu_V3.5 using HISAT2 software to obtain a rapid and precise mapping of the reads to the reference genome (Mortazavi et al. 2008 ). Based on the positional information of gene alignment on the reference genome, reads covering each gene within the start and end range were counted. Reads with alignment quality scores lower than 10, unpaired alignments, and reads aligned to multiple regions of the genome were filtered out. This analysis was conducted using the featureCounts tool within the subread software (Liao et al. 2014 ). After gene expression quantification was completed, a statistical analysis was conducted to identify genes that exhibited significant differences in expression levels across different conditions. The differential gene set was subjected to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis using the clusterProfiler software. Results The M5026 mutant displayed male sterility The male sterility M5026 mutant was isolated from the EMS-induced mutant DH line ‘FT’. Pollen grains were absent on the anther surface of the M5026 mutant (Fig. 1a, b). Pollen viability analysis revealed that ‘FT’ produced viable pollen grains. In contrast, M5026 did not produce any pollen grains (Fig. 1c). Observations of the floral organs of M5026 and ‘FT’ plants revealed that the long stamens, short stamens, pistils, petals, and buds of M5026 mutant were all shorter than those of ‘FT’ (Fig. 1d-h), while the sepals did not different (Fig. 1i). The F 1 and F 2 generations were constructed using the parental strains M5026 and ‘FT’ to conduct genetic analysis. All F 1 plants showed the male fertile phenotype. The F 2 population had 581 male fertile and 187 male sterile plants, a ratio of 3:1 (χ 2 = 0.17) (Table 1 ). A single recessive nuclear gene, Brems1 , regulated the male sterility phenotype in M5026 . Table 1 Genetic analysis in M5026. Generations Total Male fertility Male sterile Segregation ratio χ 2 χ 2 0.05 P 1 (‘FT’) 50 50 0 P 2 ( M5026 ) 50 0 50 F 1 (P 2 × P 1 ) 20 20 0 F 1 ’(P 1 × P 2 ) 20 20 0 F 2 768 581 187 3.10:1 0.17 3.841 Table 2 Candidate SNP information. Chromosome Position WT MUT Snpindex Region MutationType GeneID A10 19,022,513 C T 1 intronic / BraA10g029920.3.5C A10 19,087,850 C T 1 exonic nonsynonymous SNV BraA10g030050.3.5C A10 19,417,177 C T 1 exonic synonymous SNV BraA10g030900.3.5C Tapetum deficiency, excessive microsporocytes, and abnormal exine in the Mutant M5026 To further elucidate the cytological characteristics of the sterile mutant, we employed paraffin sectioning to examine the anthers of ‘FT’ and M5026 at different developmental stages. Based on the established categorization of Arabidopsis anther developmental stages (Sanders et al. 1999 ), we selected eight anther developmental stages of anther development in both ‘FT’ and M5026 for analysis. During stage 4 (S4) of anther development, no differences were observed between ‘FT’ (Fig. 2a) and M5026 (Fig. 2e). At stage 5 (S5), ‘FT’ developed the tapetum characterized by rectangular-shaped cells forming a layer outside the microsporocytes (Fig. 2b), however, M5026 lacked the tapetum and exhibited excess microsporocytes (Fig. 2f). At stage 6 (S6), the microsporocytes of M5026 exhibited morphological abnormalities characterized by concavity (Fig. 2g). During stages 7–12 (S7-12), the microsporocytes of ‘FT’ sequentially formed tetrads, microspores, and pollen grains (Fig. 2d, i-l). In contrast, M5026 continued to lack the tapetum, and its microsporocytes underwent continuous degradation, ultimately failing to produce pollen grains (Fig. 2h, m-p). TEM was utilized to observe further ‘FT’ and M5026 . The microspore exine structure of ‘FT’ was intact, with well-defined bacula and tectum morphology (Fig. 3a, b), while the microspore of M5026 exhibited an abnormal exine structure (Fig. 3c, d). Identification of candidate gene contributing to the M5026 male sterility phenotype. The adapted MutMap approach was employed to isolate candidate genes. Sequencing data filtering was performed to obtain high-quality data. The wild-type and F 2 mutant pool yielded 175,935,758 (95.25%) and 177,453,976 (96.23%) high quality dates, respectively. Of these, 98.98% and 99.21% could be mapped to the Brara_Chiifu_V3.5 reference genome. SNPs were detected and filtered using the GATK software, resulting in a total of 732,795 SNPs. Subsequently, SNP-index calculation was conducted, followed by plotting the distribution of the SNP-index across the chromosome (Fig. 4a). In the case of selecting regions with SNP-index values exceeding the 95th percentile as candidate regions associated with the M5026 male sterility trait, the candidate gene, Brems1 , was mapped to the region of chromosome A10 from 18,400,000 to 20,725,698. A total of 30 SNPs were found in the candidate interval, among which three have SNP-index values of 1 and were located within genes, namely SNP19,022,513, SNP19,087,850, and SNP19,417,177. SNP19,022,513 was located in an intron, SNP19,087,850 was a nonsynonymous mutation located in an exon, and SNP19,417,177 was a synonymous mutation located in an exon. To validate the causal SNP for male sterility, KASP genotyping was performed on three SNPs using two ‘FT’ and seventy F 2 male sterile individuals. The results showed that only SNP 19,022,513 co-segregated with the male sterility phenotype, with all ‘FT’ plants being C:C and all mutant phenotype plants being T:T. Recombination was identified at SNP 19,087,850, with one C:C and one C:T genotype observed among the mutant phenotype plants. Similarly, at SNP 19,022,513, recombination was also detected, revealing one C:C and two C:T genotypes in the plants exhibiting the mutant phenotype (Fig. 4b). SNP 19,022,513 was located on BraA10g029920.3.5C. Its Arabidopsis ortholog is EMS1 ( AT5G07280 ), which encodes EXCESS MICROSPOROCYTES1 (EMS1), an LRR-RLK that controls the fate of somatic and germ cells in Arabidopsis anthers. Cloning and sequence analysis of Brems1 The DNA and cDNA sequences of BraA10g029920.3.5C were cloned from M5026 and ‘FT’, respectively. The results revealed that BraA10g029920.3.5C has a full length of 3699 bp, comprising of a single exon. A non-synonymous mutation was identified at position 3439 bp, where a G to A transition occurred, changing GGG to AGG, resulting in the substitution of glycine (Gly) with arginine (Arg) (Fig. 5a). The gene sequence of BraA10g029920.3.5C cloned from our experimental materials exhibited differences compared with the reference genome Brara_Chiifu_V3.5. In the reference genome, the full length of BraA10g029920.3.5C was 4613 bp (Fig. 5b1), which was 911 bp longer than in ‘FT’ and M5026 , with the additional length located at the gene’s 5’ end of the gene (Fig. 5b2). BraA10g029920.3.5C in the reference genome contained three exons and two introns (Fig. 5b1), while those in ‘FT’ and M5026 had only one exon and no introns (Fig. 5b2). The mutation site in M5026 was located in the second intron of BraA10g029920.3.5C in the reference genome and in the exon of ‘FT’ and M5026 (Fig. 5b). The EMS1 mutation in Arabidopsis thaliana resulted in male sterility (Zhao et al. 2002 ). Hence, we presumed BraA10g029920.3.5C was the candidate gene for Brems1 . Identification of the allelic material M5073 A male sterile mutant M507 3 was induced by EMS mutagenesis in ‘FT’ (Fig. S1 a). M5026 and M5073 were crossed with ‘FT’ to generate F 1 hybrids, respectively. Allelism tests were performed by crossing these two F 1 generations. The observed ratio of male fertility to sterility in the progeny was 3:1 (χ 2 = 1.405) (Table S2 ), indicating that the male sterility traits in M5026 and M5073 were controlled by a single allelic gene. Subsequently, Mutmap method was employed to isolate candidate genes of M5073 . Regions with SNP index values exceeding the 95th percentile were selected as candidate intervals, a SNP with an SNP-index value of 1 was identified at position 19,024,244 bp on chromosome A10, located within the Brems1 gene (Fig. S1 b). Cloning of the Brems1 gene from M5073 revealed a C-to-T mutation at position 1708 bp, converting CAG to TAG, resulting in the substitution of glutamine for a stop codon and premature translation termination (Fig. S1 c). These results indicated that the mutation in Brems1 caused the male sterile phenotype observed in M5026 and M5073 , confirming its role in stamen development. Phylogenetic and structural analysis of BrEMS1 To further investigate the role of BrEMS1 in male fertility of Chinese cabbage, the structure of the BrEMS1 protein was analyzed. The BrEMS1 protein comprised 1233 amino acids, with a theoretical isoelectric point of 6.62 and a molecular weight of 133.2 kDa. It harbored a serine or threonine-specific kinase subfamily (S_TKc) along with 17 LRR structural domains (Fig. 6a). The mutation in M5026 occurred at the 1147th amino acid of the BrEMS1 protein, where Gly was substituted with Arg, located within the S_TKc domain. The SwissModel tool was employed to predict the protein structures of both BrEMS1 and Brems1, revealing that the mutation in BrEMS1 altered the three-dimensional structure of the protein (Fig. 6b). The evolutionary relationships among BrEMS1 homologous sequences across different species were elucidated through phylogenetic analysis, demonstrating its high homology among multiple crops in the Brassicaceae family (Fig. 6c). The S_TKc domain was highly conserved across multiple crops in the Brassicaceae family (Fig. 6d). Expression pattern analysis of Brems1 To elucidate the expression pattern of Brems1 , we conducted a subcellular localization and fluorescence quantitative analysis of Brems1 . BrEMS1-GFP and Brems1-GFP were introduced into tobacco and their co-localization with the cell membrane marker (Pip2a-RFP) was analyzed. The results demonstrated that the GFP fluorescence of both BrEMS1-GFP and Brems1-GFP was localized on the cell membrane (Fig. 7a), indicating that the mutation of BrEMS1 did not affect on its subcellular localization. The relative transcription levels of Brems1 in different organs of ‘FT’ and M5026 were assessed using qRT-PCR. qRT-PCR analysis revealed no significant differences in t he relative abundances of Brems1 transcripts among the flowers, leaves, roots, buds, and stems of ‘FT’ and M5026 plants. Notably, Brems1 expression was highest in buds, followed by flowers, leaves, and stems, while expression was absent in roots (Fig. 7b). Brems1 affects the transcription levels of genes associated with the tapetum and pollen tube To elucidate the molecular mechanisms underlying Brems1 regulation of male sterility, transcriptomic analysis was conducted on the buds of ‘FT’ and M5026 . ‘FT’ and M5026 obtained totals of 136,535,602 and 145,985,440 clean reads, respectively, with 88.88% and 87.99% mapped to the reference genome for gene expression analysis. A total of 3510 differentially expressed genes (DEGs) were identified (padi = 1), with 3203 upregulated DEGs and 307 downregulated DEGs in M5026 compared with ‘FT’ (Fig. 8a). DEGs associated with male sterility included genes involved in both the tapetum and pollen tube. Significantly downregulated tapetum-related genes were identified among the DEGs, including BCP1 ( BraA09g039020.3.5C , BraA08g025680.3.5C and BraA07g012900.3.5C ) and AMS ( BraA03g043930.3.5C and BraA07g004250.3.5C ) (Fig. 8a). Among the DEGs associated with pollen tube function, LRX8 ( BraA01g035230.3.5C and BraA03g039100.3.5C ), LRX9 ( BraA08g003810.3.5C ), LRX10 ( BraA07g005090.3.5C and BraA03g059520.3.5C ), and LRX11 ( BraA01g004310.3.5C , BraA08g016780.3.5C and BraA03g043440.3.5C ) exhibited significant downregulation (Fig. 8a). Brems1 belongs to the LRR-RLK family, and we identified DEGs within this family, including PRK1 ( BraA08g009650.3.5C ), PRK2A ( BraA07g006840.3.5C ), PRK3 ( BraA06g022820.3.5C ), PRK5 ( BraA05g019360.3.5C and BraA08g002990.3.5C ), and PRK8 ( BraA07g030150.3.5C ), which all exhibited significant downregulation (Fig. 8a). Mutations in Brems1 potentially affected the expression levels of genes associated with the tapetum, pollen tube and LRR-RLK family in Chinese cabbage. The KEGG pathway enrichment analysis of these DEGs revealed significant enrichment in multiple metabolic pathways, including arginine and proline metabolism (brp00330), cyanoamino acid metabolism (brp00460), arachidonic acid metabolism (brp00590), nitrogen metabolism (brp00910), amino sugar and nucleotide sugar metabolism (brp00520), galactose metabolism (brp00052), starch and sucrose metabolism (brp00500), taurine and hypotaurine metabolism (brp00430), glycerolipid metabolism (brp00561), ascorbate and aldarate metabolism (brp00053), inositol phosphate metabolism (brp00562) and glycerophospholipid metabolism (brp00564), (Fig. 8b). Therefore, the pollen sterility in M5026 might have influenced its metabolism. Discussion The utilization of male sterile lines has provided a crucial tool for breeding and producing hybrid crop varieties. Moreover, male sterile materials can also be employed to investigate genes involved in anther development. In this study, we identified two male sterile mutants, M5026 and M5073 , in Chinese cabbage and cloned the mutant gene Brems1 ( BraA10g029920.3.5C ). Brems1 encoded an LRR-RLK. According to our knowledge, this study represents the first successful cloning of the EMS1 homologous gene in the genus Brassica , contributing to elucidating the molecular mechanisms underlying male sterility in Chinese cabbage. The abnormal development of the tapetal was a primary cause of male sterility in plants (Wilson et al. 2011 ). The tapetal functioned as the innermost layer of the anther wall and directly contacted the microsporocytes to provide them with nutrients and generated enzymes that degraded sporopollenin to facilitate microspore release (Ma et al. 2015 ). Additionally, it regulated the formation of the pollen exine by secreting sporopollenin precursors (Zou et al. 2023 ). Persistent expansion of tapetum cells and abnormal aggregation of microsporocytes in the GMS line ‘AB01’ in Chinese cabbage led to male sterility (Zhou et al. 2017 ). The vacuolization anomaly of tapetum cells in the male-sterile mutant ftms in Chinese cabbage resulted in anther abortion (Tan et al. 2019 ). The programmed cell death of tapetum cells in the male-sterile mutant msm2 in Chinese cabbage was delayed, leading to an inability to produce pollen (Dong et al. 2022 ). The abnormal development of the tapetum in Chinese cabbage msm1 resulted in male sterility (Zou et al. 2023 ). The abnormal enlargement and vacuolization of tapetum cells in the Chinese cabbage mutant ftms1 resulted in infertility (Zhao et al. 2022 ). The abnormal degradation of the tapetum in the male-sterile mutant msm3 of Chinese cabbage led to complete stamen degeneration (Xu et al. 2024 ). In this study, the mutant M5026 , lacking the tapetum, resulted in complete stamen sterility. The exine is a multilayered structure precisely assembled from sporopollenin, exhibiting characteristics such as radiation protection, high-temperature resistance, and resistance to degradation (Hou et al. 2023 ), and abnormalities in the exine structure led to anther abortion. In the Arabidopsis thaliana male-sterile mutant ms1 , exine formation was aberrant with limited sporopollenin deposition (Vizcay-Barrena and Wilson 2006 ). The rice male-sterile mutant osms1 displayed aberrant exine morphology, presenting a bilayered structure with an absence of bacula (Yang et al. 2019 ). Exine thinning in the rice mutant ms7-6007 results in male sterility (An et al. 2020 ). Severe exine deficiency was observed in the Chinese cabbage male sterile mutant msm1 (Zou et al. 2023 ), and the abnormal exine of the Chinese cabbage mutant msm2 resulted in anther abortion (Dong et al. 2022 ). The tapetum provided molecules for exine formation, known as exine proteins and callose, which constitute the proteinaceous coat of pollen. The abnormal microspore exine structure in the mutant M5026 investigated in this study led to male sterility, possibly due to the absenct tapetum. Receptor protein kinases (RPKs) constituted a transmembrane protein family that played crucial roles in cellular signaling pathways in both animals and plants (Becraft and PW 1998). RPKs were typically composed of three structural domains: the extracellular domain, the transmembrane domain, and the cytosolic protein kinase domain (Walker 1994 ). Based on the presumed ligand-binding domain structures, RPKs in plants were classified into five families: wall-associated kinases (WAKs), the lectin-like RPK family, the CR4-like RPK family, the S-domain RPK family, and the LRR-RPK family (Mccarty and Chory 2000 ). LRR-RPK family genes were implicated in defense responses and diverse growth and developmental processes, including nodulation (Gresshoff et al. 2003 ), plant transpiration (Masle et al. 2005 ), organ shape and size regulation (Xu et al. 2008 ), cell fate determination and patterning in anther development (Jia et al. 2008 ), organ abscission (Kumpf et al. 2013 ), steroid hormone signaling (Santiago et al. 2013 ), nitrogen acquisition (Tabata et al. 2014 ), regulation of root and shoot meristem size (Shinohara et al. 2016 ), defense responses (Zorzatto et al. 2015 ), and pollen tube reception (Takeuchi and Higashiyama 2016 ). In this study, the causative gene Brems1 of the male-sterile mutant encoded an LRR-RLK belonging to the LRR-RLK family, resulting in defective anther development. In conclusion, this study identified two male sterile mutants through EMS-induced mutagenesis. M5026 exhibited an absence of the tapetal layer, an abundance of micropylar cells, and aberrant microspore exine structure. The LRR-RLK gene, Brems1 , homologous to the Arabidopsis male sterility gene EMS1 , was responsible for male sterility in M5026 and M5073 . A SNP G to A mutation was observed in the S_TKc domain of Brems1 in M5026 . A C-to-T SNP was identified within Brems1 of M5073 , resulting in premature translation termination. The present study provided insights into the molecular mechanisms of male sterility in Chinese cabbage. Figure legends Figure 1 Morphological characteristics of M5026 and ‘FT’. a M5026 and ‘FT’ at the flowering stage. Bar, 1 cm. b The flowers of M5026 and ‘FT’. Bar, 3 mm. c Pollen viability analysis of M5026 and ‘FT’. Bar, 50 µm. d - i The floral organs of M5026 and ‘FT’. d, pistils. e, long stamens. f, short stamens. g, petals. h, buds. i, sepals. Bars, 3mm. Figure 2 Paraffin sections during anther development in ‘FT’ ( a - d , i - l ) and M5026 ( e - h , m - p ). Bar = 50 µm. E, epidermis; PT, precursors of tapetal cells; PPT, putative precursors of tapetal cells; T, tapetal layer; MS, microsporocytes; dMS, degraded microsporocytes; Tds, tetrads; Msp, microspore; PG, pollen grain; Se, septum. Figure 3 TEM observation of microspores in ‘FT’ and M5026. a Microspore of ‘FT’. b Localized magnification of microspore in ‘FT’. c Microspore of M5026 . d Localized magnification of microspore in M5026. Bar = 1 µm. Msp, microspore; MS, microsporocytes; Ex, exine; ba, bacula; te, tectum; In, Intine; ab-Ex, abnormal exine. Figure 4 SNP-index distribution plot and KASP genotyping. a SNP-index distribution plot. The red line represented the mean SNP-index within the window, the pink line denotes the threshold line corresponding to the 99th percentile, and the orange line indicated the threshold line corresponding to the 95th percentile. b KASP genotyping of SNPs. C:C corresponds to red dots, C:T to green dots, and T:T to blue dots. Figure 5 The cloning and gene structure of BraA10g029920.3.5C . a The alignment of cloned gene sequences surrounding the mutation site in ‘FT’ and M5026. The red box highlighted the mutation site. b The gene structure of BraA10g029920.3.5C in the reference genome Brara_Chiifu_V3.5, ‘FT’, and M5026 . b1 The gene structure of BraA10g029920.3.5C in the reference genome Brara_Chiifu_V3.5. b2 The gene structure of BraA10g029920.3.5C in ‘FT’, and M5026. The red arrow highlighted the mutation site. Figure 6 Protein analysis of BrEMS1. a Protein domain analysis of BrEMS1. The red arrow indicated the mutation site. b Analysis of the three-dimensional structure of the BrEMS1 protein. b1 The overall three-dimensional structure of the BrEMS1 protein. b2 The local three-dimensional structure of BrEMS1 protein around the mutation site. b3 The overall three-dimensional structure of the Brems1 protein. b4 The local three-dimensional structure of BrEMS1 protein around the mutation site. c The phylogenetic tree of BrEMS1. d Amino acid sequence alignment of the S_TKc domain. Figure 7 Subcellular localization and qRT-PCR of ‘FT’ and M5026. a The subcellular localization of BrEMS1 and Brems1. b qRT-PCR analysis of Brems1 in different organs of ‘FT’ and M5026. Figure 8 Transcriptome analysis of buds in ‘FT’ and M5026. a Transcriptional levels of DEGs associated with the tapetum, pollen tube and LRR-RLK family. b KEGG enrichment analysis of DEGs. The red asterisk indicated metabolic pathways. Fig. S1 The allelic material M5037 of M5026 . a M5037 and ‘FT’ at the flowering stage. Bar, 1 cm. b The flowers of M5037 and ‘FT’. Bar, 3 mm. c SNP-index distribution plot. The red line represented the mean SNP-index within the window, the pink line denotes the threshold line corresponding to the 99th percentile, and the orange line indicated the threshold line corresponding to the 95th percentile. d The alignment of cloned gene sequences surrounding the mutation site in ‘FT’ and M5037. The red box highlighted the mutation site. Declarations Conflict of interest All the authors declare that they have no competing interests. Authors’ contributions CL and LW analyzed the data. CL drafted the manuscript. CL, LW and DZ participated in the creation of materials and performed the experiments. CT and ZL directed the whole study including designing experiments and revising the manuscript. Acknowledgements This work was supported by the Earmarked Fund for CARS-23-A02, the National Natural Science Foundation of China (32272736 and 31972406) and the Key Research and Development Program of Liaoning (2021JH2/10200003). Availability of data and materials The datasets supporting the conclusions of this article are included within the article and its additional fles. The Illumina RNA-Seq datasets and Mutmap data are available in the the Sequence Read Archives (SRA) of the National Center for Biotechnology Information NCBI under BioProject ID: PRJNA1105001 andPRJNA1105243. Genomic sequences and gene annotation information of B.rapa are downloaded online at http://brassicadb.cn . References Abe A, Kosugi S, Yoshida K, Natsume S, Takagi H, Kanzaki H, Matsumura H, Yoshida K, Mitsuoka C, Tamiru M (2012) Genome sequencing reveals agronomically important loci in rice using MutMap. Nat Biotechnol 30:174–178 An X, Ma B, Duan M, Dong Z, Liu R, Yuan D, Hou Q, Wu S, Zhang D, Liu D, Yu D, Zhang Y, Xie K, Zhu T, Li Z, Zhang S, Tian Y, Liu C, Li J, Yuan L, Wan X (2020) Molecular regulation of ZmMs7 required for maize male fertility and development of a dominant male-sterility system in multiple species. Proc Natl Acad Sci U S A 117:23499–23509 Becraft PW (1998) Receptor kinases in plant development. TRENDS PLANT SCI 1998,3(10):384–388 Chen L, Liu Y-G (2014) Male sterility and fertility restoration in crops. Annu Rev Plant Biol 65:579–606 Dong S, Zou J, Fang B, Zhao Y, Shi F, Song G, Huang S, Feng H (2022) Defect in BrMS1, a PHD-finger transcription factor, induces male sterility in ethyl methane sulfonate-mutagenized Chinese cabbage (Brassica rapa L. ssp. pekinensis). Front Plant Sci 13:992391 Gresshoff PM, Men AE, Laniya TS, Searle IR, Carroll B (2003) Long-Distance Signaling in Nodulation Directed by a CLAVATA1-Like Receptor Kinase. Science 299 Hang D, Yan L, Shunchun L, Lei L, Juanyan W, Zhijie Z, Qichun P, Nan L, Xiangqian Z (2021) Effect of extensin-like OsPEX1 on pollen fertility in rice. Yi Chuan 43:271–279 Hord CL, Sun YJ, Pillitteri LJ, Torii KU, Wang H, Zhang S, Ma H (2008) Regulation of Arabidopsis early anther development by the mitogen-activated protein kinases, MPK3 and MPK6, and the ERECTA and related receptor-like kinases. Mol Plant 1:645–658 Hou Q, An X, Ma B, Wu S, Wei X, Yan T, Zhou Y, Zhu T, Xie K, Zhang D, Li Z, Zhao L, Niu C, Long Y, Liu C, Zhao W, Ni F, Li J, Fu D, Yang Z, Wan X (2023) ZmMS1/ZmLBD30-orchestrated transcriptional regulatory networks precisely control pollen exine development. Mol Plant 16:1321–1338 Huang S, Liu Z, Li D, Yao R, Feng H (2016) A new method for generation and screening of Chinese cabbage mutants using isolated microspore culturing and EMS mutagenesis. Euphytica 207:23–33 Ji J, Yang L, Fang Z, Zhuang M, Zhang Y, Lv H, Liu Y, Li Z (2017) Recessive male sterility in cabbage (Brassica oleracea var. capitata) caused by loss of function of BoCYP704B1 due to the insertion of a LTR-retrotransposon. Theoretical Appl Genet 130:1441–1451 Jia G, Liu X, Owen HA, Zhao D (2008) Signaling of cell fate determination by the TPD1 small protein and EMS1 receptor kinase. Proc Natl Acad Sci U S A 105:2220–2225 Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA (2009) Circos: an information aesthetic for comparative genomics. Genome Res 19:1639–1645 Kumpf RP, Shi CL, Larrieu A, Stø IM, Butenko MA, Péret B, Riiser ES, Bennett MJ, Aalen RB (2013) Floral organ abscission peptide IDA and its HAE/HSL2 receptors control cell separation during lateral root emergence. Proc Natl Acad Sci U S A 110:5235–5240 Li H, Durbin R (2010) Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics 26:589–595 Liao Y, Smyth GK, Shi W (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30:923–930 Liu C, Song G, Wang N, Huang S, Gao Y, Fu W, Zhang M, Feng H (2021a) A single SNP in Brcer1 results in wax deficiency in Chinese cabbage (Brassica campestris L. ssp. pekinensis). Sci Hort 282:110019 Liu T, Jiang GQ, Yao XF, Liu CM (2021b) The leucine-rich repeat receptor-like kinase OsERL plays a critical role in anther lobe formation in rice. Biochem Biophys Res Commun 563:85–91 Ma Y, Kang J, Wu J, Zhu Y, Wang X (2015) Identification of tapetum-specific genes by comparing global gene expression of four different male sterile lines in Brassica oleracea. Plant Mol Biol 87:541–554 Mariani C, Beuckeleer MD, Truettner J, Leemans J, Goldberg RB (1990) Induction of male sterility in plants by a chimaeric ribonuclease gene. Nature 347:737–741 Masle J, Gilmore SR, Farquhar GD (2005) The ERECTA gene regulates plant transpiration efficiency in Arabidopsis. Nature 436:866–870 Mccarty DR, Chory J (2000) Conservation and innovation in plant signaling pathways. Cell 103:201–209 McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20:1297–1303 Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5:621–628 Sanders PM, Bui AQ, Weterings K, McIntire K, Hsu Y-C, Lee PY, Truong MT, Beals T, Goldberg R (1999) Anther developmental defects in Arabidopsis thaliana male-sterile mutants. Sex Plant Reprod 11:297–322 Santiago J, Henzler C, Hothorn M (2013) Molecular mechanism for plant steroid receptor activation by somatic embryogenesis co-receptor kinases. Science 341:889–892 Shi YL, Guo SD, Zhang R, Meng ZG, Ren MZ (2014) The role of Somatic embryogenesis receptor-like kinase 1 in controlling pollen production of the Gossypium anther. Mol Biol Rep 41:411–422 Shinohara H, Mori A, Yasue N, Sumida K, Matsubayashi Y (2016) Identification of three LRR-RKs involved in perception of root meristem growth factor in Arabidopsis. Proc Natl Acad Sci U S A 113:3897–3902 Shiu S-H, Karlowski WM, Pan R, Tzeng Y-H, Mayer KF, Li W-H (2004) Comparative analysis of the receptor-like kinase family in Arabidopsis and rice. Plant Cell 16:1220–1234 Shpak ED, Berthiaume CT, Hill EJ, Torii KU (2004) Synergistic interaction of three ERECTA-family receptor-like kinases controls Arabidopsis organ growth and flower development by promoting cell proliferation. Development 131:1491–1501 Tabata R, Sumida K, Yoshii T, Ohyama K, Shinohara H, Matsubayashi Y (2014) Perception of root-derived peptides by shoot LRR-RKs mediates systemic N-demand signaling. Science 346:343–346 Takeuchi H, Higashiyama T (2016) Tip-localized receptors control pollen tube growth and LURE sensing in Arabidopsis. Nature 531:245–248 Tan C, Liu Z, Huang S, Feng H (2019) Mapping of the male sterile mutant gene ftms in Brassica rapa L. ssp. pekinensis via BSR-Seq combined with whole-genome resequencing. Theor Appl Genet 132:355–370 Vizcay-Barrena G, Wilson ZA (2006) Altered tapetal PCD and pollen wall development in the Arabidopsis ms1 mutant. J Exp Bot 57:2709–2717 Walker JC (1994) Structure and function of the receptor-like protein kinases of higher plants. Plant Mol Biol 26:1599–1609 Wang K, Li M, Hakonarson H (2010) ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 38:e164–e164 Wilson ZA, Song J, Taylor B, Yang C (2011) The final split: the regulation of anther dehiscence. J Exp Bot 62:1633–1649 Wu Hm, Wong E, Ogdahl J, Cheung AY (2000) A pollen tube growth-promoting arabinogalactan protein from Nicotiana alata is similar to the tobacco TTS protein. Plant J 22:165–176 Xi X, Wei K, Gao B, Liu J, Liang J, Cheng F, Wang X, Wu J (2018) BrFLC5: a weak regulator of flowering time in Brassica rapa. Theoretical Appl Genet 131:2107–2116 Xu J, Liao R, Xue M, Shang S, Zhou M, Liu Z, Feng H, Huang S (2024) Mutations in BrABCG26, encoding an ATP-binding cassette transporter, are responsible for male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Theor Appl Genet 137:63 Xu SL, Rahman A, Baskin TI, Kieber JJ (2008) Two leucine-rich repeat receptor kinases mediate signaling, linking cell wall biosynthesis and ACC synthase in Arabidopsis. Plant Cell 20:3065–3079 Yamagishi H, Bhat SR (2014) Cytoplasmic male sterility in Brassicaceae crops. Breed Sci 64:38–47 Yan L, Yang M, Guo H, Yang L, Wu J, Li R, Liu P, Lian Y, Zheng X, Yan J (2013) Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol biology 20:1131–1139 Yang L, Qian X, Chen M, Fei Q, Meyers BC, Liang W, Zhang D (2016) Regulatory Role of a Receptor-Like Kinase in Specifying Anther Cell Identity. Plant Physiol 171:2085–2100 Yang S-L, Xie L-F, Mao H-Z, Puah CS, Yang W-C, Jiang L, Sundaresan V, Ye D (2003) Tapetum determinant1 is required for cell specialization in the Arabidopsis anther. Plant Cell 15:2792–2804 Yang Z, Sun L, Zhang P, Zhang Y, Yu P, Liu L, Abbas A, Xiang X, Wu W, Zhan X, Cao L, Cheng S (2019) TDR INTERACTING PROTEIN 3, encoding a PHD-finger transcription factor, regulates Ubisch bodies and pollen wall formation in rice. Plant J 99:844–861 Zhang J, Wu P, Li N, Xu X, Wang S, Chang S, Zhang Y, Wang X, Liu W, Ma Y, Manghwar H, Zhang X, Min L, Guo X (2022) A male-sterile mutant with necrosis-like dark spots on anthers was generated in cotton. Front Plant Sci 13:1102196 Zhang W, Sun Y, Timofejeva L, Chen C, Grossniklaus U, Ma H (2006) Regulation of Arabidopsis tapetum development and function by DYSFUNCTIONAL TAPETUM1 (DYT1) encoding a putative bHLH transcription factor. Development 133:3085–3095 Zhao D, Wang G, Speal B, Ma H (2002) The excess microsporocytes1 gene encodes a putative leucine-rich repeat receptor protein kinase that controls somatic and reproductive cell fates in the Arabidopsis anther. Genes Dev 16:2021–2031 Zhao Y, Huang S, Zou J, Dong S, Wang N, Feng H (2022) Mutation in BrGGL7 gene encoding a GDSL esterase / lipase causes male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Theor Appl Genet 135:3323–3335 Zhou X, Liu Z, Ji R, Feng H (2017) Comparative transcript profiling of fertile and sterile flower buds from multiple-allele-inherited male sterility in Chinese cabbage (Brassica campestris L. ssp. pekinensis). Mol Genet Genomics 292:967–990 Zorzatto C, Machado J, Lopes K, Nascimento K, Pereira W, Brustolini O, Reis P, Calil I, Deguchi M, Sachetto-Martins G, Gouveia B, Loriato V, Silva M, Silva F, Santos A, Chory J, Fontes E (2015) NIK1-mediated translation suppression functions as a plant antiviral immunity mechanism. Nature 520:679–682 Zou J, Dong S, Fang B, Zhao Y, Song G, Xin Y, Huang S, Feng H (2023) BrACOS5 mutations induced male sterility via impeding pollen exine formation in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Theor Appl Genet 136:6 Supplementary Files Fig.S1.tif Fig. S1 The allelic material M5037 of M5026 . a M5037 and ‘FT’ at the flowering stage. Bar, 1 cm. b The flowers of M5037 and ‘FT’. Bar, 3 mm. c SNP-index distribution plot. The red line represented the mean SNP-index within the window, the pink line denotes the threshold line corresponding to the 99th percentile, and the orange line indicated the threshold line corresponding to the 95th percentile. d The alignment of cloned gene sequences surrounding the mutation site in ‘FT’ and M5037. The red box highlighted the mutation site. Supplementarytables.docx Cite Share Download PDF Status: Published Journal Publication published 24 Feb, 2025 Read the published version in Theoretical and Applied Genetics → Version 1 posted Editorial decision: Major revisions 19 Dec, 2024 Reviewers agreed at journal 06 Nov, 2024 Reviewers invited by journal 10 Jul, 2024 Editor assigned by journal 30 May, 2024 First submitted to journal 28 May, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4489236","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":325354056,"identity":"6afc90c0-797b-4d59-920e-5706f4979ece","order_by":0,"name":"Chuanhong Liu","email":"","orcid":"","institution":"Shenyang Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Chuanhong","middleName":"","lastName":"Liu","suffix":""},{"id":325354057,"identity":"aec85ddd-ffe8-406b-9f1a-f2f4304023aa","order_by":1,"name":"Lin Wang","email":"","orcid":"","institution":"Shenyang Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Lin","middleName":"","lastName":"Wang","suffix":""},{"id":325354058,"identity":"7e4aa87b-fb87-472b-a1c3-e1a18f5ad1ee","order_by":2,"name":"Chong Tan","email":"","orcid":"","institution":"Shenyang Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Chong","middleName":"","lastName":"Tan","suffix":""},{"id":325354059,"identity":"c2106f70-784a-4fb3-b790-087841899bce","order_by":3,"name":"Di Zhao","email":"","orcid":"","institution":"Shenyang Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Di","middleName":"","lastName":"Zhao","suffix":""},{"id":325354060,"identity":"5ed0f414-b36c-440c-becb-d20392a9135c","order_by":4,"name":"Zhiyong Liu","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABAElEQVRIiWNgGAWjYBACPmYGBmYGAxCT+eCDD1BRCXxa2BBa2JINZxClBWQ8hMljJs1DlBZ2HuPPBQWH5cz515hJ2/w5LG9wgPngbR4GuzzcDuMxMJ5hcNjYcsazYuvctsOGGw6wJVvzMCQX49OSzGNwOHHDjcMbb+c2HGbccADswgOJDXi0HIZoOWAgbfHnsP2GA/zfCGkxbAZrOd9iJM3ABmQc4GEjoIWtGGhRurHBDWAg97alJ888zGZsOccgGacWfv7Dmz/z/LGWMzh/+OCDH3+sbfuONz+88abCDqcWKGgGxkUClAGOJgP86oGgDmjfAShjFIyCUTAKRgEaAADqQlOBhkDgDAAAAABJRU5ErkJggg==","orcid":"https://orcid.org/0000-0002-7366-6511","institution":"Shenyang Agricultural University","correspondingAuthor":true,"prefix":"","firstName":"Zhiyong","middleName":"","lastName":"Liu","suffix":""}],"badges":[],"createdAt":"2024-05-28 08:17:10","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4489236/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4489236/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s00122-025-04841-y","type":"published","date":"2025-02-24T15:57:55+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":61770269,"identity":"d15548ad-956e-4d4e-a3ee-e1f03eb01e82","added_by":"auto","created_at":"2024-08-05 11:16:23","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1069213,"visible":true,"origin":"","legend":"\u003cp\u003eMorphological characteristics of \u003cem\u003eM5026 \u003c/em\u003eand ‘FT’. \u003cstrong\u003ea\u003c/strong\u003e \u003cem\u003eM5026\u003c/em\u003e and ‘FT’ at the flowering stage. Bar, 1 cm. \u003cstrong\u003eb\u003c/strong\u003e The flowers of \u003cem\u003eM5026 \u003c/em\u003eand ‘FT’. Bar, 3 mm. \u003cstrong\u003ec \u003c/strong\u003ePollen viability analysis of \u003cem\u003eM5026 \u003c/em\u003eand ‘FT’. Bar, 50 μm. \u003cstrong\u003ed - i \u003c/strong\u003eThe floral organs of \u003cem\u003eM5026 \u003c/em\u003eand ‘FT’. d, pistils. e, long stamens. f, short stamens. g, petals. h, buds. i, sepals. Bars, 3mm.\u003c/p\u003e","description":"","filename":"Fig.1.png","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/ea57f2a96091e01f071123c4.png"},{"id":61770270,"identity":"838d82d2-17d2-4374-ae8c-f15aa5f31f81","added_by":"auto","created_at":"2024-08-05 11:16:23","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1343277,"visible":true,"origin":"","legend":"\u003cp\u003eParaffin sections during anther development in ‘FT’ (\u003cstrong\u003ea\u003c/strong\u003e-\u003cstrong\u003ed\u003c/strong\u003e, \u003cstrong\u003ei\u003c/strong\u003e-\u003cstrong\u003el\u003c/strong\u003e) and \u003cem\u003eM5026\u003c/em\u003e(\u003cstrong\u003ee\u003c/strong\u003e-\u003cstrong\u003eh\u003c/strong\u003e, \u003cstrong\u003em\u003c/strong\u003e-\u003cstrong\u003ep\u003c/strong\u003e)\u003cem\u003e. \u003c/em\u003eBar = 50 μm. E, epidermis; PT, precursors of tapetal cells; PPT, putative precursors of tapetal cells; T, tapetal layer; MS, microsporocytes; dMS, degraded microsporocytes; Tds, tetrads; Msp, microspore; PG, pollen grain; Se, septum.\u003c/p\u003e","description":"","filename":"Fig.2.png","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/bfdec4ea8d61e4d253a6ddb2.png"},{"id":61768964,"identity":"6c98f907-4480-4f3c-953b-8b7b00a9bad6","added_by":"auto","created_at":"2024-08-05 11:00:23","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":670779,"visible":true,"origin":"","legend":"\u003cp\u003eTEM observation of microspores in ‘FT’ and \u003cem\u003eM5026. \u003c/em\u003e\u003cstrong\u003ea\u003c/strong\u003e Microspore of ‘FT’. \u003cstrong\u003eb\u003c/strong\u003eLocalized magnification of microspore in ‘FT’. \u003cstrong\u003ec\u003c/strong\u003e Microspore of \u003cem\u003eM5026\u003c/em\u003e. \u003cstrong\u003ed\u003c/strong\u003e Localized magnification of microspore in \u003cem\u003eM5026. \u003c/em\u003eBar = 1 μm. Msp, microspore; MS, microsporocytes; Ex, exine; ba, bacula; te, tectum; In, Intine; ab-Ex, abnormal exine.\u003c/p\u003e","description":"","filename":"Fig.3.png","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/ce1b1e3a3935ae21865b9484.png"},{"id":61769530,"identity":"1ea221ae-a46b-4faf-9ed1-6f7471f53d95","added_by":"auto","created_at":"2024-08-05 11:08:23","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":504765,"visible":true,"origin":"","legend":"\u003cp\u003eSNP-index distribution plot and KASP genotyping. \u003cstrong\u003ea \u003c/strong\u003eSNP-index distribution plot. The red line represented the mean SNP-index within the window, the pink line denotes the threshold line corresponding to the 99th percentile, and the orange line indicated the threshold line corresponding to the 95th percentile. \u003cstrong\u003eb\u003c/strong\u003eKASP genotyping of SNPs. C:C corresponds to red dots, C:T to green dots, and T:T to blue dots.\u003c/p\u003e","description":"","filename":"Fig.4.png","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/59d93c95c9add0a1633aa647.png"},{"id":61768972,"identity":"98dfb868-0b2a-4319-8851-069734884ded","added_by":"auto","created_at":"2024-08-05 11:00:24","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":466582,"visible":true,"origin":"","legend":"\u003cp\u003eThe cloning and gene structure of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e. \u003cstrong\u003ea \u003c/strong\u003eThe alignment of cloned gene sequences surrounding the mutation site in ‘FT’ and \u003cem\u003eM5026. \u003c/em\u003eThe red box highlighted the mutation site. \u003cstrong\u003eb\u003c/strong\u003e The gene structure of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003ein the reference genome Brara_Chiifu_V3.5, ‘FT’, and \u003cem\u003eM5026\u003c/em\u003e. \u003cstrong\u003eb1\u003c/strong\u003eThe gene structure of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e in the reference genome Brara_Chiifu_V3.5. \u003cstrong\u003eb2 \u003c/strong\u003eThe gene structure of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e in ‘FT’, and \u003cem\u003eM5026. \u003c/em\u003eThe red arrow highlighted the mutation site.\u003c/p\u003e","description":"","filename":"Fig.5.png","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/6ae3c18be630aa3fa20f08b3.png"},{"id":61769533,"identity":"ad82cefb-10ed-4adf-970a-a8f270b1b4cc","added_by":"auto","created_at":"2024-08-05 11:08:24","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":472575,"visible":true,"origin":"","legend":"\u003cp\u003eProtein analysis of BrEMS1. \u003cstrong\u003ea \u003c/strong\u003eProtein domain analysis of BrEMS1. The red arrow indicated the mutation site. \u003cstrong\u003eb\u003c/strong\u003e Analysis of the three-dimensional structure of the BrEMS1 protein. \u003cstrong\u003eb1\u003c/strong\u003e The overall three-dimensional structure of the BrEMS1 protein. \u003cstrong\u003eb2\u003c/strong\u003e The local three-dimensional structure of BrEMS1 protein around the mutation site. \u003cstrong\u003eb3 \u003c/strong\u003eThe overall three-dimensional structure of the Brems1 protein. \u003cstrong\u003eb4\u003c/strong\u003e The local three-dimensional structure of BrEMS1 protein around the mutation site. \u003cstrong\u003ec \u003c/strong\u003eThe phylogenetic tree of BrEMS1. \u003cstrong\u003ed \u003c/strong\u003eAmino acid sequence alignment of the S_TKc domain.\u003c/p\u003e","description":"","filename":"Fig.6.png","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/b2b7f82e2b25f2462d6482c0.png"},{"id":61768966,"identity":"6038a3cd-5c92-44be-ae3d-13cd1c1a6d25","added_by":"auto","created_at":"2024-08-05 11:00:23","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":567104,"visible":true,"origin":"","legend":"\u003cp\u003eSubcellular localization and qRT-PCR of ‘FT’ and \u003cem\u003eM5026.\u003c/em\u003e \u003cstrong\u003ea\u003c/strong\u003e The subcellular localization of BrEMS1 and Brems1. \u003cstrong\u003eb\u003c/strong\u003e qRT-PCR analysis of \u003cem\u003eBrems1\u003c/em\u003e in different organs of ‘FT’ and \u003cem\u003eM5026.\u003c/em\u003e\u003c/p\u003e","description":"","filename":"Fig.7.png","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/0f499c441687a7ae34ee56d5.png"},{"id":61768971,"identity":"35ac9023-cd31-45c5-9a93-6b9018294a99","added_by":"auto","created_at":"2024-08-05 11:00:24","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":368317,"visible":true,"origin":"","legend":"\u003cp\u003eTranscriptome analysis of buds in ‘FT’ and \u003cem\u003eM5026. \u003c/em\u003e\u003cstrong\u003ea\u003c/strong\u003e\u003cem\u003e \u003c/em\u003eTranscriptional levels of DEGs associated with the tapetum, pollen tube and LRR-RLK family. \u003cstrong\u003eb\u003c/strong\u003e KEGG enrichment analysis of DEGs. The red asterisk indicated metabolic pathways.\u003c/p\u003e","description":"","filename":"Fig.8.png","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/86b47186e69662431835f385.png"},{"id":77622730,"identity":"7f7a6740-1f15-4d93-a193-a976b6c60b99","added_by":"auto","created_at":"2025-03-03 16:09:50","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":7441470,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/1ef97747-ea4d-478e-8c5f-c02997fc0424.pdf"},{"id":61768969,"identity":"3278b03d-8b60-4138-82ad-c14942d4be86","added_by":"auto","created_at":"2024-08-05 11:00:24","extension":"tif","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":631652,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eFig. S1 \u003c/strong\u003eThe allelic material \u003cem\u003eM5037 \u003c/em\u003eof \u003cem\u003eM5026\u003c/em\u003e. \u003cstrong\u003ea\u003c/strong\u003e \u003cem\u003eM5037\u003c/em\u003e and ‘FT’ at the flowering stage. Bar, 1 cm. \u003cstrong\u003eb\u003c/strong\u003e The flowers of \u003cem\u003eM5037 \u003c/em\u003eand ‘FT’. Bar, 3 mm. \u003cstrong\u003ec\u003c/strong\u003e SNP-index distribution plot. The red line represented the mean SNP-index within the window, the pink line denotes the threshold line corresponding to the 99th percentile, and the orange line indicated the threshold line corresponding to the 95th percentile.\u003cstrong\u003e d \u003c/strong\u003eThe alignment of cloned gene sequences surrounding the mutation site in ‘FT’ and \u003cem\u003eM5037. \u003c/em\u003eThe red box highlighted the mutation site.\u003c/p\u003e","description":"","filename":"Fig.S1.tif","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/fbf508465df6d4409ed9d541.tif"},{"id":61768963,"identity":"cb566ac7-f317-4e35-b0e1-0ae14318e0f5","added_by":"auto","created_at":"2024-08-05 11:00:23","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":23308,"visible":true,"origin":"","legend":"","description":"","filename":"Supplementarytables.docx","url":"https://assets-eu.researchsquare.com/files/rs-4489236/v1/cb8f6238c9873740c7e54235.docx"}],"financialInterests":"","formattedTitle":"Brems1 mutation induced tapetum deficiency leading to male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis)","fulltext":[{"header":"Key Message","content":"\u003cp\u003eThe mutation in \u003cem\u003eBrems1 \u003c/em\u003eresulting in male sterility in Chinese cabbage were validated through two allelic mutations.\u003c/p\u003e"},{"header":"Introduction","content":"\u003cp\u003eMale sterility in plants refers to the phenomenon in which the pistil develops normally while the stamen degenerates and fails to produce viable pollen. This is an effective method for utilizing hybrid vigor and hybrid breeding. Based on the inheritance pattern of sterility genes, male sterility can be classified into cytoplasmic male sterility and nuclear male sterility (Chen and Liu \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Cytoplasmic male sterility was typically caused by an interaction between mitochondrial or chloroplast genes and nuclear genes, and fertility could be restored by nuclear-encoded restorer genes (Yamagishi and Bhat \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Nuclear male sterility was controlled via the nuclear genome, where the molecular regulatory mechanisms of plant male sterility involved genes primarily participating in microspore development, tapetum formation, exine deposition, intine formation, and pollen wall rupture (Wilson et al. \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2011\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe tapetum, located on the innermost layer of the anther wall, directly adjacent to the microspore mother cells, played a crucial role in the nutrient supply and regulating microsporogenesis. After cellular differentiation, the tapetum transformed into polarized secretory cells, synthesizing and secreting callose enzymes via programmed cell death to release uninucleate microspores, while also promoting pollen wall formation and serving as a significant source of nutrition for microspore development (Wu et al. 2000). Plants with mutations that lacked the tapetum or experienced genetic ablation of the tapetum were unable to produce pollen grains (Mariani et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e1990\u003c/span\u003e). The \u003cem\u003eArabidopsis thaliana\u003c/em\u003e male-sterile mutants \u003cem\u003eEMS1\u003c/em\u003e (Zhao et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2002\u003c/span\u003e) and \u003cem\u003eTPD1\u003c/em\u003e (Yang et al. \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2003\u003c/span\u003e), which lacked the tapetum and resulted in male infertility. Premature or delayed degeneration tapetum could also lead to male sterility. Dysfunction of the tapetum in the \u003cem\u003eArabidopsis thaliana\u003c/em\u003e male-sterile mutant \u003cem\u003eDYT1\u003c/em\u003e resulted in abnormal pollen development (Zhang et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2006\u003c/span\u003e). Several mutants with abnormal tapetum development have been identified in Chinese cabbage, including \u003cem\u003eBrMS1\u003c/em\u003e (Dong et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), \u003cem\u003eBrACOS5\u003c/em\u003e (Zou et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), \u003cem\u003eBrGGL7\u003c/em\u003e (Zhao et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), and \u003cem\u003eBrABCG26\u003c/em\u003e (Xu et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe leucine-rich repeat receptor-like kinase (LRR-RLK) consisting of an intracellular domain with kinase activity, a transmembrane domain and extracellular domains responsible for ligand recognition (Shiu et al. \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2004\u003c/span\u003e). Some of the LRR-RLK members family participated in plant stamen development. In \u003cem\u003eArabidopsis\u003c/em\u003e, the receptor-like kinase \u003cem\u003eERECTA\u003c/em\u003e mutant exhibited abnormal flower development in the form of absent anthers (Shpak et al. \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2004\u003c/span\u003e); the triple mutant \u003cem\u003eer-105\u003c/em\u003e, \u003cem\u003eerl1-2\u003c/em\u003e, \u003cem\u003eerl2-1\u003c/em\u003e, encoding an LRR-RLK, resulted in undifferentiated anthers, leading to male sterility (Hord et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2008\u003c/span\u003e); \u003cem\u003eEMS1\u003c/em\u003e gene encoded an LRR-RLK, and its mutant exhibited male sterility due to the absence of middle layer and tapetum cells (Zhao et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2002\u003c/span\u003e). In rice, a mutation in the LRR-RLK gene \u003cem\u003eMSP1\u003c/em\u003e resulted in absent of the middle and tapetum cell layers, leading to male sterility (Yang et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2016\u003c/span\u003e); In addition, a mutation in the \u003cem\u003eERECTA\u003c/em\u003e homologous gene, \u003cem\u003eOsERL\u003c/em\u003e, resulted in significant defects in anther development (Liu et al. \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2021b\u003c/span\u003e); and knocking down the expression level of the LRR-RLK family gene \u003cem\u003eOsPEX1\u003c/em\u003e using RNAi resulted in male sterility (Hang et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). In cotton, plants with RNAi of \u003cem\u003eGhSERK1\u003c/em\u003e exhibited reduced pollen viability (Shi et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), and editing the \u003cem\u003eGhEMS1\u003c/em\u003e gene using CRISPR/Cas9 produced, male-sterile mutants (Zhang et al. \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Currently, no reports detailing the involvement of members of the LRR-RLK family in the molecular mechanisms of male fertility in Chinese cabbage have been made.\u003c/p\u003e \u003cp\u003eHere, two male sterile mutants, \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e, were obtained using ethyl methanesulfonate (EMS)-induced mutagenesis of \u0026lsquo;FT\u0026rsquo;. Cytological observations revealed that \u003cem\u003eM5026\u003c/em\u003e lacked tapetum, exhibited excessive microsporocytes production, and showed abnormal microspore exine morphology. Through MutMap, Kompetitive allele-specific PCR (KASP), and cloning sequencing, the predicted candidate gene was identified as \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e (\u003cem\u003eBrems1\u003c/em\u003e), encoding an LRR-RLK. Employing Mutmap and cloning approaches, a C-to-T SNP was identified within \u003cem\u003eBrems1\u003c/em\u003e of \u003cem\u003eM5073\u003c/em\u003e, resulting in premature translation termination. BrEMS1 and Brems1 were subcellularly localized at the cell membrane. \u003cem\u003eBrems1\u003c/em\u003e exhibited high expression levels in buds. Transcriptomic analysis indicated that the mutation in \u003cem\u003eBrems1\u003c/em\u003e reduced the expression of genes associated with trichomes, pollen tubes and the LRR-RLK family. This study cloned the \u003cem\u003eBrems1\u003c/em\u003e gene in Chinese cabbage for the first time, contributing to elucidating the molecular mechanisms underlying tapetal development and male sterility in Chinese cabbage.\u003c/p\u003e"},{"header":"Materials and methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003ePlant materials\u003c/h2\u003e \u003cp\u003eThe DH line \u0026lsquo;FT\u0026rsquo; of Chinese cabbage obtained via microspore culturing was used as the wild-type (Huang et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). The male-sterile mutants \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e were generated using EMS-induced mutagenesis of \u0026lsquo;FT\u0026rsquo;. Two parents were used to generate the F\u003csub\u003e1\u003c/sub\u003e and F\u003csub\u003e2\u003c/sub\u003e generations for genetic analysis. The F\u003csub\u003e2\u003c/sub\u003e generation was further used for gene mapping and KASP analysis.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eFloral organ observation\u003c/h2\u003e \u003cp\u003eFloral organ morphology (pistils, long stamens, short stamens, petals, sepals, and flower buds) of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e were observed at the reproductive growth stage using a stereomicroscope (Nikon SMZ800, Japan).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003ePollen viability analysis\u003c/h2\u003e \u003cp\u003eThe pollen viabilities of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e were assessed using the 2,3,5-triphenyltetrazolium chloride (TTC) staining method. Anthers were collected from flower buds and placed on a glass slide. 50 \u0026micro;L of TTC (0.2%) staining solution was added to the slide, and a coverslip was placed on top. The slides were incubated at 37\u0026deg;C for 20\u0026ndash;30 min. Pollen viability was observed using a fluorescence inverted microscope (OLYMPUS DP80, Japan).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eParaffin section observation\u003c/h2\u003e \u003cp\u003eThe buds at different developmental stages of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e were individually fixed in Formalin-Aceto-Alcohol (FAA) fixative for 12 hours, followed by dehydration using an ethanol solution gradient of 50%-100%. Subsequently, the samples were subjected to substitution with a mixed solution of ethanol and xylene (volume ratios\u0026thinsp;=\u0026thinsp;3:1, 1:1, 1:3), and different mass ratios of xylene and paraffin were then impregnated into the samples. The samples were embedded in paraffin, and sections were cut using a microtome (Leica RM2016, Germany). After staining the anther sections with safranin and fast green, these sections were observed using an optical microscope (ECLIPSE80i, Nikon).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eTEM observation\u003c/h2\u003e \u003cp\u003eBuds of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e were fixed via immersion in a 1:1 mixture of glutaraldehyde and paraformaldehyde. The fixed samples were then stored in the dark at 4\u0026deg;C. The detailed steps for transmission electron microscopy (TEM) were conducted according to the previously described protocols (Ji et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2017\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eGenetic analysis\u003c/h2\u003e \u003cp\u003eConstructing F\u003csub\u003e1\u003c/sub\u003e and F\u003csub\u003e2\u003c/sub\u003e generations were constructed using parental strains \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo; for genetic analysis. The number of male fertile and sterile plants in the F\u003csub\u003e1\u003c/sub\u003e and F\u003csub\u003e2\u003c/sub\u003e generations was observed and statistically analyzed. The chi-square (χ\u003csup\u003e2\u003c/sup\u003e) test was employed to analyze the segregation ratio in the F\u003csub\u003e2\u003c/sub\u003e population.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003eIdentification of the candidate gene\u003c/h2\u003e \u003cp\u003eCandidate genes in \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e were identified using an improved Mutmap method. (Abe et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). A pool of male sterile mutants was constructed by collecting genomic DNA from 58 F\u003csub\u003e2\u003c/sub\u003e plants exhibiting male sterility. Re-sequencing was performed using NovaSeq 6000 platform (Illumina, San Diego, California, USA) following DNA extraction from \u0026lsquo;FT\u0026rsquo; and male sterile mutant pool using a DNAsecure Plant Kit (Tiangen Biotech Ltd, Beijing, China). The sliding window method implemented in fastp (v0.20.0) was utilized for filtering raw data to obtain high-quality data (Liu et al. \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2021a\u003c/span\u003e). Subsequently, these filtered high-quality data were aligned to the reference genome using bwa (0.7.12-r1039) (Li and Durbin \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Single-nucleotide polymorphisms (SNPs) and insertion-deletions (INDELs) were detected using the Genome Analysis Toolkit (GATK) (McKenna et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). ANNOVAR was employed for functional annotation (Wang et al. \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Circos was used to visualize the mapping of mutation information onto the genome (Krzywinski et al. \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2009\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003eSNP genotype by KASP\u003c/h2\u003e \u003cp\u003eA total of 70 F\u003csub\u003e2\u003c/sub\u003e male sterile individuals and two \u0026lsquo;FT\u0026rsquo; plants were used as materials to conduct a KASP assay to detect SNP co-segregation and identify the \u003cem\u003eBrems1\u003c/em\u003e candidate gene. The thermal cycling conditions for KASP were set as described previously (Xi et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). The primers used in the KASP assay were listed in Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eCloning and sequencing of \u003cem\u003eBrems1\u003c/em\u003e\u003c/h2\u003e \u003cp\u003eThe full-length DNA and coding sequence (CDS) of the \u003cem\u003eBrems1\u003c/em\u003e were amplified from \u0026lsquo;FT\u0026rsquo;, \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e. The specific operational steps were conducted according to a previous study (Liu et al. \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2021a\u003c/span\u003e). Sanger sequencing was performed by Sangon Biotech (Shanghai, China).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003eBioinformatic characterization and phylogenetic analysis of the BrEMS1 protein\u003c/h2\u003e \u003cp\u003eThe full-length amino acid sequence of BrEMS1 was submitted to the SMART (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://smart\u003c/span\u003e\u003cspan address=\"https://smart\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e. embl- heidelberg. de/) website for protein domain analysis. The three-dimensional protein structures of BrEMS1 and Brems1 were analyzed using SwissModel. A total of 1000 bootstrap replicates were established in MEGA6.0 software, and the neighbor-joining method was employed to construct the phylogenetic tree. DNAMAN v6.0 was used to perform multiple sequence alignment of the protein sequences.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eSubcellular localization of BrEMS1 and Brems1\u003c/h2\u003e \u003cp\u003eThe CDS sequences of BrEMS1 and Brems1 were cloned from \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e, respectively, and inserted into the ProCAMV35S vector and transformed into Agrobacterium tumefaciens GV3101, and impregnated into tobacco leaf mesophyll cells at an optical density (OD) of 0.6\u0026ndash;0.8. After 24 hours of dark treatment followed by 24 hours of light treatment, fluorescence signals were observed using a confocal laser scanning microscope (Leica Microsystems, Wetzlar, Germany). Green fluorescence protein (GFP) fluorescence signals were detected at a wavelength of 496\u0026thinsp;~\u0026thinsp;540 nm.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003eQuantitative Real-Time PCR (qRT-PCR)\u003c/h2\u003e \u003cp\u003eThe total RNA isolated from flowers, leaves, roots, buds, and stems of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e plants was subsequently reverse-transcribed into cDNA. qRT-PCR was performed using Ultra SYBR Green Mix (Kangwei Century, Beijing, China) on a QuantStudio 6 PCR system. Actin was selected as the internal control. The relative expression levels were calculated using the 2\u003csup\u003e\u0026minus;ΔΔCt\u003c/sup\u003e method.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003eTranscriptome analysis\u003c/h2\u003e \u003cp\u003eBuds from \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e were individually sampled in triplicate, followed by RNA extraction. Subsequently, the RNA samples were subjected to rigorous quality control using an Agilent 2100 bioanalyzer, and libraries were constructed accordingly. The concentration of the libraries was first preliminarily quantified using a Qubit 2.0 Fluorometer, and then the insert size of the libraries was detected using an Agilent 2100 bioanalyzer. The effective concentrations of the libraries were accurately quantified via qRT-PCR to ensure the quality of the libraries. Once the libraries passed quality control, different libraries were pooled based on their effective concentrations and the desired amount of data for sequencing on the Illumina platform. These raw sequencing data were then filtered as follows: reads with adapters were removed; reads containing N (representing unknown bases) were removed; and low-quality reads (more than 50% of the bases with Qphred scores less than or equal to 5) were removed (Yan et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). These sequencing data were aligned to the reference genome Brara_Chiifu_V3.5 using HISAT2 software to obtain a rapid and precise mapping of the reads to the reference genome (Mortazavi et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2008\u003c/span\u003e). Based on the positional information of gene alignment on the reference genome, reads covering each gene within the start and end range were counted. Reads with alignment quality scores lower than 10, unpaired alignments, and reads aligned to multiple regions of the genome were filtered out. This analysis was conducted using the featureCounts tool within the subread software (Liao et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). After gene expression quantification was completed, a statistical analysis was conducted to identify genes that exhibited significant differences in expression levels across different conditions. The differential gene set was subjected to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis using the clusterProfiler software.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003eThe \u003cem\u003eM5026\u003c/em\u003e mutant displayed male sterility\u003c/h2\u003e \u003cp\u003eThe male sterility \u003cem\u003eM5026\u003c/em\u003e mutant was isolated from the EMS-induced mutant DH line \u0026lsquo;FT\u0026rsquo;. Pollen grains were absent on the anther surface of the \u003cem\u003eM5026\u003c/em\u003e mutant (Fig.\u0026nbsp;1a, b). Pollen viability analysis revealed that \u0026lsquo;FT\u0026rsquo; produced viable pollen grains. In contrast, \u003cem\u003eM5026\u003c/em\u003e did not produce any pollen grains (Fig.\u0026nbsp;1c). Observations of the floral organs of \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo; plants revealed that the long stamens, short stamens, pistils, petals, and buds of \u003cem\u003eM5026\u003c/em\u003e mutant were all shorter than those of \u0026lsquo;FT\u0026rsquo; (Fig.\u0026nbsp;1d-h), while the sepals did not different (Fig.\u0026nbsp;1i).\u003c/p\u003e \u003cp\u003eThe F\u003csub\u003e1\u003c/sub\u003e and F\u003csub\u003e2\u003c/sub\u003e generations were constructed using the parental strains \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo; to conduct genetic analysis. All F\u003csub\u003e1\u003c/sub\u003e plants showed the male fertile phenotype. The F\u003csub\u003e2\u003c/sub\u003e population had 581 male fertile and 187 male sterile plants, a ratio of 3:1 (χ\u003csup\u003e2\u003c/sup\u003e = 0.17) (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). A single recessive nuclear gene, \u003cem\u003eBrems1\u003c/em\u003e, regulated the male sterility phenotype in \u003cem\u003eM5026\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eGenetic analysis in \u003cem\u003eM5026.\u003c/em\u003e\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"7\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGenerations\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eTotal\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eMale\u0026nbsp;fertility\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMale\u0026nbsp;sterile\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eSegregation\u0026nbsp;ratio\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eχ\u003csup\u003e2\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eχ\u003csup\u003e2\u003c/sup\u003e\u003csub\u003e0.05\u003c/sub\u003e\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eP\u003csub\u003e1\u003c/sub\u003e (\u0026lsquo;FT\u0026rsquo;)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eP\u003csub\u003e2\u003c/sub\u003e (\u003cem\u003eM5026\u003c/em\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eF\u003csub\u003e1\u003c/sub\u003e (P\u003csub\u003e2\u003c/sub\u003e\u0026times;\u0026nbsp;P\u003csub\u003e1\u003c/sub\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eF\u003csub\u003e1\u003c/sub\u003e\u0026rsquo;(P\u003csub\u003e1\u003c/sub\u003e\u0026times;\u0026nbsp;P\u003csub\u003e2\u003c/sub\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e20\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eF\u003csub\u003e2\u003c/sub\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e768\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e581\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e187\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e3.10:1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.17\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e3.841\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eCandidate SNP information.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"8\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eChromosome\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePosition\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eWT\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMUT\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eSnpindex\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eRegion\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eMutationType\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c8\"\u003e \u003cp\u003eGeneID\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eA10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e19,022,513\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eT\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eintronic\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e/\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eA10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e19,087,850\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eT\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eexonic\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003enonsynonymous SNV\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cem\u003eBraA10g030050.3.5C\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eA10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e19,417,177\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eT\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eexonic\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003esynonymous SNV\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003e\u003cem\u003eBraA10g030900.3.5C\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003eTapetum deficiency, excessive microsporocytes, and abnormal exine in the Mutant \u003cem\u003eM5026\u003c/em\u003e\u003c/h2\u003e \u003cp\u003eTo further elucidate the cytological characteristics of the sterile mutant, we employed paraffin sectioning to examine the anthers of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e at different developmental stages. Based on the established categorization of \u003cem\u003eArabidopsis\u003c/em\u003e anther developmental stages (Sanders et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e1999\u003c/span\u003e), we selected eight anther developmental stages of anther development in both \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e for analysis. During stage 4 (S4) of anther development, no differences were observed between \u0026lsquo;FT\u0026rsquo; (Fig.\u0026nbsp;2a) and \u003cem\u003eM5026\u003c/em\u003e (Fig.\u0026nbsp;2e). At stage 5 (S5), \u0026lsquo;FT\u0026rsquo; developed the tapetum characterized by rectangular-shaped cells forming a layer outside the microsporocytes (Fig.\u0026nbsp;2b), however, \u003cem\u003eM5026\u003c/em\u003e lacked the tapetum and exhibited excess microsporocytes (Fig.\u0026nbsp;2f). At stage 6 (S6), the microsporocytes of \u003cem\u003eM5026\u003c/em\u003e exhibited morphological abnormalities characterized by concavity (Fig.\u0026nbsp;2g). During stages 7\u0026ndash;12 (S7-12), the microsporocytes of \u0026lsquo;FT\u0026rsquo; sequentially formed tetrads, microspores, and pollen grains (Fig.\u0026nbsp;2d, i-l). In contrast, \u003cem\u003eM5026\u003c/em\u003e continued to lack the tapetum, and its microsporocytes underwent continuous degradation, ultimately failing to produce pollen grains (Fig.\u0026nbsp;2h, m-p).\u003c/p\u003e \u003cp\u003eTEM was utilized to observe further \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e. The microspore exine structure of \u0026lsquo;FT\u0026rsquo; was intact, with well-defined bacula and tectum morphology (Fig.\u0026nbsp;3a, b), while the microspore of \u003cem\u003eM5026\u003c/em\u003e exhibited an abnormal exine structure (Fig.\u0026nbsp;3c, d).\u003c/p\u003e \u003cp\u003e \u003cb\u003eIdentification of candidate gene contributing to the\u003c/b\u003e \u003cb\u003eM5026\u003c/b\u003e \u003cb\u003emale sterility phenotype.\u003c/b\u003e\u003c/p\u003e \u003cp\u003eThe adapted MutMap approach was employed to isolate candidate genes. Sequencing data filtering was performed to obtain high-quality data. The wild-type and F\u003csub\u003e2\u003c/sub\u003e mutant pool yielded 175,935,758 (95.25%) and 177,453,976 (96.23%) high quality dates, respectively. Of these, 98.98% and 99.21% could be mapped to the Brara_Chiifu_V3.5 reference genome. SNPs were detected and filtered using the GATK software, resulting in a total of 732,795 SNPs. Subsequently, SNP-index calculation was conducted, followed by plotting the distribution of the SNP-index across the chromosome (Fig.\u0026nbsp;4a). In the case of selecting regions with SNP-index values exceeding the 95th percentile as candidate regions associated with the \u003cem\u003eM5026\u003c/em\u003e male sterility trait, the candidate gene, \u003cem\u003eBrems1\u003c/em\u003e, was mapped to the region of chromosome A10 from 18,400,000 to 20,725,698. A total of 30 SNPs were found in the candidate interval, among which three have SNP-index values of 1 and were located within genes, namely SNP19,022,513, SNP19,087,850, and SNP19,417,177. SNP19,022,513 was located in an intron, SNP19,087,850 was a nonsynonymous mutation located in an exon, and SNP19,417,177 was a synonymous mutation located in an exon.\u003c/p\u003e \u003cp\u003eTo validate the causal SNP for male sterility, KASP genotyping was performed on three SNPs using two \u0026lsquo;FT\u0026rsquo; and seventy F\u003csub\u003e2\u003c/sub\u003e male sterile individuals. The results showed that only SNP 19,022,513 co-segregated with the male sterility phenotype, with all \u0026lsquo;FT\u0026rsquo; plants being C:C and all mutant phenotype plants being T:T. Recombination was identified at SNP 19,087,850, with one C:C and one C:T genotype observed among the mutant phenotype plants. Similarly, at SNP 19,022,513, recombination was also detected, revealing one C:C and two C:T genotypes in the plants exhibiting the mutant phenotype (Fig.\u0026nbsp;4b). SNP 19,022,513 was located on \u003cem\u003eBraA10g029920.3.5C.\u003c/em\u003e Its \u003cem\u003eArabidopsis\u003c/em\u003e ortholog is \u003cem\u003eEMS1\u003c/em\u003e (\u003cem\u003eAT5G07280\u003c/em\u003e), which encodes EXCESS MICROSPOROCYTES1 (EMS1), an LRR-RLK that controls the fate of somatic and germ cells in \u003cem\u003eArabidopsis\u003c/em\u003e anthers.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section2\"\u003e \u003ch2\u003eCloning and sequence analysis of \u003cem\u003eBrems1\u003c/em\u003e\u003c/h2\u003e \u003cp\u003eThe DNA and cDNA sequences of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e were cloned from \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo;, respectively. The results revealed that \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e has a full length of 3699 bp, comprising of a single exon. A non-synonymous mutation was identified at position 3439 bp, where a G to A transition occurred, changing GGG to AGG, resulting in the substitution of glycine (Gly) with arginine (Arg) (Fig.\u0026nbsp;5a).\u003c/p\u003e \u003cp\u003eThe gene sequence of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e cloned from our experimental materials exhibited differences compared with the reference genome Brara_Chiifu_V3.5. In the reference genome, the full length of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e was 4613 bp (Fig.\u0026nbsp;5b1), which was 911 bp longer than in \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e, with the additional length located at the gene\u0026rsquo;s 5\u0026rsquo; end of the gene (Fig.\u0026nbsp;5b2). \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e in the reference genome contained three exons and two introns (Fig.\u0026nbsp;5b1), while those in \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e had only one exon and no introns (Fig.\u0026nbsp;5b2). The mutation site in \u003cem\u003eM5026\u003c/em\u003e was located in the second intron of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e in the reference genome and in the exon of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e (Fig.\u0026nbsp;5b). The \u003cem\u003eEMS1\u003c/em\u003e mutation in \u003cem\u003eArabidopsis thaliana\u003c/em\u003e resulted in male sterility (Zhao et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2002\u003c/span\u003e). Hence, we presumed \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e was the candidate gene for \u003cem\u003eBrems1\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eIdentification of the allelic material\u003c/b\u003e \u003cb\u003eM5073\u003c/b\u003e\u003c/p\u003e \u003cp\u003eA male sterile mutant \u003cem\u003eM507\u003c/em\u003e3 was induced by EMS mutagenesis in \u0026lsquo;FT\u0026rsquo; (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003ea). \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e were crossed with \u0026lsquo;FT\u0026rsquo; to generate F\u003csub\u003e1\u003c/sub\u003e hybrids, respectively. Allelism tests were performed by crossing these two F\u003csub\u003e1\u003c/sub\u003e generations. The observed ratio of male fertility to sterility in the progeny was 3:1 (χ\u003csup\u003e2\u003c/sup\u003e = 1.405) (Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e), indicating that the male sterility traits in \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e were controlled by a single allelic gene. Subsequently, Mutmap method was employed to isolate candidate genes of \u003cem\u003eM5073\u003c/em\u003e. Regions with SNP index values exceeding the 95th percentile were selected as candidate intervals, a SNP with an SNP-index value of 1 was identified at position 19,024,244 bp on chromosome A10, located within the \u003cem\u003eBrems1\u003c/em\u003e gene (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003eb). Cloning of the \u003cem\u003eBrems1\u003c/em\u003e gene from \u003cem\u003eM5073\u003c/em\u003e revealed a C-to-T mutation at position 1708 bp, converting CAG to TAG, resulting in the substitution of glutamine for a stop codon and premature translation termination (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003ec). These results indicated that the mutation in \u003cem\u003eBrems1\u003c/em\u003e caused the male sterile phenotype observed in \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e, confirming its role in stamen development.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec20\" class=\"Section2\"\u003e \u003ch2\u003ePhylogenetic and structural analysis of BrEMS1\u003c/h2\u003e \u003cp\u003eTo further investigate the role of BrEMS1 in male fertility of Chinese cabbage, the structure of the BrEMS1 protein was analyzed. The BrEMS1 protein comprised 1233 amino acids, with a theoretical isoelectric point of 6.62 and a molecular weight of 133.2 kDa. It harbored a serine or threonine-specific kinase subfamily (S_TKc) along with 17 LRR structural domains (Fig.\u0026nbsp;6a). The mutation in \u003cem\u003eM5026\u003c/em\u003e occurred at the 1147th amino acid of the BrEMS1 protein, where Gly was substituted with Arg, located within the S_TKc domain. The SwissModel tool was employed to predict the protein structures of both BrEMS1 and Brems1, revealing that the mutation in BrEMS1 altered the three-dimensional structure of the protein (Fig.\u0026nbsp;6b). The evolutionary relationships among BrEMS1 homologous sequences across different species were elucidated through phylogenetic analysis, demonstrating its high homology among multiple crops in the \u003cem\u003eBrassicaceae\u003c/em\u003e family (Fig.\u0026nbsp;6c). The S_TKc domain was highly conserved across multiple crops in the \u003cem\u003eBrassicaceae\u003c/em\u003e family (Fig.\u0026nbsp;6d).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec21\" class=\"Section2\"\u003e \u003ch2\u003eExpression pattern analysis of \u003cem\u003eBrems1\u003c/em\u003e\u003c/h2\u003e \u003cp\u003eTo elucidate the expression pattern of \u003cem\u003eBrems1\u003c/em\u003e, we conducted a subcellular localization and fluorescence quantitative analysis of \u003cem\u003eBrems1\u003c/em\u003e. BrEMS1-GFP and Brems1-GFP were introduced into tobacco and their co-localization with the cell membrane marker (Pip2a-RFP) was analyzed. The results demonstrated that the GFP fluorescence of both BrEMS1-GFP and Brems1-GFP was localized on the cell membrane (Fig.\u0026nbsp;7a), indicating that the mutation of BrEMS1 did not affect on its subcellular localization.\u003c/p\u003e \u003cp\u003eThe relative transcription levels of \u003cem\u003eBrems1\u003c/em\u003e in different organs of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e were assessed using qRT-PCR. qRT-PCR analysis revealed no significant differences in \u003cem\u003et\u003c/em\u003ehe relative abundances of \u003cem\u003eBrems1\u003c/em\u003e transcripts among the flowers, leaves, roots, buds, and stems of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e plants. Notably, \u003cem\u003eBrems1\u003c/em\u003e expression was highest in buds, followed by flowers, leaves, and stems, while expression was absent in roots (Fig.\u0026nbsp;7b).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec22\" class=\"Section2\"\u003e \u003ch2\u003e\u003cem\u003eBrems1\u003c/em\u003e affects the transcription levels of genes associated with the tapetum and pollen tube\u003c/h2\u003e \u003cp\u003eTo elucidate the molecular mechanisms underlying \u003cem\u003eBrems1\u003c/em\u003e regulation of male sterility, transcriptomic analysis was conducted on the buds of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e. \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026\u003c/em\u003e obtained totals of 136,535,602 and 145,985,440 clean reads, respectively, with 88.88% and 87.99% mapped to the reference genome for gene expression analysis. A total of 3510 differentially expressed genes (DEGs) were identified (padi\u0026thinsp;\u0026lt;\u0026thinsp;=\u0026thinsp;0.05, log2FoldChange\u0026thinsp;\u0026gt;\u0026thinsp;=\u0026thinsp;1), with 3203 upregulated DEGs and 307 downregulated DEGs in \u003cem\u003eM5026\u003c/em\u003e compared with \u0026lsquo;FT\u0026rsquo; (Fig.\u0026nbsp;8a).\u003c/p\u003e \u003cp\u003eDEGs associated with male sterility included genes involved in both the tapetum and pollen tube. Significantly downregulated tapetum-related genes were identified among the DEGs, including \u003cem\u003eBCP1\u003c/em\u003e (\u003cem\u003eBraA09g039020.3.5C\u003c/em\u003e, \u003cem\u003eBraA08g025680.3.5C\u003c/em\u003e and \u003cem\u003eBraA07g012900.3.5C\u003c/em\u003e) and \u003cem\u003eAMS\u003c/em\u003e (\u003cem\u003eBraA03g043930.3.5C\u003c/em\u003e and \u003cem\u003eBraA07g004250.3.5C\u003c/em\u003e) (Fig.\u0026nbsp;8a). Among the DEGs associated with pollen tube function, \u003cem\u003eLRX8\u003c/em\u003e (\u003cem\u003eBraA01g035230.3.5C\u003c/em\u003e and \u003cem\u003eBraA03g039100.3.5C\u003c/em\u003e), \u003cem\u003eLRX9\u003c/em\u003e (\u003cem\u003eBraA08g003810.3.5C\u003c/em\u003e), \u003cem\u003eLRX10\u003c/em\u003e (\u003cem\u003eBraA07g005090.3.5C\u003c/em\u003e and \u003cem\u003eBraA03g059520.3.5C\u003c/em\u003e), and \u003cem\u003eLRX11\u003c/em\u003e (\u003cem\u003eBraA01g004310.3.5C\u003c/em\u003e, \u003cem\u003eBraA08g016780.3.5C\u003c/em\u003e and \u003cem\u003eBraA03g043440.3.5C\u003c/em\u003e) exhibited significant downregulation (Fig.\u0026nbsp;8a). \u003cem\u003eBrems1\u003c/em\u003e belongs to the LRR-RLK family, and we identified DEGs within this family, including \u003cem\u003ePRK1\u003c/em\u003e (\u003cem\u003eBraA08g009650.3.5C\u003c/em\u003e), \u003cem\u003ePRK2A\u003c/em\u003e (\u003cem\u003eBraA07g006840.3.5C\u003c/em\u003e), \u003cem\u003ePRK3\u003c/em\u003e (\u003cem\u003eBraA06g022820.3.5C\u003c/em\u003e), \u003cem\u003ePRK5\u003c/em\u003e (\u003cem\u003eBraA05g019360.3.5C\u003c/em\u003e and \u003cem\u003eBraA08g002990.3.5C\u003c/em\u003e), and \u003cem\u003ePRK8\u003c/em\u003e (\u003cem\u003eBraA07g030150.3.5C\u003c/em\u003e), which all exhibited significant downregulation (Fig.\u0026nbsp;8a). Mutations in \u003cem\u003eBrems1\u003c/em\u003e potentially affected the expression levels of genes associated with the tapetum, pollen tube and LRR-RLK family in Chinese cabbage.\u003c/p\u003e \u003cp\u003eThe KEGG pathway enrichment analysis of these DEGs revealed significant enrichment in multiple metabolic pathways, including arginine and proline metabolism (brp00330), cyanoamino acid metabolism (brp00460), arachidonic acid metabolism (brp00590), nitrogen metabolism (brp00910), amino sugar and nucleotide sugar metabolism (brp00520), galactose metabolism (brp00052), starch and sucrose metabolism (brp00500), taurine and hypotaurine metabolism (brp00430), glycerolipid metabolism (brp00561), ascorbate and aldarate metabolism (brp00053), inositol phosphate metabolism (brp00562) and glycerophospholipid metabolism (brp00564), (Fig.\u0026nbsp;8b). Therefore, the pollen sterility in M5026 might have influenced its metabolism.\u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eThe utilization of male sterile lines has provided a crucial tool for breeding and producing hybrid crop varieties. Moreover, male sterile materials can also be employed to investigate genes involved in anther development. In this study, we identified two male sterile mutants, \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e, in Chinese cabbage and cloned the mutant gene \u003cem\u003eBrems1\u003c/em\u003e (\u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e). \u003cem\u003eBrems1\u003c/em\u003e encoded an LRR-RLK. According to our knowledge, this study represents the first successful cloning of the \u003cem\u003eEMS1\u003c/em\u003e homologous gene in the genus \u003cem\u003eBrassica\u003c/em\u003e, contributing to elucidating the molecular mechanisms underlying male sterility in Chinese cabbage.\u003c/p\u003e \u003cp\u003eThe abnormal development of the tapetal was a primary cause of male sterility in plants (Wilson et al. \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). The tapetal functioned as the innermost layer of the anther wall and directly contacted the microsporocytes to provide them with nutrients and generated enzymes that degraded sporopollenin to facilitate microspore release (Ma et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Additionally, it regulated the formation of the pollen exine by secreting sporopollenin precursors (Zou et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Persistent expansion of tapetum cells and abnormal aggregation of microsporocytes in the GMS line \u0026lsquo;AB01\u0026rsquo; in Chinese cabbage led to male sterility (Zhou et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). The vacuolization anomaly of tapetum cells in the male-sterile mutant \u003cem\u003eftms\u003c/em\u003e in Chinese cabbage resulted in anther abortion (Tan et al. \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). The programmed cell death of tapetum cells in the male-sterile mutant \u003cem\u003emsm2\u003c/em\u003e in Chinese cabbage was delayed, leading to an inability to produce pollen (Dong et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The abnormal development of the tapetum in Chinese cabbage \u003cem\u003emsm1\u003c/em\u003e resulted in male sterility (Zou et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). The abnormal enlargement and vacuolization of tapetum cells in the Chinese cabbage mutant \u003cem\u003eftms1\u003c/em\u003e resulted in infertility (Zhao et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The abnormal degradation of the tapetum in the male-sterile mutant \u003cem\u003emsm3\u003c/em\u003e of Chinese cabbage led to complete stamen degeneration (Xu et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). In this study, the mutant \u003cem\u003eM5026\u003c/em\u003e, lacking the tapetum, resulted in complete stamen sterility.\u003c/p\u003e \u003cp\u003eThe exine is a multilayered structure precisely assembled from sporopollenin, exhibiting characteristics such as radiation protection, high-temperature resistance, and resistance to degradation (Hou et al. \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), and abnormalities in the exine structure led to anther abortion. In the \u003cem\u003eArabidopsis thaliana\u003c/em\u003e male-sterile mutant \u003cem\u003ems1\u003c/em\u003e, exine formation was aberrant with limited sporopollenin deposition (Vizcay-Barrena and Wilson \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2006\u003c/span\u003e). The rice male-sterile mutant \u003cem\u003eosms1\u003c/em\u003e displayed aberrant exine morphology, presenting a bilayered structure with an absence of bacula (Yang et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Exine thinning in the rice mutant \u003cem\u003ems7-6007\u003c/em\u003e results in male sterility (An et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Severe exine deficiency was observed in the Chinese cabbage male sterile mutant \u003cem\u003emsm1\u003c/em\u003e (Zou et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), and the abnormal exine of the Chinese cabbage mutant \u003cem\u003emsm2\u003c/em\u003e resulted in anther abortion (Dong et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The tapetum provided molecules for exine formation, known as exine proteins and callose, which constitute the proteinaceous coat of pollen. The abnormal microspore exine structure in the mutant \u003cem\u003eM5026\u003c/em\u003e investigated in this study led to male sterility, possibly due to the absenct tapetum.\u003c/p\u003e \u003cp\u003eReceptor protein kinases (RPKs) constituted a transmembrane protein family that played crucial roles in cellular signaling pathways in both animals and plants (Becraft and PW 1998). RPKs were typically composed of three structural domains: the extracellular domain, the transmembrane domain, and the cytosolic protein kinase domain (Walker \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e1994\u003c/span\u003e). Based on the presumed ligand-binding domain structures, RPKs in plants were classified into five families: wall-associated kinases (WAKs), the lectin-like RPK family, the CR4-like RPK family, the S-domain RPK family, and the LRR-RPK family (Mccarty and Chory \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2000\u003c/span\u003e). LRR-RPK family genes were implicated in defense responses and diverse growth and developmental processes, including nodulation (Gresshoff et al. \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2003\u003c/span\u003e), plant transpiration (Masle et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2005\u003c/span\u003e), organ shape and size regulation (Xu et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2008\u003c/span\u003e), cell fate determination and patterning in anther development (Jia et al. \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2008\u003c/span\u003e), organ abscission (Kumpf et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2013\u003c/span\u003e), steroid hormone signaling (Santiago et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2013\u003c/span\u003e), nitrogen acquisition (Tabata et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), regulation of root and shoot meristem size (Shinohara et al. \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), defense responses (Zorzatto et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2015\u003c/span\u003e), and pollen tube reception (Takeuchi and Higashiyama \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). In this study, the causative gene \u003cem\u003eBrems1\u003c/em\u003e of the male-sterile mutant encoded an LRR-RLK belonging to the LRR-RLK family, resulting in defective anther development.\u003c/p\u003e \u003cp\u003eIn conclusion, this study identified two male sterile mutants through EMS-induced mutagenesis. \u003cem\u003eM5026\u003c/em\u003e exhibited an absence of the tapetal layer, an abundance of micropylar cells, and aberrant microspore exine structure. The LRR-RLK gene, \u003cem\u003eBrems1\u003c/em\u003e, homologous to the \u003cem\u003eArabidopsis\u003c/em\u003e male sterility gene \u003cem\u003eEMS1\u003c/em\u003e, was responsible for male sterility in \u003cem\u003eM5026\u003c/em\u003e and \u003cem\u003eM5073\u003c/em\u003e. A SNP G to A mutation was observed in the S_TKc domain of \u003cem\u003eBrems1\u003c/em\u003e in \u003cem\u003eM5026\u003c/em\u003e. A C-to-T SNP was identified within \u003cem\u003eBrems1\u003c/em\u003e of \u003cem\u003eM5073\u003c/em\u003e, resulting in premature translation termination. The present study provided insights into the molecular mechanisms of male sterility in Chinese cabbage.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure legends\u003c/b\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u0026nbsp;1\u003c/b\u003e Morphological characteristics of \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo;. \u003cb\u003ea\u003c/b\u003e \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo; at the flowering stage. Bar, 1 cm. \u003cb\u003eb\u003c/b\u003e The flowers of \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo;. Bar, 3 mm. \u003cb\u003ec\u003c/b\u003e Pollen viability analysis of \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo;. Bar, 50 \u0026micro;m. \u003cb\u003ed - i\u003c/b\u003e The floral organs of \u003cem\u003eM5026\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo;. d, pistils. e, long stamens. f, short stamens. g, petals. h, buds. i, sepals. Bars, 3mm.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u0026nbsp;2\u003c/b\u003e Paraffin sections during anther development in \u0026lsquo;FT\u0026rsquo; (\u003cb\u003ea\u003c/b\u003e-\u003cb\u003ed\u003c/b\u003e, \u003cb\u003ei\u003c/b\u003e-\u003cb\u003el\u003c/b\u003e) and \u003cem\u003eM5026\u003c/em\u003e (\u003cb\u003ee\u003c/b\u003e-\u003cb\u003eh\u003c/b\u003e, \u003cb\u003em\u003c/b\u003e-\u003cb\u003ep\u003c/b\u003e). Bar =\u0026thinsp;50 \u0026micro;m. E, epidermis; PT, precursors of tapetal cells; PPT, putative precursors of tapetal cells; T, tapetal layer; MS, microsporocytes; dMS, degraded microsporocytes; Tds, tetrads; Msp, microspore; PG, pollen grain; Se, septum.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u0026nbsp;3\u003c/b\u003e TEM observation of microspores in \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026.\u003c/em\u003e \u003cb\u003ea\u003c/b\u003e Microspore of \u0026lsquo;FT\u0026rsquo;. \u003cb\u003eb\u003c/b\u003e Localized magnification of microspore in \u0026lsquo;FT\u0026rsquo;. \u003cb\u003ec\u003c/b\u003e Microspore of \u003cem\u003eM5026\u003c/em\u003e. \u003cb\u003ed\u003c/b\u003e Localized magnification of microspore in \u003cem\u003eM5026.\u003c/em\u003e Bar =\u0026thinsp;1 \u0026micro;m. Msp, microspore; MS, microsporocytes; Ex, exine; ba, bacula; te, tectum; In, Intine; ab-Ex, abnormal exine.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u0026nbsp;4\u003c/b\u003e SNP-index distribution plot and KASP genotyping. \u003cb\u003ea\u003c/b\u003e SNP-index distribution plot. The red line represented the mean SNP-index within the window, the pink line denotes the threshold line corresponding to the 99th percentile, and the orange line indicated the threshold line corresponding to the 95th percentile. \u003cb\u003eb\u003c/b\u003e KASP genotyping of SNPs. C:C corresponds to red dots, C:T to green dots, and T:T to blue dots.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u0026nbsp;5\u003c/b\u003e The cloning and gene structure of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e. \u003cb\u003ea\u003c/b\u003e The alignment of cloned gene sequences surrounding the mutation site in \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026.\u003c/em\u003e The red box highlighted the mutation site. \u003cb\u003eb\u003c/b\u003e The gene structure of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e in the reference genome Brara_Chiifu_V3.5, \u0026lsquo;FT\u0026rsquo;, and \u003cem\u003eM5026\u003c/em\u003e. \u003cb\u003eb1\u003c/b\u003e The gene structure of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e in the reference genome Brara_Chiifu_V3.5. \u003cb\u003eb2\u003c/b\u003e The gene structure of \u003cem\u003eBraA10g029920.3.5C\u003c/em\u003e in \u0026lsquo;FT\u0026rsquo;, and \u003cem\u003eM5026.\u003c/em\u003e The red arrow highlighted the mutation site.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u0026nbsp;6\u003c/b\u003e Protein analysis of BrEMS1. \u003cb\u003ea\u003c/b\u003e Protein domain analysis of BrEMS1. The red arrow indicated the mutation site. \u003cb\u003eb\u003c/b\u003e Analysis of the three-dimensional structure of the BrEMS1 protein. \u003cb\u003eb1\u003c/b\u003e The overall three-dimensional structure of the BrEMS1 protein. \u003cb\u003eb2\u003c/b\u003e The local three-dimensional structure of BrEMS1 protein around the mutation site. \u003cb\u003eb3\u003c/b\u003e The overall three-dimensional structure of the Brems1 protein. \u003cb\u003eb4\u003c/b\u003e The local three-dimensional structure of BrEMS1 protein around the mutation site. \u003cb\u003ec\u003c/b\u003e The phylogenetic tree of BrEMS1. \u003cb\u003ed\u003c/b\u003e Amino acid sequence alignment of the S_TKc domain.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u0026nbsp;7\u003c/b\u003e Subcellular localization and qRT-PCR of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026.\u003c/em\u003e \u003cb\u003ea\u003c/b\u003e The subcellular localization of BrEMS1 and Brems1. \u003cb\u003eb\u003c/b\u003e qRT-PCR analysis of \u003cem\u003eBrems1\u003c/em\u003e in different organs of \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026.\u003c/em\u003e\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u0026nbsp;8\u003c/b\u003e Transcriptome analysis of buds in \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5026.\u003c/em\u003e \u003cb\u003ea\u003c/b\u003e Transcriptional levels of DEGs associated with the tapetum, pollen tube and LRR-RLK family. \u003cb\u003eb\u003c/b\u003e KEGG enrichment analysis of DEGs. The red asterisk indicated metabolic pathways.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e\u003c/b\u003e The allelic material \u003cem\u003eM5037\u003c/em\u003e of \u003cem\u003eM5026\u003c/em\u003e. \u003cb\u003ea\u003c/b\u003e \u003cem\u003eM5037\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo; at the flowering stage. Bar, 1 cm. \u003cb\u003eb\u003c/b\u003e The flowers of \u003cem\u003eM5037\u003c/em\u003e and \u0026lsquo;FT\u0026rsquo;. Bar, 3 mm. \u003cb\u003ec\u003c/b\u003e SNP-index distribution plot. The red line represented the mean SNP-index within the window, the pink line denotes the threshold line corresponding to the 99th percentile, and the orange line indicated the threshold line corresponding to the 95th percentile. \u003cb\u003ed\u003c/b\u003e The alignment of cloned gene sequences surrounding the mutation site in \u0026lsquo;FT\u0026rsquo; and \u003cem\u003eM5037.\u003c/em\u003e The red box highlighted the mutation site.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e \u003ch2\u003eConflict of interest\u003c/h2\u003e \u003cp\u003eAll the authors declare that they have no competing interests.\u003c/p\u003e \u003c/p\u003e\u003ch2\u003eAuthors\u0026rsquo; contributions\u003c/h2\u003e \u003cp\u003eCL and LW analyzed the data. CL drafted the manuscript. CL, LW and DZ participated in the creation of materials and performed the experiments. CT and ZL directed the whole study including designing experiments and revising the manuscript.\u003c/p\u003e\u003ch2\u003eAcknowledgements\u003c/h2\u003e \u003cp\u003eThis work was supported by the Earmarked Fund for CARS-23-A02, the National Natural Science Foundation of China (32272736 and 31972406) and the Key Research and Development Program of Liaoning (2021JH2/10200003).\u003c/p\u003e\u003ch2\u003eAvailability of data and materials\u003c/h2\u003e \u003cp\u003eThe datasets supporting the conclusions of this article are included within the article and its additional fles. The Illumina RNA-Seq datasets and Mutmap data are available in the the Sequence Read Archives (SRA) of the National Center for Biotechnology Information NCBI under BioProject ID: PRJNA1105001 andPRJNA1105243. Genomic sequences and gene annotation information of B.rapa are downloaded online at \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://brassicadb.cn\u003c/span\u003e\u003cspan address=\"http://brassicadb.cn\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eAbe A, Kosugi S, Yoshida K, Natsume S, Takagi H, Kanzaki H, Matsumura H, Yoshida K, Mitsuoka C, Tamiru M (2012) Genome sequencing reveals agronomically important loci in rice using MutMap. Nat Biotechnol 30:174\u0026ndash;178\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAn X, Ma B, Duan M, Dong Z, Liu R, Yuan D, Hou Q, Wu S, Zhang D, Liu D, Yu D, Zhang Y, Xie K, Zhu T, Li Z, Zhang S, Tian Y, Liu C, Li J, Yuan L, Wan X (2020) Molecular regulation of ZmMs7 required for maize male fertility and development of a dominant male-sterility system in multiple species. Proc Natl Acad Sci U S A 117:23499\u0026ndash;23509\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBecraft PW (1998) Receptor kinases in plant development. TRENDS PLANT SCI 1998,3(10):384\u0026ndash;388\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen L, Liu Y-G (2014) Male sterility and fertility restoration in crops. Annu Rev Plant Biol 65:579\u0026ndash;606\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDong S, Zou J, Fang B, Zhao Y, Shi F, Song G, Huang S, Feng H (2022) Defect in BrMS1, a PHD-finger transcription factor, induces male sterility in ethyl methane sulfonate-mutagenized Chinese cabbage (Brassica rapa L. ssp. pekinensis). Front Plant Sci 13:992391\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGresshoff PM, Men AE, Laniya TS, Searle IR, Carroll B (2003) Long-Distance Signaling in Nodulation Directed by a CLAVATA1-Like Receptor Kinase. Science 299\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHang D, Yan L, Shunchun L, Lei L, Juanyan W, Zhijie Z, Qichun P, Nan L, Xiangqian Z (2021) Effect of extensin-like OsPEX1 on pollen fertility in rice. Yi Chuan 43:271\u0026ndash;279\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHord CL, Sun YJ, Pillitteri LJ, Torii KU, Wang H, Zhang S, Ma H (2008) Regulation of Arabidopsis early anther development by the mitogen-activated protein kinases, MPK3 and MPK6, and the ERECTA and related receptor-like kinases. Mol Plant 1:645\u0026ndash;658\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHou Q, An X, Ma B, Wu S, Wei X, Yan T, Zhou Y, Zhu T, Xie K, Zhang D, Li Z, Zhao L, Niu C, Long Y, Liu C, Zhao W, Ni F, Li J, Fu D, Yang Z, Wan X (2023) ZmMS1/ZmLBD30-orchestrated transcriptional regulatory networks precisely control pollen exine development. Mol Plant 16:1321\u0026ndash;1338\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuang S, Liu Z, Li D, Yao R, Feng H (2016) A new method for generation and screening of Chinese cabbage mutants using isolated microspore culturing and EMS mutagenesis. Euphytica 207:23\u0026ndash;33\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJi J, Yang L, Fang Z, Zhuang M, Zhang Y, Lv H, Liu Y, Li Z (2017) Recessive male sterility in cabbage (Brassica oleracea var. capitata) caused by loss of function of BoCYP704B1 due to the insertion of a LTR-retrotransposon. Theoretical Appl Genet 130:1441\u0026ndash;1451\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJia G, Liu X, Owen HA, Zhao D (2008) Signaling of cell fate determination by the TPD1 small protein and EMS1 receptor kinase. Proc Natl Acad Sci U S A 105:2220\u0026ndash;2225\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKrzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA (2009) Circos: an information aesthetic for comparative genomics. Genome Res 19:1639\u0026ndash;1645\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKumpf RP, Shi CL, Larrieu A, St\u0026oslash; IM, Butenko MA, P\u0026eacute;ret B, Riiser ES, Bennett MJ, Aalen RB (2013) Floral organ abscission peptide IDA and its HAE/HSL2 receptors control cell separation during lateral root emergence. Proc Natl Acad Sci U S A 110:5235\u0026ndash;5240\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi H, Durbin R (2010) Fast and accurate long-read alignment with Burrows\u0026ndash;Wheeler transform. Bioinformatics 26:589\u0026ndash;595\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiao Y, Smyth GK, Shi W (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30:923\u0026ndash;930\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu C, Song G, Wang N, Huang S, Gao Y, Fu W, Zhang M, Feng H (2021a) A single SNP in Brcer1 results in wax deficiency in Chinese cabbage (Brassica campestris L. ssp. pekinensis). Sci Hort 282:110019\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu T, Jiang GQ, Yao XF, Liu CM (2021b) The leucine-rich repeat receptor-like kinase OsERL plays a critical role in anther lobe formation in rice. Biochem Biophys Res Commun 563:85\u0026ndash;91\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMa Y, Kang J, Wu J, Zhu Y, Wang X (2015) Identification of tapetum-specific genes by comparing global gene expression of four different male sterile lines in Brassica oleracea. Plant Mol Biol 87:541\u0026ndash;554\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMariani C, Beuckeleer MD, Truettner J, Leemans J, Goldberg RB (1990) Induction of male sterility in plants by a chimaeric ribonuclease gene. Nature 347:737\u0026ndash;741\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMasle J, Gilmore SR, Farquhar GD (2005) The ERECTA gene regulates plant transpiration efficiency in Arabidopsis. Nature 436:866\u0026ndash;870\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMccarty DR, Chory J (2000) Conservation and innovation in plant signaling pathways. Cell 103:201\u0026ndash;209\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20:1297\u0026ndash;1303\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5:621\u0026ndash;628\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSanders PM, Bui AQ, Weterings K, McIntire K, Hsu Y-C, Lee PY, Truong MT, Beals T, Goldberg R (1999) Anther developmental defects in Arabidopsis thaliana male-sterile mutants. Sex Plant Reprod 11:297\u0026ndash;322\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSantiago J, Henzler C, Hothorn M (2013) Molecular mechanism for plant steroid receptor activation by somatic embryogenesis co-receptor kinases. Science 341:889\u0026ndash;892\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShi YL, Guo SD, Zhang R, Meng ZG, Ren MZ (2014) The role of Somatic embryogenesis receptor-like kinase 1 in controlling pollen production of the Gossypium anther. Mol Biol Rep 41:411\u0026ndash;422\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShinohara H, Mori A, Yasue N, Sumida K, Matsubayashi Y (2016) Identification of three LRR-RKs involved in perception of root meristem growth factor in Arabidopsis. Proc Natl Acad Sci U S A 113:3897\u0026ndash;3902\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShiu S-H, Karlowski WM, Pan R, Tzeng Y-H, Mayer KF, Li W-H (2004) Comparative analysis of the receptor-like kinase family in Arabidopsis and rice. Plant Cell 16:1220\u0026ndash;1234\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShpak ED, Berthiaume CT, Hill EJ, Torii KU (2004) Synergistic interaction of three ERECTA-family receptor-like kinases controls Arabidopsis organ growth and flower development by promoting cell proliferation. Development 131:1491\u0026ndash;1501\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTabata R, Sumida K, Yoshii T, Ohyama K, Shinohara H, Matsubayashi Y (2014) Perception of root-derived peptides by shoot LRR-RKs mediates systemic N-demand signaling. Science 346:343\u0026ndash;346\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTakeuchi H, Higashiyama T (2016) Tip-localized receptors control pollen tube growth and LURE sensing in Arabidopsis. Nature 531:245\u0026ndash;248\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTan C, Liu Z, Huang S, Feng H (2019) Mapping of the male sterile mutant gene ftms in Brassica rapa L. ssp. pekinensis via BSR-Seq combined with whole-genome resequencing. Theor Appl Genet 132:355\u0026ndash;370\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVizcay-Barrena G, Wilson ZA (2006) Altered tapetal PCD and pollen wall development in the Arabidopsis ms1 mutant. J Exp Bot 57:2709\u0026ndash;2717\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWalker JC (1994) Structure and function of the receptor-like protein kinases of higher plants. Plant Mol Biol 26:1599\u0026ndash;1609\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang K, Li M, Hakonarson H (2010) ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 38:e164\u0026ndash;e164\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWilson ZA, Song J, Taylor B, Yang C (2011) The final split: the regulation of anther dehiscence. J Exp Bot 62:1633\u0026ndash;1649\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWu Hm, Wong E, Ogdahl J, Cheung AY (2000) A pollen tube growth-promoting arabinogalactan protein from Nicotiana alata is similar to the tobacco TTS protein. Plant J 22:165\u0026ndash;176\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXi X, Wei K, Gao B, Liu J, Liang J, Cheng F, Wang X, Wu J (2018) BrFLC5: a weak regulator of flowering time in Brassica rapa. Theoretical Appl Genet 131:2107\u0026ndash;2116\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXu J, Liao R, Xue M, Shang S, Zhou M, Liu Z, Feng H, Huang S (2024) Mutations in BrABCG26, encoding an ATP-binding cassette transporter, are responsible for male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Theor Appl Genet 137:63\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXu SL, Rahman A, Baskin TI, Kieber JJ (2008) Two leucine-rich repeat receptor kinases mediate signaling, linking cell wall biosynthesis and ACC synthase in Arabidopsis. Plant Cell 20:3065\u0026ndash;3079\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYamagishi H, Bhat SR (2014) Cytoplasmic male sterility in Brassicaceae crops. Breed Sci 64:38\u0026ndash;47\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYan L, Yang M, Guo H, Yang L, Wu J, Li R, Liu P, Lian Y, Zheng X, Yan J (2013) Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol biology 20:1131\u0026ndash;1139\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang L, Qian X, Chen M, Fei Q, Meyers BC, Liang W, Zhang D (2016) Regulatory Role of a Receptor-Like Kinase in Specifying Anther Cell Identity. Plant Physiol 171:2085\u0026ndash;2100\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang S-L, Xie L-F, Mao H-Z, Puah CS, Yang W-C, Jiang L, Sundaresan V, Ye D (2003) Tapetum determinant1 is required for cell specialization in the Arabidopsis anther. Plant Cell 15:2792\u0026ndash;2804\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang Z, Sun L, Zhang P, Zhang Y, Yu P, Liu L, Abbas A, Xiang X, Wu W, Zhan X, Cao L, Cheng S (2019) TDR INTERACTING PROTEIN 3, encoding a PHD-finger transcription factor, regulates Ubisch bodies and pollen wall formation in rice. Plant J 99:844\u0026ndash;861\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang J, Wu P, Li N, Xu X, Wang S, Chang S, Zhang Y, Wang X, Liu W, Ma Y, Manghwar H, Zhang X, Min L, Guo X (2022) A male-sterile mutant with necrosis-like dark spots on anthers was generated in cotton. Front Plant Sci 13:1102196\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang W, Sun Y, Timofejeva L, Chen C, Grossniklaus U, Ma H (2006) Regulation of Arabidopsis tapetum development and function by DYSFUNCTIONAL TAPETUM1 (DYT1) encoding a putative bHLH transcription factor. Development 133:3085\u0026ndash;3095\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao D, Wang G, Speal B, Ma H (2002) The excess microsporocytes1 gene encodes a putative leucine-rich repeat receptor protein kinase that controls somatic and reproductive cell fates in the Arabidopsis anther. Genes Dev 16:2021\u0026ndash;2031\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao Y, Huang S, Zou J, Dong S, Wang N, Feng H (2022) Mutation in BrGGL7 gene encoding a GDSL esterase / lipase causes male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Theor Appl Genet 135:3323\u0026ndash;3335\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhou X, Liu Z, Ji R, Feng H (2017) Comparative transcript profiling of fertile and sterile flower buds from multiple-allele-inherited male sterility in Chinese cabbage (Brassica campestris L. ssp. pekinensis). Mol Genet Genomics 292:967\u0026ndash;990\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZorzatto C, Machado J, Lopes K, Nascimento K, Pereira W, Brustolini O, Reis P, Calil I, Deguchi M, Sachetto-Martins G, Gouveia B, Loriato V, Silva M, Silva F, Santos A, Chory J, Fontes E (2015) NIK1-mediated translation suppression functions as a plant antiviral immunity mechanism. Nature 520:679\u0026ndash;682\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZou J, Dong S, Fang B, Zhao Y, Song G, Xin Y, Huang S, Feng H (2023) BrACOS5 mutations induced male sterility via impeding pollen exine formation in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Theor Appl Genet 136:6\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"theoretical-and-applied-genetics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"taag","sideBox":"Learn more about [Theoretical and Applied Genetics](https://www.springer.com/journal/122)","snPcode":"122","submissionUrl":"https://submission.nature.com/new-submission/122/3","title":"Theoretical and Applied Genetics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Chinese cabbage, male sterility, Brems1, transcriptome","lastPublishedDoi":"10.21203/rs.3.rs-4489236/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4489236/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eMale sterile lines are ideal for hybrid seed production in Chinese cabbage. Herein, the complete male sterile mutants M5026\u0026nbsp;and\u0026nbsp;M5073\u0026nbsp;were obtained through ethyl methanesulfonate (EMS) mutagenesis in the cabbage double haploid line ‘FT’. Observation of paraffin sections showed that M5026\u0026nbsp;lacked the tapetum and had excessive microsporocytes. Transmission electron microscopy (TEM) revealed abnormal exine formation in M5026\u0026nbsp;microspores. Genetic analysis revealed a single recessive nuclear gene caused the male sterility phenotype of M5026. Using Mutmap sequencing and Kompetitive allele-specific PCR (KASP) identification and gene cloning, BraA10g029920.3.5C, encoding EMS1 (Excess microsporocytes 1), a leucine-rich repeat receptor-like protein kinase (LRR-RLK), was identified as the candidate gene of M5026\u0026nbsp;and named Brems1. A nonsynonymous G-to-A mutation in an exon of the Brems1\u0026nbsp;gene in M5026\u0026nbsp;resulted in the substitution of glycine with arginine. Employing Mutmap and cloning approaches, a C-to-T SNP was identified within Brems1\u0026nbsp;of M5073, resulting in premature translation termination. Both BrEMS1 and Brems1 were subcellularly localized at the cell membrane. qRT-PCR analysis indicated \u0026nbsp;Brems1\u0026nbsp;exhibited the highest expression level in flower buds, while no expression was detected in roots. Transcriptomic analysis revealed that mutation in Brems1\u0026nbsp;reduced the expression levels of genes associated with the tapetum, pollen tube, and LRR-RLK family. These results suggested that Brems1\u0026nbsp;plays a critical role in pollen development and contributes to elucidating the molecular mechanisms underlying tapetum development and male sterility in Chinese cabbage.\u003c/p\u003e","manuscriptTitle":"Brems1 mutation induced tapetum deficiency leading to male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis)","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-08-05 11:00:19","doi":"10.21203/rs.3.rs-4489236/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revisions","date":"2024-12-19T12:49:01+00:00","index":"","fulltext":""},{"type":"reviewerAgreed","content":"","date":"2024-11-06T12:31:08+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-07-10T16:58:50+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-05-30T06:57:51+00:00","index":"","fulltext":""},{"type":"submitted","content":"Theoretical and Applied Genetics","date":"2024-05-28T04:16:39+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"theoretical-and-applied-genetics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"taag","sideBox":"Learn more about [Theoretical and Applied Genetics](https://www.springer.com/journal/122)","snPcode":"122","submissionUrl":"https://submission.nature.com/new-submission/122/3","title":"Theoretical and Applied Genetics","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"b6257fbf-1f93-46a7-bfff-27e7d56ca92c","owner":[],"postedDate":"August 5th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2025-03-03T16:05:34+00:00","versionOfRecord":{"articleIdentity":"rs-4489236","link":"https://doi.org/10.1007/s00122-025-04841-y","journal":{"identity":"theoretical-and-applied-genetics","isVorOnly":false,"title":"Theoretical and Applied Genetics"},"publishedOn":"2025-02-24 15:57:55","publishedOnDateReadable":"February 24th, 2025"},"versionCreatedAt":"2024-08-05 11:00:19","video":"","vorDoi":"10.1007/s00122-025-04841-y","vorDoiUrl":"https://doi.org/10.1007/s00122-025-04841-y","workflowStages":[]},"version":"v1","identity":"rs-4489236","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4489236","identity":"rs-4489236","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}