Next Article in Journal
Steroid-Induced Ocular Hypertension in Mice Is Differentially Reduced by Selective EP2, EP3, EP4, and IP Prostanoid Receptor Agonists
Next Article in Special Issue
ERECTA Modulates Seed Germination and Fruit Development via Auxin Signaling in Tomato
Previous Article in Journal
Antibacterial Properties of Rose Bengal Conjugated to Hyaluronic Acid
Previous Article in Special Issue
Transcriptome Analysis and Metabolic Profiling Reveal the Key Regulatory Pathways in Drought Stress Responses and Recovery in Tomatoes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 6
10.3390/ijms25063331
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

Candidate Gene Identification and Transcriptome Analysis of Tomato male sterile-30 and Functional Marker Development for ms-30 and Its Alleles, ms-33, 7B-1, and stamenless-2

by
Kai Wei
1,2,†,
Xin Li
1,†,
Xue Cao
1,‡,
Shanshan Li
1,
Li Zhang
1,2,
Feifei Lu
1,
Chang Liu
1,
Yanmei Guo
1,
Lei Liu
1,
Can Zhu
1,
Yongchen Du
1,
Junming Li
1,
Wencai Yang
2,
Zejun Huang
1,* and
Xiaoxuan Wang
1,*
1
State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
2
Department of Vegetable Science, College of Horticulture, China Agricultural University, Beijing 100193, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
‡
Current address: State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing 100700, China.
Int. J. Mol. Sci. 2024, 25(6), 3331; https://doi.org/10.3390/ijms25063331
Submission received: 15 February 2024 / Revised: 6 March 2024 / Accepted: 7 March 2024 / Published: 15 March 2024
(This article belongs to the Special Issue Advances in Tomato Breeding and Molecular Research)
Browse Figures
Review Reports Versions Notes

Abstract

:
Male sterility is a valuable trait for hybrid seed production in tomato (Solanum lycopersicum). The mutants male sterile-30 (ms-30) and ms-33 of tomato exhibit twisted stamens, exposed stigmas, and complete male sterility, thus holding potential for application in hybrid seed production. In this study, the ms-30 and ms-33 loci were fine-mapped to 53.3 kb and 111.2 kb intervals, respectively. Tomato PISTILLATA (TPI, syn. SlGLO2), a B-class MADS-box transcription factor gene, was identified as the most likely candidate gene for both loci. TPI is also the candidate gene of tomato male sterile mutant 7B-1 and sl-2. Allelism tests revealed that ms-30, ms-33, 7B-1, and sl-2 were allelic. Sequencing analysis showed sequence alterations in the TPI gene in all these mutants, with ms-30 exhibiting a transversion (G to T) that resulted in a missense mutation (S to I); ms-33 showing a transition (A to T) that led to alternative splicing, resulting in a loss of 46 amino acids in protein; and 7B-1 and sl-2 mutants showing the insertion of an approximately 4.8 kb retrotransposon. On the basis of these sequence alterations, a Kompetitive Allele Specific PCR marker, a sequencing marker, and an Insertion/Deletion marker were developed. Phenotypic analysis of the TPI gene-edited mutants and allelism tests indicated that the gene TPI is responsible for ms-30 and its alleles. Transcriptome analysis of ms-30 and quantitative RT-PCR revealed some differentially expressed genes associated with stamen and carpel development. These findings will aid in the marker-assisted selection for ms-30 and its alleles in tomato breeding and support the functional analysis of the TPI gene.

1. Introduction

