Association mapping reveals the genetic architecture of tomato response to water deficit: focus on major fruit quality traits

Mitigating Water Stress in Plants with Beneficial Bacteria: Effects on Growth and Rhizosphere Bacterial Communities

Climate change has reshaped global weather patterns and intensified extreme events, with drought and soil salinity negatively impacting the yield and quality of crop production. To mitigate the detrimental effects of drought stress, the introduction of beneficial plant growth-promoting rhizobacteria (PGPR) has proven to be a promising approach. In this study, we evaluated a synthetic microbial community (SynCom) comprising bacterial strains belonging to the species Bacillus velezensis, Pseudomonas simiae, P. salmasensis, Glutamicibacter halophytocola, and Leclercia sp., which have been demonstrated to promote tomato growth both individually and collectively. The SynCom and most of its individual bacterial strains were shown to mitigate the detrimental effects of polyethylene glycol (PEG)-induced drought stress in vitro in Arabidopsis thaliana seedlings, either by reducing alterations in xylem elements or promoting the formation of new xylem strands. In a greenhouse trial, soil drenching with the SynCom and two individual strains, B. velezensis PSE31B and P. salmasensis POE54, improved the water stress response in soilless-grown tomato plants under a 40% reduced irrigation regime. Additionally, bacterial treatments positively influenced the diversity of rhizosphere bacterial communities, with distinct changes in bacterial composition, which suggest a treatment-specific interplay between the introduced strains and the native microbiome. These findings highlight the potential of microbial consortia and individual PGPR strains as sustainable tools to improve plant resilience to abiotic stresses.

Characterisation of a major QTL for sodium accumulation in tomato grown in high salinity

Soil salinity is a serious concern for tomato culture, affecting both yield and quality parameters. Although some genes involved in tomato salt tolerance have been identified, their genetic diversity has been rarely studied. In the present study, we assessed salt tolerance-related traits at juvenile and adult stages in a large core collection and identified salt tolerance quantitative trait loci (QTLs) by genome-wide association study (GWAS). The results suggested that a major QTL is involved in leaf sodium accumulation at both physiological stages. We were able to identify the underlying candidate gene, coding for a well-known sodium transporter, called SlHKT1.2. We showed that an eQTL for the expression of this gene in roots colocalized with the above ground sodium content QTL. A polymorphism putatively responsible for its variation was identified in the gene promoter. Finally, to extend the applicability of these results, we carried out the same analysis on a test-cross panel composed of the core collection crossed with a distant line. The results indicated that the identified QTL retained its functional impact even in a hybrid genetic context: this paves the way for its use in breeding programs aimed at improving salinity tolerance in tomato cultivars.

Integration of QTL and transcriptome approaches for the identification of genes involved in tomato response to nitrogen deficiency

Optimizing plant nitrogen (N) usage and inhibiting N leaching loss in the soil-crop system is crucial to maintaining crop yield and reducing environmental pollution. This study aimed at identifying quantitative trait loci (QTLs) and differentially expressed genes (DEGs) between two N treatments in order to list candidate genes related to nitrogen-related contrasting traits in tomato varieties. We characterized a genetic diversity core-collection (CC) and a multi-parental advanced generation intercross (MAGIC) tomato population grown in a greenhouse under two nitrogen levels and assessed several N-related traits and mapped QTLs. Transcriptome response under the two N conditions was also investigated through RNA sequencing of fruit and leaves in four parents of the MAGIC population. Significant differences in response to N input reduction were observed at the phenotypic level for biomass and N-related traits. Twenty-seven QTLs were detected for three target traits (leaf N content, leaf nitrogen balance index, and petiole NO3- content), 10 and six in the low and high N condition, respectively, while 19 QTLs were identified for plasticity traits. At the transcriptome level, 4752 and 2405 DEGs were detected between the two N conditions in leaves and fruits, respectively, among which 3628 (50.6%) in leaves and 1717 (71.4%) in fruit were genotype specific. When considering all the genotypes, 1677 DEGs were shared between organs or tissues. Finally, we integrated DEG and QTL analyses to identify the most promising candidate genes. The results highlighted a complex genetic architecture of N homeostasis in tomato and novel putative genes useful for breeding tomato varieties requiring less N input.

Genome-wide identification and functional analysis of the SiCIN gene family in foxtail millet (Setaria italica L.)

Cell wall invertase (CIN) is a vital member of plant invertase (INV) and plays a key role in the breakdown of sucrose. This enzyme facilitates the hydrolysis of sucrose into glucose and fructose, which is crucial for various aspects of plant growth and development. However, the function of CIN genes in foxtail millet (Setaria italica) is less studied. In this research, we used the blast-p of NCBI and TBtools for bidirectional comparison, and a total of 13 CIN genes (named SiCINs) were identified from foxtail millet by using Arabidopsis and rice CIN sequences as reference sequences. The phylogenetic tree analysis revealed that the CIN genes can be categorized into three subfamilies: group 1, group 2, and group 3. Furthermore, upon conducting chromosomal localization analysis, it was observed that the 13 SiCINs were distributed unevenly across five chromosomes. Cis-acting elements of SiCIN genes can be classified into three categories: plant growth and development, stress response, and hormone response. The largest number of cis-acting elements were those related to light response (G-box) and the cis-acting elements related to seed-specific regulation (RY-element). qRT-PCR analysis further confirmed that the expression of SiCIN7 and SiCIN8 in the grain was higher than that in any other tissues. The overexpression of SiCIN7 in Arabidopsis improved the grain size and thousand-grain weight, suggesting that SiCIN7 could positively regulate grain development. Our findings will help to further understand the grain-filling mechanism of SiCIN and elucidate the biological mechanism underlying the grain development of SiCIN.

Genomic and phenomic predictions help capture low-effect alleles promoting seed germination in oilseed rape in addition to QTL analyses

KEY MESSAGE

