eQTL Regulating Transcript Levels Associated with Diverse Biological Processes in Tomato

A nuclear receptor HR96-related gene underlies large trans-driven differences in detoxification gene expression in a generalist herbivore

The role, magnitude, and molecular nature of trans-driven expression variation underlying the upregulation of detoxification genes in pesticide resistant arthropod populations has remained enigmatic. In this study, we performed expression quantitative trait locus (eQTL) mapping (n = 458) between a pesticide resistant and a susceptible strain of the generalist herbivore and crop pest Tetranychus urticae. We found that a single trans eQTL hotspot controlled large differences in the expression of a subset of genes in different detoxification gene families, as well as other genes associated with host plant use. As established by additional genetic approaches including RNAi gene knockdown, a duplicated gene with a nuclear hormone receptor HR96-related ligand-binding domain was identified as causal for the expression differences between strains. The presence of a large family of HR96-related genes in T. urticae may enable modular control of detoxification and host plant use genes, facilitating this species' known and rapid evolution to diverse pesticides and host plants.

Diversity of tomato leaf form provides novel insights into breeding

Tomato (Solanum lycopersicum L.) is cultivated widely globally. The crop exhibits tremendous morphological variations because of its long breeding history. Apart from the commercial tomato varieties, wild species and heirlooms are grown in certain regions of the world. Since the fruit constitutes the edible part, much of the agronomical research is focused on it. However, recent studies have indicated that leaf morphology influences fruit quality. As leaves are specialized photosynthetic organs and the vascular systems transport the photosynthetic products to sink organs, the architectural characteristics of the leaves have a strong influence on the final fruit quality. Therefore, comprehensive research focusing on both the fruit and leaf morphology is required for further tomato breeding. This review summarizes an overview of knowledge of the basic tomato leaf development, morphological diversification, and molecular mechanisms behind them and emphasizes its importance in breeding. Finally, we discuss how these findings and knowledge can be applied to future tomato breeding.

Plasticity of maternal environment-dependent expression-QTLs of tomato seeds

Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.

A comparative transcriptomics and eQTL approach identifies SlWD40 as a tomato fruit ripening regulator

Although multiple vital genes with strong effects on the tomato (Solanum lycopersicum) ripening process have been identified via the positional cloning of ripening mutants and cloning of ripening-related transcription factors (TFs), recent studies suggest that it is unlikely that we have fully characterized the gene regulatory networks underpinning this process. Here, combining comparative transcriptomics and expression QTLs, we identified 16 candidate genes involved in tomato fruit ripening and validated them through virus-induced gene silencing analysis. To further confirm the accuracy of the approach, one potential ripening regulator, SlWD40 (WD-40 repeats), was chosen for in-depth analysis. Co-expression network analysis indicated that master regulators such as RIN (ripening inhibitor) and NOR (nonripening) as well as vital TFs including FUL1 (FRUITFUL1), SlNAC4 (NAM, ATAF1,2, and CUC2 4), and AP2a (Activating enhancer binding Protein 2 alpha) strongly co-expressed with SlWD40. Furthermore, SlWD40 overexpression and RNAi lines exhibited substantially accelerated and delayed ripening phenotypes compared with the wild type, respectively. Moreover, transcriptome analysis of these transgenics revealed that expression patterns of ethylene biosynthesis genes, phytoene synthase, pectate lyase, and branched chain amino transferase 2, in SlWD40-RNAi lines were similar to those of rin and nor fruits, which further demonstrated that SlWD40 may act as an important ripening regulator in conjunction with RIN and NOR. These results are discussed in the context of current models of ripening and in terms of the use of comparative genomics and transcriptomics as an effective route for isolating causal genes underlying differences in genotypes.

eQTLs are key players in the integration of genomic and transcriptomic data for phenotype prediction

Background

Multi-omics represent a promising link between phenotypes and genome variation. Few studies yet address their integration to understand genetic architecture and improve predictability.

Results

Our study used 241 poplar genotypes, phenotyped in two common gardens, with xylem and cambium RNA sequenced at one site, yielding large phenotypic, genomic (SNP), and transcriptomic datasets. Prediction models for each trait were built separately for SNPs and transcripts, and compared to a third model integrated by concatenation of both omics. The advantage of integration varied across traits and, to understand such differences, an eQTL analysis was performed to characterize the interplay between the genome and transcriptome and classify the predicting features into cis or trans relationships. A strong, significant negative correlation was found between the change in predictability and the change in predictor ranking for trans eQTLs for traits evaluated in the site of transcriptomic sampling.

Conclusions

Consequently, beneficial integration happens when the redundancy of predictors is decreased, likely leaving the stage to other less prominent but complementary predictors. An additional gene ontology (GO) enrichment analysis appeared to corroborate such statistical output. To our knowledge, this is a novel finding delineating a promising method to explore data integration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08690-7.

Heinz-resistant tomato cultivars exhibit a lignin-based resistance to field dodder (Cuscuta campestris) parasitism

