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.

Rapid customization of Solanaceae fruit crops for urban agriculture.

Cultivation of crops in urban environments might reduce the environmental impact of food production1-4. However, lack of available land in cities and a need for rapid crop cycling, to yield quickly and continuously, mean that so far only lettuce and related 'leafy green' vegetables are cultivated in urban farms5. New fruit varieties with architectures and yields suitable for urban farming have proven difficult to breed1,5. We identified a regulator of tomato stem length (SlER) and devised a trait-stacking strategy to combine mutations for condensed shoots, rapid flowering (SP5G) and precocious growth termination (SP). Application of our strategy using one-step CRISPR-Cas9 genome editing restructured vine-like tomato plants into compact, early yielding plants suitable for urban agriculture. Field data confirmed that yields were maintained, and we demonstrated cultivation in indoor farming systems. Targeting the same stem length regulator alone in groundcherry, another Solanaceae plant, also enabled engineering to a compact stature. Our approach can expand the repertoire of crops for urban 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.

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.

Evolution of buffering in a genetic circuit controlling plant stem cell proliferation.

Precise control of plant stem cell proliferation is necessary for the continuous and reproducible development of plant organs1,2. The peptide ligand CLAVATA3 (CLV3) and its receptor protein kinase CLAVATA1 (CLV1) maintain stem cell homeostasis within a deeply conserved negative feedback circuit1,2. In Arabidopsis, CLV1 paralogs also contribute to homeostasis, by compensating for the loss of CLV1 through transcriptional upregulation3. Here, we show that compensation4,5 operates in diverse lineages for both ligands and receptors, but while the core CLV signaling module is conserved, compensation mechanisms have diversified. Transcriptional compensation between ligand paralogs operates in tomato, facilitated by an ancient gene duplication that impacted the domestication of fruit size. In contrast, we found little evidence for transcriptional compensation between ligands in Arabidopsis and maize, and receptor compensation differs between tomato and Arabidopsis. Our findings show that compensation among ligand and receptor paralogs is critical for stem cell homeostasis, but that diverse genetic mechanisms buffer conserved developmental programs.

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.

Genome-wide Analysis of Transcriptional Variability in a Large Maize-Teosinte Population.

Gene expression regulation plays an important role in controlling plant phenotypes and adaptation. Here, we report a comprehensive assessment of gene expression variation through the transcriptome analyses of a large maize-teosinte experimental population. Genome-wide mapping identified 25 660 expression quantitative trait loci (eQTL) for 17 311 genes, capturing an unprecedented range of expression variation. We found that local eQTL were more frequently mapped to adjacent genes, displaying a mode of expression piggybacking, which consequently created co-regulated gene clusters. Genes within the co-regulated gene clusters tend to have relevant functions and shared chromatin modifications. Distant eQTL formed 125 significant distant eQTL hotspots with their targets significantly enriched in specific functional categories. By integrating different sources of information, we identified putative trans- regulators for a variety of metabolic pathways. We demonstrated that the bHLH transcription factor R1 and hexokinase HEX9 might act as crucial regulators for flavonoid biosynthesis and glycolysis, respectively. Moreover, we showed that domestication or improvement has significantly affected global gene expression, with many genes targeted by selection. Of particular interest, the Bx genes for benzoxazinoid biosynthesis may have undergone coordinated cis-regulatory divergence between maize and teosinte, and a transposon insertion that inactivates Bx12 was under strong selection as maize spread into temperate environments with a distinct herbivore community.

Engineering Quantitative Trait Variation for Crop Improvement by Genome Editing.

Major advances in crop yields are needed in the coming decades. However, plant breeding is currently limited by incremental improvements in quantitative traits that often rely on laborious selection of rare naturally occurring mutations in gene-regulatory regions. Here, we demonstrate that CRISPR/Cas9 genome editing of promoters generates diverse cis-regulatory alleles that provide beneficial quantitative variation for breeding. We devised a simple genetic scheme, which exploits trans-generational heritability of Cas9 activity in heterozygous loss-of-function mutant backgrounds, to rapidly evaluate the phenotypic impact of numerous promoter variants for genes regulating three major productivity traits in tomato: fruit size, inflorescence branching, and plant architecture. Our approach allows immediate selection and fixation of novel alleles in transgene-free plants and fine manipulation of yield components. Beyond a platform to enhance variation for diverse agricultural traits, our findings provide a foundation for dissecting complex relationships between gene-regulatory changes and control of quantitative traits.

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.

The evolution of inflorescence diversity in the nightshades and heterochrony during meristem maturation.

