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

Identification and Expression Analysis of the WOX Transcription Factor Family in Foxtail Millet (Setaria italica L.)

WUSCHEL-related homeobox (WOX) transcription factors are unique to plants and play pivotal roles in plant development and stress responses. In this investigation, we acquired protein sequences of foxtail millet WOX gene family members through homologous sequence alignment and a hidden Markov model (HMM) search. Utilizing conserved domain prediction, we identified 13 foxtail millet WOX genes, which were classified into ancient, intermediate, and modern clades. Multiple sequence alignment results revealed that all WOX proteins possess a homeodomain (HD). The SiWOX genes, clustered together in the phylogenetic tree, exhibited analogous protein spatial structures, gene structures, and conserved motifs. The foxtail millet WOX genes are distributed across 7 chromosomes, featuring 3 pairs of tandem repeats: SiWOX1 and SiWOX13, SiWOX4 and SiWOX5, and SiWOX11 and SiWOX12. Collinearity analysis demonstrated that WOX genes in foxtail millet exhibit the highest collinearity with green foxtail, followed by maize. The SiWOX genes primarily harbor two categories of cis-acting regulatory elements: Stress response and plant hormone response. Notably, prominent hormones triggering responses include methyl jasmonate, abscisic acid, gibberellin, auxin, and salicylic acid. Analysis of SiWOX expression patterns and hormone responses unveiled potential functional diversity among different SiWOX genes in foxtail millet. These findings lay a solid foundation for further elucidating the functions and evolution of SiWOX genes.

Recent Advances in Tomato Gene Editing

The use of gene-editing tools, such as zinc finger nucleases, TALEN, and CRISPR/Cas, allows for the modification of physiological, morphological, and other characteristics in a wide range of crops to mitigate the negative effects of stress caused by anthropogenic climate change or biotic stresses. Importantly, these tools have the potential to improve crop resilience and increase yields in response to challenging environmental conditions. This review provides an overview of gene-editing techniques used in plants, focusing on the cultivated tomatoes. Several dozen genes that have been successfully edited with the CRISPR/Cas system were selected for inclusion to illustrate the possibilities of this technology in improving fruit yield and quality, tolerance to pathogens, or responses to drought and soil salinity, among other factors. Examples are also given of how the domestication of wild species can be accelerated using CRISPR/Cas to generate new crops that are better adapted to the new climatic situation or suited to use in indoor agriculture.

New Advances in the Study of Regulation of Tomato Flowering-Related Genes Using Biotechnological Approaches

The tomato is a convenient object for studying reproductive processes, which has become a classic. Such complex processes as flowering and fruit setting require an understanding of the fundamental principles of molecular interaction, the structures of genes and proteins, the construction of signaling pathways for transcription regulation, including the synchronous actions of cis-regulatory elements (promoter and enhancer), trans-regulatory elements (transcription factors and regulatory RNAs), and transposable elements and epigenetic regulators (DNA methylation and acetylation, chromatin structure). Here, we discuss the current state of research on tomatoes (2017-2023) devoted to studying the function of genes that regulate flowering and signal regulation systems using genome-editing technologies, RNA interference gene silencing, and gene overexpression, including heterologous expression. Although the central candidate genes for these regulatory components have been identified, a complete picture of their relationship has yet to be formed. Therefore, this review summarizes the latest achievements related to studying the processes of flowering and fruit set. This work attempts to display the gene interaction scheme to better understand the events under consideration.

Detection of mosaic and population-level structural variants with Sniffles2

Calling structural variations (SVs) is technically challenging, but using long reads remains the most accurate way to identify complex genomic alterations. Here we present Sniffles2, which improves over current methods by implementing a repeat aware clustering coupled with a fast consensus sequence and coverage-adaptive filtering. Sniffles2 is 11.8 times faster and 29% more accurate than state-of-the-art SV callers across different coverages (5-50×), sequencing technologies (ONT and HiFi) and SV types. Furthermore, Sniffles2 solves the problem of family-level to population-level SV calling to produce fully genotyped VCF files. Across 11 probands, we accurately identified causative SVs around MECP2, including highly complex alleles with three overlapping SVs. Sniffles2 also enables the detection of mosaic SVs in bulk long-read data. As a result, we identified multiple mosaic SVs in brain tissue from a patient with multiple system atrophy. The identified SV showed a remarkable diversity within the cingulate cortex, impacting both genes involved in neuron function and repetitive elements.

TBtools-II: A "one for all, all for one" bioinformatics platform for biological big-data mining

Since the official release of the stand-alone bioinformatics toolkit TBtools in 2020, its superior functionality in data analysis has been demonstrated by its widespread adoption by many thousands of users and references in more than 5000 academic articles. Now, TBtools is a commonly used tool in biological laboratories. Over the past 3 years, thanks to invaluable feedback and suggestions from numerous users, we have optimized and expanded the functionality of the toolkit, leading to the development of an upgraded version-TBtools-II. In this upgrade, we have incorporated over 100 new features, such as those for comparative genomics analysis, phylogenetic analysis, and data visualization. Meanwhile, to better meet the increasing needs of personalized data analysis, we have launched the plugin mode, which enables users to develop their own plugins and manage their selection, installation, and removal according to individual needs. To date, the plugin store has amassed over 50 plugins, with more than half of them being independently developed and contributed by TBtools users. These plugins offer a range of data analysis options including co-expression network analysis, single-cell data analysis, and bulked segregant analysis sequencing data analysis. Overall, TBtools is now transforming from a stand-alone software to a comprehensive bioinformatics platform of a vibrant and cooperative community in which users are also developers and contributors. By promoting the theme "one for all, all for one", we believe that TBtools-II will greatly benefit more biological researchers in this big-data era.

Multi-year field trials provide a massive repository of trait data on a highly diverse population of tomato and uncover novel determinants of tomato productivity

