Integrated multi-omics analyses and genome-wide association studies reveal prime candidate genes of metabolic and vegetative growth variation in canola.

Genome-wide association studies (GWAS) identified thousands of genetic loci associated with complex plant traits, including many traits of agronomical importance. However, functional interpretation of GWAS results remains challenging because of large candidate regions due to linkage disequilibrium. High-throughput omics technologies, such as genomics, transcriptomics, proteomics and metabolomics open new avenues for integrative systems biological analyses and help to nominate systems information supported (prime) candidate genes. In the present study, we capitalise on a diverse canola population with 477 spring-type lines which was previously analysed by high-throughput phenotyping of growth-related traits and by RNA sequencing and metabolite profiling for multi-omics-based hybrid performance prediction. We deepened the phenotypic data analysis, now providing 123 time-resolved image-based traits, to gain insight into the complex relations during early vegetative growth and reanalysed the transcriptome data based on the latest Darmor-bzh v10 genome assembly. Genome-wide association testing revealed 61 298 robust quantitative trait loci (QTL) including 187 metabolite QTL, 56814 expression QTL and 4297 phenotypic QTL, many clustered in pronounced hotspots. Combining information about QTL colocalisation across omics layers and correlations between omics features allowed us to discover prime candidate genes for metabolic and vegetative growth variation. Prioritised candidate genes for early biomass accumulation include A06p05760.1_BnaDAR (PIAL1), A10p16280.1_BnaDAR, C07p48260.1_BnaDAR (PRL1) and C07p48510.1_BnaDAR (CLPR4). Moreover, we observed unequal effects of the Brassica A and C subgenomes on early biomass production.

Strong prevalence of light regime-specific QTL in Arabidopsis detected using automated high-throughput phenotyping in fluctuating or constant light.

Plants have evolved and adapted under dynamic environmental conditions, particularly to fluctuating light, but plant research has often focused on constant growth conditions. To quantitatively asses the adaptation to fluctuating light, a panel of 384 natural Arabidopsis thaliana accessions was analyzed in two parallel independent experiments under fluctuating and constant light conditions in an automated high-throughput phenotyping system upgraded with supplemental LEDs. While the integrated daily photosynthetically active radiation was the same under both light regimes, plants in fluctuating light conditions accumulated significantly less biomass and had lower leaf area during their measured vegetative growth than plants in constant light. A total of 282 image-derived architectural and/or color-related traits at six common time points, and 77 photosynthesis-related traits from one common time point were used to assess their associations with genome-wide natural variation for both light regimes. Out of the 3000 significant marker-trait associations (MTAs) detected, only 183 (6.1%) were common for fluctuating and constant light conditions. The prevalence of light regime-specific QTL indicates a complex adaptation. Genes in linkage disequilibrium with fluctuating light-specific MTAs with an adjusted repeatability value >0.5 were filtered for gene ontology terms containing "photo" or "light", yielding 15 selected candidates. The candidate genes are involved in photoprotection, PSII maintenance and repair, maintenance of linear electron flow, photorespiration, phytochrome signaling, and cell wall expansion, providing a promising starting point for further investigations into the response of Arabidopsis thaliana to fluctuating light conditions.

Dynamic growth QTL action in diverse light environments: characterization of light regime-specific and stable QTL in Arabidopsis.

Plant growth is a complex process affected by a multitude of genetic and environmental factors and their interactions. To identify genetic factors influencing plant performance under different environmental conditions, vegetative growth was assessed in Arabidopsis thaliana cultivated under constant or fluctuating light intensities, using high-throughput phenotyping and genome-wide association studies. Daily automated non-invasive phenotyping of a collection of 382 Arabidopsis accessions provided growth data during developmental progression under different light regimes at high temporal resolution. Quantitative trait loci (QTL) for projected leaf area, relative growth rate, and PSII operating efficiency detected under the two light regimes were predominantly condition-specific and displayed distinct temporal activity patterns, with active phases ranging from 2 d to 9 d. Eighteen protein-coding genes and one miRNA gene were identified as potential candidate genes at 10 QTL regions consistently found under both light regimes. Expression patterns of three candidate genes affecting projected leaf area were analysed in time-series experiments in accessions with contrasting vegetative leaf growth. These observations highlight the importance of considering both environmental and temporal patterns of QTL/allele actions and emphasize the need for detailed time-resolved analyses under diverse well-defined environmental conditions to effectively unravel the complex and stage-specific contributions of genes affecting plant growth processes.

Temporal dynamics of QTL effects on vegetative growth in Arabidopsis thaliana.

