Xiaojiao Han

Integration of small RNAs, degradome and transcriptome sequencing in hyperaccumulator Sedum alfredii uncovers a complex regulatory network and provides insights into cadmium phytoremediation

Summary The hyperaccumulating ecotype of Sedum alfredii Hance is a cadmium (Cd)/zinc/lead co‐hyperaccumulating species of Crassulaceae. It is a promising phytoremediation candidate accumulating substantial heavy metal ions without obvious signs of poisoning. However, few studies have focused on the regulatory roles of miRNAs and their targets in the hyperaccumulating ecotype of S. alfredii. Here, we combined analyses of the transcriptomics, sRNAs and the degradome to generate a comprehensive resource focused on identifying key regulatory miRNA‐target circuits under Cd stress. A total of 87 721 unigenes and 356 miRNAs were identified by deep sequencing, and 79 miRNAs were differentially expressed under Cd stress. Furthermore, 754 target genes of 194 miRNAs were validated by degradome sequencing. A gene ontology (GO) enrichment analysis of differential miRNA targets revealed that auxin, redox‐related secondary metabolism and metal transport pathways responded to Cd stress. An integrated analysis uncovered 39 pairs of miRNA targets that displayed negatively correlated expression profiles. Ten miRNA‐target pairs also exhibited negative correlations according to a real‐time quantitative PCR analysis. Moreover, a coexpression regulatory network was constructed based on profiles of differentially expressed genes. Two hub genes, ARF4 (auxin response factor 4) and AAP3 (amino acid permease 3), which might play central roles in the regulation of Cd‐responsive genes, were uncovered. These results suggest that comprehensive analyses of the transcriptomics, sRNAs and the degradome provided a useful platform for investigating Cd hyperaccumulation in S. alfredii, and may provide new insights into the genetic engineering of phytoremediation.

Selection and validation of reference genes for real-time quantitative PCR in hyperaccumulating ecotype of Sedum alfredii under different heavy metals stresses.

Real-time Quantitative PCR (RT-qPCR) has become an effective method for accurate analysis of gene expression in several biological systems as well as under different experimental conditions. Although with high sensitivity, specificity and broad dynamic range, this method requires suitable reference genes for transcript normalization in order to guarantee reproducible and meaningful results. In the present study, we evaluated five traditional housekeeping genes and five novel reference genes in Hyperaccumulating ecotype of Sedum alfredii, a well known hyperaccumulator for heavy metals phytoremediation, under Cd, Pb, Zn and Cu stresses of seven different durations. The expression stability of these ten candidates were determined with three programs - geNorm, NormFinder and BestKeeper. The results showed that all the selected reference genes except for SAND could be used for RT-qPCR normalization. Among them UBC9 and TUB were ranked as the most stable candidates across all samples by three programs together. For the least stable reference genes, however, BestKeeper produced different results compared with geNorm and NormFinder. Meanwhile, the expression profiles of PCS under Cd, Pb, Zn and Cu stresses were assessed using UBC9 and TUB respectively, and similar trends were obtained from the results of the two groups. The distinct expression patterns of PCS indicated that various strategies could be taken by plants in adaption to different heavy metals stresses. This study will provide appropriate reference genes for further gene expression quantification using RT-qPCR in Hyperaccumulator S. alfredii.

Selection of suitable reference genes for quantitative real-time PCR gene expression analysis in Salix matsudana under different abiotic stresses

Salix matsudana is a deciduous, rapidly growing willow species commonly cultivated in China, which can tolerate drought, salt, and heavy metal stress conditions. Selection of suitable reference genes for quantitative real-time PCR is important for normalizing the expression of the key genes associated with various stresses. To validate suitable reference genes, we selected 11 candidate reference genes (five traditional housekeeping genes and six novel genes) and analyzed their expression stability in various samples, including different tissues and under different abiotic stress treatments. The expression of these genes was determined using five programs—geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder. The results showed that α-TUB2 (alpha-tubulin 2) and DnaJ (chaperone protein DnaJ 49) were the most stable reference genes across all the tested samples. We measured the expression profiles of the defense response gene SmCAT (catalase) using the two most stable and one least stable reference genes in all samples of S. matsudana. The relative quantification of SmCAT varied greatly according to the different reference genes. We propose that α-TUB2 and DnaJ should be the preferred reference genes for normalization and quantification of transcript levels in future gene expression studies in willow species under various abiotic stress conditions.

Validation of reference genes aiming accurate normalization of qRT-PCR data in Dendrocalamus latiflorus Munro.

