Translational genomics using Arabidopsis as a model enables the characterization of pennycress genes through forward and reverse genetics

Transcriptomic and lipidomic analysis of the differential pathway contribution to the incorporation of erucic acid to triacylglycerol during Pennycress seed maturation

Thlaspi arvense (Pennycress) is an emerging feedstock for biofuel production because of its high seed oil content enriched in erucic acid. A transcriptomic and a lipidomic study were performed to analyze the dynamics of gene expression, glycerolipid content and acyl-group distribution during seed maturation. Genes involved in fatty acid biosynthesis were expressed at the early stages of seed maturation. Genes encoding enzymes of the Kennedy pathway like diacylglycerol acyltransferase1 (TaDGAT1), lysophosphatidic acid acyltransferase (TaLPAT) or glycerol 3-phosphate acyltransferase (TaGPAT) increased their expression with maturation, coinciding with the increase in triacylglycerol species containing 22:1. Positional analysis showed that the most abundant triacylglycerol species contained 18:2 at sn-2 position in all maturation stages, suggesting no specificity of the lysophosphatidic acid acyltransferase for very long chain fatty acids. Diacylglycerol acyltransferase2 (TaDGAT2) mRNA was more abundant at the initial maturation stages, coincident with the rapid incorporation of 22:1 to triacylglycerol, suggesting a coordination between Diacylglycerol acyltransferase enzymes for triacylglycerol biosynthesis. Genes encoding the phospholipid-diacylglycerol acyltransferase (TaPDAT1), lysophosphatidylcholine acyltransferase (TaLPCAT) or phosphatidylcholine diacylglycerolcholine phosphotransferase (TaPDCT), involved in acyl-editing or phosphatidyl-choline (PC)-derived diacylglycerol (DAG) biosynthesis showed also higher expression at the early maturation stages, coinciding with a higher proportion of triacylglycerol containing C18 fatty acids. These results suggested a higher contribution of these two pathways at the early stages of seed maturation. Lipidomic analysis of the content and acyl-group distribution of diacylglycerol and phosphatidyl-choline pools was compatible with the acyl content in triacylglycerol at the different maturation stages. Our data point to a model in which a strong temporal coordination between pathways and isoforms in each pathway, both at the expression and acyl-group incorporation, contribute to high erucic triacylglycerol accumulation in Pennycress.

Arabidopsis as a model for translational research

Arabidopsis thaliana is currently the most-studied plant species on earth, with an unprecedented number of genetic, genomic, and molecular resources having been generated in this plant model. In the era of translating foundational discoveries to crops and beyond, we aimed to highlight the utility and challenges of using Arabidopsis as a reference for applied plant biology research, agricultural innovation, biotechnology, and medicine. We hope that this review will inspire the next generation of plant biologists to continue leveraging Arabidopsis as a robust and convenient experimental system to address fundamental and applied questions in biology. We aim to encourage lab and field scientists alike to take advantage of the vast Arabidopsis datasets, annotations, germplasm, constructs, methods, molecular and computational tools in our pursuit to advance understanding of plant biology and help feed the world's growing population. We envision that the power of Arabidopsis-inspired biotechnologies and foundational discoveries will continue to fuel the development of resilient, high-yielding, nutritious plants for the betterment of plant and animal health and greater environmental sustainability.

Performance and health of broiler chickens fed low erucic acid, lower fiber pennycress (CoverCressTM) grain

