Major Impacts of Widespread Structural Variation on Gene Expression and Crop Improvement in Tomato.

Structural variants (SVs) underlie important crop improvement and domestication traits. However, resolving the extent, diversity, and quantitative impact of SVs has been challenging. We used long-read nanopore sequencing to capture 238,490 SVs in 100 diverse tomato lines. This panSV genome, along with 14 new reference assemblies, revealed large-scale intermixing of diverse genotypes, as well as thousands of SVs intersecting genes and cis-regulatory regions. Hundreds of SV-gene pairs exhibit subtle and significant expression changes, which could broadly influence quantitative trait variation. By combining quantitative genetics with genome editing, we show how multiple SVs that changed gene dosage and expression levels modified fruit flavor, size, and production. In the last example, higher order epistasis among four SVs affecting three related transcription factors allowed introduction of an important harvesting trait in modern tomato. Our findings highlight the underexplored role of SVs in genotype-to-phenotype relationships and their widespread importance and utility in crop improvement.

Bypassing Negative Epistasis on Yield in Tomato Imposed by a Domestication Gene.

Selection for inflorescence architecture with improved flower production and yield is common to many domesticated crops. However, tomato inflorescences resemble wild ancestors, and breeders avoided excessive branching because of low fertility. We found branched variants carry mutations in two related transcription factors that were selected independently. One founder mutation enlarged the leaf-like organs on fruits and was selected as fruit size increased during domestication. The other mutation eliminated the flower abscission zone, providing "jointless" fruit stems that reduced fruit dropping and facilitated mechanical harvesting. Stacking both beneficial traits caused undesirable branching and sterility due to epistasis, which breeders overcame with suppressors. However, this suppression restricted the opportunity for productivity gains from weak branching. Exploiting natural and engineered alleles for multiple family members, we achieved a continuum of inflorescence complexity that allowed breeding of higher-yielding hybrids. Characterizing and neutralizing similar cases of negative epistasis could improve productivity in many agricultural organisms. VIDEO ABSTRACT.

New Horizons for Dissecting Epistasis in Crop Quantitative Trait Variation.

Uncovering the genes, variants, and interactions underlying crop diversity is a frontier in plant genetics. Phenotypic variation often does not reflect the cumulative effect of individual gene mutations. This deviation is due to epistasis, in which interactions between alleles are often unpredictable and quantitative in effect. Recent advances in genomics and genome-editing technologies are elevating the study of epistasis in crops. Using the traits and developmental pathways that were major targets in domestication and breeding, we highlight how epistasis is central in guiding the behavior of the genetic variation that shapes quantitative trait variation. We outline new strategies that illuminate how quantitative epistasis from modified gene dosage defines background dependencies. Advancing our understanding of epistasis in crops can reveal new principles and approaches to engineering targeted improvements in agriculture.

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.

Optimized sample selection for cost-efficient long-read population sequencing.

An increasingly important scenario in population genetics is when a large cohort has been genotyped using a low-resolution approach (e.g., microarrays, exome capture, short-read WGS), from which a few individuals are resequenced using a more comprehensive approach, especially long-read sequencing. The subset of individuals selected should ensure that the captured genetic diversity is fully representative and includes variants across all subpopulations. For example, human variation has historically focused on individuals with European ancestry, but this represents a small fraction of the overall diversity. Addressing this, SVCollector identifies the optimal subset of individuals for resequencing by analyzing population-level VCF files from low-resolution genotyping studies. It then computes a ranked list of samples that maximizes the total number of variants present within a subset of a given size. To solve this optimization problem, SVCollector implements a fast, greedy heuristic and an exact algorithm using integer linear programming. We apply SVCollector on simulated data, 2504 human genomes from the 1000 Genomes Project, and 3024 genomes from the 3000 Rice Genomes Project and show the rankings it computes are more representative than alternative naive strategies. When selecting an optimal subset of 100 samples in these cohorts, SVCollector identifies individuals from every subpopulation, whereas naive methods yield an unbalanced selection. Finally, we show the number of variants present in cohorts selected using this approach follows a power-law distribution that is naturally related to the population genetic concept of the allele frequency spectrum, allowing us to estimate the diversity present with increasing numbers of samples.

The Starch Granule-Associated Protein EARLY STARVATION1 Is Required for the Control of Starch Degradation in Arabidopsis thaliana Leaves.

