Next Article in Journal
Anti-Diabetic Activity of Bioactive Compound Extracted from Spondias mangifera Fruit: In-Vitro and Molecular Docking Approaches
Previous Article in Journal
Examining the Interaction between Phytophthora sojae and Soybean Cyst Nematode on Soybean (Glycine max)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electrophysiological, Morphologic, and Transcriptomic Profiling of the Ogura-CMS, DGMS and Maintainer Broccoli Lines

1
Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Ministry of Agriculture, #12 Zhong Guan Cun Nandajie Street, Beijing 100081, China
2
China Vegetable Biotechnology (Shouguang) Co., Ltd., Shouguang 262700, China
3
Institute of Agricultural Resources and Environment, Jiangsu Academy of Agricultural Sciences, Nanjing 210014, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 25 January 2022 / Revised: 3 February 2022 / Accepted: 16 February 2022 / Published: 21 February 2022
(This article belongs to the Topic Plant Breeding, Genetics and Genomics)

Abstract

:
To better serve breeding of broccoli, the electrophysiological, morphological and transcriptomic profiling of the isogenic Ogura-CMS, DGMS and their maintainer fertile lines, were carried out by scanning electron microscopy, investigation of agronomic traits and RNA-sequencing analysis. The agronomic traits of plant height, length of the largest leaf, plant spread angle, single head weight, head width and stem diameter showed stronger performance in Ogura-CMS broccoli than in DGMS line or maintainer fertile line. However, the Ogura-CMS broccoli was poorer in the seed yield and seed germination than in the DGMS line and maintainer fertile line. Additionally, the DGMS broccoli had longer maturation and flowering periods than the Ogura-CMS and maintainer fertile lines. There were obvious differences in the honey gland, happening in the male sterility and fertile lines of broccoli. Additionally, the mechanism regulating Ogura-CMS and DGMS in broccoli was investigated using florets transcriptome analyses of the Ogura-CMS, DGMS and maintainer fertile lines. As a result, a total of 2670 differentially expressed genes (DEGs) were detected, including 1054 up- and 1616 downregulated genes in the Ogura-CMS and DGMS lines compared to the maintainer fertile line. A number of functionally known genes involved in plant hormones (auxin, salicylic acid and brassinosteroid), five Mitochondrial Oxidative Phosphorylation (OXPHOS) genes of atp8, LOC106319879, LOC106324734, LOC106314622 and LOC106298585, and three upregulated genes (Lhcb1, Lhcb3 and Lhcb5) associated with the photosynthesis-antenna protein pathway, were obviously detected to be highly associated with reproductive development including flowering time, maturity and reproductive period in the Ogura-CMS and DGMS broccoli comparing to their maintainer fertile line. Our research would provide a comprehensive foundation for understanding the differences of electrophysiological, morphological and transcriptomic profiles in the Ogura-CMS, DGMS and maintainer broccoli, and as well as being beneficial to exploring the mechanism of male sterility in Brassica crops.

1. Introduction

Broccoli (Brassica oleracea L. var. italica) is a well-known vegetable grown worldwide that belongs to the cruciferous vegetable family, which includes cabbage, kale, cauliflower, Brussels sprouts, Chinese cabbage, rutabaga and turnips. Broccoli contains many nutrients and few calories [1], exhibiting anticancer [2], weight loss promoting [3], cardiovascular and cerebrovascular disease-preventing properties that are beneficial to human health [4,5]. In 2017, the area in China under cultivation for broccoli has exceeded 70,000 ha based on data from the Ministry of Agriculture [6,7].
Currently, commercial broccoli crops primarily consist of hybrid cultivars benefiting from heterosis, and the hybrid cultivars are mostly male sterility lines, especially Ogura cytoplasmic male sterility (Ogura-CMS) lines. The male sterility lines currently used in broccoli breeding were derived from a line with dominant genic male sterility (DGMS), which was discovered in cabbage, and from a line Ogura-CMS, which was discovered in radish [8,9,10]. However, there are still some disadvantages regarding bud development and seed production in male sterility lines, especially in Ogura-CMS plants; these disadvantages include bud abortion and low seed yield [8,11,12]. On the other hand, hybrid cultivars from male sterility lines have enhanced seed purity and hybrid rates [11].
Transcriptome analyses based on RNA-sequencing (RNA-Seq) technology are important methods with which to understand the differentially expressed genes (DEGs) involved in various biological processes, and these analyses have been widely performed on Brassica plants [12,13,14]. For example, a small RNA-Seq has revealed DEGs in the early development of broccoli pollen [15], and transcriptome analyses have been used to examine four transcription factors in normal and abortive buds of an Ogura-CMS line and its maintainer line [12]. In addition, the roles of pigment mechanisms in post-harvest broccoli yellowing [16], gene expression patterns associated with sulforaphane metabolism in broccoli florets [17], glucosinolate metabolism in seeds and sprouts [18] and yield heterosis in curds of broccoli have also been elucidated by RNA-Seq technology [19]. So far, few studies on the mechanism of male sterility in broccoli have been reported. To the best of our knowledge, it would be the first report of comparative transcriptome analyses of isogenic Ogura-CMS, DGMS and inbred lines in broccoli florets.
To obtain a new male sterile line of broccoli for cultivar breeding, we used multiple broccoli inbred lines with cabbage DGMS material and developed several good DGMS broccoli lines with stable agronomic performance [8]. In the same way, Ogura-CMS broccoli lines derived from CMSR3629 (Ogura-CMS) introduced by the Asgrow Seed Co. (USA) have also been developed. Among these materials, the inbred line T54 (T54S) has selfed more than 10 generations. T54S was used to develop an isogenic Ogura-CMS line (T54C) and an isogenic DGMS line (T54M).

2. Results

2.1. Investigation of Agronomic Traits and the Nectary Morphology

There were no significant differences in head characteristics among the Ogura-CMS line T54C, the DGMS line T54M and their maintainer fertile line T54S, but there were some other differences in fertility, maturity days, flowering time, production and seed characters, which are shown in and Figure 1 and Table 1. Agronomic traits of broccoli at harvest revealed that the Ogura-CMS broccoli showed stronger performance in plant height, leaf length, plant spread angle, single head weight, head width and stem diameter than the DGMS line or the maintainer fertile line. The DGMS line had longer production (days to maturation and flowering) than the Ogura-CMS and maintainer lines. In contrast, the maintainer fertile line presented a higher seed yield and seed germination rate than the DGMS and Ogura-CMS lines. The Ogura-CMS line showed the lowest seed yield and germination rate. From Figure 1, we could find more honey glands in the maintainer line than in the Ogura-CMS and DGMS lines, and most of them were open, while honey glands in the DGMS line was less significant than the other both lines.

2.2. Genomic Characteristicsand SNP Distribution

