Genome-wide identification, functional prediction, and evolutionary analysis of the R2R3-MYB superfamily in Brassica napus

Genome-Wide Analysis of the Common Fig (Ficus carica L.) R2R3-MYB Genes Reveals Their Structure, Evolution, and Roles in Fruit Color Variation

The R2R3-MYB transcription factor (TF) family is crucial for regulating plant growth, stress response, and fruit ripening. Although this TF family has been examined in a multitude of plants, the R2R3-MYB TFs in Ficus carica, a Mediterranean fruit species, have yet to be characterized. This study identified and classified 63 R2R3-MYB genes (FcMYB1 to FcMYB63) in the F. carica genome. We analyzed these genes for physicochemical properties, conserved motifs, phylogenetic relationships, gene architecture, selection pressure, and gene expression profiles and networks. The genes were classified into 29 clades, with members of the same clade showing similar exon-intron structures and motif compositions. Of the 54 orthologous gene pairs shared with mulberry (Morus notabilis), 52 evolved under negative selection, while two pairs (FcMYB55/MnMYB20 and FcMYB59/MnMYB31) experienced diversifying selection. RNA-Seq analysis showed that FcMYB26, FcMYB33, and FcMYB34 were significantly overexpressed in fig fruit peel during maturation phase III. Weighted gene co-expression network analysis (WGCNA) indicated that these genes are part of an expression module associated with the anthocyanin pathway. RT-qPCR validation confirmed these findings and revealed that the Tunisian cultivars 'Zidi' and 'Soltani' have cultivar-specific R2R3-FcMYB genes highly overexpressed during the final stage of fruit maturation and color acquisition. These genes likely influence cultivar-specific pigment synthesis. This study provides a comprehensive overview of the R2R3-MYB TF family in fig, offering a framework for selecting genes related to fruit peel color in breeding programs.

Genome-wide analysis of MYB transcription factor family and AsMYB1R subfamily contribution to ROS homeostasis regulation in Avena sativa under PEG-induced drought stress

BACKGROUND

The myeloblastosis (MYB) transcription factor (TF) family is one of the largest and most important TF families in plants, playing an important role in a life cycle and abiotic stress.

RESULTS

In this study, 268 Avena sativa MYB (AsMYB) TFs from Avena sativa were identified and named according to their order of location on the chromosomes, respectively. Phylogenetic analysis of the AsMYB and Arabidopsis MYB proteins were performed to determine their homology, the AsMYB1R proteins were classified into 5 subgroups, and the AsMYB2R proteins were classified into 34 subgroups. The conserved domains and gene structure were highly conserved among the subgroups. Eight differentially expressed AsMYB genes were screened in the transcriptome of transcriptional data and validated through RT-qPCR. Three genes in AsMYB2R subgroup, which are related to the shortened growth period, stomatal closure, and nutrient and water transport by PEG-induced drought stress, were investigated in more details. The AsMYB1R subgroup genes LHY and REV 1, together with GST, regulate ROS homeostasis to ensure ROS signal transduction and scavenge excess ROS to avoid oxidative damage.

CONCLUSION

The results of this study confirmed that the AsMYB TFs family is involved in the homeostatic regulation of ROS under drought stress. This lays the foundation for further investigating the involvement of the AsMYB TFs family in regulating A. sativa drought response mechanisms.

The R2R3 MYB gene TaMYB305 positively regulates anther and pollen development in thermo-sensitive male-sterility wheat with Aegilops kotschyi cytoplasm

MAIN CONCLUSION

The loss of TaMYB305 function down-regulated the expression of jasmonic acid synthesis pathway genes, which may disturb the jasmonic acid synthesis, resulting in abnormal pollen development and reduced fertility. The MYB family, as one of the largest transcription factor families found in plants, regulates plant development, especially the development of anthers. Therefore, it is important to identify potential MYB transcription factors associated with pollen development and to study its role in pollen development. Here, the transcripts of an R2R3 MYB gene TaMYB305 from KTM3315A, a thermo-sensitive cytoplasmic male-sterility line with Aegilops kotschyi cytoplasm (K-TCMS) wheat, was isolated. Quantitative real-time PCR (qRT-PCR) and promoter activity analysis revealed that TaMYB305 was primarily expressed in anthers. The TaMYB305 protein was localized in the nucleus, as determined by subcellular localization analysis. Our data demonstrated that silencing of TaMYB305 was related to abnormal development of stamen, including anther indehiscence and pollen abortion in KAM3315A plants. In addition, TaMYB305-silenced plants exhibited alterations in the transcriptional levels of genes involved in the synthesis of jasmonic acid (JA), indicating that TaMYB305 may regulate the expression of genes related to JA synthesis and play an important role during anther and pollen development of KTM3315A. These results provide novel insight into the function and molecular mechanism of R2R3-MYB genes in pollen development.

Systematic investigation of the R2R3-MYB gene family in Aquilaria sinensis reveals a transcriptional repressor AsMYB054 involved in 2-(2-phenylethyl)chromone biosynthesis

Trees in the genus Aquilaria produce agarwood, a valuable resin used in medicine, perfumes, and incense. 2-(2-Phenethyl)chromones (PECs) are characteristic components of agarwood; however, molecular mechanisms underlying PEC biosynthesis and regulation remain largely unknown. The R2R3-MYB transcription factors play important regulatory roles in the biosynthesis of various secondary metabolites. In this study, 101 R2R3-MYB genes in Aquilaria sinensis were systematically identified and analyzed at the genome-wide level. Transcriptomic analysis revealed that 19 R2R3-MYB genes were significantly regulated by an agarwood inducer, and showed significant correlations with PEC accumulation. Expression and evolutionary analyses revealed that AsMYB054, a subgroup 4 R2R3-MYB, was negatively correlated with PEC accumulation. AsMYB054 was located in the nucleus and functioned as a transcriptional repressor. Moreover, AsMYB054 could bind to the promoters of the PEC biosynthesis related genes AsPKS02 and AsPKS09, and inhibit their transcriptional activity. These findings suggested that AsMYB054 functions as a negative regulator of PEC biosynthesis via the inhibition of AsPKS02 and AsPKS09 in A. sinensis. Our results provide a comprehensive understanding of the R2R3-MYB subfamily in A. sinensis and lay a foundation for further functional analyses of R2R3-MYB genes in PEC biosynthesis.

Genome-wide investigation and expression profiling of LOR gene family in rapeseed under salinity and ABA stress