Cuscuta species (dodders) are agriculturally destructive, parasitic angiosperms. These parasitic plants use haustoria as physiological bridges to extract nutrients and water from hosts. Cuscuta campestris has a broad host range and wide geographical distribution. While some wild tomato relatives are resistant, cultivated tomatoes are generally susceptible to C. campestris infestations. However, some specific Heinz tomato (Solanum lycopersicum) hybrid cultivars exhibit resistance to dodders in the field, but their defense mechanism was previously unknown. Here, we discovered that the stem cortex in these resistant lines responds with local lignification upon C. campestris attachment, preventing parasite entry into the host. Lignin Induction Factor 1 (LIF1, an AP2-like transcription factor), SlMYB55, and Cuscuta R-gene for Lignin-based Resistance 1, a CC-NBS-LRR (CuRLR1) are identified as factors that confer host resistance by regulating lignification. SlWRKY16 is upregulated upon C. campestris infestation and potentially negatively regulates LIF1 function. Intriguingly, CuRLR1 may play a role in signaling or function as an intracellular receptor for receiving Cuscuta signals or effectors, thereby regulating lignification-based resistance. In summary, these four regulators control the lignin-based resistance response in specific Heinz tomato cultivars, preventing C. campestris from parasitizing resistant tomatoes. This discovery provides a foundation for investigating multilayer resistance against Cuscuta species and has potential for application in other essential crops attacked by parasitic plants.

Genetic variants associated mRNA stability in lung

Background

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs).

Results

Here, we presented a computational framework that takes advantage of recently developed methods to infer the mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3942 genes and 186,132 eQTLs for 4751 genes from 15,122,700 genetic variants for 13,476 genes on the autosomes, respectively. Interestingly, our results indicated that stQTLs were enriched in the CDS and 3’UTR regions, while eQTLs are enriched in the CDS, 3’UTR, 5’UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08405-y.

Multi-omics analysis reveals the genetic basis of rice fragrance mediated by betaine aldehyde dehydrogenase 2

INTRODUCTION

Fragrance is an important economic and quality trait in rice. The trait is controlled by the recessive gene betaine aldehyde dehydrogenase 2 (BADH2) via the production of 2-acetyl-1-pyrroline (2AP).

OBJECTIVES

Variation in BADH2 was evaluated at the population, genetic, transcriptional, and metabolic levels to obtain insights into fragrance regulation in rice.

METHODS

Whole-genome resequencing of the Korean World Rice Collection of 475 rice accessions, including 421 breeding lines and 54 wild accessions, was performed. Transcriptome analyses of a subset of 279 accessions, proteome analyses of 64 accessions, and volatile profiling of 421 breeding lines were also performed.

RESULTS

We identified over 3.1 million high-quality single nucleotide polymorphisms (SNPs) in Korean rice collection. Most SNPs were present in intergenic regions (79%), and 190,148 SNPs (6%) were located in the coding sequence, of which 53% were nonsynonymous. In total, 38 haplotypes were identified in the BADH2 coding region, including four novel haplotypes (one in cultivated and three in wild accessions). Tajima's D values suggested that BADH2 was under balancing selection in japonica rice. Furthermore, we identified 316 expression quantitative trait loci (eQTL), including 185 cis-eQTLs and 131 trans-eQTLs, involved in BADH2 regulation. A protein quantitative trait loci (pQTL) analysis revealed the presence of trans-pQTLs; 13 pQTLs were mapped 1 Mbp from the BADH2 region. Based on variable importance in projection (VIP) scores, 15 volatile compounds, including 2AP, discriminated haplotypes and were potential biomarkers for rice fragrance.

CONCLUSION

We generated a catalog of haplotypes based on a resequencing analysis of a large number of rice accessions. eQTLs and pQTLs associated with BADH2 gene expression and protein accumulation are likely involved in the regulation of 2AP variation in fragrant rice. These data improve our understanding of fragrance and provide valuable information for rice breeding.

The genetics of gene expression in a Caenorhabditis elegans multiparental recombinant inbred line population

Studying genetic variation of gene expression provides a powerful way to unravel the molecular components underlying complex traits. Expression quantitative trait locus (eQTL) studies have been performed in several different model species, yet most of these linkage studies have been based on the genetic segregation of two parental alleles. Recently, we developed a multiparental segregating population of 200 recombinant inbred lines (mpRILs) derived from four wild isolates (JU1511, JU1926, JU1931, and JU1941) in the nematode Caenorhabditis elegans. We used RNA-seq to investigate how multiple alleles affect gene expression in these mpRILs. We found 1789 genes differentially expressed between the parental lines. Transgression, expression beyond any of the parental lines in the mpRILs, was found for 7896 genes. For expression QTL mapping almost 9000 SNPs were available. By combining these SNPs and the RNA-seq profiles of the mpRILs, we detected almost 6800 eQTLs. Most trans-eQTLs (63%) co-locate in six newly identified trans-bands. The trans-eQTLs found in previous two-parental allele eQTL experiments and this study showed some overlap (17.5–46.8%), highlighting on the one hand that a large group of genes is affected by polymorphic regulators across populations and conditions, on the other hand, it shows that the mpRIL population allows identification of novel gene expression regulatory loci. Taken together, the analysis of our mpRIL population provides a more refined insight into C. elegans complex trait genetics and eQTLs in general, as well as a starting point to further test and develop advanced statistical models for detection of multiallelic eQTLs and systems genetics studying the genotype–phenotype relationship.

A combination of genome-wide and transcriptome-wide association studies reveals genetic elements leading to male sterility during high temperature stress in cotton

