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

Transcriptomic and proteomic analysis of oil body associated protein dynamics in the biofuel feedstock Pennycress (Thlaspi arvense)

Pennycress (Thlaspi arvense) is an emerging feedstock for biofuel production because of its high seed oil content enriched in erucic acid. A combination of transcriptomic and proteomic tools was used to characterize the dynamics and relative abundance of the major oil body related proteins in the Pennycress seed. Our analysis identified 21 oleosins (OLE), 6 oil body associated proteins (OBAPs), 3 SEIPINS, 3 caleosins, 3 stereolisins and 3 lipid droplet associated proteins (LDAPs) in the Pennycress genome, showing high homology with respect to Arabidopsis or rapeseed. RNA-Seq analysis on five Pennycress seed maturation stages showed that most OLE and OBAP genes increased their expression with seed maturation, coinciding with the highest accumulation of triacylglycerol. Western-blot analysis of the OLE2 protein during seed maturation confirmed this result. However, OLE5 and SEIPIN1 genes showed higher expression at the early stages of seed maturation, suggesting that both proteins could be particularly involved in the initial stages of oil body formation. Proteomic analysis on oil body enriched fractions from the YELLOW and MATURE late seed maturation stages showed that all oleosin proteins were highly abundant in oil bodies. Caleosins and stereolisins were also highly abundant. Our results indicate the existence of differential expression patterns of oil body related genes during Pennycress seed maturation, suggesting different roles of these proteins for the formation and stabilization of oil bodies in the Pennycress seed.

Propelling sustainable energy: Multi-omics analysis of pennycress FATTY ACID ELONGATION1 knockout for biofuel production

The aviation industry's growing interest in renewable jet fuel has encouraged the exploration of alternative oilseed crops. Replacing traditional fossil fuels with a sustainable, domestically sourced crop can substantially reduce carbon emissions, thus mitigating global climate instability. Pennycress (Thlaspi arvense L.) is an emerging oilseed intermediate crop that can be grown during the offseason between maize (Zea mays) and soybean (Glycine max) to produce renewable biofuel. Pennycress is being domesticated through breeding and mutagenesis, providing opportunities for trait enhancement. Here, we employed metabolic engineering strategies to improve seed oil composition and bolster the plant's economic competitiveness. FATTY ACID ELONGATION1 (FAE1) was targeted using CRISPR-Cas 9 gene editing to eliminate very long chain fatty acids (VLCFAs) from pennycress seed oil, thereby enhancing its cold flow properties. Through an integrated multiomics approach, we investigated the impact of eliminating VLCFAs in developing and mature plant embryos. Our findings revealed improved cold-germination efficiency in fae1, with seedling emergence occurring up to 3 d earlier at 10 °C. However, these alterations led to a tradeoff between storage oil content and composition. Additionally, these shifts in lipid biosynthesis were accompanied by broad metabolic changes, such as the accumulation of glucose and ADP-glucose quantities consistent with increased starch production. Furthermore, shifts to shorter FA chains triggered the upregulation of heat shock proteins, underscoring the importance of VLCFAs in stress signaling pathways. Overall, this research provides crucial insights for optimizing pennycress seed oil while preserving essential traits for biofuel applications.

Transcriptional engineering for value enhancement of oilseed crops: a forward perspective

Plant-derived oils provide 20%-35% of dietary calories and are a primary source of essential omega-6 (linoleic) and omega-3 (α-linolenic) fatty acids. While traditional breeding has significantly increased yields in key oilseed crops like soybean, sunflower, canola, peanut, and cottonseed, overall gains have plateaued over the past few decades. Oilseed crops also experience substantial yield losses in both prime and marginal agricultural areas due to biotic and abiotic stresses and shifting agro-climates. Recent genomic, transcriptomic, and metabolomics research has expanded our understanding of the genetic and physiological control of fatty acid biosynthesis and composition. Many oilseed species have inherent stress-combating mechanisms, including transcription factor regulation. Advances in genome editing tools like CRISPR/Cas9 offer precise genetic modifications, targeting transcription factors and binding sites to enhance desirable traits, such as the nutritional profile and chemical composition of fatty acids. This review explores the application of genome editing in oilseed improvement, covering recent progress, challenges, and future potential to boost yield and oil content. These advancements could play a transformative role in developing resilient, nutritious crop varieties essential for sustainable food security in a changing climate.

