Capsicum annuum S (CaS) promotes reproductive transition and is required for flower formation in pepper (Capsicum annuum)

WUSCHEL-related homeobox transcription factor SlWOX13 regulates tomato fruit ripening

Fruit ripening is a complex, genetically programmed process involving the action of critical transcription factors (TFs). Despite the established importance of WUSCHEL-related homeobox (WOX) TFs in plant development, the involvement of WOX and its underlying mechanism in the regulation of fruit ripening remain unclear. Here, we demonstrate that SlWOX13 regulates fruit ripening in tomato (Solanum lycopersicum). Overexpression of SlWOX13 accelerates fruit ripening, whereas loss-of-function mutation in SlWOX13 delays this process. Moreover, ethylene synthesis and carotenoid accumulation are significantly inhibited in slwox13 mutant fruit but accelerated in SlWOX13 transgenic fruit. Integrated analyses of RNA-seq and chromatin immunoprecipitation (ChIP)-seq identified 422 direct targets of SlWOX13, of which 243 genes are negatively regulated and 179 are positively regulated by SlWOX13. Electrophoretic mobility shift assay, RT-qPCR, dual-luciferase reporter assay, and ChIP-qPCR analyses demonstrated that SlWOX13 directly activates the expression of several genes involved in ethylene synthesis and signaling and carotenoid biosynthesis. Furthermore, SlWOX13 modulates tomato fruit ripening through key ripening-related TFs, such as RIPENING INHIBITOR (RIN), NON-RIPENING (NOR), and NAM, ATAF1, 2, and CUC2 4 (NAC4). Consequently, these effects promote fruit ripening. Taken together, these results demonstrate that SlWOX13 positively regulates tomato fruit ripening via both ethylene synthesis and signaling and by transcriptional regulation of key ripening-related TFs.

Decoding Gene Expression Signatures Underlying Vegetative to Inflorescence Meristem Transition in the Common Bean

The tropical common bean (Phaseolus vulgaris L.) is an obligatory short-day plant that requires relaxation of the photoperiod to induce flowering. Similar to other crops, photoperiod-induced floral initiation depends on the differentiation and maintenance of meristems. In this study, the global changes in transcript expression profiles were analyzed in two meristematic tissues corresponding to the vegetative and inflorescence meristems of two genotypes with different sensitivities to photoperiods. A total of 3396 differentially expressed genes (DEGs) were identified, and 1271 and 1533 were found to be up-regulated and down-regulated, respectively, whereas 592 genes showed discordant expression patterns between both genotypes. Arabidopsis homologues of DEGs were identified, and most of them were not previously involved in Arabidopsis floral transition, suggesting an evolutionary divergence of the transcriptional regulatory networks of the flowering process of both species. However, some genes belonging to the photoperiod and flower development pathways with evolutionarily conserved transcriptional profiles have been found. In addition, the flower meristem identity genes APETALA1 and LEAFY, as well as CONSTANS-LIKE 5, were identified as markers to distinguish between the vegetative and reproductive stages. Our data also indicated that the down-regulation of the photoperiodic genes seems to be directly associated with promoting floral transition under inductive short-day lengths. These findings provide valuable insight into the molecular factors that underlie meristematic development and contribute to understanding the photoperiod adaptation in the common bean.

Morphological and gene expression characterization of maf-1, a floral chili pepper mutant caused by a nonsense mutation in CaLFY

Chili peppers are important as vegetables and ornamental crops, because of the variety of fruit shapes and colors. Understanding of flower and fruit development in Capsicum is limited compared with closely related Solanaceae crops such as tomato. This study reports a novel malformed fruit mutant named malformed fruit-1 (maf-1), which was isolated from an ethyl methanesulfonate-induced mutant population of chili pepper. maf-1 exhibited homeotic changes in the floral bud, which were characterized by conversion of petals and stamens into sepal-like and carpel-like organs, respectively. In addition, the indeterminate formation of carpel-like tissue was observed. Genetic analysis demonstrated that the causative gene in maf-1 is a nonsense mutation in CaLFY. This is the first characterization of an lfy mutant in Capsicum. Unlike tomatoes, the CaLFY mutation did not affect the architecture of sympodial unit or flowering time but mainly affected the formation of flower organs. Gene expression analysis suggested that a nonsense mutation in CaLFY led to decreased expression of multiple class B genes, resulting in homeotic changes in the flower and fruit. This maf-1 mutant may provide new insights at the molecular level in understanding flower organ formation and the genetic manipulation of fruit shape in chili peppers.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01304-w.

