A complete chromosome substitution mapping panel reveals genome-wide epistasis in Arabidopsis.

Chromosome substitution lines (CSLs) are tentatively supreme resources to investigate non-allelic genetic interactions. However, the difficulty of generating such lines in most species largely yielded imperfect CSL panels, prohibiting a systematic dissection of epistasis. Here, we present the development and use of a unique and complete panel of CSLs in Arabidopsis thaliana, allowing the full factorial analysis of epistatic interactions. A first comparison of reciprocal single chromosome substitutions revealed a dependency of QTL detection on different genetic backgrounds. The subsequent analysis of the complete panel of CSLs enabled the mapping of the genetic interactors and identified multiple two- and three-way interactions for different traits. Some of the detected epistatic effects were as large as any observed main effect, illustrating the impact of epistasis on quantitative trait variation. We, therefore, have demonstrated the high power of detection and mapping of genome-wide epistasis, confirming the assumed supremacy of comprehensive CSL sets.

Enhanced recombination empowers the detection and mapping of Quantitative Trait Loci.

Modern plant breeding, such as genomic selection and gene editing, is based on the knowledge of the genetic architecture of desired traits. Quantitative trait loci (QTL) analysis, which combines high throughput phenotyping and genotyping of segregating populations, is a powerful tool to identify these genetic determinants and to decipher the underlying mechanisms. However, meiotic recombination, which shuffles genetic information between generations, is limited: Typically only one to two exchange points, called crossovers, occur between a pair of homologous chromosomes. Here we test the effect on QTL analysis of boosting recombination, by mutating the anti-crossover factors RECQ4 and FIGL1 in Arabidopsis thaliana full hybrids and lines in which a single chromosome is hybrid. We show that increasing recombination ~6-fold empowers the detection and resolution of QTLs, reaching the gene scale with only a few hundred plants. Further, enhanced recombination unmasks some secondary QTLs undetected under normal recombination. These results show the benefits of enhanced recombination to decipher the genetic bases of traits.

Genetic Mapping of Genotype-by-Ploidy Effects in Arabidopsis thaliana.

Plants can express different phenotypic responses following polyploidization, but ploidy-dependent phenotypic variation has so far not been assigned to specific genetic factors. To map such effects, segregating populations at different ploidy levels are required. The availability of an efficient haploid inducer line in Arabidopsis thaliana allows for the rapid development of large populations of segregating haploid offspring. Because Arabidopsis haploids can be self-fertilised to give rise to homozygous doubled haploids, the same genotypes can be phenotyped at both the haploid and diploid ploidy level. Here, we compared the phenotypes of recombinant haploid and diploid offspring derived from a cross between two late flowering accessions to map genotype × ploidy (G × P) interactions. Ploidy-specific quantitative trait loci (QTLs) were detected at both ploidy levels. This implies that mapping power will increase when phenotypic measurements of monoploids are included in QTL analyses. A multi-trait analysis further revealed pleiotropic effects for a number of the ploidy-specific QTLs as well as opposite effects at different ploidy levels for general QTLs. Taken together, we provide evidence of genetic variation between different Arabidopsis accessions being causal for dissimilarities in phenotypic responses to altered ploidy levels, revealing a G × P effect. Additionally, by investigating a population derived from late flowering accessions, we revealed a major vernalisation-specific QTL for variation in flowering time, countering the historical bias of research in early flowering accessions.

Drought response in Arabidopsis displays synergistic coordination between stems and leaves.

The synergy between drought-responsive traits across different organs is crucial in the whole-plant mechanism influencing drought resilience. These organ interactions, however, are poorly understood, limiting our understanding of drought response strategies at the whole-plant level. Therefore, we need more integrative studies, especially on herbaceous species that represent many important food crops but remain underexplored in their drought response. We investigated inflorescence stems and rosette leaves of six Arabidopsis thaliana genotypes with contrasting drought tolerance, and combined anatomical observations with hydraulic measurements and gene expression studies to assess differences in drought response. The soc1ful double mutant was the most drought-tolerant genotype based on its synergistic combination of low stomatal conductance, largest stomatal safety margin, more stable leaf water potential during non-watering, reduced transcript levels of drought stress marker genes, and reduced loss of chlorophyll content in leaves, in combination with stems showing the highest embolism resistance, most pronounced lignification, and thickest intervessel pit membranes. In contrast, the most sensitive Cvi ecotype shows the opposite extreme of the same set of traits. The remaining four genotypes show variations in this drought syndrome. Our results reveal that anatomical, ecophysiological, and molecular adaptations across organs are intertwined, and multiple (differentially combined) strategies can be applied to acquire a certain level of drought tolerance.

