Graph pangenome captures missing heritability and empowers tomato breeding.

Missing heritability in genome-wide association studies defines a major problem in genetic analyses of complex biological traits1,2. The solution to this problem is to identify all causal genetic variants and to measure their individual contributions3,4. Here we report a graph pangenome of tomato constructed by precisely cataloguing more than 19 million variants from 838 genomes, including 32 new reference-level genome assemblies. This graph pangenome was used for genome-wide association study analyses and heritability estimation of 20,323 gene-expression and metabolite traits. The average estimated trait heritability is 0.41 compared with 0.33 when using the single linear reference genome. This 24% increase in estimated heritability is largely due to resolving incomplete linkage disequilibrium through the inclusion of additional causal structural variants identified using the graph pangenome. Moreover, by resolving allelic and locus heterogeneity, structural variants improve the power to identify genetic factors underlying agronomically important traits leading to, for example, the identification of two new genes potentially contributing to soluble solid content. The newly identified structural variants will facilitate genetic improvement of tomato through both marker-assisted selection and genomic selection. Our study advances the understanding of the heritability of complex traits and demonstrates the power of the graph pangenome in crop breeding.

Full dynamic range proteome analysis of S. cerevisiae by targeted proteomics.

The rise of systems biology implied a growing demand for highly sensitive techniques for the fast and consistent detection and quantification of target sets of proteins across multiple samples. This is only partly achieved by classical mass spectrometry or affinity-based methods. We applied a targeted proteomics approach based on selected reaction monitoring (SRM) to detect and quantify proteins expressed to a concentration below 50 copies/cell in total S. cerevisiae digests. The detection range can be extended to single-digit copies/cell and to proteins undetected by classical methods. We illustrate the power of the technique by the consistent and fast measurement of a network of proteins spanning the entire abundance range over a growth time course of S. cerevisiae transiting through a series of metabolic phases. We therefore demonstrate the potential of SRM-based proteomics to provide assays for the measurement of any set of proteins of interest in yeast at high-throughput and quantitative accuracy.

A Solanum lycopersicoides reference genome facilitates insights into tomato specialized metabolism and immunity.

Wild relatives of tomato are a valuable source of natural variation in tomato breeding, as many can be hybridized to the cultivated species (Solanum lycopersicum). Several, including Solanum lycopersicoides, have been crossed to S. lycopersicum for the development of ordered introgression lines (ILs), facilitating breeding for desirable traits. Despite the utility of these wild relatives and their associated ILs, few finished genome sequences have been produced to aid genetic and genomic studies. Here we report a chromosome-scale genome assembly for S. lycopersicoides LA2951, which contains 37 938 predicted protein-coding genes. With the aid of this genome assembly, we have precisely delimited the boundaries of the S. lycopersicoides introgressions in a set of S. lycopersicum cv. VF36 × LA2951 ILs. We demonstrate the usefulness of the LA2951 genome by identifying several quantitative trait loci for phenolics and carotenoids, including underlying candidate genes, and by investigating the genome organization and immunity-associated function of the clustered Pto gene family. In addition, syntenic analysis of R2R3MYB genes sheds light on the identity of the Aubergine locus underlying anthocyanin production. The genome sequence and IL map provide valuable resources for studying fruit nutrient/quality traits, pathogen resistance, and environmental stress tolerance. We present a new genome resource for the wild species S. lycopersicoides, which we use to shed light on the Aubergine locus responsible for anthocyanin production. We also provide IL boundary mappings, which facilitated identifying novel carotenoid quantitative trait loci of which one was likely driven by an uncharacterized lycopene ß-cyclase whose function we demonstrate.

Metabolic networks of the Nicotiana genus in the spotlight: content, progress and outlook.

