Evolutionary divergence in embryo and seed coat development of U's Triangle Brassica species illustrated by a spatiotemporal transcriptome atlas

Domestication has altered gene expression and secondary metabolites in pea seed coat

The mature seed in legumes consists of an embryo and seed coat. In contrast to knowledge about the embryo, we know relatively little about the seed coat. We analyzed the gene expression during seed development using a panel of cultivated and wild pea genotypes. Gene co-expression analysis identified gene modules related to seed development, dormancy, and domestication. Oxidoreductase genes were found to be important components of developmental and domestication processes. Proteomic and metabolomic analysis revealed that domestication favored proteins involved in photosynthesis and protein metabolism at the expense of seed defense. Seed coats of wild peas were rich in cell wall-bound metabolites and the protective compounds predominated in their seed coats. Altogether, we have shown that domestication altered pea seed development and modified (mostly reduced) the transcripts along with the protein and metabolite composition of the seed coat, especially the content of the compounds involved in defense. We investigated dynamic profiles of selected identified phenolic and flavonoid metabolites across seed development. These compounds usually deteriorated the palatability and processing of the seeds. Our findings further provide resources to study secondary metabolism and strategies for improving the quality of legume seeds which comprise an important part of the human protein diet.

Comparative transcriptomics of seed nourishing tissues: uncovering conserved and divergent pathways in seed plants

The evolutionary and ecological success of spermatophytes is intrinsically linked to the seed habit, which provides a protective environment for the initial development of the new generation. This environment includes an ephemeral nourishing tissue that supports embryo growth. In gymnosperms this tissue originates from the asexual proliferation of the maternal megagametophyte, while in angiosperms it is a product of fertilization, and is called the endosperm. The emergence of these nourishing tissues is of profound evolutionary value, and they are also food staples for most of the world's population. Here, using Orthofinder to infer orthologue genes among newly generated and previously published datasets, we provide a comparative transcriptomic analysis of seed nourishing tissues from species of several angiosperm clades, including those of early diverging lineages, as well as of one gymnosperm. Our results show that, although the structure and composition of seed nourishing tissues has seen significant divergence along evolution, there are signatures that are conserved throughout the phylogeny. Conversely, we identified processes that are specific to species within the clades studied, and thus illustrate their functional divergence. With this, we aimed to provide a foundation for future studies on the evolutionary history of seed nourishing structures, as well as a resource for gene discovery in future functional studies.

Phylogenomic curation of Ovate Family Proteins (OFPs) in the U's Triangle of Brassica L. indicates stress-induced growth modulation

The Ovate Family Proteins (OFPs) gene family houses a class of proteins that are involved in regulating plant growth and development. To date, there is no report of the simultaneous functional characterization of this gene family in all members of U's Triangle of Brassica. Here, we retrieved a combined total of 256 OFP protein sequences and analyzed their chromosomal localization, gene structure, conserved protein motif domains, and the pattern of cis-acting regulatory elements. The abundance of light-responsive elements like G-box, MRE, and GT1 motif suggests that OFPs are sensitive to the stimuli of light. The protein-protein interaction network analysis revealed that OFP05 and its orthologous genes were involved in regulating the process of transcriptional repression through their interaction with homeodomain transcription factors like KNAT and BLH. The presence of domains like DNA binding 2 and its superfamily speculated the involvement of OFPs in regulating gene expression. The biotic and abiotic stress, and the tissue-specific expression analysis of the RNA-seq datasets revealed that some of the genes such as BjuOFP30, and BnaOFP27, BolOFP11, and BolOFP10 were highly upregulated in seed coat at the mature stage and roots under various chemical stress conditions respectively which suggests their crucial role in plant growth and development processes. Experimental validation of prominent BnaOFPs such as BnaOFP27 confirmed their involvement in regulating gene expression under salinity, heavy metal, drought, heat, and cold stress. The GO and KEGG pathway enrichment analysis also sheds light on the involvement of OFPs in regulating plant growth and development. These findings have the potential to serve as a forerunner for future studies in terms of functionally diverse analysis of the OFP gene family in Brassica and other plant species.

Systems-level proteomics and metabolomics reveals the diel molecular landscape of diverse kale cultivars