Genotypic and tissue-specific variation of Populus nigra transcriptome profiles in response to drought.

Climate change is one of the most important challenges for mankind in the far and near future. In this regard, sustainable production of woody crops on marginal land with low water availability is a major challenge to tackle. This dataset is part of an experiment, in which we exposed three genetically differentiated genotypes of Populus nigra originating from contrasting natural habitats to gradually increasing moderate drought. RNA sequencing was performed on fine roots, developing xylem and leaves of those three genotypes under control and moderate drought conditions in order to get a comprehensive dataset on the transcriptional changes at the whole plant level under water limiting conditions. This dataset has already provided insight in the transcriptional control of saccharification potential of the three Populus genotypes under drought conditions and we suggest that our data will be valuable for further in-depth analysis regarding candidate gene identification or, on a bigger scale, for meta-transcriptome analysis.

Intervessel pit membrane thickness best explains variation in embolism resistance amongst stems of Arabidopsis thaliana accessions.

BACKGROUND AND AIMS: The ability to avoid drought-induced embolisms in the xylem is one of the essential traits for plants to survive periods of water shortage. Over the past three decades, hydraulic studies have been focusing on trees, which limits our ability to understand how herbs tolerate drought. Here we investigate the embolism resistance in inflorescence stems of four Arabidopsis thaliana accessions that differ in growth form and drought response. We assess functional traits underlying the variation in embolism resistance amongst the accessions studied using detailed anatomical observations. METHODS: Vulnerability to xylem embolism was evaluated via vulnerability curves using the centrifuge technique and linked with detailed anatomical observations in stems using light microscopy and transmission electron microscopy. KEY RESULTS: The data show significant differences in stem P50, varying 2-fold from -1.58 MPa in the Cape Verde Island accession to -3.07 MPa in the woody soc1 ful double mutant. Out of all the anatomical traits measured, intervessel pit membrane thickness (TPM) best explains the differences in P50, as well as P12 and P88. The association between embolism resistance and TPM can be functionally explained by the air-seeding hypothesis. There is no evidence that the correlation between increased woodiness and increased embolism resistance is directly related to functional aspects. However, we found that increased woodiness is strongly linked to other lignification characters, explaining why mechanical stem reinforcement is indirectly related to increased embolism resistance. CONCLUSIONS: The woodier or more lignified accessions are more resistant to embolism than the herbaceous accessions, confirming the link between increased stem lignification and increased embolism resistance, as also observed in other lineages. Intervessel pit membrane thickness and, to a lesser extent, theoretical vessel implosion resistance and vessel wall thickness are the missing functional links between stem lignification and embolism resistance.

Meiotic crossover reduction by virus-induced gene silencing enables the efficient generation of chromosome substitution lines and reverse breeding in Arabidopsis thaliana.

Plant breeding applications exploiting meiotic mutant phenotypes (like the increase or decrease of crossover (CO) recombination) have been proposed over the last years. As recessive meiotic mutations in breeding lines may affect fertility or have other pleiotropic effects, transient silencing techniques may be preferred. Reverse breeding is a breeding technique that would benefit from the transient downregulation of CO formation. The technique is essentially the opposite of plant hybridization: a method to extract parental lines from a hybrid. The method can also be used to efficiently generate chromosome substitution lines (CSLs). For successful reverse breeding, the two homologous chromosome sets of a heterozygous plant must be divided over two haploid complements, which can be achieved by the suppression of meiotic CO recombination and the subsequent production of doubled haploid plants. Here we show the feasibility of transiently reducing CO formation using virus-induced gene silencing (VIGS) by targeting the meiotic gene MSH5 in a wild-type heterozygote of Arabidopsis thaliana. The application of VIGS (rather than using lengthy stable transformation) generates transgene-free offspring with the desired genetic composition: we obtained parental lines from a wild-type heterozygous F1 in two generations. In addition, we obtained 20 (of the 32 possible) CSLs in one experiment. Our results demonstrate that meiosis can be modulated at will in A. thaliana to generate CSLs and parental lines rapidly for hybrid breeding. Furthermore, we illustrate how the modification of meiosis using VIGS can open routes to develop efficient plant breeding strategies.

The Role of Transcriptional Regulation in Hybrid Vigor.