Tomato (Solanum lycopersicum) is a prominent fruit with rich genetic resources for crop improvement. By using a phenotype-guided screen of over 7900 tomato accessions from around the world, we identified new associations for complex traits such as fruit weight and total soluble solids (Brix). Here, we present the phenotypic data from several years of trials. To illustrate the power of this dataset we use two case studies. First, evaluation of color revealed allelic variation in phytoene synthase 1 that resulted in differently colored or even bicolored fruit. Secondly, in view of the negative relationship between fruit weight and Brix, we pre-selected a subset of the collection that includes high and low Brix values in each category of fruit size. Genome-wide association analysis allowed us to detect novel loci associated with total soluble solid content and fruit weight. In addition, we developed eight F2 biparental intraspecific populations. Furthermore, by taking a phenotype-guided approach we were able to isolate individuals with high Brix values that were not compromised in terms of yield. In addition, the demonstration of novel results despite the high number of previous genome-wide association studies of these traits in tomato suggests that adoption of a phenotype-guided pre-selection of germplasm may represent a useful strategy for finding target genes for breeding.

A perspective from the EU: unintended genetic changes in plants caused by NGT-their relevance for a comprehensive molecular characterisation and risk assessment

Several regions in the world are currently holding discussions in regard to the regulation of new genomic techniques (NGTs) and their application in agriculture. The European Commission, for instance, is proposing the introduction of specific regulation for NGT plants. Various questions need to be answered including e.g., the extent to which NGT-induced intended and unintended genetic modifications must be subjected to a mandatory risk assessment as part of an approval procedure. This review mostly focuses on findings in regard to unintended genetic changes that can be caused by the application of NGTs. More specifically, the review deals with the application of the nuclease CRISPR/Cas, which is currently the most important tool for developing NGT plants, and its potential to introduce double strand breaks (DSBs) at a targeted DNA sequence. For this purpose, we identified the differences in comparison to non-targeted mutagenesis methods used in conventional breeding. The review concludes that unintended genetic changes caused by NGT processes are relevant to risk assessment. Due to the technical characteristics of NGTs, the sites of the unintended changes, their genomic context and their frequency (in regard to specific sites) mean that the resulting gene combinations (intended or unintended) may be unlikely to occur with conventional methods. This, in turn, implies that the biological effects (phenotypes) can also be different and may cause risks to health and the environment. Therefore, we conclude that the assessment of intended as well as unintended genetic changes should be part of a mandatory comprehensive molecular characterisation and risk assessment of NGT plants that are meant for environmental releases or for market authorisation.

Paralog transcriptional differentiation in the D. melanogaster-specific gene family Sdic across populations and spermatogenesis stages

How recently originated gene copies become stable genomic components remains uncertain as high sequence similarity of young duplicates precludes their functional characterization. The tandem multigene family Sdic is specific to Drosophila melanogaster and has been annotated across multiple reference-quality genome assemblies. Here we show the existence of a positive correlation between Sdic copy number and total expression, plus vast intrastrain differences in mRNA abundance among paralogs, using RNA-sequencing from testis of four strains with variable paralog composition. Single cell and nucleus RNA-sequencing data expose paralog expression differentiation in meiotic cell types within testis from third instar larva and adults. Additional RNA-sequencing across synthetic strains only differing in their Y chromosomes reveal a tissue-dependent trans-regulatory effect on Sdic: upregulation in testis and downregulation in male accessory gland. By leveraging paralog-specific expression information from tissue- and cell-specific data, our results elucidate the intraspecific functional diversification of a recently expanded tandem gene family.

TM3 and STM3 Promote Flowering Together with FUL2 and MBP20, but Act Antagonistically in Inflorescence Branching in Tomato

The moment at which a plant transitions to reproductive development is paramount to its life cycle and is strictly controlled by many genes. The transcription factor SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) plays a central role in this process in Arabidopsis. However, the role of SOC1 in tomato (Solanum lycopersicum) has been sparsely studied. Here, we investigated the function of four tomato SOC1 homologs in the floral transition and inflorescence development. We thoroughly characterized the SOC1-like clade throughout the Solanaceae and selected four tomato homologs that are dynamically expressed upon the floral transition. We show that of these homologs, TOMATO MADS 3 (TM3) and SISTER OF TM3 (STM3) promote the primary and sympodial transition to flowering, while MADS-BOX PROTEIN 23 (MBP23) and MBP18 hardly contribute to flowering initiation in the indeterminate cultivar Moneyberg. Protein-protein interaction assays and whole-transcriptome analysis during reproductive meristem development revealed that TM3 and STM3 interact and share many targets with FRUITFULL (FUL) homologs, including cytokinin regulators. Furthermore, we observed that mutating TM3/STM3 affects inflorescence development, but counteracts the inflorescence-branching phenotype of ful2 mbp20. Collectively, this indicates that TM3/STM3 promote the floral transition together with FUL2/MBP20, while these transcription factors have opposite functions in inflorescence development.

The role of recombination landscape in species hybridisation and speciation

It is now well recognised that closely related species can hybridize and exchange genetic material, which may promote or oppose adaptation and speciation. In some cases, interspecific hybridisation is very common, making it surprising that species identity is preserved despite active gene exchange. The genomes of most eukaryotic species are highly heterogeneous with regard to gene density, abundance of repetitive DNA, chromatin compactisation etc, which can make certain genomic regions more prone or more resistant to introgression of genetic material from other species. Heterogeneity in local recombination rate underpins many of the observed patterns across the genome (e.g. actively recombining regions are typically gene rich and depleted for repetitive DNA) and it can strongly affect the permeability of genomic regions to interspecific introgression. The larger the region lacking recombination, the higher the chance for the presence of species incompatibility gene(s) in that region, making the entire non- or rarely recombining block impermeable to interspecific introgression. Large plant genomes tend to have highly heterogeneous recombination landscape, with recombination frequently occurring at the ends of the chromosomes and central regions lacking recombination. In this paper we review the relationship between recombination and introgression in plants and argue that large rarely recombining regions likely play a major role in preserving species identity in actively hybridising plant species.

The trans-regulatory landscape of gene networks in plants

The transcriptional effector domains of transcription factors play a key role in controlling gene expression; however, their functional nature is poorly understood, hampering our ability to explore this fundamental dimension of gene regulatory networks. To map the trans-regulatory landscape in a complex eukaryote, we systematically characterized the putative transcriptional effector domains of over 400 Arabidopsis thaliana transcription factors for their capacity to modulate transcription. We demonstrate that transcriptional effector activity can be integrated into gene regulatory networks capable of elucidating the functional dynamics underlying gene expression patterns. We further show how our characterized domains can enhance genome engineering efforts and reveal how plant transcriptional activators share regulatory features conserved across distantly related eukaryotes. Our results provide a framework to systematically characterize the regulatory role of transcription factors at a genome-scale in order to understand the transcriptional wiring of biological systems.

