Detection of major loci associated with the variation of 18 important agronomic traits between Solanum pimpinellifolium and cultivated tomatoes

A feedback loop at the THERMOSENSITIVE PARTHENOCARPY 4 locus controls tomato fruit set under heat stress

High temperatures compromise crop productivity worldwide, but breeding bottlenecks slow the delivery of climate-resilient crops. By investigating tomato fruit set under high temperatures, we discover a module comprising two linked genes, THERMOSENSITIVE PARTHENOCARPY 4a (TSP4a) and TSP4b, which encode the transcriptional regulators IAA9 and AINTEGUMENTA (ANT), respectively, to control thermosensitive parthenocarpy. TSP4a and TSP4b form a positive feedback loop upon heat stress to repress auxin signaling in ovaries. Natural TSP4a and TSP4b alleles bear regulatory-region polymorphisms and are differentially expressed to overcome the trade-off between fruit set and wider plant development. Gene editing of the TSP4a promoter and TSP4b 3' UTR in open-chromatin regions results in expression down-regulation, increased parthenocarpy without yield penalties and maintenance of fruit-sugar levels without broad auxin-related pleiotropic defects in greenhouse-grown plants. These mechanistic insights into heat-induced parthenocarpy and auxin signaling in reproductive organs demonstrate breeding utility to safeguard tomato yield under warming scenarios.

Manipulation of a strigolactone transporter in tomato confers resistance to the parasitic weed broomrape

Parasitic weeds of the Orobanchaceae family cause substantial economic losses and pose significant threats to global agriculture. However, management of such parasitism is challenging, and very few resistance genes have been cloned and characterized in depth. Here, we performed a genome-wide association study using 152 tomato accessions and identified SlABCG45 as a key gene that mediates host resistance to Phelipanche aegyptiaca by affecting the level of strigolactones (SLs) in root exudates. SLs are synthesized and released by host plants and act as germination stimulants for parasitic weeds. We found that SlABCG45 and its close homolog SlABCG44 were membrane-localized SL transporters with essential roles in exudation of SLs to the rhizosphere, resistance to Phelipanche and Orobanche, and upward transport of SLs from roots to shoots. As a predominant environmental stimulant exacerbates parasitism, phosphorus deficiency dramatically induced SlABCG45 expression and weakly induced SlABCG44 expression via the transcription factors SlNSP1 and SlNSP2. Knockout of SlABCG45 in tomato had little effect on yield traits in a broomrape-free field, but conferred increased resistance to different Phelipanche and Orobanche species, resulting in an ∼30% yield increase in a Phelipanche-infested field. Our findings reveal that targeting a single gene by genome editing can confer broad-spectrum parasite resistance in tomato, providing an effective strategy for the sustainable control of parasitic plants in agriculture.

Dissection of major QTLs and candidate genes for seedling stage salt/drought tolerance in tomato

BACKGROUND

As two of the most impactful abiotic stresses, salt and drought strongly affect tomato growth and development, especially at the seedling stage. However, dissection of the genetic basis underlying salt/drought tolerance at seedling stage in tomato remains limited in scope.

RESULTS

Here, we reported an analysis of major quantitative trait locus (QTL) and potential causal genetic variations in seedling stage salt/drought tolerance in recombinant inbred lines (n = 201) of S. pimpinellifolium and S. lycopersicum parents by whole genome resequencing. A total of 5 QTLs on chromosome 1, 3, 5, 7 and 12 for salt tolerance (ST) and 15 QTLs on chromosome 1, 3, 4, 8, 9, 10, 12 for drought tolerance (DT) were identified by linkage mapping. The proportion of phenotypic variation explained (PVE%) by these QTLs ranged from 4.91 to 15.86. Two major QTLs qST7 and qDT1-3 were detected in both two years, for which two candidate genes (methionine sulfoxide reductase SlMSRB1 and brassinosteroid insensitive 1-like receptor SlBRL1) and the potential functional variations were further analyzed. Taking advantage of the tomato population resequencing data, the frequency changes of the potential favorable QTL allele for seedling stage ST/DT during tomato breeding were explored.

CONCLUSIONS

These results will be beneficial for the exploration of salt/drought tolerance genes at seedling stages, laying a foundation for marker-assisted breeding for seedling stage salt/drought tolerance.

Genetic linkage analysis of stable QTLs in Gossypium hirsutum RIL population revealed function of GhCesA4 in fiber development

INTRODUCTION