We assessed early vegetative growth in a population of 382 accessions of Arabidopsis thaliana using automated non-invasive high-throughput phenotyping. All accessions were imaged daily from 7 d to 18 d after sowing in three independent experiments and genotyped using the Affymetrix 250k SNP array. Projected leaf area (PLA) was derived from image analysis and used to calculate relative growth rates (RGRs). In addition, initial seed size was determined. The generated datasets were used jointly for a genome-wide association study that identified 238 marker-trait associations (MTAs) individually explaining up to 8% of the total phenotypic variation. Co-localization of MTAs occurred at 33 genomic positions. At 21 of these positions, sequential co-localization of MTAs for 2-9 consecutive days was observed. The detected MTAs for PLA and RGR could be grouped according to their temporal expression patterns, emphasizing that temporal variation of MTA action can be observed even during the vegetative growth phase, a period of continuous formation and enlargement of seemingly similar rosette leaves. This indicates that causal genes may be differentially expressed in successive periods. Analyses of the temporal dynamics of biological processes are needed to gain important insight into the molecular mechanisms of growth-controlling processes in plants.

Multi-omics-based prediction of hybrid performance in canola.

KEY MESSAGE: Complementing or replacing genetic markers with transcriptomic data and use of reproducing kernel Hilbert space regression based on Gaussian kernels increases hybrid prediction accuracies for complex agronomic traits in canola. In plant breeding, hybrids gained particular importance due to heterosis, the superior performance of offspring compared to their inbred parents. Since the development of new top performing hybrids requires labour-intensive and costly breeding programmes, including testing of large numbers of experimental hybrids, the prediction of hybrid performance is of utmost interest to plant breeders. In this study, we tested the effectiveness of hybrid prediction models in spring-type oilseed rape (Brassica napus L./canola) employing different omics profiles, individually and in combination. To this end, a population of 950 F1 hybrids was evaluated for seed yield and six other agronomically relevant traits in commercial field trials at several locations throughout Europe. A subset of these hybrids was also evaluated in a climatized glasshouse regarding early biomass production. For each of the 477 parental rapeseed lines, 13,201 single nucleotide polymorphisms (SNPs), 154 primary metabolites, and 19,479 transcripts were determined and used as predictive variables. Both, SNP markers and transcripts, effectively predict hybrid performance using (genomic) best linear unbiased prediction models (gBLUP). Compared to models using pure genetic markers, models incorporating transcriptome data resulted in significantly higher prediction accuracies for five out of seven agronomic traits, indicating that transcripts carry important information beyond genomic data. Notably, reproducing kernel Hilbert space regression based on Gaussian kernels significantly exceeded the predictive abilities of gBLUP models for six of the seven agronomic traits, demonstrating its potential for implementation in future canola breeding programmes.

Strong temporal dynamics of QTL action on plant growth progression revealed through high-throughput phenotyping in canola.

A major challenge of plant biology is to unravel the genetic basis of complex traits. We took advantage of recent technical advances in high-throughput phenotyping in conjunction with genome-wide association studies to elucidate genotype-phenotype relationships at high temporal resolution. A diverse Brassica napus population from a commercial breeding programme was analysed by automated non-invasive phenotyping. Time-resolved data for early growth-related traits, including estimated biovolume, projected leaf area, early plant height and colour uniformity, were established and complemented by fresh and dry weight biomass. Genome-wide SNP array data provided the framework for genome-wide association analyses. Using time point data and relative growth rates, multiple robust main effect marker-trait associations for biomass and related traits were detected. Candidate genes involved in meristem development, cell wall modification and transcriptional regulation were detected. Our results demonstrate that early plant growth is a highly complex trait governed by several medium and many small effect loci, most of which act only during short phases. These observations highlight the importance of taking the temporal patterns of QTL/allele actions into account and emphasize the need for detailed time-resolved analyses to effectively unravel the complex and stage-specific contributions of genes affecting growth processes that operate at different developmental phases.

Genetic dissection of metabolite variation in Arabidopsis seeds: evidence for mQTL hotspots and a master regulatory locus of seed metabolism.

To gain insight into genetic factors controlling seed metabolic composition and its relationship to major seed properties, an Arabidopsis recombinant inbred line (RIL) population, derived from accessions Col-0 and C24, was studied using an MS-based metabolic profiling approach. Relative intensities of 311 polar primary metabolites were used to identify associated genomic loci and to elucidate their interactions by quantitative trait locus (QTL) mapping. A total of 786 metabolic QTLs (mQTLs) were unequally distributed across the genome, forming several hotspots. For the branched-chain amino acid leucine, mQTLs and candidate genes were elucidated in detail. Correlation studies displayed links between metabolite levels, seed protein content, and seed weight. Principal component analysis revealed a clustering of samples, with PC1 mapping to a region on the short arm of chromosome IV. The overlap of this region with mQTL hotspots indicates the presence of a potential master regulatory locus of seed metabolism. As a result of database queries, a series of candidate regulatory genes, including bZIP10, were identified within this region. Depending on the search conditions, metabolic pathway-derived candidate genes for 40-61% of tested mQTLs could be determined, providing an extensive basis for further identification and characterization of hitherto unknown genes causal for natural variation of Arabidopsis seed metabolism.

Native-sized spider silk proteins synthesized in planta via intein-based multimerization.

