The developmental transcriptome atlas of the biofuel crop Camelina sativa

Nutri-cereal tissue-specific transcriptome atlas during development: Functional integration of gene expression to identify mineral uptake pathways in little millet (Panicum sumatrense)

Little millet (Panicum sumatrense Roth ex Roem. & Schult.) is an essential minor millet of southeast Asia and Africa's temperate and subtropical regions. The plant is stress-tolerant, has a short life cycle, and has a mineral-rich nutritional profile associated with unique health benefits. We report the developmental gene expression atlas of little millet (genotype JK-8) from ten tissues representing different stages of its life cycle, starting from seed germination and vegetative growth to panicle maturation. The developmental transcriptome atlas led to the identification of 342 827 transcripts. The BUSCO analysis and comparison with the transcriptomes of related species confirm that this study presents high-quality, in-depth coverage of the little millet transcriptome. In addition, the eFP browser generated here has a user-friendly interface, allowing interactive visualizations of tissue-specific gene expression. Using these data, we identified transcripts, the orthologs of which in Arabidopsis and rice are involved in nutrient acquisition, transport, and response pathways. The comparative analysis of the expression levels of these transcripts holds great potential for enhancing the mineral content in crops, particularly zinc and iron, to address the issue of "hidden hunger" and to attain nutritional security, making it a valuable asset for translational research.

The coordinated regulation of early meiotic stages is dominated by non-coding RNAs and stage-specific transcription in wheat

Reproductive success hinges on precisely coordinated meiosis, yet our understanding of how structural rearrangements of chromatin and phase transitions during meiosis are transcriptionally regulated is limited. In crop plants, detailed analysis of the meiotic transcriptome could identify regulatory genes and epigenetic regulators that can be targeted to increase recombination rates and broaden genetic variation, as well as provide a resource for comparison among eukaryotes of different taxa to answer outstanding questions about meiosis. We conducted a meiotic stage-specific analysis of messenger RNA (mRNA), small non-coding RNA (sncRNA), and long intervening/intergenic non-coding RNA (lincRNA) in wheat (Triticum aestivum L.) and revealed novel mechanisms of meiotic transcriptional regulation and meiosis-specific transcripts. Amidst general repression of mRNA expression, significant enrichment of ncRNAs was identified during prophase I relative to vegetative cells. The core meiotic transcriptome was comprised of 9309 meiosis-specific transcripts, 48 134 previously unannotated meiotic transcripts, and many known and novel ncRNAs differentially expressed at specific stages. The abundant meiotic sncRNAs controlled the reprogramming of central metabolic pathways by targeting genes involved in photosynthesis, glycolysis, hormone biosynthesis, and cellular homeostasis, and lincRNAs enhanced the expression of nearby genes. Alternative splicing was not evident in this polyploid species, but isoforms were switched at phase transitions. The novel, stage-specific regulatory controls uncovered here challenge the conventional understanding of this crucial biological process and provide a new resource of requisite knowledge for those aiming to directly modulate meiosis to improve crop plants. The wheat meiosis transcriptome dataset can be queried for genes of interest using an eFP browser located at https://bar.utoronto.ca/efp_wheat/cgi-bin/efpWeb.cgi?dataSource=Wheat_Meiosis.

Homozygosity mapping identified loci and candidate genes responsible for freezing tolerance in Camelina sativa

Homozygosity mapping is an effective tool for detecting genomic regions responsible for a given trait when the phenotype is controlled by a limited number of dominant or co-dominant loci. Freezing tolerance is a major attribute in agricultural crops such as camelina. Previous studies indicated that freezing tolerance differences between a tolerant (Joelle) and susceptible (CO46) variety of camelina were controlled by a small number of dominant or co-dominant genes. We performed whole genome homozygosity mapping to identify markers and candidate genes responsible for freezing tolerance difference between these two genotypes. A total of 28 F3 RILs were sequenced to ∼30× coverage, and parental lines were sequenced to >30-40× coverage with Pacific Biosciences high fidelity technology and 60× coverage using Illumina whole genome sequencing. Overall, about 126k homozygous single nucleotide polymorphism markers were identified that differentiate both parents. Moreover, 617 markers were also homozygous in F3 families fixed for freezing tolerance/susceptibility. All these markers mapped to two contigs forming a contiguous stretch of chromosome 11. The homozygosity mapping detected 9 homozygous blocks among the selected markers and 22 candidate genes with strong similarity to regions in or near the homozygous blocks. Two such genes were differentially expressed during cold acclimation in camelina. The largest block contained a cold-regulated plant thionin and a putative rotamase cyclophilin 2 gene previously associated with freezing resistance in arabidopsis (Arabidopsis thaliana). The second largest block contains several cysteine-rich RLK genes and a cold-regulated receptor serine/threonine kinase gene. We hypothesize that one or more of these genes may be primarily responsible for freezing tolerance differences in camelina varieties.

Sequencing of Camelina neglecta, a diploid progenitor of the hexaploid oilseed Camelina sativa

Camelina neglecta is a diploid species from the genus Camelina, which includes the versatile oilseed Camelina sativa. These species are closely related to Arabidopsis thaliana and the economically important Brassica crop species, making this genus a useful platform to dissect traits of agronomic importance while providing a tool to study the evolution of polyploids. A highly contiguous chromosome-level genome sequence of C. neglecta with an N50 size of 29.1 Mb was generated utilizing Pacific Biosciences (PacBio, Menlo Park, CA) long-read sequencing followed by chromosome conformation phasing. Comparison of the genome with that of C. sativa shows remarkable coincidence with subgenome 1 of the hexaploid, with only one major chromosomal rearrangement separating the two. Synonymous substitution rate analysis of the predicted 34 061 genes suggested subgenome 1 of C. sativa directly descended from C. neglecta around 1.2 mya. Higher functional divergence of genes in the hexaploid as evidenced by the greater number of unique orthogroups, and differential composition of resistant gene analogs, might suggest an immediate adaptation strategy after genome merger. The absence of genome bias in gene fractionation among the subgenomes of C. sativa in comparison with C. neglecta, and the complete lack of fractionation of meiosis-specific genes attests to the neopolyploid status of C. sativa. The assembled genome will provide a tool to further study genome evolution processes in the Camelina genus and potentially allow for the identification and exploitation of novel variation for Camelina crop improvement.

Toward kingdom-wide analyses of gene expression

Gene expression data for Archaeplastida are accumulating exponentially, with more than 300 000 RNA-sequencing (RNA-seq) experiments available for hundreds of species. The gene expression data stem from thousands of experiments that capture gene expression in various organs, tissues, cell types, (a)biotic perturbations, and genotypes. Advances in software tools make it possible to process all these data in a matter of weeks on modern office computers, giving us the possibility to study gene expression in a kingdom-wide manner for the first time. We discuss how the expression data can be accessed and processed and outline analyses that take advantage of cross-species analyses, allowing us to generate powerful and robust hypotheses about gene function and evolution.

Genome-wide identification and functional analysis of Dof transcription factor family in Camelina sativa