Global warming has reduced the productivity of many field-grown crops, as the effects of high temperatures can lead to male sterility in such plants. Genetic regulation of the high temperature (HT) response in the major crop cotton is poorly understood. We determined the functionality and transcriptomes of the anthers of 218 cotton accessions grown under HT stress. By analyzing transcriptome divergence and implementing a genome-wide association study (GWAS), we identified three thermal tolerance associated loci which contained 75 protein coding genes and 27 long noncoding RNAs, and provided expression quantitative trait loci (eQTLs) for 13 132 transcripts. A transcriptome-wide association study (TWAS) confirmed six causal elements for the HT response (three genes overlapped with the GWAS results) which are involved in protein kinase activity. The most susceptible gene, GhHRK1, was confirmed to be a previously uncharacterized negative regulator of the HT response in both cotton and Arabidopsis. These functional variants provide a new understanding of the genetic basis for HT tolerance in male reproductive organs.

Cis- and Trans-Regulatory Variations in the Domestication of the Chili Pepper Fruit

The process of domestication requires the rapid transformation of the wild morphology into the cultivated forms that humans select for. This process often takes place through changes in the regulation of genes, yet, there is no definite pattern on the role of cis- and trans-acting regulatory variations in the domestication of the fruit among crops. Using allele-specific expression and network analyses, we characterized the regulatory patterns and the inheritance of gene expression in wild and cultivated accessions of chili pepper, a crop with remarkable fruit morphological variation. We propose that gene expression differences associated to the cultivated form are best explained by cis-regulatory hubs acting through trans-regulatory cascades. We show that in cultivated chili, the expression of genes associated with fruit morphology is partially recessive with respect to those in the wild relative, consistent with the hybrid fruit phenotype. Decreased expression of fruit maturation and growth genes in cultivated chili suggest that selection for loss-of-function took place in its domestication. Trans-regulatory changes underlie the majority of the genes showing regulatory divergence and had larger effect sizes on gene expression than cis-regulatory variants. Network analysis of selected cis-regulated genes, including ARP9 and MED25, indicated their interaction with many transcription factors involved in organ growth and fruit ripening. Differentially expressed genes linked to cis-regulatory variants and their interactions with downstream trans-acting genes have the potential to drive the morphological differences observed between wild and cultivated fruits and provide an attractive mechanism of morphological transformation during the domestication of the chili pepper.

Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans

The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.

Leaf shape in Populus tremula is a complex, omnigenic trait

Leaf shape is a defining feature of how we recognize and classify plant species. Although there is extensive variation in leaf shape within many species, few studies have disentangled the underlying genetic architecture. We characterized the genetic architecture of leaf shape variation in Eurasian aspen (Populus tremula L.) by performing genome-wide association study (GWAS) for physiognomy traits. To ascertain the roles of identified GWAS candidate genes within the leaf development transcriptional program, we generated RNA-Seq data that we used to perform gene co-expression network analyses from a developmental series, which is publicly available within the PlantGenIE resource. We additionally used existing gene expression measurements across the population to analyze GWAS candidate genes in the context of a population-wide co-expression network and to identify genes that were differentially expressed between groups of individuals with contrasting leaf shapes. These data were integrated with expression GWAS (eQTL) results to define a set of candidate genes associated with leaf shape variation. Our results identified no clear adaptive link to leaf shape variation and indicate that leaf shape traits are genetically complex, likely determined by numerous small-effect variations in gene expression. Genes associated with shape variation were peripheral within the population-wide co-expression network, were not highly connected within the leaf development co-expression network, and exhibited signatures of relaxed selection. As such, our results are consistent with the omnigenic model.

eQTL mapping of the 12S globulin cruciferin gene PGCRURSE5 as a novel candidate associated with starch content in potato tubers

Tuber starch content (TSC) is a very important trait in potato (Solanum tuberosum L.). This study is the first to use expression quantitative trait loci (eQTL) mapping of transcript-derived markers for TSC in potato. Thirty-four differentially expressed genes were selected by comparing the RNA-seq data of contrasting bulked segregants. For the 11 candidate genes, we determined their relative expression levels across the segregating diploid potato population using RT-qPCR. We detected 36 eQTL as candidate genes distributed on all twelve potato chromosomes, and nine of them overlapped with QTL for TSC. Peaks for two eQTL, eAGPaseS-a and ePGRCRURSE5, were close to the corresponding loci of the large subunit of ADP-glucose pyrophosphorylase (AGPaseS-a) and the 12S globulin cruciferin gene (PGCRURSE5), respectively. The eQTL peaks for AGPaseS-a and PGRCRURSE5 explained 41.0 and 28.3% of the phenotypic variation at the transcript level. We showed the association of the DNA markers for AGPaseS-a and PGRCRURSE5 with QTL for TSC, and significant correlation between the expression level of PGRCRURSE5 and TSC. We did not observe a significant correlation between the expression level of AGPaseS-a and TSC. We concluded that the cruciferin gene PGRCRURSE5 is a novel candidate involved in the regulation of starch content in potato tubers.

Xylem systems genetics analysis reveals a key regulator of lignin biosynthesis in Populus deltoides