Tomato (Solanum lycopersicum) is among the most crucial vegetables globally. According to the Food and Agriculture Organization of the United Nations statistical database 2021, approximately 189 million tons of fresh tomatoes are produced worldwide (https://www.fao.org/faostat/zh/#data/QCL). The widespread adoption of tomato F1 hybrids can be attributed to their consistently higher yields, improved quality, and enhanced resistance to diseases compared with open-pollinated varieties [1,2]. However, the traditional method of producing tomato hybrid seeds involves labor-intensive processes, such as manual emasculation and hand pollination, leading to increased costs and the risk of self-pollination impurities [3,4]. To address these challenges, researchers have explored the suitability of male sterile plants as female progenitors in the production of tomato hybrid seeds [4,5]. Different genes controlling male sterility have been examined, and their applications in plant breeding have been evaluated.
Since the first report of a male sterile mutant in tomatoes in 1915, over 55 male sterile alleles have been identified and categorized into three types: sporogenous, functional, and structural sterility [5,6]. Most tomato male sterile mutants are of the sporogenous type [7], with only a few having been cloned, such as SlMS10 and SlMS32.Both of these encode a basic helix–loop–helix (bHLH) transcription factor [8,9,10] and are homologous to DYT1 and bHLH10/89/90 in Arabidopsis, respectively [8,10]. DYT1 and bHLH10/89/90 are part of the DYT1-TDF1-AMS-MS188-MS1 pathway that regulates tapetum development and pollen formation in Arabidopsis [11]. Positional sterile (ps) is the first reported mutant of functional male sterility in tomato [12]. It produces normal pollen; however, its pollen grains remain enclosed within the anther locule due to non-splitting anthers [13]. The second reported mutant of this type is ps-2. PS-2 encodes a polygalacturonase gene (PG) [14,15], and it has been used for tomato hybrid seed production in Eastern Europe [16]. Moreover, functional male sterility can result from exserted stigma. In the exserted stigma (ex) mutant, the anthers are normal and capable of producing and releasing viable pollens, but the excessively long style prevents effective pollination [17]. Other mutants of this type include Style2.1, Stigma Exertion 3.1 (SE3.1), and SlLst (Long styles) [18,19,20]. Stamenless (sl), the first structural male sterile mutant reported in tomato, exhibited the complete transformation of stamens into carpels and the partial transformation of petals into sepals [21]. Tomato APETALA3 (TAP3) may be responsible for the sl trait [22,23]. In the sl-2 mutant, petals appear almost typical, but the stamens are twisted and distorted, bearing naked ovules [24]. The male sterile mutant 7B-1 displays a high degree of phenotypic similarity with sl-2 and an allelism test confirmed that 7B-1 and sl-2 are allelic [25]. Genetic mapping and molecular analysis indicated that the B-class MADS-box gene Tomato PISTILLATA (TPI, syn. SlGLO2) is a candidate gene of 7B-1 and sl-2; however, the specific sequence variations of TPI in 7B-1 and sl-2 mutants remain elusive [25]. Tomato ms-15 and its alleles, ms-1526 and ms-1547, producing flowers with homeotic conversion of stamens into carpels, have been fine-mapped to a 44.6 kb interval that contains the TM6 gene [26]. A recent study revealed that TM6 maintains the balance of alternative splicing of MADS-box genes involved in flower development [27].
Each type of male sterility in tomato presents both advantages and disadvantages for use in hybrid seed production. The lines ms10 and ms32 are considered beneficial for hybrid seed production [28], because they exhibit exserted stigma (accessible for pollination without emasculation), complete sterility, and stable expression of sterility regardless of environmental conditions. However, the application of these two sterile lines has some major drawbacks, such as the variation in the percentage of flowers with accessible stigmas depending on the environment and genotype [4]. The exserted stigma (ex) sterility type was initially considered advantageous for hybrid seed production due to the easy maintenance of sterile lines, absence of a need for stamen emasculation, and high hybrid seed yield. However, it presents the challenge of managing two polygenic traits (anther and style length) and the variability of phenotypic expression under different environmental conditions. ps-2 has disadvantages due to the risk of undesirable self-pollination and the necessity for stamen emasculation [4,16]. Although 7B-1 can be used for hybrid seed production [29], the absence of a functional codominant marker for 7B-1 reduces the efficiency of backcross transfer. These deficiencies in male sterile lines for tomato hybrid seed production necessitate further exploration of male sterility genes and the development of codominant molecular markers for marker-assisted selection (MAS).
Tomato male sterile mutants ms-30 and ms-33 feature shorter, twisted stamens and exerted stigmas, facilitating direct pollination without the need for emasculation [30,31]. Among these mutants, ms-33 appears highly promising as a female parent in the production of tomato hybrid seeds [32]. Linkage analysis revealed that the ms-33 locus is located on chromosome 6 [33]. The fertility of ms-33 can be partially or completely restored with the application of exogenous GA3, resulting in seeds that produce 100% male sterile plants [7,34]. Given these characteristics, ms-30 and ms-33 can be valuable for hybrid seed production. However, the specific genes responsible for ms-30 and ms-33 sterility remain unknown, and there are currently no codominant functional markers available for selecting these mutants.
The objective of this study was to characterize the genes MS-30 and MS-33, along with their sequence variations in wild-type plants and the male sterile mutants ms-30 and ms-33. Subsequently, the study aimed to develop codominant functional markers for MAS utilizing these sequence variations, and to perform an initial investigation into the role of the MS-30 gene in stamen development through RNA-seq analysis.

2. Results

2.1. Phenotypes of the ms-30 and ms-33 Mutants

The flowers of WT-30 (the wild-type corresponding to ms-30) and WT-33 exhibited normal architecture, in contrast to the stamens of ms-30 and ms-33, which displayed distinctive features, such as twisted stamens, lateral separation, reduced size, and pale color, resulting in an exerted stigma (Figure 1 and Figure S1A). Additionally, these distorted stamens were observed to transform into carpel-like structures, including carpelloid structures (CS), naked external ovules on the adaxial surface (EO), and complete transformation into carpels (TC). Notably, the occurrence of EO was more frequent than that of CS and TC in both ms-30 and ms-33 (Table S1). Due to these alterations, neither mutant was capable of producing pollens (Figure S1B), and both ms-30 and ms-33 were unable to set fruit when subjected to self-pollination using an electric vibrator. However, manual pollination with pollens from wild-type plants resulted in a successful fruit set (Figure S1C).

2.2. Genetic Analysis of ms-30 and ms-33 Loci

The F1 plants resulting from crosses between ms-30 and LA1589, as well as between ms-33 and LA1589, displayed complete fertility. Additionally, a preliminary analysis involved 200 plants in each of the F2 populations derived from these crosses. The ratio of fertile to sterile plants for both ms-30 and ms-33 conformed to the expected 3:1 segregation ratio, as confirmed by the chi-square test (Table 1). The observed segregation pattern and the statistical analysis suggested that the male sterility in both ms-30 and ms-33 is controlled by a single recessive gene.

2.3. Fine-Mapping of ms-30 and ms-33 Loci

For preliminary mapping of the ms-30 and ms-33 loci, markers located on chromosome 6 were utilized. This approach was based on prior reports indicating that ms-33 is located on chromosome 6 of tomato. Furthermore, the ms-30 and ms-33 mutants exhibited similar phenotypes, suggesting a potential linkage (Figure 1 and Figure S1; Table S1) [30,31,33,35,36]. In this preliminary mapping phase, 94 male sterile plants from each F2 population were genotyped, and their recombination rates were calculated. As anticipated, the findings confirmed that both ms-30 and ms-33 are located within the region bounded by markers HP149 and HP3661 on chromosome 6 (Figure 2A and Figure S2A). To further narrow down the locations of these loci, 76 and 39 recombinants were selected from the entire F2 populations of ms-30 and ms-33, respectively, and genotyped with additional markers. This detailed analysis allowed the ms-30 locus to be narrowed down to a 53.3 kb region, delimited by the markers CP207 and CP203 (Figure 2B). Similarly, the ms-33 locus was confined to an interval of 111.2 kb, defined by markers CP207 and HP3921 (Figure S2B). Notably, both loci shared an identical 53.3 kb region.

2.4. Candidate Gene Analysis and Allelism Testing among ms-30, ms-33, and sl-2

According to the ITAG4.0 genome annotation [37], Solyc06g059970, which encodes Tomato PISTILLATA (TPI, synonymously referred to as SlGLO2) [38,39], was identified within the fine-mapped regions of both the ms-30 and ms-33 loci (Figure 2C and Figure S2C; Table 2 and Table S2). TPI plays a pivotal role in stamen development and has been identified as the candidate gene for 7B-1 and sl-2 [25], which exhibit phenotypes similar to those of ms-30 and ms-33 (Figure 1 and Figure S1). To further ascertain the genetic relationships among ms-30, ms-33, and sl-2, allelism tests were performed. The outcomes of these tests confirmed that ms-30, ms-33, and sl-2 are allelic, indicating that they arise from mutations in the same gene (Table S3). Thus, TPI was identified as the most probable candidate gene responsible for the phenotypic traits observed in ms-30 and ms-33.

2.5. Sequence Analysis of the TPI Gene in ms-30, ms-33, 7B-1, and sl-2

To identify sequence variations within the TPI gene between wild-type plants and ms mutants, the genomic sequence of the TPI gene was sequenced in WT-30, ms-30, WT-33, and ms-33. Sequencing data analysis revealed specific mutations in each mutant. In ms-30, a single-nucleotide polymorphism (SNP) was discovered in the first exon of the TPI gene (Figure 3A), leading to a missense mutation. This mutation altered a serine residue to an isoleucine within the MADS domain of the protein, impacting its function (Figure S3). In ms-33, a different SNP was identified in the first intron of the TPI gene (Figure 3B). Amplification of cDNAs from ms-33 produced several bands, with one being smaller than the corresponding band from the wild-type (Figure S4A). Further analysis indicated that this SNP might introduce an alternative 5’ splice site in ms-33’s gene transcript. This alternative splicing resulted in a 138 bp deletion in cDNA, leading to the protein lacking a 46 amino-acid fragment (Figure S3).
Pucci et al. (2017) [25] reported that the primer pair LeGLO2-6F/R did not amplify the last exon and part of the 3’ UTR region of the TPI gene in 7B-1 and sl-2. To investigate sequence variation in this region, PCR amplification was conducted using the LeGLO2-6F/R primer pair with an extended elongation time. This method enabled the production of approximately 5000 bp long amplicons from both 7B-1 and sl-2 (Figures S4B and S5). Sequencing of these amplicons showed no sequence variations. However, a fragment of 4868 bp in the amplicon was identified as a Copia-like retrotransposon, featuring a 398 bp long terminal repeat (LTR) sequence and a 4 bp target site duplication (TSD) sequence (Figure 3C and Figure S5). Based on the SNPs in the TPI gene of ms-30 and ms-33 and the LTR retrotransposon insertion in the TPI gene of 7B-1 and sl-2, three codominant markers were developed, namely a Kompetitive Allele Specific PCR (KASP) marker, MS30KASP; a sequence marker, MS33SEQ; and an Insertion/Deletion (InDel) marker, LeGLO2-6/SL-2TE (Figure 3D–F).

2.6. Targeted Knockout of the TPI Gene by CRISPR/Cas9

To confirm that the TPI gene confers the phenotype observed in ms-30 and its alleles, the TPI gene was knocked out using the CRISPR/Cas9 system, and allelism tests were conducted between the ms-30 and TPI gene-edited mutants. Two vectors were designed to disrupt the TPI gene, with target sites on the first and last exons for vector-1 and vector-2, respectively. This resulted in two null mutants of the TPI gene. The tpicr-e1 mutant exhibited a 4 bp deletion, whereas tpicr-e7 had a 66 bp deletion (Figure 4A). The stamens of both tpicr-e1 and tpicr-e7 mutants showed curling, wrinkling, greenish coloration, and were incapable of producing viable pollen (Figure 4B–J), similar to phenotypic traits observed in the ms-30 and ms-33 mutants. Allelism tests between tpicr-e1, tpicr-e7, and ms-30 indicated that they were allelic (Table S3). These findings support the conclusion that TPI is responsible for the observed phenotypes in ms-30 and its alleles, including ms-33, 7B-1, and sl-2.

2.7. Comparative RNA-Seq Analysis of WT-30 and ms-30 Floral Organs from the Meiotic to Tetrad Stage

To further explore the role of TPI in floral organ development, RNA-seq analysis was conducted on the four floral organ whorls of young flower buds from WT-30 and ms-30 at the meiotic and tetrad stages (3–5 mm in length) [40]. RNA-Seq generated approximately 568 million paired-end reads, with each sample yielding between 20 million and 30 million reads. On average, 98% of these reads were successfully mapped to the ITAG4.0 tomato reference genome (Table S4). Principal component analysis (PCA) results showed that biological replicates clustered together, indicating good consistency within the data (Figure 5A). Notably, the stamens of ms-30 (ST_ms-30) were closely related to the carpels of WT-30 (CA_WT-30) and the carpels of ms-30 (CA_ms-30) in the PCA plot. This observation aligns with the homeotic transformation of stamens into carpels observed in the ms-30 mutant, providing molecular evidence for this phenotypic change.
Among the differentially expressed genes (DEGs) identified between WT-30 and ms-30, 96 were in sepals, 738 in petals, 5391 in stamens, and 131 in carpels. Specifically, within the stamens, 2238 genes were upregulated, whereas 3153 genes were downregulated (Figure 5B, Tables S5–S9). Regarding the ABCDE model of floral organ development, the expression of the D-class gene TAGL11 was not detected in any of the four floral organ whorls. The expression levels of other ABCDE model genes, with the exceptions of SlGLO1 and TAG1, differed significantly in the stamens between WT-30 and ms-30 (Figure S6).
Gene ontology (GO) analysis revealed that upregulated DEGs in stamens were enriched in specific biological process categories, such as cell division, regulation of transcription, DNA-templated synthesis, response to auxin, ethylene-activated signaling pathways, and response to light stimulus (Figure 5C). Conversely, DEGs downregulated in stamens were associated with categories such as response to water deprivation, response to abscisic acid, pollen sperm cell differentiation, and pollen tube growth (Figure 5D). These findings suggest that the expression of the identified DEGs is associated with key processes such as pollen development, hormonal responses, and responses to environmental stimuli.
Given the homeotic conversion of stamens into carpel-like organs observed in ms-30, our analysis focused on the DEGs related to stamen and carpel development. DEGs with homologs involved in stamen or carpel development in Arabidopsis were identified in the stamens of ms-30. Some of these genes have been previously reported to play roles in stamen or carpel development in tomato (Figure 6A,B, Tables S10 and S11). For example, genes associated with tomato tapetum development, such as SlMS10 (Solyc02g079810), Solyc03g113530 (an AtTDF1 homolog), SlAMS (Solyc08g062780), SlMS32 (Solyc01g081100), and SlPHD_MS1 (Solyc04g008420) were downregulated in the stamens of the ms-30 mutant (Figure 6C) [8,9,10,41,42]. Conversely, genes related to tomato carpel or ovule development, including SlCRCa (Solyc01g010240), SlCRCb (Solyc05g012050), SlINO (Solyc05g005240), LYRATE (Solyc05g009380), and SlMBP3 (Solyc06g064840), were upregulated in the stamens of the ms-30 mutant (Figure 6D) [43,44,45,46]. These expression differences were validated through qRT-PCR assays (Figure 6C,D). Moreover, qRT-PCR analysis showed that these genes were also differentially expressed in the ms-33, tpicr-e1, and tpicr-e7 mutants (Figures S7 and S8), further supporting their involvement in the observed phenotypic alterations.

3. Discussion

Using male sterile lines as the female parent in hybrid tomato seed production is considered a highly efficient and cost-effective method [28]. Several male sterile lines, including ms-10, ms-32, sl, 7B-1, ex, ps, and ps-2, [4,29] have been successfully applied in tomato hybrid seed production [4,29]. However, these lines have drawbacks, such as the risk of self-pollination, the need for anther emasculation, environmental dependency of phenotypic expression, reduced hybrid seed yield, or the lack of codominant markers for MAS [4]. Therefore, identifying additional genes responsible for male sterility in tomatoes, particularly those resulting in exposed stigmas and stable, complete male sterility, is crucial. Developing codominant markers for MAS of these male sterile loci is equally important. The ms-30 and ms-33 mutants are promising for hybrid tomato seed production due to their complete male sterility, exerted stigmas (Figure 1), and the ability to be manually pollinated directly without the need for hand emasculation. Both ms-30 and ms-33 exhibit single-gene recessive inheritance, with only homozygous plants displaying the male sterile phenotype, making MAS more efficient than conventional phenotypic selection for these mutants in breeding programs [47]. The ms-30 and ms-33 loci have been fine-mapped, identifying the TPI gene as the candidate gene for both loci. These mutants were found to be allelic to 7B-1 and sl-2 [25], which also involve TPI. Phenotypic analysis and allelism tests of TPI gene-edited mutants confirmed that TPI underlies the male sterile phenotype in ms-30 and its alleles, including ms-33, 7B-1, and sl-2. Sequence analysis identified variations in the TPI gene across these mutants. Based on these variations, three codominant markers—KASP marker MS30KASP, sequencing marker MS33SEQ, and InDel marker LeG-LO2-6/SL-2TE—were developed, facilitating rapid incorporation of these loci into tomato breeding. Additionally, CRISPR/Cas9-mediated targeted mutagenesis represents another approach to rapidly introduce male sterility into elite breeding lines [9,48,49,50,51,52]. The phenotypes of two TPI gene-edited mutants closely matched those of ms-30 and its alleles, suggesting that gene editing of TPI could efficiently induce male sterility in elite lines, thereby reducing the costs associated with hybrid seed production.
Carpelloid stamens are occasionally observed in angiosperms, with mutations in B-class genes identified as the primary genetic cause of this phenomenon [53]. The number of B-class genes varies across different plant species; for instance, Arabidopsis has two B-class genes (AP3 and PI), whereas tomato contains four genes (TAP3, TM6, TPI, and SlGLO1) [38,54]. This variation suggests potential functional conservation and diversification of B-class genes among plant species [54]. The role of B-class genes in stamen specification and development has been extensively studied since the 1990s [55,56,57]; however, further research is warranted. In this study, the ms-30 and ms-33 mutants displayed carpelloid stamens (Figure 1) and were found to be allelic to 7B-1 and sl-2 (Table S3). The B-class gene TPI was identified as the underlying gene for ms-30 and its alleles (Figure 2 and Figure S2). Transcriptome analysis of ms-30 was conducted to examine the function of TPI in stamen development. PCA of the transcriptome data showed that the clusters for stamens and carpels were closely related (Figure 5A), aligning with the observed phenotype of carpelloid stamens in ms-30 (Figure 1). The number of DEGs was higher in stamens than in the other three floral organs (Figure 5B). Some DEGs in stamens, specifically expressed in wild-type stamens, were downregulated in the ms-30 mutant (Figure 6A). Examples of such genes are SlMS10 (Solyc02g079810), Solyc03g113530 (AtTDF1-like1), SlAMS (Solyc08g062780), SlMS32 (Solyc01g081100), and SlPHD_MS1 (Solyc04g008420) (Figure 6C), which are key components of the DYT1-TDF1-AMS-MS188-MS1 regulatory network [8,9,10,41,42]. Conversely, DEGs that were specifically expressed in wild-type carpels and upregulated in ms-30 stamens included SlCRCa (Solyc01g010240), SlCRCb (Solyc05g012050), SlINO (Solyc05g005240), LYRATE (Solyc05g009380), and SlMBP3 (Solyc06g064840) (Figure 6D), suggesting their involvement in carpel or ovule development in tomato [43,44,45,46]. These findings suggest that TPI positively regulates the genes associated with stamen and pollen development while negatively influencing the genes involved in carpel and ovule development within stamens. Consequently, loss of function of the TPI gene leads to the low or absent expression of genes associated with stamen development and the ectopic expression of genes related to carpel and ovule development in stamens, potentially resulting in the formation of carpelloid stamens in mutants. However, further investigation is needed to elucidate the specific molecular regulatory mechanisms.

4. Materials and Methods

4.1. Plant Materials

Seeds for the wild-type (WT) plants and ms mutants (ms-30, accession number 2-455; ms-33, 2-511; sl-2, LA1801; Ailsa Craig [AC]; LA1589 and Heinz 1706) were obtained from the Tomato Genetics Resource Center (Davis, CA, USA). The male sterile plants derived from the accessions 2-455, 2-511, and LA1801 were designated as ms-30, ms-33, and sl-2, respectively. Similarly, the homozygous male fertile plants from 2-455 and 2-511 were identified as WT-30 and WT-33, respectively. AC was used for gene editing and as WT controls for the tpicr-e1 and tpicr-e7 mutants. Two F2 populations were created by crossing both ms-30 and ms-33 with the wild tomato species Solanum pimpinellifolium accession LA1589. These populations were cultivated in an open field in Shunyi District, Beijing, China, during the spring and summer of 2016, comprising 949 and 861 plants, respectively.
To confirm the allelism among ms-30, ms-33, and sl-2, the flowers of sl-2 were pollinated using pollens from the heterozygous ms-30 or ms-33 plants. Additionally, to determine whether the TPI gene is responsible for conferring the ms-30 phenotype, flowers of ms-30 were pollinated with pollens from the heterozygous TPI gene-edited mutants. Conversely, male sterile flowers of TPI gene-edited mutants were pollinated with pollens from the heterozygous ms-30 plants. The progeny from these crosses were grown in a greenhouse located in Haidian District, Beijing, China, during the springs and summers of 2019, 2020, and 2023.

4.2. Phenotypic Analysis

For phenotypic and expression analyses, WT-30, ms-30, WT-33, ms-33, and AC plants, as well as gene-edited mutants, were grown in a greenhouse in Haidian District, Beijing, China, across the springs and summers of 2019, 2020, 2021, and 2023. The flower morphology of both WT and ms mutants was examined at the anthesis stage using the previously described methods [25]. The entire flower was photographed by a camera (Canon EOS 70D, Canon Inc., Tokyo, Japan). The anther cones and dissected stamens were observed by a stereomicroscope (Carl Zeiss MicroImaging GmbH, Gottingen, Germany). Furthermore, referring to the procedures outlined by Cao et al., whole plants were photographed and pollen viability was evaluated [26]. Pollen was collected from the WT and ms plants, stained with 1% aceto-carmine, and examined using a microscope (BX51, Olympus Corporation, Tokyo, Japan). Based on crossing of the male sterile mutants ms-30 and ms-33 with Heinz 1706, we collected the fruit set data.

4.3. Molecular Marker Development and Genotyping

InDels and SNPs were identified by comparing the sequence of chromosome 6 between the tomato lines Heinz 1706 and LA1589. Whole-genome sequences for these lines were obtained from the Sol Genomics Network (SGN, https://solgenomics.net/) [58]. To evaluate polymorphisms between ms-30/ms-33 and LA1589, PCR primers were designed based on the flanking regions of these InDels and SNPs.
KASP genotyping reactions were performed using a LightCycler480I (Roche Diagnostics GmbH, Mannheim, Germany) with KASP-TF V4.0 2× Master Mix with Low Rox (Cat. No. KBS-1050-122, LGC Genomic Limited, Middlesex, United Kingdom) following the provided manual protocol (https://www.biosearchtech.com/support/education/kasp-genotyping-reagents/running-kasp-genotyping-reactions). The TPI gene fragment for ms-33 was amplified using the marker MS33SEQ and sequenced at the Beijing Genomics Institute (Beijing, China). To identify the insertion of the LTR retrotransposon in the 7B-1 and sl-2 mutants, PCR was conducted separately with primers LeGLO2-6F/R and SL-2TE-F/R. The resulting PCR products were then combined in a single tube and analyzed through electrophoresis.

4.4. Genetic Analysis and Fine-Mapping

Genetic analysis and preliminary mapping were performed using 200 plants and 94 male sterile plants selected from each F2 population. For fine-mapping of the ms-30 and ms-33 loci, all plants in the F2 populations were used. The goodness-of-fit to a 3:1 segregation ratio at the ms-30 and ms-33 loci was tested using the chi-square test with 200 individuals from each F2 population. For preliminary mapping of the ms-30 and ms-33 loci, the first 94 sterile plants in each F2 population were genotyped with 2 and 6 InDel markers, respectively. For fine-mapping of the ms-30 and ms-33 loci, 76 and 39 recombinants were selected from each F2 population using markers HP149 and HP3661, respectively. These recombinants were identified as homozygous for either the ms-30 or ms-33 mutant allele and for only one of these two markers. The recombinants were further genotyped with additional markers (Table S12), and their fertility/sterility was assessed by examining the flowers on at least three inflorescences per plant.

4.5. Gene Prediction and Sequence Polymorphism Analysis

The putative genes in the fine-mapped region containing the ms-30 and ms-33 loci were identified using the tomato gene model (ITAG release 4.0) [37] in the SGN (https://solgenomics.net).
The genomic sequence of TPI from the tomato lines WT-30, ms-30, WT-33, and ms-33 was obtained using overlapping PCR. The last exon and part of the 3′ UTR region of the TPI gene in 7B-1 and sl-2 were amplified by PCR with an elongation time of 7 min, using primers LeGLO2-6F/R [25], and 2× Phanta Max Master Mix (Dye Plus) (Cat. No. P525-03, Vazyme, Nanjing, China). The coding sequences (CDS) of the TPI gene were obtained through reverse-transcription PCR (RT-PCR). Fragments were amplified using primers LeGLO2-6F/R, and the CDS of TPI were cloned into the pEASY-Blunt Zero Cloning Vector (pEASY-Blunt Zero Cloning Kit, Cat. No. CB501-02; TransGen Biotech, Beijing, China). Both the amplified fragments and plasmids were sequenced at the Beijing Genomics Institute (Beijing, China). All primers utilized in this study are detailed in Table S12. Genomic DNA of 7B-1 was provided by Professor Huolin Shen from the College of Horticulture, China Agricultural University.

4.6. RNA Extraction, cDNA Synthesis, and qRT-PCR Analysis

The sepals, petals, stamens, and carpels of young flower buds at the meiotic and tetrad stages (3–5 mm in length) [40] were collected from WT-30, ms-30, WT-33, ms-33, AC, tpicr-e1, and tpicr-e7 plants and then immediately frozen in liquid nitrogen. Each tissue type included three biological replicates, with each replicate comprising samples from at least three plants. Total RNA extraction, cDNA synthesis, and qRT-PCR analysis were conducted following the methods previously described by Cao et al. (2019) [26] and Liu et al. (2019) [10]. All primers used for these analyses are detailed in Table S12.

4.7. RNA-Seq Analysis

A total of 24 RNA-Seq libraries, constructed with RNA from WT-30 and ms-30, were sequenced using the Illumina NovaSeq6000 platform at Berry Genomics Corporation (Beijing, China). The clean reads obtained were aligned and quantified against the tomato genome (ITAG4.0) using HISAT2 version 2.1.0 [59] and StringTie version 2.0.6 [60]. Differential gene expression analysis was conducted using the DESeq2 package, employing the MA-plot-based method in R version 3.0.3. p values were adjusted using the Benjamini–Hochberg procedure to control the false discovery rate [61]. The fold change for the sepals of WT-30 (SE_WT-30) versus sepals of ms-30 (SE_ms-30), petals of WT-30 (PE_WT-30) versus petals of ms-30 (PE_ms-30), stamens of WT-30 (ST_WT-30) versus stamens of ms-30 (ST_ms-30), and carpels of WT-30 (CA_WT-30) versus carpels of ms-30 (CA_ms-30) was determined using FPKM (fragments per kilobase of transcript sequence per million base pairs sequenced) values. The criteria for identifying DEGs were set as |log2FoldChange| ≥ 1 and the adjusted p value < 0.05. GO analysis of the DEGs was performed using their closest homologs in Arabidopsis via DAVID (Database for Annotation, Visualization, and Integrated Discovery, https://david.ncifcrf.gov/). The expression levels of certain DEGs, previously reported to be involved in the development of stamens or carpels, were validated through qRT-PCR.

4.8. Construction of Gene Editing Vectors, Plant Transformation, and Selection of Mutant Alleles

For the CRISPR/Cas9 constructs, two sgRNA binding sites per vector were identified using CRISPR-GE [62] (http://skl.scau.edu.cn/). Primers incorporating sgRNAs and Bsa I recognition sites were used to amplify the sgRNA(1)_SlU6-2t_SlU3-5p_sgRNA(2) fragments, with the pCBC-S1 vector as a template. The PCR fragments were then assembled into the binary vector pMGET through the Golden Gate cloning method, following the protocols described by Xing et al. (2014) [63] and Yang et al. (2023) [64]. These constructs were introduced into AC tomato plants via Agrobacterium tumefaciens (GV3101)-mediated cotyledon explant transformation, adhering to the method outlined by Yang et al. (2023) [64]. To confirm whether transgenic plants were successfully generated, PCR was performed using primers Cas9F and Cas9R. To identify specific types of edits, two primer pairs, TPI1 and TPI2, were used to amplify sequences including the targets of vector1 and vector2, respectively. The resulting PCR fragments were cloned into the pEASY-Blunt Zero Cloning Vector. At least 12 clones per PCR product were sequenced at the Beijing Genomics Institute. Two heterozygous T0 plants, each with a 4 bp or 66 bp deletion, were self-pollinated to produce T1 generation seeds. Homozygous mutants from the T1 generation, lacking CAS9, were identified using primers TPI1, TPI2, and Cas9. All primers used are detailed in Table S12.

4.9. Data Statistical Analysis

The mean values of floral organs length and relative expression were analyzed using IBM SPSS Statistics 20.0 software, with a t-test employed to assess significant differences among the means. The chi-square goodness-of-fit test was conducted by online data analysis platform SPSSPRO (https://www.spsspro.com).

5. Conclusions

Overall, this study demonstrated that the TPI gene (Solyc06g059970) is the causal factor for the phenotypic abnormalities observed in the mutants ms-30, ms-33, 7B-1, and sl-2. Sequence alterations in TPI across these mutants included an SNP each in ms-30 and ms-33 and an LTR retrotransposon insertion in both 7B-1 and sl-2. Leveraging these sequence variations, three codominant markers—MS30KASP, MS33SEQ, and LeGLO2-6/SL-2TE—were developed. RNA-seq and qRT-PCR analyses indicated that loss-of-function mutations in TPI altered the expression of genes involved in stamen and carpel development. These findings may be useful for the rapid development of male sterile lines through molecular marker-assisted backcrossing or CRISPR/Cas9-mediated mutagenesis targeting TPI, contributing to hybrid seed production and laying the groundwork for further functional analysis of the TPI gene.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25063331/s1. References [65,66,67,68,69] are cited in the Supplementary Material.

Author Contributions

Conceptualization, Z.H.; methodology, K.W., X.C. and S.L.; validation, formal analysis, investigation, K.W., X.L., X.C., S.L., L.Z., F.L. and C.L.; writing—original draft preparation, K.W. and Z.H.; writing—review and editing, K.W., X.L., X.C., W.Y., Z.H. and X.W.; visualization, K.W.; supervision, X.L., Y.G., L.L., C.Z., Y.D., J.L., W.Y., Z.H. and X.W.; project administration, X.W. and Z.H.; funding acquisition, X.L., X.W. and Z.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the National Key Research and Development Program of China (No. 2022YFF1003000), the National Natural Science Foundation of China (Nos. 31872949 and 31672154), and the China Agriculture Research System (No. CARS-23-A06).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The RNA sequencing datasets generated in this study have been deposited in the Sequence Read Archive (SRA) under the accession number PRJNA1036293.

Acknowledgments

We thank the Tomato Genetics Resource Center (Davis, CA, USA) for providing seed stocks and Professor Huolin Shen (China Agricultural University) for giving the tomato genome DNA of 7B-1.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bai, Y.; Lindhout, P. Domestication and breeding of tomatoes: What have we gained and what can we gain in the future? Ann. Bot. 2007, 100, 1085–1094. [Google Scholar] [CrossRef]
  2. Tamta, S.; Singh, J.P. Heterosis in tomato for growth and yield traits. Int. J. Veg. Sci. 2017, 24, 169–179. [Google Scholar] [CrossRef]
  3. Cheema, D.S.; Dhaliwal, M.S. Hybrid tomato breeding. J. New Seeds 2004, 6, 1–14. [Google Scholar] [CrossRef]
  4. Atanassova, B. Genic male sterility and its application in tomato (Lycopersicon esculentum Mill.) hybrid breeding and hybrid seed production. Acta. Hortic. 2007, 729, 45–51. [Google Scholar] [CrossRef]
  5. Cheng, Z.; Song, W.; Zhang, X. Genic male and female sterility in vegetable crops. Hortic. Res. 2023, 10, uhac232. [Google Scholar] [CrossRef] [PubMed]
  6. Kaul, M.L.H. Male Sterility in Higher Plants; Springer: Berlin/Heidelberg, Germany, 1988. [Google Scholar]
  7. Gorman, S.W.; McCormick, S. Male sterility in tomato. Crit. Rev. Plant. Sci. 1997, 16, 31–53. [Google Scholar] [CrossRef]
  8. Jeong, H.J.; Kang, J.H.; Zhao, M.; Kwon, J.K.; Choi, H.S.; Bae, J.H.; Lee, H.A.; Joung, Y.H.; Choi, D.; Kang, B.C. Tomato male sterile 1035 is essential for pollen development and meiosis in anthers. J. Exp. Bot. 2014, 65, 6693–6709. [Google Scholar] [CrossRef]
  9. Jung, Y.J.; Kim, D.H.; Lee, H.J.; Nam, K.H.; Bae, S.; Nou, I.S.; Cho, Y.G.; Kim, M.K.; Kang, K.K. Knockout of SlMS10 gene (Solyc02g079810) encoding bHLH transcription factor using CRISPR/Cas9 system confers male sterility phenotype in tomato. Plants 2020, 9, 1189. [Google Scholar] [CrossRef] [PubMed]
  10. Liu, X.; Yang, M.; Liu, X.; Wei, K.; Cao, X.; Wang, X.; Wang, X.; Guo, Y.; Du, Y.; Li, J.; et al. A putative bHLH transcription factor is a candidate gene for male sterile 32, a locus affecting pollen and tapetum development in tomato. Hortic. Res. 2019, 6, 88. [Google Scholar] [CrossRef]
  11. Zhu, E.; You, C.; Wang, S.; Cui, J.; Niu, B.; Wang, Y.; Qi, J.; Ma, H.; Chang, F. The DYT1-interacting proteins bHLH010, bHLH089 and bHLH091 are redundantly required for Arabidopsis anther development and transcriptome. Plant J. 2015, 83, 976–990. [Google Scholar] [CrossRef]
  12. Roever, W.E. A promising type of male sterility for use in hybrid tomato seed production. Science 1948, 107, 506. [Google Scholar] [CrossRef]
  13. Larson, R.; Paur, S. The description and inheritance of a functionally sterile flower mutant in tomato and its probable value in hybrid tomato seed production. Proc. Am. Soc. Hort. Sci. 1948, 52, 355–364. [Google Scholar]
  14. Gorguet, B.; Schipper, D.; van Heusden, A.W.; Lindhout, P. High-resolution fine mapping of ps-2, a mutated gene conferring functional male sterility in tomato due to non-dehiscent anthers. Theor. Appl. Genet. 2006, 113, 1437–1448. [Google Scholar] [CrossRef] [PubMed]
  15. Gorguet, B.; Schipper, D.; van Lammeren, A.; Visser, R.G.; van Heusden, A.W. ps-2, the gene responsible for functional sterility in tomato, due to non-dehiscent anthers, is the result of a mutation in a novel polygalacturonase gene. Theor. Appl. Genet. 2009, 118, 1199–1209. [Google Scholar] [CrossRef] [PubMed]
  16. Atanassova, B. Functional male sterility (ps-2) in tomato (Lycopesicon esculentum Mill.) and its application in breeding and hybrid seed production. Euphytica 1999, 107, 13–21. [Google Scholar] [CrossRef]
  17. Rick, C.M.; Robinson, J. Inherited defects of floral structure affecting fruitfulness in Lycopersicon esculentum. Am. J. Bot. 1951, 38, 639–652. [Google Scholar] [CrossRef]
  18. Chen, K.Y.; Cong, B.; Wing, R.; Vrebalov, J.; Tanksley, S.D. Changes in regulation of a transcription factor lead to autogamy in cultivated tomatoes. Science 2007, 318, 643–645. [Google Scholar] [CrossRef] [PubMed]
  19. Cheng, M.Z.; Gong, C.; Zhang, B.; Qu, W.; Qi, H.N.; Chen, X.L.; Wang, X.Y.; Zhang, Y.; Liu, J.Y.; Ding, X.D.; et al. Morphological and anatomical characteristics of exserted stigma sterility and the location and function of SlLst (Solanum lycopersicum Long styles) gene in tomato. Theor. Appl. Genet. 2021, 134, 505–518. [Google Scholar] [CrossRef] [PubMed]
  20. Shang, L.; Song, J.; Yu, H.; Wang, X.; Yu, C.; Wang, Y.; Li, F.; Lu, Y.; Wang, T.; Ouyang, B.; et al. A mutation in a C2H2-type zinc finger transcription factor contributed to the transition toward self-pollination in cultivated tomato. Plant Cell 2021, 33, 3293–3308. [Google Scholar] [CrossRef]
  21. Bishop, C.J. A stamenless male-sterile tomato. Am. J. Bot. 1954, 41, 540–542. [Google Scholar] [CrossRef]
  22. Gómez, P.; Jamilena, M.; Capel, J.; Zurita, S.; Angosto, T.; Lozano, R. Stamenless, a tomato mutant with homeotic conversions in petals and stamens. Planta 1999, 209, 172–179. [Google Scholar] [CrossRef]
  23. Quinet, M.; Bataille, G.; Dobrev, P.I.; Capel, C.; Gomez, P.; Capel, J.; Lutts, S.; Motyka, V.; Angosto, T.; Lozano, R. Transcriptional and hormonal regulation of petal and stamen development by STAMENLESS, the tomato (Solanum lycopersicum L.) orthologue to the B-class APETALA3 gene. J. Exp. Bot. 2014, 65, 2243–2256. [Google Scholar] [CrossRef]
  24. Sawhney, V.K.; Greyson, R.I. Morphogenesis of the stamenless-2 mutant in tomato. I. Comparative description of the flowers and ontogeny of stamens in the normal and mutant plants. Am. J. Bot. 1973, 60, 514–523. [Google Scholar]
  25. Pucci, A.; Picarella, M.E.; Mazzucato, A. Phenotypic, genetic and molecular characterization of 7B-1, a conditional male-sterile mutant in tomato. Theor. Appl. Genet. 2017, 130, 2361–2374. [Google Scholar] [CrossRef]
  26. Cao, X.; Liu, X.; Wang, X.; Yang, M.; van Giang, T.; Wang, J.; Liu, X.; Sun, S.; Wei, K.; Wang, X.; et al. B-class MADS-box TM6 is a candidate gene for tomato male sterile-1526. Theor. Appl. Genet. 2019, 132, 2125–2135. [Google Scholar] [CrossRef] [PubMed]
  27. Fonseca, R.; Capel, C.; Lebrón, R.; Ortiz-Atienza, A.; Yuste-Lisbona, F.J.; Angosto, T.; Capel, J.; Lozano, R. Insights into the functional role of tomato TM6 as transcriptional regulator of flower development. Hortic. Res. 2024, uhae019. [Google Scholar] [CrossRef]
  28. Atanassova, B.; Georgiev, H. Investigation of tomato male sterile lines in relation to hybrid seedproduction. Acta. Hortic. 1986, 190, 553–557. [Google Scholar] [CrossRef]
  29. Sawhney, V. Photoperiod-sensitive male-sterile mutant in tomato and its potential use in hybrid seed production. J. Hortic. Sci. Biotech. 2004, 79, 138–141. [Google Scholar] [CrossRef]
  30. Rick, C.M. Tomato Genetics Cooperative Report. 1960. No. 10. pp. 35–37. Available online: https://tgc.ifas.ufl.edu/onlinevo.htm (accessed on 15 February 2024).
  31. Rick, C.M. Tomato Genetics Cooperative Report. 1962. No. 12. pp. 39–41. Available online: https://tgc.ifas.ufl.edu/onlinevo.htm (accessed on 15 February 2024).
  32. Lona, J.L. Tomato Genetics Cooperative Report. 1966. No. 16. p. 17. Available online: https://tgc.ifas.ufl.edu/onlinevo.htm (accessed on 15 February 2024).
  33. Andrasfalvy, A. Tomato Genetics Cooperative Report. 1970. No. 20. p. 12. Available online: https://tgc.ifas.ufl.edu/onlinevo.htm (accessed on 15 February 2024).
  34. Schmidt, H.; Schmidt, V. Untersuchungen an pollensterilen, stamenless-aehnlichen mutanten von Lycopersicon esculentum Mill. II: Normalisierung von ms-15 und ms-33 mit gibberellinsaeure (GA3). Biol. Zentralbl. 1981, 10, 691–696. [Google Scholar]
  35. Tanksley, S.D.; Ganal, M.W.; Prince, J.P.; de Vicente, M.C.; Bonierbale, M.W.; Broun, P.; Fulton, T.M.; Giovannoni, J.J.; Grandillo, S.; Martin, G.B. High density molecular linkage maps of the tomato and potato genomes. Genetics 1992, 132, 1141–1160. [Google Scholar] [CrossRef] [PubMed]
  36. Weide, R.; Van Wordragen, M.F.; Lankhorst, R.K.; Verkerk, R.; Hanhart, C.; Liharska, T.; Pap, E.; Stam, P.; Zabel, P.; Koornneef, M. Integration of the classical and molecular linkage maps of tomato chromosome 6. Genetics 1993, 135, 1175–1186. [Google Scholar] [CrossRef]
  37. Hosmani, P.S.; Flores-Gonzalez, M.; Geest, H.v.d.; Maumus, F.; Bakker, L.V.; Schijlen, E.; Haarst, J.v.; Cordewener, J.; Sanchez-Perez, G.; Peters, S.; et al. An improved de novo assembly and annotation of the tomato reference genome using single-molecule sequencing, Hi-C proximity ligation and optical maps. bioRxiv 2019, 767764. [Google Scholar] [CrossRef]
  38. de Martino, G.; Pan, I.; Emmanuel, E.; Levy, A.; Irish, V.F. Functional analyses of two tomato APETALA3 genes demonstrate diversification in their roles in regulating floral development. Plant Cell 2006, 18, 1833–1845. [Google Scholar] [CrossRef] [PubMed]
  39. Mazzucato, A.; Olimpieri, I.; Siligato, F.; Picarella, M.E.; Soressi, G.P. Characterization of genes controlling stamen identity and development in a parthenocarpic tomato mutant indicates a role for the DEFICIENS ortholog in the control of fruit set. Physiol. Plantarum 2008, 132, 526–537. [Google Scholar] [CrossRef]
  40. Brukhin, V.; Hernould, M.; Gonzalez, N.; Chevalier, C.; Mouras, A. Flower development schedule in tomato Lycopersicon esculentum cv. sweet cherry. Sex. Plant Reprod. 2003, 15, 311–320. [Google Scholar] [CrossRef]
  41. Bao, H.; Ding, Y.; Yang, F.; Zhang, J.; Xie, J.; Zhao, C.; Du, K.; Zeng, Y.; Zhao, K.; Li, Z.; et al. Gene silencing, knockout and over-expression of a transcription factor ABORTED MICROSPORES (SlAMS) strongly affects pollen viability in tomato (Solanum lycopersicum). BMC Genom. 2022, 23, 346. [Google Scholar] [CrossRef]
  42. Gökdemir, G.; Seçgin, Z.; Uluisik, S.; Kavas, M. CRISPR/Cas9 knock-out of SlPHD_MS1 (Solyc04g008420) gene results in complete male sterility in tomato. Plant Growth Regul. 2022, 98, 329–341. [Google Scholar] [CrossRef]
  43. David-Schwartz, R.; Koenig, D.; Sinha, N.R. LYRATE is a key regulator of leaflet initiation and lamina outgrowth in tomato. Plant Cell 2009, 21, 3093–3104. [Google Scholar] [CrossRef]
  44. Skinner, D.J.; Brown, R.H.; Kuzoff, R.K.; Gasser, C.S. Conservation of the role of INNER NO OUTER in development of unitegmic ovules of the Solanaceae despite a divergence in protein function. BMC Plant Biol. 2016, 16, 143. [Google Scholar] [CrossRef]
  45. Zhang, J.; Wang, Y.; Naeem, M.; Zhu, M.; Li, J.; Yu, X.; Hu, Z.; Chen, G. An AGAMOUS MADS-box protein, SlMBP3, regulates the speed of placenta liquefaction and controls seed formation in tomato. J. Exp. Bot. 2019, 70, 909–924. [Google Scholar] [CrossRef]
  46. Castañeda, L.; Giménez, E.; Pineda, B.; García-Sogo, B.; Ortiz-Atienza, A.; Micol-Ponce, R.; Angosto, T.; Capel, J.; Moreno, V.; Yuste-Lisbona, F.J.; et al. Tomato CRABS CLAW paralogues interact with chromatin remodelling factors to mediate carpel development and floral determinacy. New Phytol. 2022, 234, 1059–1074. [Google Scholar] [CrossRef] [PubMed]
  47. Xu, Y.; Crouch, J.H. Marker-assisted selection in plant breeding: From publications to practice. Crop Sci. 2008, 48, 391–407. [Google Scholar] [CrossRef]
  48. Du, M.; Zhou, K.; Liu, Y.; Deng, L.; Zhang, X.; Lin, L.; Zhou, M.; Zhao, W.; Wen, C.; Xing, J.; et al. A biotechnology-based male-sterility system for hybrid seed production in tomato. Plant J. 2020, 102, 1090–1100. [Google Scholar] [CrossRef] [PubMed]
  49. Dai, X.; Han, H.; Huang, W.; Zhao, L.; Song, M.; Cao, X.; Liu, C.; Niu, X.; Lang, Z.; Ma, C.; et al. Generating novel male sterile tomatoes by editing respiratory burst oxidase homolog genes. Front. Plant Sci. 2021, 12, 817101. [Google Scholar] [CrossRef]
  50. Liu, J.; Wang, S.; Wang, H.; Luo, B.; Cai, Y.; Li, X.; Zhang, Y.; Wang, X. Rapid generation of tomato male-sterile lines with a marker use for hybrid seed production by CRISPR/Cas9 system. Mol. Breed. 2021, 41, 25. [Google Scholar] [CrossRef]
  51. Zhou, M.; Deng, L.; Yuan, G.; Zhao, W.; Ma, M.; Sun, C.; Du, M.; Li, C.; Li, C. A CRISPR-Cas9-Derived male sterility system for tomato breeding. Agronomy 2023, 13, 1785. [Google Scholar] [CrossRef]
  52. Zhou, M.; Deng, L.; Yuan, G.; Zhao, W.; Ma, M.; Sun, C.; Du, M.; Li, C.; Li, C. Rapid generation of a tomato male sterility system and its feasible application in hybrid seed production. Theor. Appl. Genet. 2023, 136, 197. [Google Scholar] [CrossRef] [PubMed]
  53. Riechmann, J.L.; Ferrándiz, C. Flower Development: Methods and Protocols; Springer Nature: New York, NY, USA, 2023. [Google Scholar]
  54. Geuten, K.; Irish, V. Hidden variability of floral homeotic B genes in Solanaceae provides a molecular basis for the evolution of novel functions. Plant Cell 2010, 22, 2562–2578. [Google Scholar] [CrossRef] [PubMed]
  55. Schwarz-Sommer, Z.; Huijser, P.; Nacken, W.; Saedler, H.; Sommer, H. Genetic control of flower development by homeotic genes in Antirrhinum majus. Science 1990, 250, 931–936. [Google Scholar] [CrossRef] [PubMed]
  56. Coen, E.S.; Meyerowitz, E.M. The war of the whorls: Genetic interactions controlling flower development. Nature 1991, 353, 31–37. [Google Scholar] [CrossRef]
  57. Theißen, G.; Saedler, H. Floral quartets. Nature 2001, 409, 469–471. [Google Scholar] [CrossRef] [PubMed]
  58. Fernandez-Pozo, N.; Menda, N.; Edwards, J.D.; Saha, S.; Tecle, I.Y.; Strickler, S.R.; Bombarely, A.; Fisher-York, T.; Pujar, A.; Foerster, H.; et al. The Sol Genomics Network (SGN)—From genotype to phenotype to breeding. Nucleic Acids Res. 2015, 43, D1036–D1041. [Google Scholar] [CrossRef] [PubMed]
  59. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar] [CrossRef]
  60. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef]
  61. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  62. Xie, X.; Ma, X.; Zhu, Q.; Zeng, D.; Li, G.; Liu, Y.G. CRISPR-GE: A convenient software toolkit for CRISPR-based genome editing. Mol. Plant 2017, 10, 1246–1249. [Google Scholar] [CrossRef] [PubMed]
  63. Xing, H.L.; Dong, L.; Wang, Z.P.; Zhang, H.Y.; Han, C.Y.; Liu, B.; Wang, X.C.; Chen, Q.J. A CRISPR/Cas9 toolkit for multiplex genome editing in plants. BMC Plant Biol. 2014, 14, 327. [Google Scholar] [CrossRef]
  64. Yang, M.X.; Liu, X.L.; Cao, C.; Wei, K.; Ning, Y.; Yang, P.; Li, S.S.; Chen, Z.Y.; Wang, X.X.; Guo, Y.M.; et al. Construction and application of a CRISPR/Cas9 system for multiplex gene editing in tomato. Acta. Hortic. Sin. 2023, 50, 1215–1229. (In Chinese) [Google Scholar] [CrossRef]
  65. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  66. Nicholas, K.B.; Nicholas, H.B. GeneDoc: A Tool for Editing and Annotating Multiple Sequence Alignments. Ver. 2.7.000. 1997. Available online: https://nrbsc.org/gfx/genedoc (accessed on 14 February 2024).
  67. Guo, X.; Hu, Z.; Yin, W.; Yu, X.; Zhu, Z.; Zhang, J.; Chen, G. The tomato floral homeotic protein FBP1-like gene, SlGLO1, plays key roles in petal and stamen development. Sci. Rep. 2016, 6, 20454. [Google Scholar] [CrossRef]
  68. Pan, I.L.; McQuinn, R.; Giovannoni, J.J.; Irish, V.F. Functional diversification of AGAMOUS lineage genes in regulating tomato flower and fruit development. J. Exp. Bot. 2010, 61, 1795–1806. [Google Scholar] [CrossRef] [PubMed]
  69. Qiu, Z.; Wang, X.; Gao, J.; Guo, Y.; Huang, Z.; Du, Y. The tomato hoffman’s anthocyaninless gene encodes a bHLH transcription factor involved in anthocyanin biosynthesis that is developmentally regulated and induced by low temperatures. PLoS ONE 2016, 11, e0151067. [Google Scholar] [CrossRef] [PubMed]
Ijms 25 03331 g001
Figure 1. Phenotypes of flowers and stamens in ms-30 and ms-33. (A) Flower, (B) anther cone, and (C) dissected stamen in WT-30. (D) Flower and (E) anther cone in ms-30. (F) Dissected stamen with carpelloid structures in ms-30. (G) Dissected stamen with carpelloid structures and external ovules in ms-30; white arrow denotes the external ovules. (H) Completely carpelloid stamen in ms-30. (I) Flower, (J) anther cone, (K) and dissected stamen in WT-33. (L) Flower and (M) anther cone in ms-33. (N) Dissected stamen with carpelloid structures in ms-33. (O) Dissected stamen with carpelloid structures and external ovules in ms-33; white arrow indicates the external ovules. (P) Completely carpelloid stamen in ms-33. The white bar indicates 0.5 cm.
Figure 1. Phenotypes of flowers and stamens in ms-30 and ms-33. (A) Flower, (B) anther cone, and (C) dissected stamen in WT-30. (D) Flower and (E) anther cone in ms-30. (F) Dissected stamen with carpelloid structures in ms-30. (G) Dissected stamen with carpelloid structures and external ovules in ms-30; white arrow denotes the external ovules. (H) Completely carpelloid stamen in ms-30. (I) Flower, (J) anther cone, (K) and dissected stamen in WT-33. (L) Flower and (M) anther cone in ms-33. (N) Dissected stamen with carpelloid structures in ms-33. (O) Dissected stamen with carpelloid structures and external ovules in ms-33; white arrow indicates the external ovules. (P) Completely carpelloid stamen in ms-33. The white bar indicates 0.5 cm.
Ijms 25 03331 g001
Ijms 25 03331 g002
Figure 2. Fine-mapping of ms-30. (A) Preliminary mapping of ms-30; PS indicates population size. (B) Fine-mapping of ms-30; N indicates the number of recombinants, and 76 recombinants were selected from 949 plants of F2 population. (C) ITAG4.0 annotated genes of ms-30. Arrows indicate the direction of transcription, and the solid arrow represents the most likely candidate gene for ms-30.
Figure 2. Fine-mapping of ms-30. (A) Preliminary mapping of ms-30; PS indicates population size. (B) Fine-mapping of ms-30; N indicates the number of recombinants, and 76 recombinants were selected from 949 plants of F2 population. (C) ITAG4.0 annotated genes of ms-30. Arrows indicate the direction of transcription, and the solid arrow represents the most likely candidate gene for ms-30.
Ijms 25 03331 g002
Ijms 25 03331 g003
Figure 3. Sequence alterations of the TPI gene in ms-30, ms-33, 7B-1, and sl-2 and the development of molecular markers (A,B) Single-nucleotide polymorphisms (SNPs) in the TPI gene of ms-30 and ms-33; black boxes indicate the exon of TPI; (C) structure of the TPI gene according to ITAG4.0 (upper) and the structure of the long terminal repeat (LTR) retrotransposon inserted in the last exon of TPI in 7B-1 and sl-2 (lower); gray boxes indicate the LTR sequence of retrotransposon, hollow blocks are two ORFs predicted on the retrotransposon, and arrows indicate the direction of transcription. “TATC” is the TSD sequences. LeGLO2-6F, LeGLO2-6R, SL-2TE-F, and SL-2TE-R are the primers used in generating the LeGLO2-6/SL-2TE marker. (D) Endpoint fluorescence scatter plot of the Kompetitive Allele Specific PCR marker MS30KASP. (E) The chromatogram of MS33SEQ; red box denotes the SNP in TPI of ms-33. (F) Agarose gel electrophoresis of PCR fragments amplified from Ailsa Craig, 7B-1, and plants from LA1801 (sl-2) using the marker LeGLO2-6/SL-2TE; F indicates the fertile plant, whereas S indicates the sterile plant.
Figure 3. Sequence alterations of the TPI gene in ms-30, ms-33, 7B-1, and sl-2 and the development of molecular markers (A,B) Single-nucleotide polymorphisms (SNPs) in the TPI gene of ms-30 and ms-33; black boxes indicate the exon of TPI; (C) structure of the TPI gene according to ITAG4.0 (upper) and the structure of the long terminal repeat (LTR) retrotransposon inserted in the last exon of TPI in 7B-1 and sl-2 (lower); gray boxes indicate the LTR sequence of retrotransposon, hollow blocks are two ORFs predicted on the retrotransposon, and arrows indicate the direction of transcription. “TATC” is the TSD sequences. LeGLO2-6F, LeGLO2-6R, SL-2TE-F, and SL-2TE-R are the primers used in generating the LeGLO2-6/SL-2TE marker. (D) Endpoint fluorescence scatter plot of the Kompetitive Allele Specific PCR marker MS30KASP. (E) The chromatogram of MS33SEQ; red box denotes the SNP in TPI of ms-33. (F) Agarose gel electrophoresis of PCR fragments amplified from Ailsa Craig, 7B-1, and plants from LA1801 (sl-2) using the marker LeGLO2-6/SL-2TE; F indicates the fertile plant, whereas S indicates the sterile plant.
Ijms 25 03331 g003
Ijms 25 03331 g004
Figure 4. Sequence alterations and phenotypes of TPI gene-edited mutants. (A) TPI mutations generated through CRISPR/Cas9 gene editing. Red letters mark the sequence of targets. Blue letters indicate the protospacer-adjacent motif. Green hyphens indicate deletions in tpicr mutants. (B–J) Representative flower, anther cone, and pollen viability images of AC, tpicr-e1, and tpicr-e7. The white bar indicates 0.5 cm (B–G), and the black bar indicates 1 mm (H–J).
Figure 4. Sequence alterations and phenotypes of TPI gene-edited mutants. (A) TPI mutations generated through CRISPR/Cas9 gene editing. Red letters mark the sequence of targets. Blue letters indicate the protospacer-adjacent motif. Green hyphens indicate deletions in tpicr mutants. (B–J) Representative flower, anther cone, and pollen viability images of AC, tpicr-e1, and tpicr-e7. The white bar indicates 0.5 cm (B–G), and the black bar indicates 1 mm (H–J).
Ijms 25 03331 g004
Ijms 25 03331 g005
Figure 5. Loss of function of TPI in ms-30 affects multiple biological processes during stamen development. (A) Principal component analysis of RNA-seq data. (B) The number of differentially expressed genes that were upregulated and downregulated in SE, PE, ST, and CA. (C) Gene ontology (GO) enrichment analysis of upregulated differentially expressed genes in ST. (D) GO enrichment analysis of downregulated differentially expressed genes in ST. SE, Sepals; PE, Petals; ST, Stamens; and CA, Carpels.
Figure 5. Loss of function of TPI in ms-30 affects multiple biological processes during stamen development. (A) Principal component analysis of RNA-seq data. (B) The number of differentially expressed genes that were upregulated and downregulated in SE, PE, ST, and CA. (C) Gene ontology (GO) enrichment analysis of upregulated differentially expressed genes in ST. (D) GO enrichment analysis of downregulated differentially expressed genes in ST. SE, Sepals; PE, Petals; ST, Stamens; and CA, Carpels.
Ijms 25 03331 g005
Ijms 25 03331 g006
Figure 6. Expression pattern of the DEGs related to stamen or carpel development between WT-30 and ms-30. (A) Heatmap of the DEGs whose homologs in Arabidopsis are related to stamen development. The genes in the red box represent those expressed mainly in the stamen of WT-30 and significantly downregulated in the stamen of ms-30. (B) Heatmap of DEGs whose homologs in Arabidopsis are related to carpel development. Genes in the red box represent those mainly expressed in the carpel and significantly upregulated in the stamen of ms-30. (C,D) FPKM and qRT-PCR analysis of the DEGs related to stamen and carpel development in four floral organs of ms-30. Asterisks indicate a significant difference (**, p < 0.01) between WT-30 and ms-30; WT-30 and ms-30 are homozygous male fertile plants and male sterile plants derived from the tomato line 2-455, respectively. SE, sepals; PE, petals; ST, stamens; and CA, carpels.
Figure 6. Expression pattern of the DEGs related to stamen or carpel development between WT-30 and ms-30. (A) Heatmap of the DEGs whose homologs in Arabidopsis are related to stamen development. The genes in the red box represent those expressed mainly in the stamen of WT-30 and significantly downregulated in the stamen of ms-30. (B) Heatmap of DEGs whose homologs in Arabidopsis are related to carpel development. Genes in the red box represent those mainly expressed in the carpel and significantly upregulated in the stamen of ms-30. (C,D) FPKM and qRT-PCR analysis of the DEGs related to stamen and carpel development in four floral organs of ms-30. Asterisks indicate a significant difference (**, p < 0.01) between WT-30 and ms-30; WT-30 and ms-30 are homozygous male fertile plants and male sterile plants derived from the tomato line 2-455, respectively. SE, sepals; PE, petals; ST, stamens; and CA, carpels.
Ijms 25 03331 g006
Table 1. Chi-square goodness-of-fit test for segregation ratios of F2 populations at ms-30 and ms-33 loci.
Table 1. Chi-square goodness-of-fit test for segregation ratios of F2 populations at ms-30 and ms-33 loci.
LocusPhenotypeObserved
Frequency		Expected
Frequency		Observed
Proportion		Expected
Proportion		Residualχ2p
ms-30Fertile1391500.6950.75−113.2270.072
Sterile61500.3050.2511
ms-33Fertile1491500.7450.75−10.0270.870
Sterile51500.2550.251
Table 2. Predicted genes in the ms-30 region.
Table 2. Predicted genes in the ms-30 region.
Gene NamePosition on SL4.0ch06Putative Function
Solyc06g05993035503650..35543738(+)sesquiterpene synthase 1
Solyc06g05996035545393..35569077(−)Histone-lysine N-methyltransferase ASHH2
Solyc06g05997035575492..35578559(−)Tomato PISTILLATA (TPI, syn. SlGLO2)
Solyc06g05998035588500..35603986(+)O-fucosyltransferase family protein
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wei, K.; Li, X.; Cao, X.; Li, S.; Zhang, L.; Lu, F.; Liu, C.; Guo, Y.; Liu, L.; Zhu, C.; et al. Candidate Gene Identification and Transcriptome Analysis of Tomato male sterile-30 and Functional Marker Development for ms-30 and Its Alleles, ms-33, 7B-1, and stamenless-2. Int. J. Mol. Sci. 2024, 25, 3331. https://doi.org/10.3390/ijms25063331

AMA Style

Wei K, Li X, Cao X, Li S, Zhang L, Lu F, Liu C, Guo Y, Liu L, Zhu C, et al. Candidate Gene Identification and Transcriptome Analysis of Tomato male sterile-30 and Functional Marker Development for ms-30 and Its Alleles, ms-33, 7B-1, and stamenless-2. International Journal of Molecular Sciences. 2024; 25(6):3331. https://doi.org/10.3390/ijms25063331

Chicago/Turabian Style

Wei, Kai, Xin Li, Xue Cao, Shanshan Li, Li Zhang, Feifei Lu, Chang Liu, Yanmei Guo, Lei Liu, Can Zhu, and et al. 2024. "Candidate Gene Identification and Transcriptome Analysis of Tomato male sterile-30 and Functional Marker Development for ms-30 and Its Alleles, ms-33, 7B-1, and stamenless-2" International Journal of Molecular Sciences 25, no. 6: 3331. https://doi.org/10.3390/ijms25063331

APA Style

Wei, K., Li, X., Cao, X., Li, S., Zhang, L., Lu, F., Liu, C., Guo, Y., Liu, L., Zhu, C., Du, Y., Li, J., Yang, W., Huang, Z., & Wang, X. (2024). Candidate Gene Identification and Transcriptome Analysis of Tomato male sterile-30 and Functional Marker Development for ms-30 and Its Alleles, ms-33, 7B-1, and stamenless-2. International Journal of Molecular Sciences, 25(6), 3331. https://doi.org/10.3390/ijms25063331

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