The 'TranSeq' 3'-end sequencing method for high-throughput transcriptomics and gene space refinement in plant genomes

Cystathionine γ-synthase expression in seeds alters metabolic and DNA methylation profiles in Arabidopsis

Arabidopsis (Arabidopsis thaliana) seeds expressing the feedback-insensitive form of cystathionine γ-synthase (AtD-CGS), the key gene of methionine (Met) synthesis, under the control of a seed-specific phaseolin promoter (SSE plants) show a significant increase in Met content. This elevation is accompanied by increased levels of other amino acids (AAs), sugars, total protein, and starch, which are important from a nutritional aspect. Here, we investigated the mechanism behind this phenomenon. Gas chromatography-mass spectrometry (GC-MS) analysis of SSE leaves, siliques, and seeds collected at 3 different developmental stages showed high levels of Met, AAs, and sugars compared to the control plants. A feeding experiment with isotope-labeled AAs showed an increased flux of AAs from nonseed tissues toward the developing seeds of SSE. Transcriptome analysis of leaves and seeds displayed changes in the status of methylation-related genes in SSE plants that were further validated by methylation-sensitive enzymes and colorimetric assay. These results suggest that SSE leaves have higher DNA methylation rates than control plants. This occurrence apparently led to accelerated senescence, together with enhanced monomer synthesis, which further resulted in increased transport of monomers from the leaves toward the seeds. The developing seeds of SSE plants, however, show reduced Met levels and methylation rates. The results provide insights into the role of Met in DNA methylation and gene expression and how Met affects the metabolic profile of the plant.

Short-term transcriptomic analysis at organ scale reveals candidate genes involved in low N responses in NUE-contrasting tomato genotypes

BACKGROUND

Understanding the complex regulatory network underlying plant nitrogen (N) responses associated with high Nitrogen Use Efficiency (NUE) is one of the main challenges for sustainable cropping systems. Nitrate (NO₃ ^-), acting as both an N source and a signal molecule, provokes very fast transcriptome reprogramming, allowing plants to adapt to its availability. These changes are genotype- and tissue-specific; thus, the comparison between contrasting genotypes is crucial to uncovering high NUE mechanisms.

METHODS

Here, we compared, for the first time, the spatio-temporal transcriptome changes in both root and shoot of two NUE contrasting tomato genotypes, Regina Ostuni (high-NUE) and UC82 (low-NUE), in response to short-term (within 24 h) low (LN) and high (HN) NO₃ ^- resupply.

RESULTS

Using time-series transcriptome data (0, 8, and 24 h), we identified 395 and 482 N-responsive genes differentially expressed (DEGs) between RO and UC82 in shoot and root, respectively. Protein kinase signaling plant hormone signal transduction, and phenylpropanoid biosynthesis were the main enriched metabolic pathways in shoot and root, respectively, and were upregulated in RO compared to UC82. Interestingly, several N transporters belonging to NRT and NPF families, such as NRT2.3, NRT2.4, NPF1.2, and NPF8.3, were found differentially expressed between RO and UC82 genotypes, which might explain the contrasting NUE performances. Transcription factors (TFs) belonging to several families, such as ERF, LOB, GLK, NFYB, ARF, Zinc-finger, and MYB, were differentially expressed between genotypes in response to LN. A complementary Weighted Gene Co-expression Network Analysis (WGCNA) allowed the identification of LN-responsive co-expression modules in RO shoot and root. The regulatory network analysis revealed candidate genes that might have key functions in short-term LN regulation. In particular, an asparagine synthetase (ASNS), a CBL-interacting serine/threonine-protein kinase 1 (CIPK1), a cytokinin riboside 5'-monophosphate phosphoribohydrolase (LOG8), a glycosyltransferase (UGT73C4), and an ERF2 were identified in the shoot, while an LRR receptor-like serine/threonine-protein kinase (FEI1) and two TFs NF-YB5 and LOB37 were identified in the root.

DISCUSSION

Our results revealed potential candidate genes that independently and/or concurrently may regulate short-term low-N response, suggesting a key role played by cytokinin and ROS balancing in early LN regulation mechanisms adopted by the N-use efficient genotype RO.

The Alteration of Tomato Chloroplast Vesiculation Positively Affects Whole-Plant Source-Sink Relations and Fruit Metabolism under Stress Conditions

