Mutants with simultaneous germline mutations were obtained in all three F5H genes and all three FAD2 genes (one to eleven mutated alleles) in order to improve the feed value of the seed meal and the fatty acid composition of the seed oil. In mutants with multiple mutated F5H alleles, sinapine in seed meal was reduced by up to 100%, accompanied by a sharp reduction in the S-monolignol content of lignin without causing lodging or stem break. A lower S-lignin monomer content in stems can contribute to improved stem degradability allowing new uses of stems. Mutants in all six FAD2 alleles showed an expected increase in MUFA from 8.7% to 74% and a reduction in PUFA from 53% to 13% in the fatty acids in seed oil. Remarkably, some full FAD2 mutants showed normal growth and seed production and not the dwarfing phenotype reported in previous studies. The relation between germline mutation allele dosage and phenotype was influenced by the still ongoing activity of the CRISPR/Cas9 system, leading to new somatic mutations in the leaves of flowering plants. The correlations between the total mutation frequency (germline plus new somatic mutations) for F5H with sinapine content, and FAD2 with fatty acid composition were higher than the correlations between germline mutation count and phenotypes. This shows the importance of quantifying both the germline mutations and somatic mutations when studying CRISPR/Cas9 effects in situations where the CRISPR/Cas9 system is not yet segregated out.

Syntenic cell wall QTLs (SQTLs) can identify genetic determinants of biomass traits in understudied species based on results from model crops. However, their effective use in plant breeding requires SQTLs to display intraspecific allelic variability and to predict causative loci in other populations/species than the ones used for SQTLs identification. In this study, genome assemblies from different accessions of Arabidopsis, rapeseed, tomato, rice, Brachypodium and maize were used to evaluate the intraspecific variability of SQTLs. In parallel, a genome-wide association study (GWAS) on cell wall quality traits was performed in miscanthus to verify the colocalization between GWAS loci and miscanthus SQTLs. Finally, an analogous approach was applied on a set of switchgrass cell wall QTLs retrieved from the literature. These analyses revealed large SQTLs intraspecific genetic variability, ranging from presence–absence gene variation to SNPs/INDELs and changes in coded proteins. Cell wall genes displaying gene dosage regulation, such as PAL and CAD, displayed presence–absence variation in Brachypodium and rapeseed, while protein INDELs were detected for the Brachypodium homologs of the rice brittle culm-like 8 locus, which may likely impact cell wall quality. Furthermore, SQTLs significantly colocalized with the miscanthus and switchgrass QTLs, with relevant cell wall genes being retained in colocalizing regions. Overall, SQTLs are useful tools to screen germplasm for relevant genes and alleles to improve biomass quality and can increase the efficiency of plant breeding in understudied biomass crops.

Abstract Plant cell walls of Poaceae and eudicots differ substantially, both in the content and composition of their components. However, the genomic and genetic basis underlying these differences is not fully resolved. In this research, we analyzed multiple genomic properties of 150 cell wall gene families across 169 angiosperm genomes. The properties analyzed include gene presence/absence, copy number, synteny, occurrence of tandem gene clusters, and phylogenetic gene diversity. Results revealed a profound genomic differentiation of cell wall genes between Poaceae and eudicots, often associated with the cell wall diversity between these plant groups. For example, overall patterns of gene copy number variation and synteny were clearly divergent between Poaceae and eudicot species. Moreover, differential Poaceae–eudicot copy number and genomic contexts were observed for all the genes within the BEL1-like HOMEODOMAIN 6 regulatory pathway, which respectively induces and represses secondary cell wall synthesis in Poaceae and eudicots. Similarly, divergent synteny, copy number, and phylogenetic gene diversification were observed for the major biosynthetic genes of xyloglucans, mannans, and xylans, potentially contributing to the differences in content and types of hemicellulosic polysaccharides differences in Poaceae and eudicot cell walls. Additionally, the Poaceae-specific tandem clusters and/or higher copy number of PHENYLALANINE AMMONIA-LYASE, CAFFEIC ACID O-METHYLTRANSFERASE, or PEROXIDASE genes may underly the higher content and larger variety of phenylpropanoid compounds observed in Poaceae cell walls. All these patterns are discussed in detail in this study, along with their evolutionary and biological relevance for cell wall (genomic) diversification between Poaceae and eudicots.