Background Dendrocalamus latiflorus Munro distributes widely in subtropical areas and plays vital roles as valuable natural resources. The transcriptome sequencing for D. latiflorus Munro has been performed and numerous genes especially those predicted to be unique to D. latiflorus Munro were revealed. qRT-PCR has become a feasible approach to uncover gene expression profiling, and the accuracy and reliability of the results obtained depends upon the proper selection of stable reference genes for accurate normalization. Therefore, a set of suitable internal controls should be validated for D. latiflorus Munro. Results In this report, twelve candidate reference genes were selected and the assessment of gene expression stability was performed in ten tissue samples and four leaf samples from seedlings and anther-regenerated plants of different ploidy. The PCR amplification efficiency was estimated, and the candidate genes were ranked according to their expression stability using three software packages: geNorm, NormFinder and Bestkeeper. GAPDH and EF1α were characterized to be the most stable genes among different tissues or in all the sample pools, while CYP showed low expression stability. RPL3 had the optimal performance among four leaf samples. The application of verified reference genes was illustrated by analyzing ferritin and laccase expression profiles among different experimental sets. The analysis revealed the biological variation in ferritin and laccase transcript expression among the tissues studied and the individual plants. Conclusions geNorm, NormFinder, and BestKeeper analyses recommended different suitable reference gene(s) for normalization according to the experimental sets. GAPDH and EF1α had the highest expression stability across different tissues and RPL3 for the other sample set. This study emphasizes the importance of validating superior reference genes for qRT-PCR analysis to accurately normalize gene expression of D. latiflorus Munro.

Overexpressing the Sedum alfredii Cu/Zn Superoxide Dismutase Increased Resistance to Oxidative Stress in Transgenic Arabidopsis

Superoxide dismutase (SOD) is a very important reactive oxygen species (ROS)-scavenging enzyme. In this study, the functions of a Cu/Zn SOD gene (SaCu/Zn SOD), from Sedum alfredii, a cadmium (Cd)/zinc/lead co-hyperaccumulator of the Crassulaceae, was characterized. The expression of SaCu/Zn SOD was induced by Cd stress. Compared with wild-type (WT) plants, overexpression of SaCu/Zn SOD gene in transgenic Arabidopsis plants enhanced the antioxidative defense capacity, including SOD and peroxidase activities. Additionally, it reduced the damage associated with the overproduction of hydrogen peroxide (H2O2) and superoxide radicals (O2•-). The influence of Cd stress on ion flux across the root surface showed that overexpressing SaCu/Zn SOD in transgenic Arabidopsis plants has greater Cd uptake capacity existed in roots. A co-expression network based on microarray data showed possible oxidative regulation in Arabidopsis after Cd-induced oxidative stress, suggesting that SaCu/Zn SOD may participate in this network and enhance ROS-scavenging capability under Cd stress. Taken together, these results suggest that overexpressing SaCu/Zn SOD increased oxidative stress resistance in transgenic Arabidopsis and provide useful information for understanding the role of SaCu/Zn SOD in response to abiotic stress.

Transcriptome analysis of immature xylem in the Chinese fir at different developmental phases

Background.Chinese fir [Cunninghamia lanceolata (Lamb.) Hook.] is one of the most important native tree species for timber production in southern China. An understanding of overall fast growing stage, stem growth stage and senescence stage cambium transcriptome variation is lacking. We used transcriptome sequencing to identify the repertoire of genes expressed during development of xylem tissue in Chinese fir, aiming to delineate the molecular mechanisms of wood formation. Results. We carried out transcriptome sequencing at three different cultivation ages (7Y, 15Y and 21Y) generating 68.71 million reads (13.88 Gbp). A total of 140,486 unigenes with a mean size of 568.64 base pairs (bp) were obtained via de novo assembly. Of these, 27,427 unigenes (19.52%) were further annotated by comparison to public protein databases. A total of 5,331 (3.79%) unigenes were mapped into 118 pathways by searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Differentially expressed genes (DEG) analysis identified 3, 16 and 5,899 DEGs from the comparison of 7Y vs. 15Y, 7Y vs. 21Y and 15Y vs. 21Y, respectively, in the immature xylem tissues, including 2,638 significantly up-regulated and 3,280 significantly down-regulated genes. Besides, five NAC transcription factors, 190 MYB transcription factors, and 34 WRKY transcription factors were identified respectively from Chinese fir transcriptome. Conclusion. Our results revealed the active transcriptional pathways and identified the DEGs at different cultivation phases of Chinese fir wood formation. This transcriptome dataset will aid in understanding and carrying out future studies on the molecular basis of Chinese fir wood formation and contribute to future artificial production and applications.

Isolation of Salt Stress-Related Genes from Aspergillus glaucus CCHA by Random Overexpression in Escherichia coli