Phenomic prediction implemented on a large diversity set can efficiently predict seed germination, capture low-effect favorable alleles that are not revealed by GWAS and identify promising genetic resources. Oilseed rape faces many challenges, especially at the beginning of its developmental cycle. Achieving rapid and uniform seed germination could help to ensure a successful establishment and therefore enabling the crop to compete with weeds and tolerate stresses during the earliest developmental stages. The polygenic nature of seed germination was highlighted in several studies, and more knowledge is needed about low- to moderate-effect underlying loci in order to enhance seed germination effectively by improving the genetic background and incorporating favorable alleles. A total of 17 QTL were detected for seed germination-related traits, for which the favorable alleles often corresponded to the most frequent alleles in the panel. Genomic and phenomic predictions methods provided moderate-to-high predictive abilities, demonstrating the ability to capture small additive and non-additive effects for seed germination. This study also showed that phenomic prediction estimated phenotypic values closer to phenotypic values than GEBV. Finally, as the predictive ability of phenomic prediction was less influenced by the genetic structure of the panel, it is worth using this prediction method to characterize genetic resources, particularly with a view to design prebreeding populations.

Two gene clusters and their positive regulator SlMYB13 that have undergone domestication-associated negative selection control phenolamide accumulation and drought tolerance in tomato

Among plant metabolites, phenolamides, which are conjugates of hydroxycinnamic acid derivatives and polyamines, play important roles in plant adaptation to abiotic and biotic stresses. However, the molecular mechanisms underlying phenolamide metabolism and regulation as well as the effects of domestication and breeding on phenolamide diversity in tomato remain largely unclear. In this study, we performed a metabolite-based genome-wide association study and identified two biosynthetic gene clusters (BGC7 and BGC11) containing 12 genes involved in phenolamide metabolism, including four biosynthesis genes (two 4CL genes, one C3H gene, and one CPA gene), seven decoration genes (five AT genes and two UGT genes), and one transport protein gene (DTX29). Using gene co-expression network analysis we further discovered that SlMYB13 positively regulates the expression of two gene clusters, thereby promoting phenolamide accumulation. Genetic and physiological analyses showed that BGC7, BGC11 and SlMYB13 enhance drought tolerance by enhancing scavenging of reactive oxygen species and increasing abscisic acid content in tomato. Natural variation analysis suggested that BGC7, BGC11 and SlMYB13 were negatively selected during tomato domestication and improvement, leading to reduced phenolamide content and drought tolerance of cultivated tomato. Collectively, our study discovers a key mechanism of phenolamide biosynthesis and regulation in tomato and reveals that crop domestication and improvement shapes metabolic diversity to affect plant environmental adaptation.

Plasticity of maternal environment-dependent expression-QTLs of tomato seeds

Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.

Salt-tolerant alternative crops as sources of quality food to mitigate the negative impact of salinity on agricultural production

An increase of abiotic stress tolerance and nutritive value of foods is currently a priority because of climate change and rising world population. Among abiotic stresses, salt stress is one of the main problems in agriculture. Mounting urbanization and industrialization, and increasing global food demand, are pressing farmers to make use of marginal lands affected by salinity and low-quality saline water. In that situation, one of the most promising approaches is searching for new sources of genetic variation like salt-tolerant alternative crops or underexploited crops. They are generally less efficient than cultivated crops in optimal conditions due to lower yield but represent an alternative in stressful growth conditions. In this review, we summarize the advances achieved in research on underexploited species differing in their genetic nature. First, we highlight advances in research on salt tolerance of traditional varieties of tomato or landraces; varieties selected and developed by smallholder farmers for adaptation to their local environments showing specific attractive fruit quality traits. We remark advances attained in screening a collection of tomato traditional varieties gathered in Spanish Southeast, a very productive region which environment is extremely stressing. Second, we explore the opportunities of exploiting the natural variation of halophytes, in particular quinoa and amaranth. The adaptation of both species in stressful growth conditions is becoming an increasingly important issue, especially for their cultivation in arid and semiarid areas prone to be affected by salinity. Here we present a project developed in Spanish Southeast, where quinoa and amaranth varieties are being adapted for their culture under abiotic stress targeting high quality grain.

Effects of limited water supply on metabolite composition in tomato fruits (Solanum lycopersicum L.) in two soils with different nutrient conditions

Effect of water supply to metabolites in tomato fruit was compared in two soils with different nutrient conditions, i.e., either limited or excess. Two types of soil nutrient condition, type A: nutrient-limited and type B: nutrient-excess, were prepared as follows; type A is a low nutrient-containing soil without a replenishment of starved nitrogen and phosphorous, type B is a high nutrient-containing soil exceeding the recommended fertilization. Soil water was adjusted either at -30 kPa (sufficient) or -80 kPa (limited). For harvested tomato fruits, we examined primary and secondary metabolites using non-targeted mass spectrometry based metabolomics. The fruit production and leaf SPAD were greatly dependent on soil nutrient levels, by contrast, the level of lycopene remained unchanged by different levels of water and nutrient supply. The perturbation of metabolites by water supply was clear in the nutrient-excess soil. In particular, limited water supply strongly decreased primary metabolites including sugars and amino acids. We demonstrated that water stress differently shifted primary metabolites of tomato fruits in two soils with different nutrient conditions via non-targeted mass spectrometry-based metabolomics. In conclusion, we suggest that the limited water supply in soils with surplus nutrient is not a recommendable way for tomato 'cv. Super Dotaerang' production if fruit nutritional quality such as sugars and amino acids is in the consideration, although there was no disadvantage in fruit yield.

Drought tolerance improvement in Solanum lycopersicum: an insight into "OMICS" approaches and genome editing

Solanum lycopersicum (tomato) is an internationally acclaimed vegetable crop that is grown worldwide. However, drought stress is one of the most critical challenges for tomato production, and it is a crucial task for agricultural biotechnology to produce drought-resistant cultivars. Although breeders have done a lot of work on the tomato to boost quality and quantity of production and enhance resistance to biotic and abiotic stresses, conventional tomato breeding approaches have been limited to improving drought tolerance because of the intricacy of drought traits. Many efforts have been made to better understand the mechanisms involved in adaptation and tolerance to drought stress in tomatoes throughout the years. "Omics" techniques, such as genomics, transcriptomics, proteomics, and metabolomics in combination with modern sequencing technologies, have tremendously aided the discovery of drought-responsive genes. In addition, the availability of biotechnological tools, such as plant transformation and the recently developed genome editing system for tomatoes, has opened up wider opportunities for validating the function of drought-responsive genes and the generation of drought-tolerant varieties. This review highlighted the recent progresses for tomatoes improvement against drought stress through "omics" and "multi-omics" technologies including genetic engineering. We have also discussed the roles of non-coding RNAs and genome editing techniques for drought stress tolerance improvement in tomatoes.