Manually curated metabolic databases residing at the Sol Genomics Network comprise two taxon-specific databases for the Solanaceae family, i.e. SolanaCyc and the genus Nicotiana, i.e. NicotianaCyc as well as six species-specific databases for Nicotiana tabacum TN90, N. tabacum K326, Nicotiana benthamiana, N. sylvestris, N. tomentosiformis and N. attenuata. New pathways were created through the extraction, examination and verification of related data from the literature and the aid of external database guided by an expert-led curation process. Here we describe the curation progress that has been achieved in these databases since the first release version 1.0 in 2016, the curation flow and the curation process using the example metabolic pathway for cholesterol in plants. The current content of our databases comprises 266 pathways and 36 superpathways in SolanaCyc and 143 pathways plus 21 superpathways in NicotianaCyc, manually curated and validated specifically for the Solanaceae family and Nicotiana genus, respectively. The curated data have been propagated to the respective Nicotiana-specific databases, which resulted in the enrichment and more accurate presentation of their metabolic networks. The quality and coverage in those databases have been compared with related external databases and discussed in terms of literature support and metabolic content.

Spatiotemporal dynamics of the tomato fruit transcriptome under prolonged water stress.

Water availability influences all aspects of plant growth and development; however, most studies of plant responses to drought have focused on vegetative organs, notably roots and leaves. Far less is known about the molecular bases of drought acclimation responses in fruits, which are complex organs with distinct tissue types. To obtain a more comprehensive picture of the molecular mechanisms governing fruit development under drought, we profiled the transcriptomes of a spectrum of fruit tissues from tomato (Solanum lycopersicum), spanning early growth through ripening and collected from plants grown under varying intensities of water stress. In addition, we compared transcriptional changes in fruit with those in leaves to highlight different and conserved transcriptome signatures in vegetative and reproductive organs. We observed extensive and diverse genetic reprogramming in different fruit tissues and leaves, each associated with a unique response to drought acclimation. These included major transcriptional shifts in the placenta of growing fruit and in the seeds of ripe fruit related to cell growth and epigenetic regulation, respectively. Changes in metabolic and hormonal pathways, such as those related to starch, carotenoids, jasmonic acid, and ethylene metabolism, were associated with distinct fruit tissues and developmental stages. Gene coexpression network analysis provided further insights into the tissue-specific regulation of distinct responses to water stress. Our data highlight the spatiotemporal specificity of drought responses in tomato fruit and indicate known and unrevealed molecular regulatory mechanisms involved in drought acclimation, during both vegetative and reproductive stages of development.

Connecting data for consumer preferences, food quality, and breeding in support of market-oriented breeding of root, tuber, and banana crops.

The 5-year project 'Breeding roots, tubers and banana products for end user preferences' (RTBfoods) focused on collecting consumers' preferences on 12 food products to guide breeding programmes. It involved multidisciplinary teams from Africa, Latin America, and Europe. Diverse data types were generated on preferred qualities of users (farmers, family and entrepreneurial processors, traders or retailers, and consumers). Country-based target product profiles were produced with a comprehensive market analysis, disaggregating gender's role and preferences, providing prioritised lists of traits for the development of new plant varieties. We describe the approach taken to create, in the roots, tubers, and banana breeding databases, a centralised and meaningful open access to sensory information on food products and genotypes. Biochemical, instrumental textural, and sensory analysis data are then directly connected to the specific plant record while user survey data, bearing personal information, were analysed, anonymised, and uploaded in a repository. Names and descriptions of food quality traits were added into the Crop Ontology for labelling data in the databases, along with the various methods of measurement used by the project. The development and application of standard operating procedures, data templates, and adapted trait ontologies improved the data quality and its format, enabling the linking of these to the plant material studied when uploaded in the breeding databases or in repositories. Some modifications to the database model were necessary to accommodate the food sensory traits and sensory panel trials. © 2023 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Large-scale genome-wide association study, using historical data, identifies conserved genetic architecture of cyanogenic glucoside content in cassava (Manihot esculenta Crantz) root.

