Exploring Camelina sativa lipid metabolism regulation by combining gene co-expression and DNA affinity purification analyses

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 laboratory and field scientists alike to take advantage of the vast Arabidopsis datasets, annotations, germplasm, constructs, methods, and 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.

Specificity landscapes of 40 R2R3-MYBs reveal how paralogs target different cis-elements by homodimeric binding

Paralogous transcription factors (TFs) frequently recognize highly similar DNA motifs. Homodimerization can help distinguish them according to their different dimeric configurations. Here, by studying R2R3-MYB TFs, we show that homodimerization can also directly change the recognized DNA motifs to distinguish between similar TFs. By high-throughput SELEX, we profiled the specificity landscape for 40 R2R3-MYBs of subfamily VIII and curated 833 motif models. The dimeric models show that homodimeric binding has evoked specificity changes for AtMYBs. Focusing on AtMYB2 as an example, we show that homodimerization has modified its specificity and allowed it to recognize additional cis-regulatory sequences that are different from the closely related CCWAA-box AtMYBs and are unique among all AtMYBs. Genomic sites described by the modified dimeric specificities of AtMYB2 are conserved in evolution and involved in AtMYB2-specific transcriptional activation. Collectively, this study provides rich data on sequence preferences of VIII R2R3-MYBs and suggests an alternative mechanism that guides closely related TFs to respective cis-regulatory sites.

Born of frustration: the emergence of Camelina sativa as a platform for lipid biotechnology

The emerging crop Camelina sativa (L.) Crantz (camelina) is a Brassicaceae oilseed with a rapidly growing reputation for the deployment of advanced lipid biotechnology and metabolic engineering. Camelina is recognized by agronomists for its traits, including yield, oil/protein content, drought tolerance, limited input requirements, plasticity, and resilience. Its utility as a platform for metabolic engineering was then quickly recognized, and biotechnologists have benefited from its short life cycle and facile genetic transformation, producing numerous transgenic interventions to modify seed lipid content and generate novel products. The desire to work with a plant that is both a model and crop has driven the expansion of research resources for camelina, including increased availability of genome and other -omics data sets. Collectively, the expansion of these resources has established camelina as an ideal plant to study the regulation of lipid metabolism and genetic improvement. Furthermore, the unique characteristics of camelina enables the design-build-test-learn cycle to be transitioned from the controlled environment to the field. Complex metabolic engineering to synthesize and accumulate high levels of novel fatty acids and modified oils in seeds can be deployed, tested, and undergo rounds of iteration in agronomically relevant environments. Engineered camelina oils are now increasingly being developed and used to sustainably supply improved nutrition, feed, biofuels, and fossil fuel replacements for high-value chemical products. In this review, we provide a summary of seed fatty acid synthesis and oil assembly in camelina, highlighting how discovery research in camelina supports the advance of metabolic engineering toward the predictive manipulation of metabolism to produce desirable bio-based products. Further examples of innovation in camelina seed lipid engineering and crop improvement are then provided, describing how technologies (e.g. genetic modification [GM], gene editing [GE], RNAi, alongside GM and GE stacking) can be applied to produce new products and denude undesirable traits. Focusing on the production of long chain polyunsaturated omega-3 fatty acids in camelina, we describe how lipid biotechnology can transition from discovery to a commercial prototype. The prospects to produce structured triacylglycerol with fatty acids in specified stereospecific positions are also discussed, alongside the future outlook for the agronomic uptake of camelina lipid biotechnology.

Unveiling the camelina MBOAT gene family: Phylogenetic insights and regulatory landscape