The INRAE Centre for Vegetable Germplasm: Geographically and Phenotypically Diverse Collections and Their Use in Genetics and Plant Breeding

The French National Research Institute for Agriculture, Food and the Environment (INRAE) conserves and distributes five vegetable collections as seeds: the aubergine* (in this article the word aubergine refers to eggplant), pepper, tomato, melon and lettuce collections, together with their wild or cultivated relatives, are conserved in Avignon, France. Accessions from the collections have geographically diverse origins, are generally well-described and fixed for traits of agronomic or scientific interest and have available passport data. In addition to currently conserving over 10,000 accessions (between 900 and 3000 accessions per crop), the centre maintains scientific collections such as core collections and bi- or multi-parental populations, which have also been genotyped with SNP markers. Each collection has its own merits and highlights, which are discussed in this review: the aubergine collection is a rich source of crop wild relatives of Solanum; the pepper, melon and lettuce collections have been screened for resistance to plant pathogens, including viruses, fungi, oomycetes and insects; and the tomato collection has been at the heart of genome-wide association studies for fruit quality traits and environmental stress tolerance.

The Bacterial Microbiome of the Tomato Fruit Is Highly Dependent on the Cultivation Approach and Correlates With Flavor Chemistry

The modes of interactions between plants and plant-associated microbiota are manifold, and secondary metabolites often play a central role in plant-microbe interactions. Abiotic and biotic (including both plant pathogens and endophytes) stress can affect the composition and concentration of secondary plant metabolites, and thus have an influence on chemical compounds that make up for the taste and aroma of fruit. While the role of microbiota in growth and health of plants is widely acknowledged, relatively little is known about the possible effect of microorganisms on the quality of fruit of plants they are colonizing. In this work, tomato (Solanum lycopersicum L.) plants of five different cultivars were grown in soil and in hydroponics to investigate the impact of the cultivation method on the flavor of fruit, and to assess whether variations in their chemical composition are attributable to shifts in bacterial microbiota. Ripe fruit were harvested and used for bacterial community analysis and for the analysis of tomato volatiles, sugars and acids, all contributing to flavor. Fruit grown in soil showed significantly higher sugar content, whereas tomatoes from plants under hydroponic conditions had significantly higher levels of organic acids. In contrast, aroma profiles of fruit were shaped by the tomato cultivars, rather than the cultivation method. In terms of bacterial communities, the cultivation method significantly defined the community composition in all cultivars, with the bacterial communities in hydroponic tomatoes being more variable that those in tomatoes grown in soil. Bacterial indicator species in soil-grown tomatoes correlated with higher concentrations of volatiles described to be perceived as "green" or "pungent." A soil-grown specific reproducibly occurring ASV (amplicon sequence variants) classified as Bacillus detected solely in "Solarino" tomatoes, which were the sweetest among all cultivars, correlated with the amount of aroma-relevant volatiles as well as of fructose and glucose in the fruit. In contrast, indicator bacterial species in hydroponic-derived tomatoes correlated with aroma compounds with "sweet" and "floral" notes and showed negative correlations with glucose concentrations in fruit. Overall, our results point toward a microbiota-related accumulation of flavor and aroma compounds in tomato fruit, which is strongly dependent on the cultivation substrate and approach.

Genome wide association mapping for agronomic, fruit quality, and root architectural traits in tomato under organic farming conditions

BACKGROUND

Opportunity and challenges of the agriculture scenario of the next decades will face increasing demand for secure food through approaches able to minimize the input to cultivations. Large panels of tomato varieties represent a valuable resource of traits of interest under sustainable cultivation systems and for genome-wide association studies (GWAS). For mapping loci controlling the variation of agronomic, fruit quality, and root architecture traits, we used a heterogeneous set of 244 traditional and improved tomato accessions grown under organic field trials. Here we report comprehensive phenotyping and GWAS using over 37,300 SNPs obtained through double digest restriction-site associated DNA (dd-RADseq).

RESULTS

A wide range of phenotypic diversity was observed in the studied collection, with highly significant differences encountered for most traits. A variable level of heritability was observed with values up to 69% for morphological traits while, among agronomic ones, fruit weight showed values above 80%. Genotype by environment analysis highlighted the strongest genotypic effect for aboveground traits compared to root architecture, suggesting that the hypogeal part of tomato plants has been a minor objective for breeding activities. GWAS was performed by a compressed mixed linear model leading to 59 significantly associated loci, allowing the identification of novel genes related to flower and fruit characteristics. Most genomic associations fell into the region surrounding SUN, OVATE, and MYB gene families. Six flower and fruit traits were associated with a single member of the SUN family (SLSUN31) on chromosome 11, in a region involved in the increase of fruit weight, locules number, and fruit fasciation. Furthermore, additional candidate genes for soluble solids content, fruit colour and shape were found near previously reported chromosomal regions, indicating the presence of synergic and multiple linked genes underlying the variation of these traits.

CONCLUSIONS

Results of this study give new hints on the genetic basis of traits in underexplored germplasm grown under organic conditions, providing a framework for the development of markers linked to candidate genes of interest to be used in genomics-assisted breeding in tomato, in particular under low-input and organic cultivation conditions.

Breeding Tomato Hybrids for Flavour: Comparison of GWAS Results Obtained on Lines and F1 Hybrids

Tomato flavour is an important goal for breeders. Volatile organic compounds (VOCs) are major determinants of tomato flavour. Although most tomato varieties for fresh market are F1 hybrids, most studies on the genetic control of flavour-related traits are performed on lines. We quantified 46 VOCs in a panel of 121 small fruited lines and in a test cross panel of 165 hybrids (the previous panel plus 44 elite cherry tomato lines crossed with a common line). High and consistent heritabilities were assessed for most VOCs in the two panels, and 65% of VOC contents were strongly correlated between lines and hybrids. Additivity was observed for most VOCs. We performed genome wide association studies (GWAS) on the two panels separately, along with a third GWAS on the test cross subset carrying only F1 hybrids corresponding to the line panel. We identified 205, 183 and 138 associations, respectively. We identified numerous overlapping associations for VOCs belonging to the same metabolic pathway within each panel; we focused on seven chromosome regions with clusters of associations simultaneously involved in several key VOCs for tomato aroma. The study highlighted the benefit of testcross panels to create tasty F1 hybrid varieties.