Antagonistic regulation of target genes by the SISTER OF TM3-JOINTLESS2 complex in tomato inflorescence branching

Inflorescence branch number is a yield-related trait controlled by cell fate determination in meristems. Two MADS-box transcription factors (TFs)-SISTER OF TM3 (STM3) and JOINTLESS 2 (J2)-have opposing regulatory roles in inflorescence branching. However, the mechanisms underlying their regulatory functions in inflorescence determinacy remain unclear. Here, we characterized the functions of these TFs in tomato (Solanum lycopersicum) floral meristem and inflorescence meristem (IM) through chromatin immunoprecipitation and sequencing analysis of their genome-wide occupancy. STM3 and J2 activate or repress the transcription of a set of common putative target genes, respectively, through recognition and binding to CArG box motifs. FRUITFULL1 (FUL1) is a shared putative target of STM3 and J2 and these TFs antagonistically regulate FUL1 in inflorescence branching. Moreover, STM3 physically interacts with J2 to mediate its cytosolic redistribution and restricts J2 repressor activity by reducing its binding to target genes. Conversely, J2 limits STM3 regulation of target genes by transcriptional repression of the STM3 promoter and reducing STM3-binding activity. Our study thus reveals an antagonistic regulatory relationship in which STM3 and J2 control tomato IM determinacy and branch number.

Tomato APETALA2 family member SlTOE1 regulates inflorescence branching by repressing SISTER OF TM3

Inflorescence architecture directly impacts yield potential in most crops. As a model of sympodial plants, tomato (Solanum lycopersicum) inflorescence exhibits highly structural plasticity. However, the genetic regulatory network of inflorescence architecture in tomato remains unclear. Here, we investigated a modulator of inflorescence branching in tomato, TARGET OF EAT1 (SlTOE1), an APETALA2 (AP2) family member found to be predominantly expressed in the floral meristem (FM) of tomato. sltoe1 knockout mutants displayed highly branched inflorescences and defective floral organs. Transcriptome analysis revealed that SISTER OF TM3 (STM3) and certain floral development-related genes were upregulated in the flower meristem of sltoe1. SlTOE1 could directly bind the promoters of STM3 and Tomato MADS-box gene 3 (TM3) to repress their transcription. Simultaneous mutation of STM3 and TM3 partially restored the inflorescence branching of the sltoe1cr mutants, suggesting that SlTOE1 regulates inflorescence development, at least in part through an SlTOE1STM3/TM3 module. Genetic analysis showed that SlTOE1 and ENHANCER OF JOINTLESS 2 (EJ2) additively regulate tomato inflorescence branching; their double mutants showed more extensive inflorescence branching. Our findings uncover a pathway controlling tomato inflorescence branching and offer deeper insight into the functions of AP2 subfamily members.

Natural variations in the PbCPK28 promoter regulate sugar content through interaction with PbTST4 and PbVHA-A1 in pear

Soluble sugars play an important role in plant growth, development and fruit quality. Pear fruits have demonstrated a considerable improvement in sugar quality during their long history of selection. However, little is known about the underlying molecular mechanisms accompanying the changes in fruit sugar content as a result of selection by horticulturists. Here, we identified a calcium-dependent protein kinase (PbCPK28), which is located on LG15 and is present within a selective sweep region, thus linked to the quantitative trait loci for soluble solids. Association analysis indicates that a single nucleotide polymorphism-13 variation (SNP13^T/C ) in the PbCPK28 regulatory region led to fructose content diversity in pear. Elevated expression of PbCPK28 resulted in significantly increased fructose levels in pear fruits. Furthermore, PbCPK28 interacts with and phosphorylates PbTST4, a proton antiporter, thereby coupling the sugar import into the vacuole with proton export. We demonstrated that residues S277 and S314 of PbTST4 are crucial for its function. Additionally, PbCPK28 interacts with and phosphorylates the vacuolar hydrogen proton pump PbVHA-A1, which could provide proton motive forces for PbTST4. We also found that the T11 and Y120 phosphorylation sites in PbVHA-A1 are essential for its function. Evolution analysis and yeast-two-hybrid results support that the CPK-TST/CPK-VHA-A regulatory network is highly conserved in plants, especially the corresponding phosphorylation sites. Together, our work identifies an agriculturally important natural variation and an important regulatory network, allowing genetic improvement of fruit sugar contents in pears through modulation of PbCPK28 expression and phosphorylation of PbTST4 and PbVHA-A1.

Reducing crop losses by gene-editing control of organ developmental physiology

Some physiological processes in reproductive organs, if not controlled, can lead to crop loss even in the absence of environmental stress. These processes may occur pre- or post- harvest, and in diverse species and include abscission processes in cereal grain, e.g., shattering and in immature fruit, e.g., preharvest drop, preharvest sprouting of cereals, and postharvest senescence in fruit. Some of the molecular mechanisms and genetic determinants underlying these processes are now better detailed, making it possible to refine them by gene editing. Here, we discuss using advanced genomics to identify genetic determinants underlying crop physiological traits. Examples of improved phenotypes developed for preharvest problems are provided, and suggestions for reducing postharvest fruit losses by gene and promoter editing were made.

Involvement of IDA-HAE Module in Natural Development of Tomato Flower Abscission

The unwanted detachment of organs such as flowers, leaves, and fruits from the main body of a plant (abscission) has significant effects on agricultural practice. Both timely and precise regulation of organ abscission from a plant is crucial as it influences the agricultural yield. The tomato (Solanum lycopersicum) has become a model system for research on organ abscission. Here, we characterized four tomato natural abscission variants named jointless (j), functionally impaired jointless (fij), functionally impaired jointless like (fij like), and normal joint (NJ), based on their cellular features within the flower abscission zones (AZ). Using eight INFLORESCENCE DEFICIENT IN ABSCISSION (SlIDA) genes and eight HAESA genes (SlHAE) identified in the genome sequence of tomato, we analyzed the pattern of gene expression during flower abscission. The AZ-specific expression for three tomato abscission polygalacturonases (SlTAPGs) in the development of flower AZ, and the progression of abscission validated our natural abscission system. Compared to that of j, fij, and fij like variants, the AZ-specific expression for SlIDA, SlIDL2, SlIDL3, SlIDL4, and SlIDL5 in the NJ largely corelated and increased with the process of abscission. Of eight SlHAE genes examined, the expression for SlHSL6 and SlHSL7 were found to be AZ-specific and increased as abscission progressed in the NJ variant. Unlike the result of gene expression obtained from natural abscission system, an in silico analysis of transcriptional binding sites uncovered that SlIDA genes (SlIDA, SlIDL6, and SlIDL7) are predominantly under the control of environmental stress, while most of the SlHSL genes are affiliated with the broader context in developmental processes and stress responses. Our result presents the potential bimodal transcriptional regulation of the tomato IDA-HAE module associated with flower abscission in tomatoes.

