Idiosyncratic and dose-dependent epistasis drives variation in tomato fruit size.

Epistasis between genes is traditionally studied with mutations that eliminate protein activity, but most natural genetic variation is in cis-regulatory DNA and influences gene expression and function quantitatively. In this study, we used natural and engineered cis-regulatory alleles in a plant stem-cell circuit to systematically evaluate epistatic relationships controlling tomato fruit size. Combining a promoter allelic series with two other loci, we collected over 30,000 phenotypic data points from 46 genotypes to quantify how allele strength transforms epistasis. We revealed a saturating dose-dependent relationship but also allele-specific idiosyncratic interactions, including between alleles driving a step change in fruit size during domestication. Our approach and findings expose an underexplored dimension of epistasis, in which cis-regulatory allelic diversity within gene regulatory networks elicits nonlinear, unpredictable interactions that shape phenotypes.

Establishing Physalis as a Solanaceae model system enables genetic reevaluation of the inflated calyx syndrome.

The highly diverse Solanaceae family contains several widely studied models and crop species. Fully exploring, appreciating, and exploiting this diversity requires additional model systems. Particularly promising are orphan fruit crops in the genus Physalis, which occupy a key evolutionary position in the Solanaceae and capture understudied variation in traits such as inflorescence complexity, fruit ripening and metabolites, disease and insect resistance, self-compatibility, and most notable, the striking inflated calyx syndrome (ICS), an evolutionary novelty found across angiosperms where sepals grow exceptionally large to encapsulate fruits in a protective husk. We recently developed transformation and genome editing in Physalis grisea (groundcherry). However, to systematically explore and unlock the potential of this and related Physalis as genetic systems, high-quality genome assemblies are needed. Here, we present chromosome-scale references for P. grisea and its close relative Physalis pruinosa and use these resources to study natural and engineered variations in floral traits. We first rapidly identified a natural structural variant in a bHLH gene that causes petal color variation. Further, and against expectations, we found that CRISPR-Cas9-targeted mutagenesis of 11 MADS-box genes, including purported essential regulators of ICS, had no effect on inflation. In a forward genetics screen, we identified huskless, which lacks ICS due to mutation of an AP2-like gene that causes sepals and petals to merge into a single whorl of mixed identity. These resources and findings elevate Physalis to a new Solanaceae model system and establish a paradigm in the search for factors driving ICS.

A conserved superlocus regulates above- and belowground root initiation.

Plants continuously form new organs in different developmental contexts in response to environmental cues. Underground lateral roots initiate from prepatterned cells in the main root, but cells can also bypass the root-shoot trajectory separation and generate shoot-borne roots through an unknown mechanism. We mapped tomato (Solanum lycopersicum) shoot-borne root development at single-cell resolution and showed that these roots initiate from phloem-associated cells through a unique transition state. This state requires the activity of a transcription factor that we named SHOOTBORNE ROOTLESS (SBRL). Evolutionary analysis reveals that SBRL's function and cis regulation are conserved in angiosperms and that it arose as an ancient duplication, with paralogs controlling wound-induced and lateral root initiation. We propose that the activation of a common transition state by context-specific regulators underlies the plasticity of plant root systems.

Dynamic evolution of small signalling peptide compensation in plant stem cell control.

Gene duplications are a hallmark of plant genome evolution and a foundation for genetic interactions that shape phenotypic diversity1-5. Compensation is a major form of paralogue interaction6-8 but how compensation relationships change as allelic variation accumulates is unknown. Here we leveraged genomics and genome editing across the Solanaceae family to capture the evolution of compensating paralogues. Mutations in the stem cell regulator CLV3 cause floral organs to overproliferate in many plants9-11. In tomato, this phenotype is partially suppressed by transcriptional upregulation of a closely related paralogue12. Tobacco lost this paralogue, resulting in no compensation and extreme clv3 phenotypes. Strikingly, the paralogues of petunia and groundcherry nearly completely suppress clv3, indicating a potent ancestral state of compensation. Cross-species transgenic complementation analyses show that this potent compensation partially degenerated in tomato due to a single amino acid change in the paralogue and cis-regulatory variation that limits its transcriptional upregulation. Our findings show how genetic interactions are remodelled following duplications and suggest that dynamic paralogue evolution is widespread over short time scales and impacts phenotypic variation from natural and engineered mutations.

Dissecting cis-regulatory control of quantitative trait variation in a plant stem cell circuit.