Upland cotton is an important allotetrapolyploid crop providing natural fibers for textile industry. Under the present high-level breeding and production conditions, further simultaneous improvement of fiber quality and yield is facing unprecedented challenges due to their complex negative correlations.

OBJECTIVES

The study was to adequately identify quantitative trait loci (QTLs) and dissect how they orchestrate the formation of fiber quality and yield.

METHODS

A high-density genetic map (HDGM) based on an intraspecific recombinant inbred line (RIL) population consisting of 231 individuals was used to identify QTLs and QTL clusters of fiber quality and yield traits. The weighted gene correlation network analysis (WGCNA) package in R software was utilized to identify WGCNA network and hub genes related to fiber development. Gene functions were verified via virus-induced gene silencing (VIGS) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 strategies.

RESULTS

An HDGM consisting of 8045 markers was constructed spanning 4943.01 cM of cotton genome. A total of 295 QTLs were identified based on multi-environmental phenotypes. Among 139 stable QTLs, including 35 newly identified ones, seventy five were of fiber quality and 64 yield traits. A total of 33 QTL clusters harboring 74 QTLs were identified. Eleven candidate hub genes were identified via WGCNA using genes in all stable QTLs and QTL clusters. The relative expression profiles of these hub genes revealed their correlations with fiber development. VIGS and CRISPR/Cas9 edition revealed that the hub gene cellulose synthase 4 (GhCesA4, GH_D07G2262) positively regulate fiber length and fiber strength formation and negatively lint percentage.

CONCLUSION

Multiple analyses demonstrate that the hub genes harbored in the QTLs orchestrate the fiber development. The hub gene GhCesA4 has opposite pleiotropic effects in regulating trait formation of fiber quality and yield. The results facilitate understanding the genetic basis of negative correlation between cotton fiber quality and yield.

Deciphering salt stress responses in Solanum pimpinellifolium through high-throughput phenotyping

Soil salinity is a major environmental stressor affecting agricultural productivity worldwide. Understanding plant responses to salt stress is crucial for developing resilient crop varieties. Wild relatives of cultivated crops, such as wild tomato, Solanum pimpinellifolium, can serve as a useful resource to further expand the resilience potential of the cultivated germplasm, S. lycopersicum. In this study, we employed high-throughput phenotyping in the greenhouse and field conditions to explore salt stress responses of a S. pimpinellifolium diversity panel. Our study revealed extensive phenotypic variations in response to salt stress, with traits such as transpiration rate, shoot mass, and ion accumulation showing significant correlations with plant performance. We found that while transpiration was a key determinant of plant performance in the greenhouse, shoot mass strongly correlated with yield under field conditions. Conversely, ion accumulation was the least influential factor under greenhouse conditions. Through a Genome Wide Association Study, we identified candidate genes not previously associated with salt stress, highlighting the power of high-throughput phenotyping in uncovering novel aspects of plant stress responses. This study contributes to our understanding of salt stress tolerance in S. pimpinellifolium and lays the groundwork for further investigations into the genetic basis of these traits, ultimately informing breeding efforts for salinity tolerance in tomato and other crops.

A highly resolved nuclear phylogeny uncovers strong phylogenetic conservatism and correlated evolution of fruit color and size in Solanum L

Mutualisms between plants and fruit-eating animals were key to the radiation of angiosperms. Still, phylogenetic uncertainties limit our understanding of fleshy-fruit evolution, as in the case of Solanum, a genus with remarkable fleshy-fruit diversity, but with unresolved phylogenetic relationships. We used 1786 nuclear genes from 247 species, including 122 newly generated transcriptomes/genomes, to reconstruct the Solanum phylogeny and examine the tempo and mode of the evolution of fruit color and size. Our analysis resolved the backbone phylogeny of Solanum, providing high support for its clades. Our results pushed back the origin of Solanum to 53.1 million years ago (Ma), with most major clades diverging between 35 and 27 Ma. Evolution of Solanum fruit color and size revealed high levels of trait conservatism, where medium-sized berries that remain green when ripe are the likely ancestral form. Our analyses revealed that fruit size and color are evolutionary correlated, where dull-colored fruits are two times larger than black/purple and red fruits. We conclude that the strong phylogenetic conservatism shown in the color and size of Solanum fruits could limit the influences of fruit-eating animals on fleshy-fruit evolution. Our findings highlight the importance of phylogenetic constraints on the diversification of fleshy-fruit functional traits.

Exploring the gene expression network involved in the heat stress response of a thermotolerant tomato genotype

BACKGROUND

