Fine mapping of restorer-of-fertility gene based on high-density genetic mapping and collinearity analysis in pepper (Capsicum annuum L.)

Characterization and Expression Analysis of Sugar Transporters through Partial Least Square Structural Equation Model (PLS-SEM) Revealed Their Role in Pepper (Capsicum annuum L.)

Pepper (Capsicum annuum L.) is one of the most important economic crops in the world. By controlling the transport and distribution of photosynthetic products between cells and organs, sugar transporters are widely involved in growth and development, environmental adaptation, and microbial interactions. The present study was carried out at the genome-wide level to systematically characterize sugar transporters. As a result, 50 MST, 3 SUT, and 29 SWEET genes were identified and classified. The expression pattern of sugar transporters in pepper was analyzed by transcriptomic data. The expression properties of sugar transporters were further explored in pepper varieties with significant differences in weight, shape, and pungency. It was shown that the pepper sugar transporter genes had obvious spatiotemporal specific expression characteristics and exhibited variety-specific expression preferences. We focus on analyzing candidate genes that may be involved in fruit development and expansion. We further explore the response of pepper sugar transporters to adversity stress using a structural equation model. Finally, we found that the MST, SUT, and SWEET families are collectively involved in balancing pepper resistance to abiotic stress by coordinating the expression strengths of different family members. Our study may contribute to the functional study of pepper sugar transporter genes and create the prospect of utilizing sugar transporter gene resources to improve pepper variety.

A 163-bp insertion in the Capana10g000198 encoding a MYB transcription factor causes male sterility in pepper (Capsicum annuum L.)

Male sterility provides an efficient approach for commercial exploitation of heterosis. Despite more than 20 genic male sterile (GMS) mutants documented in pepper (Capsicum annuum L.), only two causal genes have been successfully identified. Here, a novel spontaneous recessive GMS mutant, designated msc-3, is identified and characterized at both phenotypic and histological levels. Pollen abortion of msc-3 mutant may be due to the delayed tapetum degradation, leading to the non-degeneration of tetrads callosic wall. Then, a modified MutMap method and molecular marker linkage analysis were employed to fine mapping the msc-3 locus, which was delimited to the ~139.91-kb region harboring 10 annotated genes. Gene expression and structure variation analyses indicate the Capana10g000198, encoding a R2R3-MYB transcription factor, is the best candidate gene for the msc-3 locus. Expression profiling analysis shows the Capana10g000198 is an anther-specific gene, and a 163-bp insertion in the Capana10g000198 is highly correlated with the male sterile (MS) phenotype. Additionally, downregulation of Capana10g000198 in male fertile plants through virus-induced gene silencing resulted in male sterility. Finally, possible regulatory relationships of the msc-3 gene with the other two reported pepper GMS genes, msc-1 and msc-2, have been studied, and comparative transcriptome analysis reveals the expression of 16 GMS homologs are significantly downregulated in the MS anthers. Overall, our results reveal that Capana10g000198 is the causal gene underlying the msc-3 locus, providing important theoretical clues and basis for further in-depth study on the regulatory mechanisms of pollen development in pepper.

Identification and bioinformatic analysis of the CaCesA/Csls family members and the expression of the CaCslD1 in the flower buds of CMS/Rf system in pepper

