Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Characterization of the WRKY Gene Family Related to Anthocyanin Biosynthesis and the Regulation Mechanism under Drought Stress and Methyl Jasmonate Treatment in Lycoris radiata

Lycoris radiata, belonging to the Amaryllidaceae family, is a well-known Chinese traditional medicinal plant and susceptible to many stresses. WRKY proteins are one of the largest families of transcription factors (TFs) in plants and play significant functions in regulating physiological metabolisms and abiotic stress responses. The WRKY TF family has been identified and investigated in many medicinal plants, but its members and functions are not identified in L. radiata. In this study, a total of 31 L. radiata WRKY (LrWRKY) genes were identified based on the transcriptome-sequencing data. Next, the LrWRKYs were divided into three major clades (Group I-III) based on the WRKY domains. A motif analysis showed the members within same group shared a similar motif component, indicating a conservational function. Furthermore, subcellular localization analysis exhibited that most LrWRKYs were localized in the nucleus. The expression pattern of the LrWRKY genes differed across tissues and might be important for Lycoris growth and flower development. There were large differences among the LrWRKYs based on the transcriptional levels under drought stress and MeJA treatments. Moreover, a total of 18 anthocyanin components were characterized using an ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis and pelargonidin-3-O-glucoside-5-O-arabinoside as well as cyanidin-3-O-sambubioside were identified as the major anthocyanin aglycones responsible for the coloration of the red petals in L. radiata. We further established a gene-to-metabolite correlation network and identified LrWRKY3 and LrWRKY27 significant association with the accumulation of pelargonidin-3-O-glucoside-5-O-arabinoside in the Lycoris red petals. These results provide an important theoretical basis for further exploring the molecular basis and regulatory mechanism of WRKY TFs in anthocyanin biosynthesis and in response to drought stress and MeJA treatment.

Transcriptome-wide characterization of bHLH transcription factor genes in Lycoris radiata and functional analysis of their response to MeJA