BACKGROUND

Dof transcription factors (TFs) containing C₂-C₂ zinc finger domains are plant-specific regulatory proteins, playing crucial roles in a variety of biological processes. However, little is known about Dof in Camelina sativa, an important oil crop worldwide, with high stress tolerance. In this study, a genome-wide characterization of Dof proteins is performed to examine their basic structural characteristics, phylogenetics, expression patterns, and functions to identify the regulatory mechanism underlying lipid/oil accumulation and the candidate Dofs mediating stress resistance regulation in C. sativa.

RESULTS

Total of 103 CsDof genes unevenly distributed on 20 chromosomes were identified from the C. sativa genome, and they were classified into four groups (A, B, C and D) based on the classification of Arabidopsis Dof gene family. All of the CsDof proteins contained the highly-conserved typic CX₂C-X₂₁-CX₂C structure. Segmental duplication and purifying selection were detected for CsDof genes. 61 CsDof genes were expressed in multiple tissues, and 20 of them showed tissue-specific expression patterns, suggesting that CsDof genes functioned differentially in different tissues of C. sativa. Remarkably, a set of CsDof members were detected to be possible involved in regulation of oil/lipid biosynthesis in C. sativa. Six CsDof genes exhibited significant expression changes in seedlings under salt stress treatment.

CONCLUSIONS

The present data reveals that segmental duplication is the key force responsible for the expansion of CsDof gene family, and a strong purifying pressure plays a crucial role in CsDofs' evolution. Several CsDof TFs may mediate lipid metabolism and stress responses in C. sativa. Several CsDof TFs may mediate lipid metabolism and stress responses in C. sativa. Collectively, our findings provide a foundation for deep understanding the roles of CsDofs and genetic improvements of oil yield and salt stress tolerance in this species and the related crops.

Genetic variation and structural diversity in major seed proteins among and within Camelina species

Genetic variation in seed protein composition, seed protein gene expression and predictions of seed protein physiochemical properties were documented in C. sativa and other Camelina species. Seed protein diversity was examined in six Camelina species (C. hispida, C. laxa, C. microcarpa, C. neglecta, C. rumelica and C. sativa). Differences were observed in seed protein electrophoretic profiles, total seed protein content and amino acid composition between the species. Genes encoding major seed proteins (cruciferins, napins, oleosins and vicilins) were catalogued for C. sativa and RNA-Seq analysis established the expression patterns of these and other genes in developing seed from anthesis through to maturation. Examination of 187 C. sativa accessions revealed limited variation in seed protein electrophoretic profiles, though sufficient to group the majority into classes based on high MW protein profiles corresponding to the cruciferin region. C. sativa possessed four distinct types of cruciferins, named CsCRA, CsCRB, CsCRC and CsCRD, which corresponded to orthologues in Arabidopsis thaliana with members of each type encoded by homeologous genes on the three C. sativa sub-genomes. Total protein content and amino acid composition varied only slightly; however, RNA-Seq analysis revealed that CsCRA and CsCRB genes contributed > 95% of the cruciferin transcripts in most lines, whereas CsCRC genes were the most highly expressed cruciferin genes in others, including the type cultivar DH55. This was confirmed by proteomics analyses. Cruciferin is the most abundant seed protein and contributes the most to functionality. Modelling of the C. sativa cruciferins indicated that each type possesses different physiochemical attributes that were predicted to impart unique functional properties. As such, opportunities exist to create C. sativa cultivars with seed protein profiles tailored to specific technical applications.

Spatio-temporal transcriptome profiling and subgenome analysis in Brassica napus

Brassica napus is an important oil crop and an allotetraploid species. However, the detailed analysis of gene function and homoeologous gene expression in all tissues at different developmental stages was not explored. In this study, we performed a global transcriptome analysis of 24 vegetative and reproductive tissues at six developmental stages (totally 111 tissues). These samples were clustered into eight groups. The gene functions of silique pericarp were similar to roots, stems and leaves. In particular, glucosinolate metabolic process was associated with root and silique pericarp. Genes involved in protein phosphorylation were often associated with stamen, anther and the early developmental stage of seeds. Transcription factor (TF) genes were more specific than structural genes. A total of 17 100 genes that were preferentially expressed in one tissue (tissue-preferred genes, TPGs), including 889 TFs (5.2%), were identified in the 24 tissues. Some TPGs were identified as hub genes in the co-expression network analysis, and some TPGs in different tissues were involved in different hormone pathways. About 67.0% of the homoeologs showed balanced expression, whereas biased expression of homoeologs was associated with structural divergence. In addition, the spatiotemporal expression of homoeologs was related to the presence of transposable elements (TEs) and regulatory elements (REs); more TEs and fewer REs in the promoters resulted in divergent expression in different tissues. This study provides a valuable transcriptional map for understanding the growth and development of B. napus, for identifying important genes for future crop improvement, and for exploring gene expression patterns in the B. napus.

Molecular and archaeological evidence on the geographical origin of domestication for Camelina sativa

PREMISE

Camelina (gold-of-pleasure or false flax) is an ancient oilseed crop with emerging applications in the production of sustainable, low-input biofuels. Previous domestication hypotheses suggested a European or western Asian origin, yet little genetic evidence has existed to assess the geographical origin for this crop, and archaeological data have not been systematically surveyed.

METHODS

We utilized genotyping-by-sequencing of 185 accessions of C. sativa and its wild relatives to examine population structure within the crop species and its relationship to populations of its wild progenitor, C. microcarpa; cytotype variation was also assessed in both species. In a complementary analysis, we surveyed the archaeological literature to identify sites with archaeobotanical camelina remains and assess the timing and prevalence of usage across Europe and western Asia.

RESULTS

The majority of C. microcarpa sampled in Europe and the United States belongs to a variant cytotype (2n = 38) with a distinct evolutionary origin from that of the crop lineage (2n = 40). Populations of C. microcarpa from Transcaucasia (South Caucasus) are most closely related to C. sativa based on cytotype and population structure; in combination with archaeological insights, these data refute prior hypotheses of a European domestication origin.

CONCLUSIONS

Our findings support a Caucasus, potentially Armenian, origin of C. sativa domestication. We cannot definitively determine whether C. sativa was intentionally targeted for domestication in its own right or instead arose secondarily through selection for agricultural traits in weedy C. sativa, as originally proposed by Vavilov for this species.

Engineering a feedback inhibition-insensitive plant dihydrodipicolinate synthase to increase lysine content in Camelina sativa seeds