The cellulose synthase gene superfamily contains cellulose synthase (CesA) and cellulose synthase-like (Csl) gene families, which synthesize cellulose and hemicellulose in plant cell walls and play a crucial role in plant growth and development. However, the CesA/Csl gene family has not been reported in pepper. Therefore, the genome-wide research of the CaCesA/CaCsl gene family was conducted in pepper. In this study, a total of 39 CaCesA/CaCsls genes (10 CesAs genes and 29 Csls genes) were identified in pepper and unevenly distributed on 11 chromosomes. These CaCesA/Csls were divided into seven subfamilies (CesAs, CslAs, CslBs, CslCs, CslDs, CslEs, CslGs), and most of CaCesA/Csls genes are closely related to AtCesA/Csls genes. The cis-acting elements in the promoters of CaCesA/Csls genes are mainly related to hormone response and stress response. There are ten collinear gene pairs between the CesA/Csls gene family of pepper and Arabidopsis, and four fragment duplication gene pairs of the CaCesA/Csls genes were discovered. RNA-seq analysis shows that the majority of CaCesA/Csls are expressed in a variety of plant tissues, indicating that most CaCesA/Csls gene expression patterns are not organ-specific, and CaCslD1/D4 have the highest expression in anthers, followed by petal, ovary, and F9. RNA-seq analysis shows that most CaCesA/Csls are responsive to five hormones (IAA, GA₃, ABA, SA, and MeJA). The tissue-specific expression analysis of the CaCslD1 gene shows that the CaCslD1 gene is expressed specifically in flowers. In the flower buds IV of cytoplasmic male sterility (CMS) and its restoration of fertility (Rf) system, CaCslD1 reach the highest expression respectively. However, the relative expression level of CaCslD1 in the fertile accessions is extremely significantly higher than in the sterile accessions. This study shows an overall understanding of the CaCesA/Csls gene family and provides a new insight for understanding the function of CaCslD1 in pollen development and exploring the fertility restoration of CMS in pepper.

Fine mapping of Rf2, a minor Restorer-of-fertility (Rf) gene for cytoplasmic male sterility in chili pepper G164 (Capsicum annuum L.)

Genome re-sequencing and recombination analyses identified Capana06g000193 as a strong candidate for the minor male fertility restoration locus Rf2 in chili pepper G164 harboring two dominant male fertility restoration genes. Male fertility restoration genes of chili pepper restorer line G164 (Capsicum annuum L.) were studied using molecular marker genotypes of an F₂ population (7G) of G164 crossed with the cytoplasmic male sterility line 77013A. The ratio of sterile to fertile single plants in the F₂ population was 1:15. This result indicates that chili pepper G164 has two dominant restoration genes, which we designated as Rf1 and Rf2. An individual plant recessive for Rf1 and heterozygous for Rf2, 7G-112 (rf1rf1Rf2rf2), was identified by molecular marker selection and genetic analysis, and a single Rf2 gene-segregating population with a 3:1 ratio of fertile to sterile plants was developed from the self-pollination of male fertile individuals of 77013A and 7G-112 hybrid progeny. Bulk segregant analysis of fertile and sterile pools from the segregating populations was used to genetically map Rf2 to a 3.1-Mb region on chromosome 6. Rf2 was further narrowed to a 179.3-kb interval through recombination analysis of molecular markers and obtained the most likely candidate gene, Capana06g000193.

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

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) F₂ populations and map the Rf gene. A single dominant gene CaRfHZ conferred male fertility according to inheritance analysis. Using sterile plants from (HZ1A × HZ1C) F₂ populations and bulked segregant analysis (BSA), the CaRf_HZ 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 CaRf_HZ gene. Our results may help with hybrid pepper breeding and to elucidate the mechanism of male fertility restoration in peppers.

The Unstable Restorer-of-fertility locus in pepper (Capsicum annuum. L) is delimited to a genomic region containing PPR genes

Unstable Restorer-of-fertility (Rf_u), conferring unstable fertility restoration in the pepper CGMS system, was delimited to a genomic region near Rf and is syntenic to the PPR-like gene-rich region in tomato. The use of cytoplasmic-genic male sterility (CGMS) systems greatly increases the efficiency of hybrid seed production. Although marker development and candidate gene isolation have been performed for the Restorer-of-fertility (Rf) gene in pepper (Capsicum annuum L.), the broad use of CGMS systems has been hampered by the instability of fertility restoration among pepper accessions, especially sweet peppers, due to the widespread presence of the Unstable Restorer-of-fertility (Rf_u) locus. Therefore, to investigate the genetic factors controlling unstable fertility restoration in sweet peppers, we developed a segregation population (BC₄F₅) from crosses using a male-sterile line and an Rf_u-containing line. Segregation did not significantly deviate from a 3:1 ratio for unstable fertility restoration to sterility, indicating single dominant locus control for unstable fertility restoration in this population. Genetic mapping delimited the Rf_u locus to a 398 kb genomic region on chromosome 6, which is close to but different from the previously identified Rf-containing region. The Rf_u-containing region harbors a pentatricopeptide repeat (PPR) gene, along with 10 other candidate genes. In addition, this region is syntenic to the genomic region containing the largest number of Rf-like PPR genes in tomato. Therefore, the dynamic evolution of PPR genes might be responsible for both the restoration and instability of fertility in pepper. During genetic mapping, we developed various molecular markers, including one that co-segregated with Rf_u. These markers showed higher accuracy for genotyping than previously developed markers, pointing to their possible use in marker-assisted breeding of sweet peppers.