Six cDNA libraries from broccoli florets obtained at harvest times from the T54C, T54M and T54S lines were subjected to Illumina sequencing with two biological replicates for each sample. After the filtering of invalid reads, 159,424,568 clean reads and 47,827,370,400 clean bases (47.8 Gb) were obtained (Table S1), while the Q20 and Q30 percentages were 98.12–98.41% and 95.06–95.62%, respectively, with GC percentages of 45.08–46.80%. From Table S2, we found that eight regions of the genome contained Single-Nucleotide Polymorphism (SNP). The T54S line showed a higher percentage of SNP distribution in exonic regions (54.39%) than the T54C (52.35%) and T54M (52.79%) lines but a lower percentage in intergenic regions (15.85%) than the T54C (16.96%) and T54M (17.25%) lines. The remaining regions of SNP distribution showed no clear differences.

2.3. Gene Optimization and Correlation Test for Each Sample

A total of 5754 optimized regions involving 3872 unigenes were predicted using the gene structure optimization method, and approximately 13 mitochondrial unigenes (chrMT) were optimized. To study the reliability of the results and differences among the samples, correlation tests were performed for each sample comparison based on the log2 (FPKM) values of gene expression (Table 2 and Figure S1). All the coefficients of correlation for the pairs of samples ranged from 0.913 to 0.973, with the Pearson’s R2 test values ranging from 0.867 to 0.973 (R2 > 0.8). The coefficients of correlation among replicate T54C, T54M and T54S samples were 0.965 (R2 = 0.943), 0.960 (R2 = 0.925) and 0.966 (R2 = 0.919), respectively. The coefficients of correlation between T54C and T54M, T54C and T54S, and T54M and T54S ranged from 0.867 to 0.918, 0.867 to 0.905, and 0.894 to 0.973, respectively. Therefore, the DGMS and maintainer lines were more closely related than the Ogura-CMS and maintainer lines.

2.4. DEGs Analysis

To explore the DEGs, the gene expression variations were analyzed in three comparisons T54C vs. T54M, T54C vs. T54S and T54M vs. T54S (Figure 2). In the Ogura-CMS line compared with the maintainer line, 585 and 1073 genes were up- and downregulated, respectively. In the DGMS line compared with the maintainer line, 469 and 534 genes were up- and downregulated, respectively. Moreover, we also compared the Ogura-CMS line to the DGMS line, 505 and 1109 up- and downregulated genes were identified, respectively. The DEG cluster analysis (Figure 2) revealed that the Ogura-CMS line showed similar up- and downregulated genes in comparison with the DGMS and maintainer lines, but the T54M vs. T54S comparison presented very different numbers and types.
The KEGG pathway results were comprehensively evaluated and filtered (p < 0.01), and there were eight pathways strongly related to plant development of the Ogura-CMS line, DGMS and maintainer lines. Specifically, 7 upregulated genes related with photosynthesis-antenna proteins, 41 upregulated genes related with plant hormone signal transduction, 39 downregulated genes related with phenylpropanoid biosynthesis, 4 up- and 5 downregulated genes related with adenosine triphosphate (ATP)-binding cassette (ABC) transporters, 5 upregulated genes related with fatty acid elongation, 25 downregulated genes related with fatty acid metabolism, 15 downregulated genes related with flavonoid biosynthesis and 8 upregulated genes related with polycyclic aromatic hydrocarbon degradation pathways were distributed among different samples, which provided potential genes in regulating morphological characteristics of the Ogura-CMS line, DGMS and inbred lines (Table S3).

2.5. Genes Related to Plant Hormones

Plant hormones play important roles throughout the life span of a plant, such as auxin, gibberellin (GA), cytokinin (CK), abscisic acid (ABA), ethylene (ET), jasmonic acid (JA), salicylic acid (SA) and brassinosteroid (BR). The individual pathways of various plant hormone responses have been widely studied in the model organism Arabidopsis thaliana [20]. In this study, comparing to the maintainer line, in the plant hormone signal transduction pathway, 35 and 15 genes were upregulated in the Ogura-CMS and DGMS lines, respectively. Additionally, 49 and 17 genes were downregulated in the Ogura-CMS and DGMS lines, respectively (Figure 3). In the Ogura-CMS line T54C, these genes were associated with auxin, CK, ABA, ET, JA and SA metabolism. In the DGMS line T54M, they were obviously auxin genes (Aux/IAA and SAUR), CK-related gene (A-ARR), ABA-related gene (PP2C), ET-related gene (ERF1/2), BR-related gene (CYCD3) and JA-related gene (JAZ). Compared to the Ogura-CMS line T54C, the DGMS line T54M, there were common genes including auxin (Aux/IAA and SAUR), CK (A-ARR), ABA (PP2C), ET (ERF1/2), BR (CYCD3) and JA (JAZ). However, some others also exhibited in the DGMS line, with regard to auxin metabolism, two genes, TIR1 and GH3, were absent, and another similar gene, Aux/IAA, was present, which played a crucial role in repressing the expression levels of genes activated by auxin response factors (ARFs) in Arabidopsis [21].

2.6. Genes Related to Oxidative Phosphorylation and Photosynthesis

Mitochondrial Oxidative Phosphorylation (OXPHOS) provides Adenosine triphosphate (ATP) for driving cellular functions, which is of central importance for almost all eukaryotic cells. In plants, OXPHOS takes place in the condition of photosynthesis. Indeed, the metabolism of mitochondria and chloroplasts is tightly linked, and photosynthesis deeply affects OXPHOS [22]. In this study, 5 OXPHOS-upregulated genes of atp8, LOC106319879, LOC106324734, LOC106314622 and LOC106298585 were obviously found in the Ogura-CMS and DGMS lines, especially the atp8 gene just detected in the Ogura-CMS line with a higher relative expression level but nearly no expression in the maintainer line (Figure 4). From Figure 4, we could find the functions of these atp-related genes based on ko00190 (Oxidative phosphorylation) and ko00195 (Photosynthesis). At the same time, six genes related with photosynthesis-antenna proteins were upregulated in the Ogura-CMS line comparing to its maintainer line, while seven photosynthesis-antenna proteins and eight photosynthesis-related genes were upregulated in the DGMS line (Figure 4). According to our evaluation, three upregulated genes (Lhcb1, Lhcb3 and Lhcb5) associated with the photosynthesis-antenna protein pathway (ko00196) in the DGMS line were emphasized in this study.

2.7. Real-Time qPCR Validation of Gene Expression Patterns

To better understand and explain the different agronomic traits of the broccoli lines T54S, T54M and T54C, the number of 14 DEGs were selected and their expression was validated by qRT-PCR (Figure 5). These genes were associated with photosynthesis-antenna proteins, plant hormone signal transduction, phenylpropanoid biosynthesis, ABC transporters, fatty acid elongation, fatty acid metabolism, flavonoid biosynthesis and polycyclic aromatic hydrocarbon degradation. One additional gene associated with each term was targeted, and the transcriptome data and qRT-PCR results showed similar patterns for these selected DEGs derived from RNA-sequencing analysis.

3. Discussion

3.1. Male Sterility in Brassica Crops