Different orthology inference algorithms generate similar predicted orthogroups among Brassicaceae species

Premise

Orthology inference is crucial for comparative genomics, and multiple algorithms have been developed to identify putative orthologs for downstream analyses. Despite the abundance of proposed solutions, including publicly available benchmarks, it is difficult to assess which tool is most suitable for plant species, which commonly have complex genomic histories.

Methods

We explored the performance of four orthology inference algorithms-OrthoFinder, SonicParanoid, Broccoli, and OrthNet-on eight Brassicaceae genomes in two groups: one group comprising only diploids and another set comprising the diploids, two mesopolyploids, and one recent hexaploid genome.

Results

The composition of the orthogroups reflected the species' ploidy and genomic histories, with the diploid set having a higher proportion of identical orthogroups. While the diploid + higher ploidy set had a lower proportion of orthogroups with identical compositions, the average degree of similarity between the orthogroups was not different from the diploid set.

Discussion

Three algorithms-OrthoFinder, SonicParanoid, and Broccoli-are helpful for initial orthology predictions. Results produced using OrthNet were generally outliers but could still provide detailed information about gene colinearity. With our Brassicaceae dataset, slight discrepancies were found across the orthology inference algorithms, necessitating additional analyses such as tree inference to fine-tune results.

Discarded sequencing reads uncover natural variation in pest resistance in Thlaspi arvense

Understanding the genomic basis of natural variation in plant pest resistance is an important goal in plant science, but it usually requires large and labor-intensive phenotyping experiments. Here, we explored the possibility that non-target reads from plant DNA sequencing can serve as phenotyping proxies for addressing such questions. We used data from a whole-genome and -epigenome sequencing study of 207 natural lines of field pennycress (Thlaspi arvense) that were grown in a common environment and spontaneously colonized by aphids, mildew, and other microbes. We found that the numbers of non-target reads assigned to the pest species differed between populations, had significant SNP-based heritability, and were associated with climate of origin and baseline glucosinolate contents. Specifically, pennycress lines from cold and thermally fluctuating habitats, presumably less favorable to aphids, showed higher aphid DNA load, i.e., decreased aphid resistance. Genome-wide association analyses identified genetic variants at known defense genes but also novel genomic regions associated with variation in aphid and mildew DNA load. Moreover, we found several differentially methylated regions associated with pathogen loads, in particular differential methylation at transposons and hypomethylation in the promoter of a gene involved in stomatal closure, likely induced by pathogens. Our study provides first insights into the defense mechanisms of Thlaspi arvense, a rising crop and model species, and demonstrates that non-target whole-genome sequencing reads, usually discarded, can be leveraged to estimate intensities of plant biotic interactions. With rapidly increasing numbers of large sequencing datasets worldwide, this approach should have broad application in fundamental and applied research.

Targeted modulation of pennycress lipid droplet proteins impacts droplet morphology and seed oil content

Lipid droplets (LDs) are unusual organelles that have a phospholipid monolayer surface and a hydrophobic matrix. In oilseeds, this matrix is nearly always composed of triacylglycerols (TGs) for efficient storage of carbon and energy. Various proteins play a role in their assembly, stability and turnover, and even though the major structural oleosin proteins in seed LDs have been known for decades, the factors influencing LD formation and dynamics are still being uncovered mostly in the "model oilseed" Arabidopsis. Here we identified several key LD biogenesis proteins in the seeds of pennycress, a potential biofuel crop, that were correlated previously with seed oil content and characterized here for their participation in LD formation in transient expression assays and stable transgenics. One pennycress protein, the lipid droplet associated protein-interacting protein (LDIP), was able to functionally complement the Arabidopsis ldip mutant, emphasizing the close conservation of lipid storage among these two Brassicas. Moreover, loss-of-function ldip mutants in pennycress exhibited increased seed oil content without compromising plant growth, raising the possibility that LDIP or other LD biogenesis factors may be suitable targets for improving yields in oilseed crops more broadly.