The Brassica napus (B. napus) LOR (Lurp-One-Related) gene family is a little-known gene family characterized by a conserved LOR domain in the proteins. Limited research in Arabidopsis showed that LOR family members played important roles in Hyaloperonospora parasitica (Hpa) defense. Nevertheless, there is a paucity of research investigating the role of the LOR gene family towards their responses to abiotic stresses and hormone treatments. This study encompassed a comprehensive survey of 56 LOR genes in B. napus, which is a prominent oilseed crop that holds substantial economic significance in China, Europe, and North America. Additionally, the study evaluated the expression profiles of these genes in response to salinity and ABA stress. Phylogenetic analysis showed that 56 BnLORs could be divided into 3 subgroups (8 clades) with uneven distribution on 19 chromosomes. 37 out of 56 BnLOR members have experienced segmental duplication and 5 of them have undergone tandem repeats events with strong evidence of purifying selection. Cis-regulatory elements (CREs) analysis indicated that BnLORs involved in process such as light response, hormone response, low temperature response, heat stress response, and dehydration response. The expression pattern of BnLOR family members revealed tissue specificity. RNA-Seq and qRT-PCR were used to validate BnLOR gene expression under temperature, salinity and ABA stress, revealing that most BnLORs showed inducibility. This study enhanced our comprehension of the B. napus LOR gene family and could provide valuable information for identifying and selecting genes for stress resistant breeding.

Cloning and Functional Characterization of Cold-Inducible MYB-like 17 Transcription Factor in Rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) is an important crop for edible oil, vegetables, and biofuel. Rapeseed growth and development require a minimum temperature of ~1-3 °C. Notably, frost damage occurs during overwintering, posing a serious threat to the productivity and yield of rapeseed. MYB proteins are important transcription factors (TFs) in plants, and have been proven to be involved in the regulation of stress responses. However, the roles of the MYB TFs in rapeseed under cold stress conditions are yet to be fully elucidated. To better understand the molecular mechanisms of one MYB-like 17 gene, BnaMYBL17, in response to low temperature, the present study found that the transcript level of BnaMYBL17 is induced by cold stress. To characterize the gene's function, the 591 bp coding sequence (CDS) from rapeseed was isolated and stably transformed into rapeseed. The further functional analysis revealed significant sensitivity in BnaMYBL17 overexpression lines (BnaMYBL17-OE) after freezing stress, suggesting its involvement in freezing response. A total of 14,298 differentially expressed genes relative to freezing response were found based on transcriptomic analysis of BnaMYBL17-OE. Overall, 1321 candidate target genes were identified based on differential expression, including Phospholipases C1 (PLC1), FCS-like zinc finger 8 (FLZ8), and Kinase on the inside (KOIN). The qPCR results confirmed that the expression levels of certain genes showed fold changes ranging from two to six when compared between BnaMYBL17-OE and WT lines after exposure to freezing stress. Furthermore, verification indicated that BnaMYBL17 affects the promoter of BnaPLC1, BnaFLZ8, and BnaKOIN genes. In summary, the results suggest that BnaMYBL17 acts as a transcriptional repressor in regulating certain genes related to growth and development during freezing stress. These findings provide valuable genetic and theoretical targets for molecular breeding to enhance freezing tolerance in rapeseed.

Genome-Wide Analysis of MYB Transcription Factor Gene Superfamily Reveals BjPHL2a Involved in Modulating the Expression of BjCHI1 in Brassica juncea

Brassica juncea is an economically important vegetable and oilseed crop. The MYB transcription factor superfamily is one of the largest transcription factor families in plants, and plays crucial roles in regulating the expression of key genes involved in a variety of physiological processes. However, a systematic analysis of the MYB transcription factor genes in Brassica juncea (BjMYB) has not been performed. In this study, a total of 502 BjMYB superfamily transcription factor genes were identified, including 23 1R-MYBs, 388 R2R3-MYBs, 16 3R-MYBs, 4 4R-MYBs, 7 atypical MYBs, and 64 MYB-CCs, which is approximately 2.4-fold larger than that of AtMYBs. Phylogenetic relationship analysis revealed that the MYB-CC subfamily consists of 64 BjMYB-CC genes. The expression pattern of members of PHL2 subclade homologous genes in Brassica juncea (BjPHL2) after Botrytis cinerea infection were determined, and BjPHL2a was isolated from a yeast one-hybrid screen with the promoter of BjCHI1 as bait. BjPHL2a was found to localize mainly in the nucleus of plant cells. An EMSA assay confirmed that BjPHL2a binds to the Wbl-4 element of BjCHI1. Transiently expressed BjPHL2a activates expression of the GUS reporter system driven by a BjCHI1 mini-promoter in tobacco (Nicotiana benthamiana) leaves. Taken together, our data provide a comprehensive evaluation of BjMYBs and show that BjPHL2a, one of the members of BjMYB-CCs, functions as a transcription activator by interacting with the Wbl-4 element in the promoter of BjCHI1 for targeted gene-inducible expression.

Identification and Phylogenetic Analysis of the R2R3-MYB Subfamily in Brassica napus

The R2R3-MYB sub-family proteins are composed of most members of MYB (v-Myb avian myeloblastosis viral oncogene homolog) protein, a plant-specific transcription factor (TF) that is classified into four classes depending on the number of MYB repeats. R2R3-MYB TFs are involved in physiological and biochemical processes. However, the functions of the Brassica napus R2R3-MYB genes are still mainly unknown. In this study, 35 Brassica napus MYB (BnaMYB) genes were screened in the genome of Brassica napus, and details about their physical and chemical characteristics, evolutionary relationships, chromosome locations, gene structures, three-dimensional protein structures, cis-acting promoter elements, and gene duplications were uncovered. The BnaMYB genes have undergone segmental duplications and positive selection pressure, according to evolutionary studies. The same subfamilies have similar intron-exon patterns and motifs, according to the genes' structure and conserved motifs. Additionally, through cis-element analysis, many drought-responsive and other stress-responsive cis-elements have been found in the promoter regions of the BnaMYB genes. The expression of the BnaMYB gene displays a variety of tissue-specific patterns. Ten lignin-related genes were chosen for drought treatment. Our research screened four genes that showed significant upregulation under drought stress, and thus may be important drought-responsive genes. The findings lay a new foundation for understanding the complex mechanisms of BnaMYB in multiple developmental stages and pathways related to drought stress in rapeseed.

Genome-wide identification and expression analysis of MYB gene family under nitrogen stress in Panax notoginseng