CoverCress (low erucic acid, lower fiber pennycress) is being developed as a cover crop to be planted in the fall after corn and harvested in the spring prior to planting soybeans. Two experiments were conducted to evaluate 2 lines of the whole grain (CCWG-1: natural mutation and mutation breeding; CCWG-2: gene edited) and the whole grain pretreated with the potential palatability agent copper sulfate (CCWG-1-CuSO₄; CCWG-2-CuSO₄) as an ingredient for broilers. In Experiment 1, CCWG-1-CuSO₄ was included in the diet at 0, 4, and 6% for 41 d. Feed intake, body weight gain, feed conversion, processing characteristics, organ weights, serum thyroid, macropathology and histology data were collected. In Experiment 2, broilers were fed diets containing Control, 2% CCWG-1, 4% CCWG-1, 4% CCWG-2, and 4.35% CCWG-1-CuSO₄ for 42 d. Feed intake, body weight gain, feed conversion, organ weights, serum thyroid, blood chemistries, macropathology, and histology data were collected. In Experiment 1, feed intake and body weight were diminished with no effect on feed conversion for the birds consuming diets containing CCWG-1-CuSO₄. In Experiment 2, feed intake and body weight were lower with no difference in feed conversion in birds fed diets containing greater than 2% CoverCress grain during d 0 to 28. During d 28 to 42 no difference in feed intake, body weight and an improvement in feed conversion was observed in birds fed all of the CoverCress grain products. In both experiments no significant negative effects were observed in processing, liver, kidney, and thyroid weights, T3, T4, blood chemistries, macropathology, and histopathology between the control and any of the CoverCress grain treatments. No difference in performance was observed in birds fed the mutant (4% CCWG-1) and gene-edited (4% CCWG-2) products. Pretreating CoverCress grain with copper sulfate did not have a significant effect on improving palatability. In conclusion, CoverCress grain can be safely fed to broilers when included at a target rate of 4% in diets and with total glucosinolate levels not to exceed 4.9 µmoles g^-1.

A lipidomics platform to analyze the fatty acid compositions of non-polar and polar lipid molecular species from plant tissues: Examples from developing seeds and seedlings of pennycress (Thlaspi arvense)

The lipidome comprises the total content of molecular species of each lipid class, and is measured using the analytical techniques of lipidomics. Many liquid chromatography-mass spectrometry (LC-MS) methods have previously been described to characterize the lipidome. However, many lipidomic approaches may not fully uncover the subtleties of lipid molecular species, such as the full fatty acid (FA) composition of certain lipid classes. Here, we describe a stepwise targeted lipidomics approach to characterize the polar and non-polar lipid classes using complementary LC-MS methods. Our "polar" method measures 260 molecular species across 12 polar lipid classes, and is performed using hydrophilic interaction chromatography (HILIC) on a NH2 column to separate lipid classes by their headgroup. Our "non-polar" method measures 254 molecular species across three non-polar lipid classes, separating molecular species on their FA characteristics by reverse phase (RP) chromatography on a C30 column. Five different extraction methods were compared, with an MTBE-based extraction chosen for the final lipidomics workflow. A state-of-the-art strategy to determine and relatively quantify the FA composition of triacylglycerols is also described. This lipidomics workflow was applied to developing, mature, and germinated pennycress seeds/seedlings and found unexpected changes among several lipid molecular species. During development, diacylglycerols predominantly contained long chain length FAs, which contrasted with the very long chain FAs of triacylglycerols in mature seeds. Potential metabolic explanations are discussed. The lack of very long chain fatty acids in diacylglycerols of germinating seeds may indicate very long chain FAs, such as erucic acid, are preferentially channeled into beta-oxidation for energy production.

Genetic and environmental drivers of large-scale epigenetic variation in Thlaspi arvense

Natural plant populations often harbour substantial heritable variation in DNA methylation. However, a thorough understanding of the genetic and environmental drivers of this epigenetic variation requires large-scale and high-resolution data, which currently exist only for a few model species. Here, we studied 207 lines of the annual weed Thlaspi arvense (field pennycress), collected across a large latitudinal gradient in Europe and propagated in a common environment. By screening for variation in DNA sequence and DNA methylation using whole-genome (bisulfite) sequencing, we found significant epigenetic population structure across Europe. Average levels of DNA methylation were strongly context-dependent, with highest DNA methylation in CG context, particularly in transposable elements and in intergenic regions. Residual DNA methylation variation within all contexts was associated with genetic variants, which often co-localized with annotated methylation machinery genes but also with new candidates. Variation in DNA methylation was also significantly associated with climate of origin, with methylation levels being lower in colder regions and in more variable climates. Finally, we used variance decomposition to assess genetic versus environmental associations with differentially methylated regions (DMRs). We found that while genetic variation was generally the strongest predictor of DMRs, the strength of environmental associations increased from CG to CHG and CHH, with climate-of-origin as the strongest predictor in about one third of the CHH DMRs. In summary, our data show that natural epigenetic variation in Thlaspi arvense is significantly associated with both DNA sequence and environment of origin, and that the relative importance of the two factors strongly depends on the sequence context of DNA methylation. T. arvense is an emerging biofuel and winter cover crop; our results may hence be relevant for breeding efforts and agricultural practices in the context of rapidly changing environmental conditions.