Genome-Wide Analysis of the GW2-Like Genes in Gossypium and Functional Characterization of the Seed Size Effect of GhGW2-2D

Cotton is one of the most economically important crops worldwide. Seed size is a vital trait for plants connected with yield and germination. GW2 encodes a RING_Ubox E3 ubiquitin ligase that controls seed development by affecting cell growth. Here, are few reports on GW2-like genes in cotton, and the function of GW2 in cotton is poorly understood. In the present study, a genome-wide analysis identified 6 and 3 GW2-like genes in each of the two cultivated tetraploids (Gossypium hirsutum and G. barbadense) and each of their diploid ancestral species (G. arboreum, G. raimondii), respectively. GhGW2-2D has the same functional domain and high sequence similarity with AtDA2 in Arabidopsis. Overexpression of GhGW2-2D in Arabidopsis significantly reduced seed and seedling size, suggesting GhGW2-2D is a potential target for regulating cotton seed size. These results provided information on the genetic and molecular basis of GW2-like genes in cotton, thus establishing a foundation for functional studies of cotton seeds.

Comparative Transcriptome Analysis of the Anthers from the Cytoplasmic Male-Sterile Pepper Line HZ1A and Its Maintainer Line HZ1B

Cytoplasmic male-sterility (CMS) is important for the utilization of crop heterosis and study of the molecular mechanisms involved in CMS could improve breeding programs. In the present study, anthers of the pepper CMS line HZ1A and its maintainer line HZ1B were collected from stages S1, S2, and S3 for transcriptome sequencing. A total of 47.95 million clean reads were obtained, and the reads were assembled into 31,603 unigenes. We obtained 42 (27 up-regulated and 15 down-regulated), 691 (346 up-regulated and 345 down-regulated), and 709 (281 up-regulated and 428 down-regulated) differentially expressed genes (DEGs) in stages S1, S2, and S3, respectively. Through Gene Ontology (GO) analysis, the DEGs were found to be composed of 46 functional groups. Two GO terms involved in photosynthesis, photosynthesis (GO:0015986) and photosystem I (GO:0009522), may be related to CMS. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, oxidative phosphorylation (ko00190) and phenylpropanoid biosynthesis (ko00940) were significantly enriched in the S1 and S2 stages, respectively. Through the analysis of 104 lipid metabolism-related DEGs, four significantly enriched KEGG pathways may help to regulate male sterility during anther development. The mitochondrial genes orf470 and atp6 were identified as candidate genes of male sterility for the CMS line HZ1A. Overall, the results will provide insights into the molecular mechanisms of pepper CMS.

Development and application of KASP markers associated with Restorer-of-fertility gene in Capsicum annuum L

Fertility restoration of cytoplasmic male sterility (CMS) in Capsicum annuum is controlled by multiple alleles of Restorer-of-fertility (Rf) genes. The isolation of additional Rf genes should therefore enrich the knowledge of CMS/Rf systems and accelerate their exploitation in hybrid seed production. In this study, the fertility restorer gene CaRfm of '0601 M', a non-pungent bell pepper, was genetically mapped to a 1.2-cM region flanked by KASP markers S761 and S183. CaRfm was then physically mapped to a 128.96-Kb interval predicted from 24 recombinants with two co-segregated markers, S423 and S424. CaPPR6 encoding a pentatricopeptide repeat (PPR) protein was suggested as the most likely candidate gene for the CaRfm locus on the basis of sequence alignment as well as genotyping of tightly linked markers. In addition, molecular markers S1597 and S1609, which are immediately adjacent to CaRfm at 15.7 and 57.8-Kb respectively, were developed and applied to marker-assisted selection. The results provided friendly markers for breeding pepper restorer lines and laid the foundation for elucidating the male fertility restoration mechanism.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01109-9.

