Pollinator-assisted plant phenotyping, selection, and breeding for crop resilience to abiotic stresses

Food security is threatened by climate change, with heat and drought being the main stresses affecting crop physiology and ecosystem services, such as plant-pollinator interactions. We hypothesize that tracking and ranking pollinators' preferences for flowers under environmental pressure could be used as a marker of plant quality for agricultural breeding to increase crop stress tolerance. Despite increasing relevance of flowers as the most stress sensitive organs, phenotyping platforms aim at identifying traits of resilience by assessing the plant physiological status through remote sensing-assisted vegetative indexes, but find strong bottlenecks in quantifying flower traits and in accurate genotype-to-phenotype prediction. However, as the transport of photoassimilates from leaves (sources) to flowers (sinks) is reduced in low-resilient plants, flowers are better indicators than leaves of plant well-being. Indeed, the chemical composition and amount of pollen and nectar that flowers produce, which ultimately serve as food resources for pollinators, change in response to environmental cues. Therefore, pollinators' preferences could be used as a measure of functional source-to-sink relationships for breeding decisions. To achieve this challenging goal, we propose to develop a pollinator-assisted phenotyping and selection platform for automated quantification of Genotype × Environment × Pollinator interactions through an insect geo-positioning system. Pollinator-assisted selection can be validated by metabolic, transcriptomic, and ionomic traits, and mapping of candidate genes, linking floral and leaf traits, pollinator preferences, plant resilience, and crop productivity. This radical new approach can change the current paradigm of plant phenotyping and find new paths for crop redomestication and breeding assisted by ecological decisions.

SELF PRUNING 3C is a flowering repressor that modulates seed germination, root architecture, and drought responses

Allelic variation in the CETS (CENTRORADIALIS, TERMINAL FLOWER 1, SELF PRUNING) gene family controls agronomically important traits in many crops. CETS genes encode phosphatidylethanolamine-binding proteins that have a central role in the timing of flowering as florigenic and anti-florigenic signals. The great expansion of CETS genes in many species suggests that the functions of this family go beyond flowering induction and repression. Here, we characterized the tomato SELF PRUNING 3C (SP3C) gene, and show that besides acting as a flowering repressor it also regulates seed germination and modulates root architecture. We show that loss of SP3C function in CRISPR/Cas9-generated mutant lines increases root length and reduces root side branching relative to the wild type. Higher SP3C expression in transgenic lines promotes the opposite effects in roots, represses seed germination, and also improves tolerance to water stress in seedlings. These discoveries provide new insights into the role of SP paralogs in agronomically relevant traits, and support future exploration of the involvement of CETS genes in abiotic stress responses.

Interaction between photoperiod and variation in circadian rhythms in tomato

BACKGROUND

Many biological processes follow circadian rhythmicity and are controlled by the circadian clock. Predictable environmental changes such as seasonal variation in photoperiod can modulate circadian rhythms, allowing organisms to adjust the timing of their biological processes to the time of the year. In some crops such as rice, barley or soybean, mutations in circadian clock genes have altered photoperiod sensitivity, enhancing their cultivability in specific seasons and latitudes. However, how changes in circadian rhythms interact with the perception of photoperiod in crops remain poorly studied. In tomato, the appearance during domestication of mutations in EMPFINDLICHER IM DUNKELROTEN LICHT 1 (EID1, Solyc09g075080) and NIGHT LIGHT-INDUCIBLE AND CLOCK-REGULATED GENE 2 (LNK2, Solyc01g068560) delayed both the phase and period of its circadian rhythms. The fact that variation in period and phase are separated in tomato provides an optimal tool to study how these factors affect the perception of photoperiod.

RESULTS

Here we develop tomato near isogenic lines carrying combinations of wild alleles of EID1 and LNK2 and show that they recreate the changes in phase and period that occurred during its domestication. We perform transcriptomic profiling of these near isogenic lines under two different photoperiods, and observe that EID1, but not LNK2, has a large effect on how the tomato transcriptome responds to photoperiod. This large effect of EID1 is likely a consequence of the global phase shift elicited by this gene in tomato's circadian rhythms.

CONCLUSIONS

Our study shows that changes in phase that occurred during tomato domestication determine photoperiod perception in this species, while changes in period have little effect.

The VIL gene CRAWLING ELEPHANT controls maturation and differentiation in tomato via polycomb silencing

VERNALIZATION INSENSITIVE 3-LIKE (VIL) proteins are PHD-finger proteins that recruit the repressor complex Polycomb Repressive Complex 2 (PRC2) to the promoters of target genes. Most known VIL targets are flowering repressor genes. Here, we show that the tomato VIL gene CRAWLING ELEPHANT (CREL) promotes differentiation throughout plant development by facilitating the trimethylation of Histone H3 on lysine 27 (H3K27me3). We identified the crel mutant in a screen for suppressors of the simple-leaf phenotype of entire (e), a mutant in the AUX/IAA gene ENTIRE/SlIAA9, involved in compound-leaf development in tomato. crel mutants have increased leaf complexity, and suppress the ectopic blade growth of e mutants. In addition, crel mutants are late flowering, and have delayed and aberrant stem, root and flower development. Consistent with a role for CREL in recruiting PRC2, crel mutants show drastically reduced H3K27me3 enrichment at approximately half of the 14,789 sites enriched in wild-type plants, along with upregulation of many underlying genes. Interestingly, this reduction in H3K27me3 across the genome in crel is also associated with gains in H3K27me3 at a smaller number of sites that normally have modest levels of the mark in wild-type plants, suggesting that PRC2 activity is no longer limiting in the absence of CREL. Our results uncover a wide role for CREL in plant and organ differentiation in tomato and suggest that CREL is required for targeting PRC2 activity to, and thus silencing, a specific subset of polycomb targets.