Manihot esculenta (cassava) is a root crop originating from South America that is a major staple in the tropics, including in marginal environments. This study focused on South American and African germplasm and investigated the genetic architecture of hydrogen cyanide (HCN), a major component of root quality. HCN, representing total cyanogenic glucosides, is a plant defense component against herbivory but is also toxic for human consumption. We genotyped 3354 landraces and modern breeding lines originating from 26 Brazilian states and 1389 individuals were phenotypically characterized across multi-year trials for HCN. All plant material was subjected to high-density genotyping using genotyping by sequencing. We performed genome-wide association mapping to characterize the genetic architecture and gene mapping of HCN. Field experiments revealed strong broad- and narrow-sense trait heritability (0.82 and 0.41, respectively). Two major loci were identified, encoding for an ATPase and a MATE protein, and contributing up to 7 and 30% of the HCN concentration in roots, respectively. We developed diagnostic markers for breeding applications, validated trait architecture consistency in African germplasm and investigated further evidence for the domestication of sweet and bitter cassava. Fine genomic characterization revealed: (i) the major role played by vacuolar transporters in regulating HCN content; (ii) the co-domestication of sweet and bitter cassava major alleles are dependent upon geographical zone; and (iii) the major loci allele for high HCN in M. esculenta Crantz seems to originate from its ancestor, M. esculenta subsp. flabellifolia. Taken together, these findings expand our insights into cyanogenic glucosides in cassava roots and its glycosylated derivatives in plants.

Current status and impending progress for cassava structural genomics.

KEY MESSAGE: We demystify recent advances in genome assemblies for the heterozygous staple crop cassava (Manihot esculenta), and highlight key cassava genomic resources. Cassava, Manihot esculenta Crantz, is a crop of societal and agricultural importance in tropical regions around the world. Genomics provides a platform for accelerated improvement of cassava's nutritional and agronomic traits, as well as for illuminating aspects of cassava's history including its path towards domestication. The highly heterozygous nature of the cassava genome is widely recognized. However, the full extent and context of this heterozygosity has been difficult to reveal because of technological limitations within genome sequencing. Only recently, with several new long-read sequencing technologies coming online, has the genomics community been able to tackle some similarly difficult genomes. In light of these recent advances, we provide this review to document the current status of the cassava genome and genomic resources and provide a perspective on what to look forward to in the coming years.

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens.

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de.

Metabolic profiles of six African cultivars of cassava (Manihot esculenta Crantz) highlight bottlenecks of root yield.

Cassava is an important staple crop in sub-Saharan Africa, due to its high productivity even on nutrient poor soils. The metabolic characteristics underlying this high productivity are poorly understood including the mode of photosynthesis, reasons for the high rate of photosynthesis, the extent of source/sink limitation, the impact of environment, and the extent of variation between cultivars. Six commercial African cassava cultivars were grown in a greenhouse in Erlangen, Germany, and in the field in Ibadan, Nigeria. Source leaves, sink leaves, stems and storage roots were harvested during storage root bulking and analyzed for sugars, organic acids, amino acids, phosphorylated intermediates, minerals, starch, protein, activities of enzymes in central metabolism and yield traits. High ratios of RuBisCO:phosphoenolpyruvate carboxylase activity support a C3 mode of photosynthesis. The high rate of photosynthesis is likely to be attributed to high activities of enzymes in the Calvin-Benson cycle and pathways for sucrose and starch synthesis. Nevertheless, source limitation is indicated because root yield traits correlated with metabolic traits in leaves rather than in the stem or storage roots. This situation was especially so in greenhouse-grown plants, where irradiance will have been low. In the field, plants produced more storage roots. This was associated with higher AGPase activity and lower sucrose in the roots, indicating that feedforward loops enhanced sink capacity in the high light and low nitrogen environment in the field. Overall, these results indicated that carbon assimilation rate, the K battery, root starch synthesis, trehalose, and chlorogenic acid accumulation are potential target traits for genetic improvement.