The increase in temperatures due to the current climate change dramatically affects crop cultivation, resulting in yield losses and altered fruit quality. Tomato is one of the most extensively grown and consumed horticultural products, and although it can withstand a wide range of climatic conditions, heat stress can affect plant growth and development specially on the reproductive stage, severely influencing the final yield. In the present work, the heat stress response mechanisms of one thermotolerant genotype (E42) were investigated by exploring its regulatory gene network. This was achieved through a promoter analysis based on the identification of the heat stress elements (HSEs) mapping in the promoters, combined with a gene co-expression network analysis aimed at identifying interactions among heat-related genes.

RESULTS

Results highlighted 82 genes presenting HSEs in the promoter and belonging to one of the 52 gene networks obtained by the GCN analysis; 61 of these also interact with heat shock factors (Hsfs). Finally, a list of 13 candidate genes including two Hsfs, nine heat shock proteins (Hsps) and two GDSL esterase/lipase (GELPs) were retrieved by focusing on those E42 genes exhibiting HSEs in the promoters, interacting with Hsfs and showing variants, compared to Heinz reference genome, with HIGH and/or MODERATE impact on the translated protein. Among these, the Gene Ontology annotation analysis evidenced that only LeHsp100 (Solyc02g088610) belongs to a network specifically involved in the response to heat stress.

CONCLUSIONS

As a whole, the combination of bioinformatic analyses carried out on genomic and trascriptomic data available for tomato, together with polymorphisms detected in HS-related genes of the thermotolerant E42 allowed to determine a subset of candidate genes involved in the HS response in tomato. This study provides a novel approach in the investigation of abiotic stress response mechanisms and further studies will be conducted to validate the role of the highlighted genes.

Tomato plant response to heat stress: a focus on candidate genes for yield-related traits

Climate change and global warming represent the main threats for many agricultural crops. Tomato is one of the most extensively grown and consumed horticultural products and can survive in a wide range of climatic conditions. However, high temperatures negatively affect both vegetative growth and reproductive processes, resulting in losses of yield and fruit quality traits. Researchers have employed different parameters to evaluate the heat stress tolerance, including evaluation of leaf- (stomatal conductance, net photosynthetic rate, Fv/Fm), flower- (inflorescence number, flower number, stigma exertion), pollen-related traits (pollen germination and viability, pollen tube growth) and fruit yield per plant. Moreover, several authors have gone even further, trying to understand the plants molecular response mechanisms to this stress. The present review focused on the tomato molecular response to heat stress during the reproductive stage, since the increase of temperatures above the optimum usually occurs late in the growing tomato season. Reproductive-related traits directly affects the final yield and are regulated by several genes such as transcriptional factors, heat shock proteins, genes related to flower, flowering, pollen and fruit set, and epigenetic mechanisms involving DNA methylation, histone modification, chromatin remodelling and non-coding RNAs. We provided a detailed list of these genes and their function under high temperature conditions in defining the final yield with the aim to summarize the recent findings and pose the attention on candidate genes that could prompt on the selection and constitution of new thermotolerant tomato plant genotypes able to face this abiotic challenge.

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.

Genomic Insights into the Origin of a Thermotolerant Tomato Line and Identification of Candidate Genes for Heat Stress

Climate change represents the main problem for agricultural crops, and the constitution of heat-tolerant genotypes is an important breeder’s strategy to reduce yield losses. The aim of the present study was to investigate the whole genome of a heat-tolerant tomato genotype (E42), in order to identify candidate genes involved in its response to high temperature. E42 presented a high variability for chromosomes 1, 4, 7 and 12, and phylogenetic analysis highlighted its relationship with the wild S. pimpinellifolium species. Variants with high (18) and moderate (139) impact on protein function were retrieved from two lists of genes related to heat tolerance and reproduction. This analysis permitted us to prioritize a subset of 35 candidate gene mapping in polymorphic regions, some colocalizing in QTLs controlling flowering in tomato. Among these genes, we identified 23 HSPs, one HSF, six involved in flowering and five in pollen activity. Interestingly, one gene coded for a flowering locus T1 and mapping on chromosome 11 resides in a QTL region controlling flowering and also showed 100% identity with an S. pimpinellifolium allele. This study provides useful information on both the E42 genetic background and heat stress response, and further studies will be conducted to validate these genes.

Molecular bases of rice grain size and quality for optimized productivity