The myeloblastosis (MYB) gene family, involved in regulating many important physiological and biochemical processes, is one of the largest transcript factor superfamilies in plants. Since the identification of genome sequencing of Panax notoginseng has been completed, there was little known about the whole genome of its specific MYB gene family and the response to abiotic stresses, in consideration of the excessive application of nitrogen fertilizers in P. notoginseng. In this study, 123 PnMYB genes (MYB genes of P. notoginseng) have been identified and divided into 3 subfamilies by the phylogenetic analysis. These PnMYB genes were unevenly located on 12 chromosomes. Meanwhile, the gene structure and protein conserved domain were established by MEME Suite. The analysis of collinear relationships reflected that there were 121 homologous genes between P. notoginseng and Arabidopsis and 30 between P. notoginseng and rice. Moreover, cis-acting elements of PnMYB gene promoters were predicted which indicated that PnMYBs are involved in biotic, abiotic stress, and hormone induction. The expressions of PnMYB transcription factors in its roots, flowers, and leaves were detected by qRT-PCR and they had tissue-specific expressions and related to the growth of different tissues. Under nitrogen stress, MYB transcription factors had great feedback. Ten R2R3-MYB subfamily genes were significantly induced and indicated the possible function of protecting P. notoginseng from excess nitrogen. With further knowledge on identification of PnMYB gene related to tissue selectivity and abiotic stresses, this study laid the foundation for the functional development of PnMYB gene family and improved the cultivation of P. notoginseng.

Downregulation of Brassica napus MYB69 (BnMYB69) increases biomass growth and disease susceptibility via remodeling phytohormone, chlorophyll, shikimate and lignin levels

MYB transcription factors are major actors regulating plant development and adaptability. Brassica napus is a staple oil crop and is hampered by lodging and diseases. Here, four B. napus MYB69 (BnMYB69s) genes were cloned and functionally characterized. They were dominantly expressed in stems during lignification. BnMYB69 RNA interference (BnMYB69i) plants showed considerable changes in morphology, anatomy, metabolism and gene expression. Stem diameter, leaves, roots and total biomass were distinctly larger, but plant height was significantly reduced. Contents of lignin, cellulose and protopectin in stems were significantly reduced, accompanied with decrease in bending resistance and Sclerotinia sclerotiorum resistance. Anatomical detection observed perturbation in vascular and fiber differentiation in stems, but promotion in parenchyma growth, accompanied with changes in cell size and cell number. In shoots, contents of IAA, shikimates and proanthocyanidin were reduced, while contents of ABA, BL and leaf chlorophyll were increased. qRT-PCR revealed changes in multiple pathways of primary and secondary metabolisms. IAA treatment could recover many phenotypes and metabolisms of BnMYB69i plants. However, roots showed trends opposite to shoots in most cases, and BnMYB69i phenotypes were light-sensitive. Conclusively, BnMYB69s might be light-regulated positive regulators of shikimates-related metabolisms, and exert profound influences on various internal and external plant traits.

A Novel R2R3-MYB Transcription Factor SbMYB12 Positively Regulates Baicalin Biosynthesis in Scutellaria baicalensis Georgi

Scutellaria baicalensis Georgi is an annual herb from the Scutellaria genus that has been extensively used as a traditional medicine for over 2000 years in China. Baicalin and other flavonoids have been identified as the principal bioactive ingredients. The biosynthetic pathway of baicalin in S. baicalensis has been elucidated; however, the specific functions of R2R3-MYB TF, which regulates baicalin synthesis, has not been well characterized in S. baicalensis to date. Here, a S20 R2R3-MYB TF (SbMYB12), which encodes 263 amino acids with a length of 792 bp, was expressed in all tested tissues (mainly in leaves) and responded to exogenous hormone methyl jasmonate (MeJA) treatment. The overexpression of SbMYB12 significantly promoted the accumulation of flavonoids such as baicalin and wogonoside in S. baicalensis hairy roots. Furthermore, biochemical experiments revealed that SbMYB12 is a nuclear-localized transcription activator that binds to the SbCCL7-4, SbCHI-2, and SbF6H-1 promoters to activate their expression. These results illustrate that SbMYB12 positively regulates the generation of baicalin and wogonoside. In summary, this work revealed a novel S20 R2R3-MYB regulator and enhances our understanding of the transcriptional and regulatory mechanisms of baicalin biosynthesis, as well as sheds new light on metabolic engineering in S. baicalensis.

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.

Genome-wide identification and genetic characterization of the CaMYB family and its response to five types of heavy metal stress in hot pepper (Capsicum annuum cv. CM334)

MYB proteins play a crucial role in plant growth and development and stress responses. In this study, 160 members of the MYB gene family from the pepper genome database were used to analyze gene structures, chromosome localization, collinearity, genetic affinity and expression in response to heavy metals. The results identified R2R3-MYB members and further phylogenetically classified them into 35 subgroups based on highly conserved gene structures and motifs. Collinearity analysis showed that segmental duplication events played a crucial role in the functional expansion of the CaMYB gene family by intraspecific collinearity, and at least 12 pairs of CaMYB genes existed between species prior to the differentiation between monocots and dicots. Moreover, the upstream CaMYB genes were mainly localized to the phytohormone elements ABRE and transcription factor elements MYB and MYC. Further analysis revealed that MYB transcription factors were closely associated with a variety of abiotic stress-related proteins (e.g., MAC-complex and SKIP). Under the stress of five metal ions, Cd²⁺, Cu²⁺, Pb²⁺, Zn²⁺, and Fe³⁺, the expression levels of some CaMYB family genes were upregulated. Of these genes, pairing homologous 1 (PH-1), PH-13, and PH-15 in the roots of Capsicum annuum were upregulated to the greatest extent, indicating that these three MYB family members are particularly sensitive to these five metals. This study provides a theoretical reference for the analysis of the molecular regulatory mechanism of MYB family genes in mediating the response to heavy metals in plants. This study reveals the mode of interaction between MYB and a variety of abiotic stress proteins and clarifies the biological functions of CaMYB family members in the regulation of heavy metal stress.

Adaptation Strategies to Improve the Resistance of Oilseed Crops to Heat Stress Under a Changing Climate: An Overview

Temperature is one of the decisive environmental factors that is projected to increase by 1. 5°C over the next two decades due to climate change that may affect various agronomic characteristics, such as biomass production, phenology and physiology, and yield-contributing traits in oilseed crops. Oilseed crops such as soybean, sunflower, canola, peanut, cottonseed, coconut, palm oil, sesame, safflower, olive etc., are widely grown. Specific importance is the vulnerability of oil synthesis in these crops against the rise in climatic temperature, threatening the stability of yield and quality. The natural defense system in these crops cannot withstand the harmful impacts of heat stress, thus causing a considerable loss in seed and oil yield. Therefore, a proper understanding of underlying mechanisms of genotype-environment interactions that could affect oil synthesis pathways is a prime requirement in developing stable cultivars. Heat stress tolerance is a complex quantitative trait controlled by many genes and is challenging to study and characterize. However, heat tolerance studies to date have pointed to several sophisticated mechanisms to deal with the stress of high temperatures, including hormonal signaling pathways for sensing heat stimuli and acquiring tolerance to heat stress, maintaining membrane integrity, production of heat shock proteins (HSPs), removal of reactive oxygen species (ROS), assembly of antioxidants, accumulation of compatible solutes, modified gene expression to enable changes, intelligent agricultural technologies, and several other agronomic techniques for thriving and surviving. Manipulation of multiple genes responsible for thermo-tolerance and exploring their high expressions greatly impacts their potential application using CRISPR/Cas genome editing and OMICS technology. This review highlights the latest outcomes on the response and tolerance to heat stress at the cellular, organelle, and whole plant levels describing numerous approaches applied to enhance thermos-tolerance in oilseed crops. We are attempting to critically analyze the scattered existing approaches to temperature tolerance used in oilseeds as a whole, work toward extending studies into the field, and provide researchers and related parties with useful information to streamline their breeding programs so that they can seek new avenues and develop guidelines that will greatly enhance ongoing efforts to establish heat stress tolerance in oilseeds.