The 'Tommy Atkins' mango genome reveals candidate genes for fruit quality.

BACKGROUND: Mango, Mangifera indica L., an important tropical fruit crop, is grown for its sweet and aromatic fruits. Past improvement of this species has predominantly relied on chance seedlings derived from over 1000 cultivars in the Indian sub-continent with a large variation for fruit size, yield, biotic and abiotic stress resistance, and fruit quality among other traits. Historically, mango has been an orphan crop with very limited molecular information. Only recently have molecular and genomics-based analyses enabled the creation of linkage maps, transcriptomes, and diversity analysis of large collections. Additionally, the combined analysis of genomic and phenotypic information is poised to improve mango breeding efficiency. RESULTS: This study sequenced, de novo assembled, analyzed, and annotated the genome of the monoembryonic mango cultivar 'Tommy Atkins'. The draft genome sequence was generated using NRGene de-novo Magic on high molecular weight DNA of 'Tommy Atkins', supplemented by 10X Genomics long read sequencing to improve the initial assembly. A hybrid population between 'Tommy Atkins' x 'Kensington Pride' was used to generate phased haplotype chromosomes and a highly resolved phased SNP map. The final 'Tommy Atkins' genome assembly was a consensus sequence that included 20 pseudomolecules representing the 20 chromosomes of mango and included ~ 86% of the ~ 439 Mb haploid mango genome. Skim sequencing identified ~ 3.3 M SNPs using the 'Tommy Atkins' x 'Kensington Pride' mapping population. Repeat masking identified 26,616 genes with a median length of 3348 bp. A whole genome duplication analysis revealed an ancestral 65 MYA polyploidization event shared with Anacardium occidentale. Two regions, one on LG4 and one on LG7 containing 28 candidate genes, were associated with the commercially important fruit size characteristic in the mapping population. CONCLUSIONS: The availability of the complete 'Tommy Atkins' mango genome will aid global initiatives to study mango genetics.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication.

KEY MESSAGE: Brazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation. Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300-4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity [Formula: see text], effective population size estimate [Formula: see text]) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava's center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

Draft Genome Sequence Resource of the Citrus Stem-End Rot Fungal Pathogen Lasiodiplodia theobromae CITRA15.

Lasiodiplodia theobromae is a fungal pathogen associated with perennial tropical fruit plants worldwide. In citrus, L. theobromae causes stem-end rot (Diplodia stem-end rot), a damaging postharvest disease that is aggravated when trees are also infected with the citrus greening bacteria 'Candidatus Liberibacter asiaticus'. Due to the latent infection of L. theobromae during the preharvest stage, it becomes difficult to control the disease by chemical or physical treatment. In the current study, we sequenced and assembled strain CITRA15, the first genome of L. theobromae obtained from diseased Citrus paradise 'Flame' grapefruit in Florida, and thereby provided a genomic resource for future research on diagnostics, and postharvest and preharvest disease management of citrus and other fruit crops.

Longitudinal Transcriptomic, Proteomic, and Metabolomic Analysis of Citrus limon Response to Graft Inoculation by Candidatus Liberibacter asiaticus.

Presymptomatic detection of citrus trees infected with Candidatus Liberibacter asiaticus (CLas), the bacterial pathogen associated with Huanglongbing (HLB; citrus greening disease), is critical to controlling the spread of the disease. To test whether infected citrus trees produce systemic signals that may be used for indirect disease detection, lemon (Citrus limon) plants were graft-inoculated with either CLas-infected or control (CLas-) budwood, and leaf samples were longitudinally collected over 46 weeks and analyzed for plant changes associated with CLas infection. RNA, protein, and metabolite samples extracted from leaves were analyzed using RNA-Seq, mass spectrometry, and 1H NMR spectroscopy, respectively. Significant differences in specific transcripts, proteins, and metabolites were observed between CLas-infected and control plants as early as 2 weeks post graft (wpg). The most dramatic differences between the transcriptome and proteome of CLas-infected and control plants were observed at 10 wpg, including coordinated increases in transcripts and proteins of citrus orthologs of known plant defense genes. This integrated approach to quantifying plant molecular changes in leaves of CLas-infected plants supports the development of diagnostic technology for presymptomatic or early disease detection as part of efforts to control the spread of HLB into uninfected citrus groves.