Genetic dissection of seed characteristics in field pennycress via genome-wide association mapping studies

Field pennycress (Thlaspi arvense L.) is a new winter annual cash cover crop with high oil content and seed yield, excellent winter hardiness, early maturation, and resistance to most pests and diseases. It provides living cover on fallow croplands between summer seasons, and in doing so reduces nutrient leaching into water sources, mitigates soil erosion, and suppresses weed growth. The first ever genome-wide association study (GWAS) was conducted on a pennycress diversity panel to identify marker trait associations with important seed size and composition related traits. The entire population was phenotyped in three total environments over 2 yr, and seed area, length, width, thousand grain weight, total oil, and total protein were measured post-harvest with specialized high-throughput imaging and near-infrared spectroscopy. Basic unbiased linear prediction values were calculated for each trait. Seed size traits tended to have higher entry mean reliabilities (0.76-0.79) compared with oil content (0.51) and protein content (0.37). Genotyping-by-sequencing identified 33,606 high quality genome-wide single nucleotide polymorphism (SNPs) that were coupled with phenotypic data to perform GWAS for seed area, length, width, thousand grain weight, total oil, and total protein content. Fifty-nine total marker-trait associations were identified revealing genomic regions controlling each trait. The significant SNPs explained 0.06-0.18% of the total variance for that trait in our population. A list of candidate genes was identified based on their functional annotations and characterization in other species. Our results confirm that GWAS is an efficient strategy to identify significant marker-trait associations that can be incorporated into marker-assisted selection pipelines to accelerate pennycress breeding progress.

Improved and Highly Efficient Agrobacterium rhizogenes-Mediated Genetic Transformation Protocol: Efficient Tools for Functional Analysis of Root-Specific Resistance Genes for Solanum lycopersicum cv. Micro-Tom

Gene function analysis, molecular breeding, and the introduction of new traits in crop plants all require the development of a high-performance genetic transformation system. In numerous crops, including tomatoes, Agrobacterium-mediated genetic transformation is the preferred method. As one of our ongoing research efforts, we are in the process of mapping a broad-spectrum nematode resistance gene (Me1) in pepper. We work to transform tomato plants with candidate genes to confer resistance to nematodes in Solanaceae members. The transformation technology development is designed to produce a reproducible, rapid, and highly effective Agrobacterium-mediated genetic transformation system of Micro-Tom. In our system, a transformation efficiency of over 90% was achieved. The entire procedure, starting from the germination of seeds to the establishment of transformed plants in soil, was completed in 53 days. We confirmed the presence of the NeoR/KanR and DsRed genes in the transformed roots by polymerase chain reaction. The hairy root plants were infected with nematodes, and after 3 months, the presence of DsRed and NeoR/KanR genes was detected in the transformant roots to confirm the long-term effectiveness of the method. The presented study may facilitate root-related research and exploration of root–pathogen interactions.

Chromosome-level Thlaspi arvense genome provides new tools for translational research and for a newly domesticated cash cover crop of the cooler climates

Thlaspi arvense (field pennycress) is being domesticated as a winter annual oilseed crop capable of improving ecosystems and intensifying agricultural productivity without increasing land use. It is a selfing diploid with a short life cycle and is amenable to genetic manipulations, making it an accessible field-based model species for genetics and epigenetics. The availability of a high-quality reference genome is vital for understanding pennycress physiology and for clarifying its evolutionary history within the Brassicaceae. Here, we present a chromosome-level genome assembly of var. MN106-Ref with improved gene annotation and use it to investigate gene structure differences between two accessions (MN108 and Spring32-10) that are highly amenable to genetic transformation. We describe non-coding RNAs, pseudogenes and transposable elements, and highlight tissue-specific expression and methylation patterns. Resequencing of forty wild accessions provided insights into genome-wide genetic variation, and QTL regions were identified for a seedling colour phenotype. Altogether, these data will serve as a tool for pennycress improvement in general and for translational research across the Brassicaceae.