Despite the growing resources and tools for high-throughput characterization and analysis of genomic information, the discovery of the genetic elements that regulate complex traits remains a challenge. Systems genetics is an emerging field that aims to understand the flow of biological information that underlies complex traits from genotype to phenotype. In this study, we used a systems genetics approach to identify and evaluate regulators of the lignin biosynthesis pathway in Populus deltoides by combining genome, transcriptome, and phenotype data from a population of 268 unrelated individuals of P. deltoides The discovery of lignin regulators began with the quantitative genetic analysis of the xylem transcriptome and resulted in the detection of 6706 and 4628 significant local- and distant-eQTL associations, respectively. Among the locally regulated genes, we identified the R2R3-MYB transcription factor MYB125 (Potri.003G114100) as a putative trans-regulator of the majority of genes in the lignin biosynthesis pathway. The expression of MYB125 in a diverse population positively correlated with lignin content. Furthermore, overexpression of MYB125 in transgenic poplar resulted in increased lignin content, as well as altered expression of genes in the lignin biosynthesis pathway. Altogether, our findings indicate that MYB125 is involved in the control of a transcriptional coexpression network of lignin biosynthesis genes during secondary cell wall formation in P. deltoides.

High Photosynthetic Rates in a Solanum pennellii Chromosome 2 QTL Is Explained by Biochemical and Photochemical Changes

Enhanced photosynthesis is strictly associated with to productivity and it can be accomplished by genetic approaches through identification of genetic variation. By using a Solanum pennellii introgression lines (ILs) population, it was previously verified that, under normal (CO2), IL 2-5 and 2-6 display increased photosynthetic rates by up to 20% in comparison with their parental background (M82). However, the physiological mechanisms involved in the enhanced CO2 assimilation exhibited by these lines remained unknown, precluding their use for further biotechnological applications. Thereby, here we attempted to uncover the physiological factors involved in the upregulation of photosynthesis in ILs 2-5 and 2-6 under normal (CO2) as well as under elevated (CO2). The results provide evidence for increased biochemical capacity (higher maximum carboxylation velocity and maximum electron transport rate) in plants from IL 2-5 and 2-6, whereas the diffusive components (stomatal and mesophyll conductances) were unaltered in these ILs in comparison to M82. Our analyses revealed that the higher photosynthetic rate observed in these ILs was associated with higher levels of starch as well as total protein levels, specially increased RuBisCO content. Further analyses performed in plants under high (CO2) confirmed that biochemical properties are involved in genetic variation on chromosome 2 related to enhanced photosynthesis.

The Genetic Regulation of Alternative Splicing in Populus deltoides

Alternative splicing (AS) is a mechanism of regulation of the proteome via enabling the production of multiple mRNAs from a single gene. To date, the dynamics of AS and its effects on the protein sequences of individuals in a large and genetically unrelated population of trees have not been investigated. Here we describe the diversity of AS events within a previously genotyped population of 268 individuals of Populus deltoides and their putative downstream functional effects. Using a robust bioinformatics pipeline, the AS events and resulting transcript isoforms were discovered and quantified for each individual in the population. Analysis of the AS revealed that, as expected, most AS isoforms are conserved. However, we also identified a substantial collection of new, unannotated splice junctions and transcript isoforms. Heritability estimates for the expression of transcript isoforms showed that approximately half of the isoforms are heritable. The genetic regulators of these AS isoforms and splice junction usage were then identified using a genome-wide association analysis. The expression of AS isoforms was predominately cis regulated while splice junction usage was generally regulated in trans. Additionally, we identified 696 genes encoding alternatively spliced isoforms that changed putative protein domains relative to the longest protein coding isoform of the gene, and 859 genes exhibiting this same phenomenon relative to the most highly expressed isoform. Finally, we found that 748 genes gained or lost micro-RNA binding sites relative to the longest protein coding isoform of a given gene, while 940 gained or lost micro-RNA binding sites relative to the most highly expressed isoform. These results indicate that a significant fraction of AS events are genetically regulated and that this isoform usage can result in protein domain architecture changes.

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.

Genetic dissection of the gene coexpression network underlying photosynthesis in Populus

Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein-encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue-specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis-related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker-assisted breeding efforts to genetically improve woody plants.

Genome-wide analysis of expression quantitative trait loci (eQTLs) reveals the regulatory architecture of gene expression variation in the storage roots of sweet potato

Dissecting the genetic regulation of gene expression is critical for understanding phenotypic variation and species evolution. However, our understanding of the transcriptional variability in sweet potato remains limited. Here, we analyzed two publicly available datasets to explore the landscape of transcriptomic variations and its genetic basis in the storage roots of sweet potato. The comprehensive analysis identified a total of 724,438 high-confidence single nucleotide polymorphisms (SNPs) and 26,026 expressed genes. Expression quantitative trait locus (eQTL) analysis revealed 4408 eQTLs regulating the expression of 3646 genes, including 2261 local eQTLs and 2147 distant eQTLs. Two distant eQTL hotspots were found with target genes significantly enriched in specific functional classifications. By combining the information from regulatory network analyses, eQTLs and association mapping, we found that IbMYB1-2 acts as a master regulator and is the major gene responsible for the activation of anthocyanin biosynthesis in the storage roots of sweet potato. Our study provides the first insight into the genetic architecture of genome-wide expression variation in sweet potato and can be used to investigate the potential effects of genetic variants on key agronomic traits in sweet potato.

Trait discovery and editing in tomato

