Annotation and characterization of Cd-responsive metal transporter genes in rapeseed (Brassica napus).

In higher plants, heavy metal transporters are responsible for metal uptake, translocation and homeostasis. These metals include essential metals such as zinc (Zn) or manganese (Mn) and non-essential metals like cadmium (Cd) or lead (Pb). Although a few heavy metal transporters have been well identified in model plants (e.g. Arabidopsis and rice), little is known about their functionality in rapeseed (Brassica napus). B. napus is an important oil crop ranking the third largest sources of vegetable oil over the world. Importantly, B. napus has long been considered as a desirable candidate for phytoremediation owning to its massive dry weight productivity and moderate to high Cd accumulation. In this study, 270 metal transporter genes (MTGs) from B. napus genome were identified and annotated using bioinformatics and high-throughput sequencing. Most of the MTGs (74.8%, 202/270) were validated by RNA-sequencing (RNA-seq) the seedling libraries. Based on the sequence identity, nine superfamilies including YSL, OPT, NRAMP, COPT, ZIP, CDF/MTP, HMA, MRP and PDR have been classified. RNA-sequencing profiled 202 non-redundant MTGs from B. napus seedlings, of which, 108 MTGs were differentially expressed and 62 genes were significantly induced under Cd stress. These differentially expressed genes (DEGs) are dispersed in the rapeseed genome. Some of the genes were well confirmed by qRT-PCR. Analysis of the genomic distribution of MTGs on B. napus chromosomes revealed that their evolutional expansion was probably through localized allele duplications.

Identification of Cd-responsive RNA helicase genes and expression of a putative BnRH 24 mediated by miR158 in canola (Brassica napus).

RNA helicases play crucial roles in RNA splicing, transport, editing and degradation, protein translation initiation and siRNA-mediated gene silencing. However, knowledge about their functionality in rapeseed (Brassica napus) is rare. In the study, we identified and annotated 271 RNA helicase genes from B. napus using bioinformatics and high-throughput RNA-sequencing (RNA-seq). Three subfamilies DEAD-box, DEAH-box, or DExD/H-box have been identified. One hundred and ninety-five RNA helicases were confirmed by RNA-seq and 49 were identified to differentially respond to cadmium (Cd) stress (> 1.5 fold change, p < 0.05). As an example, we functionally specified BnaA04g26450D encoding a BnRH24 under Cd exposure. BnRH24 is a constitutive gene expressing throughout the life span. Using our previously generated degradome datasets, we found that BnRH24 can be cleaved by miR158, suggesting that BnRH24 is a target of miR158 in B. napus. The mature miR158 was induced, while BnRH24 was repressed in B. napus under Cd stress. The contrasting expression pattern of B. napus miR158 and BnRH24 under the normal and Cd would support the post-transcriptional regulation of BnRH24 by miR158. Ectopic expression of BnRH24 in Arabidopsis revealed that the transgenic lines showed more sensitivity to Cd toxicity by reducing root elongation, fresh mass production, chlorophyll accumulation and increasing oxidative products such as O2-., H2O2 and thiobarbituric acid reactive substances (TBARS), indicating that the controlling the level of BnRH24 by miR158 may be required for Cd tolerance in plants.

Genome-wide identification of Cd-responsive NRAMP transporter genes and analyzing expression of NRAMP 1 mediated by miR167 in Brassica napus.

In plants, metal transporters are responsible for metal uptake, translocation and homeostasis. These metals include essential nutrients such as zinc (Zn) and manganese (Mn) or non-essential metals like cadmium (Cd) and lead (Pb). Although a few metal transporters have been well characterized in model plants, little is known about their functionality in rapeseed (Brassica napus). In the study, 22 NRAMP transporter genes from B. napus genome were identified and annotated using bioinformatics and high-throughput RNA-sequencing (RNA-seq). Based on the sequence identity, these NRAMP transporters can be classified into 6 subfamilies. RNA-seq analysis revealed that 19 NRAMP transporters were detected and some of the genes were well confirmed by qRT-PCR. Ten NRAMP transporters (45.5%, 10/22) were found to be differentially expressed (> 2 fold change, p < 0.05) under Cd exposure. As an example, we specified expression of BnNRAMP1b under Cd exposure. BnNRAMP1b is a constitutive gene expressing throughout all development stages including seedlings, vegetative tissue, flowers and siliques. Expression of BnNRAMP1b can be strongly induced in seedlings exposed to 80, 160 and 240 µM Cd. To define whether BnNRAMP1b was specific for Cd transport, a yeast (wild-type, BY4741) system with its mutants (ycf1, zrc1, and smf1) defective in transport activity of Cd, Zn and Mn, respectively were tested. Compared to empty vectors (pYES2), cells carrying BnNRAMP1b can rescue the transport functions. As a consequence, excess Cd, Zn and Mn were taken in the cells, which led to metal toxicity, suggesting that BnNRAMP1b is responsible for transport of these metals in B. napus. Using our previously created degradome datasets, we found that BnNRAMP1b could be cleaved by miR167, suggesting that BnNRAMP1b is a target of miR167 in B. napus. The contrasting expression pattern of BnNRAMP1b and miR167 under Cd stress supported the post-transcriptional regulation of BnNRAMP1b by miR167.