Camelina sativa (camelina) is emerging as an alternative oilseed crop due to its short growing cycle, low input requirements, adaptability to less favorable growing environments and a seed oil profile suitable for biofuel and industrial applications. Camelina meal and oil are also registered for use in animal and fish feeds; however, like meals derived from most cereals and oilseeds, it is deficient in certain essential amino acids, such as lysine. In higher plants, the reaction catalyzed by dihydrodipicolinate synthase (DHDPS) is the first committed step in the biosynthesis of lysine and is subject to regulation by lysine through feedback inhibition. Here, we report enhancement of lysine content in C. sativa seed via expression of a feedback inhibition-insensitive form of DHDPS from Corynebacterium glutamicums (CgDHDPS). Two genes encoding C. sativa DHDPS were identified and the endogenous enzyme is partially insensitive to lysine inhibition. Site-directed mutagenesis was used to examine the impact of alterations, alone and in combination, present in lysine-desensitized DHDPS isoforms from Arabidopsis thaliana DHDPS (W53R), Nicotiana tabacum (N80I) and Zea mays (E84K) on C. sativa DHDPS lysine sensitivity. When introduced alone, each of the alterations decreased sensitivity to lysine; however, enzyme specific activity was also affected. There was evidence of molecular or structural interplay between residues within the C. sativa DHDPS allosteric site as coupling of the W53R mutation with the N80V mutation decreased lysine sensitivity of the latter, but not to the level with the W53R mutation alone. Furthermore, the activity and lysine sensitivity of the triple mutant (W53R/N80V/E84T) was similar to the W53R mutation alone or the C. glutamicum DHDPS. The most active and most lysine-insensitive C. sativa DHDPS variant (W53R) was not inhibited by free lysine up to 1 mM, comparable to the C. glutamicums enzyme. Seed lysine content increased 13.6 -22.6% in CgDHDPS transgenic lines and 7.6–13.2% in the mCsDHDPS lines. The high lysine-accumulating lines from this work may be used to produce superior quality animal feed with improved essential amino acid profile.

Recent Advances in Genome-wide Analyses of Plant Potassium Transporter Families

Plants require potassium (K⁺) as a macronutrient to support numerous physiological processes. Understanding how this nutrient is transported, stored, and utilized within plants is crucial for breeding crops with high K⁺ use efficiency. As K⁺ is not metabolized, cross-membrane transport becomes a rate-limiting step for efficient distribution and utilization in plants. Several K⁺ transporter families, such as KUP/HAK/KT and KEA transporters and Shaker-like and TPK channels, play dominant roles in plant K⁺ transport processes. In this review, we provide a comprehensive contemporary overview of our knowledge about these K⁺ transporter families in angiosperms, with a major focus on the genome-wide identification of K⁺ transporter families, subcellular localization, spatial expression, function and regulation. We also expanded the genome-wide search for the K⁺ transporter genes and examined their tissue-specific expression in Camelina sativa, a polyploid oil-seed crop with a potential to adapt to marginal lands for biofuel purposes and contribution to sustainable agriculture. In addition, we present new insights and emphasis on the study of K⁺ transporters in polyploids in an effort to generate crops with high K⁺ Utilization Efficiency (KUE).

CRISPR- mediated Mutation in Cinnamoyl- CoA Reductase 4 in Allohexaploid Oilseed Crop Camelina sativa, Revealed its Pivotal Role in Resistance Against Sclerotinia sclerotiorum

BACKGROUND

Sclerotinia sclerotiorum (Ss) is a broad host range necrotrophic ascomycete fungus affecting over 400 plant species. Ss causes stem rot disease on Camelina sativa (Cs) an allohexaploid crucifer species that is promoted as a low input crop and industrial oil attributes suitable as biofuel and lubricant. Histochemical and molecular studies has linked resistance to Ss in C. sativa with the cell wall lignification (Eynck et al., 2012) and reported constitutive expression of Cinnamoyl-CoA Reductase 4 (CsCCR4) gene, in the Cs resistant line CN114263. Modern breeding efforts, such as gene editing, are needed to improve commercial lines and to limit the risk of crop loss which would be substantial to producers.

OBJECTIVES

To investigate the importance of monolignol biosynthesis and the role of CsCCR4 in Camelina resistance to Ss we generated CsCCR4 knockout mutants of CN114263 Camelina line using CRISPR/Cas9-mediated gene editing.

MATERIALS AND METHODS

Thirty T1 plants were produced via floral dip transformation followed by glyphosate spraying that was used in the first step of screening procedures and were confirmed by PCR method. Transgene's T-DNA copy number variation, T-DNA CNV, in T1 and T2 progenitors were determined using digital droplet PCR (ddPCR) and the occurrence of mutation in the three copies of CsCCR4 homeologues in T1 and T2 generations were scrutinized by drop-off assay technique. To make sure that if the created mutants in T2 plants are real, TOPO TA sequencing flanking the Cas9/gRNA specific hot point of cleavage for three of them was conducted.

RESULTS

In the T1 generation, 25 plants were confirmed which had between one to nine T-DNA copies in the corresponding Camelina genome. In T2 generation the population were screened for potential mutation in CsCCR4 gene. Various types of mutations, including insertions and deletions, were demonstrated in three copies of CsCCR4. In fact, CRISPR system could have cut one, two or three copies of the gene in events numbered T2-plant 10, T2-plant 15 and T2-plant 19, respectively. The T3-plant 19 which showed mutation in all versions of CsCCR4 in previous generation had susceptibility to S. sclerotiorum invasion and was kept as real csccr4 mutant material for further investigations of Camelina-Sclerotinia interaction. Mutation in CsCCR4 had occurred through error-prone none- homologous end joining (NHEJ) nucleus DNA repair pathway. Ss challenge on the early flowering T3 generation. The T3 plants with mutation causing premature stop codon at position 217 of CsCCR4 were compromised in their resistance to Ss compared to the wildtype resistant control parent CN114263.

CONCLUSION

Using ddPCR it easily was possible to identify both the T-DNA CNV and occurrence of mutation in CsCCR4 homeologues in T1 and T2 progenitors. We illustrated that CRISPR/Cas9-mediated mutation is a decent technique that can be utilized to expedite the mutant line development which could assist to figure out the activity of a CsCCR4 gene in defense responses to the pathogens in C. sativa as prospective oilseed crop for biodiesel production.

Myxospermy Evolution in Brassicaceae: A Highly Complex and Diverse Trait with Arabidopsis as an Uncommon Model

The ability to extrude mucilage upon seed imbibition (myxospermy) occurs in several Angiosperm taxonomic groups, but its ancestral nature or evolutionary convergence origin remains misunderstood. We investigated seed mucilage evolution in the Brassicaceae family with comparison to the knowledge accumulated in Arabidopsis thaliana. The myxospermy occurrence was evaluated in 27 Brassicaceae species. Phenotyping included mucilage secretory cell morphology and topochemistry to highlight subtle myxospermy traits. In parallel, computational biology was driven on the one hundred genes constituting the so-called A. thaliana mucilage secretory cell toolbox to confront their sequence conservation to the observed phenotypes. Mucilage secretory cells show high morphology diversity; the three studied Arabidopsis species had a specific extrusion modality compared to the other studied Brassicaceae species. Orthologous genes from the A. thaliana mucilage secretory cell toolbox were mostly found in all studied species without correlation with the occurrence of myxospermy or even more sub-cellular traits. Seed mucilage may be an ancestral feature of the Brassicaceae family. It consists of highly diverse subtle traits, probably underlined by several genes not yet characterized in A. thaliana or by species-specific genes. Therefore, A. thaliana is probably not a sufficient reference for future myxospermy evo-devo studies.