Genetic diversity of tomato response to heat stress at the QTL and transcriptome levels

Tomato is a widely cultivated crop, which can grow in many environments. However, temperature above 30°C impairs its reproduction, subsequently impacting fruit yield. We assessed the impact of high-temperature stress (HS) in two tomato experimental populations, a multi-parental advanced generation intercross (MAGIC) population and a core-collection (CC) of small-fruited tomato accessions. Both populations were evaluated for 11 traits related to yield components, phenology and fruit quality in optimal and HS conditions. HS significantly impacted all traits in both populations, but a few genotypes with stable yield under HS were identified. A plasticity index was computed for each individual to measure the extent of the heat impact for each trait. Quantitative trait loci (QTL) were detected in control and HS conditions as well as for plasticity index. Linkage and genome-wide association analyses in the MAGIC and CC populations identified a total of 98 and 166 QTLs, respectively. Taking the two populations together, 69 plasticity QTLs (pQTLs) were involved in tomato heat response for 11 traits. The transcriptome changes in the ovary of six genotypes with contrasted responses to HS were studied, and 837 genes differentially expressed according to the conditions were detected. Combined with previous transcriptome studies, these results were used to propose candidate genes for HS response QTLs.

A transcriptomic, metabolomic and cellular approach to the physiological adaptation of tomato fruit to high temperature

High temperatures can negatively influence plant growth and development. Besides yield, the effects of heat stress on fruit quality traits remain poorly characterised. In tomato, insights into how fruits regulate cellular metabolism in response to heat stress could contribute to the development of heat-tolerant varieties, without detrimental effects on quality. In the present study, the changes occurring in wild type tomato fruits after exposure to transient heat stress have been elucidated at the transcriptome, cellular and metabolite level. An impact on fruit quality was evident as nutritional attributes changed in response to heat stress. Fruit carotenogenesis was affected, predominantly at the stage of phytoene formation, although altered desaturation/isomerisation arose during the transient exposure to high temperatures. Plastidial isoprenoid compounds showed subtle alterations in their distribution within chromoplast sub-compartments. Metabolite profiling suggests limited effects on primary/intermediary metabolism but lipid remodelling was evident. The heat-induced molecular signatures included the accumulation of sucrose and triacylglycerols, and a decrease in the degree of membrane lipid unsaturation, which influenced the volatile profile. Collectively, these data provide valuable insights into the underlying biochemical and molecular adaptation of fruit to heat stress and will impact on our ability to develop future climate resilient tomato varieties.

Validated MAGIC and GWAS population mapping reveals the link between vitamin E content and natural variation in chorismate metabolism in tomato

Tocochromanols constitute the different forms of vitamin E (VTE), essential components of the human diet, and display a high membrane protectant activity. By combining interval mapping and genome-wide association studies (GWAS), we unveiled the genetic determinants of tocochromanol accumulation in tomato (Solanum lycopersicum) fruits. To enhance the nutritional value of this highly consumed vegetable, we dissected the natural intraspecific variability of tocochromanols in tomato fruits and genetically engineered their biosynthetic pathway. These analyses allowed the identification of a total of 25 quantitative trait loci interspersed across the genome pinpointing the chorismate-tyrosine pathway as a regulatory hub controlling the supply of the aromatic head group for tocochromanol biosynthesis. To validate the link between the chorismate-tyrosine pathway and VTE, we engineered tomato plants to bypass the pathway at the arogenate branch point. Transgenic tomatoes showed moderate increments in tocopherols (up to approximately 20%) and a massive accumulation of tocotrienols (up to approximately 3400%). Gene expression analyses of these plants reveal a trade-off between VTE and natural variation in chorismate metabolism explained by transcriptional reprogramming of specific structural genes of the pathway. By restoring the accumulation of alpha-tocotrienols (α-t3) in fruits, the plants produced here are of high pharmacological and nutritional interest.

Traditional Tomato Varieties Improve Fruit Quality Without Affecting Fruit Yield Under Moderate Salt Stress

Identification of tomato varieties able to exhibit higher accumulation of primary and secondary metabolites in their fruits is currently a main objective in tomato breeding. One tool to improve fruit quality is to cultivate the plants under salt stress, although improvement of fruit quality is generally accompanied by productivity losses. However, it is very interesting to implement strategies aiming at enhancing fruit quality of tomato by means of growing plants in moderate salt stress that allows for a sustainable fruit yield. The traditional tomato varieties adapted to the Mediterranean environmental constraints may be very attractive plant materials to achieve this goal, given the wide range of fruit quality traits because of their genetic diversity. Here, agronomic responses and fruit quality traits, including primary and secondary metabolites, were analyzed in fruits of two Mediterranean traditional tomato varieties named "Tomate Pimiento" ("TP") and "Muchamiel Aperado" ("MA") because of the pepper and pear shape of their fruits, using as reference the commercial cultivar "Moneymaker" ("MM"). Plants were grown without salt (control) and with moderate salt stress (50 mM NaCl), which did not affect fruit yield in any variety. "TP" is of great interest because of its high soluble solids content (SSC) in control, which is even higher in salt, whereas "MA" is very attractive because of its high Brix yield index (SSC × fruit yield), used as overall fruit quality measure. Similitude between both traditional varieties were found for primary metabolism, as they significantly increased sucrose contents compared with "MM" in red ripe fruits from plants in control and, especially, salt stress conditions. The most remarkable difference was the high constitutive levels of total amino acids in "TP" fruits, including the three major free amino acids found in tomato fruit, GABA, glutamate, and glutamine, which even increased under salinity. Regarding secondary metabolites, the most interesting change induced by salinity was the increase in α-tocopherol found in red ripe fruits of both "TP" and "MA." These results reveal the interest of traditional varieties as sources of genetic variation in breeding because of their improvement of tomato fruit quality without production losses under moderate salt stress.

Harnessing High-throughput Phenotyping and Genotyping for Enhanced Drought Tolerance in Crop Plants