Genome-Wide Discovery of miRNAs with Differential Expression Patterns in Responses to Salinity in the Two Contrasting Wheat Cultivars

Salinity is a serious environmental issue. It has a substantial effect on crop yield, as many crop species are sensitive to salinity due to climate change, and it impact is continuing to increase. Plant microRNAs (miRNAs) contribute to salinity stress response in bread wheat. However, the underlying molecular mechanisms by which miRNAs confer salt tolerance in wheat are unclear. We conducted a genome-wide discovery study using Illumina high throughput sequencing and comprehensive in silico analysis to obtain insight into the underlying mechanisms by which small RNAs confer tolerance to salinity in roots of two contrasting wheat cvv., namely Suntop (salt-tolerant) and Sunmate (salt-sensitive). A total of 191 microRNAs were identified in both cultivars, consisting of 110 known miRNAs and 81 novel miRNAs; 181 miRNAs were shared between the two cultivars. The known miRNAs belonged to 35 families consisted of 23 conserved and 12 unique families. Salinity stress induced 43 and 75 miRNAs in Suntop and Sunmate, respectively. Among them, 14 and 29 known and novel miRNAs were expressed in Suntop and 37 and 38 in Sunmate. In silico analysis revealed 861 putative target mRNAs for the 75 known miRNAs and 52 putative target mRNAs for the 15 candidate novel miRNAs. Furthermore, seven miRNAs including tae-miR156, tae-miR160, tae-miR171a-b, tae-miR319, tae-miR159a-b, tae-miR9657 and novel-mir59 that regulate auxin responsive-factor, SPL, SCL6, PCF5, R2R3 MYB, and CBL-CIPK, respectively, were predicted to contribute to salt tolerance in Suntop. This information helps further our understanding of how the molecular mechanisms of salt tolerance are mediated by miRNAs and may facilitate the genetic improvement of wheat cultivars.

Editor's introduction to this issue (G&I 19:1, 2021)

Brassica napus is the third most important oilseed crop in the world; however, in Korea, it is greatly affected by cold stress, limiting seed growth and production. Plants have developed specific stress responses that are generally divided into three categories: cold-stress signaling, transcriptional/post-transcriptional regulation, and stress-response mechanisms. Large numbers of functional and regulatory proteins are involved in these processes when triggered by cold stress. Here, our objective was to investigate the different genetic factors involved in the cold-stress responses of B. napus. Consequently, we treated the Korean B. napus cultivar Naehan at the 4-week stage in cold chambers under different conditions, and RNA and cDNA were obtained. An in silico analysis included 80 cold-responsive genes downloaded from the National Center for Biotechnology Information (NCBI) database. Expression levels were assessed by reverse transcription polymerase chain reaction, and 14 cold-triggered genes were identified under cold-stress conditions. The most significant genes encoded zinc-finger proteins (33.7%), followed by MYB transcription factors (7.5%). In the future, we will select genes appropriate for improving the cold tolerance of B. napus.

Omics: The way forward to enhance abiotic stress tolerance in Brassica napus L

Plant abiotic stresses negative affects growth and development, causing a massive reduction in global agricultural production. Rapeseed (Brassica napus L.) is a major oilseed crop because of its economic value and oilseed production. However, its productivity has been reduced by many environmental adversities. Therefore, it is a prime need to grow rapeseed cultivars, which can withstand numerous abiotic stresses. To understand the various molecular and cellular mechanisms underlying the abiotic stress tolerance and improvement in rapeseed, omics approaches have been extensively employed in recent years. This review summarized the recent advancement in genomics, transcriptomics, proteomics, metabolomics, and their imploration in abiotic stress regulation in rapeseed. Some persisting bottlenecks have been highlighted, demanding proper attention to fully explore the omics tools. Further, the potential prospects of the CRISPR/Cas9 system for genome editing to assist molecular breeding in developing abiotic stress-tolerant rapeseed genotypes have also been explained. In short, the combination of integrated omics, genome editing, and speed breeding can alter rapeseed production worldwide.

Genome-Wide Comparative Analysis of R2R3 MYB Gene Family in Populus and Salix and Identification of Male Flower Bud Development-Related Genes

The MYB transcription factor (TF) family is one of the largest plant transcription factor gene family playing vital roles in plant growth and development, including defense, cell differentiation, secondary metabolism, and responses to biotic and abiotic stresses. As a model tree species of woody plants, in recent years, the identification and functional prediction of certain MYB family members in the poplar genome have been reported. However, to date, the characterization of the gene family in the genome of the poplar's sister species willow has not been done, nor are the differences and similarities between the poplar and willow genomes understood. In this study, we conducted the first genome-wide investigation of the R2R3 MYB subfamily in the willow, identifying 216 R2R3 MYB gene members, and combined with the poplar R2R3 MYB genes, performed the first comparative analysis of R2R3 MYB genes between the poplar and willow. We identified 81 and 86 pairs of R2R3 MYB paralogs in the poplar and willow, respectively. There were 17 pairs of tandem repeat genes in the willow, indicating active duplication of willow R2R3 MYB genes. A further 166 pairs of poplar and willow orthologs were identified by collinear and synonymous analysis. The findings support the duplication of R2R3 MYB genes in the ancestral species, with most of the R2R3 MYB genes being retained during the evolutionary process. The phylogenetic trees of the R2R3 MYB genes of 10 different species were drawn. The functions of the poplar and willow R2R3 MYB genes were predicted using reported functional groupings and clustering by OrthoFinder. Identified 5 subgroups in general expanded in woody species, three subgroups were predicted to be related to lignin synthesis, and we further speculate that the other two subgroups also play a role in wood formation. We analyzed the expression patterns of the GAMYB gene of subgroup 18 (S18) related to pollen development in the male flower buds of poplar and willow at different developmental stages by qRT-PCR. The results showed that the GAMYB gene was specifically expressed in the male flower bud from pollen formation to maturity, and that the expression first increased and then decreased. Both the specificity of tissue expression specificity and conservation indicated that GAMYB played an important role in pollen development in both poplar and willow and was an ideal candidate gene for the analysis of male flower development-related functions of the two species.

Genome-wide identification of R2R3-MYB family in wheat and functional characteristics of the abiotic stress responsive gene TaMYB344