The membrane-bound O-acyltransferase (MBOAT) gene family comprises enzymes responsible for transferring acyl groups to various lipid molecules. Some members of the MBOAT gene family and their functions have been extensively studied in the model plant Arabidopsis. However, research on the MBOAT gene family in camelina is still limited. In this study, 54 MBOATs were identified on 17 chromosomes and one unidentified scaffold in camelina, including seven newly identified genes. A total of 149 MBOATs were identified in 10 other species. Six subgroups of these MBOATs with different conservation were classified by phylogenetic analysis, showing diversification between monocots and dicots. Detailed analysis of the motif composition, evolutionary relationships, and gene structures of CsaMBOATs are presented. The results of the syntenic analysis suggest that the evolution of CsaMBOAT gene family is primarily driven by segmental and tandem duplications, and that there is a stronger collinearity within dicots. In addition, analysis of CsaMBOAT gene promoter cis-elements reveals a possible transcriptional regulation and tissue-specific expression, highlighting potential role in plant stress responses and hormone signaling. Furthermore, both the transcriptome and RT-qPCR data revealed the changes in the expression levels of DGAT1 during salt stress treatment. Subsequent analyses indicated that DGAT1 influenced the ratio of fatty acid fractions in the plants. Importantly, a large number of transcription factors involved in the regulation of CsaMBOAT gene expression were identified by WGCNA analysis, and the transcriptional data confirmed that the NAC032 and CAMMTA6 genes play a role upstream of DGAT1. This study not only identified the members of the MBOAT in camelina, but also provided insights and clues into its regulatory mechanisms.

Exploring genetic diversity, population structure, and subgenome differences in the allopolyploid Camelina sativa: implications for future breeding and research studies

Camelina (Camelina sativa), an allohexaploid species, is an emerging aviation biofuel crop that has been the focus of resurgent interest in recent decades. To guide future breeding and crop improvement efforts, the community requires a deeper comprehension of subgenome dominance, often noted in allopolyploid species, "alongside an understanding of the genetic diversity" and population structure of material present within breeding programs. We conducted population genetic analyses of a C. sativa diversity panel, leveraging a new genome, to estimate nucleotide diversity and population structure, and analyzed for patterns of subgenome expression dominance among different organs. Our analyses confirm that C. sativa has relatively low genetic diversity and show that the SG3 subgenome has substantially lower genetic diversity compared to the other two subgenomes. Despite the low genetic diversity, our analyses identified 13 distinct subpopulations including two distinct wild populations and others putatively representing founders in existing breeding populations. When analyzing for subgenome composition of long non-coding RNAs, which are known to play important roles in (a)biotic stress tolerance, we found that the SG3 subgenome contained significantly more lincRNAs compared to other subgenomes. Similarly, transcriptome analyses revealed that expression dominance of SG3 is not as strong as previously reported and may not be universal across all organ types. From a global analysis, SG3 "was only significant higher expressed" in flower, flower bud, and fruit organs, which is an important discovery given that the crop yield is associated with these organs. Collectively, these results will be valuable for guiding future breeding efforts in camelina.

Molecular mechanisms underlying gene regulatory variation of maize metabolic traits

Variation in gene expression levels is pervasive among individuals and races or varieties, and has substantial agronomic consequences, for example, by contributing to hybrid vigor. Gene expression level variation results from mutations in regulatory sequences (cis) and/or transcription factor (TF) activity (trans), but the mechanisms underlying cis- and/or trans-regulatory variation of complex phenotypes remain largely unknown. Here, we investigated gene expression variation mechanisms underlying the differential accumulation of the insecticidal compounds maysin and chlorogenic acid in silks of widely used maize (Zea mays) inbreds, B73 and A632. By combining transcriptomics and cistromics, we identified 1,338 silk direct targets of the maize R2R3-MYB TF Pericarp color1 (P1), consistent with it being a regulator of maysin and chlorogenic acid biosynthesis. Among these P1 targets, 464 showed allele-specific expression (ASE) between B73 and A632 silks. Allelic DNA-affinity purification sequencing identified 34 examples in which P1 allelic specific binding (ASB) correlated with cis-expression variation. From previous yeast one-hybrid studies, we identified 9 TFs potentially implicated in the control of P1 targets, with ASB to 83 out of 464 ASE genes (cis) and differential expression of 4 out of 9 TFs between B73 and A632 silks (trans). These results provide a molecular framework for understanding universal mechanisms underlying natural variation of gene expression levels, and how the regulation of metabolic diversity is established.

Prioritizing Maize Metabolic Gene Regulators through Multi-Omic Network Integration