Plants form organs continuously throughout their lives, and the number and shape of their organs is determined in a flexible manner according to the internal and external circumstances. Alongside this flexibility, plants maintain basic developmental programs to ensure proper functioning. Among the ways by which plants achieve flexible development is by tuning the pace of their maturation and differentiation, at both the plant and organ levels. One of the ways plants regulate the rate of maturation and differentiation is by changing gene expression. Here, we identified a gene that promotes plant and organ maturation and differentiation. This gene, CRAWLING ELEPHANT (CREL) acts by bringing a repressing complex to target genes. We show the importance of CREL in multiple developmental processes and in the expression of multiple genes throughout the tomato genome.

Effective Mapping by Sequencing to Isolate Causal Mutations in the Tomato Genome

Forward genetic analysis remains as one of the most powerful tools for assessing gene functions, although the identification of the causal mutation responsible for a given phenotype has been a tedious and time-consuming task until recently. Advances in deep sequencing technologies have provided new approaches for the exploitation of natural and artificially induced genetic diversity, thus accelerating the discovery of novel allelic variants. In this chapter, a mapping-by-sequencing forward genetics approach is described to identify causal mutations in tomato (Solanum lycopersicum L.), a major crop species that is also a model species for plant biology and breeding.

Functional analysis of FRIGIDA using naturally occurring variation in Arabidopsis thaliana

The FRIGIDA locus (FRI, AT4G00650) has been extensively studied in Arabidopsis thaliana because of its role creating flowering time diversity. The FRI protein regulates flowering induction by binding partner proteins on its N-terminus and C-terminus domains and creating a supercomplex that promotes transcription of the floral repressor FLOWERING LOCUS C (FLC). Despite the knowledge accumulated on FRIGIDA (FRI), the function of the highly conserved central domain of the protein is still unknown. Functional characterization of naturally occurring DNA polymorphisms can provide useful information about the role of a protein and the localization of its operative domains. For FRI, loss-of-function mutations are positively selected and widespread in nature, making them a powerful tool to study the function of the different domains of the protein. Here we explore natural sequence variation in the FRI locus in more than 1000 Arabidopsis accessions. We identify 127 mutations that alter the FRI protein, including 60 that had never been described before. We defined 103 different alleles of FRI and study their association with variation in flowering time. We confirmed these associations by cloning 22 different alleles and expressing them in a common null genetic background. Our analysis pinpoints two single amino acid changes in the central domain that render the protein non-functional. We show that these two mutations determine the stability and cellular localization of the FRI protein. In summary, our work makes use of natural variants at the FRI locus to help understanding the function of the central domain of the FRI protein.

Mild drought in the vegetative stage induces phenotypic, gene expression, and DNA methylation plasticity in Arabidopsis but no transgenerational effects

There is renewed interest in whether environmentally induced changes in phenotypes can be heritable. In plants, heritable trait variation can occur without DNA sequence mutations through epigenetic mechanisms involving DNA methylation. However, it remains unknown whether this alternative system of inheritance responds to environmental changes and if it can provide a rapid way for plants to generate adaptive heritable phenotypic variation. To assess potential transgenerational effects induced by the environment, we subjected four natural accessions of Arabidopsis thaliana together with the reference accession Col-0 to mild drought in a multi-generational experiment. As expected, plastic responses to drought were observed in each accession, as well as a number of intergenerational effects of the parental environments. However, after an intervening generation without stress, except for a very few trait-based parental effects, descendants of stressed and non-stressed plants were phenotypically indistinguishable irrespective of whether they were grown in control conditions or under water deficit. In addition, genome-wide analysis of DNA methylation and gene expression in Col-0 demonstrated that, while mild drought induced changes in the DNA methylome of exposed plants, these variants were not inherited. We conclude that mild drought stress does not induce transgenerational epigenetic effects.

Mild drought in the vegetative stage induces phenotypic, gene expression, and DNA methylation plasticity in Arabidopsis but no transgenerational effects

There is renewed interest in whether environmentally induced changes in phenotypes can be heritable. In plants, heritable trait variation can occur without DNA sequence mutations through epigenetic mechanisms involving DNA methylation. However, it remains unknown whether this alternative system of inheritance responds to environmental changes and if it can provide a rapid way for plants to generate adaptive heritable phenotypic variation. To assess potential transgenerational effects induced by the environment, we subjected four natural accessions of Arabidopsis thaliana together with the reference accession Col-0 to mild drought in a multi-generational experiment. As expected, plastic responses to drought were observed in each accession, as well as a number of intergenerational effects of the parental environments. However, after an intervening generation without stress, except for a very few trait-based parental effects, descendants of stressed and non-stressed plants were phenotypically indistinguishable irrespective of whether they were grown in control conditions or under water deficit. In addition, genome-wide analysis of DNA methylation and gene expression in Col-0 demonstrated that, while mild drought induced changes in the DNA methylome of exposed plants, these variants were not inherited. We conclude that mild drought stress does not induce transgenerational epigenetic effects.

SlHAK20: a new player in plant salt tolerance

The majority of crops remain sensitive to salt stress despite the steady increase in salt concentration in agricultural soil. In this issue of The EMBO Journal, Wang et al (2020) screen hundreds of tomato accessions to identify SlHAK20 as a gene accounting for quantitative differences in salt tolerance between accessions. SlHAK20 is a potassium transporter belonging to a poorly studied clade in a large family of transporters, and its mutation induces salt susceptibility both in tomato and rice.