Background

MYB superfamily is one of the most abundant families in plants, which plays important roles in plant growth, development, and productivity. However, to date, researches on MYBs in wheat (Triticum aestivum L.) are scattered mostly, not comprehensive.

Results

In this study, a total of 393 R2R3-MYBs and 12 R1R2R3-MYBs were identified and analyzed including gene structure, chromosomal distribution, synteny relationship, and evolutionary relationship. Then, 29 clusters tandem duplication and 8 clusters segmental duplication genes were discovered. The expression profile of the identified genes under abiotic and biotic stress was analyzed using RNA-seq data. Based on expression patterns analysis, we screened many candidate genes involved in plant response to abiotic and biotic stress. Among them, the functional characteristics of TaMYB344 were further studied. TaMYB344 was localized in the nucleus and functioned as a weak transcriptional activator. We demonstrated that TaMYB344-overexpressing transgenic tobacco plants had enhanced tolerance to drought, heat, and high salt stress.

Conclusions

In this study, 393 R2R3-MYBs and 12 R1R2R3-MYBs in wheat were systemically identified and analyzed. Differential expression analysis indicated that many R2R3-MYBs were involved in abiotic and biotic stress response. We identified a potential candidate gene TaMYB344, overexpression of which in tobacco plants enhanced drought, heat, and salt stress tolerance. These results will provide abundant molecular data for breeding new varieties of wheat in the future.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s12864-020-07175-9.

The R2R3-MYB gene family in banana (Musa acuminata): Genome-wide identification, classification and expression patterns

The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, defense responses and metabolite accumulation. To date MYB family genes have not yet been comprehensively identified in the major staple fruit crop banana. In this study, we present a comprehensive, genome-wide analysis of the MYB genes from Musa acuminata DH-Pahang (A genome). A total of 285 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Organ- and development-specific expression patterns were determined from RNA-Seq data. For 280 M. acuminata MYB genes for which expression was found in at least one of the analysed samples, a variety of expression patterns were detected. The M. acuminata R2R3-MYB genes were functionally categorised, leading to the identification of seven clades containing only M. acuminata R2R3-MYBs. The encoded proteins may have specialised functions that were acquired or expanded in Musa during genome evolution. This functional classification and expression analysis of the MYB gene family in banana establishes a solid foundation for future comprehensive functional analysis of MaMYBs and can be utilized in banana improvement programmes.

MYB43 in Oilseed Rape (Brassica napus) Positively Regulates Vascular Lignification, Plant Morphology and Yield Potential but Negatively Affects Resistance to Sclerotinia sclerotiorum

Arabidopsis thaliana MYB43 (AtMYB43) is suggested to be involved in cell wall lignification. PtrMYB152, the Populus orthologue of AtMYB43, is a transcriptional activator of lignin biosynthesis and vessel wall deposition. In this research, MYB43 genes from Brassica napus (rapeseed) and its parental species B. rapa and B. oleracea were molecularly characterized, which were dominantly expressed in stem and other vascular organs and showed responsiveness to Sclerotinia sclerotiorum infection. The BnMYB43 family was silenced by RNAi, and the transgenic rapeseed lines showed retardation in growth and development with smaller organs, reduced lodging resistance, fewer silique number and lower yield potential. The thickness of the xylem layer decreased by 28%; the numbers of sclerenchymatous cells, vessels, interfascicular fibers, sieve tubes and pith cells in the whole cross section of the stem decreased by 28%, 59%, 48%, 34% and 21% in these lines, respectively. The contents of cellulose and lignin decreased by 17.49% and 16.21% respectively, while the pectin content increased by 71.92% in stems of RNAi lines. When inoculated with S. sclerotiorum, the lesion length was drastically decreased by 52.10% in the stems of transgenic plants compared with WT, implying great increase in disease resistance. Correspondingly, changes in the gene expression patterns of lignin biosynthesis, cellulose biosynthesis, pectin biosynthesis, cell cycle, SA- and JA-signals, and defensive pathways were in accordance with above phenotypic modifications. These results show that BnMYB43, being a growth-defense trade-off participant, positively regulates vascular lignification, plant morphology and yield potential, but negatively affects resistance to S. sclerotiorum. Moreover, this lignification activator influences cell biogenesis of both lignified and non-lignified tissues of the whole vascular organ.

Transcription Factors in Plant Stress Responses: Challenges and Potential for Sugarcane Improvement

Increasing vulnerability of crops to a wide range of abiotic and biotic stresses can have a marked influence on the growth and yield of major crops, especially sugarcane (Saccharum spp.). In response to various stresses, plants have evolved a variety of complex defense systems of signal perception and transduction networks. Transcription factors (TFs) that are activated by different pathways of signal transduction and can directly or indirectly combine with cis-acting elements to modulate the transcription efficiency of target genes, which play key regulators for crop genetic improvement. Over the past decade, significant progresses have been made in deciphering the role of plant TFs as key regulators of environmental responses in particular important cereal crops; however, a limited amount of studies have focused on sugarcane. This review summarizes the potential functions of major TF families, such as WRKY, NAC, MYB and AP2/ERF, in regulating gene expression in the response of plants to abiotic and biotic stresses, which provides important clues for the engineering of stress-tolerant cultivars in sugarcane.

Comparative transcriptome analysis reveals gene expression differences between two peach cultivars under saline-alkaline stress

BACKGROUND

Saline-alkaline stress is a major abiotic stress that is harmful to plant growth worldwide. Two peach cultivars (GF677 and Maotao) display distinct phenotypes under saline-alkaline stress. The molecular mechanism explaining the differences between the two cultivars is still unclear.

RESULTS

In the present study, we systematically analysed the changes in GF677 and Maotao leaves upon saline-alkaline stress by using cytological and biochemical technologies as well as comparative transcriptome analysis. Transmission electron microscopy (TEM) observations showed that the structure of granum was dispersive in Maotao chloroplasts. The biochemical analysis revealed that POD activity and the contents of chlorophyll a and chlorophyll b, as well as iron, were notably decreased in Maotao. Comparative transcriptome analysis detected 881 genes with differential expression (including 294 upregulated and 587 downregulated) under the criteria of |log2 Ratio| ≥ 1 and FDR ≤0.01. Gene ontology (GO) analysis showed that all differentially expressed genes (DEGs) were grouped into 30 groups. MapMan annotation of DEGs showed that photosynthesis, antioxidation, ion metabolism, and WRKY TF were activated in GF677, while cell wall degradation, secondary metabolism, starch degradation, MYB TF, and bHLH TF were activated in Maotao. Several iron and stress-related TFs (ppa024966m, ppa010295m, ppa0271826m, ppa002645m, ppa010846m, ppa009439m, ppa008846m, and ppa007708m) were further discussed from a functional perspective based on the phylogenetic tree integration of other species homologues.