Elucidating gene regulatory networks is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. Here, we integrated 46 co-expression networks, 283 protein-DNA interaction (PDI) assays, and 16 million SNPs used to identify expression quantitative trait loci (eQTL) to construct TF-target networks. In total, we analyzed ∼4.6M interactions to generate four distinct types of TF-target networks: co-expression, PDI, trans -eQTL, and cis -eQTL combined with PDIs. To functionally annotate TFs based on their target genes, we implemented three different network integration strategies. We evaluated the effectiveness of each strategy through TF loss-of function mutant inspection and random network analyses. The multi-network integration allowed us to identify transcriptional regulators of several biological processes. Using the topological properties of the fully integrated network, we identified potential functionally redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

GRAPHICAL ABSTRACT

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

Functional investigation of five R2R3-MYB transcription factors associated with wood development in Eucalyptus using DAP-seq-ML

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta.

"Simple Tidy GeneCoEx": A gene co-expression analysis workflow powered by tidyverse and graph-based clustering in R

Gene co-expression analysis is an effective method to detect groups (or modules) of co-expressed genes that display similar expression patterns, which may function in the same biological processes. Here, we present "Simple Tidy GeneCoEx", a gene co-expression analysis workflow written in the R programming language. The workflow is highly customizable across multiple stages of the pipeline including gene selection, edge selection, clustering resolution, and data visualization. Powered by the tidyverse package ecosystem and network analysis functions provided by the igraph package, the workflow detects gene co-expression modules whose members are highly interconnected. Step-by-step instructions with two use case examples as well as source code are available at https://github.com/cxli233/SimpleTidy_GeneCoEx.

The coexpression of two desaturases provides an optimized reduction of saturates in camelina oil

Reducing the saturate content of vegetable oils is key to increasing their utility and adoption as a feedstock for the production of biofuels. Expression of either the FAT5 16 : 0-CoA desaturase from Caenorhabditis elegans, or an engineered cyanobacterial 16 : 0/18 : 0-glycerolipid desaturase, DES9*, in seeds of Arabidopsis (Arabidopsis thaliana) substantially lowered oil saturates. However, because pathway fluxes and regulation of oil synthesis are known to differ across species, translating this transgene technology from the model plant to crop species requires additional investigation. In the work reported here, we found that high expression of FAT5 in seeds of camelina (Camelina sativa) provided only a moderate decrease in saturates, from 12.9% of total oil fatty acids in untransformed controls to 8.6%. Expression of DES9* reduced saturates to 4.6%, but compromised seed physiology and oil content. However, the coexpression of the two desaturases together cooperatively reduced saturates to only 4.0%, less than one-third of the level in the parental line, without compromising oil yield or seedling germination and establishment. Our successful lowering of oil saturates in camelina identifies strategies that can now be integrated with genetic engineering approaches that reduce polyunsaturates to provide optimized oil composition for biofuels in camelina and other oil seed crops.

Regulatory mechanism of MeGI on sexuality in Diospyros oleifera

Dioecy system is an important strategy for maintaining genetic diversity. The transcription factor MeGI, contributes to dioecy by promoting gynoecium development in Diospyros lotus and D. kaki. However, the function of MeGI in D. oleifera has not been identified. In this study, we confirmed that MeGI, cloned from D. oleifera, repressed the androecium development in Arabidopsis thaliana. Subsequently, chromatin immunoprecipitation-sequencing (ChIP-seq), DNA affinity purification-sequencing (DAP-seq), and RNA-seq were used to uncover the gene expression response to MeGI. The results showed that the genes upregulated and downregulated in response to MeGI were mainly enriched in the circadian rhythm-related and flavonoid biosynthetic pathways, respectively. Additionally, the WRKY DNA-binding protein 28 (WRKY28) gene, which was detected by ChIP-seq, DAP-seq, and RNA-seq, was emphasized. WRKY28 has been reported to inhibit salicylic acid (SA) biosynthesis and was upregulated in MeGI-overexpressing A. thaliana flowers, suggesting that MeGI represses the SA level by increasing the expression level of WRKY28. This was confirmed that SA level was lower in D. oleifera female floral buds than male. Overall, our findings indicate that the MeGI mediates its sex control function in D. oleifera mainly by regulating genes in the circadian rhythm, SA biosynthetic, and flavonoid biosynthetic pathways.