The genetic basis of hybrid vigor in plants remains largely unsolved but strong evidence suggests that variation in transcriptional regulation can explain many aspects of this phenomenon. Natural variation in transcriptional regulation is highly abundant in virtually all species and thus a potential source of heterotic variability. Allele Specific Expression (ASE), which is tightly linked to parent of origin effects and modulated by complex interactions in cis and in trans, is generally considered to play a key role in explaining the differences between hybrids and parental lines. Here we discuss the recent developments in elucidating the role of transcriptional variation in a number of aspects of hybrid vigor, thereby bridging old paradigms and hypotheses with contemporary research in various species.

Reciprocal cybrids reveal how organellar genomes affect plant phenotypes.

Assessment of the impact of variation in chloroplast and mitochondrial DNA (collectively termed the plasmotype) on plant phenotypes is challenging due to the difficulty in separating their effect from nuclear-derived variation (the nucleotype). Haploid-inducer lines can be used as efficient plasmotype donors to generate new plasmotype-nucleotype combinations (cybrids)1. We generated a panel comprising all possible cybrids of seven Arabidopsis thaliana accessions and extensively phenotyped these lines for 1,859 phenotypes under both stable and fluctuating conditions. We show that natural variation in the plasmotype results in both additive and epistatic effects across all phenotypic categories. Plasmotypes that induce more additive phenotypic changes also cause more epistatic effects, suggesting a possible common basis for both additive and epistatic effects. On average, epistatic interactions explained twice as much of the variance in phenotypes as additive plasmotype effects. The impact of plasmotypic variation was also more pronounced under fluctuating and stressful environmental conditions. Thus, the phenotypic impact of variation in plasmotypes is the outcome of multi-level nucleotype-plasmotype-environment interactions and, as such, the plasmotype is likely to serve as a reservoir of variation that is predominantly exposed under certain conditions. The production of cybrids using haploid inducers is a rapid and precise method for assessment of the phenotypic effects of natural variation in organellar genomes. It will facilitate efficient screening of unique nucleotype-plasmotype combinations to both improve our understanding of natural variation in these combinations and identify favourable combinations to enhance plant performance.

Accelerating Climate Resilient Plant Breeding by Applying Next-Generation Artificial Intelligence.

Breeding crops for high yield and superior adaptability to new and variable climates is imperative to ensure continued food security, biomass production, and ecosystem services. Advances in genomics and phenomics are delivering insights into the complex biological mechanisms that underlie plant functions in response to environmental perturbations. However, linking genotype to phenotype remains a huge challenge and is hampering the optimal application of high-throughput genomics and phenomics to advanced breeding. Critical to success is the need to assimilate large amounts of data into biologically meaningful interpretations. Here, we present the current state of genomics and field phenomics, explore emerging approaches and challenges for multiomics big data integration by means of next-generation (Next-Gen) artificial intelligence (AI), and propose a workable path to improvement.

A genetical metabolomics approach for bioprospecting plant biosynthetic gene clusters.

OBJECTIVE: Plants produce a plethora of specialized metabolites to defend themselves against pathogens and insects, to attract pollinators and to communicate with other organisms. Many of these are also applied in the clinic and in agriculture. Genes encoding the enzymes that drive the biosynthesis of these metabolites are sometimes physically grouped on the chromosome, in regions called biosynthetic gene clusters (BGCs). Several algorithms have been developed to identify plant BGCs, but a large percentage of predicted gene clusters upon further inspection do not show coexpression or do not encode a single functional biosynthetic pathway. Hence, further prioritization is needed. RESULTS: Here, we introduce a strategy to systematically evaluate potential functions of predicted BGCs by superimposing their locations on metabolite quantitative trait loci (mQTLs). We show the feasibility of such an approach by integrating automated BGC prediction with mQTL datasets originating from a recombinant inbred line (RIL) population of Oryza sativa and a genome-wide association study (GWAS) of Arabidopsis thaliana. In these data, we identified several links for which the enzyme content of the BGCs matches well with the chemical features observed in the metabolite structure, suggesting that this method can effectively guide bioprospecting of plant BGCs.

Publisher Correction: A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms.

In the originally published Supplementary Information for this paper, the files presented as Supplementary Tables 3, 4, and 7 were duplicates of Supplementary Tables 5, 6, and 9, respectively. All Supplementary Table files are now correct online.

A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms.