Development of drought-tolerant cultivars is one of the challenging tasks for the plant breeders due to its complex inheritance and polygenic regulation. Evaluating genetic material for drought tolerance is a complex process due to its spatiotemporal interactions with environmental factors. The conventional breeding approaches are costly, lengthy, and inefficient to achieve the expected gain in drought tolerance. In this regard, genomics-assisted breeding (GAB) offers promise to develop cultivars with improved drought tolerance in a more efficient, quicker, and cost-effective manner. The success of GAB depends upon the precision in marker-trait association and estimation of genomic estimated breeding values (GEBVs), which mostly depends on coverage and precision of genotyping and phenotyping. A wide gap between the discovery and practical use of quantitative trait loci (QTL) for crop improvement has been observed for many important agronomical traits. Such a limitation could be due to the low accuracy in QTL detection, mainly resulting from low marker density and manually collected phenotypes of complex agronomic traits. Increasing marker density using the high-throughput genotyping (HTG), and accurate and precise phenotyping using high-throughput digital phenotyping (HTP) platforms can improve the precision and power of QTL detection. Therefore, both HTG and HTP can enhance the practical utility of GAB along with a faster characterization of germplasm and breeding material. In the present review, we discussed how the recent innovations in HTG and HTP would assist in the breeding of improved drought-tolerant varieties. We have also discussed strategies, tools, and analytical advances made on the HTG and HTP along with their pros and cons.

Analysis of wild tomato introgression lines elucidates the genetic basis of transcriptome and metabolome variation underlying fruit traits and pathogen response

Wild tomato species represent a rich gene pool for numerous desirable traits lost during domestication. Here, we exploited an introgression population representing wild desert-adapted species and a domesticated cultivar to establish the genetic basis of gene expression and chemical variation accompanying the transfer of wild-species-associated fruit traits. Transcriptome and metabolome analysis of 580 lines coupled to pathogen sensitivity assays resulted in the identification of genomic loci associated with levels of hundreds of transcripts and metabolites. These associations occurred in hotspots representing coordinated perturbation of metabolic pathways and ripening-related processes. Here, we identify components of the Solanum alkaloid pathway, as well as genes and metabolites involved in pathogen defense and linking fungal resistance with changes in the fruit ripening regulatory network. Our results outline a framework for understanding metabolism and pathogen resistance during tomato fruit ripening and provide insights into key fruit quality traits.

Genetic basis of phenotypic plasticity and genotype × environment interactions in a multi-parental tomato population

Deciphering the genetic basis of phenotypic plasticity and genotype × environment interactions (G×E) is of primary importance for plant breeding in the context of global climate change. Tomato (Solanum lycopersicum) is a widely cultivated crop that can grow in different geographical habitats and that displays a great capacity for expressing phenotypic plasticity. We used a multi-parental advanced generation intercross (MAGIC) tomato population to explore G×E and plasticity for multiple traits measured in a multi-environment trial (MET) comprising optimal cultural conditions together with water deficit, salinity, and heat stress over 12 environments. Substantial G×E was observed for all the traits measured. Different plasticity parameters were estimated by employing Finlay-Wilkinson and factorial regression models and these were used together with genotypic means for quantitative trait loci (QTL) mapping analyses. In addition, mixed linear models were also used to investigate the presence of QTL × environment interactions. The results highlighted a complex genetic architecture of tomato plasticity and G×E. Candidate genes that might be involved in the occurrence of G×E are proposed, paving the way for functional characterization of stress response genes in tomato and for breeding climate-adapted cultivars.

Integration of QTL, Transcriptome and Polymorphism Studies Reveals Candidate Genes for Water Stress Response in Tomato

Water deficit (WD) leads to significant phenotypic changes in crops resulting from complex stress regulation mechanisms involving responses at the physiological, biochemical and molecular levels. Tomato growth and fruit quality have been shown to be significantly affected by WD stress. Understanding the molecular mechanism underlying response to WD is crucial to develop tomato cultivars with relatively high performance under low watering conditions. Transcriptome response to WD was investigated through the RNA sequencing of fruit and leaves in eight accessions grown under two irrigation conditions, in order to get insight into the complex genetic regulation of WD response in tomato. Significant differences in genotype WD response were first observed at the phenotypic level for fruit composition and plant development traits. At the transcriptome level, a total of 14,065 differentially expressed genes (DEGs) in response to WD were detected, among which 7393 (53%) and 11,059 (79%) were genotype- and organ-specific, respectively. Water deficit induced transcriptome variations much stronger in leaves than in fruit. A significant effect of the genetic background on expression variation was observed compared to the WD effect, along with the presence of a set of genes showing a significant genotype × watering regime interaction. Integrating the DEGs with previously identified WD response quantitative trait loci (QTLs) mapped in a multi-parental population derived from the crossing of the eight genotypes narrowed the candidate gene lists to within the confidence intervals surrounding the QTLs. The results present valuable resources for further study to decipher the genetic determinants of tomato response to WD.

Fruit ripening in Lycium barbarum and Lycium ruthenicum is associated with distinct gene expression patterns

Goji berries have been used as food and medicine for millennia. Due to their high morphological similarity, fruits of two distinct species belonging to the family Solanaceae, Lycium barbarum (LB) and Lycium chinense (Chinese boxthorn), are usually marketed together as goji berries, but nearly 90% of all commercially available goji berries belong to the former species. A third closely related species, a wild perennial thorny shrub native to north‐western China, Lycium ruthenicum (LR; known as Russian box thorn, and its fruit as black wolfberry), has become a popular choice for combating soil desertification and for alleviating soil salinity/alkalinity due to its high resistance to the harsh environment of saline deserts. Despite the phylogenetic closeness of LB and LR, their fruits are very different. To identify the genes involved in these distinct phenotypes, here we studied expression patterns of 22 transcriptional regulators that may be crucial drivers of these differences during five developmental stages. BAM1 may contribute to higher sugar content in LB. High expression of BFRUCT in ripe LR is likely to be an evolutionary adaptation to fruit ripening in an arid environment. Two arogenate dehydratase paralogues, CHS and LDOX, are probably crucial elements of the mechanism by which LR accumulates much higher levels of anthocyanin. DXS2 (carotenoid accumulation in LB) and CCD4 (carotenoid degradation in ripe LR fruit) may be crucial drivers behind the much higher content of carotenoids in LB. EIL3 and ERF5 are two transcription factors that may contribute to the higher abiotic stress resilience of LR. GATA22‐like appears to have more important roles in growth than ripening in LB fruit and vice versa in LR. HAT5‐like exhibited opposite temporal patterns in two fruits: high in the 1st stage in LB and high in the 5th stage in LR. PED1 was expressed at a much lower level in LR. Finally, we hypothesise that the poorly functionally characterised SCL32 gene may play a part in the increased resistance to environmental stress of LR. We suggest that BAM1, BFRUCT, EIL3, ERF5, ADT paralogues (for functional redundancy), PED1, GATA22‐like, HAT5‐like and SCL32 warrant further functional studies.