Natural variation in HsfA2 pre-mRNA splicing is associated with changes in thermotolerance during tomato domestication

Wild relatives of crops thrive in habitats where environmental conditions can be restrictive for productivity and survival of cultivated species. The genetic basis of this variability, particularly for tolerance to high temperatures, is not well understood. We examined the capacity of wild and cultivated accessions to acclimate to rapid temperature elevations that cause heat stress (HS). We investigated genotypic variation in thermotolerance of seedlings of wild and cultivated accessions. The contribution of polymorphisms associated with thermotolerance variation was examined regarding alterations in function of the identified gene. We show that tomato germplasm underwent a progressive loss of acclimation to strong temperature elevations. Sensitivity is associated with intronic polymorphisms in the HS transcription factor HsfA2 which affect the splicing efficiency of its pre-mRNA. Intron splicing in wild species results in increased synthesis of isoform HsfA2-II, implicated in the early stress response, at the expense of HsfA2-I which is involved in establishing short-term acclimation and thermotolerance. We propose that the selection for modern HsfA2 haplotypes reduced the ability of cultivated tomatoes to rapidly acclimate to temperature elevations, but enhanced their short-term acclimation capacity. Hence, we provide evidence that alternative splicing has a central role in the definition of plant fitness plasticity to stressful conditions.

Author Correction: Duplication of a domestication locus neutralized a cryptic variant that caused a breeding barrier in tomato

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Duplication of a domestication locus neutralized a cryptic variant that caused a breeding barrier in tomato

Genome editing technologies are being widely adopted in plant breeding¹. However, a looming challenge of engineering desirable genetic variation in diverse genotypes is poor predictability of phenotypic outcomes due to unforeseen interactions with pre-existing cryptic mutations^2-4. In tomato, breeding with a classical MADS-box gene mutation that improves harvesting by eliminating fruit stem abscission frequently results in excessive inflorescence branching, flowering and reduced fertility due to interaction with a cryptic variant that causes partial mis-splicing in a homologous gene^5-8. Here, we show that a recently evolved tandem duplication carrying the second-site variant achieves a threshold of functional transcripts to suppress branching, enabling breeders to neutralize negative epistasis on yield. By dissecting the dosage mechanisms by which this structural variant restored normal flowering and fertility, we devised strategies that use CRISPR-Cas9 genome editing to predictably improve harvesting. Our findings highlight the under-appreciated impact of epistasis in targeted trait breeding and underscore the need for a deeper characterization of cryptic variation to enable the full potential of genome editing in agriculture.

Arabidopsis SME1 Regulates Plant Development and Response to Abiotic Stress by Determining Spliceosome Activity Specificity

The control of precursor-messenger RNA (pre-mRNA) splicing is emerging as an important layer of regulation in plant responses to endogenous and external cues. In eukaryotes, pre-mRNA splicing is governed by the activity of a large ribonucleoprotein machinery, the spliceosome, whose protein core is composed of the Sm ring and the related Sm-like 2-8 complex. Recently, the Arabidopsis (Arabidopsis thaliana) Sm-like 2-8 complex has been characterized. However, the role of plant Sm proteins in pre-mRNA splicing remains largely unknown. Here, we present the functional characterization of Sm protein E1 (SME1), an Arabidopsis homolog of the SME subunit of the eukaryotic Sm ring. Our results demonstrate that SME1 regulates the spliceosome activity and that this regulation is controlled by the environmental conditions. Indeed, depending on the conditions, SME1 ensures the efficiency of constitutive and alternative splicing of selected pre-mRNAs. Moreover, missplicing of most targeted pre-mRNAs leads to the generation of nonsense-mediated decay signatures, indicating that SME1 also guarantees adequate levels of the corresponding functional transcripts. In addition, we show that the selective function of SME1 in ensuring appropriate gene expression patterns through the regulation of specific pre-mRNA splicing is essential for adequate plant development and adaptation to freezing temperatures. These findings reveal that SME1 plays a critical role in plant development and interaction with the environment by providing spliceosome activity specificity.

Pivotal Roles of Cryptochromes 1a and 2 in Tomato Development and Physiology

Cryptochromes are flavin-containing blue/UVA light photoreceptors that regulate various plant light-induced physiological processes. In Arabidopsis (Arabidopsis thaliana), cryptochromes mediate de-etiolation, photoperiodic control of flowering, entrainment of the circadian clock, cotyledon opening and expansion, anthocyanin accumulation, and root growth. In tomato (Solanum lycopersicum), cryptochromes are encoded by a multigene family, comprising CRY1a, CRY1b, CRY2, and CRY3 We have previously reported the phenotypes of tomato cry1a mutants and CRY2 overexpressing plants. Here, we report the isolation by targeting induced local lesions in genomes, of a tomato cry2 knock-out mutant, its introgression in the indeterminate Moneymaker background, and the phenotypes of cry1a/cry2 single and double mutants. The cry1a/cry2 mutant showed phenotypes similar to its Arabidopsis counterpart (long hypocotyls in white and blue light), but also several additional features such as increased seed weight and internode length, enhanced hypocotyl length in red light, inhibited primary root growth under different light conditions, anticipation of flowering under long-day conditions, and alteration of the phase of circadian leaf movements. Both cry1a and cry2 control the levels of photosynthetic pigments in leaves, but cry2 has a predominant role in fruit pigmentation. Metabolites of the sterol, tocopherol, quinone, and sugar classes are differentially accumulated in cry1a and cry2 leaves and fruits. These results demonstrate a pivotal role of cryptochromes in controlling tomato development and physiology. The manipulation of these photoreceptors represents a powerful tool to influence important agronomic traits such as flowering time and fruit quality.

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.