Kale is a group of diverse Brassicaceae species that are nutritious leafy greens consumed for their abundance of vitamins and micronutrients. Typified by their curly, serrated and/or wavy leaves, kale varieties have been primarily defined based on their leaf morphology and geographic origin, despite having complex genetic backgrounds. Kale is a very promising crop for vertical farming due to its high nutritional content; however, being a non-model organism, foundational, systems-level analyses of kale are lacking. Previous studies in kale have shown that time-of-day harvesting can affect its nutritional composition. Therefore, to gain a systems-level diel understanding of kale across its wide-ranging and diverse genetic landscape, we selected nine publicly available and commercially grown kale cultivars for growth under near-sunlight LED light conditions ideal for vertical farming. We then analyzed changes in morphology, growth and nutrition using a combination of plant phenotyping, proteomics and metabolomics. As the diel molecular activities of plants drive their daily growth and development, ultimately determining their productivity as a crop, we harvested kale leaf tissue at both end-of-day (ED) and end-of-night (EN) time-points for all molecular analyses. Our results reveal that diel proteome and metabolome signatures divide the selected kale cultivars into two groups defined by their amino acid and sugar content, along with significant proteome differences involving carbon and nitrogen metabolism, mRNA splicing, protein translation and light harvesting. Together, our multi-cultivar, multi-omic analysis provides new insights into the molecular underpinnings of the diel growth and development landscape of kale, advancing our fundamental understanding of this nutritious leafy green super-food for horticulture/vertical farming applications.

Identification and characterization of the gene BraANS.A03 associated with purple leaf color in pak choi (Brassica rapa L. ssp. chinensis)

MAIN CONCLUSION

BraANS.A3 was the key gene controlling purple leaf color in pak choi, and two short fragments of promoter region in green pak choi might be interfering its normal expression. Pak choi (B. rapa L. ssp. chinensis) is an influential and important vegetable with green, yellow, or purple leaves that is cultivated worldwide. The purple leaves are rich in anthocyanins, but the underlying genetics and evolution have yet to be extensively studied. Free-hand sections of the purple leaves indicated that anthocyanins mainly accumulate throughout the adaxial and abaxial epidermal leaf cells. Segregation analyses of an F2 population of a B. rapa ssp. chinensis L. purple leaf mutant ZBC indicated that the purple trait is controlled by an incompletely dominant nuclear gene. Bulked segregant analysis (BSA) showed that the key genes controlling the trait were between 24.25 and 38.10 Mb on chromosome A03 of B. rapa. From the annotated genes, only BraA03g050560.3C, homologous to Arabidopsis AtANS, was related to the anthocyanin synthesis pathway. Genome annotation results and transcriptional sequencing analyses revealed that the BraANS.A3 gene was involved in the purple leaf trait. qRT-PCR analyses showed that BraANS.A3 was highly upregulated in ZBC but hardly expressed in the leaves of an inbred homozygous line of B. campestris ssp. chinensis L. green leaf mutant WTC, indicating that BraANS.A3 played a key role catalyzing anthocyanin synthesis in ZBC. Full-length sequence alignment of BraANS.A3 in WTC and ZBC showed that it was highly conserved in the gene region, with significant variation in the promoter region. In particular, the insertion of two short fragments of the promoter region in WTC may interfere with its normal expression. The promoter regions of ANS in six Brassica species all had multiple cis-elements involved in responses to abscisic acid, light, and stress, suggesting that ANS may be involved in multiple metabolic pathways or biological processes. Protein-protein interactions predicted that BraANS.A3 interacts with virtually all catalytic proteins in the anthocyanin synthesis pathway and has a strong relationship with Transparent Testa 8 (TT8). These results suggest that BraANS.A3 promotes anthocyanin accumulation in purple pak choi and provide new insights into the functional analysis of anthocyanin-related genes in Chinese cabbage and transcriptional regulatory networks.

Spatiotemporal transcriptomics and metabolic profiling provide insights into gene regulatory networks during lentil seed development