The halotolerant fungus Aspergillus glaucus CCHA was isolated from the surface of wild vegetation around a saltern with the salinity range being 0–31%. Here, a full-length cDNA library of A. glaucus under salt stress was constructed to identify genes related to salt tolerance, and one hundred clones were randomly selected for sequencing and bioinformatics analysis. Among these, 82 putative sequences were functionally annotated as being involved in signal transduction, osmolyte synthesis and transport, or regulation of transcription. Subsequently, the cDNA library was transformed into E. coli cells to screen for putative salt stress-related clones. Five putative positive clones were obtained from E. coli cells grown on LB agar containing 1 M NaCl, on which they showed rapid growth compared to the empty vector control line. Analysis of transgenic Arabidopsis thaliana lines overexpressing CCHA-2142 demonstrated that the gene conferred increased salt tolerance to plants as well by protecting the cellular membranes, suppressing the inhibition of chlorophyll biosynthesis. These results highlight the utility of this A. glaucus cDNA library as a tool for isolating and characterizing genes related to salt tolerance. Furthermore, the identified genes can be used for the study of the underlying biology of halotolerance.

Identification and expression analysis of salt-responsive genes using a comparative microarray approach in Salix matsudana

Salt stress exerts negative effects on plant growth, development and yields, with roots being the primary site of both perception and damage. Salix matsudana (Chinese willow) is tolerant of high salinity. However, genes associated with this trait were rarely characterized. Therefore, we first performed salt-stress treatment on S. matsudana plants, then identified differentially expressed genes by comparison of salt-treated roots and untreated controls using microarray analysis. A total of 403 salt-responsive genes were identified, of which 239 were repressed and 164 were up-regulated. Functional classification analysis revealed that these genes belonged to families encoding proteins involved in metabolism, regulation of transcription, signal transduction, hormone responses, abiotic stress responses, and other processes related to growth and development. This suggested that when S. matsudana was confronted with salt stress, coordinated adjustments are made to physiological and biochemical processes, which would then allow more resources to be allocated to protective mechanisms to avoid salt injury. The expression patterns of representative genes were further validated and the diversity of the temporal profiles indicated that a combination of several genes and the initiation of diverse pathways performed functions in S. matsudana salt tolerance. This work represents the first study employing microarrays to investigate salt tolerance in S. matsudana. The data presented herein enhance our understanding of the molecular mechanisms of S. matsudana responses to salinity stress and lay the groundwork for genetic engineering strategies to improve stress tolerance of agronomically important species.

Sedum alfredii SaNramp6 Metal Transporter Contributes to Cadmium Accumulation in Transgenic Arabidopsis thaliana

The plant natural resistance-associated macrophage protein (Nramp) family plays an important role in tolerance to heavy metal stress. However, few Nramps have been functionally characterized in the heavy metal-accumulating plant Sedum alfredii. Here, Nramp6 was cloned and identified from S. alfredii and its function analyzed in transgenic Arabidopsis thaliana. SaNramp6 cDNA contains an open reading frame of 1, 638 bp encoding 545 amino acids. SaNramp6′s expression can be induced by cadmium (Cd) stress, and, after treatment, it peaked at one week and 12 h in the roots and leaves, respectively. SaNramp6 localized to the plasma membrane in protoplasts isolated from A. thaliana, Nicotiana benthamiana lower leaf and onion (Allium cepa) epidermal cells. The heterologous expression of SaNramp6 in the Δycf1 yeast mutant increased the Cd content in yeast cells. SaNramp6 also rescued the low Cd accumulation of the A. thaliana nramp1 mutant. Transgenic A. thaliana expressing SaNramp6 exhibited high Cd accumulation levels, as determined by a statistical analysis of the Cd concentration, translocation factors and net Cd2+ fluxes under Cd stress. Thus, SaNramp6 may play a significant role in improving Cd accumulation, and the gene may be useful for the biotechnological development of transgenic plants for phytoremediation.

Functional Characterization of a Gene in Sedum alfredii Hance Resembling Rubber Elongation Factor Endowed with Functions Associated with Cadmium Tolerance.

Cadmium is a major toxic heavy-metal pollutant considering their bioaccumulation potential and persistence in the environment. The hyperaccumulating ecotype of Sedum alfredii Hance is a Zn/Cd co-hyperaccumulator inhabiting in a region of China with soils rich in Pb/Zn. Investigations into the underlying molecular regulatory mechanisms of Cd tolerance are of substantial interest. Here, library screening for genes related to cadmium tolerance identified a gene resembling the rubber elongation factor gene designated as SaREFl. The heterologous expression of SaREFl rescued the growth of a transformed Cd-sensitive strain (ycf1). Furthermore, SaREFl-expressing Arabidopsis plants were more tolerant to cadmium stress compared with wild type by measuring parameters of root length, fresh weight and physiological indexes. When under four different heavy metal treatments, we found that SaREFl responded most strongly to Cd and the root was the plant organ most sensitive to this heavy metal. Yeast two-hybrid screening of SaREFl as a bait led to the identification of five possible interacting targets in Sedum alfredii Hance. Among them, a gene annotated as prenylated Rab acceptor 1 (PRA1) domain protein was detected with a high frequency. Moreover, subcellular localization of SaREF1-GFP fusion protein revealed some patchy spots in cytosol suggesting potential association with organelles for its cellular functions. Our findings would further enrich the connotation of REF-like genes and provide theoretical assistance for the application in breeding heavy metal-tolerant plants.