Fine mapping and candidate gene analysis of the up locus determining fruit orientation in pepper (Capsicum spp.)

The up locus determining fruit orientation was fine-mapped into a region with a physical length of ~169.51 kb on chromosome P12 in pepper. Capana12g000958, encoding a developmentally regulated G protein 2, was proposed as the strongest candidate via sequence comparison and expression analysis. Fruit orientation is an important horticultural and domesticated trait, which is controlled by a single semi-dominant gene (up) in pepper. However, the gene underlying up locus has not yet been identified. In this study, the previously detected major QTL UP12.1 was firstly verified using a backcross population (n = 225) stem from the cross of BB3 (C. annuum) and its wild relative Chiltepin (C. annuum var. glabriusculum) using BB3 as the recurrent parent. Then, a large BC₁F₂ population (n = 1827) was used for recombinant screening to delimit the up locus into an interval with ~ 169.51 kb in length. Sequence comparison and expression analysis suggested that Capana12g000958, encoding a developmentally regulated G protein 2, was the most likely candidate gene for the up locus. There is no difference within the coding sequences of Capana12g000958 between BB3 and Chiltepin, while a SNP in the upstream of Capana12g000958 showed a complete correlation with the fruit orientation among a panel of 40 diverse pepper inbred lines. These findings will form a basis for gene isolation and reveal of genetic mechanism underlying the fruit orientation domestication in pepper.

Genome-wide identification and analysis of highly specific CRISPR/Cas9 editing sites in pepper (Capsicum annuum L.)

The CRISPR/Cas9 system is an efficient genome editing tool that possesses the outstanding advantages of simplicity and high efficiency. Genome-wide identification and specificity analysis of editing sites is an effective approach for mitigating the risk of off-target effects of CRISPR/Cas9 and has been applied in several plant species but has not yet been reported in pepper. In present study, we first identified genome-wide CRISPR/Cas9 editing sites based on the ‘Zunla-1’ reference genome and then evaluated the specificity of CRISPR/Cas9 editing sites through whole-genome alignment. Results showed that a total of 603,202,314 CRISPR/Cas9 editing sites, including 229,909,837 (~38.11%) NGG-PAM sites and 373,292,477 (~61.89%) NAG-PAM sites, were detectable in the pepper genome, and the systematic characterization of their composition and distribution was performed. Furthermore, 29,623,855 highly specific NGG-PAM sites were identified through whole-genome alignment analysis. There were 26,699,38 (~90.13%) highly specific NGG-PAM sites located in intergenic regions, which was 9.13 times of the number in genic regions, but the average density in genic regions was higher than that in intergenic regions. More importantly, 34,251 (~96.93%) out of 35,336 annotated genes exhibited at least one highly specific NGG-PAM site in their exons, and 90.50% of the annotated genes exhibited at least 4 highly specific NGG- PAM sites, indicating that the set of highly specific CRISPR/Cas9 editing sites identified in this study was widely applicable and conducive to the minimization of the off-target effects of CRISPR/Cas9 in pepper.

A predicted NEDD8 conjugating enzyme gene identified as a Capsicum candidate Rf gene using bulk segregant RNA sequencing

Cytoplasmic male sterility (CMS) is an important tool for producing F₁ hybrids, which can exhibit heterosis. The companion system, restorer-of-fertility (Rf), is poorly understood at the molecular level and would be valuable in producing restorer lines for hybrid seed production. The identity of the Rf gene in Capsicum (pepper) is currently unclear. In this study, using bulked segregant RNA sequencing (BSR-seq), a strong candidate Rf gene, Capana06g002866, which is annotated as a NEDD8 conjugating enzyme E2, was identified. Capana06g002866 has an ORF of 555 bp in length encoding 184 amino acids; it can be cloned from F₁ plants from the hybridization of the CMS line 8A and restorer line R1 but is not found in CMS line 8A. With qRT-PCR validation, Capana06g002866 was found to be upregulated in restorer accessions compared to sterile accessions. The relative expression in flower buds increased with the developmental stage in F₁ plants, while the expression was very low in all flower bud stages of the CMS lines. These results provide new insights into the Rf gene in pepper and will be useful for other crops utilizing the CMS system.