Next Article in Journal
Absence of Figla-like Gene Is Concordant with Femaleness in Cichlids Harboring the LG1 Sex-Determination System
Next Article in Special Issue
A Large-Scale Genomic Association Analysis Identifies the Candidate Genes Regulating Salt Tolerance in Cucumber (Cucumis sativus L.) Seedlings
Previous Article in Journal
Nano-Formulation Endows Quorum Quenching Enzyme-Antibiotic Hybrids with Improved Antibacterial and Antibiofilm Activities against Pseudomonas aeruginosa
Previous Article in Special Issue
Biosynthesis and Metabolism of Garlic Odor Compounds in Cultivated Chinese Chives (Allium tuberosum) and Wild Chinese Chives (Allium hookeri)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 14
10.3390/ijms23147633
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fine Mapping and Gene Analysis of restorer-of-fertility Gene CaRfHZ in Pepper (Capsicum annuum L.)

by
Zhixing Nie
1,2,3,†,
Yunpeng Song
1,†,
Hong Wang
2,
Jianying Chen
2,
Qingliang Niu
3,* and
Weimin Zhu
1,*
1
Shanghai Key Laboratory of Protected Horticultural Technology, Horticulture Research Institute, Shanghai Academy of Agricultural Sciences, Shanghai 201403, China
2
Vegetable Research Institute, Hangzhou Academy of Agricultural Sciences, Hangzhou 310024, China
3
School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai 200240, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(14), 7633; https://doi.org/10.3390/ijms23147633
Submission received: 7 June 2022 / Revised: 7 July 2022 / Accepted: 9 July 2022 / Published: 11 July 2022
(This article belongs to the Special Issue Vegetable Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Cytoplasmic male sterility (CMS) is a common biological phenomenon used in hybrid production of peppers (Capsicum annuum L.). Although several restorer-of-fertility (Rf) genes of pepper CMS lines have been mapped, there is no report that the Rf gene with clear gene function has been isolated. Here, pepper CMS line HZ1A and its restorer line HZ1C were used to construct (HZ1A × HZ1C) F2 populations and map the Rf gene. A single dominant gene CaRfHZ conferred male fertility according to inheritance analysis. Using sterile plants from (HZ1A × HZ1C) F2 populations and bulked segregant analysis (BSA), the CaRfHZ gene was mapped between P06gInDel-66 and P06gInDel-89 on chromosome 6. This region spans 533.81 kb, where four genes are annotated according to Zunla-1 V2.0 gene models. Based on the analysis of genomic DNA sequences, gene expressions, and protein structures, Capana06g002968 was proposed as the strongest candidate for the CaRfHZ gene. Our results may help with hybrid pepper breeding and to elucidate the mechanism of male fertility restoration in peppers.

1. Introduction

Cytoplasmic male sterility (CMS) is a common biological phenomenon observed in plants, and it is widely used for seed production of hybrid crops. Fertility restoration of CMS is critical in hybrid production, and restorer-of-fertility (Rf) gene can restore and affect the fertility of hybrids. More than ten Rf genes have been isolated and studied in major crops [1] such as rice, maize, rapeseed, wheat, and soybean. Most of the Rf genes in crops contain multiple pentatricopeptide repeat (PPR) domains, such as Rf1 [2], Rf4 [3], Rf5 [4], and Rf6 [5] in rice; Rf1 [6] and Rf3 [7] in maize; Rfp [8] in rapeseed; Rfo [9] in radish; RFL79 and RFL29a [10] in wheat; and GmPPR576 [11] in soybean. Additionally, certain Rf genes do not contain PPR domains, such as Rf2 [6] and Rf4 [12] in maize. Rf2 is similar to the mammalian aldehyde dehydrogenase, and Rf4 is a transcription factor that contains basic helix–loop–helix (bHLH) domains. Most Rf genes restore the fertility of sterile lines by affecting the transcription of sterile genes in mitochondria. For example, the rice gene Rf4 restores fertility by reducing the transcripts of the sterility gene orf352 in the cytoplasm of wild-type rice [4,13]; the rice gene Rf6 reduces the accumulation of the transcript atp6-orfH79, thereby restoring fertility [5]; the maize gene Rf3 inhibits mitochondrial gene orf77 editing and degradation, and accelerates orf355 degradation, leading to CMS-S fertility restoration [14]; the rapeseed gene Rfp restores fertility of the sterile line by processing the transcript of the mitochondrial sterility gene atp6-orf224 [8].
Pepper is one of the largest vegetable crops globally. Peterson (1958) [15] was the first to report cytoplasmic–nucleus interaction type male sterility in pepper, following which, breeders bred numerous sterile pepper lines using CMS lines [16]. Due to the different fertility segregating populations used for genetic analysis in different studies, the inheritance pattern of pepper Rf gene also has different conclusions. Wang et al. [17] used the male sterile line 77013A to cross 114 doubled haploid lines to identify the fertility of F1 offspring, and considered that fertility restoration of pepper CMS was controlled by one major and four minor genes. Wei et al. [18] analyzed the segregation of fertility using an F2 population. The fertility restoration of pepper CMS may be conferred by two pairs of additive-dominant epistatic major genes and additive-dominant polygenes. Gulyas et al. [19] and Ye et al. [20] investigated the fertility separation of F1–F4 and F2 populations constructed using male sterile and restorer lines, respectively. The results showed that a pair of dominant nuclear genes controls fertility restoration in CMS.
In early studies, owing to the limitations of molecular biology techniques and the absence of a reference genome, the Rf genes of pepper were only initially mapped. Zhang et al. [21] used the F2 population to locate the major Rf gene between two randomly amplified polymorphic DNA (RAPD) markers using bulked segregant analysis (BSA), and the genetic distance was 0.37 centimorgans (cM) and 8.12 cM. Wang et al. [17] mapped a major quantitative trait locus (QTL) of fertility restoration to chromosome 6 and four minor QTLs on chromosomes 1, 2, 3, and 5. Kim et al. [22] used BSA–amplified fragment length polymorphism (AFLP) to segregate populations by fertility, constructed a linkage map of the Rf gene, and obtained a cleaved amplified polymorphic sequence (CAPS) marker AFRF8CAPS with a genetic distance of 1.8 cM from the Rf gene. Lee et al. [23] used 205 F2 individual plants and obtained three CAPS and one sequence characterized amplified region (SCAR) markers that were closely linked to the Rf site of pepper CMS. Yang et al. [24] used the F2 population to obtain the simple sequence repeat (SSR) molecular marker AF208834 linked to the pepper Rf gene, with a genetic distance of 20.8 cM. Jo et al. [25] screened the pepper BAC library with the petunia Rf gene, cloned the Rf candidate gene PePPR1, and found three markers closely linked to the Rf gene.
Along with the in-depth research and publication of the pepper reference genome [26,27], the pepper Rf gene has been mapped to an adjacent region of chromosome 6 in different reports (Table 1). Jo et al. [28] used BSA–AFLP and comparative genomics to locate the Rf gene within the 821 kb segment on chromosome 6 and identified the candidate gene CaPPR6. Zhang et al. [29] narrowed down the region of CaPPR6 to 128.96 Kb using kompetitive allele-specific PCR (KASP) markers. Kang et al. [30] mapped an unstable Rf gene, Rfu, using the molecular markers from Jo et al. [28] and newly developed CAPS markers. Ye et al. [20] located the Rf gene between SSR markers pep43 and pep20 on chromosome 6, with a physical distance of 498.6 kb. Wu et al. [31] performed association analysis on 287 pepper lines using specific-locus amplified fragment sequencing (SLAF-seq) and genome-wide association study (GWAS) and identified two candidate Rf genes, Capana06g002967 and Capana06g002969. Cheng et al. [32] used high-density genetic mapping and collinearity analysis to map the Rf gene CaRf to a 270.10 kb region on chromosome 6, and analyzed the candidate gene Capana06g003028. Zhang et al. [33] used high-throughput sequencing combined with BSA to initially map the Rf gene in the F2 population, then mapped the Rf gene to a 148.05 kb segment on chromosome 6, and identified the candidate gene CaRf032 (CA00g82510), which has a base variation in the maintainer line leading to premature termination of transcription of the gene. Wei et al. [34] performed RNA sequencing on fertile and sterile pools constructed from the F2 population and found a candidate gene, NEDD8 (Capana06g002866), on chromosome 6. However, to date, there is no report that the Rf gene with a clear gene function of the pepper CMS line has been cloned.
In the present study, we constructed an F2 segregating population and fine-mapped the Rf gene in the pepper line HZ1C. The Rf gene was mapped to a 533.81 kb region between two insertion or deletion (InDel) markers on chromosome 6 of the Zunla-1 V2.0 reference genome [26]. According to DNA sequence alignment, gene expression, and protein structure analysis, Capana06g002968 was proposed as the strongest candidate gene for the Rf gene in pepper line HZ1C.

2. Results

2.1. The Inheritance Analysis of Fertility Restoration in CMS Line HZ1A

HZ1A shows the typical characteristics of male sterility, including no pollen dissemination around its anthers and fewer and shriveled pollen grains [35]. The fertility phenotypes were identified through visible pollen on anthers and Alexander’s staining of pollen. All (HZ1A × HZ1C) F1 plants were completely male fertile. The segregation of male-fertile to male-sterile plants in 11 F2 populations of different sizes, which had been planted in the autumn and spring of 2020 and 2021, respectively, conformed to a 3:1 segregation ratio under the χ2 criterion (p > 0.05). Overall, 1658 F2 plants were classified as 1277 male-fertile and 381 male-sterile plants, which also conformed to a 3:1 segregation ratio (χ2 = 3.50, p = 0.06) (Table 2). These results demonstrated that a single dominant gene conferred male fertility in HZ1A, which we designated as CaRfHZ.

2.2. BSA and Genetic Linkage Mapping of the CaRfHZ Gene

In total, 226 plants of the (HZ1A × HZ1C) F2 population, which had been planted in the autumn of 2020, were used to map the CaRfHZ gene according to BSA. DNA from 10 male-fertile plants was mixed to composite the male-fertile DNA bulk, and DNA from 10 male-sterile plants was mixed to composite the male-sterile DNA bulk. In total, 449 pairs of SSR markers [36] evenly distributed on pepper chromosomes were selected to detect the parents: CMS line HZ1A and restorer line HZ1C. Polymorphic markers with clear bands were detected in male-fertile and male-sterile DNA bulks. Finally, the SSR marker P06g8490 was found to be polymorphic among the DNA bulks. This implied that the CaRfHZ gene is close to P06g8490 and is located on chromosome 6. In the 12 Mb region near P06g8490, 11 polymorphic SSR markers among the parents were selected. The above 11 SSR markers and P06g8490 were amplified in 45 male-sterile plants from the (HZ1A × HZ1C) F2 population. According to the recombinant number and genetic distance between the CaRfHZ and its linked markers (Table 3), CaRfHZ was initially located between the markers P06G8229 and P06G8560 (Figure 1a). In this region, 6 SSR markers, P06g8264, P06g8490, P06g8494, P06g8497, P06g8536, and P06g8527, were co-segregated with CaRfHZ.

2.3. Fine Mapping of the CaRfHZ Gene

A larger population composed of 336 male-sterile plants out of the (HZ1A × HZ1C) F2 populations, which had been planted in the spring of 2021, was used to narrow down the region of the CaRfHZ gene. More SSR markers and newly developed InDel markers (Table 3) in the primary mapping region were selected for fine mapping. Twenty-eight markers were polymorphic in HZ1A and HZ1C, and were used to map the CaRfHZ gene. The recombinant numbers and genetic distances between the CaRfHZ gene and its linked markers are shown in Table 3. Some of these markers identified the same number of recombinant individuals, such as P06gInDel-46, P06gInDel-48, and P06gInDel-56, and these markers had the same genetic distance from the CaRfHZ gene. However, owing to image size limitations, the markers are not listed in the mapping linkage (Figure 1a,b). The figure shows five F2 recombinant individuals between the P06G8508 and P06gInDel-95 markers (Figure 1c). Finally, the CaRfHZ gene was mapped between the markers P06gInDel-66 and P06gInDel-89 with the same genetic distance of 0.15 cM (Figure 1b,c). Three markers, P06G8527, P06gInDel-79, and P06gInDel-81, were co-segregated with the CaRfHZ gene.

2.4. Analysis of the Annotation Genes

According to the pepper reference genome sequence (Zunla-1 V2.0) [26], the physical distance between the two markers, P06gInDel-66 and P06gInDel-89, is 533.81kb on Chr. 06: 215,097,259..215,631,069. In this region, there are four annotated genes (Table 4): Capana06g002965, Capana06g002967, Capana06g002968, and Capana06g002969 according to Zunla-1 V2.0 gene models.
Full-length genomic DNA sequences of the four annotated genes were sequenced in HZ1A and HZ1C, and the coding sequences were analyzed (Table 5). There were two missense variations at 129 bp and 504 bp and a synonymous variation at 436 bp in Capana06g002965. There was a missense variation of 935 bp in Capana06g002967. There were two missense variations at 20 bp and 467 bp in Capana06g002968. There were two missense variations at 196 bp and 318 bp and a synonymous variation at 144 bp in Capana06g002969. Differences in the above sites lead to differences in the amino acid sequence (Figure S1), but only the protein structure of Capana06g002969 was different. Due to the missense variation at 318 bp, the stop codon TGA of Capana06g002969 in HZ1C was changed to TGG in HZ1A, which resulted in 31 more amino acids in HZ1A.
The expression of the four annotated genes in young leaves and three developmental stages (stage 1 (S1), sporogenous tissue to meiotic stage; stage 2 (S2), tetrad to mononuclear stage; stage 3 (S3), mature pollen stage) of anthers was analyzed in HZ1A, HZ1C, and (HZ1A × HZ1C) F1 using quantitative real-time PCR (qRT-PCR) (Figure 2). The expression of Capana06g002965 in (HZ1A × HZ1C) F1 was low in young leaves but was not detected in HZ1A and HZ1C. At the S1 stage, the expression level of Capana06g002965 in HZ1A was significantly lower than that in HZ1C and (HZ1A × HZ1C) F1. Conversely, the expression level of Capana06g002965 in HZ1A was significantly higher than that in HZ1C and (HZ1A × HZ1C) F1 at the S2 stage. The expression of Capana06g002965 was not detected in HZ1A, HZ1C, or (HZ1A × HZ1C) F1. The expression of Capana06g002967 in HZ1A was not significantly different from that in HZ1C, but was higher than that in (HZ1A × HZ1C) F1 in young leaves and the S3 stage. The expression trends of Capana06g002967 were similar to Capana06g002965 at the S1 and S2 stages. The expression level of Capana06g002968 in HZ1A was significantly lower than that in HZ1C and (HZ1A × HZ1C) F1 in young leaves and the S1 stage. At the S2 and S3 stages, the expression levels of Capana06g002968 in HZ1A were significantly higher than those in HZ1C. The expression of Capana06g002969 was not detected in the young leaves of HZ1A, HZ1C, and (HZ1A × HZ1C) F1. At the S1, S2, and S3 stages, the expression level of Capana06g002969 in HZ1A was significantly lower than that in HZ1C. In summary, in HZ1A and HZ1C, four genes were differentially expressed at the S1, S2, and S3 developmental stages of anther, with the exception of Capana06g002967 at the S3 stage. The gene expression trends in HZ1C and (HZ1A × HZ1C) F1 were mostly the same.
Although there were differences in the amino acid sequences and gene expression levels of the four annotated genes, it was not possible to determine which gene was the candidate gene for CaRfHZ based on these differences.

3. Discussion

3.1. Comparison of CaRfHZ Mapping Interval and the Published Rf Gene Position in Pepper

Thus far, the seven pepper Rf genes have been fine mapped using a forward genetics strategy (Table 1). According to the gene mapping intervals, these seven genes are located in an adjacent region of approximately 7 Mb from 210–217 Mb on chromosome 6, and the mapping results overlapped in some studies [20,31]. The situation in which the mapping results are adjacent is presumed to be related to the similarity of the materials used for mapping or to the possible clustering of pepper Rf genes on chromosome 6. In the present study, the CaRfHZ gene was fine-mapped within a 533.81 kb segment in the adjacent region. The mapping interval of the CaRfHZ gene partially overlaps the mapping region in Rf [20] and is located within the mapping range of Rf (858.26 kb) [31] and NEDD8 (5.11 Mb) [34], which are related to the larger mapping range of the above research results. Similar mapping results proved the accuracy of our results in this study but were not conducive to the discovery of new pepper Rf genes. As pepper Rf genes may exist in clusters on chromosome 6, there may be multiple genes controlling fertility restoration simultaneously. This leads to inconsistent genotypes for the molecular marker Rf, making it difficult to locate the Rf gene within a shorter interval.

3.2. Prediction and Characteristics of the Candidate Gene

According to the annotation information of the Zunla-1 V2.0 genome [26], there are four annotated genes in the localization interval of the CaRfHZ gene. All four annotated genes showed sense variations and differences in expression between HZ1A and HZ1C. The candidate gene for CaRfHZ could not be determined based on the above information.
Therefore, variations between HZ1A and HZ1C were detected in nine maintainer and four restorer lines of HZ1A (Table 5), and missense variations were analyzed. In Capana06g002965, missense variations were the same as HZ1A in eight maintainer lines, but they were different from HZ1C in all four restorer lines, and the correlation between missense variations and restored traits was 61.54% in 13 lines. In Capana06g002967, missense variations were the same as HZ1C in all four restorer lines, but they were different from HZ1A in eight maintainer lines, and the correlation between missense variations and restored traits was 38.46% in 13 lines. The detection results of missense variations in Capana06g002968 and Capana06g002969 were similar; the missense variations were the same as HZ1A in eight maintainer lines and were the same as HZ1C in three restorer lines, and the correlation between missense variations and restored traits was 84.62% in 13 lines. Conclusively, the missense variations in Capana06g002965 and Capana06g002967 were not highly correlated with restored traits, while the missense variations in Capana06g002968 and Capana06g002969 were highly correlated with restored traits and may be candidate genes for CaRfHZ. Capana06g002968 was not reported as a candidate gene for pepper Rf in previous studies, but Capana06g002969 was identified as a candidate gene. There was no difference in Capana06g002969 sequence between the CMS and restorer lines in a study by Wu et al. [31], which is different from the findings in the present study.
Using the simple modular architecture research tool (SMART) (https://smart.embl.de/, accessed on 15 March 2022) [37] to analyze protein structures, it was found that Capana06g002968 in HZ1A and HZ1C both contained tetratricopeptide repeat (TPR) domains, and Capana06g002969 in HZ1A, but not HZ1C, contained an IBR domain. No reports related to plant fertility were found for the IBR domain, but there have been reports on the TPR domain. In rice, silencing the gene ORF3 containing the TPR domain can restore the spikelet fertility of heterozygous plants [38]. The C-terminal TPR and CaMbd domains of the wFKBP73 gene were critical for male fertility in transgenic rice [39]. Overexpression of AtTRP1, which contains a TRP domain, in wild-type Arabidopsis results in dwarf plants and reduced fertility [40]. High expression of TPR may affect normal plant fertility. In our study, anther development of HZ1A started abnormally from the S2 stage, and the expression of Capana06g002968 in HZ1C was lower than that in HZ1A in the S2 and S3 stages. This suggests that lower expression of Capana06g002968 may be involved in the restoration of fertility of HZ1A. Therefore, Capana06g002968 was proposed as the strongest candidate gene for CaRfHZ, but this needs to be validated in the future.

3.3. Application of the Related Markers of CaRfHZ Gene in Pepper Breeding

Since the discovery of the first pepper CMS by Peterson [15], CMS sterile lines have been used for pepper breeding. The ‘three-line’ hybrid seed production system consisting of CMS, maintainer, and restorer lines is used in pepper seed production. This method can reduce the link of manual emasculation, thereby reducing the production cost and improving the purity of the hybrid. Several excellent pepper hybrid varieties were obtained by this method in China [41,42]. Fertility stability of the hybrid is an important guarantee for the large-scale promotion of the ‘three-line’ hybrid, and the Rf gene in the restorer line plays a key role.
Based on the differential variations in Capana06g002968 (467 bp) and Capana06g002969 (318 bp) in HZ1A and HZ1C, we developed two CAPS markers (results not shown). Amplification and enzyme digestion were performed on 13 lines (Table 5), and the results were consistent with the sequencing results. The accuracy of detecting both maintainer and restorer lines was 84.62% (Capana06g002968) and 92.31% (Capana06g002969), but not 100%, which may be related to the existence of other fertility restorer genes. The accuracies of 84.62% and 92.31% were similar to those of previous studies, such as Co1Mod1-CAPS (88.0% and 92.2%, respectively) [28], and CRF3S1S, CRF-SCAR, and Co1Mod1-CAPS markers (77.2%, 79.2%, and 70.3%, respectively) [33]. The newly developed CAPS markers in this study can be used for sterile and restorer line screening, but they need validation in more maintainer and restorer lines.

4. Materials and Methods

4.1. Plant Materials

The CMS line HZ1A is a pepper line with complete sterility that is obtained by backcrossing for more than 15 generations. The F2 segregating populations (Table 2) were derived from 11 F1 hybrid plants, with HZ1A and HZIC as the female and male parents, respectively. In total, 226 F2 plants and 1442 F2 plants were planted in the autumn of 2020 and spring of 2021, respectively. Plants of F2 populations were used for inheritance analysis and mapping of the Rf gene CaRfHZ. All plant materials, consisting of HZ1A, HZ1C, maintainer lines, restorer lines, F1 plants, and F2 plants, were grown in a plastic-covered greenhouse at the Zhuanghang Experimental Station of Shanghai Academy of Agricultural Sciences (Shanghai, China).

4.2. Pollen Fertility Evaluation

Pollen fertility was evaluated visually with at least three blooming flowers and further discriminated using Alexander’s solution staining [43] under microscope, as described previously [35]. Plants that produced abundant pollen grains were considered male-fertile, whereas plants with no pollen dissemination were considered male-sterile. After Alexander’s solution staining, the pollen grains that stained magenta red were classified as male-fertile, and those stained blue–green were classified as male-sterile.

4.3. Inheritance Analysis of Fertility Restoration

Segregation of male fertility/sterility was identified based on the fertility of the plants from the (HZ1A × HZ1C) F2 population. Chi-square goodness-of-fit test was performed using the following equation: χ2 = Σ (O − E)2/E, where O is the observed frequency, and E is the expected frequency, with p = 0.05 as the threshold.

4.4. Nucleic Acid Extraction

Genomic DNA was extracted from the young leaves using a DNA extraction kit (Tiangen, Beijing, China). The anthers of HZ1A, HZ1C, and (HZ1A × HZ1C) F1 at three developmental stages: stage 1 (S1), sporogenous tissue to meiotic stage; stage 2 (S2), tetrad to mononuclear stage; and stage 3 (S3), mature pollen stage, based on an earlier study [35] were isolated from the flower buds and immediately placed into an RNA stabilization solution (Qiagen, Hilden, Germany). Total RNA was isolated from anthers using a Biospin Plant Total RNA Extraction Kit (Bioer Technology, Hangzhou, China). Reverse transcriptase was used to enrich full-length cDNA using a HiScript II One Step RT-PCR Kit (Vazyme, Nanjing, China). According to the manufacturer’s instructions, DNA and RNA were extracted, and reverse transcription was performed.

4.5. SSR and InDel Markers

SSR markers used in the present study were derived from Cheng et al. [36]. SSR markers are numbered according to their position on the chromosome, and chromosome information is added, such as P01g0001. The markers used in the CaRfHZ linkage map are shown in Table S1.
Chromosome 6 sequences of two published pepper genomes, Zunla-1 V2.0 and Chiltepin V2.0, were obtained from the China National GeneBank Sequence Archive (CNSA) (https://db.cngb.org/search/project/CNPhis0000547/, accessed on 25 October 2020). According to the initial and fine mapping of the CaRfHZ gene, 4 Mb (Chr06:195,000,000..196,000,000 and Chr06: 213,000,000..216,000,000) sequences were selected to search for InDel. The InDel variations between ‘Zunla-1’ and ‘Chiltepin’ were detected using MUMmer 3.0 software [44] (Hamburg, Germany), and the parameter of variation length ≥ 5 bp was used for screening. The 500 bp sequences before and after the InDel variation were selected, and Primer-BLAST software (Bethesda, MD, USA) was used for primer sequence design. In total, 102 InDel markers were developed, and the information is listed in Table S2. Sangon Biotech (Shanghai, China) synthesized the primers. The design and synthesis of all primers used in the present study were the same as those described above.

4.6. BSA and Linkage Analysis of CaRfHZ Gene

The CaRfHZ gene was mapped using BSA. DNA from 10 plants each from male-fertile and male-sterile plants of the (HZ1A × HZ1C) F2 population were randomly selected and mixed to composite the male-fertile and male-sterile DNA bulks. Then, 449 SSR markers [36], which covered the entire pepper genome, were tested for polymorphism in the parent HZ1A and HZ1C. Thereafter, the polymorphic markers between parents were detected in male-fertile and male-sterile DNA bulks. The polymorphic marker between DNA bulks was recognized as linked to the CaRfHZ gene.
In the present study, the male-sterile plants from the (HZ1A × HZ1C) F2 populations were used to map the CaRfHZ gene (Table 2) according to Zhang et al. [45]. The polymorphic markers between DNA bulks and the SSR and InDel markers near this marker were amplified and analyzed in male-sterile plants from the (HZ1A × HZ1C) F2 population. The recombination rate (c) between the marker and the CaRfHZ gene was calculated using the formula c = (N1 + N2/2)/N, where N is the number of all sterile plants, N1 is the number of plants with homozygous bands from the fertile parent HZ1C, and N2 is the number of plants with heterozygous bands from the two parents [45,46]. Kosambi’s formula: d = 1/4 log ((1 + 2c)/(1 − 2c)) was used to calculate the genetic distance (d) [47]. The unit of d is centimorgans (cM). MapDraw v2.1 was used to construct the genetic map [48].
PCR amplification was performed using a reaction mixture containing the following: 5 µL of 2×Hieff® PCR Master Mix (YEASEN, Shanghai, China), 1 μL of 10 µM PCR upstream primer, 1 μL of 10 µM PCR downstream primer, 20 ng template DNA, and ddH2O up to 10 µL. The PCR cycling program was as follows: pre-denaturation at 94 °C for 2 min; followed by 32 cycles of 94 °C for 30 s, primer annealing temperature for 30 s, and 72 °C for 30 s; final extension at 72 °C for 5 min; cooling to 12 °C. Amplification products were electrophoresed using 8% non-denaturing polyacrylamide gels at 150 V for 1 h. The gels were visualized using 0.1% AgNO3 solution [46].

4.7. Candidate Gene Amplification and qRT-PCR

Amplification of the candidate genes of CaRfHZ in the HZ1A and HZ1B parental lines was performed following the instructions using GoldenStar® T6 Super PCR Mix Ver.2 (Tsingke Biotechnology, Beijing, China) in a total reaction volume of 50 µL containing 45 μL of 1.1× GoldenStar® Mix Ver.2, 2 μL of 10 µM upstream primer, 2 μL of 10 µM downstream primer, and 200 ng template DNA. The PCR cycling conditions were as follows: pre-denaturation at 98 °C for 2 min; followed by 32 cycles of 98 °C for 10 s, primer annealing temperature for 15 s, and 72 °C for 40 s; final extension at 72 °C for 5 min, and cooling to 12 °C. The primers used for cloning the candidate genes of CaRfHZ are listed in Table S3. PCR products were sequenced by Sangon Biotech (Shanghai, China). The amino acid sequences of the candidate genes of CaRfHZ were aligned using BioXM 2.7.120 (Nanjing, China).
The qRT-PCR was performed using 2x Hieff UNICON® Universal Blue qPCR SYBR Green Master Mix (YEASEN, Shanghai, China) following the instructions. The PCR reactions were analyzed using the Quant Studio 5 real-time PCR system (Waltham, MA, USA). The β-actin gene from Lv et al. [49] was used as an internal control to normalize expression data. The 2−ΔΔCt method was used to calculate the relative expression levels. Three biological and three technical replicates were used for each experiment. The primers used for qRT-PCR are listed in Table S4.

4.8. Analysis of Candidate Genes

The physical positions of the molecular markers were determined with reference to the pepper reference Zunla-1 V2.0 genome sequence as described above. The annotated genes in the identified target region were obtained according to the Zunla-1 V2.0 gene models [26]. The protein structures of annotated genes were predicted by SMART (https://smart.embl.de/, accessed on 15 March 2022) [37].

5. Conclusions

To summarize, the male fertility of the pepper CMS line HZ1A could be restored by a single dominant gene, CaRfHZ. The Rf gene, CaRfHZ, was fine-mapped to a 533.81 kb region on chromosome 6 using BSA. Four annotated genes were analyzed according to Zunla-1 V2.0 gene models. Based on the comparison of genomic DNA sequences, qRT-PCR, and protein structures in HZ1A and its maintainer and restorer lines, Capana06g002968 was proposed as the strongest candidate gene for the CaRfHZ gene. Functional verification of the CaRfHZ gene needs to be performed in the future.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23147633/s1.

Author Contributions

Conceptualization, Z.N., W.Z. and Q.N.; methodology, Y.S., H.W. and J.C.; software, Z.N. and Y.S.; investigation, Z.N. and J.C.; resources, J.C.; data curation, Z.N. and Y.S.; writing—original draft preparation, Z.N.; writing—review and editing, W.Z. and Q.N.; supervision, W.Z. and Q.N.; funding acquisition, Z.N. and W.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Independent Declaration Project of Agricultural and Social Research in Hangzhou, grant number 20191203B58; the Excellent Team Building Program of Shanghai Academy of Agricultural Sciences, grant number 2022019; the National Major Vegetable Industry Technology System Construction Project, grant number CARS-25; and the Science and Technology Project of the Hangzhou Academy of Agricultural Sciences, grant number 2020HNCT-10.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, Y.J.; Zhang, D. Molecular control of male fertility for hybrid crop breeding. Trends Plant Sci. 2018, 23, 53–65. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, Z.; Zou, Y.; Li, X.; Zhang, Q.; Chen, L.; Wu, H.; Su, D.; Chen, Y.; Guo, J.; Luo, D. 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 2006, 18, 676–687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kazama, T.; Toriyama, K. A fertility restorer gene, Rf4, which is widely used for hybrid rice breeding, encodes a pentatricopeptide repeat protein. Rice 2014, 7, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hu, J.; Wang, K.; Huang, W.; Liu, G.; Gao, Y.; Wang, J.; Huang, Q.; Ji, Y.; Qin, X.; Wan, L.; et al. The rice pentatricopeptide repeat protein RF5 restores fertility in Hong-Lian cytoplasmic male-sterile lines via a complex with the glycine-rich protein GRP162. Plant Cell 2012, 24, 109–122. [Google Scholar] [CrossRef] [Green Version]
  5. Huang, W.; Yu, C.; Hu, J.; Wang, L.; Dan, Z.; Zhou, W.; He, C.; Zeng, Y.; Yao, G.; Qi, J. Pentatricopeptide-repeat family protein RF6 functions with hexokinase 6 to rescue cytoplasmic male sterility in rice. Proc. Natl. Acad. Sci. USA 2015, 112, 14984–14989. [Google Scholar] [CrossRef] [Green Version]
  6. Cui, X.; Wise, R.P.; Schnable, P.S. The rf2 nuclear restorer gene of male-sterile T-cytoplasm maize. Science 1996, 272, 1334–1336. [Google Scholar] [CrossRef] [PubMed]
  7. Zabala, G.; Gabay-Laughnan, S.; Laughnan, J.R. The nuclear gene Rf3 affects the expression of the mitochondrial chimeric sequence R implicated in S-type male sterility in maize. Genetics 1997, 147, 847–860. [Google Scholar] [CrossRef]
  8. Yang, G.; Fu, T.; Yi, B.; Hong, D. Mitochondria-targeted PPR protein restores pol cytoplasmic male sterility by reducing orf224 transcript levels in oilseed rape. Mol. Plant 2016, 9, 1082–1084. [Google Scholar]
  9. Uyttewaal, M.; Arnal, N.; Quadrado, M.; Martin-Canadell, A.; Vrielynck, N.; Hiard, S.; Gherbi, H.; Bendahmane, A.; Budar, F.; Mireau, H. Characterization of Raphanus sativus pentatricopeptide repeat proteins encoded by the fertility restorer locus for Ogura cytoplasmic male sterility. Plant Cell 2008, 20, 3331–3345. [Google Scholar] [CrossRef] [Green Version]
  10. Melonek, J.; Duarte, J.; Martin, J.; Beuf, L.; Murigneux, A.; Varenne, P.; Comadran, J.; Specel, S.; Levadoux, S.; Bernath-Levin, K.; et al. The genetic basis of cytoplasmic male sterility and fertility restoration in wheat. Nat. Commun. 2021, 12, 1036. [Google Scholar] [CrossRef]
  11. Wang, T.L.; He, T.T.; Ding, X.L.; Zhang, Q.Q.; Yang, L.S.; Nie, Z.X.; Zhao, T.J.; Gai, J.Y.; Yang, S.P. Confirmation of GmPPR576 as a fertility restorative gene of cytoplasmic male sterility in soybean. J. Exp. Bot. 2021, 72, 7729–7742. [Google Scholar] [CrossRef] [PubMed]
  12. Jaqueth, J.S.; Hou, Z.; Zheng, P.; Ren, R.; Nagel, B.A.; Cutter, G.; Niu, X.; Vollbrecht, E.; Greene, T.W.; Kumpatla, S.P. Fertility restoration of maize CMS-C altered by a single amino acid substitution within the Rf4 bHLH transcription factor. Plant J. 2020, 101, 101–111. [Google Scholar] [CrossRef] [PubMed]
  13. Tang, H.; Luo, D.; Zhou, D.; Zhang, Q.; Tian, D.; Zheng, X.; Chen, L.; Liu, Y.G. The rice restorer Rf4 for wild-abortive cytoplasmic male sterility encodes a mitochondrial-localized PPR protein that functions in reduction of WA352 transcripts. Mol. Plant 2014, 7, 1497–1500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Qin, X.; Tian, S.; Zhang, W.; Zheng, Q.; Wang, H.; Feng, Y.; Lin, Y.; Tang, J.; Wang, Y.; Yan, J.; et al. The main restorer Rf3 of maize S-type cytoplasmic male sterility encodes a PPR protein that functions in the reduction of the transcripts of orf355. Mol. Plant 2021, 14, 1961–1964. [Google Scholar] [CrossRef] [PubMed]
  15. Peterson, P. Cytoplasmically inherited male sterility in capsicum. Am. Nat. 1958, 92, 111–119. [Google Scholar] [CrossRef]
  16. Xu, F.; Yang, X.; Zhao, N.; Hu, Z.; Mackenzie, S.A.; Zhang, M.; Yang, J. Exploiting sterility and fertility variation in cytoplasmic male sterile vegetable crops. Hortic. Res. 2022, 9, uhab039. [Google Scholar] [CrossRef]
  17. Wang, L.H.; Zhang, B.X.; Lefebvre, V.; Huang, S.W.; Daubeze, A.M.; Palloix, A. QTL analysis of fertility restoration in cytoplasmic male sterile pepper. Theor. Appl. Genet. 2004, 109, 1058–1063. [Google Scholar] [CrossRef]
  18. Wei, B.; Wang, L.; Chen, L.; Zhang, R. Genetic analysis on the restoration of cytoplasmic male sterility with mixed model of major gene plus polygene in pepper. Acta Hortic. Sin. 2013, 40, 2263–2268. [Google Scholar]
  19. Gulyas, G.; Pakozdi, K.; Lee, J.S.; Hirata, Y. Analysis of fertility restoration using cytoplasmic male-sterile red pepper (Capsicum annuum L.) lines. Breed. Sci. 2006, 56, 331–334. [Google Scholar] [CrossRef] [Green Version]
  20. Ye, Q.; Ruan, M.; Wang, R.; Yao, Z.; Li, Z.; Wan, H.; Cheng, Y.; Yang, Y.; Zhou, G. Genetic analysis and molecular mapping of the restorer gene for cytoplasmic male sterility in pepper. Mol. Plant Breed. 2017, 15, 4985–4991. [Google Scholar]
  21. Zhang, B.; Huang, S.; Yang, G.; Guo, J. Two RAPD markers linked to a major fertility restoration gene in pepper. Euphytica 2000, 113, 155–161. [Google Scholar]
  22. Kim, D.S.; Kim, D.H.; Yoo, J.H.; Kim, B.D. Cleaved amplified polymorphic sequence and amplified fragment length polymorphism markers linked to the fertility restorer gene in chili pepper (Capsicum annuum L.). Mol. Cells 2006, 21, 135–140. [Google Scholar]
  23. Lee, J.; Yoon, J.B.; Park, H.G. A CAPS marker associated with the partial restoration of cytoplasmic male sterility in chili pepper (Capsicum annuum L.). Mol. Breed. 2008, 21, 95–104. [Google Scholar] [CrossRef]
  24. Yang, J.; Wang, W.; Shen, H. SSR marker mapping and marker-assisted selection of Rf gene in pepper. Chin. Cucurbits Veg. 2010, 23, 1–5. [Google Scholar]
  25. Jo, Y.D.; Kim, Y.M.; Park, M.N.; Yoo, J.H.; Park, M.; Kim, B.D.; Kang, B.C. Development and evaluation of broadly applicable markers for restorer-of-fertility in pepper. Mol. Breed. 2010, 25, 187–201. [Google Scholar] [CrossRef]
  26. Qin, C.; Yu, C.; Shen, Y.; Fang, X.; Chen, L.; Min, J.; Cheng, J.; Zhao, S.; Xu, M.; Luo, Y.; et al. Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization. Proc. Natl. Acad. Sci. USA 2014, 111, 5135–5140. [Google Scholar] [CrossRef] [Green Version]
  27. Kim, S.; Park, M.; Yeom, S.I.; Kim, Y.M.; Lee, J.M.; Lee, H.A.; Seo, E.; Choi, J.; Cheong, K.; Kim, K.T.; et al. Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nat. Genet. 2014, 46, 270–278. [Google Scholar] [CrossRef]
  28. Jo, Y.D.; Ha, Y.; Lee, J.H.; Park, M.; Bergsma, A.C.; Choi, H.I.; Goritschnig, S.; Kloosterman, B.; van Dijk, P.J.; Choi, D.; et al. Fine mapping of Restorer-of-fertility in pepper (Capsicum annuum L.) identified a candidate gene encoding a pentatricopeptide repeat (PPR)-containing protein. Theor. Appl. Genet. 2016, 129, 2003–2017. [Google Scholar] [CrossRef]
  29. Zhang, Z.H.; An, D.L.; Cao, Y.C.; Yu, H.L.; Zhu, Y.S.; Mei, Y.J.; Zhang, B.X.; Wang, L.H. Development and application of KASP markers associated with Restorer-of-fertility gene in Capsicum annuum L. Physiol. Mol. Biol. Plant 2021, 27, 2757–2765. [Google Scholar] [CrossRef]
  30. Kang, M.C.; Kang, H.J.; Jung, S.Y.; Lee, H.Y.; Kang, M.Y.; Jo, Y.D.; Kang, B.C. The Unstable Restorer-of-fertility locus in pepper (Capsicum annuum L.) is delimited to a genomic region containing PPR genes. Theor. Appl. Genet. 2022, 135, 1923–1937. [Google Scholar] [CrossRef]
  31. Wu, L.; Wang, P.; Wang, Y.; Cheng, Q.; Lu, Q.; Liu, J.; Li, T.; Ai, Y.; Yang, W.; Sun, L.; et al. Genome-wide correlation of 36 agronomic traits in the 287 pepper (Capsicum) accessions obtained from the SLAF-seq-Based GWAS. Int. J. Mol. Sci. 2019, 20, 5675. [Google Scholar] [CrossRef] [Green Version]
  32. Cheng, J.; Chen, Y.; Hu, Y.; Zhou, Z.; Hu, F.; Dong, J.; Chen, W.; Cui, J.; Wu, Z.; Hu, K. Fine mapping of restorer-of-fertility gene based on high-density genetic mapping and collinearity analysis in pepper (Capsicum annuum L.). Theor. Appl. Genet. 2020, 133, 889–902. [Google Scholar] [CrossRef]
  33. Zhang, Z.; Zhu, Y.; Cao, Y.; Yu, H.; Bai, R.; Zhao, H.; Zhang, B.; Wang, L. Fine mapping of the male fertility restoration gene CaRf032 in Capsicum annuum L. Theor. Appl. Genet. 2020, 133, 1177–1187. [Google Scholar] [CrossRef] [PubMed]
  34. Wei, B.; Bosland, P.W.; Zhang, Z.; Wang, Y.; Zhang, G.; Wang, L.; Yu, J. A predicted NEDD8 conjugating enzyme gene identified as a Capsicum candidate Rf gene using bulk segregant RNA sequencing. Hortic. Res. 2020, 7, 210. [Google Scholar] [CrossRef] [PubMed]
  35. Nie, Z.; Chen, J.; Song, Y.; Fu, H.; Wang, H.; Niu, Q.; Zhu, W. Comparative transcriptome analysis of the anthers from the cytoplasmic male-sterile pepper line HZ1A and its maintainer line HZ1B. Horticulturae 2021, 7, 580. [Google Scholar] [CrossRef]
  36. Cheng, J.; Zhao, Z.; Li, B.; Qin, C.; Wu, Z.; Trejo-Saavedra, D.L.; Luo, X.; Cui, J.; Rivera-Bustamante, R.F.; Li, S.; et al. A comprehensive characterization of simple sequence repeats in pepper genomes provides valuable resources for marker development in Capsicum. Sci. Rep. 2016, 6, 18919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Letunic, I.; Khedkar, S.; Bork, P. SMART: Recent updates, new developments and status in 2020. Nucleic Acids Res. 2021, 49, D458–D460. [Google Scholar] [CrossRef]
  38. Yu, Y.; Zhao, Z.; Shi, Y.; Tian, H.; Liu, L.; Bian, X.; Xu, Y.; Zheng, X.; Gan, L.; Shen, Y.; et al. Hybrid sterility in rice (Oryza sativa L.) involves the tetratricopeptide repeat domain containing protein. Genetics 2016, 203, 1439–1451. [Google Scholar] [CrossRef] [Green Version]
  39. Kurek, I.; Dulberger, R.; Azem, A.; Tzvi, B.B.; Sudhakar, D.; Christou, P.; Breiman, A. Deletion of the C-terminal 138 amino acids of the wheat FKBP73 abrogates calmodulin binding, dimerization and male fertility in transgenic rice. Plant Mol. Biol. 2002, 48, 369–381. [Google Scholar] [CrossRef]
  40. Lin, Z.; Ho, C.W.; Grierson, D. AtTRP1 encodes a novel TPR protein that interacts with the ethylene receptor ERS1 and modulates development in Arabidopsis. J. Exp. Bot. 2009, 60, 3697–3714. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Wang, L.; Zhang, B.; Zhang, Z.; Cao, Y.; Yu, H.; Feng, X. Status in breeding and production of Capsicum spp. in China during ‘The Thirteenth Five-year Plan’ period and future prospect. China Veg. 2021, 2, 21–29. [Google Scholar]
  42. Zhou, X.; Ma, Y.; Dai, X.; Li, X.; Yang, S. Spread and industry development of pepper in China. Acta Hortic. Sinica 2020, 47, 1715–1726. [Google Scholar]
  43. Peterson, R.; Slovin, J.P.; Chen, C. A simplified method for differential staining of aborted and non-aborted pollen grains. J. Plant Biol. 2010, 1, e13. [Google Scholar] [CrossRef] [Green Version]
  44. Kurtz, S.; Phillippy, A.; Delcher, A.L.; Smoot, M.; Shumway, M.; Antonescu, C.; Salzberg, S.L. Versatile and open software for comparing large genomes. Genome Biol. 2004, 5, R12. [Google Scholar] [CrossRef] [Green Version]
  45. Zhang, Q.; Shen, B.Z.; Dai, X.K.; Mei, M.H.; Maroof, M.A.S.; Li, Z.B. Using bulked extremes and recessive class to map genes for photoperiod-sensitive genic male sterility in rice. Proc. Natl. Acad. Sci. USA 1994, 91, 8675–8679. [Google Scholar] [CrossRef] [Green Version]
  46. Nie, Z.; Zhao, T.; Liu, M.; Dai, J.; He, T.; Lyu, D.; Zhao, J.; Yang, S.; Gai, J. Molecular mapping of a novel male-sterile gene msNJ in soybean [Glycine max (L.) Merr.]. Plant Reprod. 2019, 32, 371–380. [Google Scholar] [CrossRef]
  47. Kosambi, D.D. The estimation of map distance from recombination values. Ann. Eugen. 1943, 12, 172–175. [Google Scholar] [CrossRef]
  48. Liu, R.H.; Meng, J.L. MapDraw: A Microsoft excel macro for drawing genetic linkage maps based on given genetic linkage data. Hereditas 2003, 25, 317–321. [Google Scholar]
  49. Lv, J.H.; Liu, Z.B.; Liu, Y.H.; Ou, L.J.; Deng, M.H.; Wang, J.; Song, J.S.; Ma, Y.Q.; Chen, W.C.; Zhang, Z.Q.; et al. Comparative transcriptome analysis between cytoplasmic male-sterile line and its maintainer during the floral bud development of pepper. Hortic. Plant J. 2020, 6, 89–98. [Google Scholar] [CrossRef]
Ijms 23 07633 g001 550
Figure 1. Mapping of the pepper CaRfHZ gene. (a) Primary mapping of CaRfHZ using a small population of (HZ1A × HZ1C) F2. (b) Fine mapping of CaRfHZ using a large population of (HZ1A × HZ1C) F2. (c) Genotyping results from 5 F2 recombinants by 2 SSR markers and 12 newly developed InDel markers: ‘A’ = ‘HZ1A’ genotype; ‘B’ = ‘HZ1C’ genotype; ‘H’ = (HZ1A × HZ1C) F1 genotype. Distances are shown in cM. Due to image size limitations, the markers that co-segregated with the fertility phenotype are not all listed in the mapping linkage (a,b).
Figure 1. Mapping of the pepper CaRfHZ gene. (a) Primary mapping of CaRfHZ using a small population of (HZ1A × HZ1C) F2. (b) Fine mapping of CaRfHZ using a large population of (HZ1A × HZ1C) F2. (c) Genotyping results from 5 F2 recombinants by 2 SSR markers and 12 newly developed InDel markers: ‘A’ = ‘HZ1A’ genotype; ‘B’ = ‘HZ1C’ genotype; ‘H’ = (HZ1A × HZ1C) F1 genotype. Distances are shown in cM. Due to image size limitations, the markers that co-segregated with the fertility phenotype are not all listed in the mapping linkage (a,b).
Ijms 23 07633 g001
Ijms 23 07633 g002 550
Figure 2. Expression levels of four annotated genes in young leaves and in anthers. * and ** indicate significant difference at p < 0.05 and p < 0.01, respectively.
Figure 2. Expression levels of four annotated genes in young leaves and in anthers. * and ** indicate significant difference at p < 0.05 and p < 0.01, respectively.
Ijms 23 07633 g002
Table
Table 1. Reported Rf candidate genes on chromosome 6 in pepper.
Table 1. Reported Rf candidate genes on chromosome 6 in pepper.
Gene NameMapping MethodRangeFlanking MarkersMapping Interval
(Zunla-1 V2.0)		Candidate GenesLocations of Candidate Genes
(Zunla-1 V2.0)		References
RfCAPS, SSR, SCAR markers and BSA498.60 kbpep20~pep43Chr06: 214,671,854..
215,170,454		Capana06g002962
Capana06g002963
Capana06g002965		Chr06: 214,791,005..214,791,707
Chr06: 214,928,883..214,932,088
Chr06: 215,102,194..215,102,754		[20]
CaPPR6BSA–AFLP and comparative mapping128.96 kbS423~S424Nearby Chr06: 216,988,988CaPPR6Chr06: 214,033,734..214,035,560[28,29]
RfuCAPS, SCAR markers398.00 kb4162-SCAR
~G16-CAPS		Nearby Chr06: 214,076,189CA00g30080-[30]
RfSLAF-seq, and GWAS858.26 kbP06-247
~PW6-126		Chr06: 214,868,888..
215,727,145		Capana06g002967
Capana06g002969		Chr06: 215,172,352..215,173,332
Chr06: 215,340,394..215,340,521		[31]
CaRfConjoint analysis of recombinants and collinearity270.10 kbS3~S1Nearby Chr06: 217,065,420Capana06g003028Chr06: 217,273,835..217,280,134[32]
CaRf032Genome resequencing and recombination analysis148.05 kbS1402~S1354Chr06: 213,923,525..
214,071,576		CA00g82510Chr06: 214,033,838..214,035,589[33]
NEDD8Bulked segregant RNA sequencing5.1 MbKS18~KS22Chr06: 210,576,870..
215,685,280		Capana06g002866Chr06: 210,452,333..210,452,887[34]
“-” indicates there is no corresponding region. In some studies, CM334 v1.55 was selected as the reference genome; therefore, the flanking markers of the mapping interval were compared on the Zunla-1 V2.0 genome to obtain the corresponding positions. For the markers without corresponding positions, according to the corresponding positions of the annotated genes in the mapping interval on the Zunla-1 V2.0 genome, the approximate locations of the mapping intervals, such as CaPPR6 [28,29], Rfu [30], and CaRf [32] were speculated. The mapping interval of NEDD8 [34] in the original text is inconsistent with the candidate gene position.
Table
Table 2. Segregation patterns of male fertility and sterility in the pepper F2 population derived from the cross HZ1A × HZ1C.
Table 2. Segregation patterns of male fertility and sterility in the pepper F2 population derived from the cross HZ1A × HZ1C.
Year and SeasonLine NameNo. of PlantsExpected Ratioχ2Probability
TotalFertileSterile
2020 AutumnP202106653133:10.730.39
P202176353102.330.13
P202189765220.000.95
2021 SpringP211014134102323:10.040.84
P211017370279910.010.90
P211091129102270.930.33
P211103168134341.790.18
P211115160128321.880.17
P211129161121400.000.96
P211157218165530.020.88
P21116210275270.050.82
Total 165812773813:13.500.06
Table
Table 3. The recombinant number and genetic distance between the CaRfHZ gene and its linked markers.
Table 3. The recombinant number and genetic distance between the CaRfHZ gene and its linked markers.
MarkerPosition (bp)Recombinant IndividualsGenetic Distance (cM)
Primary mappingP06g8405203,656,182 3 HO6.71
P06g8077212,732,221 2 HO, 1 HE5.58
P06g8089204,569,609 2 HO4.46
P06g8229208,899,271 2 HE2.22
CaRfHZ
P06g8264209,767,667 00.00
P06g8490214,353,887 00.00
P06g8494214,428,701 00.00
P06g8497214,494,826 00.00
P06g8536215,672,797 00.00
P06g8527215,336,423 00.00
P06g8560215,995,316 1 HE1.11
P06g8618216,556,879 2 HO4.46
Fine mappingP06g8229208,899,271 43 HE6.43
P06g8335210,959,629 35 HE5.23
P06g8405212,732,234 24 HE3.58
P06g8434213,257,869 23 HE3.43
P06g8219208,736,376 22 HE3.28
P06gInDel-15213,906,630 20 HE2.98
P06g8528215,338,579 6 HE0.89
P06gInDel-38214,431,109 4 HE0.60
P06g8490214,353,887 4 HE0.60
P06g8499214,498,896 4 HE0.60
P06g8508214,872,305 3 HE0.45
P06gInDel-46214,885,975 2 HE0.30
P06gInDel-48214,903,798 2 HE0.30
P06gInDel-56214,950,959 2 HE0.30
P06gInDel-66215,097,259 1 HE0.15
CaRfHZ
P06g8527215,336,423 00.00
P06gInDel-79215,419,831 00.00
P06gInDel-81215,498,386 00.00
P06gInDel-89215,631,069 1 HE0.15
P06gInDel-90215,633,948 1 HE0.15
P06gInDel-91215,636,8892 HE0.30
P06gInDel-92215,669,0322 HE0.30
P06gInDel-94215,669,384 2 HE0.30
P06gInDel-95215,674,072 2 HE0.30
P06g8543215,710,422 3 HE0.45
P06gInDel-99215,741,617 3 HE0.45
P06g8549215,773,438 3 HE0.45
P06g8560215,995,331 6 HE0.89
‘HO’ indicates a recombinant individual with a homozygous band from the fertile parent HZ1C, and ‘HE’ indicates a recombinant individual with heterozygous bands from the two parents.
Table
Table 4. List of four annotation genes in the mapping region of CaRfHZ.
Table 4. List of four annotation genes in the mapping region of CaRfHZ.
Gene IDPosition (bp)Annotation Description
Capana06g002965Chr06:215,102,194..215,102,754Unknown protein
Capana06g002967Chr06:215,172,352..215,173,332CW-type zinc finger protein
Capana06g002968Chr06:215,328,827..215,334,970Tetratricopeptide repeat protein
Capana06g002969Chr06:215,340,394..215,341,013Unknown protein
Table
Table 5. The variations of four annotated genes in the CMS, maintainer, and restorer lines.
Table 5. The variations of four annotated genes in the CMS, maintainer, and restorer lines.
Line nameTypeCapana06g002965Capana06g002967Capana06g002968Capana06g002969
129 bp436 bp504 bp935 bp20 bp467 bp144 bp196 bp318 bp
HZ1ACMS lineATGGGTGTG
HZ1CRestorer lineCCAACGAGA
P21204Maintainer lineCTGGGGGGG
P21238Maintainer lineATGAGTATG
P21239Maintainer lineATGAGTATG
P21241Maintainer lineATGAGTATG
P21243Maintainer lineATGAGTATG
P21244Maintainer lineATGAGTATG
P21246Maintainer lineATGAGTATG
P21270Maintainer lineATGAGTATG
P21273Maintainer lineATGAGTATG
P21215Restorer lineATGACGAGA
P21240Restorer lineATGACGAGA
P21247Restorer lineATGACGAGA
P21274Restorer lineATGAGTATG
Shading indicates the nucleobases are consistent with the restored traits.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nie, Z.; Song, Y.; Wang, H.; Chen, J.; Niu, Q.; Zhu, W. Fine Mapping and Gene Analysis of restorer-of-fertility Gene CaRfHZ in Pepper (Capsicum annuum L.). Int. J. Mol. Sci. 2022, 23, 7633. https://doi.org/10.3390/ijms23147633

AMA Style

Nie Z, Song Y, Wang H, Chen J, Niu Q, Zhu W. Fine Mapping and Gene Analysis of restorer-of-fertility Gene CaRfHZ in Pepper (Capsicum annuum L.). International Journal of Molecular Sciences. 2022; 23(14):7633. https://doi.org/10.3390/ijms23147633

Chicago/Turabian Style

Nie, Zhixing, Yunpeng Song, Hong Wang, Jianying Chen, Qingliang Niu, and Weimin Zhu. 2022. "Fine Mapping and Gene Analysis of restorer-of-fertility Gene CaRfHZ in Pepper (Capsicum annuum L.)" International Journal of Molecular Sciences 23, no. 14: 7633. https://doi.org/10.3390/ijms23147633

APA Style

Nie, Z., Song, Y., Wang, H., Chen, J., Niu, Q., & Zhu, W. (2022). Fine Mapping and Gene Analysis of restorer-of-fertility Gene CaRfHZ in Pepper (Capsicum annuum L.). International Journal of Molecular Sciences, 23(14), 7633. https://doi.org/10.3390/ijms23147633

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