Flower Development in the Solanaceae

Flower development is the process leading from a reproductive meristem to a mature flower with fully developed floral organs. This multi-step process is complex and involves thousands of genes in intertwined regulatory pathways; navigating through the FLOR-ID website will give an impression of this complexity and of the astonishing amount of work that has been carried on the topic (Bouché et al., Nucleic Acids Res 44:D1167-D1171, 2016). Our understanding of flower development mostly comes from the model species Arabidopsis thaliana, but numerous other studies outside of Brassicaceae have helped apprehend the conservation of these mechanisms in a large evolutionary context (Moyroud and Glover, Curr Biol 27:R941-R951, 2017; Smyth, New Phytol 220:70-86, 2018; Soltis et al., Ann Bot 100:155-163, 2007). Integrating additional species and families to the research on this topic can only advance our understanding of flower development and its evolution.In this chapter, we review the contribution that the Solanaceae family has made to the comprehension of flower development. While many of the general features of flower development (i.e., the key molecular players involved in flower meristem identity, inflorescence architecture or floral organ development) are similar to Arabidopsis, our main objective in this chapter is to highlight the points of divergence and emphasize specificities of the Solanaceae. We will not discuss the large topics of flowering time regulation, inflorescence architecture and fruit development, and we will restrict ourselves to the mechanisms included in a time window after the floral transition and before the fertilization. Moreover, this review will not be exhaustive of the large amount of work carried on the topic, and the choices that we made to describe in large details some stories from the literature are based on the soundness of the functional work performed, and surely as well on our own preferences and expertise.First, we will give a brief overview of the Solanaceae family and some of its specificities. Then, our focus will be on the molecular mechanisms controlling floral organ identity, for which extended functional work in petunia led to substantial revisions to the famous ABC model. Finally, after reviewing some studies on floral organ initiation and growth, we will discuss floral organ maturation, using the examples of the inflated calyx of the Chinese lantern Physalis and petunia petal pigmentation.

Unintended Consequences of Plant Domestication

Human selection on wild populations mostly favored a common set of plant traits during domestication. This process of direct selection also altered other independent traits that were not directly perceived or desired during crop domestication and improvement. A deeper knowledge of the inadvertent and undesirable phenotypic effects and their underlying genetic causes can help design strategies to mitigate their effects and improve genetic gain in crop plants. We review different factors explaining the negative consequences of plant domestication at the phenotypic and genomic levels. We further describe the genetic causes of undesirable effects that originate from the selection of favorable alleles during plant domestication. In addition, we propose strategies that could be useful in attenuating such effects for crop improvement. With novel -omics and genome-editing tools, it is relatively approachable to understand and manipulate the genetic and biochemical mechanisms responsible for the undesirable phenotypes in domesticated plants.

Tracing the Evolution of the SEPALLATA Subfamily across Angiosperms Associated with Neo- and Sub-Functionalization for Reproductive and Agronomically Relevant Traits

SEPALLATA transcription factors (SEP TFs) have been extensively studied in angiosperms as pivotal components of virtually all the MADS-box tetrameric complex master regulators of floral organ identities. However, there are published reports that suggest that some SEP members also regulate earlier reproductive events, such as inflorescence meristem determinacy and inflorescence architecture, with potential for application in breeding programs in crops. The SEP subfamily underwent a quite complex pattern of duplications during the radiation of the angiosperms. Taking advantage of the many whole genomic sequences now available, we present a revised and expanded SEP phylogeny and link it to the known functions of previously characterized genes. This snapshot supports the evidence that the major SEP3 clade is highly specialized for the specification of the three innermost floral whorls, while its sister LOFSEP clade is functionally more versatile and has been recruited for diverse roles, such as the regulation of extra-floral bract formation and inflorescence determinacy and shape. This larger pool of angiosperm SEP genes confirms previous evidence that their evolution was driven by whole-genome duplications rather than small-scale duplication events. Our work may help to identify those SEP lineages that are the best candidates for the improvement of inflorescence traits, even in far distantly related crops.

Comparative Analysis of Environment-Responsive Alternative Splicing in the Inflorescences of Cultivated and Wild Tomato Species

Cultivated tomato (Solanum lycopersicum) is bred for fruit production in optimized environments, in contrast to harsh environments where their ancestral relatives thrive. The process of domestication and breeding has profound impacts on the phenotypic plasticity of plant development and the stress response. Notably, the alternative splicing (AS) of precursor message RNA (pre-mRNA), which is one of the major factors contributing to transcriptome complexity, is responsive to developmental cues and environmental change. To determine a possible association between AS events and phenotypic plasticity, we investigated environment-responsive AS events in the inflorescences of cultivated tomato and its ancestral relatives S. pimpinellifolium. Despite that similar AS frequencies were detected in the cultivated tomato variety Moneymaker and two S. pimpinellifolium accessions under the same growth conditions, 528 genes including splicing factors showed differential splicing in the inflorescences of plants grown in open fields and plastic greenhouses in the Moneymaker variety. In contrast, the two S. pimpinellifolium accessions, LA1589 and LA1781, had 298 and 268 genes showing differential splicing, respectively. Moreover, seven heat responsive genes showed opposite expression patterns in response to changing growth conditions between Moneymaker and its ancestral relatives. Accordingly, there were eight differentially expressed splice variants from genes involved in heat response in Moneymaker. Our results reveal distinctive features of AS events in the inflorescences between cultivated tomato and its ancestral relatives, and show that AS regulation in response to environmental changes is genotype dependent.

Transcriptome Analysis to Explore the Cause of the Formation of Different Inflorescences in Tomato

