Transcriptome Analysis and Metabolic Profiling Reveal the Key Regulatory Pathways in Drought Stress Responses and Recovery in Tomatoes.

Drought stress is a major abiotic factor affecting tomato production and fruit quality. However, the genes and metabolites associated with tomato responses to water deficiency and rehydration are poorly characterized. To identify the functional genes and key metabolic pathways underlying tomato responses to drought stress and recovery, drought-susceptible and drought-tolerant inbred lines underwent transcriptomic and metabolomic analyses. A total of 332 drought-responsive and 491 rehydration-responsive core genes were robustly differentially expressed in both genotypes. The drought-responsive and rehydration-responsive genes were mainly related to photosynthesis-antenna proteins, nitrogen metabolism, plant-pathogen interactions, and the MAPK signaling pathway. Various transcription factors, including homeobox-leucine zipper protein ATHB-12, NAC transcription factor 29, and heat stress transcription factor A-6b-like, may be vital for tomato responses to water status. Moreover, 24,30-dihydroxy-12(13)-enolupinol, caffeoyl hawthorn acid, adenosine 5'-monophosphate, and guanosine were the key metabolites identified in both genotypes under drought and recovery conditions. The combined transcriptomic and metabolomic analysis highlighted the importance of 38 genes involved in metabolic pathways, the biosynthesis of secondary metabolites, the biosynthesis of amino acids, and ABC transporters for tomato responses to water stress. Our results provide valuable clues regarding the molecular basis of drought tolerance and rehydration. The data presented herein may be relevant for genetically improving tomatoes to enhance drought tolerance.

OBV (obscure vein), a C2H2 zinc finger transcription factor, positively regulates chloroplast development and bundle sheath extension formation in tomato (Solanum lycopersicum) leaf veins.

Leaf veins play an important role in plant growth and development, and the bundle sheath (BS) is believed to greatly improve the photosynthetic efficiency of C4 plants. The OBV mutation in tomato (Solanum lycopersicum) results in dark veins and has been used widely in processing tomato varieties. However, physiological performance has difficulty explaining fitness in production. In this study, we confirmed that this mutation was caused by both the increased chlorophyll content and the absence of bundle sheath extension (BSE) in the veins. Using genome-wide association analysis and map-based cloning, we revealed that OBV encoded a C2H2L domain class transcription factor. It was localized in the nucleus and presented cell type-specific gene expression in the leaf veins. Furthermore, we verified the gene function by generating CRISPR/Cas9 knockout and overexpression mutants of the tomato gene. RNA sequencing analysis revealed that OBV was involved in regulating chloroplast development and photosynthesis, which greatly supported the change in chlorophyll content by mutation. Taken together, these findings demonstrated that OBV affected the growth and development of tomato by regulating chloroplast development in leaf veins. This study also provides a solid foundation to further decipher the mechanism of BSEs and to understand the evolution of photosynthesis in land plants.

Fine Mapping of the Ph-2 Gene Conferring Resistance to Late Blight (Phytophthora infestans) in Tomato.

Late blight is a devastating tomato disease. Breeding new varieties with multiple resistance (R) genes is highly effective for preventing late blight. The Ph-2 gene mediates resistance to Phytophthora infestans race T1 in tomato. In this study, we used an F2 population derived from a cross between Solanum lycopersicum Moboline (resistant) and LA3988 (susceptible) cultivars for the fine mapping of Ph-2. Two flanking markers, CAPS-1 and CC-Ase, mapped Ph-2 to a 141-kb genomic region containing 21 projected genes, 5 of which were identified as putative R genes. The Solyc10g085460 coding sequence varied significantly between the parents. The markers developed and candidate genes identified in this study shall be useful for the molecular breeding of tomato exhibiting increased late blight resistance and for the cloning of the Ph-2 gene.

SlGID1a Is a Putative Candidate Gene for qtph1.1, a Major-Effect Quantitative Trait Locus Controlling Tomato Plant Height.