Lentil (Lens culinaris Medik.) is a nutritious legume with seeds rich in protein, minerals and an array of diverse specialized metabolites. The formation of a seed requires regulation and tight coordination of developmental programs to form the embryo, endosperm and seed coat compartments, which determines the structure and composition of mature seed and thus its end-use quality. Understanding the molecular and cellular events and metabolic processes of seed development is essential for improving lentil yield and seed nutritional value. However, such information remains largely unknown, especially at the seed compartment level. In this study, we generated high-resolution spatiotemporal gene expression profiles in lentil embryo, seed coat and whole seeds from fertilization through maturation. Apart from anatomic differences between the embryo and seed coat, comparative transcriptomics and weighted gene co-expression network analysis revealed embryo- and seed coat-specific genes and gene modules predominant in specific tissues and stages, which highlights distinct genetic programming. Furthermore, we investigated the dynamic profiles of flavonoid, isoflavone, phytic acid and saponin in seed compartments across seed development. Coupled with transcriptome data, we identified sets of candidate genes involved in the biosynthesis of these metabolites. The global view of the transcriptional and metabolic changes of lentil seed tissues throughout development provides a valuable resource for dissecting the genetic control of secondary metabolism and development of molecular tools for improving seed nutritional quality.

Genome-wide identification and expression analysis of the anthocyanin-related genes during seed coat development in six Brassica species

Yellow seed is one favorite trait for the breeding of Brassica oilseed crops, but the performance of seed coat color is very complicated due to the involvement of various pigments. The change of seed coat color of Brassica crops is related to the specific synthesis and accumulation of anthocyanin, and the expression level of structural genes in anthocyanin synthesis pathway is specifically regulated by transcription factors. Despite some previous reports on the regulations of seed coat color from linkage marker development, gene fine-mapping and multi-omics association analysis, the trait of Brassica crops is affected by the evolutionary events such as genome triploidization, the regulatory mechanism is still largely unknown. In this study, we identified genes related to anthocyanin synthesis in six Brassica crops in U-triangle at the genome-wide level and performed collinearity analysis. A total of 1119 anthocyanin-related genes were identified, the collinear relationship of anthocyanin-related genes on subgenomic chromosomes was the best in B. napus (AACC) and the worst in B. carinata (BBCC). The comparisons of gene expressions for anthocyanin metabolic pathways in seed coats during seed development revealed differences in its metabolism among these species. Interestingly, the R2R3-MYB transcription factors MYB5 and TT2 were differentially expressed at all eight stages of seed coat development, indicating that they might be the key genes that caused the variation of the seed coat color. The expression curve and trend analyses of the seed coat development period showed that the main reason for the unexpressed copies of MYB5 and TT2 was likely gene silencing caused by gene structural variation. These results were valuable for the genetic improvement of Brassica seed coat color, and also provided new insights into gene multicopy evolution in Brassica polyploids.

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.

Identification of NPF Family Genes in Brassica rapa Reveal Their Potential Functions in Pollen Development and Response to Low Nitrate Stress

Nitrate Transporter 1/Peptide Transporter Family (NPF) genes encode membrane transporters involved in the transport of diverse substrates. However, little is known about the diversity and functions of NPFs in Brassica rapa. In this study, 85 NPFs were identified in B. rapa (BrNPFs) which comprised eight subfamilies. Gene structure and conserved motif analysis suggested that BrNFPs were conserved throughout the genus. Stress and hormone-responsive cis-acting elements and transcription factor binding sites were identified in BrNPF promoters. Syntenic analysis suggested that tandem duplication contributed to the expansion of BrNPFs in B. rapa. Transcriptomic profiling analysis indicated that BrNPF2.6, BrNPF2.15, BrNPF7.6, and BrNPF8.9 were expressed in fertile floral buds, suggesting important roles in pollen development. Thirty-nine BrNPFs were responsive to low nitrate availability in shoots or roots. BrNPF2.10, BrNPF2.19, BrNPF2.3, BrNPF5.12, BrNPF5.16, BrNPF5.8, and BrNPF6.3 were only up-regulated in roots under low nitrate conditions, indicating that they play positive roles in nitrate absorption. Furthermore, many genes were identified in contrasting genotypes that responded to vernalization and clubroot disease. Our results increase understanding of BrNPFs as candidate genes for genetic improvement studies of B. rapa to promote low nitrate availability tolerance and for generating sterile male lines based on gene editing methods.

Comprehensive Analyses of the Histone Deacetylases Tuin (HDT) Gene Family in Brassicaceae Reveals Their Roles in Stress Response