NPR1 mediates a novel regulatory pathway in cold acclimation by interacting with HSFA1 factors

NON-EXPRESSER OF PATHOGENESIS-RELATED GENES 1 (NPR1) is a master regulator of plant response to pathogens that confers immunity through a transcriptional cascade mediated by salicylic acid and TGA transcription factors. Little is known, however, about its implication in plant response to abiotic stress. Here, we provide genetic and molecular evidence supporting the fact that Arabidopsis NPR1 plays an essential role in cold acclimation by regulating cold-induced gene expression independently of salicylic acid and TGA factors. Our results demonstrate that, in response to low temperature, cytoplasmic NPR1 oligomers release monomers that translocate to the nucleus where they interact with heat shock transcription factor 1 (HSFA1) to promote the induction of HSFA1-regulated genes and cold acclimation. These findings unveil an unexpected function for NPR1 in plant response to low temperature, reveal a new regulatory pathway for cold acclimation mediated by NPR1 and HSFA1 factors, and place NPR1 as a central hub integrating cold and pathogen signalling for a better adaptation of plants to an ever-changing environment.

Prominent alterations of wild barley leaf transcriptome in response to individual and combined drought acclimation and heat shock conditions

Under field conditions, drought and heat stress typically happen simultaneously and their negative impact on the agricultural production is expected to increase worldwide under the climate change scenario. In this study, we performed RNA-sequencing analysis on leaves of wild barley (Hordeum spontaneum) originated from the northern coastal region of Egypt following individual drought acclimation (DA) and heat shock (HS) treatments and their combination (CS, combined stresses) to distinguish the unique and shared differentially expressed genes (DEG). Results indicated that the number of unique genes that were differentially expressed following HS treatment exceeded the number of those expressed following DA. In addition, the number of genes that were uniquely differentially expressed in response to CS treatment exceeded the number of those of shared responses to individual DA and HS treatments. These results indicate a better adaptation of the Mediterranean wild barley to drought conditions when compared with heat stress. It also manifests that the wild barley response to CS tends to be unique rather than common. Annotation of DEG showed that metabolic processes were the most influenced biological function in response to the applied stresses.

Profiling the onset of somatic embryogenesis in Arabidopsis

BACKGROUND

Totipotency is the ability of a cell to regenerate a whole organism. Plant somatic embryogenesis (SE) is a remarkable example of totipotency because somatic cells reverse differentiation, respond to an appropriate stimulus and initiate embryo development. Although SE is an ideal system to investigate de-differentiation and differentiation, we still lack a deep molecular understanding of the phenomenon due to experimental restraints.

RESULTS

We applied the INTACT method to specifically isolate the nuclei of those cells undergoing SE among the majority of non-embryogenic cells that make up a callus. We compared the transcriptome of embryogenic cells to the one of proliferating callus cells. Our analyses revealed that embryogenic cells are transcriptionally rather than metabolically active. Embryogenic cells shut off biochemical pathways involved in carbohydrate and lipid metabolism and activate the transcriptional machinery. Furthermore, we show how early in SE, ground tissue and leaf primordia specification are switched on before the specification of a shoot apical meristem.

CONCLUSIONS

This is the first attempt to specifically profile embryogenic cells among the different cell types that constitute plant in vitro tissue cultures. Our comparative analyses provide insights in the gene networks regulating SE and open new research avenues in the field of plant regeneration.

Multi-Omics of Tomato Glandular Trichomes Reveals Distinct Features of Central Carbon Metabolism Supporting High Productivity of Specialized Metabolites

Glandular trichomes are metabolic cell factories with the capacity to produce large quantities of secondary metabolites. Little is known about the connection between central carbon metabolism and metabolic productivity for secondary metabolites in glandular trichomes. To address this gap in our knowledge, we performed comparative metabolomics, transcriptomics, proteomics, and ¹³C-labeling of type VI glandular trichomes and leaves from a cultivated (Solanum lycopersicum LA4024) and a wild (Solanum habrochaites LA1777) tomato accession. Specific features of glandular trichomes that drive the formation of secondary metabolites could be identified. Tomato type VI trichomes are photosynthetic but acquire their carbon essentially from leaf sucrose. The energy and reducing power from photosynthesis are used to support the biosynthesis of secondary metabolites, while the comparatively reduced Calvin-Benson-Bassham cycle activity may be involved in recycling metabolic CO₂ Glandular trichomes cope with oxidative stress by producing high levels of polyunsaturated fatty acids, oxylipins, and glutathione. Finally, distinct mechanisms are present in glandular trichomes to increase the supply of precursors for the isoprenoid pathways. Particularly, the citrate-malate shuttle supplies cytosolic acetyl-CoA and plastidic glycolysis and malic enzyme support the formation of plastidic pyruvate. A model is proposed on how glandular trichomes achieve high metabolic productivity.

New Strategies and Tools in Quantitative Genetics: How to Go from the Phenotype to the Genotype