CONCLUSIONS

According to the cytological and molecular differences between the two cultivars, we suggest that the integrity of chloroplast structure and the activation of photosynthesis as well as stress-related genes are crucial for saline-alkaline resistance in GF677. The results presented in this report provide a theoretical basis for cloning saline-alkaline tolerance genes and molecular breeding to improve saline-alkaline tolerance in peach.

Genome-Wide Analysis of Phosphorus Transporter Genes in Brassica and Their Roles in Heavy Metal Stress Tolerance

Phosphorus transporter (PHT) genes encode H2PO4-/H+ co-transporters that absorb and transport inorganic nutrient elements required for plant development and growth and protect plants from heavy metal stress. However, little is known about the roles of PHTs in Brassica compared to Arabidopsis thaliana. In this study, we identified and extensively analyzed 336 PHTs from three diploid (B. rapa, B. oleracea, and B. nigra) and two allotetraploid (B. juncea and B. napus) Brassica species. We categorized the PHTs into five phylogenetic clusters (PHT1-PHT5), including 201 PHT1 homologs, 15 PHT2 homologs, 40 PHT3 homologs, 54 PHT4 homologs, and 26 PHT5 homologs, which are unevenly distributed on the corresponding chromosomes of the five Brassica species. All PHT family genes from Brassica are more closely related to Arabidopsis PHTs in the same vs. other clusters, suggesting they are highly conserved and have similar functions. Duplication and synteny analysis revealed that segmental and tandem duplications led to the expansion of the PHT gene family during the process of polyploidization and that members of this family have undergone purifying selection during evolution based on Ka/Ks values. Finally, we explored the expression profiles of BnaPHT family genes in specific tissues, at various developmental stages, and under heavy metal stress via RNA-seq analysis and qRT-PCR. BnaPHTs that were induced by heavy metal treatment might mediate the response of rapeseed to this important stress. This study represents the first genome-wide analysis of PHT family genes in Brassica species. Our findings improve our understanding of PHT family genes and provide a basis for further studies of BnaPHTs in plant tolerance to heavy metal stress.

Engineering Multiple Abiotic Stress Tolerance in Canola, Brassica napus

Impacts of climate change like global warming, drought, flooding, and other extreme events are posing severe challenges to global crop production. Contribution of Brassica napus towards the oilseed industry makes it an essential component of international trade and agroeconomics. Consequences from increasing occurrences of multiple abiotic stresses on this crop are leading to agroeconomic losses making it vital to endow B. napus crop with an ability to survive and maintain yield when faced with simultaneous exposure to multiple abiotic stresses. For an improved understanding of the stress sensing machinery, there is a need for analyzing regulatory pathways of multiple stress-responsive genes and other regulatory elements such as non-coding RNAs. However, our understanding of these pathways and their interactions in B. napus is far from complete. This review outlines the current knowledge of stress-responsive genes and their role in imparting multiple stress tolerance in B. napus. Analysis of network cross-talk through omics data mining is now making it possible to unravel the underlying complexity required for stress sensing and signaling in plants. Novel biotechnological approaches such as transgene-free genome editing and utilization of nanoparticles as gene delivery tools are also discussed. These can contribute to providing solutions for developing climate change resilient B. napus varieties with reduced regulatory limitations. The potential ability of synthetic biology to engineer and modify networks through fine-tuning of stress regulatory elements for plant responses to stress adaption is also highlighted.

Genome-wide analysis and expression profiles of the StR2R3-MYB transcription factor superfamily in potato (Solanum tuberosum L.)

MYB transcription factors comprise one of the largest families in plant kingdom, which play a variety of functions in plant developmental processes and defence responses, the R2R3-MYB members are the predominant form found in higher plants. In the present study, a total of 111 StR2R3-MYB transcription factors were identified and further phylogenetically classified into 31 subfamilies, as supported by highly conserved gene structures and motifs. Collinearity analysis showed that the segmental duplication events played a crucial role in the expansion of StR2R3-MYB gene family. Synteny analysis indicated that 37 and 13 StR2R3-MYB genes were orthologous to Arabidopsis and wheat (Triticum aestivum), respectively, and these gene pairs have evolved under strong purifying selection. RNA-seq data from different tissues and abiotic stresses revealed tissue-preferential and abiotic stress-responsive StR2R3-MYB genes. We further analyzed StR2R3-MYB genes might be involved in anthocyanin biosynthesis and drought stress by using RNA-seq data of pigmented tetraploid potato cultivars and drought-sensitive and -tolerant tetraploid potato cultivars under drought stress, respectively. Moreover, EAR motifs were found in 21 StR2R3-MYB proteins and 446 pairs of proteins were predicted to interact with 21 EAR motif-containing StR2R3-MYB proteins by constructing the interaction network with medium confidence (0.4). Additionally, Gene Ontology (GO) analysis of the 21 EAR motif-containing StR2R3-MYB proteins was performed to further investigate their functions. This work will facilitate future biologically functional studies of potato StR2R3-MYB transcription factors and enrich the knowledge of MYB superfamily genes in plant species.

Functional activation of a novel R2R3-MYB protein gene, GmMYB68, confers salt-alkali resistance in soybean (Glycine max L.)

Soil salinity significantly reduces soybean (Glycine max L.) production worldwide. Plants resistance to stress conditions is a complex characteristic regulated by multiple signaling pathways. The v-Myb avian myeloblastosis viral oncogene homolog (MYB) transcription factor (TF) plays a crucial role in plant development, secondary metabolism, and abiotic stress responses. GmMYB68-overexpression and RNA interference (RNAi) lines were established for examining the function of G. max GmMYB68 in plant responses to abiotic stresses. The predicted amino acid sequence of GmMYB68 was similar to that of R2R3-MYB proteins. Quantitative real-time PCR analysis revealed that GmMYB68 expression varied in response to abiotic stresses. GmMYB68-overexpression lines showed enhanced resistance to salt and alkali stresses and their osmotic adjustment and photosynthetic rates were also stronger than that of GmMYB68-RNAi and wild type plants. Although wild type and transgenic plants showed no significant differences in agronomic traits under normal conditions, the overexpression of GmMYB68 increased grain number and 100-grain weights under salt stress. Our study identified a valuable TF associated with stress response in soybean, as its overexpression might help improve salt and alkali tolerance in soybean and other crops.

Considerations on the Rational Design of Covalently Conjugated Cell-Penetrating Peptides (CPPs) for Intracellular Delivery of Proteins: A Guide to CPP Selection Using Glucarpidase as the Model Cargo Molecule