Changes in climate conditions can negatively affect the productivity of crop plants. They can induce chloroplast degradation (senescence), which leads to decreased source capacity, as well as decreased whole-plant carbon/nitrogen assimilation and allocation. The importance, contribution and mechanisms of action regulating source-tissue capacity under stress conditions in tomato (Solanum lycopersicum) are not well understood. We hypothesized that delaying chloroplast degradation by altering the activity of the tomato chloroplast vesiculation (CV) under stress would lead to more efficient use of carbon and nitrogen and to higher yields. Tomato CV is upregulated under stress conditions. Specific induction of CV in leaves at the fruit development stage resulted in stress-induced senescence and negatively affected fruit yield, without any positive effects on fruit quality. Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated protein 9 (CRISPR/CAS9) knockout CV plants, generated using a near-isogenic tomato line with enhanced sink capacity, exhibited stress tolerance at both the vegetative and the reproductive stages, leading to enhanced fruit quantity, quality and harvest index. Detailed metabolic and transcriptomic network analysis of sink tissue revealed that the l-glutamine and l-arginine biosynthesis pathways are associated with stress-response conditions and also identified putative novel genes involved in tomato fruit quality under stress. Our results are the first to demonstrate the feasibility of delayed stress-induced senescence as a stress-tolerance trait in a fleshy fruit crop, to highlight the involvement of the CV pathway in the regulation of source strength under stress and to identify genes and metabolic pathways involved in increased tomato sink capacity under stress conditions.

Transcriptomics reveals the molecular mechanisms of flesh colour differences in eggplant (Solanum melongena)

BACKGROUND

Fruit flesh colour is not only an important commodity attribute of eggplant but is also closely related to maturity. However, very little is known about its formation mechanism in eggplant.

RESULTS

Two inbred lines of eggplant, green 'NC7' and white 'BL', were used in this study to explain the differences in flesh colour. Transcriptome sequencing results revealed a total of 3304 differentially expressed genes (DEGs) in NC7 vs. BL. Of the DEGs obtained, 2050 were higher and 1254 were lower in BL. These DEGs were annotated to 126 pathways, where porphyrin and chlorophyll metabolism, flavonoid biosynthesis, and photosynthesis-antenna proteins play vital roles in the colour formation of eggplant flesh. At the same time, Gene Ontology (GO) enrichment significance analysis showed that a large number of unigenes involved in the formation of chloroplast structure were lower in BL, which indicated that the formation of chloroplasts in white-fleshed eggplant was blocked. This was confirmed by transmission electron microscopy (TEM), which found only leucoplasts but no chloroplasts in the flesh cells of white-fleshed eggplant. Several genes encoding ERF and bHLH transcription factors were predicted to participate in the regulation of chlorophyll biosynthetic genes.

CONCLUSIONS

The results of this study indicated that differences in the gene expression of the chlorophyll metabolic pathway were the main cause of the different flesh colour formations. These findings will increase our understanding of the genetic basis in eggplant flesh colors formation mechanism.

Uncovering Pathways Highly Correlated to NUE through a Combined Metabolomics and Transcriptomics Approach in Eggplant

Nitrogen (N) fertilization is one of the main inputs to increase crop yield and food production. However, crops utilize only 30–40% of N applied; the remainder is leached into the soil, causing environmental and health damage. In this scenario, the improvement of nitrogen-use efficiency (NUE) will be an essential strategy for sustainable agriculture. Here, we compared two pairs of NUE-contrasting eggplant (Solanum melongena L.) genotypes, employing GC-MS and UPLC-qTOF-MS-based technologies to determine the differential profiles of primary and secondary metabolites in root and shoot tissues, under N starvation as well as at short- and long-term N-limiting resupply. Firstly, differences in the primary metabolism pathways of shoots related to alanine, aspartate and glutamate; starch, sucrose and glycine; serine and threonine; and in secondary metabolites biosynthesis were detected. An integrated analysis between differentially accumulated metabolites and expressed transcripts highlighted a key role of glycine accumulation and the related glyA transcript in the N-use-efficient genotypes to cope with N-limiting stress. Interestingly, a correlation between both sucrose synthase (SUS)- and fructokinase (scrK)-transcript abundances, as well as D-glucose and D-fructose accumulation, appeared useful to distinguish the N-use-efficient genotypes. Furthermore, increased levels of L-aspartate and L-asparagine in the N-use-efficient genotypes at short-term low-N exposure were detected. Granule-bound starch synthase (WAXY) and endoglucanase (E3.2.1.4) downregulation at long-term N stress was observed. Therefore, genes and metabolites related to these pathways could be exploited to improve NUE in eggplant.