To uncover components of the mechanism that adjusts the rate of leaf starch degradation to the length of the night, we devised a screen for mutant Arabidopsis thaliana plants in which starch reserves are prematurely exhausted. The mutation in one such mutant, named early starvation1 (esv1), eliminates a previously uncharacterized protein. Starch in mutant leaves is degraded rapidly and in a nonlinear fashion, so that reserves are exhausted 2 h prior to dawn. The ESV1 protein and a similar uncharacterized Arabidopsis protein (named Like ESV1 [LESV]) are located in the chloroplast stroma and are also bound into starch granules. The region of highest similarity between the two proteins contains a series of near-repeated motifs rich in tryptophan. Both proteins are conserved throughout starch-synthesizing organisms, from angiosperms and monocots to green algae. Analysis of transgenic plants lacking or overexpressing ESV1 or LESV, and of double mutants lacking ESV1 and another protein necessary for starch degradation, leads us to propose that these proteins function in the organization of the starch granule matrix. We argue that their misexpression affects starch degradation indirectly, by altering matrix organization and, thus, accessibility of starch polymers to starch-degrading enzymes.

PROTEIN TARGETING TO STARCH is required for localising GRANULE-BOUND STARCH SYNTHASE to starch granules and for normal amylose synthesis in Arabidopsis.

The domestication of starch crops underpinned the development of human civilisation, yet we still do not fully understand how plants make starch. Starch is composed of glucose polymers that are branched (amylopectin) or linear (amylose). The amount of amylose strongly influences the physico-chemical behaviour of starchy foods during cooking and of starch mixtures in non-food manufacturing processes. The GRANULE-BOUND STARCH SYNTHASE (GBSS) is the glucosyltransferase specifically responsible for elongating amylose polymers and was the only protein known to be required for its biosynthesis. Here, we demonstrate that PROTEIN TARGETING TO STARCH (PTST) is also specifically required for amylose synthesis in Arabidopsis. PTST is a plastidial protein possessing an N-terminal coiled coil domain and a C-terminal carbohydrate binding module (CBM). We discovered that Arabidopsis ptst mutants synthesise amylose-free starch and are phenotypically similar to mutants lacking GBSS. Analysis of granule-bound proteins showed a dramatic reduction of GBSS protein in ptst mutant starch granules. Pull-down assays with recombinant proteins in vitro, as well as immunoprecipitation assays in planta, revealed that GBSS physically interacts with PTST via a coiled coil. Furthermore, we show that the CBM domain of PTST, which mediates its interaction with starch granules, is also required for correct GBSS localisation. Fluorescently tagged Arabidopsis GBSS, expressed either in tobacco or Arabidopsis leaves, required the presence of Arabidopsis PTST to localise to starch granules. Mutation of the CBM of PTST caused GBSS to remain in the plastid stroma. PTST fulfils a previously unknown function in targeting GBSS to starch. This sheds new light on the importance of targeting biosynthetic enzymes to sub-cellular sites where their action is required. Importantly, PTST represents a promising new gene target for the biotechnological modification of starch composition, as it is exclusively involved in amylose synthesis.

The Enzyme-Like Domain of Arabidopsis Nuclear ß-Amylases Is Critical for DNA Sequence Recognition and Transcriptional Activation.

Plant BZR1-BAM transcription factors contain a ß-amylase (BAM)-like domain, characteristic of proteins involved in starch breakdown. The enzyme-derived domains appear to be noncatalytic, but they determine the function of the two Arabidopsis thaliana BZR1-BAM isoforms (BAM7 and BAM8) during transcriptional initiation. Removal or swapping of the BAM domains demonstrates that the BAM7 BAM domain restricts DNA binding and transcriptional activation, while the BAM8 BAM domain allows both activities. Furthermore, we demonstrate that BAM7 and BAM8 interact on the protein level and cooperate during transcriptional regulation. Site-directed mutagenesis of residues in the BAM domain of BAM8 shows that its function as a transcriptional activator is independent of catalysis but requires an intact substrate binding site, suggesting it may bind a ligand. Microarray experiments with plants overexpressing truncated versions lacking the BAM domain indicate that the pseudo-enzymatic domain increases selectivity for the preferred cis-regulatory element BBRE (BZR1-BAM Responsive Element). Side specificity toward the G-box may allow crosstalk to other signaling networks. This work highlights the importance of the enzyme-derived domain of BZR1-BAMs, supporting their potential role as metabolic sensors.

ß-amylase-like proteins function as transcription factors in Arabidopsis, controlling shoot growth and development.