Longitudinal Transcriptomic, Proteomic, and Metabolomic Analyses of Citrus sinensis (L.) Osbeck Graft-Inoculated with "Candidatus Liberibacter asiaticus".

"Candidatus Liberibacter asiaticus" (CLas) is the bacterium associated with the citrus disease Huanglongbing (HLB). Current CLas detection methods are unreliable during presymptomatic infection, and understanding CLas pathogenicity to help develop new detection techniques is challenging because CLas has yet to be isolated in pure culture. To understand how CLas affects citrus metabolism and whether infected plants produce systemic signals that can be used to develop improved detection techniques, leaves from Washington Navel orange (Citrus sinensis (L.) Osbeck) plants were graft-inoculated with CLas and longitudinally studied using transcriptomics (RNA sequencing), proteomics (liquid chromatography-tandem mass spectrometry), and metabolomics (proton nuclear magnetic resonance). Photosynthesis gene expression and protein levels were lower in infected plants compared to controls during late infection, and lower levels of photosynthesis proteins were identified as early as 8 weeks post-grafting. These changes coordinated with higher sugar concentrations, which have been shown to accumulate during HLB. Cell wall modification and degradation gene expression and proteins were higher in infected plants during late infection. Changes in gene expression and proteins related to plant defense were observed in infected plants as early as 8 weeks post-grafting. These results reveal coordinated changes in greenhouse navel leaves during CLas infection at the transcript, protein, and metabolite levels, which can inform of biomarkers of early infection.

Molecular insight into cotton leaf curl geminivirus disease resistance in cultivated cotton (Gossypium hirsutum).

Cultivated cotton (Gossypium hirsutum) is the most important fibre crop in the world. Cotton leaf curl disease (CLCuD) is the major limiting factor and a threat to textile industry in India and Pakistan. All the local cotton cultivars exhibit moderate to no resistance against CLCuD. In this study, we evaluated an exotic cotton accession Mac7 as a resistance source to CLCuD by challenging it with viruliferous whiteflies and performing qPCR to evaluate the presence/absence and relative titre of CLCuD-associated geminiviruses/betasatellites. The results indicated that replication of pathogenicity determinant betasatellite is significantly attenuated in Mac7 and probably responsible for resistance phenotype. Afterwards, to decipher the genetic basis of CLCuD resistance in Mac7, we performed RNA sequencing on CLCuD-infested Mac7 and validated RNA-Seq data with qPCR on 24 independent genes. We performed co-expression network and pathway analysis for regulation of geminivirus/betasatellite-interacting genes. We identified nine novel modules with 52 hubs of highly connected genes in network topology within the co-expression network. Analysis of these hubs indicated the differential regulation of auxin stimulus and cellular localization pathways in response to CLCuD. We also analysed the differential regulation of geminivirus/betasatellite-interacting genes in Mac7. We further performed the functional validation of selected candidate genes via virus-induced gene silencing (VIGS). Finally, we evaluated the genomic context of resistance responsive genes and found that these genes are not specific to A or D sub-genomes of G. hirsutum. These results have important implications in understanding CLCuD resistance mechanism and developing a durable resistance in cultivated cotton.

BrAPI-an application programming interface for plant breeding applications.

