Main conclusion Genome-wide identification, classification, expression analyses, and functional characterization of GRAS genes in oil crop, Brassica napus, indicate their importance in root development and stress response.

MYB proteins are involved in diverse important biological processes in plants. Herein, we obtained the MYB superfamily from the allotetraploid Brassica napus, which contains 227 MYB-related (BnMYBR/Bn1R-MYB), 429 R2R3-MYB (Bn2R-MYB), 22 R1R2R3-MYB (Bn3R-MYB), and two R1R2R2R1/2-MYB (Bn4R-MYB) genes. Phylogenetic analysis classified the Bn2R-MYBs into 43 subfamilies, and the BnMYBRs into five subfamilies. Sequence characteristics and exon/intron structures within each subfamily of the Bn2R-MYBs and BnMYBRs were highly conserved. The whole superfamily was unevenly distributed on 19 chromosomes and underwent unbalanced expansion in B. napus. Allopolyploidy between B. oleracea and B. rapa mainly contributed to the expansion in their descendent B. napus, in which B. rapa-derived genes were more retained. Comparative phylogenetic analysis of 2R-MYB proteins from nine Brassicaceae and seven non-Brassicaceae species identified five Brassicaceae-specific subfamilies and five subfamilies that are lacking from the examined Brassicaceae species, which provided an example for the adaptive evolution of the 2R-MYB gene family alongside angiosperm diversification. Ectopic expression of four Bn2R-MYBs under the control of the viral CaMV35S and/or native promoters could rescue the lesser root hair phenotype of the Arabidopsis thaliana wer mutant plants, proving the conserved negative roles of the 2R-MYBs of the S15 subfamily in root hair development. RNA-sequencing data revealed that the Bn2R-MYBs and BnMYBRs had diverse transcript profiles in roots in response to the treatments with various hormones. Our findings provide valuable information for further functional characterizations of B. napus MYB genes.

MADS-box transcription factors are important for plant growth and development, and hundreds of MADS-box genes have been functionally characterized in plants. However, less is known about the functions of these genes in the economically important allopolyploid oil crop, Brassica napus. We identified 307 potential MADS-box genes (BnMADSs) in the B. napus genome and categorized them into type I (Mα, Mβ, and Mγ) and type II (MADS DNA-binding domain, intervening domain, keratin-like domain, and C-terminal domain [MIKC]c and MIKC*) based on phylogeny, protein motif structure, and exon-intron organization. We identified one conserved intron pattern in the MADS-box domain and seven conserved intron patterns in the K-box domain of the MIKCc genes that were previously ignored and may be associated with function. Chromosome distribution and synteny analysis revealed that hybridization between Brassica rapa and Brassica oleracea, segmental duplication, and homologous exchange (HE) in B. napus were the main BnMADSs expansion mechanisms. Promoter cis-element analyses indicated that BnMADSs may respond to various stressors (drought, heat, hormones) and light. Expression analyses showed that homologous genes in a given subfamily or sister pair are highly conserved, indicating widespread functional conservation and redundancy. Analyses of BnMADSs provide a basis for understanding their functional roles in plant development.

Auxin response factor (ARF) is a member of the plant-specific B3 DNA binding superfamily. Here, we report the results of a comprehensive analysis of ARF genes in allotetraploid Brassica napus (2n = 38, AACC). Sixty-seven ARF genes were identified in B . napus ( BnARFs ) and divided into four subfamilies (I–IV). Sixty-one BnARFs were distributed on all chromosomes except C02; the remaining were on Ann and Cnn. The full length of the BnARF proteins was highly conserved especially within each subfamily with all members sharing the N-terminal DNA binding domain (DBD) and the middle region (MR), and most contained the C-terminal dimerization domain (PBI). Twenty-one members had a glutamine-rich MR that may be an activator and the remaining were repressors. Accordingly, the intron patterns are highly conserved in each subfamily or clade, especially in DBD and PBI domains. Several members in subfamily III are potential targets for miR167. Many putative cis-elements involved in phytohormones, light signaling responses, and biotic and abiotic stress were identified in BnARF promoters, implying their possible roles. Most ARF proteins are likely to interact with auxin/indole-3-acetic acid (Aux/IAA) -related proteins, and members from different subfamilies generally shared many common interaction proteins. Whole genome-wide duplication (WGD) by hybridization between Brassica rapa and Brassica oleracea and segmental duplication led to gene expansion. Gene loss following WGD is biased with the A n -subgenome retaining more ancestral genes than the C n -subgenome. BnARFs have wide expression profiles across vegetative and reproductive organs during different developmental stages. No obvious expression bias was observed between A n - and C n -subgenomes. Most synteny-pair genes had similar expression patterns, indicating their functional redundancy. BnARFs were sensitive to exogenous IAA and 6-BA treatments especially subfamily III. The present study provides insights into the distribution, phylogeny, and evolution of ARF gene family.