The number of inflorescence branches is an important agronomic character of tomato. The meristem differentiation and development pattern of tomato inflorescence is complex and its regulation mechanism is very different from those of other model plants. Therefore, in order to explore the cause of tomato inflorescence branching, transcriptome analysis was conducted on two kinds of tomato inflorescences (single racemes and compound inflorescences). According to the transcriptome data analysis, there were many DEGs of tomato inflorescences at early, middle, and late stages. Then, GO and KEGG enrichments of DEGs were performed. DEGs are mainly enriched in metabolic pathways, biohormone signaling, and cell cycle pathways. According to previous studies, DEGs were mainly enriched in metabolic pathways, and FALSIFLORA (FA) and ANANTHA (AN) genes were the most notable of 41 DEGs related to inflorescence branching. This study not only provides a theoretical basis for understanding inflorescence branching, but also provides a new idea for the follow-up study of inflorescence.

European traditional tomatoes galore: a result of farmers' selection of a few diversity-rich loci

A comprehensive collection of 1254 tomato accessions, corresponding to European traditional and modern varieties, early domesticated varieties, and wild relatives, was analyzed by genotyping by sequencing. A continuous genetic gradient between the traditional and modern varieties was observed. European traditional tomatoes displayed very low genetic diversity, with only 298 polymorphic loci (95% threshold) out of 64 943 total variants. European traditional tomatoes could be classified into several genetic groups. Two main clusters consisting of Spanish and Italian accessions showed higher genetic diversity than the remaining varieties, suggesting that these regions might be independent secondary centers of diversity with a different history. Other varieties seem to be the result of a more recent complex pattern of migrations and hybridizations among the European regions. Several polymorphic loci were associated in a genome-wide association study with fruit morphological traits in the European traditional collection. The corresponding alleles were found to contribute to the distinctive phenotypic characteristic of the genetic varietal groups. The few highly polymorphic loci associated with morphological traits in an otherwise a low-diversity population suggests a history of balancing selection, in which tomato farmers likely maintained the morphological variation by inadvertently applying a high selective pressure within different varietal types.

Atlas of phenotypic, genotypic and geographical diversity present in the European traditional tomato

The Mediterranean basin countries are considered secondary centres of tomato diversification. However, information on phenotypic and allelic variation of local tomato materials is still limited. Here we report on the evaluation of the largest traditional tomato collection, which includes 1499 accessions from Southern Europe. Analyses of 70 traits revealed a broad range of phenotypic variability with different distributions among countries, with the culinary end use within each country being the main driver of tomato diversification. Furthermore, eight main tomato types (phenoclusters) were defined by integrating phenotypic data, country of origin, and end use. Genome-wide association study (GWAS) meta-analyses identified associations in 211 loci, 159 of which were novel. The multidimensional integration of phenoclusters and the GWAS meta-analysis identified the molecular signatures for each traditional tomato type and indicated that signatures originated from differential combinations of loci, which in some cases converged in the same tomato phenotype. Our results provide a roadmap for studying and exploiting this untapped tomato diversity.

Fine-tuning shoot meristem size to feed the world

In order to maintain food security for the world's growing population, crop yields need to be significantly improved. Domestication and crop improvement involve modification of traits such as fruit size and seed number to optimize productivity. Although these traits are selected at the mature stage, they are determined during the development of shoot meristem, a tissue that forms successive meristems and reproductive organs that make edible fruits or seeds. Therefore, the architecture of reproductive organs and yield-related traits are determined during the maturation of shoot meristem. Here, we highlight recent progress in understanding how shoot meristem size affects yield-related traits and outline the strategies to fine-tune meristem regulatory genes to meet the demands of a growing population and promote sustainable agriculture.

Genomic structural variation in tomato and its role in plant immunity

It is well known that large genomic variations can greatly impact the phenotype of an organism. Structural Variants (SVs) encompass any genomic variation larger than 30 base pairs, and include changes caused by deletions, inversions, duplications, transversions, and other genome modifications. Due to their size and complex nature, until recently, it has been difficult to truly capture these variations. Recent advances in sequencing technology and computational analyses now permit more extensive studies of SVs in plant genomes. In tomato, advances in sequencing technology have allowed researchers to sequence hundreds of genomes from tomatoes, and tomato relatives. These studies have identified SVs related to fruit size and flavor, as well as plant disease response, resistance/susceptibility, and the ability of plants to detect pathogens (immunity). In this review, we discuss the implications for genomic structural variation in plants with a focus on its role in tomato immunity. We also discuss how advances in sequencing technology have led to new discoveries of SVs in more complex genomes, the current evidence for the role of SVs in biotic and abiotic stress responses, and the outlook for genetic modification of SVs to advance plant breeding objectives.

FRUITFULL-like genes regulate flowering time and inflorescence architecture in tomato

The timing of flowering and the inflorescence architecture are critical for the reproductive success of tomato (Solanum lycopersicum), but the gene regulatory networks underlying these traits have not been fully explored. Here, we show that the tomato FRUITFULL-like (FUL-like) genes FUL2 and MADS-BOX PROTEIN 20 (MBP20) promote the vegetative-to-reproductive transition and repress inflorescence branching by inducing floral meristem (FM) maturation. FUL1 fulfils a less prominent role and appears to depend on FUL2 and MBP20 for its upregulation in the inflorescence- and floral meristems. MBP10, the fourth tomato FUL-like gene, has probably lost its function. The tomato FUL-like proteins cannot homodimerize in in vitro assays, but heterodimerize with various other MADS-domain proteins, potentially forming distinct complexes in the transition meristem and FM. Transcriptome analysis of the primary shoot meristems revealed various interesting downstream targets, including four repressors of cytokinin signaling that are upregulated during the floral transition in ful1 ful2 mbp10 mbp20 mutants. FUL2 and MBP20 can also bind in vitro to the upstream regions of these genes, thereby probably directly stimulating cell division in the meristem upon the transition to flowering. The control of inflorescence branching does not occur via the cytokinin oxidase/dehydrogenases (CKXs) but may be regulated by repression of transcription factors such as TOMATO MADS-box gene 3 (TM3) and APETALA 2b (AP2b).

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

Utilizing CRISPR-Cas in Tropical Crop Improvement: A Decision Process for Fitting Genome Engineering to Your Species