As one of the biggest plant specific transcription factor (TF) families, basic helix-loop-helix (bHLH) protein, plays significant roles in plant growth, development, and abiotic stress responses. However, there has been minimal research about the effects of methyl jasmonate (MeJA) treatment on the bHLH gene family in Lycoris radiata (L'Her.) Herb. In this study, based on transcriptome sequencing data, 50 putative L. radiata bHLH (LrbHLH) genes with complete open reading frames (ORFs), which were divided into 20 bHLH subfamilies, were identified. The protein motif analyses showed that a total of 10 conserved motifs were found in LrbHLH proteins and motif 1 and motif 2 were the most highly conserved motifs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of LrbHLH genes revealed their involvement in regulation of plant growth, jasmonic acid (JA) mediated signaling pathway, photoperiodism, and flowering. Furthermore, subcellular localization revealed that most LrbHLHs were located in the nucleus. Expression pattern analysis of LrbHLH genes in different tissues and at flower developmental stages suggested that their expression differed across lineages and might be important for plant growth and organ development in Lycoris. In addition, all LrbHLH genes exhibited specific spatial and temporal expression patterns under MeJA treatment. Moreover, protein-protein interaction (PPI) network analysis and yeast two-hybrid assay showed that numerous LrbHLHs could interact with jasmonate ZIM (zinc-finger inflorescence meristem) domain (JAZ) proteins. This research provides a theoretical basis for further investigation of LrbHLHs to find their functions and insights for their regulatory mechanisms involved in JA signaling pathway.

Comprehensive Analysis of Arabinogalactan Protein-Encoding Genes Reveals the Involvement of Three BrFLA Genes in Pollen Germination in Brassica rapa

Arabinogalactan proteins (AGPs) are a superfamily of hydroxyproline-rich glycoproteins that are massively glycosylated, widely implicated in plant growth and development. No comprehensive analysis of the AGP gene family has been performed in Chinese cabbage (Brassica rapa ssp. chinensis). Here, we identified a total of 293 putative AGP-encoding genes in B. rapa, including 25 classical AGPs, three lysine-rich AGPs, 30 AG-peptides, 36 fasciclin-like AGPs (FLAs), 59 phytocyanin-like AGPs, 33 xylogen-like AGPs, 102 other chimeric AGPs, two non-classical AGPs and three AGP/extensin hybrids. Their protein structures, phylogenetic relationships, chromosomal location and gene duplication status were comprehensively analyzed. Based on RNA sequencing data, we found that 73 AGP genes were differentially expressed in the floral buds of the sterile and fertile plants at least at one developmental stage in B. rapa, suggesting a potential role of AGPs in male reproductive development. We further characterized BrFLA2, BrFLA28 and BrFLA32, three FLA members especially expressed in anthers, pollen grains and pollen tubes. BrFLA2, BrFLA28 and BrFLA32 are indispensable for the proper timing of pollen germination under high relative humidity. Our study greatly extends the repertoire of AGPs in B. rapa and reveals a role for three members of the FLA subfamily in pollen germination.

A Genome-Wide Survey of MATE Transporters in Brassicaceae and Unveiling Their Expression Profiles under Abiotic Stress in Rapeseed

The multidrug and toxic compound extrusion (MATE) protein family is important in the export of toxins and other substrates, but detailed information on this family in the Brassicaceae has not yet been reported compared to Arabidopsis thaliana. In this study, we identified 57, 124, 81, 85, 130, and 79 MATE genes in A. thaliana, Brassica napus, Brassica oleracea, Brassica rapa, Brassica juncea, and Brassica nigra, respectively, which were unevenly distributed on chromosomes owing to both tandem and segmental duplication events. Phylogenetic analysis showed that these genes could be classified into four subgroups, shared high similarity and conservation within each group, and have evolved mainly through purifying selection. Furthermore, numerous B. napusMATE genes showed differential expression between tissues and developmental stages and between plants treated with heavy metals or hormones and untreated control plants. This differential expression was especially pronounced for the Group 2 and 3 BnaMATE genes, indicating that they may play important roles in stress tolerance and hormone induction. Our results provide a valuable foundation for the functional dissection of the different BnaMATE homologs in B. napus and its parental lines, as well as for the breeding of more stress-tolerant B. napus genotypes.

Comparative analysis of basic helix-loop-helix gene family among Brassica oleracea, Brassica rapa, and Brassica napus

Background

The basic helix–loop–helix (bHLH) is the second largest gene family in the plant, some members play important roles in pistil development and response to drought, waterlogging, cold stress and salt stress. The bHLH gene family has been identified in many species, except for Brassica oleracea and B. napus thus far. This study aims to identify the bHLH family members in B. oleracea, B. rapa and B. napus, and elucidate the expression, duplication, phylogeny and evolution characters of them.

Result

A total of 268 bHLH genes in B. oleracea, 440 genes in B. napus, and 251 genes in B. rapa, including 21 new bHLH members, have been identified. Subsequently, the analyses of the phylogenetic trees, conserved motifs and gene structures showed that the members in the same subfamily were highly conserved. Most Ka/Ks values of homologous gene were < 1, which indicated that these genes suffered from strong purifying selection for retention. The retention rates of BrabHLH and BolbHLH genes were 51.6 and 55.1%, respectively. The comparative expression patterns between B. rapa and B. napus showed that they had similar expression patterns in the root and contrasting patterns in the stems, leaves, and reproductive tissues. In addition, there were 41 and 30 differential expression bHLH genes under the treatments of ABA and JA, respectively, and the number of down regulation genes was significantly more than up regulation genes.

Conclusion

In the present study, we identified and performed the comparative genomics analysis of bHLH gene family among B. oleracea, B. rapa and B. napus, and also investigated their diversity. The expression patterns between B. rapa and B. napus shows that they have the similar expression pattern in the root and opposite patterns in the stems, leaves, and reproduction tissues. Further analysis demonstrated that some bHLH gene members may play crucial roles under the abiotic and biotic stress conditions. This is the first to report on the bHLH gene family analysis in B. oleracea and B. napus, which can offer useful information on the functional analysis of the bHLH gene in plants.

Identification, Evolution, and Expression Profiling of Histone Lysine Methylation Moderators in Brassica rapa

Histone modifications, such as methylation and demethylation, are vital for regulating chromatin structure, thus affecting its expression patterns. The objective of this study is to understand the phylogenetic relationships, genomic organization, diversification of motif modules, gene duplications, co-regulatory network analysis, and expression dynamics of histone lysine methyltransferases and histone demethylase in Brassica rapa. We identified 60 SET (HKMTases), 53 JmjC, and 4 LSD (HDMases) genes in B. rapa. The domain composition analysis subcategorized them into seven and nine subgroups, respectively. Duplication analysis for paralogous pairs of SET and JmjC (eight and nine pairs, respectively) exhibited variation. Interestingly, three pairs of SET exhibited Ka/Ks > 1.00 values, signifying positive selection, whereas the remaining underwent purifying selection with values less than 1.00. Furthermore, RT-PCR validation analysis and RNA-sequence data acquired on six different tissues (i.e., leaf, stem, callus, silique, flower, and root) revealed dynamic expression patterns. This comprehensive study on the abundance, classification, co-regulatory network analysis, gene duplication, and responses to heat and cold stress of SET and JmjC provides insights into the structure and diversification of these family members in B. rapa. This study will be helpful to reveal functions of these putative SET and JmjC genes in B. rapa.

Genome-Wide Identification and Characterization of Pectin Methylesterase Inhibitor Genes in Brassica oleracea

The activities of pectin methylesterases (PMEs) are regulated by pectin methylesterase inhibitors (PMEIs), which consequently control the pectin methylesterification status. However, the role of PMEI genes in Brassica oleracea, an economically important vegetable crop, is poorly understood. In this study, 95 B. oleracea PMEI (BoPMEI) genes were identified. A total of 77 syntenic ortholog pairs and 10 tandemly duplicated clusters were detected, suggesting that the expansion of BoPMEI genes was mainly attributed to whole-genome triplication (WGT) and tandem duplication (TD). During diploidization after WGT, BoPMEI genes were preferentially retained in accordance with the gene balance hypothesis. Most homologous gene pairs experienced purifying selection with ω (Ka/Ks) ratios lower than 1 in evolution. Five stamen-specific BoPMEI genes were identified by expression pattern analysis. By combining the analyses of expression and evolution, we speculated that nonfunctionalization, subfunctionalization, neofunctionalization, and functional conservation can occur in the long evolutionary process. This work provides insights into the characterization of PMEI genes in B. oleracea and contributes to the further functional studies of BoPMEI genes.

Genomic, Molecular Evolution, and Expression Analysis of Genes Encoding Putative Classical AGPs, Lysine-Rich AGPs, and AG Peptides in Brassica rapa

Arabinogalactan proteins (AGPs) belong to a class of Pro/Hyp-rich glycoproteins and are some of the most complex types of macromolecules found in plants. In the economically important plant species, Brassica rapa, only chimeric AGPs have been identified to date. This has significantly limited our understanding of the functional roles of AGPs in this plant. In this study, 64 AGPs were identified in the genome of B. rapa, including 33 classical AGPs, 28 AG peptides and three lys-rich AGPs. Syntenic gene analysis between B. rapa and A. thaliana suggested that the whole genome triplication event dominated the expansion of the AGP gene family in B. rapa. This resulted in a high retained proportion of the AGP family in the B. rapa genome, especially in the least fractionated subgenome. Phylogenetic and motif analysis classified the classical AGPs into six clades and three orphan genes, and the AG peptides into three clades and five orphan genes. Classical AGPs has a faster rate of molecular evolution than AG peptides revealed by estimation of molecular evolution rates. However, no significant differences were observed between classical AGPs and lys-rich AGPs. Under control conditions and in response to phytohormones treatment, a complete expression profiling experiment has identified five anther-specific AGPs and quite a number of AGPs responding to abscisic acid, methyl jasmonate and/or gibberellin. In this study, we presented a bioinformatics approach to identify important types of AGPs. Moreover, the association between their function and their protein structure, as well as the evolution and the expression of AGP genes were investigated, which might provide fundamental information for revealing the roles of AGPs in B. rapa.

Comparative paleogenomics of crucifers: ancestral genomic blocks revisited

A decade ago the concept of the Ancestral Crucifer Karyotype (ACK) and the definition of 24 conserved genomic blocks was presented. Subsequently, 35 cytogenetic reconstructions and/or draft genome sequences of crucifer species (members of the Brassicaceae family) have been analyzed in the context of this system; placing crucifers at the forefront of plant phylogenomics. In this review, we highlight how the ACK and genomic blocks have facilitated and guided genomic analysis of crucifers in the last 10 years and provide an update of this robust model.

Diversification and evolution of the SDG gene family in Brassica rapa after the whole genome triplication

Histone lysine methylation, controlled by the SET Domain Group (SDG) gene family, is part of the histone code that regulates chromatin function and epigenetic control of gene expression. Analyzing the SDG gene family in Brassica rapa for their gene structure, domain architecture, subcellular localization, rate of molecular evolution and gene expression pattern revealed common occurrences of subfunctionalization and neofunctionalization in BrSDGs. In comparison with Arabidopsis thaliana, the BrSDG gene family was found to be more divergent than AtSDGs, which might partly explain the rich variety of morphotypes in B. rapa. In addition, a new evolutionary pattern of the four main groups of SDGs was presented, in which the Trx group and the SUVR subgroup evolved faster than the E(z), Ash groups and the SUVH subgroup. These differences in evolutionary rate among the four main groups of SDGs are perhaps due to the complexity and variability of the regions that bind with biomacromolecules, which guide SDGs to their target loci.

Genomic and Post-Translational Modification Analysis of Leucine-Rich-Repeat Receptor-Like Kinases in Brassica rapa

Among several receptor-like kinases (RLKs), leucine-rich-repeat receptor-like kinases (LRR-RLKs) are a major group of genes that play crucial roles in growth, development and stress responses in plant systems. Given that they have several functional roles, it is important to investigate their roles in Brassica rapa. In the present study, 303 LRR-RLKs were identified in the genome of B. rapa and comparative phylogenetic analysis of 1213 combined LRR-RLKs of B. rapa, Arabidopsis thaliana, Oryza sativa and Populus trichocarpa helped us to categorize the gene family into 15 subfamilies based on their sequence and structural similarities. The chromosome localizations of 293 genes allowed the prediction of duplicates, and motif conservation and intron/exon patterns showed differences among the B. rapa LRR-RLK (BrLRR-RLK) genes. Additionally, computational function annotation and expression analysis was used to predict their possible functional roles in the plant system. Biochemical results for 11 selected genes showed variations in phosphorylation activity. Interestingly, BrBAK1 showed strong auto-phosphorylation and trans-phosphorylation on its tyrosine and threonine residues compared with AtBAK1 in previous studies. The AtBAK1 receptor kinase is involved in plant growth and development, plant innate immunity, and programmed cell death, and our results suggest that BrBAK1 might also be involved in the same functions. Another interesting result was that BrBAK1, BrBRI1, BrPEPR1 and BrPEPR2 showed activity with both anti-phosphotyrosine and anti-phosphothreonine antibodies, indicating that they might have dual-specificity kinase activity. This study provides comprehensive results for the BrLRR-RLKs, revealing expansion of the gene family through gene duplications, structural similarities and variations among the genes, and potential functional roles according to gene ontology, transcriptome profiling and biochemical analysis.

Resolution of Brassicaceae Phylogeny Using Nuclear Genes Uncovers Nested Radiations and Supports Convergent Morphological Evolution

Brassicaceae is one of the most diverse and economically valuable angiosperm families with widely cultivated vegetable crops and scientifically important model plants, such as Arabidopsis thaliana. The evolutionary history, ecological, morphological, and genetic diversity, and abundant resources and knowledge of Brassicaceae make it an excellent model family for evolutionary studies. Recent phylogenetic analyses of the family revealed three major lineages (I, II, and III), but relationships among and within these lineages remain largely unclear. Here, we present a highly supported phylogeny with six major clades using nuclear markers from newly sequenced transcriptomes of 32 Brassicaceae species and large data sets from additional taxa for a total of 55 species spanning 29 out of 51 tribes. Clade A consisting of Lineage I and Macropodium nivale is sister to combined Clade B (with Lineage II and others) and a new Clade C. The ABC clade is sister to Clade D with species previously weakly associated with Lineage II and Clade E (Lineage III) is sister to the ABCD clade. Clade F (the tribe Aethionemeae) is sister to the remainder of the entire family. Molecular clock estimation reveals an early radiation of major clades near or shortly after the Eocene–Oligocene boundary and subsequent nested divergences of several tribes of the previously polytomous Expanded Lineage II. Reconstruction of ancestral morphological states during the Brassicaceae evolution indicates prevalent parallel (convergent) evolution of several traits over deep times across the entire family. These results form a foundation for future evolutionary analyses of structures and functions across Brassicaceae.

Did androgen-binding protein paralogs undergo neo- and/or Subfunctionalization as the Abp gene region expanded in the mouse genome?

The Androgen-binding protein (Abp) region of the mouse genome contains 30 Abpa genes encoding alpha subunits and 34 Abpbg genes encoding betagamma subunits, their products forming dimers composed of an alpha and a betagamma subunit. We endeavored to determine how many Abp genes are expressed as proteins in tears and saliva, and as transcripts in the exocrine glands producing them. Using standard PCR, we amplified Abp transcripts from cDNA libraries of C57BL/6 mice and found fifteen Abp gene transcripts in the lacrimal gland and five in the submandibular gland. Proteomic analyses identified proteins corresponding to eleven of the lacrimal gland transcripts, all of them different from the three salivary ABPs reported previously. Our qPCR results showed that five of the six transcripts that lacked corresponding proteins are expressed at very low levels compared to those transcripts with proteins. We found 1) no overlap in the repertoires of expressed Abp paralogs in lacrimal gland/tears and salivary glands/saliva; 2) substantial sex-limited expression of lacrimal gland/tear expressed-paralogs in males but no sex-limited expression in females; and 3) that the lacrimal gland/tear expressed-paralogs are found exclusively in ancestral clades 1, 2 and 3 of the five clades described previously while the salivary glands/saliva expressed-paralogs are found only in clade 5. The number of instances of extremely low levels of transcription without corresponding protein production in paralogs specific to tears and saliva suggested the role of subfunctionalization, a derived condition wherein genes that may have been expressed highly in both glands ancestrally were down-regulated subsequent to duplication. Thus, evidence for subfunctionalization can be seen in our data and we argue that the partitioning of paralog expression between lacrimal and salivary glands that we report here occurred as the result of adaptive evolution.