Male sterility lines, including Ogura-CMS and DGMS (especially Ogura-CMS), have been widely used in Brassica crops [8,12,23]. The Ogura-CMS line contributes significantly to hybrid seed production in Brassica crops. Therefore, some studies focus on Ogura-CMS, although the mechanism of hybrid seed production has not been elucidated. Previous transcriptome and proteome analyses of the Ogura-CMS line have provided basic information about the DEGs and differential abundant proteins (DAPs) between Ogura-CMS lines and their maintainer lines [12,13,14,24]. However, few studies focus on broccoli and the mechanism of male sterility in the Ogura-CMS, DGMS and maintainer lines. In addition, significant differences of some traits mentioned in this study, especially the maturity periods, flowering time and seed yield, do exist in practice which play an important role in broccoli breeding. Moreover, honey yield and the numbers of bees gathering honey in flowers are usually different in the male sterility and fertile crops [25]. However, there is still no comprehensive study on the mechanism of male sterility including the Ogura-CMS and DGMS lines and their maintainer line in broccoli. In this study, some important agronomic traits were noted and compared among the Ogura-CMS and DGMS lines and their maintainer line based on a two-year investigation. In contrast to previous studies, this study might elucidate the reasons of the stronger performance of plant height, largest leaf length, plant spread angle, single head weight, head width and stem diameter in the Ogura-CMS broccoli than in the DGMS line, and why the DGMS broccoli had longer maturation and flowering periods than the Ogura-CMS and inbred lines. However, the Ogura-CMS broccoli was poorer in seed yield and germination rate than the maintainer line.

3.2. Plant Hormone Signal Transduction and ABC Transporters Affected the Growth Performance and Reproductive Traits

As is known to all, plant hormones play important roles in plant growth and development. In this study, the plant hormone signal transduction and ABC transporters pathways were both significantly enriched for DEGs in the Ogura-CMS line compared to the DGMS and maintainer lines, respectively. To better understand the morphological characterization difference in three broccoli lines, we explored the possible molecular mechanism associated with DEGs of plant hormones.
Auxin, the polar distribution of the plant hormone, is unique and the primary cause for the establishment of local auxin maxima and minima that trigger a plethora of physiological and developmental programs [26]. The generation, control and flexibility of these gradients require the coordinated action of more than three major auxin transporter families including the PIN-FORMED (PIN), the AUXIN1-RESISTANT1 (AUX1)/LIKE AUX1 (AUX1/LAX) and the ABC transporters of the B family (ABCBs) [27,28,29] (Figure 6). Additionally, the subgroup of auxin-transporting ABCBs functions as primary active (ATP-dependent) auxin pumps could transport against steep auxin gradients [30,31]. Thus, ATP is quite necessary for the auxin influx and efflux transporters, and photosynthesis is the initial and powerful basis for the growth and development of plants, especially in plant hormone signal transduction and ABC transporters which could strongly affect all kinds of reproductive and vegetative organs in broccoli and the other eukaryotes [32] (Figure 6). In our study, a number of plant hormone genes of auxin (Aux/IAA and SAUR), CK (A-ARR), ABA (PP2C), ET (ERF1/2), BR (CYCD3) and JA (JAZ) had been detected and highly upregulated in the Ogura-CMS and DGMS lines compared to the maintainer lines, repeatedly. Meanwhile, five OXPHOS-upregulated genes of atp8, LOC106319879, LOC106324734, LOC106314622 and LOC106298585 were also obviously found in the Ogura-CMS and DGMS lines (Figure 6). So, we declared that those functional genes associated with plant hormones and OXPHOS might provide a good explanation for the differences in morphological and transcriptomic characterization in the Ogura-CMS, DGMS and maintainer broccoli lines, and as well as some other Brassica crops.
In addition, with regard to auxin metabolism, the TIR1, GH3 and SAUR families were all obviously upregulated, which could promote cell enlargement during plant growth promotion. The SCF (TIR1) E3 ubiquitin ligase complex is involved in the auxin-mediated signaling pathway that regulates root and hypocotyl growth, lateral root formation, cell elongation and gravitropism. The auxin-responsive GH3 family is essential for plant growth and development in rice, and several genes of this family have been proved to have indoleacetic acid (IAA) amido synthetase biochemical activity [33,34]. The SAUR gene family was initially defined as a set of auxin-inducible genes regulating development, especially in hypocotyls, and many SAUR genes act by promoting cell expansion [35,36]. Additionally, the NPR1 gene associated with SA metabolism was highly upregulated in the Ogura-CMS line rather than the DGMS and maintainer lines, so we inferred that the NPR1 gene might be a domain regulator in the SA-mediated systemic acquired resistance (SAR) pathway and play an essential role in plant immunity [37]. In Arabidopsis, NPR1 can interact with transcription factors to induce the expression of pathogenesis-related (PR) genes and strongly promote defense responses [38], which was consistent with its agronomic characters in this study.

3.3. The DGMS Line Responds to Plant Hormone Signal Transduction

Notably, Arabidopsis plants overexpressing CYCD3-1 show extensive leaf curling, disorganized meristems, increased leaf numbers, late flowering and delayed senescence [39]. In our study, the CYCD3 gene associated with BR metabolism was upregulated in the DGMS broccoli line, which might be an important find for explaining longer maturation and flowering periods than the Ogura-CMS and maintainer fertile lines. Additionally, the J-related gene (JAZ) was present in the DGMS line, while two JA-related genes (JAZ and MYC2) were present in the Ogura-CMS line, which might be act as reliable evidence to explaining the better seed yields in the DGMS line. No SA gene upregulation was found in the DGMS line. Therefore, we found that the BR gene CYCD3 might function in late flowering and delayed senescence of broccoli during the growth and developmental stages consistent with the investigation data obtained over two ecotypes.

4. Materials and Methods

4.1. Plant Materials and Investigation of Agronomic Traits

The broccoli (Brassica oleracea var. italica) maintainer fertile T54S (inbred line F10) has selfed more than 10 generations, and the Ogura-CMS T54C derived from CMSR3629 (Ogura-CMS) introduced by the Asgrow Seed Co. (USA) have also been developed [40,41]. Additionally, the DGMS line T54M was initially derived from the cabbage line 79-399-3 [8,42]. All the materials were bred at the Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences (IVF-CAAS). In the autumn of 2019 and 2020, three replicates with the plant spacing 50 × 45 cm2, were successively grown in the same experimental greenhouse on August 12 at IVF-CAAS (Beijing, China) (40°15′ N, 116°83′ E). When plant harvest began on October 10, three types of bud samples were typically collected from the middle and around the head (diameter 14–15 cm) of T54S, T54M and T54C and labeled with two biological replicates. The isolated buds were immediately frozen in liquid nitrogen and stored at −80 °C for use.

4.2. Scanning Electron Microscopy