Chia (Salvia hispanica) Gene Expression Atlas Elucidates Dynamic Spatio-Temporal Changes Associated With Plant Growth and Development

Chia (Salvia hispanica L.), now a popular superfood and a pseudocereal, is one of the richest sources of dietary nutrients such as protein, fiber, and polyunsaturated fatty acids (PUFAs). At present, the genomic and genetic information available in the public domain for this crop are scanty, which hinders an understanding of its growth and development and genetic improvement. We report an RNA-sequencing (RNA-Seq)-based comprehensive transcriptome atlas of Chia sampled from 13 tissue types covering vegetative and reproductive growth stages. We used ~355 million high-quality reads of total ~394 million raw reads from transcriptome sequencing to generate de novo reference transcriptome assembly and the tissue-specific transcript assemblies. After the quality assessment of the merged assemblies and implementing redundancy reduction methods, 82,663 reference transcripts were identified. About 65,587 of 82,663 transcripts were translated into 99,307 peptides, and we were successful in assigning InterPro annotations to 45,209 peptides and gene ontology (GO) terms to 32,638 peptides. The assembled transcriptome is estimated to have the complete sequence information for ~86% of the genes found in the Chia genome. Furthermore, the analysis of 53,200 differentially expressed transcripts (DETs) revealed their distinct expression patterns in Chia's vegetative and reproductive tissues; tissue-specific networks and developmental stage-specific networks of transcription factors (TFs); and the regulation of the expression of enzyme-coding genes associated with important metabolic pathways. In addition, we identified 2,411 simple sequence repeats (SSRs) as potential genetic markers from the transcripts. Overall, this study provides a comprehensive transcriptome atlas, and SSRs, contributing to building essential genomic resources to support basic research, genome annotation, functional genomics, and molecular breeding of Chia.

Functional characterization of an novel acyl-CoA:diacylglycerol acyltransferase 3-3 (CsDGAT3-3) gene from Camelina sativa

Diacylglycerol acyltransferases (DGAT) catalyze the final committed step of de novo biosynthesis of triacylglycerol (TAG) in plant seeds. This study was to functionally characterize DGAT3 genes in Camelina sativa, an important oil crops accumulating high levels of unsaturated fatty acids (UFAs) in seeds. Three camelina DGAT3 genes (CsDGAT3-1, CsDGAT3-2 and CsDGAT3-3) were identified, and the encoded proteins were predicted to be cytosolic-soluble proteins present as a homodimer containing the 2Fe-2S domain. They had divergent expression patterns in various tissues, suggesting that they may function in tissue-specific manner with CsDGAT3-1 in roots, CsDGAT3-2 in flowers and young seedlings, and CsDGAT3-3 in developing seeds. Functional complementation assay in yeast demonstrated that CsDGAT3-3 restored TAG synthesis. TAG content and UFAs, particularly eicosenoic acid (EA, 20:1n-9) were largely increased by adding exogenous UFAs in the yeast medium. Further heterogeneously transient expression in N. benthamiana leaves and seed-specific expression in tobacco seeds indicated that CsDGAT3-3 significantly enhanced oil and UFA accumulation with much higher level of EA. Overall, CsDGAT3-3 exhibited a strong abilty catalyzing TAG synthesis and high substrate preference for UFAs, especially for 20:1n-9. The present data provide new insights for further understanding oil biosynthesis mechanism in camelina seeds, indicating that CsDGAT3-3 may have practical applications for increasing both oil yield and quality.

CamRegBase: a gene regulation database for the biofuel crop, Camelina sativa