Cis-regulatory mutations underlie important crop domestication and improvement traits1,2. However, limited allelic diversity has hindered functional dissection of the large number of cis-regulatory elements and their potential interactions, thereby precluding a deeper understanding of how cis-regulatory variation impacts traits quantitatively. Here, we engineered over 60 promoter alleles in two tomato fruit size genes3,4 to characterize cis-regulatory sequences and study their functional relationships. We found that targeted mutations in conserved promoter sequences of SlCLV3, a repressor of stem cell proliferation5,6, have a weak impact on fruit locule number. Pairwise combinations of these mutations mildly enhance this phenotype, revealing additive and synergistic relationships between conserved regions and further suggesting even higher-order cis-regulatory interactions within the SlCLV3 promoter. In contrast, SlWUS, a positive regulator of stem cell proliferation repressed by SlCLV3 (refs. 5,6), is more tolerant to promoter perturbations. Our results show that complex interplay among cis-regulatory variants can shape quantitative variation, and suggest that empirical dissections of this hidden complexity can guide promoter engineering to predictably modify crop traits.

Conserved pleiotropy of an ancient plant homeobox gene uncovered by cis-regulatory dissection.

Divergence of gene function is a hallmark of evolution, but assessing functional divergence over deep time is not trivial. The few alleles available for cross-species studies often fail to expose the entire functional spectrum of genes, potentially obscuring deeply conserved pleiotropic roles. Here, we explore the functional divergence of WUSCHEL HOMEOBOX9 (WOX9), suggested to have species-specific roles in embryo and inflorescence development. Using a cis-regulatory editing drive system, we generate a comprehensive allelic series in tomato, which revealed hidden pleiotropic roles for WOX9. Analysis of accessible chromatin and conserved cis-regulatory sequences identifies the regions responsible for this pleiotropic activity, the functions of which are conserved in groundcherry, a tomato relative. Mimicking these alleles in Arabidopsis, distantly related to tomato and groundcherry, reveals new inflorescence phenotypes, exposing a deeply conserved pleiotropy. We suggest that targeted cis-regulatory mutations can uncover conserved gene functions and reduce undesirable effects in crop improvement.

Tomato HAIRY MERISTEM genes are involved in meristem maintenance and compound leaf morphogenesis.

The HAIRY MERISTEM (HAM) genes function in meristem maintenance but play minor roles in the morphogenesis of a simple leaf that is determinate. Here, we functionally analyzed HAM genes in tomato and uncovered their involvement in compound leaf morphogenesis. Tomato encodes three HAM homologs, of which SlHAM and SlHAM2 (SlHAMs) are guided for cleavage by microRNA171 and are abundant in the shoot and floral meristems as well as in the compound leaf primordia. We found that SlHAMs silencing led to overproliferation of cells in the periphery of the meristems where SlHAM is localized. As in meristems, leaf-specific silencing of SlHAMs provoked overproliferation of meristematic cells in the organogenic compound leaf rachis. We further demonstrate that the meristematic cell overproliferation in both meristems and leaves was in part due to the misexpression of the stem cell regulator WUSCHEL, previously shown to be induced by cytokinin. Strikingly, reduction of cytokinin levels in SlHAMs-silenced leaves completely suppressed the overproliferation phenotype, suggesting a regulatory link between SlHAMs and cytokinin, a key hormone found to promote indeterminacy in meristems and leaves. Taken together, our data provide evidence that in addition to their conserved function in meristem maintenance, SlHAMs are also required for the proper morphogenesis of the compound leaf.

The tomato NAC transcription factor SlNAM2 is involved in flower-boundary morphogenesis.

Being composed of several whorls of distinct floral organs, the flower is one of the most complex organs in the plant. As such, the formation and maintenance of boundaries that separate the meristem from the floral organ primordium and adjacent organs are critical for its normal development. In Arabidopsis, the miR164-regulated NAM genes play key roles in floral-boundary specification. By contrast, much less is known about floral-boundary establishment in the model crop tomato. It was found that the miR164-regulated NAM gene GOBLET is expressed in the floral meristem-organ boundaries and its loss-of-function mutant produces flowers with fused organs, indicating its requirement for tomato floral-boundary formation. It was found here that sly-miR164 targets the transcripts of three additional uncharacterized NAM genes in developing flowers. It is shown that, after floral-boundary initiation, the NAM gene Solyc03g115850 (SlNAM2) is expressed as stripes that mark the boundaries between sepals and between different floral whorls. Furthermore, ectopic accumulation of SlNAM2-encoding transcripts caused various growth-suppression and extraorgan phenotypes typically observed in plants over-expressing known boundary genes. Flower-specific silencing of sly-miR164-targeted NAM genes (AP1>>MIR164) caused defects in the separation of sepals and floral whorls indicating abnormal boundary specification. However, supplementing these NAM-deficient flowers with miR164-resistant SlNAM2 suppressed their fusion phenotypes and completely restored floral boundaries. Together, our results strongly suggest that SlNAM2 participates in the establishment of tomato flower whorl and sepal boundaries.