Access of proteins to their intracellular targets is limited by a hydrophobic barrier called the cellular membrane. Conjugation with cell-penetrating peptides (CPPs) has been shown to improve protein transduction into the cells. This conjugation can be either covalent or non-covalent, each with its unique pros and cons. The CPP-protein covalent conjugation may result in undesirable structural and functional alterations in the target protein. Therefore, we propose a systematic approach to evaluate different CPPs for covalent conjugations. This guide is presented using the carboxypeptidase G2 (CPG2) enzyme as the target protein. Seventy CPPs -out of 1155- with the highest probability of uptake efficiency were selected. These peptides were then conjugated to the N- or C-terminus of CPG2. Translational efficacy of the conjugates, robustness and thermodynamic properties of the chimera, aggregation possibility, folding rate, backbone flexibility, and aspects of in vivo administration such as protease susceptibility were predicted. The effect of the position of conjugation was evaluated using unpaired t-test (p < 0.05). It was concluded that N-terminal conjugation resulted in higher quality constructs. Seventeen CPP-CPG2/CPG2-CPP constructs were identified as the most promising. Based on this study, the bioinformatics workflow that is presented may be universally applied to any CPP-protein conjugate design.

MYB Gene Family in Potato (Solanum tuberosum L.): Genome-Wide Identification of Hormone-Responsive Reveals Their Potential Functions in Growth and Development

As an important nongrain crop, the growth and yield of potato (Solanum tuberosum L.) is often affected by an unfavorable external environment in the process of cultivation. The MYB family is one of the largest and most important gene families, participating in the regulation of plant growth and development and response to abiotic stresses. Several MYB genes in potato that regulate anthocyanin synthesis and participate in abiotic stress responses have been identified. To identify all Solanum tuberosum L. MYB (StMYB) genes involved in hormone or stress responses to potentially regulate potato growth and development, we identified the MYB gene family at the genome-wide level. In this work, 158 StMYB genes were found in the potato genome. According to the amino acid sequence of the MYB domain and gene structure, the StMYB genes were divided into R2R3-MYB and R1R2R3-MYB families, and the R2R3-MYB family was divided into 20 subgroups (SGs). The expression of 21 StMYB genes from different SGs in roots, stems, leaves, flowers, shoots, stolons, young tubers, and mature tubers was determined by quantitative real-time polymerase chain reaction (qRT-PCR). The expression patterns of StMYB genes in potatoes treated with abscisic acid (ABA), indole-3-acetic acid (IAA), gibberellin acid 3 (GA3), NaCl, mannitol, and heat were also measured. We have identified several potential candidate genes that regulate the synthesis of potato flavonoids or participate in hormone or stress responses. This work provides a comprehensive understanding of the MYB family in potato and will lay a foundation for the future investigation of the potential functions of StMYB genes in the growth and development of potato.

Characterization of the Rosa roxburghii Tratt transcriptome and analysis of MYB genes

Rosa roxburghii Tratt (Rosaceae) has a fruit that is flavorful, economically valuable, and highly nutritious, providing numerous health benefits. Myeloblastosis (MYB) proteins play key roles in the development and fruit quality of R. roxburghii. However, there is little available genomic and transcriptomic information for R. roxburghii. Here, a normalized cDNA library was constructed from five tissues, including the stem, leaf, flower, young fruit, and mature fruit, using the Illumina HiSeq 3000 platform. De novo assembly was performed, and 470.66 million clean reads were obtained. In total, 63,727 unigenes, with an average GC content of 42.08%, were discovered, 60,406 of which were annotated. In addition, 9,354 unigenes were assigned to Gene Ontology categories, and 20,202 unigenes were assigned to 25 Eukaryotic Ortholog Groups. Additionally, 19,508 unigenes were classified into 140 pathways of the Kyoto Encyclopedia of Genes and Genomes database. Based on the transcriptome, 163 unigenes associated with MYBs were detected. Among these genes, 75 genes were significantly expressed in the various tissues, including 10 R1 MYB, 42 R2R3 MYB, one R1R2R3 MYB, three R4 MYB and 19 atypical MYB-like proteins. The expression levels of the 12 MYB genes randomly selected for quantitative real-time PCR analysis corroborated the RNA sequencing results. A total of 37,545 microsatellites were detected, with an average expressed sequence tag–simple sequence repeat frequency of 0.59 (37,545/63,727). This transcriptome data improves our understanding of the role of MYB in R. roxburghii and will be valuable for identifying genes of interest.

A Novel Sweetpotato Transcription Factor Gene IbMYB116 Enhances Drought Tolerance in Transgenic Arabidopsis

Several members of the MYB transcription factor family have been found to regulate growth, developmental processes, metabolism, and biotic and abiotic stress responses in plants. However, the role of MYB116 in plants is still unclear. In this study, a MYB transcription factor gene IbMYB116 was cloned and characterized from the sweetpotato [Ipomoea batatas (L.) Lam.] line Xushu55-2, a line that is considered to be drought resistant. We show here that IbMYB116 is a nuclear protein and that it possesses a transactivation domain at the C terminus. This gene exhibited a high expression level in the leaf tissues of Xushu55-2 and was strongly induced by PEG6000 and methyl-jasmonate (MeJA). The IbMYB116-overexpressing Arabidopsis plants showed significantly enhanced drought tolerance, increased MeJA content, and a decreased H₂O₂ level under drought stress. The overexpression of IbMYB116 in Arabidopsis systematically upregulated jasmonic acid (JA) biosynthesis genes and activated the JA signaling pathway as well as reactive oxygen species (ROS)-scavenging system genes under drought stress conditions. The overall results suggest that the IbMYB116 gene might enhance drought tolerance by activating a ROS-scavenging system through the JA signaling pathway in transgenic Arabidopsis. These findings reveal, for the first time, the crucial role of IbMYB116 in the drought tolerance of plants.

Identification and characterization of a sterically robust phenylalanine ammonia-lyase among 481 natural isoforms through association of in silico and in vitro studies

The enzyme phenylalanine ammonia lyase (PAL) is of special importance for the treatment of phenylketonuria patients. The aim of this study was to find a stable recombinant PAL with suitable kinetic properties among all natural PAL producing species using in silico and experimental approaches. To find such a stable PAL among 481 natural isoforms, 48,000 of 3-D models were predicted using the Modeller 9.10 program and evaluated by Ramachandran plot. Correlation analysis between Ramachandran plot and the energy of different thermodynamic components indicated that this plot could be an appropriate tool to predict protein stability. Hence, PAL6 from Lotus japonicus (LjPAL6) was selected as a stable isoform. Molecular dynamic (MD) simulation for 50 ns and docking has been conducted for LjPAL6-phenylalanine complex. The best PAL-phenylalanine frame was selected by re-docking with l-phenylalanine (L-Phe) and root-mean-square deviation (RMSD) value. MD simulation showed that the complex has a good stability, depicted by the low RMSD value, binding free energy and hydrogen bindings. Docking results showed that LjPAL6 has a high affinity toward l-Phe according to the low level of binding free energy. By overexpressing Ljpal6 in E. coli BL21, a total of 33.5 mg/l of protein was obtained, which has been increased to 83.7 mg/l via the optimization of LjPAL6 production using response surface methodology. The optimal pH and temperature were 8.5 and 50 °C, respectively. LjPAL6 showed a specific activity of 42 nkat/mg protein, with K_m, K_cat and K_cat/K_m values of 0.483 mM, 7 S^-1 and 14.5 S^-1 mM^-1 for l-phe, respectively. In conclusion, finding models with the most reasonable stereo-chemical quality and lowest numbers of steric clashes would result in easier folding. Hence, in silico analyses of bulk data from natural origin will lead one to find an optimal model for in vitro studies and drug design.