Establishing Lupinus mutabilis as a protein and oil crop requires improved varieties adapted to EU climates. The genetic regulation of strategic breeding traits, including plant architecture, growing cycle length and yield, is unknown. This study aimed to identify associations between 16 669 single nucleotide polymorphisms (SNPs) and 9 agronomic traits on a panel of 223 L. mutabilis accessions, grown in four environments, by applying a genome wide association study (GWAS). Seven environment-specific QTLs linked to vegetative yield, plant height, pods number and flowering time, were identified as major effect QTLs, being able to capture 6 to 20% of the phenotypic variation observed in these traits. Furthermore, two QTLs across environments were identified for flowering time on chromosome 8. The genes FAF, GAMYB and LNK, regulating major pathways involved in flowering and growth habit, as well as GA30X1, BIM1, Dr1, HDA15, HAT3, interacting with these pathways in response to hormonal and environmental cues, were prosed as candidate genes. These results are pivotal to accelerate the development of L. mutabilis varieties adapted to European cropping conditions by using marker-assisted selection (MAS), as well as to provide a framework for further functional studies on plant development and phenology in this species.

Translational genomics can enable a quicker improvement of orphan crops toward novel agricultural applications, including the advancement of orphan biomass species for cultivation on marginal lands. In this sense, cell wall quality is a preeminent breeding target. However, tools to efficiently project genetic data on target traits across large sets of species are currently missing. This study aimed at closing this gap by developing a strategy to project a set of cell wall QTLs across a large group of plants by using genome synteny. This strategy is suited for large-scale analyses and detected 362 syntenic cell wall QTLs (SQTLs) across 74 angiosperms, including several (orphan) biomass species. SQTLs analyses revealed that they span large portions of the initial cell wall QTLs and are extensively conserved across diverse species. Moreover, numerous QTLs cell wall genes were conserved through SQTLs, including genes displaying allelic variation associated with cell wall composition. Functional analyses showed that highly conserved genes of SQTLs include important cell wall transcription factors and genes involved in the remodeling of cell wall polymers. For some of these gene families, SQTLs indicated the presence of differentially conserved genomic contexts for different gene members, highlighting their utility as a tool to pinpoint gene targets that maximize the likelihood of functional gene conservation. Overall, the results of this study can facilitate “universal” approaches for breeding (orphan) biomass crops, while the strategy for QTLs translation can be applied to other sets of traits and species, helping to unlock the potential of orphan species.

The Cellulose synthase superfamily synthesizes cellulose and different hemicellulosic polysaccharides in plant cell walls. While much has been discovered about the evolution and function of these genes, their genomic architecture and relationship with gene (sub-)functionalization and evolution remains unclear. By using 242 genomes covering plant evolution from green algae to eudicots, we performed a large-scale analysis of synteny, phylogenetic, and functional data of the CesA superfamily. Results revealed considerable gene copy number variation across species and gene families, and also two patterns – singletons vs. tandem arrays – in chromosomic gene arrangement. Synteny analysis revealed exceptional conservation of gene architecture across species, but also lineage-specific patterns across gene (sub-)families. Synteny patterns correlated with gene sub-functionalization into primary and secondary CesAs and distinct CslD functional isoforms. Furthermore, a genomic context shift of a group of cotton secondary CesAs was associated with peculiar properties of cotton fiber synthesis. Finally, phylogenetics suggested that primary CesA sequences appeared before the secondary CesAs, while phylogenomic analyses unveiled the genomic trace of the CslD duplication that initiated the CslF family. Our results describe in detail the genomic architecture of the CesA superfamily in plants, highlighting its crucial relevance for gene diversification and sub-functionalization, and for understanding their evolution.