Domestication and engineering of pennycress (Thlaspi arvense L.): challenges and opportunities for sustainable bio-based feedstocks

MAIN CONCLUSION

Pennycress, as an emerging oilseed crop with high oil content, faces challenges but offers potential for sustainable bioproducts; ongoing research aims to enhance its traits and quality. Pennycress (Thlaspi arvense L.) is an emerging oilseed crop with many advantages, such as high seed oil (27-39%) and monounsaturated fatty acid (55.6%) content, making it an attractive candidate to produce sustainable bioproducts. However, several challenges are associated with domesticating pennycress, including high silicle shatter, which reduces seed yield during harvest, non-uniformed germination rate and high contents of erucic acid and glucosinolates, which have adverse health effects on humans and animals. Pennycress, which can be easily and rapidly transformed using the floral dip method under vacuum, can achieve trait improvements. Ongoing research for pennycress domestication using mutation breeding, including ethylmethylsulfonate treatment and genome editing, aims to improve its quality. Pennycress can be used as an excellent platform for producing industrially important fatty acids such as hydroxy and epoxy fatty acids and docosahexaenoic acid. In conclusion, pennycress is a promising oilseed crop with multiple advantages and potential applications. Continuous improvements in quality and engineering for producing high-value bio-based feedstocks in pennycress will establish it as a sustainable and economically valuable crop.

High-resolution methylome analysis uncovers stress-responsive genomic hotspots and drought-sensitive transposable element superfamilies in the clonal Lombardy poplar

DNA methylation is environment-sensitive and can mediate stress responses. In trees, changes in the environment might cumulatively shape the methylome landscape over time. However, because high-resolution methylome studies usually focus on single environmental cues, the stress-specificity and long-term stability of methylation responses remain unclear. Here, we studied the methylome plasticity of a Populus nigra cv. 'Italica' clone widely distributed across Europe. Adult trees from different geographic locations were clonally propagated in a common garden experiment and exposed to cold, heat, drought, herbivory, rust infection, and salicylic acid treatments. Whole-genome bisulfite sequencing revealed stress-induced and naturally occurring DNA methylation variants. In CG/CHG contexts, the same genomic regions were often affected by multiple stresses, suggesting a generic methylome response. Moreover, these variants showed striking overlap with naturally occurring methylation variants between trees from different locations. Drought treatment triggered CHH hypermethylation of transposable elements, affecting entire superfamilies near drought-responsive genes. Thus, we revealed genomic hotspots of methylation change that are not stress-specific and that contribute to natural DNA methylation variation, and identified stress-specific hypermethylation of entire transposon superfamilies with possible functional consequences. Our results underscore the importance of studying multiple stressors in a single experiment for recognizing general versus stress-specific methylome responses.

Current status of community resources and priorities for weed genomics research

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

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.

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

The CLAVATA3/ESR-related peptide family in the biofuel crop pennycress

CLAVATA3/ESR-related (CLE) peptides perform a variety of important functions in plant development and historically have been targeted during the domestication of existing crops. Pennycress (Thlaspi arvense) is an emerging biofuel crop currently undergoing domestication that offers novel monetary and environmental incentives as a winter cover crop during an otherwise fallow period of the corn/soybean farming rotation. Here we report the characterization of the CLE gene family in pennycress through homology comparison of the CLE motif with other dicot species by conducting a homology comparison and maximum likelihood phylogenetic analysis supplemented with manual annotation. Twenty-seven pennycress CLE genes were identified, and their expression analyzed through transcriptome profiling and RT-qPCR. Our study provides a genome-wide analysis of the CLE gene family in pennycress and carries significant value for accelerating the domestication of this crop through identification of potential key developmental regulatory genes.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