Exploring transcription factors reveals crucial members and regulatory networks involved in different abiotic stresses in Brassica napus L

BACKGROUND

Brassica napus (B. napus) encompasses diverse transcription factors (TFs), but thorough identification and characterization of TF families, as well as their transcriptional responsiveness to multifarious stresses are still not clear.

RESULTS

Totally 2167 TFs belonging to five families were genome-widely identified in B. napus, including 518 BnAP2/EREBPs, 252 BnbZIPs, 721 BnMYBs, 398 BnNACs and 278 BnWRKYs, which contained some novel members in comparison with existing results. Sub-genome distributions of BnAP2/EREBPs and BnMYBs indicated that the two families might have suffered from duplication and divergence during evolution. Synteny analysis revealed strong co-linearity between B. napus and its two ancestors, although chromosomal rearrangements have occurred and 85 TFs were lost. About 7.6% and 9.4% TFs of the five families in B. napus were novel genes and conserved genes, which both showed preference on the C sub-genome. RNA-Seq revealed that more than 80% TFs were abiotic stress inducible and 315 crucial differentially expressed genes (DEGs) were screened out. Network analysis revealed that the 315 DEGs are highly co-expressed. The homologous gene network in A. thaliana revealed that a considerable amount of TFs could trigger the differential expression of targeted genes, resulting in a complex clustered network with clusters of genes responsible for targeted stress responsiveness.

CONCLUSIONS

We identified and characterized five TF families in B. napus. Some crucial members and regulatory networks involved in different abiotic stresses have been explored. The investigations deepen our understanding of TFs for stress tolerance in B. napus.

Genome-wide survey, characterization, and expression analysis of RING finger protein genes in Brassica oleracea and their syntenic comparison to Brassica rapa and Arabidopsis thaliana

The ubiquitin-mediated post-translational regulatory pathway regulates a broad range of cell functions in all eukaryotes. It requires the involvement of a large number of E3 ligases, of which more than one third belongs to the RING protein family as in Arabidopsis thaliana. In this study, a total of 756 RING domains in 734 predicted proteins were identified in Brassica oleracea. Their encoding genes were characterized by RING domain type, additional domain, and expression pattern, and mapped on the nine chromosomes of B. oleracea. Comparison of these results with B. rapa and A. thaliana revealed some common as well as species-specific features. Our results showed that the differential gene loss following the whole genome triplication has largely contributed to the RING protein gene number variation among these species, although other factors such as tandem duplication, RING domain loss, or modification had also contributed to this variation. Analysis of RNA-seq data showed that these RING protein genes were functionally diversified and involved in all the stages of plant growth and development, and that the triplicated members were also diverged in expression with one member often more dominantly expressed over the two others in the majority of cases. Our study lays the foundation for further functional determination of each RING protein gene among species of the genus Brassica.

Genome-wide identification and classification of MYB superfamily genes in peach

The MYB transcription factor superfamily is one of the largest superfamilies modulating various biological processes in plants. Over the past few decades, many MYB superfamily genes have been identified and characterized in some plant species. However, genes belonging to the MYB superfamily in peach (Prunus persica) have not been comprehensively identified and characterized although the genome sequences of peach were released several years ago. In total, this study yielded a set of 256 MYB superfamily genes that was divided into five subfamilies: the R2R3-MYB (2R-MYB), R1R2R3-MYB (3R-MYB), MYB-related (1R-MYB), 4R-MYB, and Atypical-MYB subfamilies. These subfamilies contained 128, 4, 109, 1, and 14 members, respectively. The 128 R2R3-MYB subfamily genes in peach were further clustered into 35 groups, and the 109 MYB-related subfamily genes were further clustered into 6 groups: the CCA1-like, CPC-like, TBP-like, I-box-binding-like, R-R-type, and Peach-specific groups. The motif compositions and exon/intron structures within each group within the R2R3-MYB or MYB-related subfamily in peach were highly conserved. The logo sequences of the R2 and R3 repeats of R2R3-MYB subfamily members were highly conserved with those in these repeats of several other plant species. Except for 48 novel peach-specific MYB genes, the remaining 208 out of 256 MYB genes in peach were conserved with the corresponding 198 MYB genes in A. thaliana. Additionally, the 256 MYB genes unevenly distributed on chromosomes 1 to 8 of the peach genome. Eighty-one orthologous pairs of peach/A. thaliana MYB genes were identified among 256 MYB genes in peach and 198 MYB genes in A. thaliana in this study. In addition, 146 pairs of paralogous MYB genes were identified on the eight chromosomes of peach. The expression levels of some of the 51 MYB genes selected for qRT-PCR analysis decreased or increased with red-fleshed fruit development, while the expression patterns of some genes followed no clear rules over the five developmental stages of fruits. This study laid the foundation for further functional analysis of MYB superfamily genes in peach and enriched the knowledge of MYB superfamily genes in plant species.

Genome-wide analysis of basic helix-loop-helix superfamily members in peach

The basic helix-loop-helix (bHLH) transcription factor family is a superfamily found in all eukaryotes that plays important roles in regulating growth and development. Over the past several decades, many bHLH superfamily genes have been identified and characterized in herbaceous and woody plants. However, the genes belonging to the bHLH superfamily in peach (Prunus persica) have not yet been comprehensively identified and characterized. Here, we identified 95 members of the bHLH superfamily in the peach genome, and these genes were classified into 19 subfamilies based on a phylogenetic comparison with bHLH proteins from Arabidopsis. The members within each subfamily were highly conserved according to the analysis of motif compositions and exon/intron organizations. The 95 bHLH genes were unevenly distributed on chromosomes 1 to 8 of the peach genome. We identified 57 pairs of bHLH members that were orthologous between peach and Arabidopsis. Additionally, 48 pairs of paralogous bHLH genes were identified on the eight chromosomes of the peach genome. Coupled with relative expression analysis of bHLH genes in red-fleshed peach fruit at five developmental stages, we identified several bHLH genes that might be involved in fruit development and anthocyanin biosynthesis. This study provides insight into the molecular mechanisms through which these genes are involved in the regulation of biological and biochemical processes in peach and lays the foundation for further studies on these genes.