Tomato (Solanum lycopersicum), which is used for both processing and fresh markets, is a major crop species that is the top ranked vegetable produced over the world. Tomato is also a model species for research in genetics, fruit development and disease resistance. Genetic resources available in public repositories comprise the 12 wild related species and thousands of landraces, modern cultivars and mutants. In addition, high quality genome sequences are available for cultivated tomato and for several wild relatives, hundreds of accessions have been sequenced, and databases gathering sequence data together with genetic and phenotypic data are accessible to the tomato community. Major breeding goals are productivity, resistance to biotic and abiotic stresses, and fruit sensorial and nutritional quality. New traits, including resistance to various biotic and abiotic stresses and root architecture, are increasingly being studied. Several major mutations and quantitative trait loci (QTLs) underlying traits of interest in tomato have been uncovered to date and, thanks to new populations and advances in sequencing technologies, the pace of trait discovery has considerably accelerated. In recent years, clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing (GE) already proved its remarkable efficiency in tomato for engineering favorable alleles and for creating new genetic diversity by gene disruption, gene replacement, and precise base editing. Here, we provide insight into the major tomato traits and underlying causal genetic variations discovered so far and review the existing genetic resources and most recent strategies for trait discovery in tomato. Furthermore, we explore the opportunities offered by CRISPR/Cas9 and their exploitation for trait editing in tomato.

Regulation of Root Angle and Gravitropism

Regulation of plant root angle is critical for obtaining nutrients and water and is an important trait for plant breeding. A plant’s final, long-term root angle is the net result of a complex series of decisions made by a root tip in response to changes in nutrient availability, impediments, the gravity vector and other stimuli. When a root tip is displaced from the gravity vector, the short-term process of gravitropism results in rapid reorientation of the root toward the vertical. Here, we explore both short- and long-term regulation of root growth angle, using natural variation in tomato to identify shared and separate genetic features of the two responses. Mapping of expression quantitative trait loci mapping and leveraging natural variation between and within species including Arabidopsis suggest a role for PURPLE ACID PHOSPHATASE 27 and CELL DIVISION CYCLE 73 in determining root angle.

Network Analysis Reveals a Role for Salicylic Acid Pathway Components in Shade Avoidance

Plants have sophisticated mechanisms for sensing neighbor shade. To maximize their ability to compete for light, plants respond to shade through enhanced elongation and physiological changes. The shade avoidance response affects many different organs and growth stages, yet the signaling pathways underlying this response have mostly been studied in seedlings. We assayed transcriptome changes in response to shade across a 2-d time course in the wild type and 12 Arabidopsis (Arabidopsis thaliana) mutants. The resulting temporal map of transcriptional responses to shade defines early and late responses in adult plants, enabling us to determine connections between key signaling genes and downstream responses. We found a pervasive and unexpectedly strong connection between shade avoidance and genes related to salicylic acid, suggesting salicylic acid signaling to be an important shade avoidance growth regulator. We tested this connection and found that several mutants disrupting salicylic acid levels or signaling were defective in shade avoidance. The effect of these mutations on shade avoidance was specific to petiole elongation; neither hypocotyl nor flowering time responses were altered, thereby defining important stage-specific differences in the downstream shade avoidance signaling pathway. Shade treatment did not change salicylic acid levels, indicating that the mediation of shade avoidance by salicylic acid is not dependent on the modulation of salicylic acid levels. These results demonstrate that salicylic acid pathway genes also are key components of petiole shade avoidance.

Allele-specific expression and genetic determinants of transcriptomic variations in response to mild water deficit in tomato

Characterizing the natural diversity of gene expression across environments is an important step in understanding how genotype-by-environment interactions shape phenotypes. Here, we analyzed the impact of water deficit onto gene expression levels in tomato at the genome-wide scale. We sequenced the transcriptome of growing leaves and fruit pericarps at cell expansion stage in a cherry and a large fruited accession and their F₁ hybrid grown under two watering regimes. Gene expression levels were steadily affected by the genotype and the watering regime. Whereas phenotypes showed mostly additive inheritance, ~80% of the genes displayed non-additive inheritance. By comparing allele-specific expression (ASE) in the F₁ hybrid to the allelic expression in both parental lines, respectively, 3005 genes in leaf and 2857 genes in fruit deviated from 1:1 ratio independently of the watering regime. Among these genes, ~55% were controlled by cis factors, ~25% by trans factors and ~20% by a combination of both types of factors. A total of 328 genes in leaf and 113 in fruit exhibited significant ASE-by-watering regime interaction, among which ~80% presented trans-by-watering regime interaction, suggesting a response to water deficit mediated through a majority of trans-acting loci in tomato. We cross-validated the expression levels of 274 transcripts in fruit and leaves of 124 recombinant inbred lines (RILs) and identified 163 expression quantitative trait loci (eQTLs) mostly confirming the divergences identified by ASE. Combining phenotypic and expression data, we observed a complex network of variation between genes encoding enzymes involved in the sugar metabolism.

Quantitative Trait Loci Analysis Identifies a Prominent Gene Involved in the Production of Fatty Acid-Derived Flavor Volatiles in Tomato

To gain insight into the genetic regulation of lipid metabolism in tomato, we conducted metabolic trait loci (mQTL) analysis following the lipidomic profiling of fruit pericarp and leaf tissue of the Solanum pennellii introgression lines (IL). To enhance mapping resolution for selected fruit-specific mQTL, we profiled the lipids in a subset of independently derived S. pennellii backcross inbred lines, as well as in a near-isogenic sub-IL population. We identified a putative lecithin:cholesterol acyltransferase that controls the levels of several lipids, and two members of the class III lipase family, LIP1 and LIP2, that were associated with decreased levels of diacylglycerols (DAGs) and triacylglycerols (TAGs). Lipases of this class cleave fatty acids from the glycerol backbone of acylglycerols. The released fatty acids serve as precursors of flavor volatiles. We show that LIP1 expression correlates with fatty acid-derived volatile levels. We further confirm the function of LIP1 in TAG and DAG breakdown and volatile synthesis using transgenic plants. Taken together, our study extensively characterized the genetic architecture of lipophilic compounds in tomato and demonstrated at molecular level that release of free fatty acids from the glycerol backbone can have a major impact on downstream volatile synthesis.