MOTIVATION: Modern genomic breeding methods rely heavily on very large amounts of phenotyping and genotyping data, presenting new challenges in effective data management and integration. Recently, the size and complexity of datasets have increased significantly, with the result that data are often stored on multiple systems. As analyses of interest increasingly require aggregation of datasets from diverse sources, data exchange between disparate systems becomes a challenge. RESULTS: To facilitate interoperability among breeding applications, we present the public plant Breeding Application Programming Interface (BrAPI). BrAPI is a standardized web service API specification. The development of BrAPI is a collaborative, community-based initiative involving a growing global community of over a hundred participants representing several dozen institutions and companies. Development of such a standard is recognized as critical to a number of important large breeding system initiatives as a foundational technology. The focus of the first version of the API is on providing services for connecting systems and retrieving basic breeding data including germplasm, study, observation, and marker data. A number of BrAPI-enabled applications, termed BrAPPs, have been written, that take advantage of the emerging support of BrAPI by many databases. AVAILABILITY AND IMPLEMENTATION: More information on BrAPI, including links to the specification, test suites, BrAPPs, and sample implementations is available at https://brapi.org/. The BrAPI specification and the developer tools are provided as free and open source.

A quick guide for student-driven community genome annotation.

High quality gene models are necessary to expand the molecular and genetic tools available for a target organism, but these are available for only a handful of model organisms that have undergone extensive curation and experimental validation over the course of many years. The majority of gene models present in biological databases today have been identified in draft genome assemblies using automated annotation pipelines that are frequently based on orthologs from distantly related model organisms and usually have minor or major errors. Manual curation is time consuming and often requires substantial expertise, but is instrumental in improving gene model structure and identification. Manual annotation may seem to be a daunting and cost-prohibitive task for small research communities but involving undergraduates in community genome annotation consortiums can be mutually beneficial for both education and improved genomic resources. We outline a workflow for efficient manual annotation driven by a team of primarily undergraduate annotators. This model can be scaled to large teams and includes quality control processes through incremental evaluation. Moreover, it gives students an opportunity to increase their understanding of genome biology and to participate in scientific research in collaboration with peers and senior researchers at multiple institutions.

De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus.

The Cucurbita genus (squashes, pumpkins and gourds) includes important domesticated species such as C. pepo, C. maxima and C. moschata. In this study, we present a high-quality draft of the zucchini (C. pepo) genome. The assembly has a size of 263 Mb, a scaffold N50 of 1.8 Mb and 34 240 gene models. It includes 92% of the conserved BUSCO core gene set, and it is estimated to cover 93.0% of the genome. The genome is organized in 20 pseudomolecules that represent 81.4% of the assembly, and it is integrated with a genetic map of 7718 SNPs. Despite the small genome size, three independent lines of evidence support that the C. pepo genome is the result of a whole-genome duplication: the topology of the gene family phylogenies, the karyotype organization and the distribution of 4DTv distances. Additionally, 40 transcriptomes of 12 species of the genus were assembled and analysed together with all the other published genomes of the Cucurbitaceae family. The duplication was detected in all the Cucurbita species analysed, including C. maxima and C. moschata, but not in the more distant cucurbits belonging to the Cucumis and Citrullus genera, and it is likely to have occurred 30 ± 4 Mya in the ancestral species that gave rise to the genus.

The Tomato Expression Atlas.

SUMMARY: With the development of new high-throughput DNA sequencing technologies and decreasing costs, large gene expression datasets are being generated at an accelerating rate, but can be complex to visualize. New, more interactive and intuitive tools are needed to visualize the spatiotemporal context of expression data and help elucidate gene function. Using tomato fruit as a model, we have developed the Tomato Expression Atlas to facilitate effective data analysis, allowing the simultaneous visualization of groups of genes at a cell/tissue level of resolution within an organ, enhancing hypothesis development and testing in addition to candidate gene identification. This atlas can be adapted to different types of expression data from diverse multicellular species. AVAILABILITY AND IMPLEMENTATION: The Tomato Expression Atlas is available at http://tea.solgenomics.net/ . Source code is available at https://github.com/solgenomics/Tea . CONTACT: jr286@cornell.edu or lam87@cornell.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.