Quantitative genetics has a long history in plants: It has been used to study specific biological processes, identify the factors important for trait evolution, and breed new crop varieties. These classical approaches to quantitative trait locus mapping have naturally improved with technology. In this review, we show how quantitative genetics has evolved recently in plants and how new developments in phenotyping, population generation, sequencing, gene manipulation, and statistics are rejuvenating both the classical linkage mapping approaches (for example, through nested association mapping) as well as the more recently developed genome-wide association studies. These strategies are complementary in most instances, and indeed, one is often used to confirm the results of the other. Despite significant advances, an emerging trend is that the outcome and efficiency of the different approaches depend greatly on the genetic architecture of the trait in the genetic material under study.

New Strategies and Tools in Quantitative Genetics: How to Go from the Phenotype to the Genotype

Quantitative genetics has a long history in plants: It has been used to study specific biological processes, identify the factors important for trait evolution, and breed new crop varieties. These classical approaches to quantitative trait locus mapping have naturally improved with technology. In this review, we show how quantitative genetics has evolved recently in plants and how new developments in phenotyping, population generation, sequencing, gene manipulation, and statistics are rejuvenating both the classical linkage mapping approaches (for example, through nested association mapping) as well as the more recently developed genome-wide association studies. These strategies are complementary in most instances, and indeed, one is often used to confirm the results of the other. Despite significant advances, an emerging trend is that the outcome and efficiency of the different approaches depend greatly on the genetic architecture of the trait in the genetic material under study.

A footprint of desiccation tolerance in the genome of Xerophyta viscosa

Desiccation tolerance is common in seeds and various other organisms, but only a few angiosperm species possess vegetative desiccation tolerance. These 'resurrection species' may serve as ideal models for the ultimate design of crops with enhanced drought tolerance. To understand the molecular and genetic mechanisms enabling vegetative desiccation tolerance, we produced a high-quality whole-genome sequence for the resurrection plant Xerophyta viscosa and assessed transcriptome changes during its dehydration. Data revealed induction of transcripts typically associated with desiccation tolerance in seeds and involvement of orthologues of ABI3 and ABI5, both key regulators of seed maturation. Dehydration resulted in both increased, but predominantly reduced, transcript abundance of genomic 'clusters of desiccation-associated genes' (CoDAGs), reflecting the cessation of growth that allows for the expression of desiccation tolerance. Vegetative desiccation tolerance in X. viscosa was found to be uncoupled from drought-induced senescence. We provide strong support for the hypothesis that vegetative desiccation tolerance arose by redirection of genetic information from desiccation-tolerant seeds.

Alternative splicing enhances transcriptome complexity in desiccating seeds

Before being dispersed in the environment, mature seeds need to be dehydrated. The survival of seeds after dispersal depends on their low hydration in combination with high desiccation tolerance. These characteristics are established during seed maturation. Some key seed maturation genes have been reported to be regulated by alternative splicing (AS). However, so far AS was described only for single genes and a comprehensive analysis of AS during seed maturation has been lacking. We investigated gene expression and AS during Arabidopsis thaliana seed development at a global level, before and after desiccation. Bioinformatics tools were developed to identify differentially spliced regions within genes. Our data suggest the importance and shows the peculiar features of AS during seed desiccation. We identified AS in 34% of genes that are expressed at both timepoints before and after desiccation. Most of these AS transcript variants had not been found before in other tissues. Among the AS genes some seed master regulators could be found. Interestingly, 6% of all expressed transcripts were not transcriptionally regulated during desiccation, but only modified by AS. We propose that AS should be more routinely taken into account in the analysis of transcriptomic data to prevent overlooking potentially important regulators.

CLAUSA Is a MYB Transcription Factor That Promotes Leaf Differentiation by Attenuating Cytokinin Signaling

Leaf morphogenesis and differentiation are highly flexible processes, resulting in a large diversity of leaf forms. The development of compound leaves involves an extended morphogenesis stage compared with that of simple leaves, and the tomato (Solanum lycopersicum) mutant clausa (clau) exposes a potential for extended morphogenesis in tomato leaves. Here, we report that the CLAU gene encodes a MYB transcription factor that has evolved a unique role in compound-leaf species to promote an exit from the morphogenetic phase of tomato leaf development. We show that CLAU attenuates cytokinin signaling, and that clau plants have increased cytokinin sensitivity. The results suggest that flexible leaf patterning involves a coordinated interplay between transcription factors and hormones.

Targeted and Untargeted Approaches Unravel Novel Candidate Genes and Diagnostic SNPs for Quantitative Resistance of the Potato (Solanum tuberosum L.) to Phytophthora infestans Causing the Late Blight Disease

The oomycete Phytophthora infestans causes late blight of potato, which can completely destroy the crop. Therefore, for the past 160 years, late blight has been the most important potato disease worldwide. The identification of cultivars with high and durable field resistance to P. infestans is an objective of most potato breeding programs. This type of resistance is polygenic and therefore quantitative. Its evaluation requires multi-year and location trials. Furthermore, quantitative resistance to late blight correlates with late plant maturity, a negative agricultural trait. Knowledge of the molecular genetic basis of quantitative resistance to late blight not compromised by late maturity is very limited. It is however essential for developing diagnostic DNA markers that facilitate the efficient combination of superior resistance alleles in improved cultivars. We used association genetics in a population of 184 tetraploid potato cultivars in order to identify single nucleotide polymorphisms (SNPs) that are associated with maturity corrected resistance (MCR) to late blight. The population was genotyped for almost 9000 SNPs from three different sources. The first source was candidate genes specifically selected for their function in the jasmonate pathway. The second source was novel candidate genes selected based on comparative transcript profiling (RNA-Seq) of groups of genotypes with contrasting levels of quantitative resistance to P. infestans. The third source was the first generation 8.3k SolCAP SNP genotyping array available in potato for genome wide association studies (GWAS). Twenty seven SNPs from all three sources showed robust association with MCR. Some of those were located in genes that are strong candidates for directly controlling quantitative resistance, based on functional annotation. Most important were: a lipoxygenase (jasmonate pathway), a 3-hydroxy-3-methylglutaryl coenzyme A reductase (mevalonate pathway), a P450 protein (terpene biosynthesis), a transcription factor and a homolog of a major gene for resistance to P. infestans from the wild potato species Solanum venturii. The candidate gene approach and GWAS complemented each other as they identified different genes. The results of this study provide new insight in the molecular genetic basis of quantitative resistance in potato and a toolbox of diagnostic SNP markers for breeding applications.