Identification of Genetic Factors Controlling the Formation of Multiple Flowers Per Node in Pepper (Capsicum spp.)

Flower production provides the foundation for crop yield and increased profits. Capsicum annuum is a pepper species with a sympodial shoot structure with solitary flowers. By contrast, C. chinense produces multiple flowers per node. C. annuum accounts for 80% of pepper production worldwide. The identification of C. chinense genes that control multiple flowers and their transfer into C. annuum may open the way to increasing fruit yield. In this study, we dissected the genetic factors were dissected controlling the multiple-flower-per-node trait in Capsicum. 85 recombinant inbred lines (RILs) between the contrasting C. annuum 'TF68' and C. chinense 'Habanero' accessions were phenotyped and genotyped. Quantitative Trait Loci (QTL) analysis identified four novel QTLs on chromosomes 1, 2, 7, and 11 that accounted for 65% of the total phenotypic variation. Genome-wide association study was also performed on a panel of 276 genotyped and phenotyped C. annuum accessions, which revealed 28 regions significantly associated with the multiple-flower trait, of which three overlapped the identified QTLs. Five candidate genes involved in the development of the shoot and flower meristems were identified and these genes could cause multiple flowers per node in pepper. These results contribute to our understanding of multiple flower formation in Capsicum and will be useful to develop high-yielding cultivars.

The APETALA2 homolog CaFFN regulates flowering time in pepper

Flowering time is an important agronomic trait that contributes to fitness in plants. However, the genetic basis of flowering time has not been extensively studied in pepper. To understand the genetics underlying flowering time, we constructed an F₂ population by crossing a spontaneous early flowering mutant and a late-flowering pepper line. Using bulked segregant RNA-seq, a major locus controlling flowering time in this population was mapped to the end of chromosome 2. An APETALA2 (AP2) homolog (CaFFN) cosegregated with flowering time in 297 individuals of the F₂ population. A comparison between the parents revealed a naturally occurring rare SNP (SNP2T > C) that resulted in the loss of a start codon in CaFFN in the early flowering mutant. Transgenic Nicotiana benthamiana plants with high CaFFN expression exhibited a delay in flowering time and floral patterning defects. On the other hand, pepper plants with CaFFN silencing flowered early. Therefore, the CaFFN gene acts as a flowering repressor in pepper. CaFFN may function as a transcriptional activator to activate the expression of CaAGL15 and miR156e and as a transcriptional repressor to repress the expression of CaAG, CaAP1, CaSEP3, CaSOC1, and miR172b based on a qRT-PCR assay. Direct activation of CaAGL15 by CaFFN was detected using yeast one-hybrid and dual-luciferase reporter assays, consistent with the hypothesis that CaFFN regulates flowering time. Moreover, the CaFFN gene association analysis revealed a significant association with flowering time in a natural pepper population, indicating that the CaFFN gene has a broad effect on flowering time in pepper. Finally, the phylogeny, evolutionary expansion and expression patterns of CaFFN/AP2 homologs were analyzed to provide valuable insight into CaFFN. This study increases our understanding of the involvement of CaFFN in controlling flowering time in pepper, thus making CaFFN a target gene for breeding early maturing pepper.

Conserved pleiotropy of an ancient plant homeobox gene uncovered by cis-regulatory dissection