A Myb transcription factor, PgMyb308-like, enhances the level of shikimate, aromatic amino acids, and lignins, but represses the synthesis of flavonoids and hydrolyzable tannins, in pomegranate (Punica granatum L.)

Pomegranate fruit peels are highly abundant in metabolites derived from the shikimate pathway, such as hydrolyzable tannins (HTs) and flavonoids. These metabolites are beneficial to human health (commercial juice is enriched with peel metabolites), and also protect the fruit from environmental stresses. To understand the transcriptional control of shikimate pathway-related metabolites in pomegranate, we cloned and characterized a subgroup S4 R2R3 Myb transcription factor, PgMyb308-like. Overexpressing PgMyb308-like in pomegranate hairy roots increased the accumulation of shikimate, aromatic amino acids, isoferulic acid, and total lignins, but led to reduced gallic acid and its downstream products HTs, as well as multiple flavonoids. Changes in these metabolites are supported by the increased expression of 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase and shikimate dehydrogenase 1 (PgSDH1) (the SDH isoform associated with shikimate biosynthesis), and the reduced expression of PgSDH4 (the SDH isoform suggested to produce gallic acid). Transcriptome analysis of PgMyb308-like-overexpressing hairy roots further revealed reprogramming of cell wall-related genes, while overexpression of PgMyb308-like in Arabidopsis thaliana plants uncovered its distinct role in a different genetic and metabolic background. These results together suggest that PgMyb308-like activates genes in the shikimate pathway and lignin biosynthesis, but suppresses those involved in the production of HTs and flavonoids.

Boosting the power of transcriptomics by developing an efficient gene expression profiling approach

Recent advances in plant genomics are scaling up gene expression profiling from the individual level to the population level, making transcriptomics a more powerful tool while deciphering the genome function. This study developed an efficient 3'RNA-seq method, Simplified Poly(A) Anchored Sequencing (SiPAS), to perform large-scale experiments of gene expression quantification. Aside from being cost-effective, by conducting a comprehensive performance assessment of SiPAS in hexaploid wheat, we demonstrated that SiPAS is highly sensitive, accurate, and reproducible while quantifying gene expression. Our method is anticipated to boost studies of population transcriptomics in plants and improve our understanding of genome biology.

A High-Throughput 3'-Tag RNA Sequencing for Large-Scale Time-Series Transcriptome Studies

RNA sequencing (RNA-seq) has proven invaluable for exploring gene expression variation under complex environmental cues. However, the cost of standard RNA-seq (e.g., Illumina TruSeq or NEBNext) remains a barrier for high-throughput applications. 3'-Tag RNA-seq (3'-TagSeq) is a cost-effective solution that permits large-scale experiments. Unlike standard RNA-seq, which generates sequencing libraries for full-length mRNAs, 3'-TagSeq only generates a single fragment from the 3' end of each transcript (a tag read) and quantifies gene expression by tag abundance. Consequently, 3'-TagSeq requires lower sequencing depth (~5 million reads per sample) than standard RNA-seq (~30 million reads per sample), which reduces costs and allows increased technical and biological replication in experiments. Because 3'-TagSeq is considerably cheaper than standard RNA-seq while exhibiting comparable accuracy and reproducibility, researchers focusing on gene expression levels in large or extensive time-series experiments might find 3'-TagSeq to be superior to standard RNA-seq. In this chapter, we describe 3'-TagSeq sequencing library preparation and provide example bioinformatics and statistical analyses of gene expression data.

MP3RNA-seq: Massively parallel 3' end RNA sequencing for high-throughput gene expression profiling and genotyping

Transcriptome deep sequencing (RNA-seq) has become a routine method for global gene expression profiling. However, its application to large-scale experiments remains limited by cost and labor constraints. Here we describe a massively parallel 3' end RNA-seq (MP3RNA-seq) method that introduces unique sample barcodes during reverse transcription to permit sample pooling immediately following this initial step. MP3RNA-seq allows for handling of hundreds of samples in a single experiment, at a cost of about $6 per sample for library construction and sequencing. MP3RNA-seq is effective for not only high-throughput gene expression profiling, but also genotyping. To demonstrate its utility, we applied MP3RNA-seq to 477 double haploid lines of maize. We identified 19,429 genes expressed in at least 50% of the lines and 35,836 high-quality single nucleotide polymorphisms for genotyping analysis. Armed with these data, we performed expression and agronomic trait quantitative trait locus (QTL) mapping and identified 25,797 expression QTLs for 15,335 genes and 21 QTLs for plant height, ear height, and relative ear height. We conclude that MP3RNA-seq is highly reproducible, accurate, and sensitive for high-throughput gene expression profiling and genotyping, and should be generally applicable to most eukaryotic species.

