Pectins are fundamental polysaccharides in the plant primary cell wall. Polygalacturonases (PGs) and pectin methylesterases (PMEs), major components of the pectin remodeling and disassembly network, are involved in cell separation processes during many stages of plant development. A comprehensive study of these genes in plants could shed light on the evolution patterns of their structural development. In this study, we conducted whole-genome annotation, molecular evolution and gene expression analyses of PGs and PMEs in Brassica rapa and 8 other plant species. A total of 100 PGs and 110 PMEs were identified in B. rapa; they primarily diverged from 12–18 MYA and PMEs were retained more than PGs. Along with another 305 PGs and 348 PMEs in the 8 species, two different expansion or evolution types were discovered: a new branch of class A PGs appeared after the split of gymnosperms and angiosperms, which led to the rapid expansion of PGs; the pro domain was obtained or lost in the proPMEs through comprehensive analyses among PME genes. In addition, the PGs and PMEs exhibit diverged expression patterns. These findings will lead to novel insight regarding functional divergence and conservation in the gene families and provide more support for molecular evolution analyses.

Ascorbic acid (AsA) is an important antioxidant in plants and an essential vitamin for humans. Extending the study of AsA-related genes from Arabidopsis thaliana to Brassica rapa could shed light on the evolution of AsA in plants and inform crop breeding. In this study, we conducted whole-genome annotation, molecular-evolution and gene-expression analyses of all known AsA-related genes in B. rapa. The nucleobase–ascorbate transporter (NAT) gene family and AsA l-galactose pathway genes were also compared among plant species. Four important insights gained are that: 1) 102 AsA-related gene were identified in B. rapa and they mainly diverged 12–18 Ma accompanied by the Brassica-specific genome triplication event; 2) during their evolution, these AsA-related genes were preferentially retained, consistent with the gene dosage hypothesis; 3) the putative proteins were highly conserved, but their expression patterns varied; and 4) although the number of AsA-related genes is higher in B. rapa than in A. thaliana, the AsA contents and the numbers of expressed genes in leaves of both species are similar, the genes that are not generally expressed may serve as substitutes during emergencies. In summary, this study provides genome-wide insights into evolutionary history and mechanisms of AsA-related genes following whole-genome triplication in B. rapa.

Ascorbic acid (AsA) participates in diverse biological processes, is regulated by multiple factors and is a potent antioxidant and cellular reductant. The D-mannose/L-galactose pathway is a major plant AsA biosynthetic pathway that is highly connected within biosynthetic networks, and generally conserved across plants. Previous work has shown that, although most genes of this pathway are expressed under standard growth conditions in Brassica rapa, some paralogs of these genes are not. We hypothesize that regulatory evolution in duplicate AsA pathway genes has occurred as an adaptation to environmental stressors, and that gene retention has been influenced by polyploidation events in Brassicas. To test these hypotheses, we explored the conservation of these genes in Brassicas and their expression patterns divergence in B. rapa. Similar retention and a high degree of gene sequence similarity were identified in B. rapa (A genome), Brassica oleracea (C genome) and Brassica napus (AC genome). However, the number of genes that encode the same type of enzymes varied among the three plant species. With the exception of GMP, which has nine genes, there were one to four genes that encoded the other enzymes. Moreover, we found that expression patterns divergence widely exists among these genes. i) VTC2 and VTC5 are paralogous genes, but only VTC5 is influenced by FLC. ii) Under light treatment, PMI1 co-regulates the AsA pool size with other D-Man/L-Gal pathway genes, whereas PMI2 is regulated only by darkness. iii) Under NaCl, Cu2+, MeJA and wounding stresses, most of the paralogs exhibit different expression patterns. Additionally, GME and GPP are the key regulatory enzymes that limit AsA biosynthesis in response to these treatments. In conclusion, our data support that the conservative and divergent expression patterns of D-Man/L-Gal pathway genes not only avoid AsA biosynthesis network instability but also allow B. rapa to better adapt to complex environments.