Plant height is an important agronomic trait in crops. Several genes underlying tomato (Solanum lycopersicum) plant height mutants have been cloned. However, few quantitative trait genes for plant height have been identified in tomato. In this study, seven quantitative trait loci (QTLs) controlling plant height were identified in tomato. Of which, qtph1.1 (QTL for tomato plant height 1.1), qtph3.1 and qtph12.1 were major QTLs and explained 15, 16, and 12% of phenotypic variation (R2), respectively. The qtph1.1 was further mapped to an 18.9-kb interval on chromosome 1. Based on the annotated tomato genome (version SL2.50, annotation ITAG2.40), Solyc01g098390 encoding GA receptor SlGID1a was the putative candidate gene. The SlGID1a gene underlying the qtph1.1 locus contained a single nucleotide polymorphism (SNP) that resulted in an amino acid alteration in protein sequence. The near-isogenic line containing the qtph1.1 locus (NIL-qtph1.1) exhibited shorter internode length and cell length than the wild type (NIL-WT). The dwarf phenotype of NIL-qtph1.1 could not be rescued by exogenous GA3 treatment. Transcriptome analysis and real-time quantitative reverse transcription PCR (qPCR) showed that several genes related to biosynthesis and signaling of GA and auxin were differentially expressed in stems between NIL-qtph1.1 and NIL-WT. These findings might pave the road for understanding the molecular regulation mechanism of tomato plant height.

Exploration of resistance to Phelipanche aegyptiaca in tomato.

BACKGROUND: Cultivated tomatoes are highly susceptible to the destructive parasite Phelipanche aegyptiaca. Wild relatives show the potential resistance for genetic improvement. However, their genetic and molecular mechanisms are still unknown. RESULTS: Among 50 wild tomato accessions were evaluated for resistance to P. aegyptiaca, most of the wild relatives exhibited varying degrees of resistance compared to the cultivars. Solanum pennellii LA0716 performed the most promising and solid resistance with very low infection by the broomrape. The resistance involved in LA0716 was further confirmed by cytological analysis, and explored by employing a permanent introgression line (IL) population. Thirteen putative quantitative trait loci (QTLs) conferring the different resistance traits were identified. They are located on chromosomes 1, 2, 3, 4, 6, 8 and 9. The most attractive QTLs are positioned in IL6-2 and overlap with IL6-3. Specially, IL6-2 showed the highest and most consistent resistance for multiple traits and explained the major phenotypic variation of LA0716. Analysis of candidate genes involved in these regions showed that Beta (Solyc06g074240) and P450 (Solyc06g073570, Solyc06g074180 and Solyc06g074420) genes are substantially related to the strigolactone (SL) pathway. Transcript analysis further demonstrated that both Solyc06g073570 and Solyc06g074180 might play an important role in the reduction of P. aegyptiaca infection. CONCLUSION: Germplasms resistant to P. aegyptiaca were found in wild tomato species. QTLs conferring P. aegyptiaca tolerance in LA0716 were identified. IL6-2 is identified as a prospective line possessing the major QTLs. The candidate genes would provide the availability to assist the introgression of the resistance in future breeding programmes. © 2020 Society of Chemical Industry.

A putative bHLH transcription factor is a candidate gene for male sterile 32, a locus affecting pollen and tapetum development in tomato.

The tomato (Solanum lycopersicum) male sterile 32 (ms32) mutant has been used in hybrid seed breeding programs largely because it produces no pollen and has exserted stigmas. In this study, histological examination of anthers revealed dysfunctional pollen and tapetum development in the ms32 mutant. The ms32 locus was fine mapped to a 28.5 kb interval that encoded four putative genes. Solyc01g081100, a homolog of Arabidopsis bHLH10/89/90 and rice EAT1, was proposed to be the candidate gene of MS32 because it contained a single nucleotide polymorphism (SNP) that led to the formation of a premature stop codon. A codominant derived cleaved amplified polymorphic sequence (dCAPS) marker, MS32D, was developed based on the SNP. Real-time quantitative reverse-transcription PCR showed that most of the genes, which were proposed to be involved in pollen and tapetum development in tomato, were downregulated in the ms32 mutant. These findings may aid in marker-assisted selection of ms32 in hybrid breeding programs and facilitate studies on the regulatory mechanisms of pollen and tapetum development in tomato.

