The real time quantitative reverse transcription PCR (qRT-PCR) is becoming increasingly important to gain insight into function of genes. Given the increased sensitivity, ease and reproducibility of qRT-PCR, the requirement of suitable reference genes for normalization has become important and stringent. It is now known that the expression of internal control genes in living organism vary considerably during developmental stages and under different experimental conditions. For economically important Brassica crops, only a couple of reference genes are reported till date. In this study, expression stability of 12 candidate reference genes including ACT2, ELFA, GAPDH, TUA, UBQ9 (traditional housekeeping genes), ACP, CAC, SNF, TIPS-41, TMD, TSB and ZNF (new candidate reference genes), in a diverse set of 49 tissue samples representing different developmental stages, stress and hormone treated conditions and cultivars of Brassica juncea has been validated. For the normalization of vegetative stages the ELFA, ACT2, CAC and TIPS-41 combination would be appropriate whereas TIPS-41 along with CAC would be suitable for normalization of reproductive stages. A combination of GAPDH, TUA, TIPS-41 and CAC were identified as the most suitable reference genes for total developmental stages. In various stress and hormone treated samples, UBQ9 and TIPS-41 had the most stable expression. Across five cultivars of B. juncea, the expression of CAC and TIPS-41 did not vary significantly and were identified as the most stably expressed reference genes. This study provides comprehensive information that the new reference genes selected herein performed better than the traditional housekeeping genes. The selection of most suitable reference genes depends on the experimental conditions, and is tissue and cultivar-specific. Further, to attain accuracy in the results more than one reference genes are necessary for normalization.

The globally cultivated Brassica species possess diverse aliphatic glucosinolates, which are important for plant defense and animal nutrition. The committed step in the side chain elongation of methionine-derived aliphatic glucosinolates is catalyzed by methylthioalkylmalate synthase, which likely evolved from the isopropylmalate synthases of leucine biosynthesis. However, the molecular basis for the evolution of methylthioalkylmalate synthase and its generation of natural product diversity in Brassica is poorly understood. Here, we show that Brassica genomes encode multiple methylthioalkylmalate synthases that have differences in expression profiles and 2-oxo substrate preferences, which account for the diversity of aliphatic glucosinolates across Brassica accessions. Analysis of the 2.1 Å resolution x-ray crystal structure of Brassica juncea methylthioalkylmalate synthase identified key active site residues responsible for controlling the specificity for different 2-oxo substrates and the determinants of side chain length in aliphatic glucosinolates. Overall, these results provide the evolutionary and biochemical foundation for the diversification of glucosinolate profiles across globally cultivated Brassica species, which could be used with ongoing breeding strategies toward the manipulation of beneficial glucosinolate compounds for animal health and plant protection.

SummaryBrassica juncea (Indian mustard), a globally important oilseed crop, contains relatively high amount of seed glucosinolates ranging from 80 to 120 lmol/g dry weight (DW). One of the major breeding objectives in oilseed Brassicas is to improve the seed-meal quality through the development of low-seed-glucosinolate lines (<30 lmol/g DW), as high amounts of certain seed glucosinolates are known to be anti-nutritional and reduce the meal palatability. Here, we report the development of transgenic B. juncea lines having seed glucosinolates as low as 11.26 lmol/g DW, through RNAi-based targeted suppression of BjMYB28, a R2R3-MYB transcription factor family gene involved in aliphatic glucosinolate biosynthesis. Targeted silencing of BjMYB28 homologs provided significant reduction in the anti-nutritional aliphatic glucosinolates fractions, without altering the desirable nonaliphatic glucosinolate pool, both in leaves and seeds of transgenic plants. Molecular characterization of single-copy, low glucosinolate homozygous lines confirmed significant down-regulation of BjMYB28 homologs vis-a-vis enhanced accumulation of BjMYB28-specific siRNA pool. Consequently, these low glucosinolate lines also showed significant suppression of genes involved in aliphatic glucosinolate biosynthesis. The low glucosinolate trait was stable in subsequent generations of the transgenic lines with no visible off-target effects on plant growth and development. Various seed quality parameters including fatty acid composition, oil content, protein content and seed weight of the low glucosinolate lines also remained unaltered, when tested under containment conditions in the field. Our results indicate that targeted silencing of a key glucosinolate transcriptional regulator MYB28 has huge potential for reducing the glucosinolates content and improving the seed-meal quality of oilseed Brassica crops.