The accomplishment of further optimization of crop productivity in grain yield and quality is a great challenge. Grain size is one of the crucial determinants of rice yield and quality; all of these traits are typical quantitative traits controlled by multiple genes. Research advances have revealed several molecular and developmental pathways that govern these traits of agronomical importance. This review provides a comprehensive summary of these pathways, including those mediated by G-protein, the ubiquitin-proteasome system, mitogen-activated protein kinase, phytohormone, transcriptional regulators, and storage product biosynthesis and accumulation. We also generalize the excellent precedents for rice variety improvement of grain size and quality, which utilize newly developed gene editing and conventional gene pyramiding capabilities. In addition, we discuss the rational and accurate breeding strategies, with the aim of better applying molecular design to breed high-yield and superior-quality varieties.

Targeting a gene regulatory element enhances rice grain yield by decoupling panicle number and size

Crop genetic improvement requires balancing complex tradeoffs caused by gene pleiotropy and linkage drags, as exemplified by IPA1 (Ideal Plant Architecture 1), a typical pleiotropic gene in rice that increases grains per panicle but reduces tillers. In this study, we identified a 54-base pair cis-regulatory region in IPA1 via a tiling-deletion-based CRISPR-Cas9 screen that, when deleted, resolves the tradeoff between grains per panicle and tiller number, leading to substantially enhanced grain yield per plant. Mechanistic studies revealed that the deleted fragment is a target site for the transcription factor An-1 to repress IPA1 expression in panicles and roots. Targeting gene regulatory regions should help dissect tradeoff effects and provide a rich source of targets for breeding complementary beneficial traits.

SISTER OF TM3 activates FRUITFULL1 to regulate inflorescence branching in tomato

Selection for favorable inflorescence architecture to improve yield is one of the crucial targets in crop breeding. Different tomato varieties require distinct inflorescence-branching structures to enhance productivity. While a few important genes for tomato inflorescence-branching development have been identified, the regulatory mechanism underlying inflorescence branching is still unclear. Here, we confirmed that SISTER OF TM3 (STM3), a homolog of Arabidopsis SOC1, is a major positive regulatory factor of tomato inflorescence architecture by map-based cloning. High expression levels of STM3 underlie the highly inflorescence-branching phenotype in ST024. STM3 is expressed in both vegetative and reproductive meristematic tissues and in leaf primordia and leaves, indicative of its function in flowering time and inflorescence-branching development. Transcriptome analysis shows that several floral development-related genes are affected by STM3 mutation. Among them, FRUITFULL1 (FUL1) is downregulated in stm3cr mutants, and its promoter is bound by STM3 by ChIP-qPCR analysis. EMSA and dual-luciferase reporter assays further confirmed that STM3 could directly bind the promoter region to activate FUL1 expression. Mutation of FUL1 could partially restore inflorescence-branching phenotypes caused by high STM3 expression in ST024. Our findings provide insights into the molecular and genetic mechanisms underlying inflorescence development in tomato.

Tomato genomic prediction for good performance under high-temperature and identification of loci involved in thermotolerance response

Many studies showed that few degrees above tomato optimum growth temperature threshold can lead to serious loss in production. Therefore, the development of innovative strategies to obtain tomato cultivars with improved yield under high temperature conditions is a main goal both for basic genetic studies and breeding activities. In this paper, a F4 segregating population was phenotypically evaluated for quantitative and qualitative traits under heat stress conditions. Moreover, a genotyping by sequencing (GBS) approach has been employed for building up genomic selection (GS) models both for yield and soluble solid content (SCC). Several parameters, including training population size, composition and marker quality were tested to predict genotype performance under heat stress conditions. A good prediction accuracy for the two analyzed traits (0.729 for yield production and 0.715 for SCC) was obtained. The predicted models improved the genetic gain of selection in the next breeding cycles, suggesting that GS approach is a promising strategy to accelerate breeding for heat tolerance in tomato. Finally, the annotation of SNPs located in gene body regions combined with QTL analysis allowed the identification of five candidates putatively involved in high temperatures response, and the building up of a GS model based on calibrated panel of SNP markers.

Two zinc-finger roteins control the initiation and elongation of long stalk trichomes in tomato