Mapping of quantitative trait loci for antioxidant molecules in tomato fruit: Carotenoids, vitamins C and E, glutathione and phenolic acids

The nutritional value of a crop lies not only in its protein, lipid, and sugar content but also involves compounds such as the antioxidants lycopene, β-carotene and vitamin C. In the present study, wild tomato Solanum pimpinellifolium LA 1589 was assessed for its potential to improve antioxidant content. This wild species was found to be a good source of alleles for increasing β-carotene, lycopene, vitamin C and vitamin E contents in cultivated tomato. Characterization of an LA 1589 interspecific inbred backcross line (IBL) mapping population revealed many individuals with transgressive segregation for the antioxidants confirming the usefulness of this wild species for breeding of these traits. Molecular markers were used to identify QTLs for the metabolites in the IBL population. In total, 64 QTLs were identified for the antioxidants and their locations were compared to the map positions of previously identified QTLs for confirmation. Four (57 %) of the carotenoid QTLs, four (36 %) of the vitamin QTLs, and 11 (25 %) of the phenolic acid QTLs were supported by previous studies. Furthermore, several potential candidate genes were identified for vitamins C and E and phenolic acids loci. These candidate genes might be used as markers in breeding programs to increase tomato's antioxidant content.

Responses of water accumulation and solute metabolism in tomato fruit to water scarcity and implications for main fruit quality variables

Fruit is important for human health, and applying deficit irrigation in fruit production is a strategy to regulate fruit quality and support environmental sustainability. Responses of different fruit quality variables to deficit irrigation have been widely documented, and much progress has been made in understanding the mechanisms of these responses. We review the effects of water shortage on fruit water accumulation considering water transport from the parent plant into the fruit determined by hydraulic properties of the pathway (including xylem water transport and transmembrane water transport regulated by aquaporins) and the driving force for water movement. We discuss water relations and solute metabolism that affect the main fruit quality variables (e.g. size, flavour, nutrition, and firmness) at the cellular level under water shortage. We also summarize the most recent advances in the understanding of responses of the main fruit quality variables to water shortage, considering the effects of variety, the severity of water deficit imposed, and the developmental stage of the fruit. We finally identify knowledge gaps and suggest avenues for future research. This review provides new insights into the stress physiology of fleshy fruit, which will be beneficial for the sustainable production of high-quality fruit under deficit irrigation.

Tomato Fruit Development and Metabolism

Tomato (Solanum lycopersicum L.) belongs to the Solanaceae family and is the second most important fruit or vegetable crop next to potato (Solanum tuberosum L.). It is cultivated for fresh fruit and processed products. Tomatoes contain many health-promoting compounds including vitamins, carotenoids, and phenolic compounds. In addition to its economic and nutritional importance, tomatoes have become the model for the study of fleshy fruit development. Tomato is a climacteric fruit and dramatic metabolic changes occur during its fruit development. In this review, we provide an overview of our current understanding of tomato fruit metabolism. We begin by detailing the genetic and hormonal control of fruit development and ripening, after which we document the primary metabolism of tomato fruits, with a special focus on sugar, organic acid, and amino acid metabolism. Links between primary and secondary metabolic pathways are further highlighted by the importance of pigments, flavonoids, and volatiles for tomato fruit quality. Finally, as tomato plants are sensitive to several abiotic stresses, we briefly summarize the effects of adverse environmental conditions on tomato fruit metabolism and quality.

Genome Wide Association Study of Karnal Bunt Resistance in a Wheat Germplasm Collection from Afghanistan

Karnal bunt disease of wheat, caused by the fungus Neovossia indica, is one of the most important challenges to the grain industry as it affects the grain quality and also restricts the international movement of infected grain. It is a seed-, soil- and airborne disease with limited effect of chemical control. Currently, this disease is contained through the deployment of host resistance but further improvement is limited as only a few genotypes have been found to carry partial resistance. To identify genomic regions responsible for resistance in a set of 339 wheat accessions, genome-wide association study (GWAS) was undertaken using the DArTSeq^® technology, in which 18 genomic regions for Karnal bunt resistance were identified, explaining 5–20% of the phenotypic variation. The identified quantitative trait loci (QTL) on chromosome 2BL showed consistently significant effects across all four experiments, whereas another QTL on 5BL was significant in three experiments. Additional QTLs were mapped on chromosomes 1DL, 2DL, 4AL, 5AS, 6BL, 6BS, 7BS and 7DL that have not been mapped previously, and on chromosomes 4B, 5AL, 5BL and 6BS, which have been reported in previous studies. Germplasm with less than 1% Karnal bunt infection have been identified and can be used for resistance breeding. The SNP markers linked to the genomic regions conferring resistance to Karnal bunt could be used to improve Karnal bunt resistance through marker-assisted selection.

Unravelling the Complex Genetics of Karnal Bunt (Tilletia indica) Resistance in Common Wheat (Triticum aestivum) by Genetic Linkage and Genome-Wide Association Analyses