Background: The basic helix-loop-helix (bHLH) gene family is one of the largest transcription factor families in plants and is functionally characterized in diverse species. However, less is known about its functions in the economically important allopolyploid oil crop, Brassica napus. Results: We identified 602 potential bHLHs in the B. napus genome (BnabHLHs) and categorized them into 35 subfamilies, including seven newly separated subfamilies, based on phylogeny, protein structure, and exon-intron organization analysis. The intron insertion patterns of this gene family were analyzed and a total of eight types were identified in the bHLH regions of BnabHLHs. Chromosome distribution and synteny analyses revealed that hybridization between Brassica rapa and Brassica oleracea was the main expansion mechanism for BnabHLHs. Expression analyses showed that BnabHLHs were widely in different plant tissues and formed seven main patterns, suggesting they may participate in various aspects of B. napus development. Furthermore, when roots were treated with five different hormones (IAA, auxin; GA3, gibberellin; 6-BA, cytokinin; ABA, abscisic acid and ACC, ethylene), the expression profiles of BnabHLHs changed significantly, with many showing increased expression. The induction of five candidate BnabHLHs was confirmed following the five hormone treatments via qRT-PCR. Up to 246 BnabHLHs from nine subfamilies were predicted to have potential roles relating to root development through the joint analysis of their expression profiles and homolog function. Conclusion: The 602 BnabHLHs identified from B. napus were classified into 35 subfamilies, and those members from the same subfamily generally had similar sequence motifs. Overall, we found that BnabHLHs may be widely involved in root development in B. napus. Moreover, this study provides important insights into the potential functions of the BnabHLHs super gene family and thus will be useful in future gene function research.

Background: The basic helix-loop-helix (bHLH) gene family is one of the largest transcription factor families in plants and is functionally characterized in diverse species. However, less is known about its functions in the economically important allopolyploid oil crop, Brassica napus . Results : We identified 602 potential bHLHs in the B. napus genome ( BnabHLHs ) and categorized them into 35 subfamilies, including seven newly separated subfamilies, based on phylogeny, protein structure, and exon-intron organization analysis. The intron insertion patterns of this gene family were analyzed and a total of

No abstract

Nitrogen (N) is an essential factor for crop yield. Here, we characterized 605 genes from 25 gene families that form the complex gene networks of N utilization pathway in Brassica napus. We found unequal gene distribution between the An- and Cn-sub-genomes, and that genes derived from Brassica rapa were more retained. Transcriptome analysis indicated that N utilization pathway gene activity shifted in a spatio-temporal manner in B. napus. A low N (LN) stress RNA-seq of B. napus seedling leaves and roots was generated, which proved that most N utilization related genes were sensitive to LN stress, thereby forming co-expression network modules. Nine candidate genes in N utilization pathway were confirmed to be significantly induced under N deficiency conditions in B. napus roots, indicating their potential roles in LN stress response process. Analyses of 22 representative species confirmed that the N utilization gene networks were widely present in plants ranging from Chlorophyta to angiosperms with a rapid expansion trend. Consistent with B. napus, the genes in this pathway commonly showed a wide and conserved expression profile in response to N stress in other plants. The network, genes, and gene-regulatory modules identified here represent resources that may enhance the N utilization efficiency or the LN tolerance of B. napus.