Natural variation and improved genome annotation of the emerging biofuel crop field pennycress (Thlaspi arvense)

The Brassicaceae family comprises more than 3,700 species with a diversity of phenotypic characteristics, including seed oil content and composition. Recently, the global interest in Thlaspi arvense L. (pennycress) has grown as the seed oil composition makes it a suitable source for biodiesel and aviation fuel production. However, many wild traits of this species need to be domesticated to make pennycress ideal for cultivation. Molecular breeding and engineering efforts require the availability of an accurate genome sequence of the species. Here, we describe pennycress genome annotation improvements, using a combination of long- and short-read transcriptome data obtained from RNA derived from embryos of 22 accessions, in addition to public genome and gene expression information. Our analysis identified 27,213 protein-coding genes, as well as on average 6,188 biallelic SNPs. In addition, we used the identified SNPs to evaluate the population structure of our accessions. The data from this analysis support that the accession Ames 32872, originally from Armenia, is highly divergent from the other accessions, while the accessions originating from Canada and the United States cluster together. When we evaluated the likely signatures of natural selection from alternative SNPs, we found 7 candidate genes under likely recent positive selection. These genes are enriched with functions related to amino acid metabolism and lipid biosynthesis and highlight possible future targets for crop improvement efforts in pennycress.

Optimisation of root traits to provide enhanced ecosystem services in agricultural systems: A focus on cover crops

Roots are the interface between the plant and the soil and play a central role in multiple ecosystem processes. With intensification of agricultural practices, rhizosphere processes are being disrupted and are causing degradation of the physical, chemical and biotic properties of soil. However, cover crops, a group of plants that provide ecosystem services, can be utilised during fallow periods or used as an intercrop to restore soil health. The effectiveness of ecosystem services provided by cover crops varies widely as very little breeding has occurred in these species. Improvement of ecosystem service performance is rarely considered as a breeding trait due to the complexities and challenges of belowground evaluation. Advancements in root phenotyping and genetic tools are critical in accelerating ecosystem service improvement in cover crops. In this study, we provide an overview of the range of belowground ecosystem services provided by cover crop roots: (1) soil structural remediation, (2) capture of soil resources and (3) maintenance of the rhizosphere and building of organic matter content. Based on the ecosystem services described, we outline current and promising phenotyping technologies and breeding strategies in cover crops that can enhance agricultural sustainability through improvement of root traits.

Quantitative trait locus mapping combined with variant and transcriptome analyses identifies a cluster of gene candidates underlying the variation in leaf wax between upland and lowland switchgrass ecotypes

Mapping combined with expression and variant analyses in switchgrass, a crop with complex genetics, identified a cluster of candidate genes for leaf wax in a fast-evolving region of chromosome 7K. Switchgrass (Panicum virgatum L.) is a promising warm-season candidate energy crop. It occurs in two ecotypes, upland and lowland, which vary in a number of phenotypic traits, including leaf glaucousness. To initiate trait mapping, two F₂ mapping populations were developed by crossing two different F₁ sibs derived from a cross between the tetraploid lowland genotype AP13 and the tetraploid upland genotype VS16, and high-density linkage maps were generated. Quantitative trait locus (QTL) analyses of visually scored leaf glaucousness and of hydrophobicity of the abaxial leaf surface measured using a drop shape analyzer identified highly significant colocalizing QTL on chromosome 7K (Chr07K). Using a multipronged approach, we identified a cluster of genes including Pavir.7KG077009, which encodes a Type III polyketide synthase-like protein, and Pavir.7KG013754 and Pavir.7KG030500, two highly similar genes that encode putative acyl-acyl carrier protein (ACP) thioesterases, as strong candidates underlying the QTL. The lack of homoeologs for any of the three genes on Chr07N, the relatively low level of identity with other switchgrass KCS proteins and thioesterases, as well as the organization of the surrounding region suggest that Pavir.7KG077009 and Pavir.7KG013754/Pavir.7KG030500 were duplicated into a fast-evolving chromosome region, which led to their neofunctionalization. Furthermore, sequence analyses showed all three genes to be absent in the two upland compared to the two lowland accessions analyzed. This study provides an example of and practical guide for trait mapping and candidate gene identification in a complex genetic system by combining QTL mapping, transcriptomics and variant analysis.