Karnal bunt caused by Tilletia indica Mitra [syn. Neovossia indica (Mitra) Mundkur] is a significant biosecurity concern for wheat-exporting countries that are free of the disease. It is a seed-, soil-and air-borne disease with no effective chemical control measures. The current study used data from multi-year field experiments of two bi-parental populations and a genome-wide association (GWA) mapping panel to unravel the genetic basis for resistance in common wheat. Broad-sense heritability for Karnal bunt resistance in the populations varied from 0.52 in the WH542×HD29 population, to 0.61 in the WH542×W485 cross and 0.71 in a GWAS panel. Quantitative trait locus (QTL) analysis with seven years of phenotypic data identified a major locus on chromosome 3B (R² = 27.8%) and a minor locus on chromosome 1A (R² = 12.2%), in the WH542×HD29 population, with both parents contributing the high-value alleles. A major locus (R² = 27.8%) and seven minor loci (R² = 4.4–15.8%) were detected in the WH542×W485 population. GWA mapping validated QTL regions in the bi-parent populations, but also identified novel loci not previously associated with Karnal bunt resistance. Meta-QTL analysis aligned the results from this study with those reported in wheat over the last two decades. Two major clusters were detected, the first on chromosome 4B, which clustered with Qkb.ksu-4B, QKb.cimmyt-4BL, Qkb.cim-4BL, and the second on chromosome 3B, which clustered with Qkb.cnl-3B, QKb.cimmyt-3BS and Qkb.cim-3BS1. The results provide definitive chromosomal assignments for QTL/genes controlling Karnal bunt resistance in common wheat, and will be useful in pre-emptive breeding against the pathogen in wheat-producing areas that are free of the disease.

Genome-wide association multi-locus and multi-variate linear mixed models reveal two linked loci with major effects on partial resistance of apricot to bacterial canker

BACKGROUND

Diseases caused by Pseudomonas syringae (Ps) are recognized as the most damaging factors in fruit trees with a significant economic and sanitary impact on crops. Among them, bacterial canker of apricot is exceedingly difficult to control due to a lack of efficient prophylactic measures. Several sources of partial resistance have been identified among genetic resources but the underlying genetic pattern has not been elucidated thus far. In this study, we phenotyped bacterial canker susceptibility in an apricot core-collection of 73 accessions over 4 years by measuring canker and superficial browning lengths issued from artificial inoculations in the orchard. In order to investigate the genetic architecture of partial resistance, we performed a genome-wide association study using best linear unbiased predictors on genetic (G) and genetic x year (G × Y) interaction effects extracted from linear mixed models. Using a set of 63,236 single-nucleotide polymorphism markers genotyped in the germplasm over the whole genome, multi-locus and multi-variate mixed models aimed at mapping the resistance while controlling for relatedness between individuals.

RESULTS

We detected 11 significant associations over 7 candidate loci linked to disease resistance under the two most severe years. Colocalizations between G and G × Y terms indicated a modulation on allelic effect depending on environmental conditions. Among the candidate loci, two loci on chromosomes 5 and 6 had a high impact on both canker length and superficial browning, explaining 41 and 26% of the total phenotypic variance, respectively. We found unexpected long-range linkage disequilibrium (LD) between these two markers revealing an inter-chromosomal LD block linking the two underlying genes. This result supports the hypothesis of a co-adaptation effect due to selection through population demography. Candidate genes annotations suggest a functional pathway involving abscisic acid, a hormone mainly known for mediating abiotic stress responses but also reported as a potential factor in plant-pathogen interactions.

CONCLUSIONS

Our study contributed to the first detailed characterization of the genetic determinants of partial resistance to bacterial canker in a Rosaceae species. It provided tools for fruit tree breeding by identifying progenitors with favorable haplotypes and by providing major-effect markers for a marker-assisted selection strategy.

Trait discovery and editing in tomato

Tomato (Solanum lycopersicum), which is used for both processing and fresh markets, is a major crop species that is the top ranked vegetable produced over the world. Tomato is also a model species for research in genetics, fruit development and disease resistance. Genetic resources available in public repositories comprise the 12 wild related species and thousands of landraces, modern cultivars and mutants. In addition, high quality genome sequences are available for cultivated tomato and for several wild relatives, hundreds of accessions have been sequenced, and databases gathering sequence data together with genetic and phenotypic data are accessible to the tomato community. Major breeding goals are productivity, resistance to biotic and abiotic stresses, and fruit sensorial and nutritional quality. New traits, including resistance to various biotic and abiotic stresses and root architecture, are increasingly being studied. Several major mutations and quantitative trait loci (QTLs) underlying traits of interest in tomato have been uncovered to date and, thanks to new populations and advances in sequencing technologies, the pace of trait discovery has considerably accelerated. In recent years, clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing (GE) already proved its remarkable efficiency in tomato for engineering favorable alleles and for creating new genetic diversity by gene disruption, gene replacement, and precise base editing. Here, we provide insight into the major tomato traits and underlying causal genetic variations discovered so far and review the existing genetic resources and most recent strategies for trait discovery in tomato. Furthermore, we explore the opportunities offered by CRISPR/Cas9 and their exploitation for trait editing in tomato.

Allele-specific expression and genetic determinants of transcriptomic variations in response to mild water deficit in tomato

Characterizing the natural diversity of gene expression across environments is an important step in understanding how genotype-by-environment interactions shape phenotypes. Here, we analyzed the impact of water deficit onto gene expression levels in tomato at the genome-wide scale. We sequenced the transcriptome of growing leaves and fruit pericarps at cell expansion stage in a cherry and a large fruited accession and their F₁ hybrid grown under two watering regimes. Gene expression levels were steadily affected by the genotype and the watering regime. Whereas phenotypes showed mostly additive inheritance, ~80% of the genes displayed non-additive inheritance. By comparing allele-specific expression (ASE) in the F₁ hybrid to the allelic expression in both parental lines, respectively, 3005 genes in leaf and 2857 genes in fruit deviated from 1:1 ratio independently of the watering regime. Among these genes, ~55% were controlled by cis factors, ~25% by trans factors and ~20% by a combination of both types of factors. A total of 328 genes in leaf and 113 in fruit exhibited significant ASE-by-watering regime interaction, among which ~80% presented trans-by-watering regime interaction, suggesting a response to water deficit mediated through a majority of trans-acting loci in tomato. We cross-validated the expression levels of 274 transcripts in fruit and leaves of 124 recombinant inbred lines (RILs) and identified 163 expression quantitative trait loci (eQTLs) mostly confirming the divergences identified by ASE. Combining phenotypic and expression data, we observed a complex network of variation between genes encoding enzymes involved in the sugar metabolism.

Validation of QTLs for plant ideotype, earliness and growth habit traits in pigeonpea (Cajanus cajan Millsp.)