The Plantarray 3.0 phenotyping platform ® was used to monitor the growth and water use of the quinoa varieties Pasto and selRiobamba under salinity (0–300 mM NaCl). Salinity reduced the cumulative transpiration of both varieties by 60% at 200 mM NaCl and by 75 and 82% at 300 mM NaCl for selRiobamba and Pasto, respectively. Stomatal conductance was reduced by salinity, but at 200 mM NaCl Pasto showed a lower reduction (15%) than selRiobamba (35%), along with decreased specific leaf area. Diurnal changes in water use parameters indicate that under salt stress, daily transpiration in quinoa is less responsive to changes in light irradiance, and stomatal conductance is modulated to maximize CO 2 uptake and minimize water loss following the changes in VPD (vapor pressure deficit). These changes might contribute to the enhanced water use efficiency of both varieties under salt stress. The mechanistic crop model LINTUL was used to integrate physiological responses into the radiation use efficiency of the plants (RUE), which was more reduced in Pasto than selRiobamba under salinity. By the end of the experiment (eleven weeks after sowing, six weeks after stress), the growth of Pasto was significantly lower than selRiobamba, fresh biomass was 50 and 35% reduced at 200 mM and 70 and 50% reduced at 300 mM NaCl for Pasto and selRiobamba, respectively. We argue that contrasting water management strategies can at least partly explain the differences in salt tolerance between Pasto and selRiobamba. Pasto adopted a “conservative-growth” strategy, saving water at the expense of growth, while selRiobamba used an “acquisitive-growth” strategy, maximizing growth in spite of the stress. The implementation of high-resolution phenotyping could help to dissect these complex growth traits that might be novel breeding targets for abiotic stress tolerance.

Quinoa is a nutritious seed crop with a great potential to grow in saline soils. Here, we studied ion concentrations in quinoa tissues throughout the life cycle of the plant, and linked ion dynamics to responses in growth parameters, seed yield and efficiency of photosynthesis under salinity (0-400 mM NaCl). Ion dynamics changed from high ion exclusion (> 99 %, root contents lower than root medium and low accumulation of ions in the leaves) before flowering, to a build-up of ions during seed filling. This indicates a change in strategy in maintaining the necessary gradient of water potential from the root medium to the leaves. K+ concentrations in leaves also increased by more than 100 % in response to prolonged severe salt stress, which may point to a role of this ion in leaf osmotic adjustment. Accumulation of ions in epidermal bladder cells did not contribute substantially to Na+-exclusion as it was less than 6 % of the total Na+ taken up in leaves. Growth under salt stress was mostly impaired by anatomical adaptations (reduced SLA), while initial light use efficiency (Fv/Fm) and NAR were not affected. The variety Pasto showed a "survival strategy" to high salinity with higher ion exclusion and a higher reduction in transpiration than the other varieties, at the expense of lower biomass and seed yield.

Hemp (Cannabis sativa L.) is a bast-fiber crop with a great potential in the emerging bio-based economy. Yet, hemp breeding for fiber quality is restricted and that is mainly due to the limited knowledge of the genetic architecture of its fiber quality. A panel of 123 hemp accessions, with large phenotypic variability, was used to study the genetic basis of seven cell wall and bast fiber traits relevant to fiber quality. These traits showed large genetic variance components and high values of broad sense heritability in this hemp panel, as concluded from the phenotypic evaluation across three test locations with contrasting environments. The hemp panel was genotyped using restriction site associated DNA sequencing (RAD-seq). Subsequently, a large set (> 600,000) of selected genome-wide single nucleotide polymorphism (SNP) markers was used for a genome-wide association study (GWAS) approach to get insights into quantitative trait loci (QTLs) controlling fiber quality traits. In absence of a complete hemp genome sequence, identification of QTLs was based on the following characteristics: (i) association level to traits, (ii) fraction of explained trait variance, (iii) collinearity between QTLs, and (iv) detection across different environments. Using this approach, 16 QTLs were identified across locations for different fiber quality traits, including contents of glucose, glucuronic acid, mannose, xylose, lignin, and bast fiber content. Among them, six were found across the three environments. The genetic markers composing the QTLs that are common across locations are valuable tools to develop novel genotypes of hemp with improved fiber quality. Underneath the QTLs, 12 candidate genes were identified which are likely to be involved in the biosynthesis and modification of monosaccharides, polysaccharides, and lignin. These candidate genes were suggested to play an important role in determining fiber quality in hemp. This study provides new insights into the genetic architecture of fiber traits, identifies QTLs and candidate genes that form the basis for molecular breeding for high fiber quality hemp cultivars.

