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.

Abstract A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis -regulatory sequences. It remains unclear how such drastic cis -regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis- regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ∼125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 ( CLV3 ) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis -regulatory sequences. Our results demonstrate remarkable malleability in cis -regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis -regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis -regulation to effectively engineer trait variation from conserved productivity genes in crops. Author summary We investigated the evolution of cis -regulatory elements (CREs) and their interactions in the regulation of a plant stem cell regulator gene, CLAVATA3 (CLV3) , in Arabidopsis and tomato. Despite diverging ∼125 million years ago, the function and expression of CLV3 is conserved in these species; however, cis -regulatory sequences upstream and downstream have drastically diverged, preventing identification of conserved non-coding sequences between them. We used CRISPR-Cas9 to engineer dozens of mutations within the cis -regulatory regions of Arabidopsis and tomato CLV3. In tomato, our results show that tomato CLV3 function primarily relies on interactions among CREs in the 5’ non-coding region, unlike Arabidopsis CLV3 , which depends on a more balanced distribution of functional CREs between the 5’ and 3’ regions. Therefore, despite a high degree of functional conservation, our study demonstrates divergent regulatory strategies between two distantly related CLV3 orthologs, with substantial alterations in regulatory sequences, their spatial arrangement, and their relative effects on CLV3 regulation. These results suggest that regulatory regions are not only extremely robust to mutagenesis, but also that the sequences underlying this robustness can be lineage-specific for conserved genes, due to the complex and often redundant interactions among CREs that ensure proper gene function amidst large-scale sequence turnover.

ABSTRACT Calcium (Ca 2+ )-dependent protein kinases (CDPKs or CPKs) are a unique family of Ca 2+ -sensor/kinase-effector proteins with diverse functions in plants. In Arabidopsis thaliana , CPK28 contributes to immune homeostasis by promoting degradation of the key immune signaling receptor-like cytoplasmic kinase BOTRYTIS-INDUCED KINASE 1 (BIK1), and additionally functions in vegetative-to-reproductive stage transition. How CPK28 controls these seemingly disparate pathways is unknown. Here, we identify a single phosphorylation site in the kinase domain of CPK28 (Ser318) that is differentially required for its function in immune homeostasis and stem elongation. We show that CPK28 undergoes intra- and inter-molecular auto-phosphorylation on Ser318 and can additionally be trans-phosphorylated on this residue by BIK1. Analysis of several other phosphorylation sites demonstrates that Ser318 phosphorylation is uniquely required to prime CPK28 for Ca 2+ activation at physiological concentrations of Ca 2+ , possibly through stabilization of the Ca 2+ -bound active state as indicated by intrinsic fluorescence experiments. Together, our data indicate that phosphorylation of Ser318 is required for the activation of CPK28 at low intracellular [Ca 2+ ] to prevent initiation of an immune response in the absence of infection. By comparison, phosphorylation of Ser318 is not required for stem-elongation, indicating pathway specific requirements for phosphorylation-based Ca 2+ -sensitivity priming. We additionally provide evidence for a conserved function for Ser318 phosphorylation in related group IV CDPKs which holds promise for biotechnological applications by generating CDPK alleles that enhance resistance to microbial pathogens without consequences to yield.