Pigeonpea productivity is greatly constrained by poor plant ideotype of existing Indian cultivars. Enhancing pigeonpea yield demands a renewed focus on restructuring the ideal plant type by using more efficient approaches like genomic tools. Therefore, the present study aims to identify and validate a set of QTLs/gene(s) presumably associated with various plant ideotype traits in pigeonpea. A total of 133 pigeonpea germplasms were evaluated along with four checks in the augmented design for various ideotype traits i.e. initiation of flowering (IF), days to 50% flowering (DFF), days to maturity (DM), plant height (PH), primary branches (PB), seeds per pod (SP) and pod length (PL). We observed significant genetic diversity in the germplasm lines for these traits. The genetic control of IF, DFF, DM and PH renders these traits suitable for detection of marker trait associations. By using residual maximum likelihood algorithm, we obtained appropriate variance-covariance structures for modeling heterogeneity, correlation of genetic effects and non-genetic residual effects. The estimates of genetic correlations indicated a strong association among earliness traits. The best linear unbiased prediction values were calculated for individual traits, and association analysis was performed in a panel of 95 diverse genotypes with 19 genic SSRs. Out of five QTL-flanking SSRs used here for validation, only ASSR295 could show significant association with FDR and Bonferroni corrections, and accounted for 15.4% IF, 14.2% DFF and 16.2% DM of phenotypic variance (PV). Remaining SSR markers (ASSR1486, ASSR206 and ASSR408) could not qualify false discovery rate (FDR) and Bonferroni criteria, hence declared as false positives. Additionally, we identified two highly significant SSR markers, ASSR8 and ASSR390 on LG 1 and LG 2, respectively. The SSR marker ASSR8 explained up to 22 and 11% PV for earliness traits and PB respectively, whereas ASSR390 controlled up to 17% PV for earliness traits. The validation and identification of new QTLs in pigeonpea across diverse genetic backgrounds brightens the prospects for marker-assisted selection to improve yield gains in pigeonpea.

Putting primary metabolism into perspective to obtain better fruits

Background

One of the key goals of fruit biology is to understand the factors that influence fruit growth and quality, ultimately with a view to manipulating them for improvement of fruit traits.

Scope

Primary metabolism, which is not only essential for growth but is also a major component of fruit quality, is an obvious target for improvement. However, metabolism is a moving target that undergoes marked changes throughout fruit growth and ripening.

Conclusions

Agricultural practice and breeding have successfully improved fruit metabolic traits, but both face the complexity of the interplay between development, metabolism and the environment. Thus, more fundamental knowledge is needed to identify further strategies for the manipulation of fruit metabolism. Nearly two decades of post-genomics approaches involving transcriptomics, proteomics and/or metabolomics have generated a lot of information about the behaviour of fruit metabolic networks. Today, the emergence of modelling tools is providing the opportunity to turn this information into a mechanistic understanding of fruits, and ultimately to design better fruits. Since high-quality data are a key requirement in modelling, a range of must-have parameters and variables is proposed.

Water Deficit and Salinity Stress Reveal Many Specific QTL for Plant Growth and Fruit Quality Traits in Tomato

Quality is a key trait in plant breeding, especially for fruit and vegetables. Quality involves several polygenic components, often influenced by environmental conditions with variable levels of genotype × environment interaction that must be considered in breeding strategies aiming to improve quality. In order to assess the impact of water deficit and salinity on tomato fruit quality, we evaluated a multi-parent advanced generation intercross (MAGIC) tomato population in contrasted environmental conditions over 2 years, one year in control vs. drought condition and the other in control vs. salt condition. Overall 250 individual lines from the MAGIC population-derived from eight parental lines covering a large diversity in cultivated tomato-were used to identify QTL in both experiments for fruit quality and yield component traits (fruit weight, number of fruit, Soluble Solid Content, firmness), phenology traits (time to flower and ripe) and a vegetative trait, leaf length. All the traits showed a large genotype variation (33-86% of total phenotypic variation) in both experiments and high heritability whatever the year or treatment. Significant genotype × treatment interactions were detected for five of the seven traits over the 2 years of experiments. QTL were mapped using 1,345 SNP markers. A total of 54 QTL were found among which 15 revealed genotype × environment interactions and 65% (35 QTL) were treatment specific. Confidence intervals of the QTL were projected on the genome physical map and allowed identifying regions carrying QTL co-localizations, suggesting pleiotropic regulation. We then applied a strategy for candidate gene detection based on the high resolution mapping offered by the MAGIC population, the allelic effect of each parental line at the QTL and the sequence information of the eight parental lines.

Profiling of drought-responsive microRNA and mRNA in tomato using high-throughput sequencing

BACKGROUND

Abiotic stresses cause severe loss of crop production. Among them, drought is one of the most frequent environmental stresses, which limits crop growth, development and productivity. Plant drought tolerance is fine-tuned by a complex gene regulatory network. Understanding the molecular regulation of this polygenic trait is crucial for the eventual success to improve plant yield and quality. Recent studies have demonstrated that microRNAs play critical roles in plant drought tolerance. However, little is known about the microRNA in drought response of the model plant tomato. Here, we described the profiling of drought-responsive microRNA and mRNA in tomato using high-throughput next-generation sequencing.

RESULTS

Drought stress was applied on the seedlings of M82, a drought-sensitive cultivated tomato genotype, and IL9-1, a drought-tolerant introgression line derived from the stress-resistant wild species Solanum pennellii LA0716 and M82. Under drought, IL9-1 performed superior than M82 regarding survival rate, H2O2 elimination and leaf turgor maintenance. A total of four small RNA and eight mRNA libraries were constructed and sequenced using Illumina sequencing technology. 105 conserved and 179 novel microRNAs were identified, among them, 54 and 98 were differentially expressed upon drought stress, respectively. The majority of the differentially-expressed conserved microRNAs was up-regulated in IL9-1 whereas down-regulated in M82. Under drought stress, 2714 and 1161 genes were found to be differentially expressed in M82 and IL9-1, respectively, and many of their homologues are involved in plant stress, such as genes encoding transcription factor and protein kinase. Various pathways involved in abiotic stress were revealed by Gene Ontology and pathway analysis. The mRNA sequencing results indicated that most of the target genes were regulated by their corresponding microRNAs, which suggested that microRNAs may play essential roles in the drought tolerance of tomato.

CONCLUSION

In this study, numerous microRNAs and mRNAs involved in the drought response of tomato were identified using high-throughput sequencing, which will provide new insights into the complex regulatory network of plant adaption to drought stress. This work will also help to exploit new players functioning in plant drought-stress tolerance.