This study deals with approaches for a social-ecological friendly European bioeconomy based on biomass from industrial crops cultivated on marginal agricultural land. The selected crops to be investigated are: Biomass sorghum, camelina, cardoon, castor, crambe, Ethiopian mustard, giant reed, hemp, lupin, miscanthus, pennycress, poplar, reed canary grass, safflower, Siberian elm, switchgrass, tall wheatgrass, wild sugarcane, and willow. The research question focused on the overall crop growth suitability under low-input management. The study assessed: (i) How the growth suitability of industrial crops can be defined under the given natural constraints of European marginal agricultural lands; and (ii) which agricultural practices are required for marginal agricultural land low-input systems (MALLIS). For the growth-suitability analysis, available thresholds and growth requirements of the selected industrial crops were defined. The marginal agricultural land was categorized according to the agro-ecological zone (AEZ) concept in combination with the marginality constraints, so-called ‘marginal agro-ecological zones’ (M-AEZ). It was found that both large marginal agricultural areas and numerous agricultural practices are available for industrial crop cultivation on European marginal agricultural lands. These results help to further describe the suitability of industrial crops for the development of social-ecologically friendly MALLIS in Europe.

Increased soil salinity is a significant agricultural problem that decreases yields for common agricultural crops. Its dynamics require cost and labour effective measurement techniques and widely acknowledged methods are not present yet. We investigated the potential of Unmanned Aerial Vehicle (UAV) remote sensing to measure salt stress in quinoa plants. Three different UAV sensors were used: a WIRIS thermal camera, a Rikola hyperspectral camera and a Riegl VUX-SYS Light Detection and Ranging (LiDAR) scanner. Several vegetation indices, canopy temperature and LiDAR measured plant height were derived from the remote sensing data and their relation with ground measured parameters like salt treatment, stomatal conductance and actual plant height is analysed. The results show that widely used multispectral vegetation indices are not efficient in discriminating between salt affected and control quinoa plants. The hyperspectral Physiological Reflectance Index (PRI) performed best and showed a clear distinction between salt affected and treated plants. This distinction is also visible for LiDAR measured plant height, where salt treated plants were on average 10 cm shorter than control plants. Canopy temperature was significantly affected, though detection of this required an additional step in analysis - Normalised Difference Vegetation Index (NDVI) clustering. This step assured temperature comparison for equally vegetated pixels. Data combination of all three sensors in a Multiple Linear Regression model increased the prediction power and for the whole dataset R-2 reached 0.46, with some subgroups reaching an R-2 of 0.64. We conclude that UAV borne remote sensing is useful for measuring salt stress in plants and a combination of multiple measurement techniques is advised to increase the accuracy.

Chenopodium quinoa (quinoa) is a highly nutritious grain identified as an important crop to improve world food security. Unfortunately, few resources are available to facilitate its genetic improvement. Here we report the assembly of a high-quality, chromosome-scale reference genome sequence for quinoa, which was produced using single-molecule real-time sequencing in combination with optical, chromosome-contact and genetic maps. We also report the sequencing of two diploids from the ancestral gene pools of quinoa, which enables the identification of sub-genomes in quinoa, and reduced-coverage genome sequences for 22 other samples of the allotetraploid goosefoot complex. The genome sequence facilitated the identification of the transcription factor likely to control the production of anti-nutritional triterpenoid saponins found in quinoa seeds, including a mutation that appears to cause alternative splicing and a premature stop codon in sweet quinoa strains. These genomic resources are an important first step towards the genetic improvement of quinoa.

Camelina (Camelina sativa L. Crantz) is considered a relatively new oilseed Brassicacea in both Europe and North America, even though its history as a crop dates back to the Bronze Age. Camelina has recently received renewed interest from both the scientific community and bio-based industries around the world. The main attractive features of this species are: drought and frost tolerance, disease and pest resistance, a unique seed oil composition with high levels of n-3 fatty acids, a considerably high seed oil content, and satisfactory seed yields, in particular under low-input management and in limiting environments. Aiming at evaluating the feasible introduction of recently released camelina breeding lines under different environmental conditions and their productive potential a multi-location trial was set up. The agronomic performance of nine improved genotypes of camelina was evaluated in a wide range of environments in Europe (Greece, Italy, Poland) and in five locations across Canada, in two consecutive growing seasons (2015 and 2016). Sowing time was optimized for each location according to the different climatic conditions. Camelina proved to be a highly adaptable species, reaching seed yields of about 1MgDMha−1 under the most limiting conditions (i.e., low precipitation, poor soil quality, extremely high temperature at flowering). Growing environments characterized by mild temperatures and adequate rainfall (>170mm, during the growing season) resulted in higher average seed yields. The length of the growing cycle varied greatly between different locations (80–110d), but the cumulative thermal time was quite stable (∼1200 GDD, growing degree days). The advanced breeding line 787–08, which possesses up to 30% larger seed compared to the mean seed size of all other test entries, proved to be the most promising genotype across all locations in Europe and Canada, combining high seed yields (1.1–2.7MgDMha−1) with improved yield stability. To the best of our knowledge, for the first time, camelina lines with improved oil composition (i.e., increased oleic and α-linolenic and lower linoleic acid contents) for feed, food and industrial applications were identified (789–02 and 887).