Plant glandular trichomes are epidermal secretory structures that are important for plant resistance to pests. Although several regulatory genes have been characterized in trichome development, the molecular mechanisms conferring glandular trichome morphogenesis are unclear. We observed the differences in trichomes in cultivated tomato cv. 'Moneymaker' (MM) and the wild species Solanum pimpinellifolium PI365967 (PP), and used a recombinant inbred line (RIL) population to identify the genes that control trichome development in tomato. We found that the genomic variations in two genes, H and SH, contribute to the trichome differences between MM and PP. H and SH encode two paralogous C2H2 zinc-finger proteins that function redundantly in regulating trichome formation. Loss-of-function h/sh double mutants exhibited a significantly decreased number of Type I trichomes and complete loss of long stalk trichomes. Molecular and genetic analyses further indicate that H and SH act upstream of ZFP5. Overexpression of ZFP5 partially restored the trichome defects in NIL-h^PPsh^PP. Moreover, H and SH expression is induced by high temperatures, and their mutations inhibit the elongation of trichomes that reduce the plant repellent to whiteflies. Our findings confirm that H and SH are two vital transcription factors controlling initiation and elongation of Type I and III multicellular trichomes in tomato.

Genome-wide association study reveals the genetic architecture of 27 agronomic traits in tomato

Tomato (Solanum lycopersicum) is a highly valuable fruit crop, and yield is one of the most important agronomic traits. However, the genetic architecture underlying tomato yield-related traits has not been fully addressed. Based on ∼4.4 million single nucleotide polymorphisms obtained from 605 diverse accessions, we performed a comprehensive genome-wide association study for 27 agronomic traits in tomato. A total of 239 significant associations corresponding to 129 loci, harboring many previously reported and additional genes related to vegetative and reproductive development, were identified, and these loci explained an average of ∼8.8% of the phenotypic variance. A total of 51 loci associated with 25 traits have been under selection during tomato domestication and improvement. Furthermore, a candidate gene, Sl-ACTIVATED MALATE TRANSPORTER15, that encodes an aluminum-activated malate transporter was functionally characterized and shown to act as a pivotal regulator of leaf stomata formation, thereby affecting photosynthesis and drought resistance. This study provides valuable information for tomato genetic research and breeding.

Tomato fruit quality traits and metabolite content are affected by reciprocal crosses and heterosis

Heterosis occurs when the F1s outperform their parental lines for a trait. Reciprocal hybrids are obtained by changing the cross direction of parental genotypes. Both biological phenomena could affect the external and internal attributes of fleshy fruits. This work aimed to detect reciprocal effects and heterosis in tomato (Solanum lycopersicum) fruit quality traits and metabolite content. Twelve agronomic traits and 28 metabolites identified and estimated by 1H-NMR were evaluated in five cultivars grown in two environments. Given that the genotype component was more important than the phenotype, the traits were evaluated following a full diallel mating design among those cultivars, in a greenhouse. Hybrids showed a higher phenotypic diversity than parental lines. Interestingly, the metabolites, mainly amino acids, displayed more reciprocal effects and heterosis. Agronomic traits were more influenced by general combining ability (GCA) and metabolites by specific combining ability (SCA). Furthermore, the genetic distance between parental lines was not causally related to the occurrence of reciprocal effects or heterosis. Hybrids with heterosis and a high content of metabolites linked to tomato flavour and nutritious components were obtained. Our results highlight the impact of selecting a cultivar as male or female in a cross to enhance the variability of fruit attributes through hybrids as well as the possibility to exploit heterosis for fruit composition.

Genetic dissection of an allotetraploid interspecific CSSLs guides interspecific genetics and breeding in cotton

Background

The low genetic diversity of Upland cotton limits the potential for genetic improvement. Making full use of the genetic resources of Sea-island cotton will facilitate genetic improvement of widely cultivated Upland cotton varieties. The chromosome segments substitution lines (CSSLs) provide an ideal strategy for mapping quantitative trait loci (QTL) in interspecific hybridization.

Results

In this study, a CSSL population was developed by PCR-based markers assisted selection (MAS), derived from the crossing and backcrossing of Gossypium hirsutum (Gh) and G. barbadense (Gb), firstly. Then, by whole genome re-sequencing, 11,653,661 high-quality single nucleotide polymorphisms (SNPs) were identified which ultimately constructed 1211 recombination chromosome introgression segments from Gb. The sequencing-based physical map provided more accurate introgressions than the PCR-based markers. By exploiting CSSLs with mutant morphological traits, the genes responding for leaf shape and fuzz-less mutation in the Gb were identified. Based on a high-resolution recombination bin map to uncover genetic loci determining the phenotypic variance between Gh and Gb, 64 QTLs were identified for 14 agronomic traits with an interval length of 158 kb to 27 Mb. Surprisingly, multiple alleles of Gb showed extremely high value in enhancing cottonseed oil content (SOC).

Conclusions