Adopting modern gene-editing technologies for trait improvement in agriculture requires important workflow developments, yet these developments are not often discussed. Using tropical crop systems as a case study, we describe a workflow broken down into discrete processes with specific steps and decision points that allow for the practical application of the CRISPR-Cas gene editing platform in a crop of interest. While we present the steps of developing genome-edited plants as sequential, in practice parts can be done in parallel, which are discussed in this perspective. The main processes include 1) understanding the genetic basis of the trait along with having the crop’s genome sequence, 2) testing and optimization of the editing reagents, development of efficient 3) tissue culture and 4) transformation methods, and 5) screening methods to identify edited events with commercial potential. Our goal in this perspective is to help any lab that wishes to implement this powerful, easy-to-use tool in their pipeline, thus aiming to democratize the technology.

SISTER OF TM3 activates FRUITFULL1 to regulate inflorescence branching in tomato

Selection for favorable inflorescence architecture to improve yield is one of the crucial targets in crop breeding. Different tomato varieties require distinct inflorescence-branching structures to enhance productivity. While a few important genes for tomato inflorescence-branching development have been identified, the regulatory mechanism underlying inflorescence branching is still unclear. Here, we confirmed that SISTER OF TM3 (STM3), a homolog of Arabidopsis SOC1, is a major positive regulatory factor of tomato inflorescence architecture by map-based cloning. High expression levels of STM3 underlie the highly inflorescence-branching phenotype in ST024. STM3 is expressed in both vegetative and reproductive meristematic tissues and in leaf primordia and leaves, indicative of its function in flowering time and inflorescence-branching development. Transcriptome analysis shows that several floral development-related genes are affected by STM3 mutation. Among them, FRUITFULL1 (FUL1) is downregulated in stm3cr mutants, and its promoter is bound by STM3 by ChIP-qPCR analysis. EMSA and dual-luciferase reporter assays further confirmed that STM3 could directly bind the promoter region to activate FUL1 expression. Mutation of FUL1 could partially restore inflorescence-branching phenotypes caused by high STM3 expression in ST024. Our findings provide insights into the molecular and genetic mechanisms underlying inflorescence development in tomato.

Leaf form diversification in an ornamental heirloom tomato results from alterations in two different HOMEOBOX genes

Domesticated plants display diverse phenotypic traits. However, the influence of breeding effort on this phenotypic diversity remains unknown. Here, we demonstrate that a single nucleotide deletion in the homeobox motif of BIPINNATA, a BEL-LIKE HOMEODOMAIN gene, led to a highly complex leaf phenotype in an heirloom tomato (Solanum lycopersicum), Silvery Fir Tree (SiFT), which is used as a landscaping and ornamental plant. A comparative gene network analysis revealed that repression of SOLANIFOLIA, the ortholog of WUSCHEL RELATED HOMEOBOX 1, caused the narrow leaflet phenotype seen in SiFT. Comparative genomics indicated that the bip mutation in SiFT likely arose de novo and is unique to SiFT and not introgressed from other tomato genomes. These results provide new insights into the natural variation in phenotypic traits introduced into crops during improvement processes after domestication and establish homeobox genes as evolutionary hotspots.

Leaf form diversification in an ornamental heirloom tomato results from alterations in two different HOMEOBOX genes

Domesticated plants display diverse phenotypic traits. However, the influence of breeding effort on this phenotypic diversity remains unknown. Here, we demonstrate that a single nucleotide deletion in the homeobox motif of BIPINNATA, a BEL-LIKE HOMEODOMAIN gene, led to a highly complex leaf phenotype in an heirloom tomato (Solanum lycopersicum), Silvery Fir Tree (SiFT), which is used as a landscaping and ornamental plant. A comparative gene network analysis revealed that repression of SOLANIFOLIA, the ortholog of WUSCHEL RELATED HOMEOBOX 1, caused the narrow leaflet phenotype seen in SiFT. Comparative genomics indicated that the bip mutation in SiFT likely arose de novo and is unique to SiFT and not introgressed from other tomato genomes. These results provide new insights into the natural variation in phenotypic traits introduced into crops during improvement processes after domestication and establish homeobox genes as evolutionary hotspots.

Meristem transitions and plant architecture-learning from domestication for crop breeding

Genetic networks that regulate meristem transitions were recurrent targets of selection during crop domestication and allow fine-tuning of plant architecture for improved crop productivity.

Domestication shapes recombination patterns in tomato

Meiotic recombination is a biological process of key importance in breeding, to generate genetic diversity and develop novel or agronomically relevant haplotypes. In crop tomato, recombination is curtailed as manifested by linkage disequilibrium decay over a longer distance and reduced diversity compared with wild relatives. Here, we compared domesticated and wild populations of tomato and found an overall conserved recombination landscape, with local changes in effective recombination rate in specific genomic regions. We also studied the dynamics of recombination hotspots resulting from domestication and found that loss of such hotspots is associated with selective sweeps, most notably in the pericentromeric heterochromatin. We detected footprints of genetic changes and structural variants, among them associated with transposable elements, linked with hotspot divergence during domestication, likely causing fine-scale alterations to recombination patterns and resulting in linkage drag.

GWAS for genetics of complex quantitative traits: Genome to pangenome and SNPs to SVs and k-mers

The development of improved methods for genome-wide association studies (GWAS) for genetics of quantitative traits has been an active area of research during the last 25 years. This activity initially started with the use of mixed linear model (MLM), which was variously modified. During the last decade, however, with the availability of high throughput next generation sequencing (NGS) technology, development and use of pangenomes and novel markers including structural variations (SVs) and k-mers for GWAS has taken over as a new thrust area of research. Pangenomes and SVs are now available in humans, livestock, and a number of plant species, so that these resources along with k-mers are being used in GWAS for exploring additional genetic variation that was hitherto not available for analysis. These developments have resulted in significant improvement in GWAS methodology for detection of marker-trait associations (MTAs) that are relevant to human healthcare and crop improvement.

Natural Variation in Crops: Realized Understanding, Continuing Promise

Crops feed the world's population and shape human civilization. The improvement of crop productivity has been ongoing for almost 10,000 years and has evolved from an experience-based to a knowledge-driven practice over the past three decades. Natural alleles and their reshuffling are long-standing genetic changes that affect how crops respond to various environmental conditions and agricultural practices. Decoding the genetic basis of natural variation is central to understanding crop evolution and, in turn, improving crop breeding. Here, we review current advances in the approaches used to map the causal alleles of natural variation, provide refined insights into the genetics and evolution of natural variation, and outline how this knowledge promises to drive the development of sustainable agriculture under the dome of emerging technologies.