Divergence of gene function is a hallmark of evolution, but assessing functional divergence over deep time is not trivial. The few alleles available for cross-species studies often fail to expose the entire functional spectrum of genes, potentially obscuring deeply conserved pleiotropic roles. Here, we explore the functional divergence of WUSCHEL HOMEOBOX9 (WOX9), suggested to have species-specific roles in embryo and inflorescence development. Using a cis-regulatory editing drive system, we generate a comprehensive allelic series in tomato, which revealed hidden pleiotropic roles for WOX9. Analysis of accessible chromatin and conserved cis-regulatory sequences identifies the regions responsible for this pleiotropic activity, the functions of which are conserved in groundcherry, a tomato relative. Mimicking these alleles in Arabidopsis, distantly related to tomato and groundcherry, reveals new inflorescence phenotypes, exposing a deeply conserved pleiotropy. We suggest that targeted cis-regulatory mutations can uncover conserved gene functions and reduce undesirable effects in crop improvement.

CaFT-LIKE is a flowering promoter in pepper and functions as florigen in tomato

We identified a pepper late-flowering mutant that is disrupted in the sequence of CaFT-LIKE, the ortholog of tomato SINGLE FLOWER TRUSS (SFT). Heterologous expression in tomato indicated that CaFT-LIKE has a conserved function as a flowering promoter and can rescue the wild-type phenotype of the tomato sft mutant. CaFT-LIKE confers a graft-transmissible signal for flowering initiation in tomato, implicating its function as a florigen. To test the relationship between CaFT-LIKE and FASCICULATE (FA), the ortholog of tomato SELF PRUNING (SP), we constructed the double mutant Caft-like fa. The phenotype of Caft-like fa resembled that of Caft-like, indicating epistasis of Caft-like over fa in controlling flowering time and sympodial shoot structure. To examine the association between the expression pattern of flowering genes and natural variation in flowering time, the expression levels of CaFT-LIKE and the flowering repressor CaAP2 were determined in a panel of early-flowering cultivars and late-flowering landraces and wild accessions. Strong positive and negative correlations between flowering time and expression levels of CaAP2 and CaFT-LIKE, respectively, were observed, indicating that high-expression alleles of CaFT-LIKE and low-expression alleles of CaAP2 were selected for early flowering during pepper domestication and breeding.

High-density genetic map construction and QTL mapping of first flower node in pepper (Capsicum annuum L.)

BACKGROUND

First flower node (FFN) is an important trait for evaluating fruit earliness in pepper (Capsicum annuum L.). The trait is controlled by quantitative trait loci (QTL); however, studies have been limited on QTL mapping and genes contributing to the trait.

RESULTS

In this study, we developed a high density genetic map using specific-locus amplified fragment sequencing (SLAF-seq), a high-throughput strategy for de novo single nucleotide polymorphism discovery, based on 146 recombinant inbred lines (RILs) derived from an intraspecific cross between PM702 and FS871. The map contained 9328 SLAF markers on 12 linkage groups (LGs), and spanned a total genetic distance of 2009.69 centimorgan (cM) with an average distance of 0.22 cM. The sequencing depth for the map was 72.39-fold in the male parent, 57.04-fold in the female parent, and 15.65-fold in offspring. Using the genetic map, two major QTLs, named Ffn2.1 and Ffn2.2, identified on LG02 were strongly associated with FFN, with a phenotypic variance explanation of 28.62 and 19.56%, respectively. On the basis of the current annotation of C. annuum cv. Criollo de Morelos (CM334), 59 candidate genes were found within the Ffn2.1 and Ffn2.2 region, but only 3 of 59 genes were differentially expressed according to the RNA-seq results. Eventually we identified one gene associated with the FFN based on the function through GO, KEGG, and Swiss-prot analysis.

CONCLUSIONS

Our research showed that the construction of high-density genetic map using SLAF-seq is a valuable tool for fine QTL mapping. The map we constructed is by far the most saturated complete genetic map of pepper, and using it we conducted fine QTL mapping for the important trait, FFN. QTLs and candidate genes obtained in this study lay a good foundation for the further research on FFN-related genes and other genetic applications in pepper.