Analysis of Circadian Leaf Movements

The circadian clock is a molecular timekeeper that controls a wide variety of biological processes. In plants, clock outputs range from the molecular level, with rhythmic gene expression and metabolite content, to physiological processes such as stomatal conductance or leaf movements. Any of these outputs can be used as markers to monitor the state of the circadian clock. In the model plant Arabidopsis thaliana, much of the current knowledge about the clock has been gained from time course experiments profiling expression of endogenous genes or reporter constructs regulated by the circadian clock. Since these methods require labor-intensive sample preparation or transformation, monitoring leaf movements is an interesting alternative, especially in non-model species and for natural variation studies. Technological improvements both in digital photography and image analysis allow cheap and easy monitoring of circadian leaf movements. In this chapter we present a protocol that uses an autonomous point and shoot camera and free software to monitor circadian leaf movements in tomato.

Nonsense-mediated mRNA decay modulates immune receptor levels to regulate plant antibacterial defense

Nonsense-mediated mRNA decay (NMD) is a conserved eukaryotic RNA surveillance mechanism that degrades aberrant mRNAs. NMD impairment in Arabidopsis is linked to constitutive immune response activation and enhanced antibacterial resistance, but the underlying mechanisms are unknown. Here we show that NMD contributes to innate immunity in Arabidopsis by controlling the turnover of numerous TIR domain-containing, nucleotide-binding, leucine-rich repeat (TNL) immune receptor-encoding mRNAs. Autoimmunity resulting from NMD impairment depends on TNL signaling pathway components and can be triggered through deregulation of a single TNL gene, RPS6. Bacterial infection of plants causes host-programmed inhibition of NMD, leading to stabilization of NMD-regulated TNL transcripts. Conversely, constitutive NMD activity prevents TNL stabilization and impairs plant defense, demonstrating that host-regulated NMD contributes to disease resistance. Thus, NMD shapes plant innate immunity by controlling the threshold for activation of TNL resistance pathways.

The relationship between flowering time and growth responses to drought in the Arabidopsis Landsberg erecta x Antwerp-1 population

Limited water availability is one of the most prominent abiotic constraints to plant survival and reproduction. Thus, plants have evolved different strategies to cope with water deficit, including modification of their growth and timing of developmental events such as flowering. In this work, we explore the link between flowering time and growth responses to moderate drought stress in Arabidopsis thaliana using natural variation for these traits found in the Landsberg erecta x Antwerp-1 recombinant inbred line population. We developed and phenotyped near isogenic lines containing different allelic combinations at three interacting quantitative trait loci (QTL) affecting both flowering time and growth in response to water deficit. We used these lines to confirm additive and epistatic effects of the three QTL and observed a strong association between late flowering and reduced sensitivity to drought. Analyses of growth responses to drought over time revealed that late flowering plants were able to recover their growth in the second half of their vegetative development. In contrast, early flowering, a common drought escape strategy that ensures plant survival under severe water deficit, was associated with strongly impaired plant fitness. The results presented here indicate that late flowering may be advantageous under continuous mild water deficit as it allows stress acclimatization over time.

Functional analysis of the Landsberg erecta allele of FRIGIDA

BACKGROUND

Most of the natural variation in flowering time in Arabidopsis thaliana can be attributed to allelic variation at the gene FRIGIDA (FRI, AT4G00650), which activates expression of the floral repressor FLOWERING LOCUS C (FLC, AT5G10140). Usually, late-flowering accessions carry functional FRI alleles (FRI-wt), whereas early flowering accessions contain non-functional alleles. The two most frequent alleles found in early flowering accessions are the ones present in the commonly used lab strains Columbia (FRI-Col) and Landsberg erecta (FRI-Ler), which contain a premature stop codon and a deletion of the start codon respectively.

RESULTS

Analysis of flowering time data from various Arabidopsis natural accessions indicated that the FRI-Ler allele retains some functionality. We generated transgenic lines carrying the FRI-Col or FRI-Ler allele in order to compare their effect on flowering time, vernalization response and FLC expression in the same genetic background. We characterize their modes of regulation through allele-specific expression and their relevance in nature through re-analysis of published datasets. We demonstrate that the FRI-Ler allele induces FLC expression, delays flowering time and confers sensitivity to vernalization in contrast to the true null FRI-Col allele. Nevertheless, the FRI-Ler allele revealed a weaker effect when compared to the fully functional FRI-wt allele, mainly due to reduced expression.

CONCLUSIONS

The present study defines for the first time the existence of a new class of Arabidopsis accessions with an intermediate phenotype between slow and rapid cycling types. Although using available data from a common garden experiment we cannot observe fitness differences between accessions carrying the FRI-Ler or the FRI-Col allele, the phenotypic changes observed in the lab suggest that variation in these alleles could play a role in adaptation to specific natural environments.