A systems genetics approach to deciphering the effect of dosage variation on leaf morphology in Populus

Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

Structural variation at the maize WUSCHEL1 locus alters stem cell organization in inflorescences

Structural variation in plant genomes is a significant driver of phenotypic variability in traits important for the domestication and productivity of crop species. Among these are traits that depend on functional meristems, populations of stem cells maintained by the CLAVATA-WUSCHEL (CLV-WUS) negative feedback-loop that controls the expression of the WUS homeobox transcription factor. WUS function and impact on maize development and yield remain largely unexplored. Here we show that the maize dominant Barren inflorescence3 (Bif3) mutant harbors a tandem duplicated copy of the ZmWUS1 gene, ZmWUS1-B, whose novel promoter enhances transcription in a ring-like pattern. Overexpression of ZmWUS1-B is due to multimerized binding sites for type-B RESPONSE REGULATORs (RRs), key transcription factors in cytokinin signaling. Hypersensitivity to cytokinin causes stem cell overproliferation and major rearrangements of Bif3 inflorescence meristems, leading to the formation of ball-shaped ears and severely affecting productivity. These findings establish ZmWUS1 as an essential meristem size regulator in maize and highlight the striking effect of cis-regulatory variation on a key developmental program.

The WUSCHEL transcription factor promotes plant stem cell proliferation. Here the authors show that the maize Bif3 mutant contains a duplication of the ZmWUS1 locus leading to cytokinin hypersensitivity and overproliferation at the shoot meristem demonstrating the role of WUSCHEL in maize and how structural variation can impact plant morphology.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

Coordination of differentiation rate and local patterning in compound-leaf development

The variability in leaf form in nature is immense. Leaf patterning occurs by differential growth, taking place during a limited window of morphogenetic activity at the leaf marginal meristem. While many regulators have been implicated in the designation of the morphogenetic window and in leaf patterning, how these effectors interact to generate a particular form is still not well understood. We investigated the interaction among different effectors of tomato (Solanum lycopersicum) compound-leaf development, using genetic and molecular analyses. Mutations in the tomato auxin response factor SlARF5/SlMP, which normally promotes leaflet formation, suppressed the increased leaf complexity of mutants with extended morphogenetic window. Impaired activity of the NAC/CUC transcription factor GOBLET (GOB), which specifies leaflet boundaries, also reduced leaf complexity in these backgrounds. Analysis of genetic interactions showed that the patterning factors SlMP, GOB and the MYB transcription factor LYRATE (LYR) coordinately regulate leaf patterning by modulating in parallel different aspects of leaflet formation and shaping. This work places an array of developmental regulators in a morphogenetic context. It reveals how organ-level differentiation rate and local growth are coordinated to sculpture an organ. These concepts are applicable to the coordination of pattering and differentiation in other species and developmental processes.

CRISPR-Mediated Engineering across the Central Dogma in Plant Biology for Basic Research and Crop Improvement

The central dogma (CD) of molecular biology is the transfer of genetic information from DNA to RNA to protein. Major CD processes governing genetic flow include the cell cycle, DNA replication, chromosome packaging, epigenetic changes, transcription, posttranscriptional alterations, translation, and posttranslational modifications. The CD processes are tightly regulated in plants to maintain genetic integrity throughout the life cycle and to pass genetic materials to next generation. Engineering of various CD processes involved in gene regulation will accelerate crop improvement to feed the growing world population. CRISPR technology enables programmable editing of CD processes to alter DNA, RNA, or protein, which would have been impossible in the past. Here, an overview of recent advancements in CRISPR tool development and CRISPR-based CD modulations that expedite basic and applied plant research is provided. Furthermore, CRISPR applications in major thriving areas of research, such as gene discovery (allele mining and cryptic gene activation), introgression (de novo domestication and haploid induction), and application of desired traits beneficial to farmers or consumers (biotic/abiotic stress-resilient crops, plant cell factories, and delayed senescence), are described. Finally, the global regulatory policies, challenges, and prospects for CRISPR-mediated crop improvement are discussed.

Genome of Solanum pimpinellifolium provides insights into structural variants during tomato breeding

Solanum pimpinellifolium (SP) is the wild progenitor of cultivated tomato. Because of its remarkable stress tolerance and intense flavor, SP has been used as an important germplasm donor in modern tomato breeding. Here, we present a high-quality chromosome-scale genome sequence of SP LA2093. Genome comparison identifies more than 92,000 structural variants (SVs) between LA2093 and the modern cultivar, Heinz 1706. Genotyping these SVs in ~600 representative tomato accessions identifies alleles under selection during tomato domestication, improvement and modern breeding, and discovers numerous SVs overlapping genes known to regulate important breeding traits such as fruit weight and lycopene content. Expression quantitative trait locus (eQTL) analysis detects hotspots harboring master regulators controlling important fruit quality traits, including cuticular wax accumulation and flavonoid biosynthesis, and SVs contributing to these complex regulatory networks. The LA2093 genome sequence and the identified SVs provide rich resources for future research and biodiversity-based breeding.

Solanum pimpinellifolium (SP) is the progenitor of cultivated tomato and an important germplasm. Here, the authors assemble SP genome, identify structural variants (SVs) by comparing with modern cultivar, reveal SVs associated with important breeding traits, and detect SVs harboring master regulators of fruit quality traits.

Multi-gene metabolic engineering of tomato plants results in increased fruit yield up to 23%

The capacity to assimilate carbon and nitrogen, to transport the resultant sugars and amino acids to sink tissues, and to convert the incoming sugars and amino acids into storage compounds in the sink tissues, are key determinants of crop yield. Given that all of these processes have the potential to co-limit growth, multiple genetic interventions in source and sink tissues, plus transport processes may be necessary to reach the full yield potential of a crop. We used biolistic combinatorial co-transformation (up to 20 transgenes) for increasing C and N flows with the purpose of increasing tomato fruit yield. We observed an increased fruit yield of up to 23%. To better explore the reconfiguration of metabolic networks in these transformants, we generated a dataset encompassing physiological parameters, gene expression and metabolite profiling on plants grown under glasshouse or polytunnel conditions. A Sparse Partial Least Squares regression model was able to explain the combination of genes that contributed to increased fruit yield. This combinatorial study of multiple transgenes targeting primary metabolism thus offers opportunities to probe the genetic basis of metabolic and phenotypic variation, providing insight into the difficulties in choosing the correct combination of targets for engineering increased fruit yield.