Transcriptomics reveal new insights into molecular regulation of nitrogen use efficiency in Solanum melongena

Nitrogen-use efficiency (NUE) is a complex trait of great interest in breeding programs because through its improvement, high crop yields can be maintained whilst N supply is reduced. In this study, we report a transcriptomic analysis of four NUE-contrasting eggplant (Solanum melongena) genotypes following short- and long-term exposure to low N, to identify key genes related to NUE in the roots and shoots. The differentially expressed genes in the high-NUE genotypes are involved in the light-harvesting complex and receptor, a ferredoxin-NADP reductase, a catalase and WRKY33. These genes were then used as bait for a co-expression gene network analysis in order to identify genes with the same trends in expression. This showed that up-regulation of WRKY33 triggered higher expression of a cluster of 21 genes and also of other genes, many of which were related to N-metabolism, that were able to improve both nitrogen uptake efficiency and nitrogen utilization efficiency, the two components of NUE. We also conducted an independent de novo experiment to validate the significantly higher expression of WRKY33 and its gene cluster in the high-NUE genotypes. Finally, examination of an Arabidopsis transgenic 35S::AtWRKY33 overexpression line showed that it had a bigger root system and was more efficient at taking up N from the soil, confirming the pivotal role of WRKY33 for NUE improvement.

Induced defense response in red mango fruit against Colletotrichum gloeosporioides

Mango fruit exposed to sunlight develops red skin and are more resistant to biotic and abiotic stresses. Here we show that harvested red mango fruit that was exposed to sunlight at the orchard is more resistant than green fruit to Colletotrichum gloeosporioides. LCMS analysis showed high amounts of antifungal compounds, as glycosylated flavonols, glycosylated anthocyanins, and mangiferin in red vs. green mango skin, correlated with higher antioxidant and lower ROS. However, also the green side of red mango fruit that has low levels of flavonoids was resistant, indicated induced resistance. Transcriptomes of red and green fruit inoculated on their red and green sides with C. gloeosporioides were analyzed. Overall, in red fruit skin, 2,187 genes were upregulated in response to C. gloeosporioides. On the green side of red mango, upregulation of 22 transcription factors and 33 signaling-related transcripts indicated induced resistance. The RNA-Seq analysis suggests that resistance of the whole red fruit involved upregulation of ethylene, brassinosteroid, and phenylpropanoid pathways. To conclude, red fruit resistance to fungal pathogen was related to both flavonoid toxicity and primed resistance of fruit that was exposed to light at the orchard.

Innovative transcriptome-based genotyping highlights environmentally responsive genes for phenology, growth and yield in a non-model grain legume

The narrow-leafed lupin, Lupinus angustifolius L., is a grain legume crop, cultivated both as a green manure and as a source of protein for animal feed and human food production. During its domestication process, numerous agronomic traits were improved, however, only two trait-related genes were identified hitherto, both by linkage mapping. Genome-wide association studies (GWAS), exploiting genomic sequencing, did not select any novel candidate gene. In the present study, an innovative method of 3'-end reduced representation transcriptomic profiling, a massive analysis of cDNA ends, has been used for genotyping of 126 L. angustifolius lines surveyed by field phenotyping. Significant genotype × environment interactions were identified for all phenology and yield traits analysed. Principal component analysis of population structure evidenced European domestication bottlenecks, visualized by clustering of breeding materials and cultivars. GWAS provided contribution towards deciphering vernalization pathway in legumes, and, apart from highlighting known domestication loci (Ku/Julius and mol), designated novel candidate genes for L. angustifolius traits. Early phenology was associated with genes from vernalization, cold-responsiveness and phosphatidylinositol signalling pathways whereas high yield with genes controlling photosynthesis performance and abiotic stress (drought or heat) tolerance. PCR-based toolbox was developed and validated to enable tracking desired alleles in marker-assisted selection. Narrow-leafed lupin was genotyped with an innovative method of transcriptome profiling and phenotyped for phenology, growth and yield traits in field. Early phenology was found associated with genes from cold-response, vernalization and phosphatidylinositol signalling pathways, whereas high yield with genes running photosystem II and drought or heat stress response. Key loci were supplied with PCR-based toolbox for marker-assisted selection.