The completely open flowers derived from T54S, T54M and T54C, were dissected and dehydrated through a series of increasing ethanol solutions. Additionally, the materials were then critically point-dried with solvent-substituted liquid carbon dioxide and coated with a thin layer of gold palladium. Photomicrographs were obtained using a JEOL 5800 LV at 20 kV (JEOL USA, Peabody, MA, USA). At maturity, the nectaries of broccoli at bolting stage were cut off and pre-fixed in 1% glutaraldehyde and 4% formaldehyde in phosphate buffer (pH 7.2) for 4 h. The nectaries were then post-fixed in 1.5% OsO4 at 4 °C in the same buffer for 4 h, and then, the samples were dehydrated in ascending graded series of acetone and embedded in Spurr’s resin. Ultrathin sections were made using a Reichert ultramicrotome and stained with uranyl acetate and lead citrate [43,44]. The equipment of JEOL-JEM 1200 Ex II transmission electron microscopy (TEM) was used for identifying all the samples at 85 kV.
In addition, the agronomic traits of the broccoli lines T54S, T54M and 54C were investigated during the harvest, flowering, seed ripening and seed germination stages. All the data were analyzed with SPSS 19.0 software (IBM Co., Ltd., New York, NY, USA), and one-way ANOVA was carried out on the agronomic trait data at p < 0.05.

4.3. RNA Extraction and Illumina Sequencing

Total RNA was extracted using an RNA Pure Plant Kit (Tiangen Co., Ltd., Beijing, China) following the manufacturer’s instructions. The integrity of the RNA was assessed with an RNA Nano 6000 Assay Kit on an Agilent Bioanalyzer 2100 platform (Agilent Technologies Inc., Santa Clara, CA, USA). High-quality RNA from each sample was used for cDNA library construction and RNA-Seq on an Illumina HiSeqTM 2500 platform (Berry Genomics Co., Beijing, China). To obtain clean high-quality reads, adapter sequences, low-quality reads (in which >50% bases had Q-values ≤3) and unknown sequences (“N” reads > 10% unknown nucleotides) were removed from the raw reads. The Q20, Q30 and GC contents were then calculated based on the clean reads.

4.4. Quality Control and Transcriptome Analysis

The clean reads from each sample were mapped to the reference genome (ftp://brassicadb.org/Brassica_oleracea/) with TopHat2 (Version 2.0.3.12) [45]. The gene expression levels were normalized using the fragments per kilobase of transcript per million mapped reads (FPKM) method. The edgeR package (Version 2.15.2) was used to identify differentially expressed genes (DEGs) between two samples. Each gene was quantified by calculating the FPKM value with RSEM (V1.2.15). Thresholds of an FDR < 0.05 and a log2 |fold change| >1 were set to select the DEGs. The DEGs were then analyzed for enrichment of Gene Ontology functions and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [17]. GO terms or pathways with p values <0.05 were defined as significantly enriched for the DEGs. Correlations among the expression levels of DEGs were analyzed with Pearson correlation tests [46].

4.5. Quantitative qRT-PCR Analysis

Quantitative real-time (qRT-PCR) analyses were performed to validate the results from the DEGs with the three biological replicates. These are the results of the study: fourteen DEGs were chosen and tested, and six of them were specifically upregulated, and eight of them were specifically downregulated in T54M vs. T54S samples, T54C vs. T54M samples and T54C vs. T54S samples. Total RNA was extracted from the head buds of each sample at harvest time using an RNA Pure Plant Kit (Tiangen Co., Ltd., Beijing, China). First-strand cDNA was synthesized using a PrimeScriptTM 1st Strand cDNA Synthesis Kit (Takara, Dalian, China). qRT-PCR was carried out according to the SYBR PrimeScript RT-PCR Kit manufacturer specifications (Takara, Dalian, China) on an ABI Prism® 7900HT Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). Three biological replicates (with three technical replicates for each biological replicate) were analyzed for each gene, and the relative expression level was estimated by the 2-ΔΔCt method [47,48]. The β-actin gene was used as an internal control, and the specific primers were designed using Primer Premier 6 (Premier Co., Palo Alto, PA, Canada) (Table 3) [47].

4.6. Statistical Analysis

All data were analyzed by analysis of variance (ANOVA) using SPSS Statistics version 19.0 (SPSS, Inc., Chicago, IL, USA) and were expressed as the mean ± standard errors. One-way ANOVA and Tukey’s multiple-range test were used to evaluate the significant differences (p < 0.05). The statistical analyses and charts were produced using the IBM SPSS software and GraphPad Prism 8.2.0 (GraphPad Software Inc., San Diego, CA, USA), respectively.

5. Conclusions

In this study, we firstly significant differences of reproductive traits did exist in three isogenic lines. A total of 5 OXPHOS genes might play role in reproductive traits of the Ogura-CMS broccoli. Meanwhile, 3 photosynthesis genes were found to be highly relative with flowering time in the DGMS broccoli. Our research would provide a comprehensive foundation for understanding the differences of electrophysiological, morphological and transcriptomic profiles in the Ogura-CMS, DGMS and maintainer broccoli, and as well as beneficial to exploring the mechanism of male sterility in Brassica crops.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/plants11040561/s1, Table S1: Summary of the transcriptome sequencing data, Table S2: The KEGG pathways enriched for the DEGs, Table S3: SNP distribution in genomic regions of the Ogura-CMS, DGMS and maintainer broccoli lines, Figure S1: Correlation analysis of gene expression levels between replicate samples of broccoli.

Author Contributions

