Genes underlying environmental adaptability tend to be over-retained in polyploid plant species. Zinc deficiency (ZnD) and iron deficiency (FeD), excess Zn (ZnE) and cadmium exposure (CdE) are major environmental problems for crop cultivation, but little is known about the differential expression of duplicated genes upon these stress conditions. Applying Tag-Seq technology to leaves of Brassica rapa grown under FeD, ZnD, ZnE or CdE conditions, with normal conditions as a control, we examined global gene expression changes and compared the expression patterns of multiple paralogs. We identified 812, 543, 331 and 447 differentially expressed genes under FeD, ZnD, ZnE and CdE conditions, respectively, in B.rapa leaves. Genes involved in regulatory networks centered on the transcription factors bHLH038 or bHLH100 were differentially expressed under (ZnE-induced) FeD. Further analysis revealed that genes associated with Zn, Fe and Cd responses tended to be over-retained in the B.rapa genome. Most of these multiple-copy genes showed the same direction of expression change under stress conditions. We conclude that the duplicated genes involved in trace element responses in B.rapa are functionally redundant, making the regulatory network more complex in B.rapa than in Arabidopsis thaliana.