Crambe abyssinica (crambe) is a new industrial oil crop that can grow on saline soil and tolerates salty water irrigation. Genetically engineered crambe in which the seed‐oil composition is manipulated for more erucic acid and less polyunsaturated fatty acid (PUFA) would be highly beneficial to industry. In this research, lysophosphatidic acid acyltransferase 2 RNA interference (CaLPAT2‐RNAi) was introduced into the crambe genome to manipulate its oil composition. The result showed in comparison with wild type, CaLPAT2‐RNAi could significantly reduce linoleic and linolenic acid content, simultaneously increasing erucic acid content. Systematic metabolism engineering was then carried out to further study CaLPAT2‐RNAi, combined with the overexpression of Brassica napus fatty acid elongase (BnFAE), Limnanthes douglasii LPAT (LdLPAT), and RNAi of endogenous fatty acid desaturase 2 (CaFAD2‐RNAi). Oil composition analysis on the tranformants' seeds showed that (a) with CaFAD2‐RNAi, PUFA content could be dramatically decreased, in comparison with BnFAE + LdLPAT + CaFAD2‐RNAi, and BnFAE + LdLPAT + CaFAD2‐RNAi + CaLPAT2‐RNAi seeds showed lower linolenic acid content; (b) BnFAE + LdLPAT + CaFAD2‐RNAi + CaLPAT2‐RNAi could increase the erucic acid content in crambe seed oil from less than 66.6% to 71.6%, whereas the highest erucic acid content of BnFAE + LdLPAT + CaFAD2‐RNAi was 79.2%; (c) although the four‐gene combination could not increase the erucic acid content of seed oil to a higher level than the others, it led to increased carbon resource deposited into C22:1 and C18:1 moieties and lower PUFA. Summarily, the present research indicates that suppression of LPAT2 is a new, promising strategy for seed‐oil biosynthesis pathway engineering, which would increase the value of crambe oil.

Nature 542, 307–312 (2017); doi:10.1038/nature21370 The Acknowledgements section of this Article should have included the following sentence: “Metabolite imaging was conducted at Metabolomics Australia (School of BioSciences, The University of Melbourne, Australia), a NCRIS initiative under Bioplatforms Australia Pty Ltd.