Network types and their application in natural variation studies in plants

We are in the age of data-driven biology. Not even a decade after the invention of high-throughput sequencing technologies, there are methods that accurately monitor DNA polymorphisms, transcription profiles, methylation states, transcription factor binding sites, chromatin compactness, nucleosome positions, dynamic histone marks, and so on. We are starting to generate comparable amounts of protein or metabolite data. A key issue is how are we going to make sense of all this information. Network analysis is the most promising method to integrate, query and display large amounts of data for human interpretation. This review shortly summarizes the basic types of networks, their properties and limitations. In addition, I introduce the application of networks to the study of the molecular mechanisms behind natural phenotypic variation.

Comparative transcriptomics reveals patterns of selection in domesticated and wild tomato

Although applied over extremely short timescales, artificial selection has dramatically altered the form, physiology, and life history of cultivated plants. We have used RNAseq to define both gene sequence and expression divergence between cultivated tomato and five related wild species. Based on sequence differences, we detect footprints of positive selection in over 50 genes. We also document thousands of shifts in gene-expression level, many of which resulted from changes in selection pressure. These rapidly evolving genes are commonly associated with environmental response and stress tolerance. The importance of environmental inputs during evolution of gene expression is further highlighted by large-scale alteration of the light response coexpression network between wild and cultivated accessions. Human manipulation of the genome has heavily impacted the tomato transcriptome through directed admixture and by indirectly favoring nonsynonymous over synonymous substitutions. Taken together, our results shed light on the pervasive effects artificial and natural selection have had on the transcriptomes of tomato and its wild relatives.

A quantitative genetic basis for leaf morphology in a set of precisely defined tomato introgression lines

Introgression lines (ILs), in which genetic material from wild tomato species is introgressed into a domesticated background, have been used extensively in tomato (Solanum lycopersicum) improvement. Here, we genotype an IL population derived from the wild desert tomato Solanum pennellii at ultrahigh density, providing the exact gene content harbored by each line. To take advantage of this information, we determine IL phenotypes for a suite of vegetative traits, ranging from leaf complexity, shape, and size to cellular traits, such as stomatal density and epidermal cell phenotypes. Elliptical Fourier descriptors on leaflet outlines provide a global analysis of highly heritable, intricate aspects of leaf morphology. We also demonstrate constraints between leaflet size and leaf complexity, pavement cell size, and stomatal density and show independent segregation of traits previously assumed to be genetically coregulated. Meta-analysis of previously measured traits in the ILs shows an unexpected relationship between leaf morphology and fruit sugar levels, which RNA-Seq data suggest may be attributable to genetically coregulated changes in fruit morphology or the impact of leaf shape on photosynthesis. Together, our results both improve upon the utility of an important genetic resource and attest to a complex, genetic basis for differences in leaf morphology between natural populations.

Fine genetic mapping of RXopJ4, a bacterial spot disease resistance locus from Solanum pennellii LA716

The RXopJ4 resistance locus from the wild accession Solanum pennellii (Sp) LA716 confers resistance to bacterial spot disease of tomato (S. lycopersicum, Sl) caused by Xanthomonas perforans (Xp). RXopJ4 resistance depends on recognition of the pathogen type III effector protein XopJ4. We used a collection of Sp introgression lines (ILs) to narrow the RXopJ4 locus to a 4.2-Mb segment on the long arm of chromosome 6, encompassed by the ILs 6-2 and 6-2-2. We then adapted or developed a collection of 14 molecular markers to map on a segregating F(2) population from a cross between the susceptible parent Sl FL8000 and the resistant parent RXopJ4 8000 OC(7). In the F(2) population, a 190-kb segment between the markers J350 and J352 cosegregated with resistance. This fine mapping will enable both the identification of candidate genes and the detection of resistant plants using cosegregating markers. The RXopJ4 resistance gene(s), in combination with other recently characterized genes and a quantitative trait locus (QTL) for bacterial spot disease resistance, will likely be an effective tool for the development of durable resistance in cultivated tomato.

Native environment modulates leaf size and response to simulated foliar shade across wild tomato species

The laminae of leaves optimize photosynthetic rates by serving as a platform for both light capture and gas exchange, while minimizing water losses associated with thermoregulation and transpiration. Many have speculated that plants maximize photosynthetic output and minimize associated costs through leaf size, complexity, and shape, but a unifying theory linking the plethora of observed leaf forms with the environment remains elusive. Additionally, the leaf itself is a plastic structure, responsive to its surroundings, further complicating the relationship. Despite extensive knowledge of the genetic mechanisms underlying angiosperm leaf development, little is known about how phenotypic plasticity and selective pressures converge to create the diversity of leaf shapes and sizes across lineages. Here, we use wild tomato accessions, collected from locales with diverse levels of foliar shade, temperature, and precipitation, as a model to assay the extent of shade avoidance in leaf traits and the degree to which these leaf traits correlate with environmental factors. We find that leaf size is correlated with measures of foliar shade across the wild tomato species sampled and that leaf size and serration correlate in a species-dependent fashion with temperature and precipitation. We use far-red induced changes in leaf length as a proxy measure of the shade avoidance response, and find that shade avoidance in leaves negatively correlates with the level of foliar shade recorded at the point of origin of an accession. The direction and magnitude of these correlations varies across the leaf series, suggesting that heterochronic and/or ontogenic programs are a mechanism by which selective pressures can alter leaf size and form. This study highlights the value of wild tomato accessions for studies of both morphological and light-regulated development of compound leaves, and promises to be useful in the future identification of genes regulating potentially adaptive plastic leaf traits.

