Design of a tobacco exon array with application to investigate the differential cadmium accumulation property in two tobacco varieties
- Florian Martin†1,
- Lucien Bovet†1,
- Audrey Cordier1,
- Mario Stanke2,
- Irfan Gunduz3,
- Manuel C Peitsch1 and
- Nikolai V Ivanov1Email author
© Martin et al; licensee BioMed Central Ltd. 2012
Received: 26 August 2012
Accepted: 23 November 2012
Published: 28 November 2012
For decades the tobacco plant has served as a model organism in plant biology to answer fundamental biological questions in the areas of plant development, physiology, and genetics. Due to the lack of sufficient coverage of genomic sequences, however, none of the expressed sequence tag (EST)-based chips developed to date cover gene expression from the whole genome. The availability of Tobacco Genome Initiative (TGI) sequences provides a useful resource to build a whole genome exon array, even if the assembled sequences are highly fragmented. Here, the design of a Tobacco Exon Array is reported and an application to improve the understanding of genes regulated by cadmium (Cd) in tobacco is described.
From the analysis and annotation of the 1,271,256 Nicotiana tabacum fasta and quality files from methyl filtered genomic survey sequences (GSS) obtained from the TGI and ~56,000 ESTs available in public databases, an exon array with 272,342 probesets was designed (four probes per exon) and tested on two selected tobacco varieties.
Two tobacco varieties out of 45 accumulating low and high cadmium in leaf were identified based on the GGE biplot analysis, which is analysis of the genotype main effect (G) plus analysis of the genotype by environment interaction (GE) of eight field trials (four fields over two years) showing reproducibility across the trials. The selected varieties were grown under greenhouse conditions in two different soils and subjected to exon array analyses using root and leaf tissues to understand the genetic make-up of the Cd accumulation.
An Affymetrix Exon Array was developed to cover a large (~90%) proportion of the tobacco gene space. The Tobacco Exon Array will be available for research use through Affymetrix array catalogue. As a proof of the exon array usability, we have demonstrated that the Tobacco Exon Array is a valuable tool for studying Cd accumulation in tobacco leaves. Data from field and greenhouse experiments supported by gene expression studies strongly suggested that the difference in leaf Cd accumulation between the two specific tobacco cultivars is dependent solely on genetic factors and genetic variability rather than on the environment.
KeywordsTobacco plant Nicotiana tabacum Exon array Gene expression Cadmium accumulation
Tobacco (Nicotiana tabacum L.) is a species in the large family of the Solanaceae and is important as an agronomic crop, since more than six million tons of tobacco are produced per year throughout the world. For decades the tobacco plant has served as a model organism in plant biology and has helped to answer fundamental biological questions in the areas of plant development, physiology, and genetics. It has made its scientific reputation due to the ease of gene transformation and genetic manipulation. Tobacco is an allopolyploid species (2n=4x=48) and shares its basic chromosome number of x = 12 with many other Solanaceae species, such as tomato, potato, pepper, and eggplant. It is most likely the result of a tetraploidization event [1, 2] involving Nicotiana sylvestris (S-genome) and a species closely related to modern day Nicotiana tomentosiformis (T-genome). Considering both genomes together, tobacco is at the high end of genome sizes (4.5 Gbp) in the Solanaceae  and contains a large proportion of repetitive sequences [4, 5].
The North Carolina State University (NCSU) Tobacco Genome Initiative (TGI) was started in 2002 in cooperation with Philip Morris USA to gather genetic information on N. tabacum by means of sequencing gene-rich regions of genomic DNA and cDNA libraries of Hicks Broadleaf variety. The TGI website (http://www.pngg.org/tgi/) contains related project information and links for data download. TGI has been leveraged to build a dense tobacco genetic map  and microsatellite marker kits for variety identification . Many phenotypes of importance in tobacco, such as heavy metal accumulation, nutrient deficiency, and yield, are thought to be transcriptionally controlled. Therefore, microarray technology is a suitable tool to study genetic variation and environmental effects with the objective to improve varieties of crops.
Expressed sequence tag (EST)-based tobacco microarrays have been used successfully in tobacco plant research [8, 9]. Cui et al. compared the gene expression of 2,831 selected tobacco genes between the trichomes and the leaves with removed trichomes. Trichomes predominantly expressed genes involved in the second metabolic processes, defence responses, and metabolism regulation . A Tobacco Expression Atlas (TobEA)  was constructed through systematic measurements of gene expression across different tobacco samples. To achieve this, a custom-built Affymetrix GeneChip™ was designed from the cDNA sequences originating from multiple tissues (seeds, roots, leaves, flowers, etc.) of several tobacco varieties. However, due to the lack of sufficient coverage of genomic sequences before the TGI data were released, none of these arrays were intended to cover genome-wide gene expression.
There is evidence that in many plant genomes, including tobacco, the repetitive regions are heavily methylated and the gene regions are undermethylated . Consequently, a methyl filtration approach was used in the TGI to reduce the complexity of tobacco genome sequencing and at the same time to cover a large portion of the coding regions in the genome.
The challenge was to design a whole genome functional chip that can be used to provide a more consistent and comprehensive picture than an EST-based microarray. Because the assembly of unmethylated tobacco genome sequences yields highly fragmented contigs covering a large set of functional genes representing gene-rich portion (~13%) of the genome, this exon array is a feasible alternative to a full genome gene array. In this study we present the application of the Tobacco Exon Array to measure differential gene expression with the objective to improve our understanding of cadmium (Cd) accumulation by the tobacco plant.
Cd accumulation in crop plants such as tobacco can lead to human exposure to this carcinogenic metal. Therefore, there is a considerable interest to find strategies to produce plants with low Cd content. Tobacco tends to sequester higher concentrations of Cd in the leaves than in the roots. The Cd concentration in tobacco leaves usually ranges from 0.1 to 5 μg/g dry weight; the Cd concentration in lower leaves is usually more elevated compared to higher and mid/upper-stalk position leaves . In addition, large differences in leaf Cd content exist among tobacco varieties and between growing regions , which suggests the predominant role of both genetic background and environmental conditions.
In this study, to better understand the genotype-environment interactions, eight field trials over a period of two crop years were performed to identify two variety candidates that differentially accumulate Cd in leaves, independently of the environment. The Tobacco Exon Array was then used to unravel the difference in transcription between those two varieties.
Design of tobacco exon array
Due to the large size and complexity of the tobacco genome, a methyl filtration approach [10, 13, 14] was used within the TGI to enrich the genomic clone library with gene-rich genomic DNA sequences . At the time of the study (June 2007), 1,271,256 N. tabacum fasta and quality files from methyl filtered genomic survey sequences (GSS) were obtained from the TGI. The cleaning procedure with seqclean removed 13,010 reads matching known vectors, including the ones used by TGI, and 146 reads originating from E. coli contamination. Further quality control scripts removed 23,883 mitochondrial and plastid reads, 245 short reads, and 257 reads containing undetermined bases, resulting in 1,233,715 reads. Simple repeats were identified and masked in 114,897 (9.0%) reads. Complex repeats were identified and masked in 156,899 (12.7%) reads. In the absence of tobacco plant-specific repeats, repeat libraries containing retrotransposons and DNA transposons of Brassica, maize, rice, barley, wheat, and Solanum family from mips-REdat  and SOL Genomics Network  were used. The final set contained 1,172,176 reads adjusted for cleaning and trimming. Due to lack of sufficient genome coverage (<1X) of the data (even over the euchromatin part of the genome estimated to be ~500 Mbp based on the tomato genome), the tobacco genome could only be partially assembled. Celera assembler produced 387,927 singletons and 183,198 contigs with the N50 = 990 (contigs only) and 4.28 reads per contig, resulting in ~0.13X coverage of the tobacco genome.
As described in the Methods section, the evidence from various sources for exons (ab initio predictions, transcript sequences, homology within related species, conservation at genome level) and against exons (repeats) was integrated towards a set of ~525,000 consensus exonic regions and passed as hints to AUGUSTUS to predict exact exon boundaries. Each exon was assigned a score based on the number of supporting evidence sources. The exons were further filtered to remove those with a lower degree of evidence, resulting in 374,907 exon candidates.
At the same time, about 56,000 contigs of ESTs were used to identify transcript fragments not already included in the genomic exonic regions. In order to determine the strand, an open reading frame (ORF) of at least 150 bp was required. EST contigs that appeared to be redundant with the exonic regions from the genome were removed. The resulting 312,053 exonic regions were scored by supporting evidence (0–7) so they could later be prioritized after further filtering by Affymetrix and after determining the exact size of the exon array. Finally, 272,342 exons were included to make 1,089,368 probes (four probes per exon) selected by Affymetrix.
Field trial results
Cadmium concentrations for V21 (flue-cured) and V5 (burley) grown in field trials
P-value 2 (Corrected data)
To perform transcriptomic analysis, V21 and V5 were also grown under controlled conditions in the greenhouse to identify putative gene candidates involved in Cd accumulation, since low variability were found in the field.
Greenhouse experiment with V21 and V5
Cadmium concentrations in leaf and leaf weight for V21 (flue-cured) and V5 (burley) grown in the greenhouse n two different soils
Leaf Weight (g)
Quality control (QC) for the tobacco exon array
Gene expression analysis (set of tools to analyze exon array)
Once the chips meet the quality standard of the QC metrics, the probesets with high DABG P-values were filtered out (threshold was set to 0.1). The experimental design led to the definition of a linear model (see Methods) that was fitted for each probeset using LIMMA  R package. For each coefficient of the linear model, P-values were adjusted by Benjamini-Hochberg false discovery rate (fdr) correction . Finally the probesets of interest were selected based on coefficients and fdr P-value thresholds (0.05). Because many of probesets had high coefficients, we also put a threshold on the coefficient values (3 on the log2-scale, corresponding to an 8-fold change) to be able to interpret the results.
The highest differentially expressed exons (absolute coefficient values above 5) are shown in the Additional file 1: Tables S1 and S2. For these transcripts, the number of differentially expressed exons in leaf is actually well-distributed between V5 and V21; 29 exons are highly expressed in V21 and 32 in V5. This suggests that the two varieties have each a specific set of expressed exons in leaf contributing to the known physiological differences observed between flue-cured tobacco (V21) and burley tobacco (V5), i.e., color of the stems and leaf thickness. However, such a gene expression pattern may contribute to the differential Cd sequestration profiles within the leaves of the two varieties as shown in Tables 1 and 2. Interestingly, the only exon-related gene that could play a specific role in Cd accumulation corresponds to a MATE transporter as identified by the probe NtPMIa1g52548e2_s_st (5’-CAGTACGACAATTCTAGGGTGGGTATTCATGATTTCTCTTGGCTTCAATGCAGCAGCAAG-3’). MATE transporters have already been reported to transport metals such as aluminium in barley and Arabidopsis [20, 21], and to be regulated upon Cd stress  or involved in Cd detoxification and tolerance . In order to support the biological interpretation of the experiment, gene set enrichment analysis was performed for each contrast (e.g., Variety effect and Soil effect) using biological function-associated gene sets as a priori knowledge (see Methods). Thirty-six significant gene sets (fdr<0.01) were identified for the variety effect in leaves (see Additional file 1: Table S3). Among them, ABC transporters, known to be involved in Cd transport in plants [24–26] may play a role in the differential Cd accumulation and/or sequestration in leaves of both species. Besides ABC transporters, six other gene sets are related to transcription factors and seven gene sets to other type of transporters, both may contribute to leaf Cd sequestration. Other generic gene sets linked to light photosynthesis reactions and sugar pathways are certainly more specific markers for burley and flue-cured tobacco.
In roots, the difference between V5 and V21 was more apparent: 31 exons were highly expressed in V21, but only 14 exons in V5. In this set of genes, no probeset annotations were found to be directly linked to Cd uptake, sequestration, or translocation to the shoot. In V21, more genes likely contributing to higher active root Cd sequestration were activated than in V5. In addition, only 5 significant gene sets (fdr<0.01) were found in roots (root variety effect, see Additional file 1: Table S3), suggesting that only a few specific root functions are different between V5 and V21. In conclusion, microarray analysis confirmed that the Cd accumulation trait is linked to a constitutive differential transcriptional gene expression of multiple genes explained by general genetic variability between V5 and V21 and not due to soil properties.
The analysis of the set of root genes responding to the soil effect shows that 132 gene sets (see Additional file 1: Table S3) play a role by activating or deactivating several root functions, depending on the soil composition. Among these gene sets, the major gene families involved in Cd transport, sequestration, root-to-shoot translocation, and Cd cell responses are represented, namely “ABC transporters” [27, 28], “metal”, “cell wall general”, “cell wall lignification” , and “callus formation” , “enzyme ascorbate glutathione cycle” , “stress_GST”, “formation of GST-complexes” , ”Ca2+ cation antiporter” , “Multidrug Resistance exporter” (vacuolar Cd sequestration via MRP transporter) ; “Pleiotropic Drug Resistance protein” [26, 35], and “P-type ATPase” (Cd root-to-shoot translocation) , and “Cd sequestration” . According to the literature, these sets of genes cover the main pathway driving Cd from root uptake to leaf storage. Interestingly, although the transcription of these gene sets is affected by the soil environment, no major differences in Cd were found in the leaf between field trials and greenhouse experiments (compare relative decrease between Table 1 and Table 2). This suggests that plants adapt to the soil environmental conditions to maintain fitness and optimal growth conditions. The resulting effect is visualized by Cd content which is not affected or only slightly affected by the environmental conditions.
Altogether, these data indicate that only a few Cd-induced genes are involved in the differential accumulation of Cd between the flue-cured variety V21 and the Burley variety V5. Thus, a global genetic variability between these two tobacco species likely exists, including differential constitutive expression of certain genes. Nevertheless, we cannot exclude the possibility that Cd variation may also result from post-translational regulation involving specific proteins in one or the other variety.
Semi-quantitative RT-PCR validation
While the whole tobacco genome sequence is still not completed, the TGI has provided the scientific community with valuable sequence information that has created opportunities to develop plant biology tools, such as the Tobacco Exon Array. The Tobacco Exon Array was built to include a large set of experimentally confirmed and predicted genes within the tobacco gene space reaching an estimated coverage of more than 90%. Many applications of exon arrays have been developed for other plant species . Our experience with the Tobacco Exon Array demonstrated that it is a powerful tool, as shown in the present study on Cd accumulation.
Across all internal experiments, only a few chips were discarded after examination of the QC metrics. The DABG calls were used to eliminate probesets prior to statistical analysis and few expression data were contradicted by RT-PCR validation. It should be stressed, however, that due to the fragmented nature of the draft genome of N. tabacum, the Tobacco Exon Array was analyzed as a gene chip disregarding alternative splicing events, even though this type of analysis can be performed for well annotated multi-exon genes [42–44].
The experimental goal was to find suitable tobacco varieties exhibiting reduced Cd translocation from root to shoot and ultimately to discover genes involved in leaf Cd accumulation. The field experiments allowed us to highlight two tobacco varieties differentially accumulating Cd that were stable across the environments (avoiding the GxE effect), thereby enabling greenhouse experiments to identify gene transcripts involved in Cd accumulation using Tobacco Exon Array.
An Affymetrix Tobacco Exon Array was developed based on the current genome and EST sequence data from the TGI to cover a large proportion of the tobacco gene space. The advantages of the exon array design include (i) representation of the genes not yet found in the currently available EST libraries, (ii) ability to investigate alternative splicing in the tobacco plant, and (iii) equal probe coverage of each exon of the gene. The Tobacco Exon Array will be available for research use through Affymetrix array catalogue. The experimental data described in this work have been released to Gene Expression Omnibus (GEO) under accession number GSE42319 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE42319) and will also be publicly available through GeneVestigator (https://www.genevestigator.com/gv/plant.jsp) or SGN (http://solgenomics.net) web sites.
As a proof of the exon array usability, we have demonstrated in this study that the Tobacco Exon Array is a valuable tool for studying stresses in tobacco, in particular Cd accumulation in leaves. Data from field and greenhouse experiments supported by gene expression studies strongly suggest that the difference in leaf Cd accumulation between the two specific tobacco cultivars is dependent solely on genetic factors and genetic variability.
Tobacco genomic read assembly and gene prediction
Fasta and quality files generated using phred (Phil Green, http://www.phrap.org) from N. tabacum methyl filtered GSS reads were provided by North Carolina State University as a part of the TGI . A cleaning script seqclean (http://sourceforge.net/projects/seqclean/) with “-A -L” options was run to remove reads matching either to known vectors included in the UniVec collection (http://www.ncbi.nlm.nih.gov/VecScreen/UniVec.html) or to contamination with E. coli K12 genome (NC_000913). Reads corresponding to N. tabacum mitochondrial genome (NC_006581) or to N. tabacum plastid genome (NC_001879), reads shorter than 100 bases, and reads with more than 3% of undetermined bases were also removed. Lower-case masking of simple and complex repeats was performed using Tandem Repeat Finder  and RepeatMasker , respectively. The library of complex repeats consisted of mips-REdat  and repeats from SOL Genomics Network . Sixty bases at each 5’- and 3’-end were trimmed to make sure that they did not contain lower-case masked sequences. Tobacco genome assembly was performed on the set of cleaned and trimmed reads using Celera assembler (http://sourceforge.net/projects/wgs-assembler/?source=directory) with optimized parameters.
Genes were predicted in the genomic contigs with two gene prediction programs: AUGUSTUS  and FGENESH . AUGUSTUS was trained on a set of bona fide gene structures of known tobacco genes available in the GenBank. The tobacco-trained version of AUGUSTUS included models for the untranslated regions (5’ and 3’, spliced and unspliced) and the promoter region of a gene. Gene prediction with AUGUSTUS was performed using “hints” based on evidence from tobacco ESTs and cDNAs, homology at the protein level, genomic conservation with four plants (tomato, Arabidopsis thaliana, Medicago trancatula and rice), and cDNA from the tomato and potato. FGENESH was trained for tobacco by SoftBerry Inc. and run as an ab initio gene finder with default parameters. Further iterations of the tobacco genome assembly were also masked with a new N. tabacum repeat library constructed ab initio by running RepeatScout  on the assembled genome.
In order to measure genomic conservation, a 5-way alignment of the genomes of tobacco, tomato, A. thaliana, Medicago truncatula, and rice was performed using the multiple sequence aligning programs TBA  and BLASTZ . A conservation landscape (see Figure 2) was computed using the phylogenetic Hidden Markov Model phastcons . Tobacco ESTs were obtained from several libraries as a part of the TGI.
Tobacco exon array design
The tobacco ESTs were assembled into 55,520 contigs, 8,451 of which were filtered out by the criteria for transposable elements and repeats. An additional 26,882 EST contigs were filtered because they matched a consensus exon candidate from genomic DNA (min 98% identity); in the case of duplicate representation we gave precedence to exon candidates represented in the genomic assembly. In the remaining 20,187 EST contigs we searched for the longest ORF and kept the 16,258 EST contigs with an ORF of minimum length 150. For these, we inferred a putative strand from the ORF.
Exon candidates from genomic DNA and from EST contigs were joined and the genomic candidates were cleaned for redundancies (98% threshold). This resulted in a set of 312,053 exon candidates, 12,925 of which were represented by ESTs, but were not included in the genome assembly.
Field trials and GGE analysis
Varieties were attributed to subplots using a completely randomized block design. The experimental unit was a pooled sample of 20 plants. The objective of the field trials was to study the variety (genetic make-up) effect on Cd uptake. Forty-five varieties were grown in four environments (two countries, Poland and the Philippines, over two years). In each country, two fields in which burley and flue-cured cultivation practices were used.
Fields were designed by completely randomized block design, with six Blocks and 45 Rows. Plants in each block were sampled at medium stalk position and a pooled sample was created for each block and analyzed by Inductively Coupled Plasma Mass Spectrometry (ICP-MS) for Cd content in leaves. The data were modeled as Cd~E+G+GxE+e, with E being a random effect (G=Genotype, E=Environment=Field x Year). Differences in soil composition and irrigation are known to impact Cd bioavailability creating a confounding with the genotype effect under study, therefore, prior to the GGE analysis , the data were corrected by Generalized Additive Models (GAM)  for the spatial effect within the field (X~1+Variety+s(Rows,Columns)+e, where s is a cubic regression spline) and then scaled by Environment. Modeling the spatial effect showed its importance in the majority of the fields (see Additional file 2: Figure S1 for an example).
Plant material and growth (V21 vs. V5 Experiment)
Tobacco seeds of V5 and V21 varieties were germinated in floating trays and grown for three weeks before being transferred to 5 liter pots. Two soils were prepared, one a peaty soil (minus Cd soil, referred to as Soil 1) and the other a soil mixture containing half of the same peaty soil and half of an agricultural clay soil from Poland (pH = 5.1) containing more bioavailable Cd (plus Cd soil, referred to as Soil 2). To ensure that there were no toxic Cd effects and to allow the monitoring of differential Cd accumulation in both tobacco varieties, the soils did not exceed 0.12 mg/kg total Cd. In each soil, six plants were planted and placed on a table according to a randomized pot design. Plants were watered daily on an automatic drip with a full fertilizer (93.03 mg N l-1, 49.94 mg P2O5 l-1, 172.27 mg K2O l-1, 13.67 mg Mg l-1, 62.03 mg Ca l-1, 27.62 mg S l-1 and microelements: 0.472 mg Fe-EDTA l-1, 0.309 mg Mn l-1, 0.147 mg Zn l-1, 0.122 mg B l-1, 0.027 mg Cu l-1, 0.027 mg Mo l-1). The illumination cycle was 14/10 day/night with 24°C during the day and 20°C during the night. After nine weeks of growth and just before flowering, a pool of leaf disks at leaf stalk positions 5, 6 and, 7 and one lateral root were collected simultaneously from each plant. Full leaves 5, 6, and 7 were then collected and dried at 65°C in an oven for three days.
Preparation of RNA samples
Approximately 100 μg of plant tissue was placed in a 1.5 ml Eppendorf tube (1/3 in volume), and snap frozen in liquid nitrogen. RNA extraction was performed by two methods: Trizol® for microarray analysis (see below) and Qiagen RNeasy Plant mini kit for semi-quantitative RT-PCR.
In case of RNA extraction for microarray experiment, 400 μl of Trizol® reagent (Invitrogen) was added to the frozen plant tissue and the sample was ground to homogeneity before an additional 600 μl of Trizol® was added. Samples were vortexed for 15 s, then incubated at for 5 min room temperature. 200 μl of chloroform was added and the tubes were gently mixed before centrifugation (12,000 × g) at 4°C for 15 min. The aqueous phase was transferred to an Eppendorf tube containing 500 μl of isopropyl alcohol. After 15 min incubation on ice, RNA was pelleted by centrifugation at 4°C for 10 min. The RNA pellet was finally washed with 1 ml of 70% ethanol, and centrifuged for five minutes a 7,500 × g. After air-drying, RNA was resuspended in 50 μl of nucleotide-free H2O. RNA concentration was analyzed by measuring optical density (OD) at 260 nm. RNA quality control was determined by OD260/OD280 and using Agilent 2100 Bioanalyzer before cDNA preparation.
Frozen samples (−80°C) packed in dry ice were sent to DNA-Vision (Charleroi, Belgium). Total RNA isolation using a Trizol® method (Invitrogen 155596–018), microarray hybridization and quality checked by Agilent 2100 Bioanalyzer were executed at DNA-Vision. Affymetrix hybridizations were performed using Affymetrix kits with catalog numbers 900652 and 900454; probe labeling was checked as suggested by the manufacturer.
Microarray data analysis
An in-house QC pipeline was developed to assess the quality of the gene-expression data. In addition to the standard quality metrics suggested by Affymetrix http://media.affymetrix.com/support/technical/whitepapers/exon_gene_arrays_qa_whitepaper.pdf, an in-house QC metrics including probe-level models, Normalized Unscaled Standard Error (NUSE) and Relative Log Expression (RLE) plots , and the analysis of DABG results http://media.affymetrix.com/support/technical/whitepapers/exon_background_correction_whitepaper.pdf were used. As the exon array design had no mismatch probes, summarization was performed using Robust Multiarray Average (RMA) method . A total of 272,342 probeset expression values were generated, and DABG P-values were computed to assess the significance of the signal obtained for each probeset. This involved the background probes that are spread over the chip. These random probes have a varying GC content http://media.affymetrix.com/support/technical/whitepapers/exon_background_correction_whitepaper.pdf. The QC pipeline involves a combination of Affymetrix Power Tools (APT) http://www.affymetrix.com/partners_programs/programs/developer/tools/powertools.affx and Bioconductor packages, for which the Tobacco Exon Array (TobArray520623F) cdf environment was created. Once the expression values were available, differential gene expression analysis was performed using moderated t-statistics in linear model LIMMA .
In addition, gene sets were defined by first annotating probesets by homology to A. thaliana genes and using A. thaliana gene sets [http://raetschlab.org//suppl/kirmes/a.thaliana-geneset-FASTAs.zip/view.html], and mean-rank gene-set enrichment analysis was performed.
Cd and Zn chemical analysis
Cd and Zn metal analyses were performed by ALS Group Czech Republic (Prague, http://www.alsglobal.com/). Prior to the metal determination, the samples were homogenized and then digested in organic matrices according to standard operating procedure (SOP) from ALS Group. Cd and Zn concentrations were determined by Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) or ICP-MS according to the manufacturer’s instructions.
Semi-quantitative RT-PCR method
Total RNA was treated with DNAse (RQ1, RNAse free, Promega Catalys, Wallisellen, Switzerland), followed by M-MLV reverse transcriptase (RNAse H minus, point mutant, Promega, Catalys) according to the manufacturer’s recommendation and stored at −20°C. An aliquot of the cDNAs diluted at 1/10 was used in the PCR reaction. After 2 min denaturation at 95°C, 30 PCR cycles (95°C for 30 s, 54°C for 30 s and 72°C for 30 s) were run. The PCR reactions were performed in a final volume of 25 μl containing 1 μM of both forward and reverse primers, 1 U Hot Start Go Taq DNA polymerase 2x Mix (Promega, Wallisellen, Switzerland). The amplified PCR products have a size between 200 and 300 bp. As internal control, tubulin A6 (house-keeping gene) transcripts were amplified using forward 5’- ATTTGTTGACTGGTGCCCAAC and reverse 5’- TCTTCATCGTCAACTTCAGCA primers. The three other transcripts subjected to semi-quantitative RT-PCR were NtPMIa1g222869e1 (gi|48761132) corresponding to a metallophosphatase using the forward 5’-ACAATGCTAGCTTCGGTTATGG and reverse 5’-GGCTCTGCTTCTGTTTTTGTCT primers, NtPMIa1g11112e1_s matching with the protein At4g24110 using the forward 5’-TGTGGCAGCAAATATTTCAAAG and reverse 5’-TCGGATCTTGGAGTCGTTAATC primers and NtPMIa1g68186e1 corresponding to Q9ZS31 (NL27 in S. tuberosum) and the Nicotiana glutinosa virus resistance (N) gene using the forward 5’-GCAGACTGTATTCGGCATATTG and reverse 5’-TGCTATTGTGGTTTTACCCATTC primers. PCR products were separated by 1.5% agarose gel electrophoresis and visualized by staining with GelRed Nucleic Acid Stain (Chemie Brunswig, Basel, Switzerland).
Availability of supporting data
The accession number from Gene Expression Omnibus (GEO) GSE42319 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE42319) provides access to all CEL files and meta-data sheet for the experiments described in the manuscript. In addition, GEO platform GPL16290 accession number (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL16290) contains all other files (probeset annotations, Affymetrix CDF, and Affymetrix 1lq) necessary to process CEL files generated from Tobacco Exon Array experiments.
This work is dedicated to the memory of Dr. Nicolas Lugon-Moulin.
The authors would like to acknowledge Dr. Paolo Donini, Dr. Joanne Ross, and Dr. Ferruccio Gadani for supporting the Cd research at PMI and thank Dr. Ferruccio Gadani for being the principal initiator of the TGI. Finally, we would like to thank Lynda Conroy-Schneider for the extensive review of the manuscript and Dr. Sam Ansari for submission of the gene expression data to the Gene Expression Omnibus (GEO).
- Clarkson JJ, Lim KY, Kovarik A, Chase MW, Knapp S, Leitch AR: Long-term genome diploidization in allopolyploid Nicotiana section Repandae (Solanaceae). New Phytol. 2005, 168 (1): 241-252.View ArticlePubMedGoogle Scholar
- Lim KY, Matyasek R, Kovarik A, Leitch AR: Genome evolution in allotetraploid Nicotiana. Biol J Linn Soc. 2004, 82: 599-606.View ArticleGoogle Scholar
- Arumuganathan K, Earle E: Nuclear DNA content of some important plant species. Plant Mol Biol Rep. 1991, 9 (3): 208-218.View ArticleGoogle Scholar
- Kenton A, Parokonny AS, Gleba YY, Bennett MD: Characterization of the Nicotiana tabacum L. Genome by molecular cytogenetics. Mol Gen Genet. 1993, 240 (2): 159-169.View ArticlePubMedGoogle Scholar
- Zimmerman JL, Goldberg RB: DNA sequence organization in tobacco. Chromosoma. 1977, 59 (3): 227-252.View ArticleGoogle Scholar
- Bindler G, Plieske J, Bakaher N, Gunduz I, Ivanov N, Van der Hoeven R, Ganal M, Donini P: A high density genetic map of tobacco (Nicotiana tabacum L.) obtained from large scale microsatellite marker development. Theor Appl Genet. 2011, 123 (2): 219-230.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin F: An application of kernel methods to variety identification based on SSR markers genetic fingerprinting. BMC Bioinformatics. 2011, 12: 177-PubMed CentralView ArticlePubMedGoogle Scholar
- Cui H, Zhang ST, Yang HJ, Ji H, Wang XJ: Gene expression profile analysis of tobacco leaf trichomes. BMC Plant Biol. 2011, 11: 76-PubMed CentralView ArticlePubMedGoogle Scholar
- Edwards KD, Bombarely A, Story GW, Allen F, Mueller LA, Coates SA, Jones L: TobEA: an atlas of tobacco gene expression from seed to senescence. BMC Genomics. 2010, 11: 142-PubMed CentralView ArticlePubMedGoogle Scholar
- Rabinowicz PD, Citek R, Budiman MA, Nunberg A, Bedell JA, Lakey N, O'Shaughnessy AL, Nascimento LU, McCombie WR, Martienssen RA: Differential methylation of genes and repeats in land plants. Genome Res. 2005, 15 (10): 1431-1440.PubMed CentralView ArticlePubMedGoogle Scholar
- Lugon-Moulin N, Zhang M, Gadani F, Rossi L, Koller D, Krauss M, Wagner GJ: Critical review of the science and options for reducing cadmium in tobacco (Nicotiana tabacum L.) and other plants. Adv Agron. 2004, 83: 111-180.View ArticleGoogle Scholar
- Lugon-Moulin N, Martin F, Krauss MR, Ramey PB, Rossi L: Cadmium concentration in tobacco (Nicotiana tabacum L.) from different countries and its relationship with other elements. Chemosphere. 2006, 63 (7): 1074-1086.View ArticlePubMedGoogle Scholar
- Palmer LE, Rabinowicz PD, O'Shaughnessy AL, Balija VS, Nascimento LU, Dike S, de la Bastide M, Martienssen RA, McCombie WR: Maize genome sequencing by methylation filtration. Science. 2003, 302 (5653): 2115-2117.View ArticlePubMedGoogle Scholar
- Whitelaw CA, Barbazuk WB, Pertea G, Chan AP, Cheung F, Lee Y, Zheng L, van Heeringen S, Karamycheva S, Bennetzen JL, et al: Enrichment of gene-coding sequences in maize by genome filtration. Science. 2003, 302 (5653): 2118-2120.View ArticlePubMedGoogle Scholar
- Tobacco Genome Initiative. [http://www.pngg.org/tgi/],
- Spannagl M, Noubibou O, Haase D, Yang L, Gundlach H, Hindemitt T, Klee K, Haberer G, Schoof H, Mayer KF: MIPSPlantsDB--plant database resource for integrative and comparative plant genome research. Nucleic Acids Res. 2007, 35 (Database issue): D834-D840.PubMed CentralView ArticlePubMedGoogle Scholar
- Bombarely A, Menda N, Tecle IY, Buels RM, Strickler S, Fischer-York T, Pujar A, Leto J, Gosselin J, Mueller LA: The Sol genomics network (solgenomics.net): growing tomatoes using Perl. Nucleic Acids Res. 2011, 39 (Database issue): D1149-D1155.PubMed CentralView ArticlePubMedGoogle Scholar
- Smyth GK, Michaud J, Scott HS: Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics. 2005, 21 (9): 2067-2075.View ArticlePubMedGoogle Scholar
- Reiner A, Yekutieli D, Benjamini Y: Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003, 19 (3): 368-375.View ArticlePubMedGoogle Scholar
- Furukawa J, Yamaji N, Wang H, Mitani N, Murata Y, Sato K, Katsuhara M, Takeda K, Ma JF: An aluminum-activated citrate transporter in barley. Plant Cell Physiol. 2007, 48 (8): 1081-1091.View ArticlePubMedGoogle Scholar
- Liu J, Magalhaes JV, Shaff J, Kochian LV: Aluminum-activated citrate and malate transporters from the MATE and ALMT families function independently to confer Arabidopsis aluminum tolerance. Plant J. 2009, 57 (3): 389-399.View ArticlePubMedGoogle Scholar
- Ogawa I, Nakanishi H, Mori S, Nishizawa N: Time course analysis of gene regulation under cadmium stress in rice. Plant Soil. 2009, 325 (1): 97-108.View ArticleGoogle Scholar
- Li L, He Z, Pandey GK, Tsuchiya T, Luan S: Functional cloning and characterization of a plant efflux carrier for multidrug and heavy metal detoxification. J Biol Chem. 2002, 277 (7): 5360-5368.View ArticlePubMedGoogle Scholar
- Mendoza-Cozatl DG, Zhai Z, Jobe TO, Akmakjian GZ, Song WY, Limbo O, Russell MR, Kozlovskyy VI, Martinoia E, Vatamaniuk OK, et al: Tonoplast-localized Abc2 transporter mediates phytochelatin accumulation in vacuoles and confers cadmium tolerance. J Biol Chem. 2010, 285 (52): 40416-40426.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim DY, Bovet L, Kushnir S, Noh EW, Martinoia E, Lee Y: AtATM3 Is involved in heavy metal resistance in Arabidopsis. Plant Physiol. 2006, 140 (3): 922-932.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim DY, Bovet L, Maeshima M, Martinoia E, Lee Y: The ABC transporter AtPDR8 is a cadmium extrusion pump conferring heavy metal resistance. Plant J. 2007, 50 (2): 207-218.View ArticlePubMedGoogle Scholar
- Rea PA: Plant ATP-binding cassette transporters. Annual Review of Plant Biology. 2007, 58 (1): 347-375.View ArticlePubMedGoogle Scholar
- Verrier PJ, Bird D, Burla B, Dassa E, Forestier C, Geisler M, Klein M, Kolukisaoglu Ü, Lee Y, Martinoia E, et al: Plant ABC proteins – a unified nomenclature and updated inventory. Trends Plant Sci. 2008, 13 (4): 151-159.View ArticlePubMedGoogle Scholar
- Finger-Teixeira A, Lucio Ferrarese ML, Ricardo Soares A, da Silva D, Ferrarese-Filho O: Cadmium-induced lignification restricts soybean root growth. Ecotoxicol Environ Saf. 2010, 73 (8): 1959-1964.View ArticlePubMedGoogle Scholar
- Ueki S, Citovsky V: Identification of an interactor of cadmium ion-induced glycine-rich protein involved in regulation of callose levels in plant vasculature. Proc Natl Acad Sci U S A. 2005, 102 (34): 12089-12094.PubMed CentralView ArticlePubMedGoogle Scholar
- Lyubenova L, Nehnevajova E, Herzig R, Schröder P: Response of antioxidant enzymes in Nicotiana tabacum clones during phytoextraction of heavy metals. Environ Sci Pollut Res. 2009, 16 (5): 573-581.View ArticleGoogle Scholar
- Dixit P, Mukherjee PK, Ramachandran V, Eapen S: Glutathione transferase from Trichoderma virens enhances cadmium tolerance without enhancing its accumulation in transgenic Nicotiana tabacum. PLoS One. 2011, 6 (1): e16360-PubMed CentralView ArticlePubMedGoogle Scholar
- Korenkov V, Hirschi K, Crutchfield J, Wagner G: Enhancing tonoplast Cd/H antiport activity increases Cd, Zn, and Mn tolerance, and impacts root/shoot Cd partitioning in Nicotiana tabacum L. Planta. 2007, 226 (6): 1379-1387.View ArticlePubMedGoogle Scholar
- Bovet L, Eggmann T, Meylan-Bettex M, Polier J, Kammer P, Marin E, Feller U, Martinoia E: Transcript levels of AtMRPs after cadmium treatment: induction of AtMRP3. Plant Cell Environ. 2003, 26 (3): 371-381.View ArticleGoogle Scholar
- Gaillard S, Jacquet H, Vavasseur A, Leonhardt N, Forestier C: AtMRP6/AtABCC6, An ATP-binding cassette transporter gene expressed during early steps of seedling development and up-regulated by cadmium in Arabidopsis thaliana. BMC Plant Biology. 2008, 8 (1): 22-PubMed CentralView ArticlePubMedGoogle Scholar
- Wong CKE, Jarvis RS, Sherson SM, Cobbett CS: Functional analysis of the heavy metal binding domains of the Zn/Cd-transporting ATPase, HMA2, in Arabidopsis thaliana. New Phytol. 2009, 181 (1): 79-88.View ArticlePubMedGoogle Scholar
- Ueno D, Yamaji N, Kono I, Huang CF, Ando T, Yano M, Ma JF: Gene limiting cadmium accumulation in rice. Proc Natl Acad Sci. 2010, 107 (38): 16500-16505.PubMed CentralView ArticlePubMedGoogle Scholar
- Mustroph A, Lee SC, Oosumi T, Zanetti ME, Yang H, Ma K, Yaghoubi-Masihi A, Fukao T, Bailey-Serres J: Cross-kingdom comparison of transcriptomic adjustments to low-oxygen stress highlights conserved and plant-specific responses. Plant Physiol. 2010, 152 (3): 1484-1500.PubMed CentralView ArticlePubMedGoogle Scholar
- van Dongen JT, Fröhlich A, Ramírez-Aguilar SJ, Schauer N, Fernie AR, Erban A, Kopka J, Clark J, Langer A, Geigenberger P: Transcript and metabolite profiling of the adaptive response to mild decreases in oxygen concentration in the roots of Arabidopsis plants. Ann Bot. 2009, 103 (2): 269-280.PubMed CentralView ArticlePubMedGoogle Scholar
- Moorhead GB, De Wever V, Templeton G, Kerk D: Evolution of protein phosphatases in plants and animals. Biochem J. 2009, 417 (2): 401-409.View ArticlePubMedGoogle Scholar
- Love CG, O'Lochlainn S, Graham NS, Bowen HC, May ST, White PJ, Broadley MR, Hammond JP, King GJ: A brassica exon array for whole-transcript gene expression profiling. PLoS One. 2010, 5 (9): e12812-[http://dx.doi.org/10.1371/journal.pone.0012812]PubMed CentralView ArticlePubMedGoogle Scholar
- Bemmo A, Benovoy D, Kwan T, Gaffney D, Jensen R, Majewski J: Gene expression and isoform variation analysis using Affymetrix exon arrays. BMC Genomics. 2008, 9 (1): 529-PubMed CentralView ArticlePubMedGoogle Scholar
- Whistler T, Chiang C-F, Lonergan W, Hollier M, Unger E: Implementation of exon arrays: alternative splicing during T-cell proliferation as determined by whole genome analysis. BMC Genomics. 2010, 11 (1): 496-PubMed CentralPubMedGoogle Scholar
- Ha K, Coulombe-Huntington J, Majewski J: Comparison of Affymetrix gene array with the exon array shows potential application for detection of transcript isoform variation. BMC Genomics. 2009, 10 (1): 519-PubMed CentralView ArticlePubMedGoogle Scholar
- Liu X, Peng K, Wang A, Lian C, Shen Z: Cadmium accumulation and distribution in populations of Phytolacca americana L. and the role of transpiration. Chemosphere. 2010, 78 (9): 1136-1141.View ArticlePubMedGoogle Scholar
- Uraguchi S, Mori S, Kuramata M, Kawasaki A, Arao T, Ishikawa S: Root-to-shoot Cd translocation via the xylem is the major process determining shoot and grain cadmium accumulation in rice. J Exp Bot. 2009, 60 (9): 2677-2688.PubMed CentralView ArticlePubMedGoogle Scholar
- Lux A, Martinka M, Vaculík M, White PJ: Root responses to cadmium in the rhizosphere: a review. J Exp Bot. 2011, 62 (1): 21-37.View ArticlePubMedGoogle Scholar
- Hanikenne M, Talke IN, Haydon MJ, Lanz C, Nolte A, Motte P, Kroymann J, Weigel D, Kramer U: Evolution of metal hyperaccumulation required cis-regulatory changes and triplication of HMA4. Nature. 2008, 453 (7193): 391-395.View ArticlePubMedGoogle Scholar
- Marmiroli M, Visioli G, Maestri E, Marmiroli N: Correlating SNP genotype with the phenotypic response to exposure to cadmium in Populus spp. Environ Sci Technol. 2011, 45 (10): 4497-4505.View ArticlePubMedGoogle Scholar
- Studart-Guimarães C, Fait A, Nunes-Nesi A, Carrari F, Usadel B, Fernie AR: Reduced expression of succinyl-coenzyme A ligase can be compensated for by up-regulation of the γ-aminobutyrate shunt in illuminated tomato leaves. Plant Physiol. 2007, 145 (3): 626-639.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsolais F, Pajak A, Yin F, Taylor M, Gabriel M, Merino DM, Ma V, Kameka A, Vijayan P, Pham H, et al: Proteomic analysis of common bean seed with storage protein deficiency reveals up-regulation of sulfur-rich proteins and starch and raffinose metabolic enzymes, and down-regulation of the secretory pathway. Journal of Proteomics. 2010, 73 (8): 1587-1600.View ArticlePubMedGoogle Scholar
- Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27 (2): 573-580.PubMed CentralView ArticlePubMedGoogle Scholar
- Smit A, Hubley R, Green P: RepeatMasker open-3.0 (1996-2010). [http://www.repeatmasker.org]
- Stanke M, Diekhans M, Baertsch R, Haussler D: Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008, 24 (5): 637-644.View ArticlePubMedGoogle Scholar
- Salamov AA, Solovyev VV: Ab initio gene finding in drosophila genomic DNA. Genome Res. 2000, 10 (4): 516-522.PubMed CentralView ArticlePubMedGoogle Scholar
- Price AL, Jones NC, Pevzner PA: De novo identification of repeat families in large genomes. Bioinformatics. 2005, 21 (Suppl 1): i351-358.View ArticlePubMedGoogle Scholar
- Blanchette M, Kent WJ, Riemer C, Elnitski L, Smit AF, Roskin KM, Baertsch R, Rosenbloom K, Clawson H, Green ED, et al: Aligning multiple genomic sequences with the threaded blockset aligner. Genome Res. 2004, 14 (4): 708-715.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwartz S, Kent WJ, Smit A, Zhang Z, Baertsch R, Hardison RC, Haussler D, Miller W: Human-mouse alignments with BLASTZ. Genome Res. 2003, 13 (1): 103-107.PubMed CentralView ArticlePubMedGoogle Scholar
- Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, et al: Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005, 15 (8): 1034-1050.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan W, Kang MS: Biplot analysis of multi-environment trial data. GGE biplot analysis: A Graphical Tool for Breeders, Geneticists, and Agronomists. Edited by: Yan W, Kang MS. 2002, Boca Raton: CRC Press, 63-99.View ArticleGoogle Scholar
- Wood SN: Stable and efficient multiple smoothing parameter estimation for generalized additive models. J Am Stat Assoc. 2004, 99 (467): 673-686.View ArticleGoogle Scholar
- Heber S, Sick B: Quality assessment of Affymetrix GeneChip data. OMICS. 2006, 10 (3): 358-368.View ArticlePubMedGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264.View ArticlePubMedGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (1): Article 3 (Online) http://dx.doi.org/10.2202/1544-6115.1027Google Scholar