The developmental outcomes of P0-mediated ARGONAUTE destabilization in tomato.

The plant protein ARGONAUTE1 (AGO1) functions in multiple RNA-silencing pathways, including those of microRNAs, key regulators of growth and development. Genetic analysis of ago1 mutants with informative defects has provided valuable insights into AGO1's biological functions. Tomato encodes two AGO1 homologs (SlAGO1s), but mutants have not been described to date. To analyze SlAGO1s' involvement in development, we confirmed that both undergo decay in the presence of the Polerovirus silencing suppressor P0 and produce a transgenic responder line (OP:P0HA) that, upon transactivation, expresses P0 C-terminally fused to a hemagglutinin (HA) tag (P0HA) and destabilizes SlAGO1s at the site of expression. By crossing OP:P0HA with a battery of driver lines, constitutive as well as organ- and stage-specific SlAGO1 downregulation was induced in the F1 progeny. Activated plants exhibited various developmental phenotypes that partially overlapped with those of Arabidopsis ago1 mutants. Plants that constitutively expressed P0HA had reduced SlAGO1 levels and increased accumulation of miRNA targets, indicating compromised SlAGO1-mediated silencing. Consistent with this, they exhibited pleiotropic morphological defects and their growth was arrested post-germination. Transactivation of P0HA in young leaf and floral organ primordia dramatically modified corresponding organ morphology, including the radialization of leaflets, petals and anthers, suggesting that SlAGO1s' activities are required for normal lateral organ development and polarity. Overall, our results suggest that the OP:P0HA responder line can serve as a valuable tool to suppress SlAGO1 silencing pathways in tomato. The suppression of additional SlAGOs by P0HA and its contribution to the observed phenotypes awaits investigation.

Identification and characterization of a novel miR159 target not related to MYB in tomato.

MicroRNA 159 (miR159) is a highly conserved miRNA with roles in flowering under short days, anther development and seed germination via repression of GAMYB-like genes. In tomato, the function of miR159 (Sl-miR159) is currently unknown and target transcripts have not been experimentally validated. Here, we identified and characterized a new miR159 target gene (SGN-U567133) in Solanum lycopersicum (tomato) that is not related to MYB. SGN-U567133 is predominantly expressed in flowers and encodes a nuclear-localized protein that contains a unique NOZZLE-like domain at its N terminus. In tomato, SGN-U567133 represents a small gene family and orthologs have been identified in other plant species, all containing a conserved miR159 target site in their coding sequence. Accordingly, 5'-RACE cleavage assay supported miRNA-mediated cleavage of SGN-U567133 transcripts in vivo. Moreover, the SGN-U567133 transcript accumulated in P19-HA-expressing tomato leaves in which miRNA-mediated cleavage is inhibited. In addition, transgenic tomato plants expressing a miR159-resistant form of SGN-U567133 accumulated higher levels of the SGN-U567133 transcript and exhibited defects in leaf and flower development. Together, our results suggest that SGN-U567133 represents a novel class of miR159 targets in plants and raise the possibility that its post-transcriptional regulation by Sl-miR159 is essential for normal tomato development.

Transgenic expression of tomato bushy stunt virus silencing suppressor P19 via the pOp/LhG4 transactivation system induces viral-like symptoms in tomato.

During natural infection, the Tomato bushy stunt virus (TBSV) silencing suppressor protein P19 is expressed at high levels, which are required for optimum viral pathogenicity and silencing suppression. To date, expression of P19 in transgenic host plants has failed to achieve comparable expression levels and thus has provided only limited information on its in planta effects. To obtain high P19 expression and study its effects on host plant development in the absence of virus infection, we generated HA-tagged P19 (P19HA)-transgenic tomato reporter plants using the pOp/LhG4 transactivation system, which separates transformation from transgene expression. Upon reporter plant activation with a strong constitutive promoter, the transactivated F1 plants expressed high levels of a functional P19HA protein and displayed multiple abnormal phenotypes, some of which were highly reminiscent of the symptoms described previously for TBSV-infected tomato. Moreover, phenotype severity correlated with P19HA expression level, amount of bound miRNA/miRNA* duplexes, and accumulation of miRNA target transcripts. Together our results demonstrate that the tomato miRNA pathway is markedly compromised by P19, in particular when this protein is relatively abundant, as occurs during natural infection. We suggest that such interference with endogenous silencing may be responsible for at least some of the symptoms characteristic of TBSV-infected tomato.