Analysis of wild tomato introgression lines elucidates the genetic basis of transcriptome and metabolome variation underlying fruit traits and pathogen response

Wild tomato species represent a rich gene pool for numerous desirable traits lost during domestication. Here, we exploited an introgression population representing wild desert-adapted species and a domesticated cultivar to establish the genetic basis of gene expression and chemical variation accompanying the transfer of wild-species-associated fruit traits. Transcriptome and metabolome analysis of 580 lines coupled to pathogen sensitivity assays resulted in the identification of genomic loci associated with levels of hundreds of transcripts and metabolites. These associations occurred in hotspots representing coordinated perturbation of metabolic pathways and ripening-related processes. Here, we identify components of the Solanum alkaloid pathway, as well as genes and metabolites involved in pathogen defense and linking fungal resistance with changes in the fruit ripening regulatory network. Our results outline a framework for understanding metabolism and pathogen resistance during tomato fruit ripening and provide insights into key fruit quality traits.

Transcriptome profiling at osmotic and ionic phases of salt stress response in bread wheat uncovers trait-specific candidate genes

BACKGROUND

Bread wheat is one of the most important crops for the human diet, but the increasing soil salinization is causing yield reductions worldwide. Improving salt stress tolerance in wheat requires the elucidation of the mechanistic basis of plant response to this abiotic stress factor. Although several studies have been performed to analyze wheat adaptation to salt stress, there are still some gaps to fully understand the molecular mechanisms from initial signal perception to the onset of responsive tolerance pathways. The main objective of this study is to exploit the dynamic salt stress transcriptome in underlying QTL regions to uncover candidate genes controlling salt stress tolerance in bread wheat. The massive analysis of 3'-ends sequencing protocol was used to analyze leave samples at osmotic and ionic phases. Afterward, stress-responsive genes overlapping QTL for salt stress-related traits in two mapping populations were identified.

RESULTS

Among the over-represented salt-responsive gene categories, the early up-regulation of calcium-binding and cell wall synthesis genes found in the tolerant genotype are presumably strategies to cope with the salt-related osmotic stress. On the other hand, the down-regulation of photosynthesis-related and calcium-binding genes, and the increased oxidative stress response in the susceptible genotype are linked with the greater photosynthesis inhibition at the osmotic phase. The specific up-regulation of some ABC transporters and Na+/Ca2+ exchangers in the tolerant genotype at the ionic stage indicates their involvement in mechanisms of sodium exclusion and homeostasis. Moreover, genes related to protein synthesis and breakdown were identified at both stress phases. Based on the linkage disequilibrium blocks, salt-responsive genes within QTL intervals were identified as potential components operating in pathways leading to salt stress tolerance. Furthermore, this study conferred evidence of novel regions with transcription in bread wheat.

CONCLUSION

The dynamic transcriptome analysis allowed the comparison of osmotic and ionic phases of the salt stress response and gave insights into key molecular mechanisms involved in the salt stress adaptation of contrasting bread wheat genotypes. The leveraging of the highly contiguous chromosome-level reference genome sequence assembly facilitated the QTL dissection by targeting novel candidate genes for salt tolerance.

SUBERMAN regulates developmental suberization of the Arabidopsis root endodermis

Root endodermis, the innermost cortical layer surrounding the root vasculature, serves as the foremost barrier to water, solutes, and nutrients taken up from soil. Endodermis barrier functionality is achieved via its hydrophobic coating of lignified Casparian strips and the suberin lamellae; nonetheless the regulatory mechanisms underlying endodermis suberization are still elusive. Here, we discovered that the Arabidopsis SUBERMAN (SUB) transcription factor controls the establishment of the root suberin lamellae. Transient expression of SUB in Nicotiana benthamiana leaves resulted in the induction of heterologous suberin genes, the accumulation of suberin-type monomers, and consequent deposition of suberin-like lamellae. We demonstrate that SUB exerts its regulatory roles by transactivating promoters of suberin genes. In Arabidopsis, SUB is expressed in patchy and continuous suberization root endodermal cells, and thus roots with higher or lower expression of SUB display altered suberin polymer deposition patterns and modified composition. While these changes did not interfere with Casparian strip formation they had a substantial effect on root uptake capacity, resulting in varied root and leaf ionomic phenotypes. Gene expression profiling revealed that SUB function impacts transcriptional networks associated with suberin, phenylpropanoids, lignin, and cuticular lipid biosynthesis, as well as root transport activities, hormone signalling, and cell wall modification. Our findings highlight SUB as a regulator of root endodermis suberization during normal development, and its characterization is thus a key step towards dissecting the molecular mechanisms partaking in root endodermal barrier functionalities.