Technologies enabling rapid crop improvements for sustainable agriculture: example pennycress (Thlaspi arvense L.)

Growing concerns over food insecurity and ecosystems health related to population growth and climate change have challenged scientists to develop new crops, employing revolutionary technologies in combination with traditional methods. In this review, we discuss the domestication of the oilseed-producing cover crop pennycress, which along with the development of other new crops and improvements to farming practices can provide sustainable solutions to address malnutrition and environmental impacts of production agriculture. We highlight some of the new technologies such as bioinformatics-enabled next-generation sequencing and CRISPR genome editing in combination with traditional mutation breeding that has accelerated pennycress development as a new crop and a potential model system. Furthermore, we provide a brief overview of the technologies that can be integrated for improving pennycress and other crops and the status of pennycress development using these technologies.

Open Inquiry versus Broadly Relevant Short-Term Research Experiences for Non-Biology Majors

Undergraduate student participation in course-based research experiences results in many positive outcomes, but there is a lack of evidence demonstrating which elements of a research experience are necessary, especially for non-biology majors. Broad relevance is one element that can be logistically challenging to incorporate into research experiences in large-enrollment courses. We investigated the impacts of broad relevance in a short-term research experience in an introductory biology course for non-majors. Students either participated in an open-inquiry research experience (OI-RE), where they developed their own research question, or a broadly relevant research experience (BR-RE), where they investigated a question assigned to them that was relevant to an ongoing research project. We found a significant association between the type of research project experienced and students' preference for an experience, with half of the students in the OI-RE group and nearly all students in the BR-RE group preferring a broadly relevant research experience. However, since science confidence increased over the course for both groups, these findings indicate that while students who participated in a BR-RE valued it, broadly relevant research experiences may not be necessary for positive outcomes for non-majors.

CRISPR/Cas9-Induced fad2 and rod1 Mutations Stacked With fae1 Confer High Oleic Acid Seed Oil in Pennycress (Thlaspi arvense L.)

Pennycress (Thlaspi arvense L.) is being domesticated as an oilseed cash cover crop to be grown in the off-season throughout temperate regions of the world. With its diploid genome and ease of directed mutagenesis using molecular approaches, pennycress seed oil composition can be rapidly tailored for a plethora of food, feed, oleochemical and fuel uses. Here, we utilized Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 technology to produce knockout mutations in the FATTY ACID DESATURASE2 (FAD2) and REDUCED OLEATE DESATURATION1 (ROD1) genes to increase oleic acid content. High oleic acid (18:1) oil is valued for its oxidative stability that is superior to the polyunsaturated fatty acids (PUFAs) linoleic (18:2) and linolenic (18:3), and better cold flow properties than the very long chain fatty acid (VLCFA) erucic (22:1). When combined with a FATTY ACID ELONGATION1 (fae1) knockout mutation, fad2 fae1 and rod1 fae1 double mutants produced ∼90% and ∼60% oleic acid in seed oil, respectively, with PUFAs in fad2 fae1 as well as fad2 single mutants reduced to less than 5%. MALDI-MS spatial imaging analyses of phosphatidylcholine (PC) and triacylglycerol (TAG) molecular species in wild-type pennycress embryo sections from mature seeds revealed that erucic acid is highly enriched in cotyledons which serve as storage organs, suggestive of a role in providing energy for the germinating seedling. In contrast, PUFA-containing TAGs are enriched in the embryonic axis, which may be utilized for cellular membrane expansion during seed germination and seedling emergence. Under standard growth chamber conditions, rod1 fae1 plants grew like wild type whereas fad2 single and fad2 fae1 double mutant plants exhibited delayed growth and overall reduced heights and seed yields, suggesting that reducing PUFAs below a threshold in pennycress had negative physiological effects. Taken together, our results suggest that combinatorial knockout of ROD1 and FAE1 may be a viable route to commercially increase oleic acid content in pennycress seed oil whereas mutations in FAD2 will likely require at least partial function to avoid fitness trade-offs.