Weed genomics: yielding insights into the genetics of weedy traits for crop improvement

Weeds cause tremendous economic and ecological damage worldwide. The number of genomes established for weed species has sharply increased during the recent decade, with some 26 weed species having been sequenced and de novo genomes assembled. These genomes range from 270 Mb (Barbarea vulgaris) to almost 4.4 Gb (Aegilops tauschii). Importantly, chromosome-level assemblies are now available for 17 of these 26 species, and genomic investigations on weed populations have been conducted in at least 12 species. The resulting genomic data have greatly facilitated studies of weed management and biology, especially origin and evolution. Available weed genomes have indeed revealed valuable weed-derived genetic materials for crop improvement. In this review, we summarize the recent progress made in weed genomics and provide a perspective for further exploitation in this emerging field.

Sinapis genomes provide insights into whole-genome triplication and divergence patterns within tribe Brassiceae

Sinapis alba and Sinapis arvensis are mustard crops within the Brassiceae tribe of the Brassicaceae family, and represent an important genetic resource for crop improvement. We performed the de novo assembly of Brassica nigra, S. alba, and S. arvensis, and conducted comparative genomics to investigate the pattern of genomic evolution since an ancient whole-genome triplication event. Both Sinapis species retained evidence of the Brassiceae whole-genome triplication approximately 20.5 million years ago (Mya), with subgenome dominance observed in gene density, gene expression, and selective constraint. While S. alba diverged from the ancestor of Brassica and Raphanus at approximately 12.5 Mya, the divergence time of S. arvensis and B. nigra was approximately 6.5 Mya. S. arvensis and B. nigra had greater collinearity compared with their relationship to either Brassica rapa or Brassica oleracea. Two chromosomes of S. alba (Sal03 and Sal08) were completely collinear with two ancestral chromosomes proposed in the Ancestral Crucifer Karyotype (ACK) genomic block model, the first time this has been observed in the Brassiceae. These results are consistent with S. alba representing a relatively ancient lineage of the species evolved from the common ancestor of tribe Brassiceae, and suggest that the phylogeny of the Brassica and Sinapis genera requires some revision. Our study provides new insights into the genome evolution and phylogenetic relationships of Brassiceae and provides genomic information for genetic improvement of these plants.

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.

DNA methylation dynamics during stress response in woodland strawberry (Fragaria vesca)

Environmental stresses can result in a wide range of physiological and molecular responses in plants. These responses can also impact epigenetic information in genomes, especially at the level of DNA methylation (5-methylcytosine). DNA methylation is the hallmark heritable epigenetic modification and plays a key role in silencing transposable elements (TEs). Although DNA methylation is an essential epigenetic mechanism, fundamental aspects of its contribution to stress responses and adaptation remain obscure. We investigated epigenome dynamics of wild strawberry (Fragaria vesca) in response to variable ecologically relevant environmental conditions at the DNA methylation level. F. vesca methylome responded with great plasticity to ecologically relevant abiotic and hormonal stresses. Thermal stress resulted in substantial genome-wide loss of DNA methylation. Notably, all tested stress conditions resulted in marked hot spots of differential DNA methylation near centromeric or pericentromeric regions, particularly in the non-symmetrical DNA methylation context. Additionally, we identified differentially methylated regions (DMRs) within promoter regions of transcription factor (TF) superfamilies involved in plant stress-response and assessed the effects of these changes on gene expression. These findings improve our understanding on stress-response at the epigenome level by highlighting the correlation between DNA methylation, TEs and gene expression regulation in plants subjected to a broad range of environmental stresses.

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.

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.