Methodology for Constructing a Knowledgebase for Plant Gene Regulation Information

The amount of biological data is growing at a rapid pace as many high-throughput omics technologies and data pipelines are developed. This is resulting in the growth of databases for DNA and protein sequences, gene expression, protein accumulation, structural, and localization information. The diversity and multi-omics nature of such bioinformatic data requires well-designed databases for flexible organization and presentation. Besides general-purpose online bioinformatic databases, users need narrowly focused online databases to quickly access a meaningful collection of related data for their research. Here, we describe the methodology used to implement a plant gene regulatory knowledgebase, with data, query, and tool features, as well as the ability to expand to accommodate future datasets. We exemplify this methodology for the GRASSIUS knowledgebase, but it is applicable to developing and updating similar plant gene regulatory knowledgebases. GRASSIUS organizes and presents gene regulatory data from grass species with a central focus on maize (Zea mays). The main class of data presented include not only the families of transcription factors (TFs) and co-regulators (CRs) but also protein-DNA interaction data, where available.

TDTHub, a web server tool for the analysis of transcription factor binding sites in plants

Transcriptional regulation underlies most developmental programs and physiological responses to environmental changes in plants. Transcription factors (TFs) play a key role in the regulation of gene expression by binding specifically to short DNA sequences in the regulatory regions of genes: the TF binding sites (TFBSs). In recent years, several bioinformatic tools have been developed to detect TFBSs in candidate genes, either by de novo prediction or by directly mapping experimentally known TFBSs. However, most of these tools contain information for only a few species or require multi-step procedures, and are not always intuitive for non-experienced researchers. Here we present TFBS-Discovery Tool Hub (TDTHub), a web server for quick and intuitive studies of transcriptional regulation in plants. TDTHub uses pre-computed TFBSs in 40 plant species and allows the choice of two mapping algorithms, providing a higher versatility. Besides the main TFBS enrichment tool, TDTHub includes additional tools to assist in the analysis and visualization of data. In order to demonstrate the effectiveness of TDTHub, we analyzed the transcriptional regulation of the anthocyanin biosynthesis pathway. We also analyzed the transcriptional cascades in response to jasmonate and wounding in Arabidopsis and tomato (Solanum lycopersicum), respectively. In these studies, TDTHub helped to verify the most relevant TF nodes and to propose new ones with a prominent role in these pathways. TDTHub is available at http://acrab.cnb.csic.es/TDTHub/, and it will be periodically upgraded and expanded for new species and gene annotations.

De novo Transcriptome Analysis Revealed the Putative Pathway Genes Involved in Biosynthesis of Moracins in Morus alba L

Moracins, a kind of 2-phenyl-benzofuran compound from Moraceae, serve as phytoalexins with antimicrobial, anti-inflammatory, antitumor, and antidiabetes activities and respond to biotic and abiotic stresses, while their biosynthetic pathway and regulatory mechanism remain unclear. Here, we report a de novo transcriptome sequencing for different tissues of seedlings, as well as leaves under different stresses, in M. alba L. A total of 88 282 unigenes were assembled with an average length of 937 bp, and 82.2% of them were annotated. On the basis of the differential expression analysis and enzymatic activity assays in vitro, moracins were traced to the phenylpropanoid pathway, and a putative biosynthetic pathway of moracins was proposed. Unigenes coding key enzymes in the pathway were identified and their expression levels were verified by real-time quantitative reverse transcription PCR (qRT-PCR). Particularly, a p-coumaroyl CoA 2'-hydroxylase was presumed to be involved in the biosynthesis of stilbenes and deoxychalcones in mulberry. Additionally, the transcription factors that might participate in the regulation of moracin biosynthesis were obtained by coexpression analysis. These results shed light on the putative biosynthetic pathway of moracins, providing a basis for further investigation in functional characterization and transcriptional regulation of moracin biosynthesis in mulberry.