Glucosinolates are Capparales-specific secondary metabolites that have immense potential in human health and agriculture. Unlike Arabidopsis thaliana, our knowledge about glucosinolate regulators in the Brassica crops is sparse. In the current study, four MYB28 homologues were identified (BjuMYB28-1,-2,-3,-4) from the polyploid Brassica juncea, and the effects of allopolyploidization on the divergence of gene sequence, structure, function, and expression were assessed. The deduced protein sequences of the four BjuMYB28 genes showed 76.1–83.1% identity with the Arabidopsis MYB28. Phylogenetic analysis revealed that the four BjuMYB28 proteins have evolved via the hybridization and duplication processes forming the B. juncea genome (AABB) from B. rapa (AA) and B. nigra (BB), while retaining high levels of sequence conservation. Mutant complementation and over-expression studies in A. thaliana showed that all four BjuMYB28 genes encode functional MYB28 proteins and resulted in similar aliphatic glucosinolate composition and content. Detailed expression analysis using qRT-PCR assays and promoter-GUS lines revealed that the BjuMYB28 genes have both tissue- and cell-specific expression partitioning in B. juncea. The two B-genome origin BjuMYB28 genes had more abundant transcripts during the early stages of plant development than the A-genome origin genes. However, with the onset of the reproductive phase, expression levels of all four BjuMYB28 increased significantly, which may be necessary for producing and maintaining high amounts of aliphatic glucosinolates during the later stages of plant development. Taken together, our results suggest that the four MYB28 genes are differentially expressed and regulated in B. juncea to play discrete though overlapping roles in controlling aliphatic glucosinolate biosynthesis.

Glucosinolates are amino acids derived secondary metabolites, invariably present in Brassicales, which have huge health and agricultural benefits. Sulphoraphane, the breakdown product of glucosinolate glucoraphanin is known to posses anti-cancer properties. AOP (2-oxoglutarate-dependent dioxygenases) or GSL-ALK enzyme catalyzes the conversion of desirable glucoraphanin to deleterious gluconapin and progoitrin, which are present in very high amounts in most of the cultivable Brassica species including Brassica juncea. In this study we showed that B. juncea encodes four functional homologs of GSL-ALK gene and constitutive silencing of GSL-ALK homologs resulted in accumulation of glucoraphanin up to 43.11 μmoles g−1 DW in the seeds with a concomitant reduction in the anti-nutritional glucosinolates. Glucoraphanin content was found remarkably high in leaves as well as sprouts of the transgenic lines. Transcript quantification of high glucoraphanin lines confirmed significant down-regulation of GSL-ALK homologs. Growth and other seed quality parameters of the transgenic lines did not show drastic difference, compared to the untransformed control. High glucoraphanin lines also showed higher resistance towards stem rot pathogen Sclerotinia sclerotiorum. Our results suggest that metabolic engineering of GSL-ALK has huge potential for enriching glucoraphanin content, and improve the oil quality and vegetable value of Brassica crops.

No abstract

14-3-3s are highly conserved, multigene family proteins that have been implicated in modulating various biological processes. The presence of inherent polyploidy and genome complexity has limited the identification and characterization of 14-3-3 proteins from globally important Brassica crops. Through data mining of Brassica rapa, the model Brassica genome, we identified 21 members encoding 14-3-3 proteins namely, BraA.GRF14.a to BraA.GRF14.u. Phylogenetic analysis indicated that B. rapa contains both ε (epsilon) and non-ε 14-3-3 isoforms, having distinct intron-exon structural organization patterns. The non-ε isoforms showed lower divergence rate (Ks < 0.45) compared to ε protein isoforms (Ks > 0.48), suggesting class-specific divergence pattern. Synteny analysis revealed that mesohexaploid B. rapa genome has retained 1–5 orthologs of each Arabidopsis 14-3-3 gene, interspersed across its three fragmented sub-genomes. qRT-PCR analysis showed that 14 of the 21 BraA.GRF14 were expressed, wherein a higher abundance of non-ε transcripts was observed compared to the ε genes, indicating class-specific transcriptional bias. The BraA.GRF14 genes showed distinct expression pattern during plant developmental stages and in response to abiotic stress, phytohormone treatments, and nutrient deprivation conditions. Together, the distinct expression pattern and differential regulation of BraA.GRF14 genes indicated the occurrence of functional divergence of B. rapa 14-3-3 proteins during plant development and stress responses.

GTR1 and GTR2 transporters are components of the source to sink translocation network of glucosinolates, which are major defence metabolites in the Brassicaceae. These transporters can be genetically manipulated for reduction of seed-glucosinolates without inhibiting glucosinolate biosynthesis, thereby maintaining the inherent defence potential of plants. However, the different roles of GTRs in influencing tissue-specific distribution of glucosinolates in agriculturally important Brassica crops are yet unknown. Here, we report functional characterization of two groups of glucosinolate transporters (GTR1 and GTR2) from Brassica juncea based on gene expression data, biochemical analysis, gene-complementation studies in GTR-deficient mutants and RNAi-based knockdown followed by insect feeding experiments. Although both GTRs showed ubiquitous expression patterns and broad substrate specificity, the single-gene knockdown lines displayed different phenotypes. The GTR2-knockdown plants showed a significant reduction of glucosinolates in seeds and a higher accumulation in leaves and pods, while the GTR1-knockdown plants displayed a smaller reduction of glucosinolates in seeds and significantly lower glucosinolate levels in leaves. Consequently, knockdown of GTR2 resulted in higher resistance towards the generalist pest, Spodoptera litura. Overall, our study highlights the distinctive roles of B. juncea GTRs in tissue-specific accumulation of glucosinolates and the potential for manipulating GTR2 for enhanced nutrition and plant defence.