Plants contain ß-amylase-like proteins (BAMs; enzymes usually associated with starch breakdown) present in the nucleus rather than targeted to the chloroplast. They possess BRASSINAZOLE RESISTANT1 (BZR1)-type DNA binding domains--also found in transcription factors mediating brassinosteroid (BR) responses. The two Arabidopsis thaliana BZR1-BAM proteins (BAM7 and BAM8) bind a cis-regulatory element that both contains a G box and resembles a BR-responsive element. In protoplast transactivation assays, these BZR1-BAMs activate gene expression. Structural modeling suggests that the BAM domain's glucan binding cleft is intact, but the recombinant proteins are at least 1000 times less active than chloroplastic ß-amylases. Deregulation of BZR1-BAMs (the bam7bam8 double mutant and BAM8-overexpressing plants) causes altered leaf growth and development. Of the genes upregulated in plants overexpressing BAM8 and downregulated in bam7bam8 plants, many carry the cis-regulatory element in their promoters. Many genes that respond to BRs are inversely regulated by BZR1-BAMs. We propose a role for BZR1-BAMs in controlling plant growth and development through crosstalk with BR signaling. Furthermore, we speculate that BZR1-BAMs may transmit metabolic signals by binding a ligand in their BAM domain, although diurnal changes in the concentration of maltose, a candidate ligand produced by chloroplastic ß-amylases, do not influence their transcription factor function.

Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing.

Advancing crop genomics requires efficient genetic systems enabled by high-quality personalized genome assemblies. Here, we introduce RagTag, a toolset for automating assembly scaffolding and patching, and we establish chromosome-scale reference genomes for the widely used tomato genotype M82 along with Sweet-100, a new rapid-cycling genotype that we developed to accelerate functional genomics and genome editing in tomato. This work outlines strategies to rapidly expand genetic systems and genomic resources in other plant species.

Non-invasive topology analysis of membrane proteins in the secretory pathway.

We present a novel method to experimentally visualize in vivo the topology of transmembrane proteins residing in the endoplasmic reticulum (ER) membrane or passing through the secretory pathway on their way to their final destination. This approach, so-called redox-based topology analysis (ReTA), is based on fusion of transmembrane proteins with redox-sensitive GFP (roGFP) and ratiometric imaging. The ratio images provide direct information on the orientation of roGFP relative to the membrane as the roGFP fluorescence alters with changes in the glutathione redox potential across the ER membrane. As proof of concept, we produced binary read-outs using oxidized roGFP inside the ER lumen and reduced roGFP on the cytosolic side of the membrane for both N- and C-terminal fusions of single and multi-spanning membrane proteins. Further, successive deletion of hydrophobic domains from the C-terminus of the K/HDEL receptor ERD2 resulted in alternating localization of roGFP and a topology model for AtERD2 with six transmembrane domains.

RaGOO: fast and accurate reference-guided scaffolding of draft genomes.

We present RaGOO, a reference-guided contig ordering and orienting tool that leverages the speed and sensitivity of Minimap2 to accurately achieve chromosome-scale assemblies in minutes. After the pseudomolecules are constructed, RaGOO identifies structural variants, including those spanning sequencing gaps. We show that RaGOO accurately orders and orients 3 de novo tomato genome assemblies, including the widely used M82 reference cultivar. We then demonstrate the scalability and utility of RaGOO with a pan-genome analysis of 103 Arabidopsis thaliana accessions by examining the structural variants detected in the newly assembled pseudomolecules. RaGOO is available open source at https://github.com/malonge/RaGOO .

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 breeding1. 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 mutations2-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 gene5-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.

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.

Rapid improvement of domestication traits in an orphan crop by genome editing.

Genome editing holds great promise for increasing crop productivity, and there is particular interest in advancing breeding in orphan crops, which are often burdened by undesirable characteristics resembling wild relatives. We developed genomic resources and efficient transformation in the orphan Solanaceae crop 'groundcherry' (Physalis pruinosa) and used clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein-9 nuclease (Cas9) (CRISPR-Cas9) to mutate orthologues of tomato domestication and improvement genes that control plant architecture, flower production and fruit size, thereby improving these major productivity traits. Thus, translating knowledge from model crops enables rapid creation of targeted allelic diversity and novel breeding germplasm in distantly related orphan crops.

Variation in the flowering gene SELF PRUNING 5G promotes day-neutrality and early yield in tomato.

Plants evolved so that their flowering is triggered by seasonal changes in day length. However, day-length sensitivity in crops limits their geographical range of cultivation, and thus modification of the photoperiod response was critical for their domestication. Here we show that loss of day-length-sensitive flowering in tomato was driven by the florigen paralog and flowering repressor SELF-PRUNING 5G (SP5G). SP5G expression is induced to high levels during long days in wild species, but not in cultivated tomato because of cis-regulatory variation. CRISPR/Cas9-engineered mutations in SP5G cause rapid flowering and enhance the compact determinate growth habit of field tomatoes, resulting in a quick burst of flower production that translates to an early yield. Our findings suggest that pre-existing variation in SP5G facilitated the expansion of cultivated tomato beyond its origin near the equator in South America, and they provide a compelling demonstration of the power of gene editing to rapidly improve yield traits in crop breeding.