Identification and Characterization of EI (Elongated Internode) Gene in Tomato (Solanum lycopersicum).

Internode length is an important agronomic trait affecting plant architecture and crop yield. However, few genes for internode elongation have been identified in tomato. In this study, we characterized an elongated internode inbred line P502, which is a natural mutant of the tomato cultivar 05T606. The mutant P502 exhibits longer internode and higher bioactive GA concentration compared with wild-type 05T606. Genetic analysis suggested that the elongated internode trait is controlled by quantitative trait loci (QTL). Then, we identified a major QTL on chromosome 2 based on molecular markers and bulked segregant analysis (BSA). The locus was designated as EI (Elongated Internode), which explained 73.6% genetic variance. The EI was further mapped to a 75.8-kb region containing 10 genes in the reference Heinz 1706 genome. One single nucleotide polymorphism (SNP) in the coding region of solyc02g080120.1 was identified, which encodes gibberellin 2-beta-dioxygenase 7 (SlGA2ox7). SlGA2ox7, orthologous to AtGA2ox7 and AtGA2ox8, is involved in the regulation of GA degradation. Overexpression of the wild EI gene in mutant P502 caused a dwarf phenotype with a shortened internode. The difference of EI expression levels was not significant in the P502 and wild-type, but the expression levels of GA biosynthetic genes including CPS, KO, KAO, GA20ox1, GA20ox2, GA20ox4, GA3ox1, GA2ox1, GA2ox2, GA2ox4, and GA2ox5, were upregulated in mutant P502. Our results may provide a better understanding of the genetics underlying the internode elongation and valuable information to improve plant architecture of the tomato.

Normal and Abortive Buds Transcriptomic Profiling of Broccoli ogu Cytoplasmic Male Sterile Line and Its Maintainer.

Bud abortion is the main factor affecting hybrid seeds' yield during broccoli cross breeding when using ogura cytoplasmic male sterile (ogu CMS) lines. However, the genes associated with bud abortion are poorly understood. We applied RNA sequencing to analyze the transcriptomes of normal and abortive buds of broccoli maintainer and ogu CMS lines. Functional analysis showed that among the 54,753 annotated unigenes obtained, 74 and 21 differentially expressed genes in common were upregulated and downregulated in ogu CMS abortive buds compared with ogu CMS normal buds, maintainer normal, and abortive buds, respectively. Nineteen of the common differentially expressed genes were enriched by GO terms associated with glycosyl hydrolases, reactive oxygen species scavenging, inhibitor, and protein degradation. Ethylene-responsive transcription factor 115 and transcriptional factor basic helix-loop-helix 137 were significantly upregulated; transcription factors DUO1 and PosF21/RF2a/BZIP34 were downregulated in ogu CMS abortive buds compared with the other groups. Genes related to polygalacturonase metabolism, glycosyl hydrolases, oxidation reduction process, phenylalanine metabolism, and phenylpropanoid biosynthesis were significantly changed in ogu CMS abortive buds. Our results increase our understanding of bud abortion, provide a valuable resource for further functional characterization of ogu CMS during bud abortion, and will aid in future cross breeding of Brassica crops.

QTL-seq for rapid identification of candidate genes for flowering time in broccoli × cabbage.