Variation of DNA methylation patterns associated with gene expression in rice (Oryza sativa) exposed to cadmium.

We report genome-wide single-base resolution maps of methylated cytosines and transcriptome change in Cd-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between Cd-exposed and Cd-free rice genomes. There are 2320 non-redundant differentially methylated regions detected in the genome. RNA sequencing revealed 2092 DNA methylation-modified genes differentially expressed under Cd exposure. More genes were found hypermethylated than those hypomethylated in CG, CHH and CHG (where H is A, C or T) contexts in upstream, gene body and downstream regions. Many of the genes were involved in stress response, metal transport and transcription factors. Most of the DNA methylation-modified genes were transcriptionally altered under Cd stress. A subset of loss of function mutants defective in DNA methylation and histone modification activities was used to identify transcript abundance of selected genes. Compared with wide type, mutation of MET1 and DRM2 resulted in general lower transcript levels of the genes under Cd stress. Transcripts of OsIRO2, OsPR1b and Os09g02214 in drm2 were significantly reduced. A commonly used DNA methylation inhibitor 5-azacytidine was employed to investigate whether DNA demethylation affected physiological consequences. 5-azacytidine provision decreased general DNA methylation levels of selected genes, but promoted growth of rice seedlings and Cd accumulation in rice plant.

DEAD-box RNA helicase 6 regulates drought and abscisic acid stress responses in rapeseed (Brassica napus).

DEAD-box RNA helicase is a major subfamily of RNA helicases with vital roles played in plant growth, development, and plant-environment interactions. RNA helicase 6 in rapeseed (Brassica napus) (BnRH6) is a member of DEAD-box RNA helicase. While previous research has demonstrated the role of BnRH6 in salt stress regulation, the involvement of BnRH6 in drought stress adaptation remains unknown. This report described a function of BnRH6 in drought stress response. BnRH6 was sufficiently induced by osmotic stress. Transgenic Brassica napus and Arabidopsis thaliana (Arabidopsis) overexpressing BnRH6 (OE) showed a drought tolerance phenotype, characterized by improved plant growth, increased survival rates, reduced water loss, leaf chlorosis and oxidative stress. Furthermore, BnRH6 was also induced by exogenous abscisic acid (ABA). BnRH6 overexpression plants exhibited ABA hypersensitivity with lagging seed germination, growth stunt and diminished stomatal opening in the presence of ABA, suggesting the involvement of ABA signal. Assessment of several well-identified drought stress responsive genes such as Calcium-dependent Protein Kinase 14 (BnCDPK14), Enhanced Response to ABA1 (BnERA1) and ABA Insensitive 1 (BnABI1) revealed that their expressions were accordingly changed in BnOE plants, possibly with interplays between BnRH6 and those genes. Our data highlighted the functional role of BnRH6, which plays a positive role in active regulation of drought stress response likely through an ABA-dependent manner in rapeseed plants.

Identification of genomic ATP binding cassette (ABC) transporter genes and Cd-responsive ABCs in Brassica napus.

The plant ATP binding cassette (ABC) transporters are one of the integral membrane proteins responsible for uptake and allocation of a wide range of metabolites and xenobiotics including heavy metals (e.g. zinc, manganese and cadmium). They play multiple roles in plant growth, development and environmental stress responses. Although the ABC transporters have been identified in model plants such as Arabidopsis and rice, they are have not been annotated and identified in rapeseed (Brassica napus) and also, little is known about functionality of these metal transporters. B. napus is an important oil crop ranking the third largest source of vegetable oil worldwide. Importantly, it is long considered as a desirable candidate for phytoremediation owning to its massive dry weight productivity and moderate Cd accumulation. In this study, we identified 314 ABC protein genes from B. napus using bioinformatics and high-throughput sequencing. Eight subfamilies including ABCA-G and ABCI have been categorized. The ABCG proteins constitute the largest subfamily with 116 members, and the ABCB and ABCC subfamilies ranks second and third with 69 and 47 members, respectively. Analyses of ABCs in B. napus genome reveal that their evolutional expansion was through localized allele duplications. Most of the ABC genes (74.2%, 233/314) were validated by RNA-sequencing rapeseed seedlings. Among the 233 profiled BnaABCs, 132 genes were differentially expressed (>1.5 fold change, p < 0.05) and 84 genes were significantly induced under Cd stress. Analyses of specific cis-elements in the upstream of eight representative genes show diverse motifs, which potentially respond to environmental stress, hormone responsiveness and other development signals.