Identification of major loci and genomic regions controlling acid and volatile content in tomato fruit: implications for flavor improvement.

Plant metabolites are important to world food security due to their roles in crop yield and nutritional quality. Here we report the metabolic profile of 300 tomato accessions (Solanum lycopersicum and related wild species) by quantifying 60 primary and secondary metabolites, including volatile organic compounds, over a period of 2 yr. Metabolite content and genetic inheritance of metabolites varied broadly, both within and between different genetic groups. Using genotype information gained from 10 000 single nucleotide polymorphism markers, we performed a metabolite genome-wide association mapping (GWAS) study. We identified 79 associations influencing 13 primary and 19 secondary metabolites with large effects at high resolution. Four genome regions were detected, highlighting clusters of associations controlling the variation of several metabolites. Local linkage disequilibrium analysis and allele mining identified possible candidate genes which may modulate the content of metabolites that are of significant importance for human diet and fruit consumption. We precisely characterized two associations involved in fruit acidity and phenylpropanoid volatile production. Taken together, this study reveals complex and distinct metabolite regulation in tomato subspecies and demonstrates that GWAS is a powerful tool for gene-metabolite annotation and identification, pathways elucidation, and further crop improvement.

Use of modern tomato breeding germplasm for deciphering the genetic control of agronomical traits by Genome Wide Association study.

KEY MESSAGE: A panel of 300 tomato accessions including breeding materials was built and characterized with >11,000 SNP. A population structure in six subgroups was identified. Strong heterogeneity in linkage disequilibrium and recombination landscape among groups and chromosomes was shown. GWAS identified several associations for fruit weight, earliness and plant growth. Genome-wide association studies (GWAS) have become a method of choice in quantitative trait dissection. First limited to highly polymorphic and outcrossing species, it is now applied in horticultural crops, notably in tomato. Until now GWAS in tomato has been performed on panels of heirloom and wild accessions. Using modern breeding materials would be of direct interest for breeding purpose. To implement GWAS on a large panel of 300 tomato accessions including 168 breeding lines, this study assessed the genetic diversity and linkage disequilibrium decay and revealed the population structure and performed GWA experiment. Genetic diversity and population structure analyses were based on molecular markers (>11,000 SNP) covering the whole genome. Six genetic subgroups were revealed and associated to traits of agronomical interest, such as fruit weight and disease resistance. Estimates of linkage disequilibrium highlighted the heterogeneity of its decay among genetic subgroups. Haplotype definition allowed a fine characterization of the groups and their recombination landscape revealing the patterns of admixture along the genome. Selection footprints showed results in congruence with introgressions. Taken together, all these elements refined our knowledge of the genetic material included in this panel and allowed the identification of several associations for fruit weight, plant growth and earliness, deciphering the genetic architecture of these complex traits and identifying several new loci useful for tomato breeding.

Genes involved in floral meristem in tomato exhibit drastically reduced genetic diversity and signature of selection.

BACKGROUND: Domestication and selection of crops have notably reshaped fruit morphology. With its large phenotypic diversity, tomato (Solanum lycopersicum) illustrates this evolutive trend. Genes involved in flower meristem development are known to regulate also fruit morphology. To decipher the genetic variation underlying tomato fruit morphology, we assessed the nucleotide diversity and selection footprints of candidate genes involved in flower and fruit development and performed genome-wide association studies. RESULTS: Thirty candidate genes were selected according to their similarity with genes involved in meristem development or their known causal function in Arabidopsis thaliana. In tomato, these genes and flanking regions were sequenced in a core collection of 96 accessions (including cultivated, cherry-type and wild relative accessions) maximizing the molecular diversity, using the Roche 454 technology. A total amount of 17 Mb was sequenced allowing the discovery of 6,106 single nucleotide polymorphisms (SNPs). The annotation of the 30 gene regions identified 231 exons carrying 517 SNPs. Subsequently, the nucleotide diversity (π) and the neutral evolution of each region were compared against genome-wide values within the collection, using a SNP array carrying 7,667 SNPs mainly distributed in coding sequences.About half of the genes revealed footprints of selection and polymorphisms putatively involved in fruit size variation by showing negative Tajima's D and nucleotide diversity reduction in cultivated tomato compared to its wild relative. Among the candidates, FW2.2 and BAM1 sequences revealed selection footprints within their promoter regions suggesting their potential involvement in their regulation. Two associations co-localized with previously identified loci: LC (locule number) and Ovate (fruit shape). CONCLUSION: Compared to whole genome genotypic data, a drastic reduction of nucleotide diversity was shown for several candidate genes. Strong selection patterns were identified in 15 candidates highlighting the critical role of meristem maintenance genes as well as the impact of domestication on candidates. The study highlighted a set of polymorphisms putatively important in the evolution of these genes.