Next generation quantitative genetics in plants

Most characteristics in living organisms show continuous variation, which suggests that they are controlled by multiple genes. Quantitative trait loci (QTL) analysis can identify the genes underlying continuous traits by establishing associations between genetic markers and observed phenotypic variation in a segregating population. The new high-throughput sequencing (HTS) technologies greatly facilitate QTL analysis by providing genetic markers at genome-wide resolution in any species without previous knowledge of its genome. In addition HTS serves to quantify molecular phenotypes, which aids to identify the loci responsible for QTLs and to understand the mechanisms underlying diversity. The constant improvements in price, experimental protocols, computational pipelines, and statistical frameworks are making feasible the use of HTS for any research group interested in quantitative genetics. In this review I discuss the application of HTS for molecular marker discovery, population genotyping, and expression profiling in QTL analysis.

Network analysis identifies ELF3 as a QTL for the shade avoidance response in Arabidopsis

Quantitative Trait Loci (QTL) analyses in immortal populations are a powerful method for exploring the genetic mechanisms that control interactions of organisms with their environment. However, QTL analyses frequently do not culminate in the identification of a causal gene due to the large chromosomal regions often underlying QTLs. A reasonable approach to inform the process of causal gene identification is to incorporate additional genome-wide information, which is becoming increasingly accessible. In this work, we perform QTL analysis of the shade avoidance response in the Bayreuth-0 (Bay-0, CS954) x Shahdara (Sha, CS929) recombinant inbred line population of Arabidopsis. We take advantage of the complex pleiotropic nature of this trait to perform network analysis using co-expression, eQTL and functional classification from publicly available datasets to help us find good candidate genes for our strongest QTL, SAR2. This novel network analysis detected EARLY FLOWERING 3 (ELF3; AT2G25930) as the most likely candidate gene affecting the shade avoidance response in our population. Further genetic and transgenic experiments confirmed ELF3 as the causative gene for SAR2. The Bay-0 and Sha alleles of ELF3 differentially regulate developmental time and circadian clock period length in Arabidopsis, and the extent of this regulation is dependent on the light environment. This is the first time that ELF3 has been implicated in the shade avoidance response and that different natural alleles of this gene are shown to have phenotypic effects. In summary, we show that development of networks to inform candidate gene identification for QTLs is a promising technique that can significantly accelerate the process of QTL cloning.

Sequence diversity in three tomato species: SNPs, markers, and molecular evolution

BACKGROUND

Tomato species are of significant agricultural and ecological interest, with cultivated tomato being among the most common vegetable crops grown. Wild tomato species are native to diverse habitats in South America and show great morphological and ecological diversity that has proven useful in breeding programs. However, relatively little is known about nucleotide diversity between tomato species. Until recently limited sequence information was available for tomato, preventing genome-wide evolutionary analyses. Now, an extensive collection of tomato expressed sequence tags (ESTs) is available at the SOL Genomics Network (SGN). This database holds sequences from several species, annotated with quality values, assembled into unigenes, and tested for homology against other genomes. Despite the importance of polymorphism detection for breeding and natural variation studies, such analyses in tomato have mostly been restricted to cultivated accessions. Importantly, previous polymorphisms surveys mostly ignored the linked meta-information, limiting functional and evolutionary analyses. The current data in SGN is thus an under-exploited resource. Here we describe a cross-species analysis taking full-advantage of available information.

RESULTS

We mined 20,000 interspecific polymorphisms between Solanum lycopersicum and S. habrochaites or S. pennellii and 28,800 intraspecific polymorphisms within S. lycopersicum. Using the available meta-information we classified genes into functional categories and obtained estimations of single nucleotide polymorphisms (SNP) quality, position in the gene, and effect on the encoded proteins, allowing us to perform evolutionary analyses. Finally, we developed a set of more than 10,000 between-species molecular markers optimized by sequence quality and predicted intron position. Experimental validation of 491 of these molecular markers resulted in confirmation of 413 polymorphisms.

CONCLUSION

We present a new analysis of the extensive tomato EST sequences available that represents the most comprehensive survey of sequence diversity across Solanum species to date. These SNPs, plus thousands of molecular makers designed to detect the polymorphisms are available to the community via a website. Evolutionary analyses on these polymorphism uncovered sets of genes potentially important for the evolution and domestication of tomato; interestingly these sets were enriched for genes involved in response to the environment.

Quantitative genetic analysis of flowering time in tomato

Artificial selection of cultivated tomato (Solanum lycopersicum L.) has resulted in the generation of early-flowering, day-length-insensitive cultivars, despite its close relationship to other Solanum species that need more time and specific photoperiods to flower. To investigate the genetic mechanisms controlling flowering time in tomato and related species, we performed a quantitative trait locus (QTL) analysis for flowering time in an F2 mapping population derived from S. lycopersicum and its late-flowering wild relative S. chmielewskii. Flowering time was scored as the number of days from sowing to the opening of the first flower (days to flowering), and as the number of leaves under the first inflorescence (leaf number). QTL analyses detected 2 QTLs affecting days to flowering, which explained 55.3% of the total phenotypic variance, and 6 QTLs for leaf number, accounting for 66.7% of the corresponding phenotypic variance. Four of the leaf number QTLs had not previously been detected for this trait in tomato. Colocation of some QTLs with flowering-time genes included in the genetic map suggests PHYB2, FALSIFLORA, and a tomato FLC-like sequence as candidate genes that might have been targets of selection during the domestication of tomato.