- 1State Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology/Innovative Academy of Seed Design, Chinese Academy of Sciences, Beijing, China
- 2University of Chinese Academy of Sciences, Beijing, China
- 3Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, China
- 4Institute of Crop Research, Xinjiang Academy of Agri-Reclamation Sciences, Shihezi, China
Cytoplasmic male sterility (CMS) plays an essential role in hybrid seeds production. In wheat, orf279 was reported as a CMS gene of AL-type male sterile line (AL18A), but its sterility mechanism is still unclear. Therefore, transcriptomic and proteomic analyses of the anthers of AL18A and its maintainer line (AL18B) were performed to interpret the sterility mechanism. Results showed that the electron transport chain and ROS scavenging enzyme expression levels changed in the early stages of the anther development. Biological processes, i.e., fatty acid synthesis, lipid transport, and polysaccharide metabolism, were abnormal, resulting in pollen abortion in AL18A. In addition, we identified several critical regulatory genes related to anther development through combined analysis of transcriptome and proteome. Most of the genes were enzymes or transcription factors, and 63 were partially homologous to the reported genic male sterile (GMS) genes. This study provides a new perspective of the sterility mechanism of AL18A and lays a foundation to study the functional genes of anther development.
Introduction
Heterosis refers to the phenomenon that offspring are superior to the parents in terms of viability, stress resistance, growth potential, and other traits (Chen and Liu, 2014). Heterosis can be utilized to develop new varieties with high yield, good quality, and stress resistance. However, it is challenging to utilize heterosis in self-pollinated crops. Therefore making it necessary to develop male sterile lines (Berlan, 2018). The cytoplasmic male sterile (CMS) line is one of the leading male sterility systems used to produce hybrid seeds. The CMS genes exist in mitochondria while the restorer genes (RF) exist in the nucleus. The CMS/RF system is produced due to abnormal between the nucleus and the cytoplasm (Tang et al., 2017). Different types of CMS lines have been reported in wheat, e.g., T-CMS, K-CMS, V-CMS, AL-CMS, etc. (Gupta et al., 2019). Studies on the male sterility mechanism of wheat can accelerate the utilization of heterosis in wheat.
The molecular mechanisms of cytoplasmic male sterility can be classified into four types: cytotoxicity model, energy deficit model, abnormal PCD model, and retrograde regulation model (Chen and Liu, 2014). In the cytotoxicity model, CMS genes encode toxic proteins, i.e., URF13, ORF138, or ORF79, inhibiting cell growth and development (Korth and Levings, 1993; Duroc et al., 2005; Wang et al., 2006). In the energy deficit model, CMS genes primarily encode transmembrane proteins in the inner membrane of mitochondria and the chimeric part of the mitochondrial genes, affecting the electron transport chain and ATP production leading to energy deficiency and pollen abortion (Sabar et al., 2003; Logan, 2006; Duroc et al., 2009). Timely PCD in the tapetum is crucial to pollen development (Zhang et al., 2016; Yu and Zhang, 2019). In abnormal PCD models, tapetum cells advance or delay PCD, leading to abnormal nutrient absorption of microspores and pollen abortion, like WA352 (Reape and McCabe, 2010; Falasca et al., 2013; Luo et al., 2013). In the retrograde regulation model, CMS genes in mitochondria regulate nuclear gene expression through various retrograde signals, thereby regulating pollen fertility, like orf355 (Xiao et al., 2020). Wan et al. classified the reported genic male sterility genes into four categories: transcriptional regulation, lipid metabolism, polysaccharide metabolism, and others (Wan et al., 2019). The main component of pollen exine, sporopollenin, is synthesized by lipid and phenylpropane metabolism in the tapetum cells (Zhang et al., 2016; Wan et al., 2020).
At present, many genes related to anthers development have been reported in rice and Arabidopsis (Wan et al., 2019). Nevertheless, the mechanism of wheat anther development is still unknown. The advancement in biotechnology, transcriptome, proteome, metabolome, and other omics methods has collectively been utilized in anthers development studies (Zheng et al., 2014; Xing et al., 2018; Cheng et al., 2019; Ning et al., 2019; Zhang et al., 2021), accelerated the research of anther development. The transcriptome analysis of the Aegilops kotschyi thermo-sensitive cytoplasmic male sterile and fertile lines showed that genes related to phenylpropane biosynthesis and jasmonic acid biosynthesis pathway, regulated by MYB transcription factors, were up-regulated in fertile anthers during fertility conversion (Ye et al., 2017). Yang et al. (2018) compared the transcriptome of the sterile line (EP26A) and the maintainer line (EP26) and identified 1716 DEGs (differentially expressed genes) at different anther development stages. GO (Gene Ontology), and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis showed that these genes were mainly involved in redox, carbohydrate, and amino acid metabolism. Proteomic analysis of the anthers of the Yunnan early rice A (ZD-CMS) and its maintainer line revealed that the down-regulation of the 45 DEPs (differentially expressed proteins) involved in the carbohydrate metabolism and stress response in male sterile line. The ROS content in the anther mitochondria of ZD-CMS was increased, and the ultrastructural analysis revealed the destruction of mitochondria (Yan et al., 2014). Proteome and transcriptome analysis of flower buds from the sterile line, Shanan2A and its maintainer Shanan2B, of Brassica napus revealed the involvement of DEGs in carbon metabolism, lipid and flavonoid metabolism, mitochondrial electron transport and ATP synthesis pathways. The DEPs are mainly related to carbohydrate metabolism, energy metabolism, and genetic information processing pathways (Ning et al., 2019). In wheat AL-CMS type sterile line (AL18A), orf279 was reported as the CMS gene (Hao et al., 2021; Melonek et al., 2021), but the sterility mechanism needs further study.
In this study, transcriptomics and proteomics were used to analyze variations between AL18A and AL18B at different stages of anthers development, thus explored the biological processes and metabolic pathways related to the AL18A pollen abortion and identified the genes related to anthers development. This study laid a foundation for the sterility mechanism analysis of wheat AL18A cytoplasmic male sterile line.
Materials and Methods
Plant Material
AL18A, a stable wheat AL-CMS line, was developed from a cross between AL-CMS line 781A and Xindong18 (recurrent parent). The maintainer line Xindong18 will be mentioned as AL18B, hereafter. The plant material was obtained from the Xinjiang Academy of Agri-Reclamation Sciences and planted in the Institute of Genetics and Developmental Biology experimental farm (2018). Anthers of AL18A and AL18B collected at the tetrad, uni-nucleate, bi-nucleate, and tri-nucleate stages were immediately frozen in liquid nitrogen and stored at −80°C, The developmental stage of pollen was determined using the procedure described byBrowne et al. (2018).
Transcriptomic Sequencing
Samples for RNA-seq analysis were prepared from the wheat lines AL18A and AL18B at four different stages of anther development with three biological replications, which made 2×4×3 = 24 libraries in total. Total RNA was extracted using RNEasy Plant Mini Kit (Qiangen, Germany). RNA quality was evaluated on 1% agarose gel, and the purity was checked on a Nano-Photometer spectrophotometer (IMPLEN, CA, United States). RNA concentration was measured using the Qubit RNA Assay Kit and a Qubit 2.0 Fluorometer (Life Technologies, CA, United States), and RNA integrity was assessed using the RNA Nano 6000 Assay Kit on a Bioanalyzer 2100 system (Agilent Technologies, CA, United States). RNA quality testing is done by the cooperative company (BGI, China). The total amount of RNA in each sample was greater than 4 μg, OD260/OD280 ranged from 1.8 to 2.2, OD260/OD230 ranged from 2 to 2.4, and RIN value was greater than 6.3. Ribosomal RNA was depleted using the EpicentreRibo-Zero Gold Kit (Illumina, San Diego), and after purification, depleted rRNA was fragmented into small pieces using fragmentation buffer. The fragmented RNA was reverse transcribed to create the cDNA libraries according to the mRNA-Seq sample preparation kit (Illumina), and paired-end sequencing (2 × 150 bp) was performed on Illumina Hiseq2500 sequencer at BGI company. Trimmomatic (V0.32) (Bolger et al., 2014) was used for data quality filtering, and the bases with Phred quality score mass above 20 were retained. After removing the smaller reads (<25) and one-end reads, the filtered data were aligned to the Chinese Spring wheat reference genome (IWGSC RefSeq V1.0) by HISAT2 (Trapnell et al., 2009). The htseq-count was used to estimate the expression value of each transcript, and the differentially expressed genes (DEGs) were analyzed with DEGseq software (Robinson et al., 2010). P-value < 0.05 and |log2 (fold change) | >1 were used as the threshold values for significantly different expression levels.
DIA Quantification Proteomics
Samples for DIA analysis were prepared from the wheat lines AL18A and AL18B at two different stages (uni-nuclear stage and bi-nuclear stage) of anther development with three biological replications. Each sample was collected with more than 100 mg anther for protein extraction (Imin et al., 2001).The protein samples were sequentially digested with trypsin in a 1:40 ratio for 12 h (4 h + 8 h) at 37°C. Enzymatic peptides were desalted using a Strata X column and vacuumed to dryness. The trypsin digestion was analyzed by LC-MS/MS technique using an Ultimate 3000 nanoLCsystem (Thermo Scientific) coupled to an Orbitrap Fusion Lumos mass spectrometer (Thermo Scientific) operated in either DDA (Data-dependent Acquisition) or DIA (Data-independent Acquisition) acquisition mode at a gradient time of 180 min. Andromeda search engine within MaxQuant (Cox and Mann, 2008) was used for DDA data identification, and the results were used for spectral library construction.
Reference protein database (https://urgi.versailles.inra.fr/download/iwgsc/IWGSC RefSeq_Annotations/v1.0/). For large-scale DIA data, mProphet algorithm was used to complete analytical quality control, thus obtaining a large number of reliable quantitative results. MSstats (Choi et al., 2014) was used to evaluate significant differences in proteins or peptides from different samples. Differentially expressed proteins, having a fold change difference of >2 at the adjusted P-value of <0.05, were considered significant.
Bioinformatics Analysis
Genes identified through the transcriptomics and proteomics were annotated using the following databases: NCBI non-redundant protein sequences (Nr); Kyoto Encyclopedia of Genes and Genomes (KEGG); Gene Ontology (GO); Pfam (a database of protein families); Swiss-Prot (a manually annotated and reviewed protein sequence database); and Eukaryotic orthologous groups (KOGs). Bingo (V2.44) was used for GO enrichment analysis, GO terms with P-value < 0.05 were defined as significantly enriched in DEGs/DEPs, and KEGG functional enrichment analysis of the DEGs/DEPs was executed by KOBAS software. Protein-protein interaction (PPI) was observed using the STRING database (http://string-db.org). Bi-direction clustering of genes and samples was done in the R package.
Results
Alterations in the Electron Transport Chain and ROS Level Can Affect Pollen Fertility of AL18A
More than 12 GB of high-quality RNA-seq data were generated from each sample using the Illumina sequencing platform. The reads were mapped to IWGSCWGA v1.0 with a unique mapped ratio of 69.43–83.32% (Supplementary Table S1). To gain insight into the transcriptome dynamic of anthers development, we performed principal component analysis for the four-time series of AL18A and AL18B (Figure 1A). The plant background and timing of tetrad, uni-nuclear, bi-nuclear, and tri-nuclear stages for anther development are in line. 1504, 916, 2973, and 3889 DEGs between the AL18A and AL18B were obtained at the tetrad, uni-nuclear, bi-nuclear, and the tri-nuclear stages of anther development, respectively (Figure 1B). The number of DEGs gradually increased from the uni-nuclear stage to the tri-nuclear stage, this result is consistent with our previous discovery that AL18A anthers have delayed PCD in the tapetal layer at the uni-nuclear stage (Hao et al., 2021). It is suggested that the early abnormal development of the anther of AL18A was near the uni-nuclear stage. In addition, 153 DEGs were present at all four developmental stages (Supplementary Figure S1; Supplementary Table S2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/762332/fgene-12-762332-HTML/image_m/fgene-12-762332-g001.jpg)
FIGURE 1. Differentially expressed genes and GO function analysis in AL18A at different developmental stages. (A) Principal component analysis (PCA) of RNA-seq data for four stages of wheat anther development. B0, B1, B2, and B3 represent the tetrad, uni-nucleate, bi-nucleate and tri-nucleate stages of anther development of AL18B; A0, A1, A2, and A3 represent the tetrad, uni-nucleate, bi-nucleate and tri-nucleate stages of anther development of AL18A, respectively. (B) The number of DEGs between AL18A and AL18B at different anther developmental stages. DEGs: differentially expressed genes. (C) GO function enrichment analysis of DEGs between AL18A and AL18B at different anther development stages.
To study the molecular mechanism of the pollen abortion in the AL18A line, GO enrichment analysis was performed on the DEGs at four stages of anther development (Figure 1C; Supplementary Table S3). Results showed that the L-phenylalanine metabolic process, fatty acid biosynthetic process, cell wall modification, and the polysaccharide catabolic process were abnormal in the late anther development of AL18A. Fatty acid biosynthetic process and L-phenylalanine metabolic process are known to be involved in pollen wall development (Wan et al., 2019). These DEGs involved in biological processes are consistent with previous reports and provide strong evidence for the pollen abortion of AL18A. In addition, many DEGs participate in electron transport chains and photosynthesis during the early anther development of AL18A. This result is consistent with the genetic background of AL18A and AL18B, because they have the same nucleus but different cytoplasm. Combined with our previous studies that orf279 is the CMS gene of AL18A, the first 96 amino acids of ORF279 and ATP8 were identical (Hao et al., 2021), therefore, orf279 may lead to abortion of AL18A by affecting the electron transport chain.
Most of the CMS genes interfere with the mitochondrial electron transport chain, resulting in energy deficiency or ROS level abnormality leading to advance or delay of tapetal PCD and eventually cause pollen abortion (Chen and Liu, 2014). We found that 31 of the oxidative phosphorylation-related genes were abnormally regulated at the uni-nuclear stage in AL18A anthers (Figure 2A; Supplementary Table S4) and classified into four clusters. Cluster one included 12 NADH dehydrogenase-related genes, cluster two included four cytochrome bc1 complex related genes, cluster three included six cytochrome c complex related genes, and cluster four included nine ATP synthases. In addition, most of the 26 photosynthesis related genes were up-regulated at uni-nuclear stage in AL18A anthers (Figure 2B; Supplementary Table S5). These genes are associated with photosystem I, photosystem II, cytochrome b6f complex, and ATP synthase. We can see that these DEGs are mainly involved in the electron transport chain. Studies have shown that electron transport chain abnormalities can affect ROS levels (Van Aken and Van Breusegem, 2015). We detected that the ROS scavenging enzyme expression changes at different stages of anther development in AL18A and AL18B (Figure 2C; Supplementary Table S6). 32 ROS scavenging enzyme related genes were abnormally regulated in AL18A at anther development stage, and 30 of them encode peroxidases. 7, 8, 20, and 17 ROS scavenging enzymes were differentially expressed at the tetrad, uni-nuclear, bi-nuclear, and tri-nuclear stages. We found that six of the peroxidases were up-regulated, and one was down-regulated in the tetrad stage (Supplementary Table S6). Previous studies have shown that during anther development, the process of PCD in tapetum is regulated by ROS levels (Yu and Zhang, 2019). Differential expression in the ROS scavenging enzymes suggests changes in ROS levels in the AL18A, this is associated with our previous discovery of delayed PCD in the tapetal cells of AL18A (Hao et al., 2021).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/762332/fgene-12-762332-HTML/image_m/fgene-12-762332-g002.jpg)
FIGURE 2. Expression analysis of DEGs involved in electron transport chain and ROS scavenging enzymes. (A) Expression analysis of DEGs in mitochondrial electron transport chains of AL18A and AL18B in uni-nucleate stage. (B) Expression analysis of DEGs in chloroplasts electron transport chains of AL18A and AL18B in uni-nucleate stage. (C) Expression analysis of DEGs in ROS scavenging enzymes of AL18A and AL18B in different anther development stages.
Proteomic Analysis of AL18A and AL18B
To better understand the sterility mechanisms underlying AL18A, we performed the proteomic analysis of AL18A and AL18B. Principal component analysis showed that the proteomic dynamic of anthers development is in line with the plant background and the timing of uni-nuclear and bi-nuclear stages (Figure 3A). The number of DEPs between AL18A and AL18B gradually increased from the uni-nuclear stage (656) to the bi-nuclear stage (1322) (Figure 3B), and these results were consistent with the transcriptome analysis. Eukaryotic orthologous groups (KOGs) analysis was conducted to understand the potential functions of the DEPs globally (Figure 3C). At the uni-nuclear and bi-nuclear stages, most of the DEPs were categorised into post-translational modification and chaperones; translation, ribosomal structure and biogenesis; carbohydrate transport and metabolism; amino acid transport and metabolism; energy production and conversion; lipid transport and metabolism; and signal transduction mechanisms, suggested that the male sterility in AL18A might be related to the disruption of anther cells metabolism at multiple levels.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/762332/fgene-12-762332-HTML/image_m/fgene-12-762332-g003.jpg)
FIGURE 3. Differentially expressed proteins (DEPs) and functional classification of DEPs from anther proteome between AL18A and AL18B. (A) Principal component analysis of proteomic data for two stages of wheat anther development. (B) The number of DEPs between AL18A and AL18B at different anther developmental stages. (C) Functional classification of DEPs at two stages of anther development. All proteins aligned with the KOG database were divided into 24 clusters.
We performed a conjoint analysis in AL18A and AL18B to explore the correlation between the DEGs and DEPs. At the uni-nuclear and bi-nuclear stages, only 1.8% (12/656) and 10.8% (143/1322) of the DEPs overlapped with the DEGs, a total of 146 genes were differentially expressed at both RNA and protein levels (Figures 4A,B; Supplementary Table S7), respectively, indicating a poor correlation between the DEPs and DEGs species accumulation patterns as described in previous studies (Han et al., 2018). These results also infer that proteomics can generate information different from transcriptomics. Therefore, GO analysis was performed on DEPs identified only by the proteome in the uni-nuclear and bi-nuclear stages (Figure 4C; Supplementary Table S8). The result indicated that at uni-nuclear and bi-nuclear stages, the DEPs were involved in fatty acid biosynthesis, lipid transport, aromatic amino acid family metabolism, and DNA catabolism. In addition, polysaccharide catabolic process was also abnormal in uni-nuclear stage of AL18A. Previous studies have shown that abnormal lipid and polysaccharide metabolism often leads to male sterility (Wan et al., 2019); therefore, the DEPs involved in the biological processes leading to AL18A abortion might be essential for anthers fertility. Further analysis of these DEPs would significantly improve our understanding of anthers development.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/762332/fgene-12-762332-HTML/image_m/fgene-12-762332-g004.jpg)
FIGURE 4. Analysis of DEPs detected only by proteomics. (A) Venn diagrams of DEPs and DEGs in uni-nucleate stage. (B) Venn diagrams of DEPs and DEGs in bi-nucleate stage. (C) GO function enrichment analysis of DEPs detected only by proteome in uni-nucleate and bi-nucleate stages.
We further investigated the expression pattern of ROS scavenging enzymes in DEPs. The protein abundance of 30 ROS scavengers between AL18A and AL18B was significantly different (Figure 5; Supplementary Table S9). In AL18A anthers, 12 down-regulated and three up-regulated ROS scavenging enzymes were identified at the uni-nuclear stage; whereas, 18 down-regulated and ten up-regulated enzymes were observed at the bi-nucleated stage. These results were consistent with the transcriptome analysis that the redox homeostasis of AL18A is disrupted, which could explain our earlier discovery that AL18A’s tapetum PCD was delayed (Hao et al., 2021).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/762332/fgene-12-762332-HTML/image_m/fgene-12-762332-g005.jpg)
FIGURE 5. Expression analysis of DEPs in ROS scavenging enzymes of AL18A and AL18B at two anther development stages.
Screening of Important Genes Related to Anthers Development
Many GMS genes have been reported in rice, Arabidopsis, and maize, and most have conserved functions (Wan et al., 2019). However, few wheat GMS genes have been identified and functionally characterized so far. We performed the protein BLAST and the protein interaction analysis of DEPs, and subsequently detected 18 putative GMS orthologs in wheat and their 268 associated protein pairs interacting at the binuclear stage (Figure 6; Supplementary Tables S10, S11).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/762332/fgene-12-762332-HTML/image_m/fgene-12-762332-g006.jpg)
FIGURE 6. Protein interaction networks of 18 putative GMS orthologs genes in wheat. These 18 proteins were DEPs in bi-nucleate stage. The strength of protein interaction is represented by edge color, and the up-down relationship of DEPs is represented by node color. In addition, the putative GMS homologous genes were labeled in the outer part of the figure, and the gene annotations of some of their interacting proteins were labeled in the inner part of the figure.
These results will be useful for the analysis of anther development mechanism. RBOHE in Arabidopsis thaliana encodes NADPH oxidase gene (Xie et al., 2014), TraesCS4B01G327800.2 is a putative homologue of RBOHE in wheat, which prediction interacts with catalase (TraesCS4B01G325800.1). OsGT1 in rice encodes a glycosyltransferase (Moon et al., 2013), TraesCS4A01G446700.1 is a putative homologue of OsGT1 in wheat, which prediction interacts with trehalose-6-phosphate synthase (TraesCS5B01G202200.1, TraesCS3A01G289300.1, TraesCS3D01G289100.1). ABCG26 in Arabidopsis thaliana encodes a member of the ATP-binding cassette (ABC) transporter superfamily (Chang et al., 2016), TraesCS4A01G372100.1 is a putative homologue of ABCG26 in wheat, which prediction interacts with ABC transporter G family member (TraesCS2A01G374100.1, TraesCS3A01G354500.1). CYP704B1 in Arabidopsis thaliana encodes a long-chain fatty acid omega-hydroxylase, is essential for exine development (Dobritsa et al., 2009), TraesCS4B01G341800.1 is a putative homologue of CYP704B1 in wheat, which prediction interacts with 7-dehydrocholesterol reductase (TraesCS3A01G488600.1), methyltransferase (TraesCS6A01G269700.1, TraesCS6B01G297000.1), ferredoxin NADP reductase (TraesCS6D01G012200.1). ACOS5 in Arabidopsis thaliana encodes an enzyme that participates in a conserved and ancient biochemical pathway required for sporopollen in monomer biosynthesis (de Azevedo Souza et al., 2009), TraesCS6A01G151700.1 is a putative homologue of ACOS5 in wheat, which prediction interacts with cytochrome P450 family cinnamate 4-hydroxylase (TraesCS3A01G136100.1, TraesCS3B01G154000.1, TraesCS3A01G342900.1, TraesCS3B01G375100.1), protochlorophyllide reductase (TraesCS1A01G171000.1, TraesCS1D01G168700.1, TraesCS2D01G563600.1), Phenylalanine ammonia-lyase (TraesCS6B01G258400.1), Peroxisomal membrane protein (TraesCS4A01G442900.1). The differential expression of these putative homologous of GMS proteins provided evidence for pollen abortion of AL18A, and these interesting interactions may be related to anthers development and need to be further explored.
Previous studies have shown that biological processes such as redox homeostasis, lipid metabolism and polysaccharide metabolism are involved in tapetum PCD and pollen wall development (Wan et al., 2019). Similar results were obtained in our study. In addition, we found that trehalose biosynthesis process and sexual reproduction are also important for anthers development. We studied six important biological processes involved in anthers development (Figure 7A; Supplementary Tables S12, S13). A total of 46 cell wall organization or biogenesis related DEPs/DEGs, 24 encoding pectinesterase and nine encoding xyloglucan endotransglucosylase/hydrolase (Supplementary Table S14). 72 lipid metabolic process related DEPs/DEGs were detected, 14 encode 3-ketoacyl-CoA synthase, five encode GDSL esterase/lipase, 13 encode lipid transfer protein, seven encode Acyl-desaturase. Ten genes are both detected in DEPs and DEGs, five (TraesCS2A01G279600.1, TraesCS2D01G278500.1, TraesCS4B01G274700.1, TraesCS6A01G104300.1, TraesCS6B01G422900.1) of them are maize orthologs of ms30 (An et al., 2019), and PKSA/LAP6and PKSB/LAP5 in Arabidopsis (Kim et al., 2010). In the process of lipid metabolic process, we identified 15 DEGs that were putative homologue of the reported GMS (Supplementary Table S15). 19 L-phenylalanine metabolic process related DEPs/DEGs were detected, 14 encode phenylalanine ammonialyase, three encode 4-hydroxyphenylpyruvate dioxygenase and two encode glutathione S-transferase (Supplementary Table S16). 38 polysaccharide metabolic process related DEPs/DEGs were detected, two encode Alpha-1,4-glucan-protein synthase, 18 encode beta-amylase, four encode callose synthase. Eight (TraesCS2A01G122600.1, TraesCS3B01G333100.1, TraesCS4A01G166900.1, TraesCS7A01G146200.1, TraesCS7B01G048700.1, TraesCS7B01G243200.1, TraesCS7B01G480300.1, TraesCS7D01G554900.1) of them are rice orthologs of OsUAM3 (Sumiyoshi et al., 2015), and Cals 5 (Dong et al., 2005) and TPD1 (Jia et al., 2008) in Arabidopsis (Supplementary Table S17). In addition, we detected 12 DEPs/DEGs associated with sexual reproduction, 11 encode expansin protein (Supplementary Table S18). Nine DEPs/DEGs associated with trehalose biosynthetic process, seven encode trehalose-6-phosphate synthase, two (TraesCS2D01G168100.1, TraesCS5A01G190000.1) encode trehalose 6-phosphate phosphatase (Supplementary Table S19).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/762332/fgene-12-762332-HTML/image_m/fgene-12-762332-g007.jpg)
FIGURE 7. Screening of genes related to anthers development. (A) Identified of DEPs and DEGs in six important biological processes of anther development. (B) Expression analysis of 122 transcription factors in DEGs.
Anther development is one of the most critical events in the reproductive phase of plant development, involving hundreds of genes and proteins. Transcriptional factors (TFs) are the key connection points in regulatory networks associated with tapetal function and pollen development (Jiang et al., 2021). By comparing AL18A and AL18B, we analyzed the expression pattern of TFs in DEGs during the anther development (Figure 7B; Supplementary Table S20). As a result, we found 8, 7, 62, and 72 differentially expressed TFs at the tetrad, uni-nuclear, bi-nuclear, and tri-nuclear stages, respectively. The number of differentially expressed TFs increased as the anther developed from the uni-nuclear to the tri-nuclear stage, indicating that the anther abortion might have initiated at the uni-nuclear stage. In total, 122 differentially expressed TFs were identified, 40 of them were orthologs of ARF17 (Yang et al., 2013), CDM1 (Lu et al., 2014), DYT1 (Zhang et al., 2006), TDF1 (Zhu et al., 2008), and TGA9 (Murmu et al., 2010) from Arabidopsis, and Ocl4 (Vernoud et al., 2009) from maize (Supplementary Table S20). Although we have identified the GMS genes orthologs in wheat, whether these genes play similar roles in anthers development remains to be investigated.
Discussion
Redox Homeostasis is Important in AL18A Anther Development
Our research showed the abnormal electron transport chain at the uni-nuclear stage and significantly different expression patterns of ROS scavenging enzymes in AL18A compared to the AL18B (Figures 2, 5). Studies have shown that disturbance of the electron transport chain can affect redox homeostasis and the redox homeostasis in tapetum cells can regulate tapetum growth and PCD (Yu and Zhang, 2019). By interacting with COX11, a rice CMS gene WA352 affects the removal of reactive oxygen species in tapetum cells and leads to the premature tapetum PCD, resulting in pollen abortion (Luo et al., 2013). As the CMS gene of AL18A is orf279 and its function was annotated as an ATP synthase subunit 8 gene (Hao et al., 2021; Melonek et al., 2021). We speculate that the CMS gene may affect the homeostasis of reactive oxygen species by interfering with the electron transport chain, affecting the fertility of AL18A. More experiments are needed to confirm this hypothesis in the future.
Genes Related to Anthers Development in Wheat
Studying genes by map-based cloning is time-consuming, and sometimes cloning of homologous genes is fast and effective. By analyzing the putative orthologs of GMS genes in Arabidopsis and rice, a putative regulatory network for anther and pollen development in maize was constructed (Wan et al., 2019; Jiang et al., 2021). In our study, DEGs/DEPs are involved in fatty acid biosynthesis process, lipid transport, L-phenylalanine metabolic process and polysaccharide metabolism in AL18A’s anthers development (Figure 7). The results are similar to previous studies in rice, Arabidopsis thaliana and maize (Wan et al., 2019). Meanwhile, we identified a number of GMS genes that are partially homologous to other species in wheat. The reverse genetics approach can be used to study these partially homologous GMS genes which are found in wheat.
A Possible Working Model Related to Male Sterility in AL18A
Based on previous reports and this study, we constructed a proposed working model of impaired anther development in the AL18A (Figure 8). The orf279 of AL18A may changes ATP synthesis and ROS levels by affecting the mitochondrial electron transport chain of tapetum cells, which leads to delayed tapetum PCD (Hao et al., 2021). The aberrant tapetum PCD causes abnormality in cell wall organization, lipid metabolic process, L-phenylalanine metabolic process, polysaccharide metabolic process, trehalose biosynthetic process, and sexual reproduction. The aberration in the biological process causes pollen abortion in AL18A.
Data Availability Statement
The data presented in the study are deposited in the NCBI repository, accession number SAMN19979141–SAMN19979164.
Author Contributions
WY and AZ conceived and designed the study. JS, DL, XL, YN, and XT prepared the materials. MH and TL performed experiments and wrote the initial manuscript. WY and MH performed data analysis. MS, WY, and AZ review and editing the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Key R and D Program of China (grant number, 2016YFD0100804).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.762332/full#supplementary-material
References
An, X., Dong, Z., Tian, Y., Xie, K., Wu, S., Zhu, T., et al. (2019). ZmMs30 Encoding a Novel GDSL Lipase Is Essential for Male Fertility and Valuable for Hybrid Breeding in maize. Mol. Plant 12 (3), 343–359. doi:10.1016/j.molp.2019.01.011
Berlan, J.-P. (2018). Hybrid Corn and the Unsettled Question of Heterosis. J. Genet. 97 (5), 1075–1082. doi:10.1007/s12041-018-1037-2
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a Flexible Trimmer for Illumina Sequence Data. Bioinformatics 30 (15), 2114–2120. doi:10.1093/bioinformatics/btu170
Browne, R. G., Iacuone, S., Li, S. F., Dolferus, R., and Parish, R. W. (2018). Anther Morphological Development and Stage Determination in Triticum aestivum. Front. Plant Sci. 9, 228. doi:10.3389/fpls.2018.00228
Chang, Z., Chen, Z., Yan, W., Xie, G., Lu, J., Wang, N., et al. (2016). An ABC Transporter, OsABCG26, Is Required for Anther Cuticle and Pollen Exine Formation and Pollen-Pistil Interactions in rice. Plant Sci. 253, 21–30. doi:10.1016/j.plantsci.2016.09.006
Chen, L., and Liu, Y.-G. (2014). Male Sterility and Fertility Restoration in Crops. Annu. Rev. Plant Biol. 65, 579–606. doi:10.1146/annurev-arplant-050213-040119
Cheng, Q., Li, T., Ai, Y., Lu, Q., Wang, Y., Sun, L., et al. (2019). Complementary Transcriptomic and Proteomic Analysis Reveals a Complex Network Regulating Pollen Abortion in GMS (Msc-1) Pepper (Capsicum Annuum L). Int. J. Mol. Sci. 20 (7), 1789. doi:10.3390/ijms20071789
Choi, M., Chang, C.-Y., Clough, T., Broudy, D., Killeen, T., MacLean, B., et al. (2014). MSstats: an R Package for Statistical Analysis of Quantitative Mass Spectrometry-Based Proteomic Experiments. Bioinformatics 30 (17), 2524–2526. doi:10.1093/bioinformatics/btu305
Cox, J., and Mann, M. (2008). MaxQuant Enables High Peptide Identification Rates, Individualized p.p.b.-range Mass Accuracies and Proteome-wide Protein Quantification. Nat. Biotechnol. 26 (12), 1367–1372. doi:10.1038/nbt.1511
de Azevedo Souza, C., Kim, S. S., Koch, S., Kienow, L., Schneider, K., McKim, S. M., et al. (2009). A Novel Fatty Acyl-CoA Synthetase Is Required for Pollen Development and Sporopollenin Biosynthesis inArabidopsis. Plant Cell 21 (2), 507–525. doi:10.1105/tpc.108.062513
Dobritsa, A. A., Shrestha, J., Morant, M., Pinot, F., Matsuno, M., Swanson, R., et al. (2009). CYP704B1 Is a Long-Chain Fatty Acidω-Hydroxylase Essential for Sporopollenin Synthesis in Pollen of Arabidopsis. Arabidopsis. Plant Physiol. 151 (2), 574–589. doi:10.1104/pp.109.144469
Dong, X., Hong, Z., Sivaramakrishnan, M., Mahfouz, M., and Verma, D. P. S. (2005). Callose Synthase (CalS5) Is Required for Exine Formation during Microgametogenesis and for Pollen Viability in Arabidopsis. Plant J. 42 (3), 315–328. doi:10.1111/j.1365-313x.2005.02379.x
Duroc, Y., Hiard, S., Vrielynck, N., Ragu, S., and Budar, F. (2009). The Ogura Sterility-Inducing Protein Forms a Large Complex without Interfering with the Oxidative Phosphorylation Components in Rapeseed Mitochondria. Plant Mol. Biol. 70 (1-2), 123–137. doi:10.1007/s11103-009-9461-6
Duroc, Y., Gaillard, C., Hiard, S., Defrance, M.-C., Pelletier, G., and Budar, F. (2005). Biochemical and Functional Characterization of ORF138, a Mitochondrial Protein Responsible for Ogura Cytoplasmic Male Sterility in Brassiceae. Biochimie 87 (12), 1089–1100. doi:10.1016/j.biochi.2005.05.009
Falasca, G., D'Angeli, S., Biasi, R., Fattorini, L., Matteucci, M., Canini, A., et al. (2013). Tapetum and Middle Layer Control Male Fertility in Actinidia Deliciosa. Ann. Bot. 112 (6), 1045–1055. doi:10.1093/aob/mct173
Gupta, P. K., Balyan, H. S., Gahlaut, V., Saripalli, G., Pal, B., Basnet, B. R., et al. (2019). Hybrid Wheat: Past, Present and Future. Theor. Appl. Genet. 132 (9), 2463–2483. doi:10.1007/s00122-019-03397-y
Han, F., Zhang, X., Yang, L., Zhuang, M., Zhang, Y., Li, Z., et al. (2018). iTRAQ-Based Proteomic Analysis of Ogura-CMS Cabbage and its Maintainer Line. Int. J. Mol. Sci. 19 (10), 3180. doi:10.3390/ijms19103180
Hao, M., Yang, W., Lu, W., Sun, L., Shoaib, M., Sun, J., et al. (2021). Characterization of the Mitochondrial Genome of a Wheat AL-type Male Sterility Line and the Candidate CMS Gene. Int. J. Mol. Sci. 22 (12), 6388. doi:10.3390/ijms22126388
Imin, N., Kerim, T., Weinman, J. J., and Rolfe, B. G. (2001). Characterisation of rice Anther Proteins Expressed at the Young Microspore Stage. Proteomics 1 (9), 1149–1161. doi:10.1002/1615-9861(200109)1:9<1149::AID-PROT1149>3.0.CO;2-R
Jia, G., Liu, X., Owen, H. A., and Zhao, D. (2008). Signaling of Cell Fate Determination by the TPD1 Small Protein and EMS1 Receptor Kinase. Proc. Natl. Acad. Sci. 105 (6), 2220–2225. doi:10.1073/pnas.0708795105
Jiang, Y., An, X., Li, Z., Yan, T., Zhu, T., Xie, K., et al. (2021). CRISPR/Cas9‐based Discovery of maize Transcription Factors Regulating Male Sterility and Their Functional Conservation in Plants. Plant Biotechnol. J. 19, 1769–1784. doi:10.1111/pbi.13590
Kim, S. S., Grienenberger, E., Lallemand, B., Colpitts, C. C., Kim, S. Y., Souza, C. d. A., et al. (2010). LAP6/Polyketide Synthase AandLAP5/Polyketide Synthase Bencode Hydroxyalkyl α-Pyrone Synthases Required for Pollen Development and Sporopollenin Biosynthesis inArabidopsis Thaliana. Plant Cell 22 (12), 4045–4066. doi:10.1105/tpc.110.080028
Korth, K. L., and Levings, C. S. (1993). Baculovirus Expression of the maize Mitochondrial Protein URF13 Confers Insecticidal Activity in Cell Cultures and Larvae. Proc. Natl. Acad. Sci. 90 (8), 3388–3392. doi:10.1073/pnas.90.8.3388
Logan, D. C. (2006). The Mitochondrial Compartment. J. Exp. Bot. 57 (6), 1225–1243. doi:10.1093/jxb/erj151
Lu, P., Chai, M., Yang, J., Ning, G., Wang, G., and Ma, H. (2014). The Arabidopsis Callose Defective Microsporel Gene Is Required for Male Fertility through Regulating Callose Metabolism during Microsporogenesis. Plant Physiol. 164 (4), 1893–1904. doi:10.1104/pp.113.233387
Luo, D., Xu, H., Liu, Z., Guo, J., Li, H., Chen, L., et al. (2013). A Detrimental Mitochondrial-Nuclear Interaction Causes Cytoplasmic Male Sterility in rice. Nat. Genet. 45 (5), 573–577. doi:10.1038/ng.2570
Melonek, J., Duarte, J., Martin, J., Beuf, L., Murigneux, A., Varenne, P., et al. (2021). The Genetic Basis of Cytoplasmic Male Sterility and Fertility Restoration in Wheat. Nat. Commun. 12 (1), 1036. doi:10.1038/s41467-021-21225-0
Moon, S., Kim, S.-R., Zhao, G., Yi, J., Yoo, Y., Jin, P., et al. (2013). Rice GLYCOSYLTRANSFERASE1 Encodes a Glycosyltransferase Essential for Pollen Wall Formation. Plant Physiol. 161 (2), 663–675. doi:10.1104/pp.112.210948
Murmu, J., Bush, M. J., DeLong, C., Li, S., Xu, M., Khan, M., et al. (2010). Arabidopsis Basic Leucine-Zipper Transcription Factors TGA9 and TGA10 Interact with floral Glutaredoxins ROXY1 and ROXY2 and Are Redundantly Required for Anther Development. Plant Physiol. 154 (3), 1492–1504. doi:10.1104/pp.110.159111
Ning, L., Wang, H., Li, D., Lin, Z., Li, Y., Zhao, W., et al. (2019). Transcriptomic and Proteomic Analysis of Shaan2A Cytoplasmic Male Sterility and its Maintainer Line in Brassica Napus. Front. Plant Sci. 10, 252. doi:10.3389/fpls.2019.00252
Reape, T. J., and McCabe, P. F. (2010). Apoptotic-like Regulation of Programmed Cell Death in Plants. Apoptosis 15 (3), 249–256. doi:10.1007/s10495-009-0447-2
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26 (1), 139–140. doi:10.1093/bioinformatics/btp616
Sabar, M., Gagliardi, D., Balk, J., and Leaver, C. J. (2003). ORFB Is a Subunit of F 1 F O ‐ATP Synthase: Insight into the Basis of Cytoplasmic Male Sterility in sunflower. EMBO Rep. 4 (4), 381–386. doi:10.1038/sj.embor.embor800
Sumiyoshi, M., Inamura, T., Nakamura, A., Aohara, T., Ishii, T., Satoh, S., et al. (2015). UDP-arabinopyranose Mutase 3 Is Required for Pollen wall Morphogenesis in rice (Oryza Sativa). Plant Cel Physiol. 56 (2), 232–241. doi:10.1093/pcp/pcu132
Tang, H., Xie, Y., Liu, Y.-G., and Chen, L. (2017). Advances in Understanding the Molecular Mechanisms of Cytoplasmic Male Sterility and Restoration in rice. Plant Reprod. 30 (4), 179–184. doi:10.1007/s00497-017-0308-z
Trapnell, C., Pachter, L., and Salzberg, S. L. (2009). TopHat: Discovering Splice Junctions with RNA-Seq. Bioinformatics 25 (9), 1105–1111. doi:10.1093/bioinformatics/btp120
Van Aken, O., and Van Breusegem, F. (2015). Licensed to Kill: Mitochondria, Chloroplasts, and Cell Death. Trends Plant Sci. 20 (11), 754–766. doi:10.1016/j.tplants.2015.08.002
Vernoud, V., Laigle, G., Rozier, F., Meeley, R. B., Perez, P., and Rogowsky, P. M. (2009). The HD-ZIP IV Transcription Factor OCL4 Is Necessary for Trichome Patterning and Anther Development in maize. Plant J. 59 (6), 883–894. doi:10.1111/j.1365-313x.2009.03916.x
Wan, X., Wu, S., Li, Z., An, X., and Tian, Y. (2020). Lipid Metabolism: Critical Roles in Male Fertility and Other Aspects of Reproductive Development in Plants. Mol. Plant 13 (7), 955–983. doi:10.1016/j.molp.2020.05.009
Wan, X., Wu, S., Li, Z., Dong, Z., An, X., Ma, B., et al. (2019). Maize Genic Male-Sterility Genes and Their Applications in Hybrid Breeding: Progress and Perspectives. Mol. Plant 12 (3), 321–342. doi:10.1016/j.molp.2019.01.014
Wang, Z., Zou, Y., Li, X., Zhang, Q., Chen, L., Wu, H., et al. (2006). Cytoplasmic Male Sterility of rice with Boro II Cytoplasm Is Caused by a Cytotoxic Peptide and Is Restored by Two Related PPR Motif Genes via Distinct Modes of mRNA Silencing. Plant Cell 18 (3), 676–687. doi:10.1105/tpc.105.038240
Xiao, S., Zang, J., Pei, Y., Liu, J., Liu, J., Song, W., et al. (2020). Activation of Mitochondrial Orf355 Gene Expression by a Nuclear-Encoded DREB Transcription Factor Causes Cytoplasmic Male Sterility in Maize. Mol. Plant 13 (9), 1270–1283. doi:10.1016/j.molp.2020.07.002
Xie, H.-T., Wan, Z.-Y., Li, S., and Zhang, Y. (2014). Spatiotemporal Production of Reactive Oxygen Species by NADPH Oxidase Is Critical for Tapetal Programmed Cell Death and Pollen Development in Arabidopsis. Plant Cell 26 (5), 2007–2023. doi:10.1105/tpc.114.125427
Xing, M., Sun, C., Li, H., Hu, S., Lei, L., and Kang, J. (2018). Integrated Analysis of Transcriptome and Proteome Changes Related to the Ogura Cytoplasmic Male Sterility in Cabbage. PLoS One 13 (3), e0193462. doi:10.1371/journal.pone.0193462
Yan, J., Tian, H., Wang, S., Shao, J., Zheng, Y., Zhang, H., et al. (2014). Pollen Developmental Defects in ZD-CMS rice Line Explored by Cytological, Molecular and Proteomic Approaches. J. Proteomics 108, 110–123. doi:10.1016/j.jprot.2014.05.014
Yang, J., Tian, L., Sun, M.-X., Huang, X.-Y., Zhu, J., Guan, Y.-F., et al. (2013). Auxin Response Factor17 Is Essential for Pollen Wall Pattern Formation in Arabidopsis. Arabidopsis. Plant Physiol. 162 (2), 720–731. doi:10.1104/pp.113.214940
Yang, Y., Bao, S., Zhou, X., Liu, J., and Zhuang, Y. (2018). The Key Genes and Pathways Related to Male Sterility of Eggplant Revealed by Comparative Transcriptome Analysis. BMC Plant Biol. 18 (1), 209. doi:10.1186/s12870-018-1430-2
Ye, J., Duan, Y., Hu, G., Geng, X., Zhang, G., Yan, P., et al. (2017). Identification of Candidate Genes and Biosynthesis Pathways Related to Fertility Conversion by Wheat KTM3315A Transcriptome Profiling. Front. Plant Sci. 8, 449. doi:10.3389/fpls.2017.00449
Yu, J., and Zhang, D. (2019). Molecular Control of Redox Homoeostasis in Specifying the Cell Identity of Tapetal and Microsporocyte Cells in rice. Rice (N Y) 12, 42. doi:10.1186/s12284-019-0300-3
Zhang, D., Shi, J., and Yang, X. (2016). Role of Lipid Metabolism in Plant Pollen Exine Development. Subcell Biochem. 86, 315–337. doi:10.1007/978-3-319-25979-6_13
Zhang, W., Sun, Y., Timofejeva, L., Chen, C., Grossniklaus, U., and Ma, H. (2006). Regulation of Arabidopsis Tapetum Development and Function by DYSFUNCTIONAL TAPETUM1 (DYT1) Encoding a Putative bHLH Transcription Factor. Development 133 (16), 3085–3095. doi:10.1242/dev.02463
Zhang, Y., Song, Q., Zhang, L., Li, Z., Wang, C., and Zhang, G. (2021). Comparative Proteomic Analysis of Developmental Changes in P-type Cytoplasmic Male Sterile and Maintainer Anthers in Wheat. Int. J. Mol. Sci. 22 (4), 2012. doi:10.3390/ijms22042012
Zheng, B.-B., Fang, Y.-N., Pan, Z.-Y., Sun, L., Deng, X.-X., Grosser, J. W., et al. (2014). iTRAQ-Based Quantitative Proteomics Analysis Revealed Alterations of Carbohydrate Metabolism Pathways and Mitochondrial Proteins in a Male Sterile Cybrid Pummelo. J. Proteome Res. 13 (6), 2998–3015. doi:10.1021/pr500126g
Keywords: AL-type cytoplasmic male sterility, AL18A, transport, anthers development genes, wheat
Citation: Hao M, Yang W, Li T, Shoaib M, Sun J, Liu D, Li X, Nie Y, Tian X and Zhang A (2021) Combined Transcriptome and Proteome Analysis of Anthers of AL-type Cytoplasmic Male Sterile Line and Its Maintainer Line Reveals New Insights into Mechanism of Male Sterility in Common Wheat. Front. Genet. 12:762332. doi: 10.3389/fgene.2021.762332
Received: 28 August 2021; Accepted: 03 December 2021;
Published: 17 December 2021.
Edited by:
Karthikeyan Adhimoolam, Jeju National University, South KoreaReviewed by:
Rajeev Ranjan, Purdue University, United StatesPallavi Sinha, IRRI South Asia Hub, India
Copyright © 2021 Hao, Yang, Li, Shoaib, Sun, Liu, Li, Nie, Tian and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenlong Yang, wlyang@genetics.ac.cn; Aimin Zhang, amzhang@genetics.ac.cn
†These authors have contributed equally to this work