This study provides guidance for studying interspecific inheritance, especially breeding researchers, for future studies using the traditional PCR-based molecular markers and high-throughput re-sequencing technology in the study of CSSLs. Available resources include candidate position for controlling cotton quality and quantitative traits, and excellent breeding materials. Collectively, our results provide insights into the genetic effects of Gb alleles on the Gh, and provide guidance for the utilization of Gb alleles in interspecific breeding.

High-Throughput Genotyping of Resilient Tomato Landraces to Detect Candidate Genes Involved in the Response to High Temperatures

The selection of tolerant varieties is a powerful strategy to ensure highly stable yield under elevated temperatures. In this paper, we report the phenotypic and genotypic characterization of 10 tomato landraces to identify the best performing under high temperatures. The phenotyping of five yield-related traits allowed us to select one genotype that exhibits highly stable yield performances in different environmental conditions. Moreover, a Genotyping-by-Sequencing approach allowed us to explore the genetic variability of the tested genotypes. The high and stable yielding landrace E42 was the most polymorphic one, with ~49% and ~47% private SNPs and InDels, respectively. The effect of 26,113 mutations on proteins’ structure was investigated and it was discovered that 37 had a high impact on the structure of 34 proteins of which some are putatively involved in responses to high temperatures. Additionally, 129 polymorphic sequences aligned against tomato wild species genomes revealed the presence in the genotype E42 of several introgressed regions deriving from S. pimpinellifolium. The position on the tomato map of genes affected by moderate and high impact mutations was also compared with that of known markers/QTLs (Quantitative Trait Loci) associated with reproductive and yield-related traits. The candidate genes/QTLs regulating heat tolerance in the selected landrace E42 could be further investigated to better understand the genetic mechanisms controlling traits for high and stable yield trait under high temperatures.

Dissection of complex traits of tomato in the post-genome era

We present the main advances of dissection of complex traits in tomato by omics, the genes identified to control complex traits and the application of CRISPR/Cas9 in tomato breeding. Complex traits are believed to be under the control of multiple genes, each with different effects and interaction with environmental factors. Advance development of sequencing and molecular technologies has enabled the recognition of the genomic structure of most organisms and the identification of a nearly limitless number of markers that have made it to accelerate the speed of QTL identification and gene cloning. Meanwhile, multiomics have been used to identify the genetic variations among different tomato species, determine the expression profiles of genes in different tissues and at distinct developmental stages, and detect metabolites in different pathways and processes. The combination of these data facilitates to reveal mechanism underlying complex traits. Moreover, mutants generated by mutagens and genome editing provide relatively rich genetic variation for deciphering the complex traits and exploiting them in tomato breeding. In this article, we present the main advances of complex trait dissection in tomato by omics since the release of the tomato genome sequence in 2012. We provide further insight into some tomato complex traits because of the causal genetic variations discovered so far and explore the utilization of CRISPR/Cas9 for the modification of tomato complex traits.

Benchmarking variant identification tools for plant diversity discovery

BACKGROUND

The ability to accurately and comprehensively identify genomic variations is critical for plant studies utilizing high-throughput sequencing. Most bioinformatics tools for processing next-generation sequencing data were originally developed and tested in human studies, raising questions as to their efficacy for plant research. A detailed evaluation of the entire variant calling pipeline, including alignment, variant calling, variant filtering, and imputation was performed on different programs using both simulated and real plant genomic datasets.

RESULTS

A comparison of SOAP2, Bowtie2, and BWA-MEM found that BWA-MEM was consistently able to align the most reads with high accuracy, whereas Bowtie2 had the highest overall accuracy. Comparative results of GATK HaplotypCaller versus SAMtools mpileup indicated that the choice of variant caller affected precision and recall differentially depending on the levels of diversity, sequence coverage and genome complexity. A cross-reference experiment of S. lycopersicum and S. pennellii reference genomes revealed the inadequacy of single reference genome for variant discovery that includes distantly-related plant individuals. Machine-learning-based variant filtering strategy outperformed the traditional hard-cutoff strategy resulting in higher number of true positive variants and fewer false positive variants. A 2-step imputation method, which utilized a set of high-confidence SNPs as the reference panel, showed up to 60% higher accuracy than direct LD-based imputation.

CONCLUSIONS

Programs in the variant discovery pipeline have different performance on plant genomic dataset. Choice of the programs is subjected to the goal of the study and available resources. This study serves as an important guiding information for plant biologists utilizing next-generation sequencing data for diversity characterization and crop improvement.