Functional analysis of β-ketoacyl-CoA synthase from biofuel feedstock Thlaspi arvense reveals differences in the triacylglycerol biosynthetic pathway among Brassicaceae

Differences in FAE1 enzyme affinity for the acyl-CoA substrates, as well as the balance between the different pathways involved in their incorporation to triacylglycerol might be determinant of the different composition of the seed oil in Brassicaceae. Brassicaceae present a great heterogeneity of seed oil and fatty acid composition, accumulating Very Long Chain Fatty Acids with industrial applications. However, the molecular determinants of these differences remain elusive. We have studied the β-ketoacyl-CoA synthase from the high erucic feedstock Thlaspi arvense (Pennycress). Functional characterization of the Pennycress FAE1 enzyme was performed in two Arabidopsis backgrounds; Col-0, with less than 2.5% of erucic acid in its seed oil and the fae1-1 mutant, deficient in FAE1 activity, that did not accumulate erucic acid. Seed-specific expression of the Pennycress FAE1 gene in Col-0 resulted in a 3 to fourfold increase of erucic acid content in the seed oil. This increase was concomitant with a decrease of eicosenoic acid levels without changes in oleic ones. Interestingly, only small changes in eicosenoic and erucic acid levels occurred when the Pennycress FAE1 gene was expressed in the fae1-1 mutant, with high levels of oleic acid available for elongation, suggesting that the Pennycress FAE1 enzyme showed higher affinity for eicosenoic acid substrates, than for oleic ones in Arabidopsis. Erucic acid was incorporated to triacylglycerol in the transgenic lines without significant changes in their levels in the diacylglycerol fraction, suggesting that erucic acid was preferentially incorporated to triacylglycerol via DGAT1. Expression analysis of FAE1, AtDGAT1, AtLPCAT1 and AtPDAT1 genes in the transgenic lines further supported this conclusion. Differences in FAE1 affinity for the oleic and eicosenoic substrates among Brassicaceae, as well as their incorporation to triacylglycerol might explain the differences in composition of their seed oil.

Pennycress, carbon wise: labeling experiments reveal how pennycress seeds efficiently incorporate carbon into biomass

This article comments on: Tsogtbaatar E, Cocuron J-C, Alonso AP. 2020. Non-conventional pathways enable pennycress (Thlaspi arvense L.) embryos to achieve high biosynthetic efficiency. Journal of Experimental Botany 71, 3037–3051.

Genetic Diversity of Field Pennycress (Thlaspi arvense) Reveals Untapped Variability and Paths Toward Selection for Domestication

Evaluation of genetic diversity within wild populations is an essential process for improvement and domestication of new crop species. This process involves evaluation of population structure and individual accessions based on genetic markers, growth habits, and geographic collection area. In this study, accessions of field pennycress were analyzed to identify population structure and variation in germplasm available for breeding. A total of 9157 genome-wide single nucleotide polymorphisms (SNPs) were identified among the 121 accessions analyzed, and linkage disequilibrium based pruning resulted in 3497 SNPs. Bayesian cluster analysis was implemented in STRUCTURE v2.3.4 to identify four population groups. These groups were confirmed based on principal components analysis and geographic origins. Pairwise diversity among accessions was evaluated and revealed considerable genetic variation. Notably, a subset of accessions from Armenia with exceptional genetic variation was identified. This survey is the first to report significant genetic diversity among pennycress accessions and explain some of the phenotypic differences previously observed in the germplasm. Understanding variation in pennycress accessions will be a crucial step for selection, breeding, and domestication of a new cash cover crop for cold climates.