High-resolution mapping of a fruit firmness-related quantitative trait locus in tomato reveals epistatic interactions associated with a complex combinatorial locus.

Fruit firmness in tomato (Solanum lycopersicum) is determined by a number of factors including cell wall structure, turgor, and cuticle properties. Firmness is a complex polygenic trait involving the coregulation of many genes and has proved especially challenging to unravel. In this study, a quantitative trait locus (QTL) for fruit firmness was mapped to tomato chromosome 2 using the Zamir Solanum pennellii interspecific introgression lines (ILs) and fine-mapped in a population consisting of 7,500 F2 and F3 lines from IL 2-3 and IL 2-4. This firmness QTL contained five distinct subpeaks, Fir(s.p.)QTL2.1 to Fir(s.p.)QTL2.5, and an effect on a distal region of IL 2-4 that was nonoverlapping with IL 2-3. All these effects were located within an 8.6-Mb region. Using genetic markers, each subpeak within this combinatorial locus was mapped to a physical location within the genome, and an ethylene response factor (ERF) underlying Fir(s.p.)QTL2.2 and a region containing three pectin methylesterase (PME) genes underlying Fir(s.p.)QTL2.5 were nominated as QTL candidate genes. Statistical models used to explain the observed variability between lines indicated that these candidates and the nonoverlapping portion of IL 2-4 were sufficient to account for the majority of the fruit firmness effects. Quantitative reverse transcription-polymerase chain reaction was used to quantify the expression of each candidate gene. ERF showed increased expression associated with soft fruit texture in the mapping population. In contrast, PME expression was tightly linked with firm fruit texture. Analysis of a range of recombinant lines revealed evidence for an epistatic interaction that was associated with this combinatorial locus.

Diploid/polyploid syntenic shuttle mapping and haplotype-specific chromosome walking toward a rust resistance gene (Bru1) in highly polyploid sugarcane (2n approximately 12x approximately 115).

The genome of modern sugarcane cultivars is highly polyploid (approximately 12x), aneuploid, of interspecific origin, and contains 10 Gb of DNA. Its size and complexity represent a major challenge for the isolation of agronomically important genes. Here we report on the first attempt to isolate a gene from sugarcane by map-based cloning, targeting a durable major rust resistance gene (Bru1). We describe the genomic strategies that we have developed to overcome constraints associated with high polyploidy in the successive steps of map-based cloning approaches, including diploid/polyploid syntenic shuttle mapping with two model diploid species (sorghum and rice) and haplotype-specific chromosome walking. Their applications allowed us (i) to develop a high-resolution map including markers at 0.28 and 0.14 cM on both sides and 13 markers cosegregating with Bru1 and (ii) to develop a physical map of the target haplotype that still includes two gaps at this stage due to the discovery of an insertion specific to this haplotype. These approaches will pave the way for the development of future map-based cloning approaches for sugarcane and other complex polyploid species.

Sugarcane genomics: depicting the complex genome of an important tropical crop.