Chromosome-Scale Assembly of the Bread Wheat Genome Reveals Thousands of Additional Gene Copies

Bread wheat (Triticum aestivum) is a major food crop and an important plant system for agricultural genetics research. However, due to the complexity and size of its allohexaploid genome, genomic resources are limited compared to other major crops. The IWGSC recently published a reference genome and associated annotation (IWGSC CS v1.0, Chinese Spring) that has been widely adopted and utilized by the wheat community. Although this reference assembly represents all three wheat subgenomes at chromosome-scale, it was derived from short reads, and thus is missing a substantial portion of the expected 16 Gbp of genomic sequence. We earlier published an independent wheat assembly (Triticum_aestivum_3.1, Chinese Spring) that came much closer in length to the expected genome size, although it was only a contig-level assembly lacking gene annotations. Here, we describe a reference-guided effort to scaffold those contigs into chromosome-length pseudomolecules, add in any missing sequence that was unique to the IWGSC CS v1.0 assembly, and annotate the resulting pseudomolecules with genes. Our updated assembly, Triticum_aestivum_4.0, contains 15.07 Gbp of nongap sequence anchored to chromosomes, which is 1.2 Gbps more than the previous reference assembly. It includes 108,639 genes unambiguously localized to chromosomes, including over 2000 genes that were previously unplaced. We also discovered >5700 additional gene copies, facilitating the accurate annotation of functional gene duplications including at the Ppd-B1 photoperiod response locus.

Genetic dissection of the auxin response network

The expansion of gene families during evolution, which can generate functional overlap or specialization among their members, is a characteristic feature of signalling pathways in complex organisms. For example, families of transcriptional activators and repressors mediate responses to the plant hormone auxin. Although these regulators were identified more than 20 years ago, their overlapping functions and compensating negative feedbacks have hampered their functional analyses. Studies using loss-of-function approaches in basal land plants and gain-of-function approaches in angiosperms have in part overcome these issues but have still left an incomplete understanding. Here, we propose that renewed emphasis on genetic analysis of multiple mutants and species will shed light on the role of gene families in auxin response. Combining loss-of-function mutations in auxin-response activators and repressors can unravel complex outputs enabled by expanded gene families, such as fine-tuned developmental outcomes and robustness. Similar approaches and concepts may help to analyse other regulatory pathways whose components are also encoded by large gene families.

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.

Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression

Melon exhibits substantial natural variation especially in fruit ripening physiology, including both climacteric (ethylene-producing) and non-climacteric types. However, genomic mechanisms underlying such variation are not yet fully understood. Here, we report an Oxford Nanopore-based high-grade genome reference in the semi-climacteric cultivar Harukei-3 (378 Mb + 33,829 protein-coding genes), with an update of tissue-wide RNA-seq atlas in the Melonet-DB database. Comparison between Harukei-3 and DHL92, the first published melon genome, enabled identification of 24,758 one-to-one orthologue gene pairs, whereas others were candidates of copy number variation or presence/absence polymorphisms (PAPs). Further comparison based on 10 melon genome assemblies identified genome-wide PAPs of 415 retrotransposon Gag-like sequences. Of these, 160 showed fruit ripening-inducible expression, with 59.4% of the neighboring genes showing similar expression patterns (r > 0.8). Our results suggest that retrotransposons contributed to the modification of gene expression during diversification of melon genomes, and may affect fruit ripening-inducible gene expression.

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.

A genome-wide survey of copy number variations reveals an asymmetric evolution of duplicated genes in rice

Background

Copy number variations (CNVs) are an important type of structural variations in the genome that usually affect gene expression levels by gene dosage effect. Understanding CNVs as part of genome evolution may provide insights into the genetic basis of important agricultural traits and contribute to the crop breeding in the future. While available methods to detect CNVs utilizing next-generation sequencing technology have helped shed light on prevalence and effects of CNVs, the complexity of crop genomes poses a major challenge and requires development of additional tools.

Results

Here, we generated genomic and transcriptomic data of 93 rice (Oryza sativa L.) accessions and developed a comprehensive pipeline to call CNVs in this large-scale dataset. We analyzed the correlation between CNVs and gene expression levels and found that approximately 13% of the identified genes showed a significant correlation between their expression levels and copy numbers. Further analysis showed that about 36% of duplicate pairs were involved in pseudogenetic events while only 5% of them showed functional differentiation. Moreover, the offspring copy mainly contributed to the expression levels and seemed more likely to become a pseudogene, whereas the parent copy tended to maintain the function of ancestral gene.

Conclusion

We provide a high-accuracy CNV dataset that will contribute to functional genomics studies and molecular breeding in rice. We also showed that gene dosage effect of CNVs in rice is not exponential or linear. Our work demonstrates that the evolution of duplicated genes is asymmetric in both expression levels and gene fates, shedding a new insight into the evolution of duplicated genes.

Dissection of complex traits of tomato in the post-genome era

We present the main advances of dissection of complex traits in tomato by omics, the genes identified to control complex traits and the application of CRISPR/Cas9 in tomato breeding. Complex traits are believed to be under the control of multiple genes, each with different effects and interaction with environmental factors. Advance development of sequencing and molecular technologies has enabled the recognition of the genomic structure of most organisms and the identification of a nearly limitless number of markers that have made it to accelerate the speed of QTL identification and gene cloning. Meanwhile, multiomics have been used to identify the genetic variations among different tomato species, determine the expression profiles of genes in different tissues and at distinct developmental stages, and detect metabolites in different pathways and processes. The combination of these data facilitates to reveal mechanism underlying complex traits. Moreover, mutants generated by mutagens and genome editing provide relatively rich genetic variation for deciphering the complex traits and exploiting them in tomato breeding. In this article, we present the main advances of complex trait dissection in tomato by omics since the release of the tomato genome sequence in 2012. We provide further insight into some tomato complex traits because of the causal genetic variations discovered so far and explore the utilization of CRISPR/Cas9 for the modification of tomato complex traits.