We report a computational approach (implemented in MS-DIAL 3.0; http://prime.psc.riken.jp/) for metabolite structure characterization using fully 13C-labeled and non-labeled plants and LC-MS/MS. Our approach facilitates carbon number determination and metabolite classification for unknown molecules. Applying our method to 31 tissues from 12 plant species, we assigned 1,092 structures and 344 formulae to 3,604 carbon-determined metabolite ions, 69 of which were found to represent structures currently not listed in metabolome databases.

Converging phenomics and genomics to study natural variation in plant photosynthetic efficiency.

In recent years developments in plant phenomic approaches and facilities have gradually caught up with genomic approaches. An opportunity lies ahead to dissect complex, quantitative traits when both genotype and phenotype can be assessed at a high level of detail. This is especially true for the study of natural variation in photosynthetic efficiency, for which forward genetics studies have yielded only a little progress in our understanding of the genetic layout of the trait. High-throughput phenotyping, primarily from chlorophyll fluorescence imaging, should help to dissect the genetics of photosynthesis at the different levels of both plant physiology and development. Specific emphasis should be directed towards understanding the acclimation of the photosynthetic machinery in fluctuating environments, which may be crucial for the identification of genetic variation for relevant traits in food crops. Facilities should preferably be designed to accommodate phenotyping of photosynthesis-related traits in such environments. The use of forward genetics to study the genetic architecture of photosynthesis is likely to lead to the discovery of novel traits and/or genes that may be targeted in breeding or bio-engineering approaches to improve crop photosynthetic efficiency. In the near future, big data approaches will play a pivotal role in data processing and streamlining the phenotype-to-gene identification pipeline.

Epigenetic mapping of the Arabidopsis metabolome reveals mediators of the epigenotype-phenotype map.

Identifying the sources of natural variation underlying metabolic differences between plants will enable a better understanding of plant metabolism and provide insights into the regulatory networks that govern plant growth and morphology. So far, however, the contribution of epigenetic variation to metabolic diversity has been largely ignored. In the present study, we utilized a panel of Arabidopsis thaliana epigenetic recombinant inbred lines (epiRILs) to assess the impact of epigenetic variation on the metabolic composition. Thirty epigenetic QTL (QTLepi) were detected, which partly overlap with QTLepi linked to growth and morphology. In an effort to identify causal candidate genes in the QTLepi regions and their putative trans-targets, we performed in silico small RNA and qPCR analyses. Differentially expressed genes were further studied by phenotypic and metabolic analyses of knockout mutants. Three genes were detected that recapitulated the detected QTLepi effects, providing evidence for epigenetic regulation in cis and in trans These results indicate that epigenetic mechanisms impact metabolic diversity, possibly via small RNAs, and thus aid in further disentangling the complex epigenotype-phenotype map.

Assessment of heterosis in two Arabidopsis thaliana common-reference mapping populations.

Hybrid vigour, or heterosis, has been of tremendous importance in agriculture for the improvement of both crops and livestock. Notwithstanding large efforts to study the phenomenon of heterosis in the last decades, the identification of common molecular mechanisms underlying hybrid vigour remain rare. Here, we conducted a systematic survey of the degree of heterosis in Arabidopsis thaliana hybrids. For this purpose, two overlapping Arabidopsis hybrid populations were generated by crossing a large collection of naturally occurring accessions to two common reference lines. In these Arabidopsis hybrid populations the range of heterosis for several developmental and yield related traits was examined, and the relationship between them was studied. The traits under study were projected leaf area at 17 days after sowing, flowering time, height of the main inflorescence, number of side branches from the main stem or from the rosette base, total seed yield, seed weight, seed size and the estimated number of seeds per plant. Predominantly positive heterosis was observed for leaf area and height of the main inflorescence, whereas mainly negative heterosis was observed for rosette branching. For the other traits both positive and negative heterosis was observed in roughly equal amounts. For flowering time and seed size only low levels of heterosis were detected. In general the observed heterosis levels were highly trait specific. Furthermore, no correlation was observed between heterosis levels and the genetic distance between the parental lines. Since all selected lines were a part of the Arabidopsis genome wide association (GWA) mapping panel, a genetic mapping approach was applied to identify possible regions harbouring genetic factors causal for heterosis, with separate calculations for additive and dominance effects. Our study showed that the genetic mechanisms underlying heterosis were highly trait specific in our hybrid populations and greatly depended on the genetic background, confirming the elusive character of heterosis.

QTLepi Mapping in Arabidopsis thaliana.

While DNA sequence variation is known to be a major driver of phenotypic divergence, epigenetic variation has long been disregarded. One reason for that was the lack of suitable tools. The creation of epigenetically divergent but otherwise largely isogenic Arabidopsis populations has now alleviated some of these constraints. Epigenetic recombinant inbred line (epiRIL) populations allow for examining the effects of epigenetic variation on phenotypes. In addition, epiRILs enabled the development of epigenetic quantitative trait locus (QTLepi) mapping, an approach to identify causal epigenetic factors. Here, we describe the successive steps of QTLepi mapping in a broad sense, from the creation of epigenetically divergent populations to the identification of causal genes underlying particular phenotypes in Arabidopsis.

Parental DNA Methylation States Are Associated with Heterosis in Epigenetic Hybrids.

Despite the importance and wide exploitation of heterosis in commercial crop breeding, the molecular mechanisms behind this phenomenon are not completely understood. Recent studies have implicated changes in DNA methylation and small RNAs in hybrid performance; however, it remains unclear whether epigenetic changes are a cause or a consequence of heterosis. Here, we analyze a large panel of over 500 Arabidopsis (Arabidopsis thaliana) epigenetic hybrid plants (epiHybrids), which we derived from near-isogenic but epigenetically divergent parents. This proof-of-principle experimental system allowed us to quantify the contribution of parental methylation differences to heterosis. We measured traits such as leaf area, growth rate, flowering time, main stem branching, rosette branching, and final plant height and observed several strong positive and negative heterotic phenotypes among the epiHybrids. Using an epigenetic quantitative trait locus mapping approach, we were able to identify specific differentially methylated regions in the parental genomes that are associated with hybrid performance. Sequencing of methylomes, transcriptomes, and genomes of selected parent-epiHybrid combinations further showed that these parental differentially methylated regions most likely mediate the remodeling of methylation and transcriptional states at specific loci in the hybrids. Taken together, our data suggest that locus-specific epigenetic divergence between the parental lines can directly or indirectly trigger heterosis in Arabidopsis hybrids independent of genetic changes. These results add to a growing body of evidence that points to epigenetic factors as one of the key determinants of hybrid performance.

Genes and gene clusters related to genotype and drought-induced variation in saccharification potential, lignin content and wood anatomical traits in Populus nigra.

Wood is a renewable resource that can be employed for the production of second generation biofuels by enzymatic saccharification and subsequent fermentation. Knowledge on how the saccharification potential is affected by genotype-related variation of wood traits and drought is scarce. Here, we used three Populus nigra L. genotypes from habitats differing in water availability to (i) investigate the relationships between wood anatomy, lignin content and saccharification and (ii) identify genes and co-expressed gene clusters related to genotype and drought-induced variation in wood traits and saccharification potential. The three poplar genotypes differed in wood anatomy, lignin content and saccharification potential. Drought resulted in reduced cambial activity, decreased vessel and fiber lumina, and increased the saccharification potential. The saccharification potential was unrelated to lignin content as well as to most wood anatomical traits. RNA sequencing of the developing xylem revealed that 1.5% of the analyzed genes were differentially expressed in response to drought, while 67% differed among the genotypes. Weighted gene correlation network analysis identified modules of co-expressed genes correlated with saccharification potential. These modules were enriched in gene ontology terms related to cell wall polysaccharide biosynthesis and modification and vesicle transport, but not to lignin biosynthesis. Among the most strongly saccharification-correlated genes, those with regulatory functions, especially kinases, were prominent. We further identified transcription factors whose transcript abundances differed among genotypes, and which were co-regulated with genes for biosynthesis and modifications of hemicelluloses and pectin. Overall, our study suggests that the regulation of pectin and hemicellulose metabolism is a promising target for improving wood quality of second generation bioenergy crops. The causal relationship of the identified genes and pathways with saccharification potential needs to be validated in further experiments.

Genetic Components of Root Architecture Remodeling in Response to Salt Stress.

Salinity of the soil is highly detrimental to plant growth. Plants respond by a redistribution of root mass between main and lateral roots, yet the genetic machinery underlying this process is still largely unknown. Here, we describe the natural variation among 347 Arabidopsis thaliana accessions in root system architecture (RSA) and identify the traits with highest natural variation in their response to salt. Salt-induced changes in RSA were associated with 100 genetic loci using genome-wide association studies. Two candidate loci associated with lateral root development were validated and further investigated. Changes in CYP79B2 expression in salt stress positively correlated with lateral root development in accessions, and cyp79b2 cyp79b3 double mutants developed fewer and shorter lateral roots under salt stress, but not in control conditions. By contrast, high HKT1 expression in the root repressed lateral root development, which could be partially rescued by addition of potassium. The collected data and multivariate analysis of multiple RSA traits, available through the Salt_NV_Root App, capture root responses to salinity. Together, our results provide a better understanding of effective RSA remodeling responses, and the genetic components involved, for plant performance in stress conditions.