Heritable temporal gene expression patterns correlate with metabolomic seed content in developing hexaploid oat seed

Oat ranks sixth in world cereal production and has a higher content of health-promoting compounds compared with other cereals. However, there is neither a robust oat reference genome nor transcriptome. Using deeply sequenced full-length mRNA libraries of oat cultivar Ogle-C, a de novo high-quality and comprehensive oat seed transcriptome was assembled. With this reference transcriptome and QuantSeq 3' mRNA sequencing, gene expression was quantified during seed development from 22 diverse lines across six time points. Transcript expression showed higher correlations between adjacent time points. Based on differentially expressed genes, we identified 22 major temporal co-expression (TCoE) patterns of gene expression and revealed enriched gene ontology biological processes. Within each TCoE set, highly correlated transcripts, putatively commonly affected by genetic background, were clustered and termed genetic co-expression (GCoE) sets. Seventeen of the 22 TCoE sets had GCoE sets with median heritabilities higher than 0.50, and these heritability estimates were much higher than that estimated from permutation analysis, with no divergence observed in cluster sizes between permutation and non-permutation analyses. Linear regression between 634 metabolites from mature seeds and the PC1 score of each of the GCoE sets showed significantly lower p-values than permutation analysis. Temporal expression patterns of oat avenanthramides and lipid biosynthetic genes were concordant with previous studies of avenanthramide biosynthetic enzyme activity and lipid accumulation. This study expands our understanding of physiological processes that occur during oat seed maturation and provides plant breeders the means to change oat seed composition through targeted manipulation of key pathways.

Glycosylated flavonoids: fruit's concealed antifungal arsenal

Fruit defense against pathogens relies on induced and preformed mechanisms. The present contribution evaluated performed resistance of red and green mango fruit against the fungal pathogen Colletotrichum gloeosporioides and identified the main active antifungal components. High-performance liquid chromatography analysis of nonhydrolyzed mango peel extracts identified major anthocyanin peaks of glycosylated cyanidin and methylcyanidin, and flavonol peaks of glycosylated quercetin and kaempferol, which were more abundant on the 'red side' of red mango fruit. Organic extracts of red vs green mango peel were more efficient in inhibiting C. gloeosporioides. Transcriptome analysis of the mango-C. gloeosporioides interaction showed increased expression of glucosidase genes related to both fungal pathogenicity and host defense. Glucosidase treatment of organic peel extract increased its antifungal activity. Additionally, quercetin and cyanidin had significantly higher antifungal activity than their glycosylated derivatives. Peel extract volatiles treated with glucosidase had antifungal activity. GCMS analysis identified 15 volatiles after glucosidase hydrolysis, seven of them present only in red fruit. These results suggest that the fruit obtains a concealed arsenal of glycosylated flavonoids in its peel when they are hydrolyzed by β-glucosidase that is induced in both fungus and host during infection process, become more toxic to the fungal pathogen, inhibiting decay development.

MapMan4: A Refined Protein Classification and Annotation Framework Applicable to Multi-Omics Data Analysis

Genome sequences from over 200 plant species have already been published, with this number expected to increase rapidly due to advances in sequencing technologies. Once a new genome has been assembled and the genes identified, the functional annotation of their putative translational products, proteins, using ontologies is of key importance as it places the sequencing data in a biological context. Furthermore, to keep pace with rapid production of genome sequences, this functional annotation process must be fully automated. Here we present a redesigned and significantly enhanced MapMan4 framework, together with a revised version of the associated online Mercator annotation tool. Compared with the original MapMan, the new ontology has been expanded almost threefold and enforces stricter assignment rules. This framework was then incorporated into Mercator4, which has been upgraded to reflect current knowledge across the land plant group, providing protein annotations for all embryophytes with a comparably high quality. The annotation process has been optimized to allow a plant genome to be annotated in a matter of minutes. The output results continue to be compatible with the established MapMan desktop application.