Genetic mapping identifies loci that influence tomato resistance against Colorado potato beetles

The Colorado potato beetle (CPB; Leptinotarsa decemlineata Say), the most economically important insect pest on potato (Solanum tuberosum L.), also feeds on other Solanaceae, including cultivated tomato (Solanum lycopersicum L.). We used tomato genetic mapping populations to investigate natural variation in CPB resistance. CPB bioassays with 74 tomato lines carrying introgressions of Solanum pennellii in S. lycopersicum cv. M82 identified introgressions from S. pennellii on chromosomes 1 and 6 conferring CPB susceptibility, whereas introgressions on chromosomes 1, 8 and 10 conferred higher resistance. Mapping of CPB resistance using 113 recombinant inbred lines derived from a cross between S. lycopersicum cv UC-204B and Solanum galapagense identified significant quantitative trait loci on chromosomes 6 and 8. In each case, the S. galapagense alleles were associated with lower leaf damage and reduced larval growth. Results of both genetic mapping approaches converged on the same region of chromosome 6, which may have important functions in tomato defense against CPB herbivory. Although genetic mapping identified quantitative trait loci encompassing known genes for tomato acyl sugar and glycoalkaloid biosynthesis, experiments with acyl sugar near-isogenic lines and transgenic GAME9 glycoalkaloid-deficient and overproducing lines showed no significant effect of these otherwise insect-defensive metabolites on CPB performance.

Linkage mapping of yeast cross protection connects gene expression variation to a higher-order organismal trait

Gene expression variation is extensive in nature, and is hypothesized to play a major role in shaping phenotypic diversity. However, connecting differences in gene expression across individuals to higher-order organismal traits is not trivial. In many cases, gene expression variation may be evolutionarily neutral, and in other cases expression variation may only affect phenotype under specific conditions. To understand connections between gene expression variation and stress defense phenotypes, we have been leveraging extensive natural variation in the gene expression response to acute ethanol in laboratory and wild Saccharomyces cerevisiae strains. Previous work found that the genetic architecture underlying these expression differences included dozens of “hotspot” loci that affected many transcripts in trans. In the present study, we provide new evidence that one of these expression QTL hotspot loci affects natural variation in one particular stress defense phenotype—ethanol-induced cross protection against severe doses of H₂O₂. A major causative polymorphism is in the heme-activated transcription factor Hap1p, which we show directly impacts cross protection, but not the basal H₂O₂ resistance of unstressed cells. This provides further support that distinct cellular mechanisms underlie basal and acquired stress resistance. We also show that Hap1p-dependent cross protection relies on novel regulation of cytosolic catalase T (Ctt1p) during ethanol stress in a wild oak strain. Because ethanol accumulation precedes aerobic respiration and accompanying reactive oxygen species formation, wild strains with the ability to anticipate impending oxidative stress would likely be at an advantage. This study highlights how strategically chosen traits that better correlate with gene expression changes can improve our power to identify novel connections between gene expression variation and higher-order organismal phenotypes.

A major goal in genetics is to understand how individuals with different genetic makeups respond to their environment. Understanding these “gene-environment interactions” is important for the development of personalized medicine. For example, gene-environment interactions can explain why some people are more sensitive to certain drugs or are more likely to get certain cancers. While the underlying causes of gene-environment interactions are unclear, one possibility is that differences in gene expression across individuals are responsible. In this study, we examined that possibility using baker’s yeast as a model. We were interested in a phenomenon called acquired stress resistance, where cells exposed to a mild dose of one stress can become resistant to an otherwise lethal dose of severe stress. This response is observed in diverse organisms ranging from bacteria to humans, though the specific mechanisms governing acquisition of higher stress resistance are poorly understood. To understand the differences between yeast strains with and without the ability to acquire further stress resistance, we employed genetic mapping. We found that part of the variation in acquired stress resistance was due to sequence differences in a key regulatory protein, thus providing new insight into how different individuals respond to acute environmental change.

Verticillium dahliae Disease Resistance and the Regulatory Pathway for Maturity and Tuberization in Potato

Kleb. is a pathogenic fungus causing wilting, chlorosis, and early dying in potato ( L.). Genetic mapping of resistance to was done using a diploid population of potato. The major quantitative trait locus (QTL) for resistance was found on chromosome 5. The gene, controlling earliness of maturity and tuberization, was mapped within the interval. Another QTL on chromosome 9 co-localized with the wilt resistance gene marker. Epistasis analysis indicated that the loci on chromosomes 5 and 9 had a highly significant interaction, and that functioned downstream of The alleles were sequenced and found to encode StCDF1.1 and StCDF1.3. Interaction between the resistance allele and the was demonstrated, but not for Genome-wide expression QTL (eQTL) analysis was performed and genes with eQTL at the and loci were both found to have similar functions involving the chloroplast, including photosynthesis, which declines in both maturity and wilt. Among the gene ontology (GO) terms that were specific to genes with eQTL at the , but not the locus, were those associated with fungal defense. These results suggest that controls fungal defense and reduces early dying in wilt through affecting genetic pathway controlling tuberization timing.