Histone deacetylases tuin (HDT) is a plant-specific protein subfamily of histone deacetylation enzymes (HDAC) which has a variety of functions in plant development, hormone signaling and stress response. Although the HDT family's genes have been studied in many plant species, they have not been characterized in Brassicaceae. In this study, 14, 8 and 10 HDT genes were identified in Brassica napus, Brassica rapa and Brassica oleracea, respectively. According to phylogenetic analysis, the HDTs were divided into four groups: HDT1(HD2A), HDT2(HD2B), HDT3(HD2C) and HDT4(HD2D). There was an expansion of HDT2 orthologous genes in Brassicaceae. Most of the HDT genes were intron-rich and conserved in gene structure, and they coded for proteins with a nucleoplasmin-like (NPL) domain. Expression analysis showed that B. napus, B. rapa, and B. oleracea HDT genes were expressed in different organs at different developmental stages, while different HDT subgroups were specifically expressed in specific organs and tissues. Interestingly, most of the Bna/Br/BoHDT2 members were expressed in flowers, buds and siliques, suggesting they have an important role in the development of reproductive organs in Brassicaceae. Expression of BnaHDT was induced by various hormones, such as ABA and ethylene treatment, and some subgroups of genes were responsive to heat treatment. The expression of most HDT members was strongly induced by cold stress and freezing stress after non-cold acclimation, while it was slightly induced after cold acclimation. In this study, the HDT gene family of Brassicaceae was analyzed for the first time, which helps in understanding the function of BnaHDT in regulating plant responses to abiotic stresses, especially freezing stresses.

Spatiotemporal Transcriptomic Atlas of Developing Embryos and Vegetative Tissues in Flax

Flax (Linum usitatissimum L.) is an important multipurpose crop widely grown for oil and fiber. Despite recent advances in genomics, detailed gene activities during the important reproductive phase of its development are not well defined. In this study, we employed high-throughput RNA-sequencing methods to generate in-depth transcriptome profiles of flax tissues with emphasis on the reproductive phases of five key stages of embryogenesis (globular embryo, heart embryo, torpedo embryo, cotyledon embryo, and mature embryo), mature seed, and vegetative tissues viz. ovary, anther, and root. These datasets were used to establish the co-expression networks covering 36 gene modules based on the expression patterns for each gene through weighted gene co-expression network analysis (WGCNA). Functional interrogation with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) of dominantly expressed genetic modules in tissues revealed pathways involved in the development of different tissues. Moreover, the essential genes in embryo development and synthesis of storage reserves were identified based on their dynamic expression patterns. Together, this comprehensive dataset for developing embryos, mature seeds and vegetative tissues provides new insights into molecular mechanisms of seed development with potential for flax crop improvement.

Genome-wide identification of R2R3-MYB gene family and association with anthocyanin biosynthesis in Brassica species

Brassica species include important oil crops and vegetables in the world. The R2R3-MYB gene participates in a variety of plant functions, including the activation or inhibition of anthocyanin biosynthesis. Although previous studies have reported its phylogenetic relationships, gene structures, and expression patterns in Arabidopsis, the number and sequence variation of this gene family in Brassica crops and its involvement in the natural quantitative variation in anthocyanin biosynthesis regulation are still largely unknown. In this study, by using whole genome sequences and comprehensive genome-wide comparative analysis among the six cultivated Brassica species, 2120 R2R3-MYB genes were identified in six Brassica species, in total These R2R3-MYB genes were phylogenetically clustered into 12 groups. The R2R3-MYB family between A and C subgenomes showed better collinearity than between B and C and between A and B. From comparing transcriptional changes of five Brassica species with the purple and green leaves for the detection of the R2R3-MYB genes associated with anthocyanin biosynthesis, 7 R2R3-MYB genes were co-differentially expressed. The promoter and structure analysis of these genes showed that some variations between non-coding region, but they were highly conserved at the protein level and spatial structure. Co-expression analysis of anthocyanin-related genes and R2R3-MYBs indicated that MYB90 was strongly co-expressed with TT8, and they were co-expressed with structural genes F3H, LDOX, ANS and UF3GT at the same time. These results further clarified the roles of the R2R3-MYBs for leaf coloration in Brasica species, which provided new insights into the functions of the R2R3-MYB gene family in Brasica species.