Construction of a high density genetic map of an interspecific cross of Capsicum chinense and Capsicum annuum and QTL analysis of floral traits

The yield of pepper plants (Capsicum spp.) is their most important trait and is affected by the flower number and flowering time. Capsicum annuum produces a single flower per node and has an early flowering habit. By contrast, Capsicum chinense yields multiple flowers per node and has a late flowering character. However, the genetic mechanism underlying the control of these floral traits remains largely unknown. In this study, 150 F₂ populations from an interspecific cross between the inbred lines 740 (C. chinense) and CA1 (C. annuum) and their parents were used to construct a molecular genetic linkage map using the specific length amplified fragment sequencing (SLAF-seq) technique. This linkage map, spanning 1,586.78 cM in length, contained 9,038 markers on 12 chromosomes, with a mean marker distance of 0.18 cM. Phenotypic data on the flowering time and flower number per node were collected over multiple years, and QTL analysis identified 6 QTLs for the flowering time and flower number per node by composite interval mapping (CIM) and genome-wide composite interval mapping (GCIM) methods at least in two environments. The candidate genes within the major QTL were predicted. In the major flowering time QTL, the candidate gene Capana02g000700, which encodes the homeotic protein APETALA2, was identified. Quantitative reverse-transcription PCR (qRT-PCR) analysis indicated that its expression level in 740 was higher than that in CA1. Gene expression analysis indicated that the expression of Capana02g000700 was significantly upregulated in flowers, and many floral development-related genes were found to be coexpressed with Capana02g000700, supporting the function of this gene in association with flowering time in C. chinense and C. annuum species.

CaVIL1, a plant homeodomain gene that promotes flowering in pepper

CaVIL1 is a homolog of VIL1, a regulator of vernalization response in Arabidopsis and acts as a flowering promoter in pepper which does not respond to vernalization and photoperiod. As part of our goal to study the genetic and molecular basis of transition to flowering in pepper, we isolated the late-flowering mutant E-2698. Aside from late flowering, multiple pleiotropic alterations of the shoot structure, such as enlarged and distorted leaves, weak apical dominance, and reduced angle of the lateral branches were observed, indicating a broad role for the mutated gene in pepper development. Genetic mapping and sequence analyses revealed that the disrupted gene in E-2698 is the pepper homolog of VERNALIZATION INSENSITIVE 3-LIKE 1 (VIL1) that acts as a regulator of vernalization in Arabidopsis through chromatin modification. The pepper gene, CaVIL1, contains a plant homeodomain motif associated with chromatin modification and a VERNALIZATION INSENSITIVE 3-interacting domain that is truncated in E-2698 and in two other allelic mutants. Because pepper flowering does not respond to vernalization, we postulate that CaVIL1 regulates flowering time via chromatin modification of unknown targets. Expression analysis indicated that CaVIL1 activates the flowering promoter CaFLOWERING LOCUS T and represses the flowering repressor CaAPETALA2. Furthermore, CaVIL1 represses several genes from the FLOWERING LOCUS C (FLC)-LIKE clade that are clustered together in the pepper genome. This indicates their possible involvement in flowering regulation in this species. Our results show that CaVIL1 is a major regulator of flowering and interacts with other flowering promoters and repressors, as well as with FLC-LIKE genes whose function in flowering regulation is not yet known in pepper.

Deregulation of MADS-box transcription factor genes in a mutant defective in the WUSCHEL-LIKE HOMEOBOX gene EVERGREEN of Petunia hybrida