Integrated QTL and eQTL Mapping Provides Insights and Candidate Genes for Fatty Acid Composition, Flowering Time, and Growth Traits in a F2 Population of a Novel Synthetic Allopolyploid Brassica napus

Brassica napus (B. napus, AACC), is an economically important allotetraploid crop species that resulted from hybridization between two diploid species, Brassica rapa (AA) and Brassica olereacea (CC). We have created one new synthetic B. napus genotype Da-Ae (AACC) and one introgression line Da-Ol-1 (AACC), which were used to generate an F₂ mapping population. Plants in this F₂ mapping population varied in fatty acid content, flowering time, and growth-related traits. Using quantitative trait locus (QTL) mapping, we aimed to determine if Da-Ae and Da-Ol-1 provided novel genetic variation beyond what has already been found in B. napus. Making use of the genotyping information generated from RNA-seq data of these two lines and their F₂ mapping population of 166 plants, we constructed a genetic map consisting of 2,021 single nucleotide polymorphism markers that spans 2,929 cM across 19 linkage groups. Besides the known major QTL identified, our high resolution genetic map facilitated the identification of several new QTL contributing to the different fatty acid levels, flowering time, and growth-related trait values. These new QTL probably represent novel genetic variation that existed in our new synthetic B. napus strain. By conducting genome-wide expression variation analysis in our F₂ mapping population, genetic regions that potentially regulate many genes across the genome were revealed. A FLOWERING LOCUS C gene homolog, which was identified as a candidate regulating flowering time and multiple growth-related traits, was found underlying one of these regions. Integrated QTL and expression QTL analyses also helped us identified candidate causative genes associated with various biological traits through expression level change and/or possible protein function modification.

Construction of a High-Density Genetic Map from RNA-Seq Data for an Arabidopsis Bay-0 × Shahdara RIL Population

High-density genetic maps are essential for high resolution mapping of quantitative traits. Here, we present a new genetic map for an Arabidopsis Bayreuth × Shahdara recombinant inbred line (RIL) population, built on RNA-seq data. RNA-seq analysis on 160 RILs of this population identified 30,049 single-nucleotide polymorphisms (SNPs) covering the whole genome. Based on a 100-kbp window SNP binning method, 1059 bin-markers were identified, physically anchored on the genome. The total length of the RNA-seq genetic map spans 471.70 centimorgans (cM) with an average marker distance of 0.45 cM and a maximum marker distance of 4.81 cM. This high resolution genotyping revealed new recombination breakpoints in the population. To highlight the advantages of such high-density map, we compared it to two publicly available genetic maps for the same population, comprising 69 PCR-based markers and 497 gene expression markers derived from microarray data, respectively. In this study, we show that SNP markers can effectively be derived from RNA-seq data. The new RNA-seq map closes many existing gaps in marker coverage, saturating the previously available genetic maps. Quantitative trait locus (QTL) analysis for published phenotypes using the available genetic maps showed increased QTL mapping resolution and reduced QTL confidence interval using the RNA-seq map. The new high-density map is a valuable resource that facilitates the identification of candidate genes and map-based cloning approaches.

Time-course expression QTL-atlas of the global transcriptional response of wheat to Fusarium graminearum

Fusarium head blight is a devastating disease of small grain cereals such as bread wheat (Triticum aestivum). The pathogen switches from a biotrophic to a nectrotrophic lifestyle in course of disease development forcing its host to adapt its defence strategies. Using a genetical genomics approach, we illustrate genome-wide reconfigurations of genetic control over transcript abundances between two decisive time points after inoculation with the causative pathogen Fusarium graminearum. Whole transcriptome measurements have been recorded for 163 lines of a wheat doubled haploid population segregating for several resistance genes yielding 15 552 at 30 h and 15 888 eQTL at 50 h after inoculation. The genetic map saturated with transcript abundance-derived markers identified of a novel QTL on chromosome 6A, besides the previously reported QTL Fhb1 and Qfhs.ifa-5A. We find a highly different distribution of eQTL between time points with about 40% of eQTL being unique for the respective assessed time points. But also for more than 20% of genes governed by eQTL at either time point, genetic control changes in time. These changes are reflected in the dynamic compositions of three major regulatory hotspots on chromosomes 2B, 4A and 5A. In particular, control of defence-related biological mechanisms concentrated in the hotspot at 4A shift to hotspot 2B as the disease progresses. Hotspots do not colocalize with phenotypic QTL, and within their intervals no higher than expected number of eQTL was detected. Thus, resistance conferred by either QTL is mediated by few or single genes.

De Novo Assembly of a New Solanum pennellii Accession Using Nanopore Sequencing

Updates in nanopore technology have made it possible to obtain gigabases of sequence data. Prior to this, nanopore sequencing technology was mainly used to analyze microbial samples. Here, we describe the generation of a comprehensive nanopore sequencing data set with a median read length of 11,979 bp for a self-compatible accession of the wild tomato species Solanum pennellii We describe the assembly of its genome to a contig N50 of 2.5 MB. The assembly pipeline comprised initial read correction with Canu and assembly with SMARTdenovo. The resulting raw nanopore-based de novo genome is structurally highly similar to that of the reference S. pennellii LA716 accession but has a high error rate and was rich in homopolymer deletions. After polishing the assembly with Illumina reads, we obtained an error rate of <0.02% when assessed versus the same Illumina data. We obtained a gene completeness of 96.53%, slightly surpassing that of the reference S. pennellii Taken together, our data indicate that such long read sequencing data can be used to affordably sequence and assemble gigabase-sized plant genomes.