Camelina is an annual oilseed plant from the Brassicaceae family that is gaining momentum as a biofuel winter cover crop. However, a significant limitation in further enhancing its utility as a producer of oils that can be used as biofuels, jet fuels or bio-based products is the absence of a repository for all the gene expression and regulatory information that is being rapidly generated by the community. Here, we provide CamRegBase (https://camregbase.org/) as a one-stop resource to access Camelina information on gene expression and co-expression, transcription factors, lipid associated genes and genome-wide orthologs in the close-relative reference plant Arabidopsis. We envision this as a resource of curated information for users, as well as a repository of new gene regulation information.

Identification of the rhizospheric microbe and metabolites that led by the continuous cropping of ramie (Boehmeria nivea L. Gaud)

Continuous cropping lowers the production and quality of ramie (Boehmeria nivea L. Gaud). This study aimed to reveal the metagenomic and metabolomic changes between the healthy- and obstacle-plant after a long period of continuous cropping. After 10 years of continuous cropping, ramie planted in some portions of the land exhibited weak growth and low yield (Obstacle-group), whereas, ramie planted in the other portion of the land grew healthy (Health-group). We collected rhizosphere soil and root samples from which measurements of soil chemical and plant physiochemical properties were taken. All samples were subjected to non-targeted gas chromatograph-mass spectrometer (GS/MS) metabolome analysis. Further, metagenomics was performed to analyze the functional genes in rhizospheric soil organisms. Based on the findings, ramie in Obstacle-group were characterized by shorter plant height, smaller stem diameter, and lower fiber production than that in Health-group. Besides, the Obstacle-group showed a lower relative abundance of Rhizobiaceae, Lysobacter antibioticus, and Bradyrhizobium japonicum, but a higher relative abundance of Azospirillum lipoferum and A. brasilense compared to the Health-group. Metabolomic analysis results implicated cysteinylglycine (Cys-Gly), uracil, malonate, and glycerol as the key differential metabolites between the Health- and Obstacle-group. Notably, this work revealed that bacteria such as Rhizobia potentially synthesize IAA and are likely to reduce the biotic stress of ramie. L. antibioticus also exerts a positive effect on plants in the fight against biotic stress and is mediated by metabolites including orthophosphate, uracil, and Cys-Gly, which may serve as markers for disease risk. These bacterial effects can play a key role in plant resistance to biotic stress via metabolic and methionine metabolism pathways.

Genome-wide identification and functional characterization of the Camelina sativa WRKY gene family in response to abiotic stress

Background

WRKY transcription factors are a superfamily of regulators involved in diverse biological processes and stress responses in plants. However, there is limited knowledge about the WRKY family in camelina (Camelina sativa), an important Brassicaceae oil crop with strong tolerance for various stresses. Here, a genome-wide characterization of WRKY proteins is performed to examine their gene structures, phylogenetics, expression, conserved motif organizations, and functional annotation to identify candidate WRKYs that mediate stress resistance regulation in camelinas.

Results

A total of 242 CsWRKY proteins encoded by 224 gene loci distributed unevenly over the chromosomes were identified, and they were classified into three groups by phylogenetic analysis according to their WRKY domains and zinc finger motifs. The 15 CsWRKY gene loci generated 33 spliced variants. Orthologous WRKY gene pairs were identified, with 173 pairs in the C. sativa and Arabidopsis genomes as well as 282 pairs in the C. sativa and B. napus genomes, respectively. A total of 137 segmental duplication events were observed, but there was no tandem duplication in the camelina genome. Ten major conserved motifs were examined, with WRKYGQK being the most conserved, and several variants were present in many CsWRKYs. Expression analysis revealed that 50% more CsWRKY genes were expressed constitutively, and a set of them displayed tissue-specific expression. Notably, 11 CsWRKY genes exhibited significant expression changes in seedlings under cold, salt, and drought stresses, showing a preferentially inducible expression pattern in response to the stress.

Conclusions

The present article describes a detailed analysis of the CsWRKY gene family and its expression profiles in 12 tissues and under several stress conditions. Segmental duplication is the major force underlying the broad expansion of this gene family, and a strong purifying pressure occurred for CsWRKY proteins during their evolution. CsWRKY proteins play important roles in plant development, with differential functions in different tissues. Exceptionally, eleven CsWRKYs, particularly five alternative spliced isoforms, were found to be the possible key players in mediating plant responses to various stresses. Overall, our results provide a foundation for understanding the roles of CsWRKYs and the precise mechanism through which CsWRKYs regulate high stress resistance as well as the development of stress tolerance cultivars among Cruciferae crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07189-3.

Genome-Wide Analysis of the Five Phosphate Transporter Families in Camelina sativa and Their Expressions in Response to Low-P

Phosphate transporters (PHTs) play pivotal roles in phosphate (Pi) acquisition from the soil and distribution throughout a plant. However, there is no comprehensive genomic analysis of the PHT families in Camelina sativa, an emerging oilseed crop. In this study, we identified 73 CsPHT members belonging to the five major PHT families. A whole-genome triplication event was the major driving force for CsPHT expansion, with three homoeologs for each Arabidopsis ortholog. In addition, tandem gene duplications on chromosome 11, 18 and 20 further enlarged the CsPHT1 family beyond the ploidy norm. Phylogenetic analysis showed clustering of the CsPHT1 and CsPHT4 family members into four distinct groups, while CsPHT3s and CsPHT5s were clustered into two distinct groups. Promoter analysis revealed widespread cis-elements for low-P response (P1BS) specifically in CsPHT1s, consistent with their function in Pi acquisition and translocation. In silico RNA-seq analysis revealed more ubiquitous expression of several CsPHT1 genes in various tissues, whereas CsPHT2s and CsPHT4s displayed preferential expression in leaves. While several CsPHT3s were expressed in germinating seeds, most CsPHT5s were expressed in floral and seed organs. Suneson, a popular Camelina variety, displayed better tolerance to low-P than another variety, CS-CROO, which could be attributed to the higher expression of several CsPHT1/3/4/5 family genes in shoots and roots. This study represents the first effort in characterizing CsPHT transporters in Camelina, a promising polyploid oilseed crop that is highly tolerant to abiotic stress and low-nutrient status, and may populate marginal soils for biofuel production.

Novel genes associated with folic acid-mediated metabolism in mouse: A bioinformatics study

Folic acid plays an essential role in the central nervous system and cancer. This study aimed to screen genes related to folic acid metabolism. Datasets (GSE80587, GSE65267 and GSE116299) correlated to folic acid were screened in the Gene Expression Omnibus. Weighed gene co-expression network analysis was performed to identify modules associated with sample traits of folic acid and organs (brain, prostate and kidney). Functional enrichment analysis was performed for the eigengenes in modules that were significantly correlated with sample traits. Accordingly, the hub genes and key nodes in the modules were identified using the protein interaction network. A total of 17,252 genes in three datasets were identified. One module, which included 97 genes that were highly correlated with sample traits (including folic acid treatment [cor = -0.57, P = 3e-04] and kidney [cor = -0.68, p = 4e-06]), was screened out. Hub genes, including tetratricopeptide repeat protein 38 (Ttc38) and miR-185, as well as those (including Sema3A, Insl3, Dll1, Msh4 and Snai1) associated with "neuropilin binding", "regulation of reproductive process" and "vitamin D metabolic process", were identified. Genes, including Ttc38, Sema3A, Insl3, Dll1, Msh4 and Snai1, were the novel factors that may be associated with the development of the kidneys and related to folic acid treatment.

Assessing Diversity in the Camelina Genus Provides Insights into the Genome Structure of Camelina sativa

Camelina sativa (L.) Crantz an oilseed crop of the Brassicaceae family is gaining attention due to its potential as a source of high value oil for food, feed or fuel. The hexaploid domesticated C. sativa has limited genetic diversity, encouraging the exploration of related species for novel allelic variation for traits of interest. The current study utilized genotyping by sequencing to characterize 193 Camelina accessions belonging to seven different species collected primarily from the Ukrainian-Russian region and Eastern Europe. Population analyses among Camelina accessions with a 2n = 40 karyotype identified three subpopulations, two composed of domesticated C. sativa and one of C. microcarpa species. Winter type Camelina lines were identified as admixtures of C. sativa and C. microcarpa. Eighteen genotypes of related C. microcarpa unexpectedly shared only two subgenomes with C. sativa, suggesting a novel or cryptic sub-species of C. microcarpa with 19 haploid chromosomes. One C. microcarpa accession (2n = 26) was found to comprise the first two subgenomes of C. sativa suggesting a tetraploid structure. The defined chromosome series among C. microcarpa germplasm, including the newly designated C. neglecta diploid née C. microcarpa, suggested an evolutionary trajectory for the formation of the C. sativa hexaploid genome and re-defined the underlying subgenome structure of the reference genome.

Tung Tree (Vernicia fordii) Genome Provides A Resource for Understanding Genome Evolution and Improved Oil Production

Tung tree (Vernicia fordii) is an economically important woody oil plant that produces tung oil rich in eleostearic acid. Here, we report a high-quality chromosome-scale genome sequence of tung tree. The genome sequence was assembled by combining Illumina short reads, Pacific Biosciences single-molecule real-time long reads, and Hi-C sequencing data. The size of tung tree genome is 1.12 Gb, with 28,422 predicted genes and over 73% repeat sequences. The V. fordii underwent an ancient genome triplication event shared by core eudicots but no further whole-genome duplication in the subsequent ca. 34.55 million years of evolutionary history of the tung tree lineage. Insertion time analysis revealed that repeat-driven genome expansion might have arisen as a result of long-standing long terminal repeat retrotransposon bursts and lack of efficient DNA deletion mechanisms. The genome harbors 88 resistance genes encoding nucleotide-binding sites; 17 of these genes may be involved in early-infection stage of Fusarium wilt resistance. Further, 651 oil-related genes were identified, 88 of which are predicted to be directly involved in tung oil biosynthesis. Relatively few phosphoenolpyruvate carboxykinase genes, and synergistic effects between transcription factors and oil biosynthesis-related genes might contribute to the high oil content of tung seed. The tung tree genome constitutes a valuable resource for understanding genome evolution, as well as for molecular breeding and genetic improvements for oil production.

Interspecific Variation in the Unsaturation Level of Seed Oils Were Associated With the Expression Pattern Shifts of Duplicated Desaturase Genes and the Potential Role of Other Regulatory Genes

Seed oils are of great economic importance both for human consumption and industrial applications. The nutritional quality and industrial value of seed oils are mostly determined by their fatty acid profiles, especially the relative proportions of unsaturated fatty acids. Tree peony seed oils have recently been recognized as novel edible oils enriched in α-linolenic acid (ALA). However, congeneric species, such as Paeonia ostii and P. ludlowii, showed marked variation in the relative proportions of different unsaturated fatty acids. By comparing the dynamics of fatty acid accumulation and the time-course gene expression patterns between P. ostii and P. ludlowii, we identified genes that were differentially expressed between two species in developing seeds, and showed congruent patterns of variation between expression levels and phenotypes. In addition to the well-known desaturase and acyltransferase genes associated with fatty acid desaturation, among them were some genes that were conservatively co-expressed with the desaturation pathway genes across phylogenetically distant ALA-rich species, including Camelina sativa and Perilla frutescens. Go enrichment analysis revealed that these genes were mainly involved in transcriptional regulation, protein post-translational modification and hormone biosynthesis and response, suggesting that the fatty acid synthesis and desaturation pathway might be subject to multiple levels of regulation.

Mapping quantitative trait loci for seed traits in Camelina sativa

Genetic dissection of oil content and seed size in Camelina sativa was conducted by QTL mapping using a SNP-based linkage map and a recombinant inbred population. Camelina (Camelina sativa L. Crantz) is an oilseed crop that has great potential to provide sustainable feedstock for biofuel production and to improve dryland agriculture. A major breeding objective for camelina is to increase seed size and oil content. Understanding the genetics behind variations of seed size and associated traits such as oil content would help breeders develop varieties of increased oil yield that are more robust, easier to plant and harvest, and better for oil processing. In this study, we developed a recombinant inbred population derived from the two camelina accessions, Suneson and Pryzeth, with contrasting traits, especially seed size and oil content. Using 189 lines, a genetic map was constructed containing 2376 single nucleotide polymorphism markers spanning 2034.6 cM of 20 linkage groups with an average density of 1.5 cM per locus. Field trials were conducted for 2 years (2017 and 2018) in two environments (dryland and irrigated) in Bozeman, Montana. The results revealed important correlations of seed size with other associated traits such as oil content, pod size and seed number per pod. Significant QTLs were also discovered for these traits. The results of this study are the first step to isolate genes controlling seed development and oil accumulation and to develop advanced varieties of camelina better adapted to modern agriculture by marker-assisted breeding.

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

CRISPR/Cas9 editing of three CRUCIFERIN C homoeologues alters the seed protein profile in Camelina sativa

BACKGROUND

The oilseed Camelina sativa is grown for a range of applications, including for biofuel, biolubricants, and as a source of omega-3 fatty acids for the aquaculture feed industry. The seed meal co-product is used as a source of protein for animal feed; however, the low value of the meal hinders profitability and more widespread application of camelina. The nutritional quality of the seed meal is largely determined by the abundance of specific seed storage proteins and their amino acid composition. Manipulation of seed storage proteins has been shown to be an effective means for either adjustment of nutritional content of seeds or for enhancing accumulation of high-value recombinant proteins in seeds.

RESULTS

CRISPR/Cas9 gene editing technology was used to generate deletions in the first exon of the three homoeologous genes encoding the seed storage protein CRUCIFERIN C (CsCRUC), creating an identical premature stop-codon in each and resulting in a CsCRUC knockout line. The mutant alleles were detected by applying a droplet digital PCR drop-off assay. The quantitative nature of this technique is particularly valuable when applied to polyploid species because it can accurately determine the number of mutated alleles in a gene family. Loss of CRUC protein did not alter total seed protein content; however, the abundance of other cruciferin isoforms and other seed storage proteins was altered. Consequently, seed amino acid content was significantly changed with an increase in the proportion of alanine, cysteine and proline, and decrease of isoleucine, tyrosine and valine. CsCRUC knockout seeds did not have changed total oil content, but the fatty acid profile was significantly altered with increased relative abundance of all saturated fatty acids.

CONCLUSIONS

This study demonstrates the plasticity of the camelina seed proteome and establishes a CRUC-devoid line, providing a framework for modifying camelina seed protein composition. The results also illustrate a possible link between the composition of the seed proteome and fatty acid profile.

Genes associated with chloroplasts and hormone-signaling, and transcription factors other than CBFs are associated with differential survival after low temperature treatments of Camelina sativa biotypes

Winter annual biotypes of Camelina sativa regularly survive after winter conditions experienced in northern regions of the U.S., whereas summer annual biotypes do not. To determine potential molecular mechanisms associated with these biotype differences in survival after low temperature treatments, we examined genetic and transcript variations in both a winter- (Joelle) and a summer- (CO46) biotype. It was determined that as few as one or two dominant genes may control differential survival after low temperature treatments. Of the 1797 genes that were differentially expressed in response to cold in both the winter and summer biotypes many COR genes were identified, indicating that the CBF regulon is functional in both. However, only 153 and 76 genes from Joelle and CO46, respectively, were either differentially expressed or not expressed at all in one biotype versus the other following cold acclimation. We hypothesize that these 229 genes play a significant role in, or are primarily responsive to, differences in survival after freezing between these two biotypes. Promoter analysis provided few clues as to the regulation or these genes; however, genes that were down-regulated specifically in the winter biotype Joelle were enriched with the sequence TGGCCCTCGCTCAC, which is over-represented among genes associated with chloroplasts in Arabidopsis. Additionally, several genes involved in auxin signaling were down-regulated specifically in Joelle. A transcription factor with strong similarity to MYB47, known to be up-regulated by salt, drought, and jasmonic acid, but not cold in Arabidopsis, was essentially off in the freezing sensitive biotype CO46, but was cold-induced in the winter biotype Joelle. Several other transcription factors genes including three with similarity to WRKY70, that may be involved in SA/JA-dependent responses, a HOMEOBOX 6 gene involved in ABA signaling, and two others (NUCLEAR FACTOR Y and CONSTANS-like 2) known to be implicated in photoperiodic flowering were also differentially expressed between the two biotypes.

Mutagenesis of FAD2 genes in peanut with CRISPR/Cas9 based gene editing

Background

Increasing the content of oleic acid in peanut seeds is one of the major goals in peanut breeding due to consumer and industry benefits, such as anti-oxidation and long shelf-life. Homeologous ahFAD2A and ahFAD2B genes encode fatty acid desaturases, which are the key enzymes for converting oleic acid to linoleic acid that oxidizes readily. To date, all high oleic acid peanut varieties result from natural mutations occurred in both genes. A method to induce mutations in the genes of other elite cultivars could speed introgression of this valuable trait. The gene-editing approach utilizing CRISPR/Cas9 technology was employed to induce de novo mutations in the ahFAD2 genes using peanut protoplasts and hairy root cultures as models.

Results

The hot spot of natural mutation in these genes was selected as the target region. Appropriate sgRNAs were designed and cloned into a CRISPR/Cas9 expression plasmid. As a result of CRISPR/Cas9 activity, three mutations were identified - G448A in ahFAD2A, and 441_442insA and G451T in ahFAD2B. The G448A and 441_442insA mutations are the same as those seen in existing high oleate varieties and the G451T is new mutation. Because natural mutations appear more often in the ahFAD2A gene than in the ahFAD2B gene in subspecies A. hypogaea var. hypogaea, the mutations induced in ahFAD2B by gene editing may be useful in developing high oleate lines with many genetic backgrounds after validation of oleic acid content in the transformed lines. The appearance of the G448A mutation in ahFAD2A is a further benefit for high oleic acid oil content.

Conclusions

Overall, these results showed that mutations were, for the first time, induced by CRISPR-based gene editing approach in peanut. This research demonstrated the potential application of gene editing for mutagenesis in peanut and suggested that CRISPR/Cas9 technology may be useful in the peanut breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12896-019-0516-8) contains supplementary material, which is available to authorized users.