One of the most remarkable manifestations of plant evolution is the diversity for floral branching systems. These "inflorescences" arise from stem cell populations in shoot meristems that mature gradually to reproductive states in response to environmental and endogenous signals. The morphology of the shoot meristem maturation process is conserved across distantly related plants, raising the question of how diverse inflorescence architectures arise from seemingly common maturation programs. In tomato and related nightshades (Solanaceae), inflorescences range from solitary flowers to highly branched structures bearing hundreds of flowers. Since reproductive barriers between even closely related Solanaceae have precluded a genetic dissection, we captured and compared meristem maturation transcriptomes from five domesticated and wild species reflecting the evolutionary continuum of inflorescence complexity. We find these divergent species share hundreds of dynamically expressed genes, enriched for transcription factors. Meristem stages are defined by distinct molecular states and point to modified maturation schedules underlying architectural variation. These modified schedules are marked by a peak of transcriptome expression divergence during the reproductive transition, driven by heterochronic shifts of dynamic genes, including transcriptional regulators with known roles in flowering. Thus, evolutionary diversity in Solanaceae inflorescence complexity is determined by subtle modifications of transcriptional programs during a critical transitional window of meristem maturation, which we propose underlies similar cases of plant architectural variation. More broadly, our findings parallel the recently described transcriptome "inverse hourglass" model for animal embryogenesis, suggesting both plant and animal morphological variation is guided by a mid-development period of transcriptome divergence.

The role of cis regulatory evolution in maize domestication.

Gene expression differences between divergent lineages caused by modification of cis regulatory elements are thought to be important in evolution. We assayed genome-wide cis and trans regulatory differences between maize and its wild progenitor, teosinte, using deep RNA sequencing in F1 hybrid and parent inbred lines for three tissue types (ear, leaf and stem). Pervasive regulatory variation was observed with approximately 70% of ∼17,000 genes showing evidence of regulatory divergence between maize and teosinte. However, many fewer genes (1,079 genes) show consistent cis differences with all sampled maize and teosinte lines. For ∼70% of these 1,079 genes, the cis differences are specific to a single tissue. The number of genes with cis regulatory differences is greatest for ear tissue, which underwent a drastic transformation in form during domestication. As expected from the domestication bottleneck, maize possesses less cis regulatory variation than teosinte with this deficit greatest for genes showing maize-teosinte cis regulatory divergence, suggesting selection on cis regulatory differences during domestication. Consistent with selection on cis regulatory elements, genes with cis effects correlated strongly with genes under positive selection during maize domestication and improvement, while genes with trans regulatory effects did not. We observed a directional bias such that genes with cis differences showed higher expression of the maize allele more often than the teosinte allele, suggesting domestication favored up-regulation of gene expression. Finally, this work documents the cis and trans regulatory changes between maize and teosinte in over 17,000 genes for three tissues.

Genetic dissection of a genomic region with pleiotropic effects on domestication traits in maize reveals multiple linked QTL.

The domesticated crop maize and its wild progenitor, teosinte, have been used in numerous experiments to investigate the nature of divergent morphologies. This study examines a poorly understood region on the fifth chromosome of maize associated with a number of traits under selection during domestication, using a quantitative trait locus (QTL) mapping population specific to the fifth chromosome. In contrast with other major domestication loci in maize where large-effect, highly pleiotropic, single genes are responsible for phenotypic effects, our study found the region on chromosome five fractionates into multiple-QTL regions, none with singularly large effects. The smallest 1.5-LOD support interval for a QTL contained 54 genes, one of which was a MADS MIKC(C) transcription factor, a family of proteins implicated in many developmental programs. We also used simulated trait data sets to investigate the power of our mapping population to identify QTL for which there is a single underlying causal gene. This analysis showed that while QTL for traits controlled by single genes can be accurately mapped, our population design can detect no more than ∼4.5 QTL per trait even when there are 100 causal genes. Thus when a trait is controlled by ≥5 genes in the simulated data, the number of detected QTL can represent a simplification of the underlying causative factors. Our results show how a QTL region with effects on several domestication traits may be due to multiple linked QTL of small effect as opposed to a single gene with large and pleiotropic effects.

Defining the Role of prolamin-box binding factor1 Gene During Maize Domestication.

The prolamin-box binding factor1 (pbf1) gene encodes a transcription factor that controls the expression of seed storage protein (zein) genes in maize. Prior studies show that pbf1 underwent selection during maize domestication although how it affected trait change during domestication is unknown. To assay how pbf1 affects phenotypic differences between maize and teosinte, we compared nearly isogenic lines (NILs) that differ for a maize versus teosinte allele of pbf1 Kernel weight for the teosinte NIL (162mg) is slightly but significantly greater than that for the maize NIL (156mg). RNAseq data for developing kernels show that the teosinte allele of pbf1 is expressed at about twice the level of the maize allele. However, RNA and protein assays showed no difference in zein profile between the two NILs. The lower expression for the maize pbf1 allele suggests that selection may have favored this change; however, how reduced pbf1 expression alters phenotype remains unknown. One possibility is that pbf1 regulates genes other than zeins and thereby is a domestication trait. The observed drop in seed weight associated with the maize allele of pbf1 is counterintuitive but could represent a negative pleiotropic effect of selection on some other aspect of kernel composition.