Angiosperm inflorescences develop in two fundamentally different ways. In monopodial plants, for example in Arabidopsis thaliana, the flowers are initiated as lateral appendages of a central indeterminate inflorescence meristem. In sympodial plants, flowers arise by terminal differentiation of the inflorescence meristem, while further inflorescence development proceeds from new sympodial meristems that are generated at the flank of the terminal flower. We have used the sympodial model species Petunia hybrida to investigate inflorescence development. Here, we describe a mutant, bonsai (bns), which is defective in flower formation, inflorescence branching, and control of meristem size. Detailed microscopic analysis revealed that bns meristems retain vegetative charateristics including spiral phyllotaxis. Consistent with a block in flower formation, bns mutants exhibit a deregulated expression of various MADS-box genes. Molecular analysis revealed that the bns mutant carries a transposon insertion in the previously described EVERGREEN (EVG) gene, which belongs to the WUSCHEL-LIKE HOMEOBOX (WOX) transcription factor gene family. EVG falls in the WOX9 subfamily, which has diverse developmental functions in angiosperms. The comparison of WOX9 orthologues in five model species for flowering shows that these genes play functionally divergent roles in monopodial and sympodial plants, indicating that the WOX9 regulatory node may have played an important role in the evolution of shoot architecture.

Candidate genes for first flower node identified in pepper using combined SLAF-seq and BSA

First flower node (FFN) is an important trait for evaluating fruit earliness in pepper (Capsicum annuum L.), but the genetic mechanisms that control FFN are still poorly understood. In the present study, we developed 249 F₂ plants derived from an intraspecific cross between the inbred pepper lines Z4 and Z5. Thirty plants with the highest FFN and 30 plants with the lowest FFN were chosen and their DNAs were pooled according to phenotype to construct two bulked DNA pools. Specific-locus amplified fragment sequencing (SLAF-seq) was combined with bulked segregant analysis (BSA) to identify candidate regions related to FFN. According to our genetic analysis, the FFN trait is quantitatively inherited. A total of 106,848 high-quality single nucleotide polymorphism (SNP) markers were obtained, and 393 high-quality SNP markers associated with FFN were detected. Ten candidate regions within an interval of 3.98 Mb on chromosome 12 harboring 23 candidate genes were identified as closely correlated with FFN. Five genes (CA12g15130, CA12g15160, CA12g15370, CA12g15360, and CA12g15390) are predicted based on their annotations to be associated with expression of the FFN trait. The present study demonstrates an efficient genetic mapping strategy and lays a good foundation for molecular marker-assisted breeding using SNP markers linked to FFN and for cloning and functional analysis of the key genes controlling FFN.

The evolution of inflorescence diversity in the nightshades and heterochrony during meristem maturation

One of the most remarkable manifestations of plant evolution is the diversity for floral branching systems. These “inflorescences” arise from stem cell populations in shoot meristems that mature gradually to reproductive states in response to environmental and endogenous signals. The morphology of the shoot meristem maturation process is conserved across distantly related plants, raising the question of how diverse inflorescence architectures arise from seemingly common maturation programs. In tomato and related nightshades (Solanaceae), inflorescences range from solitary flowers to highly branched structures bearing hundreds of flowers. Since reproductive barriers between even closely related Solanaceae have precluded a genetic dissection, we captured and compared meristem maturation transcriptomes from five domesticated and wild species reflecting the evolutionary continuum of inflorescence complexity. We find these divergent species share hundreds of dynamically expressed genes, enriched for transcription factors. Meristem stages are defined by distinct molecular states and point to modified maturation schedules underlying architectural variation. These modified schedules are marked by a peak of transcriptome expression divergence during the reproductive transition, driven by heterochronic shifts of dynamic genes, including transcriptional regulators with known roles in flowering. Thus, evolutionary diversity in Solanaceae inflorescence complexity is determined by subtle modifications of transcriptional programs during a critical transitional window of meristem maturation, which we propose underlies similar cases of plant architectural variation. More broadly, our findings parallel the recently described transcriptome “inverse hourglass” model for animal embryogenesis, suggesting both plant and animal morphological variation is guided by a mid-development period of transcriptome divergence.

Development of a SNP array and its application to genetic mapping and diversity assessment in pepper (Capsicum spp.)