Assessment of Antinutritional Compounds and Chemotaxonomic Relationships between Camelina sativa and Its Wild Relatives

We compared the secondary metabolite composition in seeds of Camelina sativa and its wild relatives to identify potential germplasm with reduced levels of antinutritional compounds. Twenty Camelina accessions, from five different species, were analyzed by liquid chromatography mass spectrometry and subjected to principal component analysis, which revealed that Camelina spp. separated into distinct chemotaxonomic groups. Three major glucosinolates (GSs) were identified in our study, namely, 9-methylsulfinylnonyl GS (GS9), 10-methylsulfinyldecyl GS (GS10), and 11-methylsulfinylundecyl GS (GS11). While there were differences in total GS levels, species-specific patterns for GS9 and GS11 were noted. Sinapine content ranged between 1.4 and 5.6 mg/g FW, with the lowest levels observed in  C. laxa and C. sativa. Lignin levels were also lowest in C. sativa, with most accessions containing less than 6 mg/g FW. Our results show that wild Camelina spp. have distinct metabolomes, and based on their levels of major antinutritionals, some could be incorporated into breeding programs with C. sativa.

Unbiased subgenome evolution following a recent whole-genome duplication in pear (Pyrus bretschneideri Rehd.)

Genome fractionation (also known as diploidization) frequently occurs following paleopolyploidization events. Biased fractionation between subgenomes has been found in some paleo-allopolyploids, while this phenomenon is absent in paleo-autopolyploids. Pear (Pyrus bretschneideri Rehd.) experienced a recent whole-genome duplication (WGD, ~30 million years ago); however, the evolutionary fate of the two subgenomes derived from this WGD event is not clear. In this study, we identified the two paleo-subgenomes in pear using peach (Prunus persica) as an outgroup and investigated differences in the gene loss rate, evolutionary rate, gene expression level, and DNA methylation level between these two subgenomes. Fractionation bias was not found between the two pear subgenomes, which evolved at similar evolutionary rates. The DNA methylation level of the two subgenomes showed little bias, and we found no expression dominance between the subgenomes. However, we found that singleton genes and homeologous genes within each subgenome showed divergent evolutionary patterns of selective constraints, expression and epigenetic modification. These results provide insights into subgenome evolution following paleopolyploidization in pear.