KEY MESSAGE: A major QTL controlling early flowering in broccoli × cabbage was identified by marker analysis and next-generation sequencing, corresponding to GRF6 gene conditioning flowering time in Arabidopsis. Flowering is an important agronomic trait for hybrid production in broccoli and cabbage, but the genetic mechanism underlying this process is unknown. In this study, segregation analysis with BC1P1, BC1P2, F2, and F2:3 populations derived from a cross between two inbred lines "195" (late-flowering) and "93219" (early flowering) suggested that flowering time is a quantitative trait. Next, employing a next-generation sequencing-based whole-genome QTL-seq strategy, we identified a major genomic region harboring a robust flowering time QTL using an F2 mapping population, designated Ef2.1 on cabbage chromosome 2 for early flowering. Ef2.1 was further validated by indel (insertion or deletion) marker-based classical QTL mapping, explaining 51.5% (LOD = 37.67) and 54.0% (LOD = 40.5) of the phenotypic variation in F2 and F2:3 populations, respectively. Combined QTL-seq and classical QTL analysis narrowed down Ef1.1 to a 228-kb genomic region containing 29 genes. A cabbage gene, Bol024659, was identified in this region, which is a homolog of GRF6, a major gene regulating flowering in Arabidopsis, and was designated BolGRF6. qRT-PCR study of the expression level of BolGRF6 revealed significantly higher expression in the early flowering genotypes. Taken together, our results provide support for BolGRF6 as a possible candidate gene for early flowering in the broccoli line 93219. The identified candidate genomic regions and genes may be useful for molecular breeding to improve broccoli and cabbage flowering times.

Detection of the Diversity of Cytoplasmic Male Sterility Sources in Broccoli (Brassica Oleracea var. Italica) Using Mitochondrial Markers.

Broccoli (Brassica oleracea var. italica) is an important commercial vegetable crop. As part of an efficient pollination system, cytoplasmic male sterility (CMS) has been widely used for broccoli hybrid production. Identifying the original sources of CMS in broccoli accessions has become an important part of broccoli breeding. In this study, the diversity of the CMS sources of 39 broccoli accessions, including 19 CMS lines and 20 hybrids, were analyzed using mitochondrial markers. All CMS accessions contained the ogu orf138-related DNA fragment and the key genes of nap CMS, pol CMS, and tour CMS were not detected. The 39 CMS accessions were divided into five groups using six orf138-related and two simple sequence repeat markers. We observed that ogu CMS R3 constituted 79.49% of the CMS sources. CMS6 and CMS26 were differentiated from the other accessions using a specific primer. CMS32 was distinguished from the other accessions based on a 78-nucleotide deletion at the same locus as the orf138-related sequence. When the coefficient was about 0.90, five CMS accessions (13CMS6, 13CMS23, 13CMS24, 13CMS37, and 13CMS39) exhibiting abnormal floral organs with poor seed setting were grouped together. The polymerase chain reaction amplification profiles for these five accessions differed from those of the other accessions. We identified eight useful molecular markers that can be used to detect CMS types during broccoli breeding. Our data also provide important information relevant to future studies on the possible origins and molecular mechanisms of CMS in broccoli.

Organelle Simple Sequence Repeat Markers Help to Distinguish Carpelloid Stamen and Normal Cytoplasmic Male Sterile Sources in Broccoli.

We previously discovered carpelloid stamens when breeding cytoplasmic male sterile lines in broccoli (Brassica oleracea var. italica). In this study, hybrids and multiple backcrosses were produced from different cytoplasmic male sterile carpelloid stamen sources and maintainer lines. Carpelloid stamens caused dysplasia of the flower structure and led to hooked or coiled siliques with poor seed setting, which were inherited in a maternal fashion. Using four distinct carpelloid stamens and twelve distinct normal stamens from cytoplasmic male sterile sources and one maintainer, we used 21 mitochondrial simple sequence repeat (mtSSR) primers and 32 chloroplast SSR primers to identify a mitochondrial marker, mtSSR2, that can differentiate between the cytoplasm of carpelloid and normal stamens. Thereafter, mtSSR2 was used to identify another 34 broccoli accessions, with an accuracy rate of 100%. Analysis of the polymorphic sequences revealed that the mtSSR2 open reading frame of carpelloid stamen sterile sources had a deletion of 51 bases (encoding 18 amino acids) compared with normal stamen materials. The open reading frame is located in the coding region of orf125 and orf108 of the mitochondrial genomes in Brassica crops and had the highest similarity with Raphanus sativus and Brassica carinata. The current study has not only identified a useful molecular marker to detect the cytoplasm of carpelloid stamens during broccoli breeding, but it also provides evidence that the mitochondrial genome is maternally inherited and provides a basis for studying the effect of the cytoplasm on flower organ development in plants.