In the past few years, approaches such as molecular cytogenetics and the use of molecular markers have permitted significant advances in the establishment of the evolutionary origin and genome structure of sugarcane, an important polyploid crop. The availability of new resources, such as a bacterial artificial chromosome library and a huge collection of expressed sequence tags, has opened the gateway to promising functional analyses on a genomic scale.

Orthologous comparison in a gene-rich region among grasses reveals stability in the sugarcane polyploid genome.

Modern sugarcane (Saccharum spp.) is an important grass that contributes 60% of the raw sugar produced worldwide and has a high biofuel production potential. It was created about a century ago through hybridization of two highly polyploid species, namely S. officinarum and S. spontaneum. We investigated genome dynamics in this highly polyploid context by analyzing two homoeologous sequences (97 and 126 kb) in a region that has already been studied in several cereals. Our findings indicate that the two Saccharum species diverged 1.5-2 million years ago from one another and 8-9 million years ago from sorghum. The two sugarcane homoeologous haplotypes show perfect colinearity as well as high gene structure conservation. Apart from the insertion of a few retrotransposable elements, high homology was also observed for the non-transcribed regions. Relative to sorghum, the sugarcane sequences displayed colinearity, with the exception of two genes present only in sorghum, and striking homology in most non-coding parts of the genome. The gene distribution highlighted high synteny and colinearity with rice, and partial colinearity with each homoeologous maize region, which became perfect when the sequences were combined. The haplotypes observed in sugarcane may thus closely represent the ancestral Andropogoneae haplotype. This analysis of sugarcane haplotype organization at the sequence level suggests that the high ploidy in sugarcane did not induce generalized reshaping of its genome, thus challenging the idea that polyploidy quickly induces generalized rearrangement of genomes. These results also confirm the view that sorghum is the model of choice for sugarcane.

Evaluation of monocot and eudicot divergence using the sugarcane transcriptome.

Over 40,000 sugarcane (Saccharum officinarum) consensus sequences assembled from 237,954 expressed sequence tags were compared with the protein and DNA sequences from other angiosperms, including the genomes of Arabidopsis and rice (Oryza sativa). Approximately two-thirds of the sugarcane transcriptome have similar sequences in Arabidopsis. These sequences may represent a core set of proteins or protein domains that are conserved among monocots and eudicots and probably encode for essential angiosperm functions. The remaining sequences represent putative monocot-specific genetic material, one-half of which were found only in sugarcane. These monocot-specific cDNAs represent either novelties or, in many cases, fast-evolving sequences that diverged substantially from their eudicot homologs. The wide comparative genome analysis presented here provides information on the evolutionary changes that underlie the divergence of monocots and eudicots. Our comparative analysis also led to the identification of several not yet annotated putative genes and possible gene loss events in Arabidopsis.

Analysis and functional annotation of an expressed sequence tag collection for tropical crop sugarcane.

To contribute to our understanding of the genome complexity of sugarcane, we undertook a large-scale expressed sequence tag (EST) program. More than 260,000 cDNA clones were partially sequenced from 26 standard cDNA libraries generated from different sugarcane tissues. After the processing of the sequences, 237,954 high-quality ESTs were identified. These ESTs were assembled into 43,141 putative transcripts. Of the assembled sequences, 35.6% presented no matches with existing sequences in public databases. A global analysis of the whole SUCEST data set indicated that 14,409 assembled sequences (33% of the total) contained at least one cDNA clone with a full-length insert. Annotation of the 43,141 assembled sequences associated almost 50% of the putative identified sugarcane genes with protein metabolism, cellular communication/signal transduction, bioenergetics, and stress responses. Inspection of the translated assembled sequences for conserved protein domains revealed 40,821 amino acid sequences with 1415 Pfam domains. Reassembling the consensus sequences of the 43,141 transcripts revealed a 22% redundancy in the first assembling. This indicated that possibly 33,620 unique genes had been identified and indicated that >90% of the sugarcane expressed genes were tagged.