Camelina sativa, an oilseed at the nexus between model system and commercial crop

The rapid assessment of metabolic engineering strategies in plants is aided by crops that provide simple, high throughput transformation systems, a sequenced genome, and the ability to evaluate the resulting plants in field trials. Camelina sativa provides all of these attributes in a robust oilseed platform. The ability to perform field evaluation of Camelina is a useful, and in some studies essential benefit that allows researchers to evaluate how traits perform outside the strictly controlled conditions of a greenhouse. In the field the plants are subjected to higher light intensities, seasonal diurnal variations in temperature and light, competition for nutrients, and watering regimes dictated by natural weather patterns, all which may affect trait performance. There are difficulties associated with the use of Camelina. The current genetic resources available for Camelina pale in comparison to those developed for the model plant Arabidopsis thaliana; however, the sequence similarity of the Arabidopsis and Camelina genomes often allows the use of Arabidopsis as a reference when additional information is needed. Camelina's genome, an allohexaploid, is more complex than other model crops, but the diploid inheritance of its three subgenomes is straightforward. The need to navigate three copies of each gene in genome editing or mutagenesis experiments adds some complexity but also provides advantages for gene dosage experiments. The ability to quickly engineer Camelina with novel traits, advance generations, and bulk up homozygous lines for small-scale field tests in less than a year, in our opinion, far outweighs the complexities associated with the crop.

AgriSeqDB: an online RNA-Seq database for functional studies of agriculturally relevant plant species

BACKGROUND

The genome-wide expression profile of genes in different tissues/cell types and developmental stages is a vital component of many functional genomic studies. Transcriptome data obtained by RNA-sequencing (RNA-Seq) is often deposited in public databases that are made available via data portals. Data visualization is one of the first steps in assessment and hypothesis generation. However, these databases do not typically include visualization tools and establishing one is not trivial for users who are not computational experts. This, as well as the various formats in which data is commonly deposited, makes the processes of data access, sharing and utility more difficult. Our goal was to provide a simple and user-friendly repository that meets these needs for data-sets from major agricultural crops.

DESCRIPTION