Genetic Architecture and Molecular Networks Underlying Leaf Thickness in Desert-Adapted Tomato Solanum pennellii

Thicker leaves allow plants to grow in water-limited conditions. However, our understanding of the genetic underpinnings of this highly functional leaf shape trait is poor. We used a custom-built confocal profilometer to directly measure leaf thickness in a set of introgression lines (ILs) derived from the desert tomato Solanum pennellii and identified quantitative trait loci. We report evidence of a complex genetic architecture of this trait and roles for both genetic and environmental factors. Several ILs with thick leaves have dramatically elongated palisade mesophyll cells and, in some cases, increased leaf ploidy. We characterized the thick IL2-5 and IL4-3 in detail and found increased mesophyll cell size and leaf ploidy levels, suggesting that endoreduplication underpins leaf thickness in tomato. Next, we queried the transcriptomes and inferred dynamic Bayesian networks of gene expression across early leaf ontogeny in these lines to compare the molecular networks that pattern leaf thickness. We show that thick ILs share S. pennellii-like expression profiles for putative regulators of cell shape and meristem determinacy as well as a general signature of cell cycle-related gene expression. However, our network data suggest that leaf thickness in these two lines is patterned at least partially by distinct mechanisms. Consistent with this hypothesis, double homozygote lines combining introgression segments from these two ILs show additive phenotypes, including thick leaves, higher ploidy levels, and larger palisade mesophyll cells. Collectively, these data establish a framework of genetic, anatomical, and molecular mechanisms that pattern leaf thickness in desert-adapted tomato.

eQTLs Regulating Transcript Variations Associated with Rapid Internode Elongation in Deepwater Rice

To avoid low oxygen, oxygen deficiency or oxygen deprivation, deepwater rice cultivated in flood planes can develop elongated internodes in response to submergence. Knowledge of the gene regulatory networks underlying rapid internode elongation is important for an understanding of the evolution and adaptation of major crops in response to flooding. To elucidate the genetic and molecular basis controlling their deepwater response we used microarrays and performed expression quantitative trait loci (eQTL) and phenotypic QTL (phQTL) analyses of internode samples of 85 recombinant inbred line (RIL) populations of non-deepwater (Taichung 65)- and deepwater rice (Bhadua). After evaluating the phenotypic response of the RILs exposed to submergence, confirming the genotypes of the populations, and generating 188 genetic markers, we identified 10,047 significant eQTLs comprised of 2,902 cis-eQTLs and 7,145 trans-eQTLs and three significant eQTL hotspots on chromosomes 1, 4, and 12 that affect the expression of many genes. The hotspots on chromosomes 1 and 4 located at different position from phQTLs detected in this study and other previous studies. We then regarded the eQTL hotspots as key regulatory points to infer causal regulatory networks of deepwater response including rapid internode elongation. Our results suggest that the downstream regulation of the eQTL hotspots on chromosomes 1 and 4 is independent, and that the target genes are partially regulated by SNORKEL1 and SNORKEL2 genes (SK1/2), key ethylene response factors. Subsequent bioinformatic analyses, including gene ontology-based annotation and functional enrichment analysis and promoter enrichment analysis, contribute to enhance our understanding of SK1/2-dependent and independent pathways. One remarkable observation is that the functional categories related to photosynthesis and light signaling are significantly over-represented in the candidate target genes of SK1/2. The combined results of these investigations together with genetical genomics approaches using structured populations with a deepwater response are also discussed in the context of current molecular models concerning the rapid internode elongation in deepwater rice. This study provides new insights into the underlying genetic architecture of gene expression regulating the response to flooding in deepwater rice and will be an important community resource for analyses on the genetic basis of deepwater responses.

A New Advanced Backcross Tomato Population Enables High Resolution Leaf QTL Mapping and Gene Identification

Quantitative Trait Loci (QTL) mapping is a powerful technique for dissecting the genetic basis of traits and species differences. Established tomato mapping populations between domesticated tomato (Solanum lycopersicum) and its more distant interfertile relatives typically follow a near isogenic line (NIL) design, such as the S. pennellii Introgression Line (IL) population, with a single wild introgression per line in an otherwise domesticated genetic background. Here, we report on a new advanced backcross QTL mapping resource for tomato, derived from a cross between the M82 tomato cultivar and S. pennellii This so-called Backcrossed Inbred Line (BIL) population is comprised of a mix of BC2 and BC3 lines, with domesticated tomato as the recurrent parent. The BIL population is complementary to the existing S. pennellii IL population, with which it shares parents. Using the BILs, we mapped traits for leaf complexity, leaflet shape, and flowering time. We demonstrate the utility of the BILs for fine-mapping QTL, particularly QTL initially mapped in the ILs, by fine-mapping several QTL to single or few candidate genes. Moreover, we confirm the value of a backcrossed population with multiple introgressions per line, such as the BILs, for epistatic QTL mapping. Our work was further enabled by the development of our own statistical inference and visualization tools, namely a heterogeneous hidden Markov model for genotyping the lines, and by using state-of-the-art sparse regression techniques for QTL mapping.