The synthesis of native-sized proteins is a pre-requisite for exploiting the potential of spider silk as a bio-based material. The unique properties of spider silk, such as extraordinary tensile strength and elasticity, result from the highly repetitive nature of spider silk protein motifs. The present report describes the combination of spider silk flagelliform protein (FLAG) production in the endoplasmic reticulum of tobacco plant leaf cells with an intein-based posttranslational protein fusion technology. The repeated ligation of FLAG monomers resulted in the formation of large multimers. This method avoids the need for highly repetitive transgenes, which may result in a higher genetic and transcriptional stability. Here we show, for the first time, the production of synthetic, high molecular weight spider silk proteins larger than 250 kDa based on the assembly of protein monomers via intein-mediated trans-splicing in planta. The resulting multimeric structures form microfibers, thereby demonstrating their great potential as a biomaterial.

Natural plant growth and development achieved in the IPK PhenoSphere by dynamic environment simulation.

In plant science, the suboptimal match of growing conditions hampers the transfer of knowledge from controlled environments in glasshouses or climate chambers to field environments. Here we present the PhenoSphere, a plant cultivation infrastructure designed to simulate field-like environments in a reproducible manner. To benchmark the PhenoSphere, the effects on plant growth of weather conditions of a single maize growing season and of an averaged season over three years are compared to those of a standard glasshouse and of four years of field trials. The single season simulation proves superior to the glasshouse and the averaged season in the PhenoSphere: The simulated weather regime of the single season triggers plant growth and development progression very similar to that observed in the field. Hence, the PhenoSphere enables detailed analyses of performance-related trait expression and causal biological mechanisms in plant populations exposed to weather conditions of current and anticipated future climate scenarios.

Integrated phenotyping of root and shoot growth dynamics in maize reveals specific interaction patterns in inbreds and hybrids and in response to drought.

In recent years, various automated methods for plant phenotyping addressing roots or shoots have been developed and corresponding platforms have been established to meet the diverse requirements of plant research and breeding. However, most platforms are only either able to phenotype shoots or roots of plants but not both simultaneously. This substantially limits the opportunities offered by a joint assessment of the growth and development dynamics of both organ systems, which are highly interdependent. In order to overcome these limitations, a root phenotyping installation was integrated into an existing automated non-invasive high-throughput shoot phenotyping platform. Thus, the amended platform is now capable of conducting high-throughput phenotyping at the whole-plant level, and it was used to assess the vegetative root and shoot growth dynamics of five maize inbred lines and four hybrids thereof, as well as the responses of five inbred lines to progressive drought stress. The results showed that hybrid vigour (heterosis) occurred simultaneously in roots and shoots and was detectable as early as 4 days after transplanting (4 DAT; i.e., 8 days after seed imbibition) for estimated plant height (EPH), total root length (TRL), and total root volume (TRV). On the other hand, growth dynamics responses to progressive drought were different in roots and shoots. While TRV was significantly reduced 10 days after the onset of the water deficit treatment, the estimated shoot biovolume was significantly reduced about 6 days later, and EPH showed a significant decrease even 2 days later (8 days later than TRV) compared with the control treatment. In contrast to TRV, TRL initially increased in the water deficit period and decreased much later (not earlier than 16 days after the start of the water deficit treatment) compared with the well-watered plants. This may indicate an initial response of the plants to water deficit by forming longer but thinner roots before growth was inhibited by the overall water deficit. The magnitude and the dynamics of the responses were genotype-dependent, as well as under the influence of the water consumption, which was related to plant size.

Genome-wide association study in two-row spring barley landraces identifies QTL associated with plantlets root system architecture traits in well-watered and osmotic stress conditions.

Water availability is undoubtedly one of the most important environmental factors affecting crop production. Drought causes a gradual deprivation of water in the soil from top to deep layers and can occur at diverse stages of plant development. Roots are the first organs that perceive water deficit in soil and their adaptive development contributes to drought adaptation. Domestication has contributed to a bottleneck in genetic diversity. Wild species or landraces represent a pool of genetic diversity that has not been exploited yet in breeding program. In this study, we used a collection of 230 two-row spring barley landraces to detect phenotypic variation in root system plasticity in response to drought and to identify new quantitative trait loci (QTL) involved in root system architecture under diverse growth conditions. For this purpose, young seedlings grown for 21 days in pouches under control and osmotic-stress conditions were phenotyped and genotyped using the barley 50k iSelect SNP array, and genome-wide association studies (GWAS) were conducted using three different GWAS methods (MLM GAPIT, FarmCPU, and BLINK) to detect genotype/phenotype associations. In total, 276 significant marker-trait associations (MTAs; p-value (FDR)< 0.05) were identified for root (14 and 12 traits under osmotic-stress and control conditions, respectively) and for three shoot traits under both conditions. In total, 52 QTL (multi-trait or identified by at least two different GWAS approaches) were investigated to identify genes representing promising candidates with a role in root development and adaptation to drought stress.