AgriSeqDB ( https://expression.latrobe.edu.au/agriseqdb ) is a database for viewing, analysing and interpreting developmental and tissue/cell-specific transcriptome data from several species, including major agricultural crops such as wheat, rice, maize, barley and tomato. The disparate manner in which public transcriptome data is often warehoused and the challenge of visualizing raw data are both major hurdles to data reuse. The popular eFP browser does an excellent job of presenting transcriptome data in an easily interpretable view, but previous implementation has been mostly on a case-by-case basis. Here we present an integrated visualisation database of transcriptome data-sets from six species that did not previously have public-facing visualisations. We combine the eFP browser, for gene-by-gene investigation, with the Degust browser, which enables visualisation of all transcripts across multiple samples. The two visualisation interfaces launch from the same point, enabling users to easily switch between analysis modes. The tools allow users, even those without bioinformatics expertise, to mine into data-sets and understand the behaviour of transcripts of interest across samples and time. We have also incorporated an additional graphic download option to simplify incorporation into presentations or publications.

CONCLUSION

Powered by eFP and Degust browsers, AgriSeqDB is a quick and easy-to-use platform for data analysis and visualization in five crops and Arabidopsis. Furthermore, it provides a tool that makes it easy for researchers to share their data-sets, promoting research collaborations and data-set reuse.

Expression of FLOWERING LOCUS C and a frameshift mutation of this gene on chromosome 20 differentiate a summer and winter annual biotype of Camelina sativa

The nature of the vegetative to reproductive transition in the shoot apical meristem of Camelina sativa summer annual cultivar CO46 and winter annual cultivar Joelle was confirmed by treating seedlings with or without 8 weeks of vernalization. True to their life cycle classification, Joelle required a vernalization treatment to induce bolting and flowering, whereas CO46 did not. In this study, whole genome sequence, RNAseq, and resequencing of PCR-amplified transcripts for a key floral repressor were used to better understand factors involved in the flowering habit of summer and winter biotypes at the molecular level. Analysis of transcriptome data indicated that abundance for one of the three genes encoding the floral repressor FLOWERING LOCUS C (FLC; Csa20 g015400) was 16-fold greater in Joelle compared to CO46 prior to vernalization. Abundance of this transcript decreased only slightly in CO46 postvernalization, compared to a substantial decrease in Joelle. The results observed in the winter annual biotype Joelle are consistent with repression of FLC by vernalization. Further characterization of FLC at both the genome and transcriptome levels identified a one base deletion in the 5th exon coding for a keratin-binding domain in chromosome 20 of CO46 and Joelle. The one base deletion detected in chromosome 20 FLC is predicted to result in a frameshift that would produce a nonfunctional protein. Analysis of whole genome sequence indicated that the one base deletion in chromosome 20 FLC occurred at a greater ratio in the summer biotype CO46 (2:1) compared to the winter biotype Joelle (1:4); similar trends were also observed for RNAseq and cDNA transcripts mapping to chromosome 20 FLC of CO46 and Joelle.

Different Modes of Gene Duplication Show Divergent Evolutionary Patterns and Contribute Differently to the Expansion of Gene Families Involved in Important Fruit Traits in Pear (Pyrus bretschneideri)

Pear is an important fruit crop of the Rosaceae family and has experienced two rounds of ancient whole-genome duplications (WGDs). However, whether different types of gene duplications evolved differently after duplication remains unclear in the pear genome. In this study, we identified the different modes of gene duplication in pear. Duplicate genes derived from WGD, tandem, proximal, retrotransposed, DNA-based transposed or dispersed duplications differ in genomic distribution, gene features, selection pressure, expression divergence, regulatory divergence and biological roles. Widespread sequence, expression and regulatory divergence have occurred between duplicate genes over the 30-45 million years of evolution after the recent genome duplication in pear. The retrotransposed genes show relatively higher expression and regulatory divergence than other gene duplication modes. In contrast, WGD genes underwent a slower sequence divergence and may be influenced by abundant gene conversion events. Moreover, the different classes of duplicate genes exhibited biased functional roles. We also investigated the evolution and expansion patterns of the gene families involved in sugar and organic acid metabolism pathways, which are closely related to the fruit quality and taste in pear. Single-gene duplications largely account for the extensive expansion of gene families involved in the sorbitol metabolism pathway in pear. Gene family expansion was also detected in the sucrose metabolism pathway and tricarboxylic acid cycle pathways. Thus, this study provides insights into the evolutionary fates of duplicated genes.

Mutagenesis of the FAE1 genes significantly changes fatty acid composition in seeds of Camelina sativa

Camelina sativa is a re-emerging low-input oilseed crop that has great potentials. It is necessary to ameliorate camelina oils for optimized fatty acid composition that can meet different application requirements. Camelina seed contains significant amounts of C20-C24 very long-chain fatty acids (VLCFAs) that may not be desirable. We demonstrated that these VLCFAs can be effectively reduced by deactivating the Fatty Acid Elongase1 (FAE1) in camelina. The allohexaploid camelina contains three alleles of FAE1 genes. Ethyl methanesulfonate (EMS) induced mutation at the FAE1-B gene caused over 60% reduction of VLCFAs in seed. Homozygous knockout mutants were successfully created in a single generation by simultaneously targeting three FAE1 alleles using the CRISPR technology with an egg cell-specific Cas9 expression. VLCFAs were reduced to less than 2% of total fatty acids compared to over 22% in the wild type, and the C18 unsaturated fatty acids were concomitantly increased. The fae1 mutants were indistinguishable from wild type in seed physiology and plant growth. This study demonstrated that the CRISPR/Cas9 technology can be effectively applied to the polyploid crop camelina to rapidly obtain desired traits such as optimal fatty acid composition in its seed oil. Knocking out FAE1 also provides a means to increase the levels of oleic acid or α-linolenic acid in camelina oils that are desirable for industrial or food/feed uses.

Simultaneous Targeting of Multiple Gene Homeologs to Alter Seed Oil Production in Camelina sativa

The ability to transform Camelina sativa easily with biosynthetic enzymes derived from other plants has made this oil seed crop an ideal platform for the production of unusual lipids valuable for different applications. However, in addition to expressing transgenic enzymes, the suppression of endogenous enzyme activity to reduce competition for common substrates or cofactors is also required to enhance the production of target compounds. As camelina possesses a relatively undifferentiated hexaploid genome, up to three gene homeologs can code for any particular enzymatic activity, complicating efforts to alter endogenous biosynthetic pathways. New genome editing technologies, such as that offered by the CRISPR/Cas (clustered regularly interspaced short palindromic repeats/CRISPR-associated protein) system, offer the capability to introduce mutations into specifically targeted genomic sites. Here, by using a carefully designed guide RNA identical to all three homeologs, we demonstrate the ability of the CRISPR/Cas genome editing system to introduce mutations in all three CsDGAT1 or CsPDAT1 homeologous genes important for triacylglycerol (TAG) synthesis in developing seeds. Sequence analysis from transgenic T1 plants revealed that each CsDGAT1 or each CsPDAT1 homeolog was altered by multiple mutations, resulting in a genetic mosaic in the plants. Interestingly, seed harvested from both CsDGAT1- and CsPDAT1-targeted lines was often shrunken and wrinkled. Further, lipid analysis revealed that many lines produced seed with reduced oil content and altered fatty acid composition, consistent with the role of the targeted genes in seed oil biosynthesis. The CRISPR/Cas system therefore represents a useful method to alter endogenous biosynthetic pathways efficiently in polyploid species such as camelina.

Selective gene dosage by CRISPR-Cas9 genome editing in hexaploid Camelina sativa

In many plant species, gene dosage is an important cause of phenotype variation. Engineering gene dosage, particularly in polyploid genomes, would provide an efficient tool for plant breeding. The hexaploid oilseed crop Camelina sativa, which has three closely related expressed subgenomes, is an ideal species for investigation of the possibility of creating a large collection of combinatorial mutants. Selective, targeted mutagenesis of the three delta-12-desaturase (FAD2) genes was achieved by CRISPR-Cas9 gene editing, leading to reduced levels of polyunsaturated fatty acids and increased accumulation of oleic acid in the oil. Analysis of mutations over four generations demonstrated the presence of a large variety of heritable mutations in the three isologous CsFAD2 genes. The different combinations of single, double and triple mutants in the T3 generation were isolated, and the complete loss-of-function mutants revealed the importance of delta-12-desaturation for Camelina development. Combinatorial association of different alleles for the three FAD2 loci provided a large diversity of Camelina lines with various lipid profiles, ranging from 10% to 62% oleic acid accumulation in the oil. The different allelic combinations allowed an unbiased analysis of gene dosage and function in this hexaploid species, but also provided a unique source of genetic variability for plant breeding.