The life cycle based greenhouse gas (GHG) balances of Fatty Acid Methyl Esters (FAME also called “Biodiesel”) from various resources have been set in the Renewable Energy Directive (RED). Due to technology and scientific progress there are various options to improve the GHG balances of FAME. In this Supporting Action 10 most interesting options were assessed: 1) “Biomethanol”: Substitution of fossil methanol with biomethanol; 2) “Bioethanol”: Substitution of fossil methanol with bioethanol; 3) “CHP residues”: Use of residues and co-products in an CHP plant; 4) “New plant species”: Examination of new plants for vegetable oils, that could increase the biomass weight without any detrimental effect on the oil seed; 5) “Bioplastics and biochemicals”: Production of bioplastics and biochemicals from process residues; 6) “Advanced agriculture”: Advanced agricultural practices in terms of N2O emissions and soil carbon accumulation; 7) “Organic residues”: Use of organic versus mineral fertilizer for feedstock cultivation; 8) “FAME as fuel”: Use of FAME in machinery for cultivation, transportation and distribution; 9) “Retrofitting multi feedstock”: Retrofitting of single feedstock plants for blending fatty residues; and 10) “Green electricity”: Use of renewable electricity produced in a PV plant on site. The assessment approach started with the GHG standard values of the RED and the corresponding background data documented in BioGrace. For the most relevant FAME production possibilities in Europe, characterized by the feedstock (rapeseed, sunflower, palm oil, soybean, used cooking oil, animal fat) and FAME production capacity (50 - 200 kt/a), the technical and economic data of “Best Available Technology in 2015” (BAT) were used as starting point to assess the improvement options. Based on the calculation of GHG emissions (g CO2-eq/MJ) and production cost (€/tFAME) an overall assessment (incl SWOT-Analyses and Stakeholder involvement) of the options was made and summarized in “Fact Sheets”. A significant GHG reduction compared to the RED values in processing is possible, if best available technology (BAT) is applied. The GHG emissions of cultivation compared to RED are higher due to improved data on the correlation between fertilizer input and yields. The assessed GHG improvements options show that the potential to reduce emissions is relatively large in agriculture cultivation, but a relatively low in processing. The production cost analysis shows that revenues from co-produced animal feed and oil yield per hectare have a strong influence on total production costs, e.g. mainly animal feed from soybeans. The total FAME production cost of BAT are 280 – 1,000 €/tFAME, including revenues from co-products. Cost ranges arise due to different feedstock and capacities. The greenhouse gas analysis of the improvement options results in a GHG reduction potential of 0 - 37 g CO2-eq/MJ compared to BAT. The greenhouse gas mitigation costs of improvement options range between -260 and +1,000 €/t CO2-eq. Options with negative greenhouse gas mitigation costs generate economic benefits compared to the base case. Summing up the assessment one can conclude that the future FAME production has several options to further improve its GHG balance thus contributing substantially to a more sustainable transportation sector.

Feedstocks for industrial applications ranging from polymers to lubricants are largely derived from petroleum, a non-renewable resource. Vegetable oils with fatty acid structures and storage forms tailored for specific industrial uses offer renewable and potentially sustainable sources of petrochemical-type functionalities. A wide array of industrial vegetable oils can be generated through biotechnology, but will likely require non-commodity oilseed platforms dedicated to specialty oil production for commercial acceptance. Here we show the feasibility of three Brassicaceae oilseeds crambe, camelina, and carinata, none of which are widely cultivated for food use, as hosts for complex metabolic engineering of wax esters for lubricant applications. Lines producing wax esters >20% of total seed oil were generated for each crop and further improved for high temperature oxidative stability by down-regulation of fatty acid polyunsaturation. Field cultivation of optimized wax ester-producing crambe demonstrated commercial utility of these engineered crops and a path for sustainable production of other industrial oils in dedicated specialty oilseeds.

BACKGROUND:Economical cultivation of the oilseed crop Jatropha curcas is currently hampered in part due to the non-availability of purpose-bred cultivars. Although genetic maps and genome sequence data exist for this crop, marker-assisted breeding has not yet been implemented due to a lack of available marker-trait association studies. To identify the location of beneficial alleles for use in plant breeding, we performed quantitative trait loci (QTL) analysis for a number of agronomic traits in two biparental mapping populations.RESULTS:The mapping populations segregated for a range of traits contributing to oil yield, including plant height, stem diameter, number of branches, total seeds per plant, 100-seed weight, seed oil content and fatty acid composition. QTL were detected for each of these traits and often over multiple years, with some variation in the phenotypic variance explained between different years. In one of the mapping populations where we recorded vegetative traits, we also observed co-localization of QTL for stem diameter and plant height, which were both overdominant, suggesting a possible locus conferring a pleotropic heterosis effect. By using a candidate gene approach and integrating physical mapping data from a recent high-quality release of the Jatropha genome, we were also able to position a large number of genes involved in the biosynthesis of storage lipids onto the genetic map. By comparing the position of these genes with QTL, we were able to detect a number of genes potentially underlying seed traits, including phosphatidate phosphatase genes.CONCLUSIONS:The QTL we have identified will serve as a useful starting point in the creation of new varieties of J. curcas with improved agronomic performance for seed and oil productivity. Our ability to physically map a significant proportion of the Jatropha genome sequence onto our genetic map could also prove useful in identifying the genes underlying particular traits, allowing more controlled and precise introgression of desirable alleles and permitting the pyramiding or stacking of multiple QTL.