Conceptualization, Z.L.; writing—original draft preparation, Z.L.; writing—review and editing, L.S., Y.L., F.H. and W.L.; supervision, Z.L.; project administration, Z.L.; funding acquisition, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the National Natural Science Foundation of China (Grant No. 32172580), the Central Public-interest Scientific Institution Basal Research Fund (Grant No. IVF-BRF2021003), the China Agriculture Research System (Grant No. CARS-23-05A), the Agricultural Science and Technology Innovation Program (Grant No. ASTIP), and the State Key Laboratory of Vegetable Germplasm Innovation (Grant No. SKL-VGI).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bai, Y.; Wang, X.; Zhao, S.; Ma, C.; Cui, J.; Zheng, Y. Sulforaphane Protects against Cardiovascular Disease via Nrf2 Activation. Oxidative Med. Cell. Longev. 2015, 2015, 407580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Li, Z.; Liu, Y.; Yuan, S.; Han, F.; Fang, Z.; Yang, L.; Zhuang, M.; Zhang, Y.; Lv, H.; Wang, Y.; et al. Fine mapping of the major QTLs for biochemical variation of sulforaphane in broccoli florets using a DH population. Sci. Rep. 2021, 11, 9004. [Google Scholar] [CrossRef] [PubMed]
  3. Kensler, T.W.; Egner, P.A.; Agyeman, A.S.; Visvanathan, K.; Groopman, J.D.; Chen, J.G.; Chen, T.Y.; Fahey, J.W.; Talalay, P. Keap1-nrf2 signaling: A target for cancer prevention by sulforaphane. Top. Curr. Chem. 2013, 329, 163–177. [Google Scholar] [PubMed] [Green Version]
  4. Appendino, G.; Bardelli, A. Broccoli, PTEN deletion and prostate cancer: Where is the link? Mol. Cancer 2010, 9, 1–3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Li, Z.; Zheng, S.; Liu, Y.; Fang, Z.; Yang, L.; Zhuang, M.; Zhang, Y.; Lv, H.; Wang, Y.; Xu, D. Characterization of glucosinolates in 80 broccoli genotypes and different organs using UHPLC-Triple-TOF-MS method. Food Chem. 2021, 334, 127519. [Google Scholar] [CrossRef]
  6. Li, Z.S.; Mei, Y.J.; Liu, Y.M.; Fang, Z.Y.; Yang, L.M.; Zhuang, M.; Zhang, Y.Y.; Lv, H.H. The evolution of genetic diversity of broccoli cultivars in China since 1980. Sci. Hortic. Amst. 2019, 250, 69–80. [Google Scholar] [CrossRef]
  7. Huang, J.; Liu, Y.; Han, F.; Fang, Z.; Yang, L.; Zhuang, M.; Zhang, Y.; Lv, H.; Wang, Y.; Ji, J.; et al. Genetic diversity and population structure analysis of 161 broccoli cultivars based on SNP markers. Hortic. Plant J. 2021, 7, 423–433. [Google Scholar] [CrossRef]
  8. Fang, Z.Y.; Sun, P.T.; Liu, Y.M.; Yang, L.M.; Wang, X.W.; Hou, A.F.; Bian, C.S. A male sterile line with dominant gene (Ms) in cabbage (Brassica oleracea var. capitata) and its utilization for hybrid seed production. Euphytica 1997, 97, 265–268. [Google Scholar] [CrossRef]
  9. Kaminski, P.; Podwyszynska, M.; Starzycki, M.; Starzycka-Korbas, E. Interspecific hybridisation of cytoplasmic male-sterile rapeseed with Ogura cytoplasm and Brassica rapa var. pekinensis as a method to obtain male-sterile Chinese cabbage inbred lines. Euphytica 2016, 208, 519–534. [Google Scholar] [CrossRef] [Green Version]
  10. Kirti, P.B.; Banga, S.S.; Prakash, S.; Chopra, V.L. Transfer of Ogu Cytoplasmic Male-Sterility to Brassica-Juncea and Improvement of the Male-Sterile Line through Somatic-Cell Fusion. Appl. Genet. 1995, 91, 517–521. [Google Scholar] [CrossRef]
  11. Chamola, R.; Balyan, H.S.; Bhat, S.R. Transfer of cytoplasmic male sterility from alloplasmic Brassica juncea and B. napus to cauliflower (B. oleracea var. botrytis) through interspecific hybridization and embryo culture. Indian J. Genet. Plant Breed. 2013, 73, 203–210. [Google Scholar] [CrossRef]
  12. Shu, J.S.; Zhang, L.L.; Liu, Y.M.; Li, Z.S.; Fang, Z.Y.; Yang, L.M.; Zhuang, M.; Zhang, Y.Y.; Lv, H.H. Normal and Abortive Buds Transcriptomic Profiling of Broccoli ogu Cytoplasmic Male Sterile Line and Its Maintainer. Int. J. Mol. Sci. 2018, 19, 2501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Han, F.Q.; Zhang, X.L.; Yang, L.M.; Zhuang, M.; Zhang, Y.Y.; Li, Z.S.; Fang, Z.Y.; Lv, H.H. iTRAQ-Based Proteomic Analysis of Ogura-CMS Cabbage and Its Maintainer Line. Int. J. Mol. Sci. 2018, 19, 3180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Xing, M.M.; Sun, C.; Li, H.L.; Hu, S.L.; Lei, L.; Kang, J.G. Integrated analysis of transcriptome and proteome changes related to the Ogura cytoplasmic male sterility in cabbage. PLoS ONE 2018, 13, e0193462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Li, H.; Wang, Y.; Wu, M.; Li, L.H.; Jin, C.; Zhang, Q.L.; Chen, C.B.; Song, W.Q.; Wang, C.G. Small RNA Sequencing Reveals Differential miRNA Expression in the Early Development of Broccoli (Brassica oleracea var. italica) Pollen. Front. Plant Sci. 2017, 8, 404. [Google Scholar] [CrossRef] [Green Version]
  16. Luo, F.; Cai, J.H.; Kong, X.M.; Zhou, Q.; Zhou, X.; Zhao, Y.B.; Ji, S.J. Transcriptome profiling reveals the roles of pigment mechanisms in postharvest broccoli yellowing. Hortic. Res.-Engl. 2019, 6, 1–14. [Google Scholar] [CrossRef] [Green Version]
  17. Li, Z.S.; Liu, Y.M.; Li, L.Y.; Fang, Z.Y.; Yang, L.M.; Zhuang, M.; Zhang, Y.Y.; Lv, H.H. Transcriptome reveals the gene expression patterns of sulforaphane metabolism in broccoli florets. PLoS ONE 2019, 14, e0213902. [Google Scholar] [CrossRef]
  18. Ku, K.M.; Becker, T.M.; Juvik, J.A. Transcriptome and Metabolome Analyses of Glucosinolates in Two Broccoli Cultivars Following Jasmonate Treatment for the Induction of Glucosinolate Defense to Trichoplusia ni (Hubner). Int. J. Mol. Sci. 2016, 17, 1135. [Google Scholar] [CrossRef] [Green Version]
  19. Li, H.; Yuan, J.Y.; Wu, M.; Han, Z.P.; Li, L.H.; Jiang, H.M.; Jia, Y.L.; Han, X.; Liu, M.; Sun, D.L.; et al. Transcriptome and DNA methylome reveal insights into yield heterosis in the curds of broccoli (Brassica oleracea L var. italic). Bmc Plant Biol. 2018, 18, 1–17. [Google Scholar] [CrossRef] [Green Version]
  20. Jiang, Z.; Guo, H. A comparative genomic analysis of plant hormone related genes in different species. J. Genet. Genom. 2010, 37, 219–230. [Google Scholar] [CrossRef]
  21. Muto, H.; Nagao, I.; Demura, T.; Fukuda, H.; Kinjo, M.; Yamamoto, K.T. Fluorescence cross-correlation analyses of the molecular interaction between an Aux/IAA protein, MSG2/IAA19, and protein-protein interaction domains of auxin response factors of Arabidopsis expressed in HeLa cells. Plant Cell Physiol. 2006, 47, 1095–1101. [Google Scholar] [CrossRef] [PubMed]
  22. Braun, H.P. The Oxidative Phosphorylation system of the mitochondria in plants. Mitochondrion 2020, 53, 66–75. [Google Scholar] [CrossRef] [PubMed]
  23. Wu, J.Y.; Shen, J.R.; Mao, X.Z.; Liu, K.D.; Wei, L.P.; Liu, P.W.; Yang, G.S. Isolation and analysis of differentially expressed genes in dominant genic male sterility (DGMS) Brassica napus L. using subtractive PCR and cDNA microarray. Plant Sci. 2007, 172, 204–211. [Google Scholar] [CrossRef]
  24. Wei, X.C.; Zhang, X.H.; Yao, Q.J.; Yuan, Y.X.; Li, X.X.; Wei, F.; Zhao, Y.Y.; Zhang, Q.; Wang, Z.Y.; Jiang, W.S.; et al. The miRNAs and their regulatory networks responsible for pollen abortion in Ogura-CMS Chinese cabbage revealed by high-throughput sequencing of miRNAs, degradomes, and transcriptomes. Front. Plant Sci. 2015, 6, 894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Abdelhafiz, A.T.; Muhamad, J.A. Midcycle pericoital intravaginal bee honey and royal jelly for male factor infertility. Int. J. Gynaecol. Obs. 2008, 101, 146–149. [Google Scholar] [CrossRef] [PubMed]
  26. Hao, P.; Xia, J.; Liu, J.; Di Donato, M.; Pakula, K.; Bailly, A.; Jasinski, M.; Geisler, M. Auxin-transporting ABC transporters are defined by a conserved D/E-P motif regulated by a prolylisomerase. J. Biol. Chem. 2020, 295, 13094–13105. [Google Scholar] [CrossRef] [PubMed]
  27. Naramoto, S. Polar transport in plants mediated by membrane transporters: Focus on mechanisms of polar auxin transport. Curr. Opin. Plant Biol. 2017, 40, 8–14. [Google Scholar] [CrossRef]
  28. Müller, K.; Hošek, P.; Laňková, M.; Vosolsobě, S.; Malínská, K.; Čarná, M.; Fílová, M.; Dobrev, P.I.; Helusová, M.; Hoyerová, K.; et al. Transcription of specific auxin efflux and influx carriers drives auxin homeostasis in tobacco cells. Plant J. Cell Mol. Biol. 2019, 100, 627–640. [Google Scholar] [CrossRef]
  29. Lee, S.H.; Cho, H.-T. PINOID positively regulates auxin efflux in Arabidopsis root hair cells and tobacco cells. Plant Cell 2006, 18, 1604–1616. [Google Scholar] [CrossRef] [Green Version]
  30. Anfang, M.; Shani, E. Transport mechanisms of plant hormones. Curr Opin Plant Biol 2021, 63, 102055. [Google Scholar] [CrossRef]
  31. Grones, P.; Friml, J. Auxin transporters and binding proteins at a glance. J. Cell Sci. 2015, 128, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Skokan, R.; Medvecká, E.; Viaene, T.; Vosolsobě, S.; Zwiewka, M.; Müller, K.; Skůpa, P.; Karady, M.; Zhang, Y.; Janacek, D.P.; et al. PIN-driven auxin transport emerged early in streptophyte evolution. Nat. Plants 2019, 5, 1114–1119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Du, H.; Wu, N.; Fu, J.; Wang, S.P.; Li, X.H.; Xiao, J.H.; Xiong, L.Z. A GH3 family member, OsGH3-2, modulates auxin and abscisic acid levels and differentially affects drought and cold tolerance in rice. J. Exp. Bot. 2012, 63, 6467–6480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Yang, J.H.; Han, S.J.; Yoon, E.K.; Lee, W.S. Evidence of an auxin signal pathway, microRNA167-ARF8-GH3, and its response to exogenous auxin in cultured rice cells. Nucleic Acids Res. 2006, 34, 1892–1899. [Google Scholar] [CrossRef]
  35. Xu, N.F.; Hagen, G.; Guilfoyle, T. Multiple auxin response modules in the soybean SAUR 15A promoter. Plant Sci. 1997, 126, 193–201. [Google Scholar] [CrossRef]
  36. Spartz, A.K.; Ren, H.; Park, M.Y.; Grandt, K.N.; Lee, S.H.; Murphy, A.S.; Sussman, M.R.; Overvoorde, P.J.; Gray, W.M. SAUR Inhibition of PP2C-D Phosphatases Activates Plasma Membrane H+-ATPases to Promote Cell Expansion in Arabidopsis. Plant Cell 2014, 26, 2129–2142. [Google Scholar] [CrossRef] [Green Version]
  37. Nic-Matos, G.; Narvaez, M.; Peraza-Echeverria, S.; Saenz, L.; Oropeza, C. Molecular cloning of two novel NPR1 homologue genes in coconut palm and analysis of their expression in response to the plant defense hormone salicylic acid. Genes Genom. 2017, 39, 1007–1019. [Google Scholar] [CrossRef]
  38. Chen, J.; Mohan, R.; Zhang, Y.; Li, M.; Chen, H.; Palmer, I.A.; Chang, M.; Qi, G.; Spoel, S.H.; Mengiste, T.; et al. NPR1 Promotes Its Own and Target Gene Expression in Plant Defense by Recruiting CDK8. Plant Physiol. 2019, 181, 289–304. [Google Scholar] [CrossRef]
  39. Menges, M.; Samland, A.K.; Planchais, S.; Murray, J.A.H. The D-type cyclin CYCD3;1 is limiting for the G1-to-S-phase transition in Arabidopsis. Plant Cell 2006, 18, 893–906. [Google Scholar] [CrossRef] [Green Version]
  40. Dong, X.; Kim, W.K.; Lim, Y.P.; Kim, Y.K.; Hur, Y. Ogura-CMS in Chinese cabbage (Brassica rapa ssp. pekinensis) causes delayed expression of many nuclear genes. Plant Sci. 2013, 199–200, 7–17. [Google Scholar] [CrossRef]
  41. Yu, H.L.; Fang, Z.Y.; Liu, Y.M.; Yang, L.M.; Zhuang, M.; Lv, H.H.; Li, Z.S.; Han, F.Q.; Liu, X.P.; Zhang, Y.Y. Development of a novel allele-specific Rfo marker and creation of Ogura CMS fertility-restored interspecific hybrids in Brassica oleracea. Appl. Genet. 2016, 129, 1625–1637. [Google Scholar] [CrossRef] [PubMed]
  42. Chen, C.; Zhuang, M.; Fang, Z.Y.; Wang, Q.B.; Zhang, Y.Y.; Liu, Y.M.; Yang, L.M.; Cheng, F. A Co-Dominant Marker BoE332 Applied to Marker-Assisted Selection of Homozygous Male-Sterile Plants in Cabbage (Brassica oleracea var. capitata L.). J. Integr. Agric. 2013, 12, 596–602. [Google Scholar] [CrossRef]
  43. Solís, S.M.; Zini, L.M.; González, V.V.; Ferrucci, M.S. Floral nectaries in Sapindaceae s.s.: Morphological and structural diversity, and their systematic implications. Protoplasma 2017, 254, 2169–2188. [Google Scholar] [CrossRef] [PubMed]
  44. Zini, L.M.; Solís, S.M.; Ferrucci, M.S. Anatomical and developmental studies on floral nectaries in Cardiospermum species: An approach to the evolutionary trend in Paullinieae. Plant Syst. Evol. 2014, 300, 1515–1523. [Google Scholar] [CrossRef]
  45. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [Green Version]
  46. Brysbaert, G.; Lorgouilloux, K.; Vranken, W.F.; Lensink, M.F. RINspector: A Cytoscape app for centrality analyses and DynaMine flexibility prediction. Bioinformatics 2018, 34, 294–296. [Google Scholar] [CrossRef]
  47. Li, Z.S.; Liu, Y.M.; Fang, Z.Y.; Yang, L.M.; Zhuang, M.; Zhang, Y.Y.; Zhao, W.; Sun, P.T. Variation of Sulforaphane Levels in Broccoli (Brassica Oleracea Var. Italica) during Flower Development and the Role of Gene Aop2. J. Liq. Chromatogr. Relat. Technol. 2014, 37, 1199–1211. [Google Scholar] [CrossRef]
  48. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Figure 1. Agronomic traits of the flower, stamen and pistil, and the nectary morphology of the Ogura-CMS line T54C, the DGMS line T54M and their maintainer line T54S. The red arrows show the honey gland.