The development and application of single nucleotide polymorphisms (SNPs) is in its infancy for pepper. Here, a set of 15,000 SNPs were chosen from the resequencing data to develop an array for pepper with 12,720 loci being ultimately synthesized. Of these, 8,199 (~64.46%) SNPs were found to be scorable and covered ~81.18% of the whole genome. With this array, a high-density interspecific genetic map with 5,569 SNPs was constructed using 297 F2 individuals, and genetic diversity of a panel of 399 pepper elite/landrace lines was successfully characterized. Based on the genetic map, one major QTL, named Up12.1, was detected for the fruit orientation trait. A total of 65 protein-coding genes were predicted within this QTL region based on the current annotation of the Zunla-1 genome. In summary, the thousands of well-validated SNP markers, high-density genetic map and genetic diversity information will be useful for molecular genetics and innovative breeding in pepper. Furthermore, the mapping results lay foundation for isolating the genes underlying variation in fruit orientation of Capsicum.

CaAP2 transcription factor is a candidate gene for a flowering repressor and a candidate for controlling natural variation of flowering time in Capsicum annuum

The APETALA2 transcription factor homolog CaAP2 is a candidate gene for a flowering repressor in pepper, as revealed by induced-mutation phenotype, and a candidate underlying a major QTL controlling natural variation in flowering time. To decipher the genetic control of transition to flowering in pepper (Capsicum spp.) and determine the extent of gene function conservation compared to model species, we isolated and characterized several ethyl methanesulfonate (EMS)-induced mutants that vary in their flowering time compared to the wild type. In the present study, we report on the isolation of an early-flowering mutant that flowers after four leaves on the primary stem compared to nine leaves in the wild-type 'Maor'. By genetic mapping and sequencing of putative candidate genes linked to the mutant phenotype, we identified a member of the APETALA2 (AP2) transcription factor family, CaAP2, which was disrupted in the early-flowering mutant. CaAP2 is a likely ortholog of AP2 that functions as a repressor of flowering in Arabidopsis. To test whether CaAP2 has an effect on controlling natural variation in the transition to flowering in pepper, we performed QTL mapping for flowering time in a cross between early and late-flowering C. annuum accessions. We identified a major QTL in a region of chromosome 2 in which CaAP2 was the most significant marker, explaining 52 % of the phenotypic variation of the trait. Sequence comparison of the CaAP2 open reading frames in the two parents used for QTL mapping did not reveal significant variation. In contrast, significant differences in expression level of CaAP2 were detected between near-isogenic lines that differ for the flowering time QTL, supporting the putative function of CaAP2 as a major repressor of flowering in pepper.

Construction of an interspecific genetic map based on InDel and SSR for mapping the QTLs affecting the initiation of flower primordia in pepper (Capsicum spp.)

Re-sequencing permits the mining of genome-wide variations on a large scale and provides excellent resources for the research community. To accelerate the development and application of molecular markers and identify the QTLs affecting the flowering time-related trait in pepper, a total of 1,038 pairs of InDel and 674 SSR primers from different sources were used for genetic mapping using the F2 population (n = 154) derived from a cross between BA3 (C. annuum) and YNXML (C. frutescens). Of these, a total of 224 simple PCR-based markers, including 129 InDels and 95 SSRs, were validated and integrated into a map, which was designated as the BY map. The BY map consisted of 13 linkage groups (LGs) and spanned a total genetic distance of 1,249.77 cM with an average marker distance of 5.60 cM. Comparative analysis of the genetic and physical map based on the anchored markers showed that the BY map covered nearly the whole pepper genome. Based on the BY map, one major and five minor QTLs affecting the number of leaves on the primary axis (Nle) were detected on chromosomes P2, P7, P10 and P11 in 2012. The major QTL on P2 was confirmed based on another subset of the same F2 population (n = 147) in 2014 with selective genotyping of markers from the BY map. With the accomplishment of pepper whole genome sequencing and annotations (release 2.0), 153 candidate genes were predicted to embed in the Nle2.2 region, of which 12 important flowering related genes were obtained. The InDel/SSR-based interspecific genetic map, QTLs and candidate genes obtained by the present study will be useful for the downstream isolation of flowering time-related gene and other genetic applications for pepper.