Figure 1. Agronomic traits of the flower, stamen and pistil, and the nectary morphology of the Ogura-CMS line T54C, the DGMS line T54M and their maintainer line T54S. The red arrows show the honey gland.
Plants 11 00561 g001
Figure 2. The DEGs identified in the Ogura-CMS line, DGMS and maintainer lines. (a) Clustering heat map of DEGs in T54C and T54S, T54C and T54M, and T54M and T54S, green and red colors indicated up- and downregulation of the genes, respectively. (b) Volcano plot of DEGs in comparisons of T54M vs. T54S, T54C vs. T54S, T54C vs. T54M, green and red colors indicated up- and downregulation of the genes, respectively. (c) The DEGs among broccoli lines of T54C vs. T54S, T54C vs. T54m, and T54M vs. T54S.
Figure 2. The DEGs identified in the Ogura-CMS line, DGMS and maintainer lines. (a) Clustering heat map of DEGs in T54C and T54S, T54C and T54M, and T54M and T54S, green and red colors indicated up- and downregulation of the genes, respectively. (b) Volcano plot of DEGs in comparisons of T54M vs. T54S, T54C vs. T54S, T54C vs. T54M, green and red colors indicated up- and downregulation of the genes, respectively. (c) The DEGs among broccoli lines of T54C vs. T54S, T54C vs. T54m, and T54M vs. T54S.
Plants 11 00561 g002
Figure 3. KEGG pathway analysis of the Ogura-CMS (a) and DGMS (b) lines vs. the maintainer line T54S. The red star stands for plant hormone signal transduction pathway and the green stars stands for photosynthesis-antenna proteins and photosynthesis pathways.
Figure 3. KEGG pathway analysis of the Ogura-CMS (a) and DGMS (b) lines vs. the maintainer line T54S. The red star stands for plant hormone signal transduction pathway and the green stars stands for photosynthesis-antenna proteins and photosynthesis pathways.
Plants 11 00561 g003
Figure 4. The oxidative phosphorylation pathway (a) and five upregulated DEGs with relative expression (FPKM) (b). The red frame in oxidative phosphorylation pathway represents upregulated genes (a) and red statements in (b) indicate the functional annotation. (*) indicated significant difference at p < 0.05 (b).
Figure 4. The oxidative phosphorylation pathway (a) and five upregulated DEGs with relative expression (FPKM) (b). The red frame in oxidative phosphorylation pathway represents upregulated genes (a) and red statements in (b) indicate the functional annotation. (*) indicated significant difference at p < 0.05 (b).
Plants 11 00561 g004
Figure 5. Verification of the selected DEGs carried out by qRT-PCR. Fourteen DEGs expression in T54S, T54C and T54M lines were selected for qRT-PCR validation. The relative expression level of each gene was shown as the FPKM between different two samples in the RNA-Seq data (red bar) and qRT-PCR results (blue line).
Figure 5. Verification of the selected DEGs carried out by qRT-PCR. Fourteen DEGs expression in T54S, T54C and T54M lines were selected for qRT-PCR validation. The relative expression level of each gene was shown as the FPKM between different two samples in the RNA-Seq data (red bar) and qRT-PCR results (blue line).
Plants 11 00561 g005
Figure 6. The role of photosynthesis, auxin transporters and plant hormones in reproductive and vegetative organs of broccoli. The genes in stuffed red frames surrounded by green dotted line were from pathways of OXPHOS (5), photosynthesis-antenna protein (3) and plant hormones (1). The orange box indicates plant hormones and the dummy arrow shows their direction of transport, while the black line indicates their synthetic organs or tissues. The stuffed olive frames indicate the functions of plant hormones in different organs and tissues. The red round stands for apical point of broccoli and green round in florets stands for one bud.
Figure 6. The role of photosynthesis, auxin transporters and plant hormones in reproductive and vegetative organs of broccoli. The genes in stuffed red frames surrounded by green dotted line were from pathways of OXPHOS (5), photosynthesis-antenna protein (3) and plant hormones (1). The orange box indicates plant hormones and the dummy arrow shows their direction of transport, while the black line indicates their synthetic organs or tissues. The stuffed olive frames indicate the functions of plant hormones in different organs and tissues. The red round stands for apical point of broccoli and green round in florets stands for one bud.
Plants 11 00561 g006
Table 1. Investigation on agronomic traits of broccoli lines T54S, T54M and T54C in the autumn of 2019 and 2020.
Table 1. Investigation on agronomic traits of broccoli lines T54S, T54M and T54C in the autumn of 2019 and 2020.
Trait Name2019 2020
T54ST54MT54CT54ST54MT54C
Plant height (cm) *72.48 ± 4.25 b70.33 ± 4.29 c73.27 ± 3.14 a68.67 ± 3.65 a65.39 ± 3.74 c68.18 ± 2.85 b
Length of the largest leaf (cm) *60.37 ± 2.49 b58.56 ± 3.87 c61.48 ± 3.04 a58.88 ± 3.07 b56.49 ± 2.28 c59.71 ± 3.63 a
Width of the largest leaf (cm)25.28 ± 2.79 a24.17 ± 1.89 b25.08 ± 2.12 a24.22 ± 1.75 b24.01 ± 2.01 b25.53 ± 2.28 a
Plant spread angle (cm) *91.04 ± 5.85 b86.34 ± 6.15 c92.86 ± 6.19 a88.18 ± 5.35 b86.27 ± 6.05 c90.84 ± 7.56 a
Single head weight (kg) *0.24 ± 0.04 b0.23 ± 0.07 c0.26 ± 0.04 a0.22 ± 0.05 b0.21 ± 0.08 c0.24 ± 0.05 a
Head height (cm)15.24 ± 1.16 a15.26 ± 2.06 a15.36 ± 2.18 a15.11 ± 1.89 a14.84 ± 1.78 a15.26 ± 2.24 a
Head width (cm) *15.37 ± 1.55 b15.22 ± 2.08 b16.04 ± 2.57 a15.08 ± 1.87 a14.39 ± 2.06 b15.53 ± 3.15 a
Stem diameter (cm) *5.12 ± 0.18 a4.98 ± 0.27 b5.16 ± 0.42 a4.87 ± 0.38 a4.35 ± 0.18 b4.96 ± 0.49 a
Stem hollow/solidSolidSolidSolidSolidSolidSolid
Leaf number10 ± 2 a9 ± 2 a10 ± 2 a9 ± 1 a10 ± 3 a9 ± 2 a
Leaf colorGreenGreenGreenGreenGreenGreen
Leaf waxinessHighHighHighHighHighHigh
Flower colorYellowYellowYellowYellowYellowYellow
Days to maturity (day)*105 ± 3 b109 ± 4 a106 ± 4 b126 ± 3 b134 ± 4 a128 ± 4 b
Days to flowering (day) *126 ± 3 b129 ± 3 a126 ± 4 b143 ± 2 b147 ± 3 a142 ± 2 b
Seeds yield per plant (g) *4.13 ± 0.24 a3.52 ± 0.58 b2.41 ± 0.27 c2.86 ± 0.29 a2.21 ± 0.56 b1.74 ± 0.36 c
Seeds germination rate (%) *86.58 ± 3.77 a83.04 ± 4.18 b71.76 ± 2.47 c85.03 ± 2.61 a79.59 ± 2.92 b71.46 ± 2.88 c
Note: T54C, T54M and T54S represented samples from the Ogura-CMS, DGMS and maintainer lines, respectively. All the data were shown as the mean ± standard deviation (n = 3). Days to maturity represented the number of days from sowing to head maturity, and days to flowering represented the number of days from sowing to flowering. The asterisk (*) and minuscule show that a significant difference was identified among the T54S, T54M and T54C groups in the same season based on one-way ANOVA at p < 0.05.
Table 2. Correlation analysis of gene expression levels between replicate samples of broccoli.
Table 2. Correlation analysis of gene expression levels between replicate samples of broccoli.
Serial NumberSample NameLinear EquationPearson Test (R2)
1T54C-1 vs. T54C-2_FPKMY = 0.965*X + 0.0250.943
2T54C-1 vs. T54M-1_FPKMY = 0.939*X − 0.0530.918
3T54C-1 vs. T54M-2_FPKMY = 0.915*X + 0.0460.875
4T54C-1 vs. T54S-1_FPKMY = 0.913*X + 0.1770.867
5T54C-1 vs. T54S-2_FPKMY = 0.937*X + 0.1500.890
6T54C-2 vs. T54M-1_FPKMY = 0.945*X − 0.0730.918
7T54C-2 vs. T54M-2_FPKMY = 0.917*X + 0.0280.867
8T54C-2 vs. T54S-1_FPKMY = 0.918*X + 0.1580.867
9T54C-2 vs. T54S-2_FPKMY = 0.946*X + 0.1300.905
10T54M-1 vs. T54M-2_FPKMY = 0.960*X + 0.0990.925
11T54M-1 vs. T54S-1_FPKMY = 0.947*X + 0.2310.896
12T54M-1 vs. T54S-2_FPKMY = 0.973*X + 0.2050.973
13T54M-2 vs. T54S-1_FPKMY = 0.947*X + 0.1410.894
14T54M-2 vs. T54S-2_FPKMY = 0.967*X + 0.0130.917
15T54S-1 vs. T54S-2_FPKMY = 0.966*X − 0.0130.919
Note: T54C, T54M and T54S represented the samples of the Ogura-CMS, DGMS and maintainer lines, respectively. Y and X represented the values of log2 (FPKM) between two samples, and * indicates that correlation was highly significant between two samples at p < 0.05.
Table 3. Genes targeted by the specific primers for qRT-PCR analysis.
Table 3. Genes targeted by the specific primers for qRT-PCR analysis.
GeneFunction DescriptionChromosomeForward Primer 5′-3′Reverse Primer 5′-3′Ta (°C)
LOC106318096Chlorophyll a-b binding protein 13, chloroplasticC9TCTCCACCAAACCAGCAAAGTCAAGACCGTTGATGCGGAA59
LOC106313340Auxin-responsive protein IAA4C9CTTGGGTTTTCGGGGGAAGACCGTTTAAGCCTCTCACTGGT60
LOC106316002Ethylene receptor 2 (BO-ETR2)C1TGGGTTTAAGCATTTTCATGGGAGAAAGGGTGGGGACCGTAAG60
LOC106296937Cinnamoyl-CoA reductase 2-likeC6AACGGAGCCAAGTTCGTGATCGAAAAGAAGAGGGCATTCAGC60
LOC106299021Cytochrome P450 98A8-likeC6TGGTGGGCATCCAACATACCTTGTCGCTAACGAGCCACTT60
LOC106307625ABC transporter B family member 2C1GGAGGCGTCCAGTGATCCTCCAACGAGTACTTGGCGAC60
LOC106305741ABC transporter B family member 19C7AGATCATCCGGGGAGGTGAATCATCAAGCATCAAATCATAAGCGT60
LOC1063122873-ketoacyl-CoA synthase 1 (KCS)C8GCCGTCTCTATCGGCAATGACACTCCACATATCCGGTCCA60
LOC106304303Palmitoyl-protein thioesterase 1C7TTCGGGTTTTACCCGGATGGTACGCAAAGGCTCCTTGGTT60
LOC1063256113-oxoacyl-[acyl-carrier-protein] synthase II, chloroplastic-likeC2TGGCTGCCTCTTCCTGTTACCCGGAGTTAGTTGCTCGGTT59
LOC106328481Chalcone synthase 3-likeC3CCAAGCTCCTTGGTCTTCGTACGATGCTGGGAAGTATGTGT59
LOC106298325Chalcone--flavonone isomerase-likeC6TTGGCTACCTACCTCCTCCCTCAAACGAACCCGTGACGAT60
LOC106306241Cytochrome P450 81D11C7ACGAATCTGCGAAGGTGGAGTGCCTGCATCTGCGGAATAA61
LOC106298990Cytochrome P450 71B5C6CCGACCTAAGACGGTAGGGACCGGAGATCCGGTCAATGAG60
β-actinBeta-actin ATCTGGCATCACACTTTCTACATCTCTTTGCTCATACGGTCT55
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Z.; Song, L.; Liu, Y.; Han, F.; Liu, W. Electrophysiological, Morphologic, and Transcriptomic Profiling of the Ogura-CMS, DGMS and Maintainer Broccoli Lines. Plants 2022, 11, 561. https://0-doi-org.brum.beds.ac.uk/10.3390/plants11040561

AMA Style

Li Z, Song L, Liu Y, Han F, Liu W. Electrophysiological, Morphologic, and Transcriptomic Profiling of the Ogura-CMS, DGMS and Maintainer Broccoli Lines. Plants. 2022; 11(4):561. https://0-doi-org.brum.beds.ac.uk/10.3390/plants11040561

Chicago/Turabian Style

Li, Zhansheng, Lixiao Song, Yumei Liu, Fengqing Han, and Wei Liu. 2022. "Electrophysiological, Morphologic, and Transcriptomic Profiling of the Ogura-CMS, DGMS and Maintainer Broccoli Lines" Plants 11, no. 4: 561. https://0-doi-org.brum.beds.ac.uk/10.3390/plants11040561

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop