Skip to main content

Large-scale analysis of antisense transcription in wheat using the Affymetrix GeneChip Wheat Genome Array

Abstract

Background

Natural antisense transcripts (NATs) are transcripts of the opposite DNA strand to the sense-strand either at the same locus (cis-encoded) or a different locus (trans-encoded). They can affect gene expression at multiple stages including transcription, RNA processing and transport, and translation. NATs give rise to sense-antisense transcript pairs and the number of these identified has escalated greatly with the availability of DNA sequencing resources and public databases. Traditionally, NATs were identified by the alignment of full-length cDNAs or expressed sequence tags to genome sequences, but an alternative method for large-scale detection of sense-antisense transcript pairs involves the use of microarrays. In this study we developed a novel protocol to assay sense- and antisense-strand transcription on the 55 K Affymetrix GeneChip Wheat Genome Array, which is a 3' in vitro transcription (3'IVT) expression array. We selected five different tissue types for assay to enable maximum discovery, and used the 'Chinese Spring' wheat genotype because most of the wheat GeneChip probe sequences were based on its genomic sequence. This study is the first report of using a 3'IVT expression array to discover the expression of natural sense-antisense transcript pairs, and may be considered as proof-of-concept.

Results

By using alternative target preparation schemes, both the sense- and antisense-strand derived transcripts were labeled and hybridized to the Wheat GeneChip. Quality assurance verified that successful hybridization did occur in the antisense-strand assay. A stringent threshold for positive hybridization was applied, which resulted in the identification of 110 sense-antisense transcript pairs, as well as 80 potentially antisense-specific transcripts. Strand-specific RT-PCR validated the microarray observations, and showed that antisense transcription is likely to be tissue specific. For the annotated sense-antisense transcript pairs, analysis of the gene ontology terms showed a significant over-representation of transcripts involved in energy production. These included several representations of ATP synthase, photosystem proteins and RUBISCO, which indicated that photosynthesis is likely to be regulated by antisense transcripts.

Conclusion

This study demonstrated the novel use of an adapted labeling protocol and a 3'IVT GeneChip array for large-scale identification of antisense transcription in wheat. The results show that antisense transcription is relatively abundant in wheat, and may affect the expression of valuable agronomic phenotypes. Future work should select potentially interesting transcript pairs for further functional characterization to determine biological activity.

Background

Natural antisense transcripts (NATs) are defined as transcripts of the opposite DNA strand to the sense-strand either at the same locus (cis-encoded) or a different locus (trans-encoded). The first NATs were detected in viruses, followed by prokaryotes and then eukaryotes. For an excellent review of current NAT knowledge, please refer to Lapidot and Pilpel [1]. NATs usually possess a negative regulatory effect and can affect gene expression at multiple stages including transcription, RNA processing and transport, and translation [2, 3]. Thus, NATs may be involved in the regulation of varying biological functions such as the adaptation to stresses and development. NATs are involved in RNA interference [4, 5], methylation [6] and genomic imprinting [7]. NATs give rise to sense-antisense transcript pairs that were once considered as rare, but the number identified has escalated greatly with the availability of DNA sequencing resources and public databases. For example, 22% of annotated genes in the fruit fly genome are reported to overlap as transcript pairs [8], and more than 20% of human transcripts may form sense-antisense transcript pairs [9]. In plants, few sense-antisense transcript pairs had been reported until recent large-scale studies in rice [10, 11] and A. thaliana [12, 13]. In the rice study, full-length cDNA data revealed that approximately 7% of transcripts formed sense-antisense transcript pairs [10]. In these plant studies, the alignment of full-length cDNAs and expressed sequence tags (ESTs) to the genome sequence was used to identify the sense-antisense transcript pairs, which is limited to the detection of cis-encoded pairs. In wheat, antisense transcripts have been discovered from serial analysis of gene expression (SAGE) tags of developing grain [14], where it was reported that 25.7% of forward (sense) tags had a matching reverse (antisense) tag, which indicated widespread antisense transcription in wheat.

An alternative method for large-scale discovery of sense-antisense transcript pairs involves the use of microarrays. In the first study of this type, Yelin et al. [15] used a strand-specific oligonucleotide probe array to detect antisense transcription in human cell lines. A study in mouse using a custom oligonucleotide array to assay the expression of 1,947 known sense-antisense transcript pairs has also been reported [16]. However, these studies required prior knowledge of the sense-antisense transcript pairs to enable the design of strand specific probes. To overcome this, Werner et al. [17] took advantage of the approximately 25% of incorrectly orientated probes on the Affymetrix GeneChip U74A and U74B 3'in vitro transcription (3'IVT) mouse arrays to detect novel antisense transcription in mouse brain and kidney tissues. The results showed that the commercial expression arrays were sensitive enough to detect antisense transcription, but because it cannot be assumed that current commercial arrays contain incorrectly orientated probes, this type of study could not be repeated. Subsequently, Ge et al. [18] developed a method called 'Antisense Transcriptome analysis using Exon array (ATE)' that used an altered target synthesis and labeling method that allowed both sense- and antisense-strand transcription to be assayed on Affymetrix Whole-Transcript Expression arrays (ie. Exon and Gene arrays). This protocol was successful but cannot be applied to the numerous Affymetrix (3'IVT) expression arrays, because these arrays are constructed with probes of the opposite strand to the Whole-Transcript Expression arrays, thus they use a different target labeling procedure altogether.

In the current study, we sought to develop a protocol that could be used to assay sense- and antisense-strand transcription on the Affymetrix GeneChip Wheat Genome array, which is a 3'IVT expression array. The 3'IVT expression arrays rely on in vitro transcription of double-stranded cDNA to both amplify and label the target cRNA before hybridization. The wheat array currently provides the most comprehensive coverage of the wheat genome for a microarray and is a commonly used resource for transcript expression studies [19, 20] and hybridization-based DNA marker discovery [21]. This study is the first report of using a 3'IVT expression array to discover the expression of natural sense-antisense transcript pairs without relying on the presence of incorrectly oriented probes, and may be considered as proof-of-concept. By using alternative target preparation schemes, both the sense- and antisense-strand derived transcripts were labeled and hybridized to the Wheat Genome Array. To enable maximum discovery we selected five different tissue types for assay and used 'Chinese Spring' wheat genotype, since most of the GeneChip probe sequences were based on its genomic sequence. The functional annotation of detected wheat sense-antisense transcript pairs is discussed, as well as the performance and validation of the technique.

Results

Target preparation

Total RNA was extracted from five 'Chinese Spring' tissue types (germinated seed, shoot, flag leaf, spike pre-anthesis and spike post-anthesis; see materials and methods). In addition to maximizing discovery, these tissue types were also selected to align with the predominant tissues used for wheat EST sequencing efforts, including the International Wheat Genome Sequencing Consortium (IWGSC). All samples were of excellent quality as assessed by gel electrophoresis and spectrophotometry. The total RNA samples were mixed at equal concentrations before target preparation. The assay of sense-strand transcription followed the regular scheme as for all Affymetrix 3'IVT GeneChips (see materials and methods). However, to assay antisense-strand transcription, the Affymetrix Whole Transcript (WT) Sense Target Labeling Assay was used, which was designed specifically for use on Whole-Transcript Expression arrays. The WT target preparation method resulted in labeling the opposite strand to the 3'IVT assay and was therefore used in this study to assess antisense-strand transcription (Figure 1). The mixed total RNA sample was used as starting material for both the 3'IVT and WT target preparation, and each hybridization was carried out once according to standard Affymetrix protocol for the Wheat Genome Array.

Figure 1
figure 1

Principles of the two target preparation methods used to assay both sense- and antisense-strand transcription. A 5' (head-to-head) overlapping sense-antisense transcript pair is used as an example. The standard Affymetrix 3' in vitro transcription (3'IVT) assay was used to detect sense-strand transcription, while a modified Affymetrix Whole Transcript (WT) assay was used to detect antisense-strand transcription.

Data analysis

Following hybridization and scanning, CEL files were analyzed to identify probe sets that showed successful hybridization in each of the 3'IVT and WT assay. The quality control metrics from the affyQCreport package [22] of Bioconductor [23] showed that data was of high-quality for both assays, but as expected the WT assay resulted in a lower percentage of detected transcripts (16.27%) than the 3'IVT assay (47.95%) using Affymetrix PMA (Present/Marginal/Absent) calls. Subsequently the relationship of array distributions (Figure 2) showed a skewing towards the 3'IVT assay, but it is clear that successful hybridization did occur in the WT assay. Figure 2 also showed that the spiked-in hybridization controls from Affymetrix (bioB, bioC, bioD and creX) produced similar signals from both assays although signals were slightly higher in the WT assay. The MAS 5.0 PMA calls and Robust Multi-array Average (RMA) summarized expression values were used to determine successful hybridization for each probe set. Because of differences in the two labeling methods, including starting amount of RNA and RNA amplification, the expression values of each array could not be validly compared. However, the PMA calls in combination with the expression values were used to determine positive hybridization to a particular probe set in each assay. This provided a qualitative measure of expression rather than quantitative, but for the purposes of this study which was to detect natural sense-antisense transcript pairs, this measure was satisfactory.

Figure 2
figure 2

Comparison of the 3'IVT (sense signal) and WT (antisense signal) array distributions. The summarized log2 expression value for all probe sets (n = 61,127). Red spots indicate probe sets with absent MAS 5.0 PMA calls (P < 0.01) in both assays, black spots indicate probe sets called as present in one assay only, blue spots indicate probe sets with present calls in both assays, and green spots indicate the replicated spiked-in hybridization controls (bioB, bioC, bioD and creX). The black box shows the probe sets with expression values greater than the spiked-in hybridization controls in both assays, which were considered as positively hybridizing for this study.

Identifying sense-antisense transcript pairs

To determine a confident positive threshold for expression value in both assays, the expression values of spiked-in hybridization controls (bioB, bioC, bioD and creX) were used. Because these controls are spiked-in immediately before hybridization, they were expected to behave in the same way in both assays. The bioB control is spiked-in at the detection limit, while the others are spiked-in at staggered concentrations after bioB. Thus, the log2 expression value of bioB was considered the threshold for positive hybridization in each assay (Figure 2). Because the bioB probe set is replicated three times on the wheat GeneChip, the log2 expression value of the lowest individual probe set was used. For the 3'IVT assay this value was 8.76, and for the WT assay it was 9.86. These cut-off values were used in combination with the MAS 5.0 PMA calls and corresponding probability (p) values to detect successful hybridization in each assay. In both assays, a probe set must have firstly been called 'Present' with Wilcoxon rank sum test p- value < 0.01, and the RMA summarized log2 expression value must have been greater than 8.76 in the 3'IVT assay, and greater than 9.86 in the WT assay. This threshold cutoff identified 110 probe sets as positively hybridizing in both the 3'IVT and WT assays. In addition to the 110 probe sets, 8940 probe sets uniquely hybridized in the 3'IVT assay and 80 uniquely hybridized to the WT assay (potentially antisense-specific transcripts). Because the aim of this study was detect transcript pairs transcribed from both strands, we mainly focused on probe sets detected in both assays. These stringent detection criteria ensured that the probe sets left were highly expressed in both assays and could more reliably be considered as sense-antisense transcript pairs. In fact, the 80 antisense-specific probe sets could not necessarily be classified as antisense transcripts, because these may represent incorrectly orientated probes. Also, because the probes for a given transcript do not cover the entire sequence, there is a possibility for bias during hybridization. However, to form the basis of future studies these 80 probe sets were also given some attention.

Annotation of probe sets

Each of the 110 candidate sense-antisense transcript pair probe sets were functionally annotated using HarvEST (Affymetrix Wheat1 Chip version 1.52). Gene Ontology (GO) was based on the TIGR rice genome annotation such that if a unigene possessed a significant (<1e-10) BLASTx match to rice, as identified in HarvEST, the corresponding GO terms for the rice protein were used, if available. Of the 110 probe sets 76 could be annotated (see Additional file 1), of which 46 (59%) were classified as involved in energy production ('Energy'), including several representations of ATP synthase, photosystem proteins and RUBISCO. To determine the significance of overrepresentation of the number of energy-related transcripts identified, a hypergeometric test of selected energy-related terms in the HarvEST annotated transcript description were used (see 'Methods'). For the transcripts identified as present in one or both hybridizations, energy-related terms were identified in 1831 of the 24578 transcripts (7.4%) that possessed a transcript description. Using the same search terms, 46 of the 76 annotated probe sets identified in this study possessed energy-related terms. Subsequently, energy-related transcripts were found to be significantly over-represented in this study with a p-value of 4.88 × 10-37. The diversity of the annotated probe sets is summarized in Table 1.

Table 1 Summary of the 76 successfully detected and annotated sense-antisense transcript pairs in both the 3'IVT and WT assay.

The 80 potential antisense-specific probe sets were also annotated as described for the transcript pairs. Of the 80 probe sets only 31 could be annotated (see Additional file 2), of which 10 (32%) were classified as involved in energy production ('Energy') including several representations of RUBISCO. Nine (28%) were involved in transcription ('Transcription') and included several DNA-directed RNA polymerase transcripts. However, the majority of antisense-specific transcripts were of unknown function.

Strand-specific transcription validation

Ten probe sets selected to represent a range of functional categories were validated for sense- and antisense-strand transcription using strand-specific reverse transcription-PCR (RT-PCR). An example of the electrophoresis results is shown in Figure 3. Sense-strand transcription was detected for all 10 targets sets in each tissue except for the target RUBISCO activase in the 'Germinated seed' tissue (Table 2). In fact, the 'Germinated seed' tissue was most different to the other tissues and showed the least amount of antisense-strand transcription for the 10 targets. The 'Shoot', 'Flag leaf', 'Spike pre-anthesis' and 'Spike post-anthesis' tissues all showed the same pattern of sense- and antisense-strand transcription. These results indicate that antisense-strand transcription is likely specific to certain tissues and/or developmental stages, although not to a great extent in the 10 target transcripts analyzed in this study. Only one of the 10 targets (10%) did not show any antisense-strand transcription in any tissue, thus was not in agreement with the microarray results. However, this could be due to the position of the RT-PCR primer for amplifying the antisense-strand transcript. Because antisense-strand transcripts may not necessarily span the full-length of their complementary sense-strand transcript, the RT-PCR primer may have been targeted to a missing region in the antisense-strand transcript, thus the RT-PCR would fail.

Figure 3
figure 3

Strand-specific RT-PCR for target Carbonic anhydrase (Ta.5227.1). Sense-strand transcription was detected only for germinated seed tissue, and both sense- and antisense-strand transcription for all other tissues. M: DNA marker, 1: Germinated seed antisense, 2: Germinated seed sense, 3: Shoot antisense, 4: Shoot sense, 5: Flag leaf antisense, 6: Flag leaf sense, 7: Spike pre-anthesis antisense, 8: Spike pre-anthesis sense, 9: Spike post-anthesis antisense, 10: Spike post-anthesis sense, N: Negative controls.

Table 2 Validation of selected sense-antisense transcript pairs with strand-specific reverse transcription PCR (RT-PCR).

Discussion

This study reports on the first use of an Affymetrix GeneChip 3'IVT expression array for discovering both sense- and antisense-strand transcription. Through the adaptation of the Affymetrix WT assay, the antisense transcribed strand was successfully labeled and hybridized to the Wheat Genome Array, which allowed for the detection of natural sense-antisense transcript pairs. To our knowledge, the Wheat Genome Array does not contain any probes for known sense-antisense transcript pairs, thus the data from the hybridizations could not be standardized and/or normalized to a known sense-antisense transcript pair. Subsequently, a highly stringent data acceptance threshold was applied, based on PMA call and expression value cutoffs. This increased the confidence in detecting true antisense transcription. It is important to recognize the limitations of this study, which stem from the 'closed' nature of microarray systems. Because the Wheat Genome Array contains only known transcript sequences, the study is clearly limited to detection of transcript pairs that are present on the array. Further, the probes for each transcript are biased to the 3' end of transcripts and do not span the entire gene. Thus, because antisense-strand transcripts commonly have a different splice structure they may not be detected. Subsequently the 110 candidate sense-antisense transcript pairs and the 80 potentially antisense-specific transcripts that were identified are likely to under-represent the number of true transcript pairs. In future studies, custom microarrays containing probes for sense and antisense transcripts would be useful as different target preparation assays would not be required, but because we aimed to obtain a broad representation of the extent of antisense transcription we chose to use the most comprehensive Wheat Genome Array.

The function of antisense-strand transcription is widely believed to regulate the expression of sense-strand transcripts at either transcription, mRNA maturation or translation [2]. In fact, Lapidot and Pilpel [1] reviewed the literature and postulated four mechanisms of action; transcriptional interference, RNA masking, double-stranded RNA (ds-RNA)-dependent mechanisms, and chromatin remodeling. The ds-RNA mechanisms would likely be the result of RNA-dependent RNA polymerases, which generate ds-RNA that are the precursors of short interfering RNA (siRNA). The timing of sense- and antisense-strand transcription is also important; for example, if the sense-strand is transcribed first up to a certain level followed by transcription of the antisense-strand, the biological result would be delayed inhibition of the sense-strand gene expression. Conversely, if the antisense-strand was transcribed first, this would result in pre-inhibition of sense-strand gene expression up to a threshold. Differences in the half-life of the sense-and antisense-strand transcripts, as well as tissue-specificity and potential light and/or diurnal transcript regulation [24] would also affect these scenarios. In the present study the timing of transcription and relative level of sense- and antisense-transcripts could not be determined because a single time-point was used for RNA extraction in each tissue, and the design of the assay did not allow valid comparisons between the 3'IVT and WT results to estimate transcript levels. Thus the mode of action of the detected sense-antisense transcript pairs would require further study.

An important observation in this study was the functional annotation of the sense-antisense transcript pairs, which indicated a significant over-representation of those involved in energy production, particularly photosynthesis. Additionally, many transcripts for ribosomal proteins involved in protein synthesis were identified. The abundance of antisense transcripts for these common plant processes may indicate that they are negatively regulated by antisense transcripts. Alternatively, the antisense transcripts could possibly be the result of ectopic expression. There is little data on large-scale antisense transcription profiling in plants to compare these results with, but a study in rice of leaf and seed tissue using Serial Analysis of Gene Expression (SAGE) identified sense-antisense transcript pairs and also found that the most abundant pairs were annotated as involved in energy production, including RUBISCO and a Photosystem I protein [11]. The similarity between studies shows that transcripts involved in photosynthesis are likely to be controlled by antisense transcripts in plants. An appealing explanation is the possibility for diurnal regulation of photosynthesis through antisense regulation. Although this study did not span a time-course required to demonstrate diurnal regulation, the results warrant further exploration of this hypothesis.

The results of the strand-specific RT-PCR also showed that antisense transcription is likely to be tissue-specific. Only one of the RT-PCR results was not in complete agreement with the microarray result, which could be due to truncated antisense transcripts where the priming sites were absent. In their microarray study of human cell antisense transcription, Ge et al. [18] found that 26% of the RT-PCR results were not consistent with microarray observations. In this study we also identified 80 transcripts as potentially antisense-specific, although further studies would be needed to confirm this because of the possibility for incorrectly oriented probes or strand bias during hybridization. The majority of these transcripts were annotated as unknown, but of those that were there was again a trend towards function in photosynthesis. A high percentage were also functionally involved in controlling transcription, including transcripts with homology to DNA-directed RNA polymerase, which indicates that gene expression in wheat may be regulated by antisense transcripts at the transcriptional level.

A recent study in wheat involving SAGE of developing grain also identified antisense transcripts [14], where the most abundant functional categories aside from unknown tags were associated with storage and reproduction. The abundance of these functional categories was due to the sampling of developing grain tissue, while the abundance of energy-related transcripts in our study is most likely due to the selection of photosynthetic tissues. For this reason, these two studies complement each other well. As in our study, Poole et al. [14] found that most antisense tags were of unknown function and that many transcripts were highly expressed in both sense and antisense, which may suggest a function of the antisense transcript for mediating alternative polyadenylation rather than down-regulation of the sense transcript, although there is no evidence for this at this stage. One other similarity to our study was the identification by Poole et al. [14] of antisense transcripts related to transcription, such as nucleotide binding proteins, which the authors suggest may enable the control of multiple pathways that require large scale changes during development. Other than these similarities, the results of our study differ from Poole et al. [14], which again is likely due to the complementary tissues analyzed.

This study was exploratory and revealed that the method was successful in identifying sense-antisense transcript pairs using the commercial Wheat Genome Array. The next step from this study is to select potentially interesting antisense transcripts for further study. There were several transcript pairs belonging to functional categories including 'Cell death' and 'Transcription' that may be involved in the regulation of important biological processes, and the antisense-specific transcripts related to transcription are also of interest. An understanding of the role of antisense transcription as it relates to gene expression may be important for the expression of certain phenotypes of interest. Additionally, knowledge of natural antisense transcripts may also be important for altering gene expression through transgenic studies in plants. The abundance of antisense-strand transcripts in plants is supported by recent studies using 'open' transcriptomics systems including SAGE [11, 14] and Massively Parallel Signature Sequencing (MPSS) [12]. With the advent of RNA-Seq (RNA sequencing), which is high-throughput transcriptome sequencing method [25] that incorporates the use of next-generation sequence-by-synthesis technologies, the future will see a greatly enhanced discovery and understanding of antisense-strand transcription in plants.

Conclusion

This study demonstrated the novel use of an adapted labeling protocol and a 3'IVT Affymetrix GeneChip microarray for large-scale identification of antisense transcription in wheat, a crop of great economic importance. The results show that antisense transcription is relatively abundant in wheat, and may affect the expression of valuable agronomic phenotypes. Strand-specific RT-PCR validated the microarray observations, and showed that antisense transcription is likely to be tissue specific. Most of the identified sense-antisense transcript pairs were annotated as genes involved in energy production, indicating that photosynthesis is likely to be under regulation by antisense transcripts.

Methods

Plant material and RNA extraction

The spring wheat genotype 'Chinese Spring' was selected for this study because the majority of GeneChip Wheat Genome Array (Affymetrix, Santa Clara, California, USA) probe sequences were based on its DNA sequence. Five tissue types were selected for this study; i.Germinated seed (germinated on wetted filter paper in a petri-dish in the dark for two days, radicle and plumule emerged), ii.Shoot (Feekes 1.2, emergence with second leaf unfolded), iii.Flag leaf (Feekes 9.0, flag leaf stage with ligule emerged), iv.Spike pre-anthesis (Feekes 10.5, heading complete), and v.Spike post-anthesis (Feekes 10.5.3, flowering complete to base of spike). For each tissue type except 'Germinated seed', six seeds were planted in 6-inch round pots using a potting mix (6 peat moss: 4 vermiculite with lime: 3 sand: 3 commercial potting mix: 2 perlite: 1.7 g/L lime: 3.3 g/L Osmocote: 2.2 g/L ammonium nitrate). Seedlings were grown to the appropriate growth stages for tissue collection in a greenhouse with a stable temperature of 25 ± 2°C and a 16 h light/8 h dark cycle. All tissue samples were immediately frozen in liquid nitrogen, and total RNA was extracted from 1.0 g of each pooled tissue type using the Trizol® Plus RNA Purification Kit (Invitrogen, Carlsbad, CA) with an on-column DNase treatment. Purified total RNA samples were quantified with a NanoDrop® ND-1000 (NanoDrop, Wilmington, DE) spectrophotometer, and satisfactory purity was indicated by A260:280 ratios of 1.9–2.1 in 10 mM Tris-Hcl (pH 7.5). Integrity of total RNA samples was assessed by denaturing formaldehyde gel electrophoresis, where the presence of sharp 28S and 18S ribosomal RNA bands at an intensity ratio of ~2:1 (28S:18S) indicated good integrity. The five total RNA samples were then mixed at equal concentrations before analysis.

Wheat Genome Array

The GeneChip Wheat Genome Array (Affymetrix, Santa Clara, California, USA) is a 3'IVT array that includes 61,127 probe sets representing 55,052 transcripts for all 21 wheat chromosomes in the genome. 59,356 probes sets represent modern hexaploid (A, B and D genomes) bread wheat (T. aestivum) and are derived from the public content of the T. aestivum UniGene Build #38 (April 24, 2004). 1,215 probe sets are derived from ESTs of a diploid near relative of the A genome (T. monococcum), a further 539 represent ESTs of the tetraploid (A and B genomes) macaroni wheat species T. turgidum, and five are from ESTs of a diploid near relative of the D genome known as Aegilops tauschii. Probe sets consisted of pairs of 11 perfect match (PM) and mismatch (MM) 25-mer oligonucleotides designed from the 3' end of exemplar sequences, with nucleotide 13 as the MM. Array annotation information is available on the NetAffx data analysis center http://www.affymetrix.com.

Sense-strand transcription analysis

The 3'IVT Wheat Genome Array detects sense strand transcription by generating antisense-orientated labeled complementary RNA (cRNA) from the original RNA sample that is then hybridized to the probes that are designed to hybridize to the antisense-orientated labeled cRNA. Although the system generates antisense-orientated labeled cRNA, the assayed strand for transcription is the sense strand. Thus, to assay sense strand transcription in the mixed Chinese Spring total RNA sample, the regular 3'IVT Affymetrix protocol was carried out http://www.affymetrix.com. Briefly, double-stranded cDNA was generated from mRNA using a T7-oligo(dT) primer. The double-stranded cDNA was cleaned up and used as template for in vitro transcription (IVT) in the presence of T7 RNA Polymerase and a biotinylated nucleotide analog/ribonucleotide mix for cRNA amplification and biotin labeling. The biotinylated cRNA target was then cleaned up, fragmented, and hybridized to the Wheat Genome Array. All hybridizations and data acquisition was performed at the Genomics Core Facility at Washington State University (Pullman, Washington, USA) according to standard Affymetrix protocols http://www.biotechnology.wsu.edu/Core_Laboratories.aspx#.

Antisense-strand transcription analysis

To assay antisense strand transcription, the Affymetrix Whole Transcript (WT) Sense Target Labeling Assay http://www.affymetrix.com was used. This WT assay is intended for use on Affymetrix Whole Transcript expression arrays, which contain probes designed to hybridize with sense-orientated labeled cDNA. Because the probes of the 3'IVT Wheat Genome Array are designed to hybridize to antisense-orientated labeled cRNA derived from sense strand transcription, the sense-orientated labeled cDNA generated by the WT assay will not hybridize to the array unless it was derived from antisense transcription. Thus to discover antisense transcription in the mixed Chinese Spring total RNA sample, the target was prepared using the WT assay but hybridized to the 3'IVT Wheat Genome Array. Briefly, double-stranded cDNA was synthesized with random hexamers coupled to the T7 promoter, followed by IVT amplification with T7 RNA polymerase to produce cRNA. A second cycle cDNA synthesis was then performed using random primers for reverse transcription, which converted the cRNA into single-stranded cDNA in the same orientation as the original mRNA (sense-orientation). The single-stranded cDNA was then cleaned up, fragmented, and hybridized to the Wheat Genome Array according to standard Affymetrix protocol.

Data analysis

Using GeneChip Operating Software (GCOS) v.1.4 (Affymetrix, Santa Clara, California, USA), image quality control was performed by inspecting raw intensity (DAT) files for scratches/smears and uniform performance of the B2 oligo around the border of each image. Data quality control was from raw data in CEL files using the affyQCreport package [22] of Bioconductor [23], which provided Affymetrix recommended quality metrics, per array intensity distributions, between array comparisons, and other diagnostic plots for each hybridization. The Bioconductor [23] package affy [26] was used to read in the raw Affymetrix 'CEL' files, which were pre-processed using Robust Multi-array Average (RMA) [27, 28]. Pre-processing was modified so that only expression value summarization was applied. Background correction and normalization were omitted because the arrays were hybridized using different labeling assays. PMA (present, marginal and absent) calls were calculated for each probe set using a Wilcoxon rank sum test from MAS 5.0 [29]. Only those probe sets with a Wilcoxon p-value < 0.01 were considered 'present'.

Probe sets called as present were also required to possess a summarized log2 expression value greater than the bioB spiked-in hybridization control (>8.76 in the 3'IVT assay and >9.86 in the WT assay). Probe sets meeting these criteria were annotated using HarvEST (Affymetrix Wheat1 Chip version 1.52), which identified the corresponding unigene for each probe set and provided the current best BLASTX hit from the non-redundant (nr) database of NCBI, as well as the best BLASTX hits from rice and Arabidopsis thaliana TIGR databases http://www.tigr.org/plantProjects.shtml. A database hit <1e-10 was considered as significant, otherwise the unigene was annotated as 'no homology'. Unigenes were assigned to functional categories based on Munich Information Center for Protein Sequences (MIPS; http://mips.gsf.de/) classifications. All minimum information about microarray experiments (MIAME) guidelines were observed and GeneChip data was deposited into WheatPLEX [30] accession TA21, as well as NCBI's Gene Expression Omnibus [31] accession number GSE12528. For gene ontology (GO), the rice locus matching each probe set in the HarvEST output provided the most comprehensive annotation set. To assess the significance of energy-related transcripts, common terms in the rice transcript description for energy-related transcripts were selected and used as search terms across the annotation of transcripts found to be present in one or both assays. The search terms used were: photosystem, ribosomal protein, chloroplast, chlorophyll, cp12, oxygen-evolving, carbonic anhydrase, ATP synthase, ribulose, cytochrome, NADH. Each probe set was inspected as to whether or not it contained one or more the 'energy' search terms. To assess whether 'energy' related transcripts were significantly overrepresented in the identified sense-antisense transcript pairs than was expected by random chance, we performed a hypergeometric test.

Strand-specific transcription validation

Ten probe sets that were found to be transcribed on both strands were selected for validation using strand-specific reverse transcription-PCR (RT-PCR) (Table 2). The strand-specificity of the Qiagen One Step RT-PCR kit (Qiagen, Valencia, California, USA) has been confirmed in previous studies [32], thus was selected for use in this study. DNase-treated total RNA from each tissue type was used to determine potential tissue specificity of transcription. Unigene sequences for each of the 10 probe sets were identified using HarvEST (Affymetrix Wheat1 Chip version 1.52), and primer pairs were designed using Vector NTI (v. 10.3.0, Invitrogen Corporation). Strand-specificity was achieved by selective use of primers in the reverse transcription step, where the reverse primer was used to detect sense transcripts, and the forward primer to detect antisense transcripts. PCR reactions following reverse transcription were carried out in the presence of both forward and reverse primers, with the following cycling parameters; i. 50°C for 60 min (reverse transcription), ii. 95°C for 15 min (activate polymerase and deactivate RT enzymes), iii. 4°C for 5 min (added missing primer/s for PCR at this point), iv. 94°C for 30 s, 60°C for 30 s, 72°C for 45 s (PCR cycling repeated 35 times), and v. 72°C for 10 min (final extension). All PCR products were visualized on 1.5% agarose gels.

References

  1. Lapidot M, Pilpel Y: Genome-wide natural antisense transcription: coupling its regulation to its different regulatory mechanisms. EMBO Rep. 2006, 7 (12): 1216-1222. 10.1038/sj.embor.7400857.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Terryn N, Rouze P: The sense of naturally transcribed antisense RNAs in plants. Trends in Plant Science. 2000, 5 (9): 394-396. 10.1016/S1360-1385(00)01696-4.

    Article  CAS  PubMed  Google Scholar 

  3. Vanhee-Brossollet C, Vaquero C: Do natural antisense transcripts make sense in eukaryotes?. Gene. 1998, 211 (1): 1-9. 10.1016/S0378-1119(98)00093-6.

    Article  CAS  PubMed  Google Scholar 

  4. Bosher JM, Labouesse M: RNA interference: genetic wand and genetic watchdog. Nat Cell Biol. 2000, 2 (2): E31-E36. 10.1038/35000102.

    Article  CAS  PubMed  Google Scholar 

  5. Lu C, Jeong DH, Kulkarni K, Pillay M, Nobuta K, German R, Thatcher SR, Maher C, Zhang L, Ware D, et al: Genome-wide analysis for discovery of rice microRNAs reveals natural antisense microRNAs (nat-miRNAs). Proc Natl Acad Sci USA. 2008, 105 (12): 4951-4956. 10.1073/pnas.0708743105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Tufarelli C, Stanley JAS, Garrick D, Sharpe JA, Ayyub H, Wood WG, Higgs DR: Transcription of antisense RNA leading to gene silencing and methylation as a novel cause of human genetic disease. Nature Genetics. 2003, 34 (2): 157-165. 10.1038/ng1157.

    Article  CAS  PubMed  Google Scholar 

  7. Reik W, Walter J: Genomic imprinting: Parental influence on the genome. Nat Rev Genet. 2001, 2 (1): 21-32. 10.1038/35047554.

    Article  CAS  PubMed  Google Scholar 

  8. Misra S, Crosby M, Mungall C, Matthews B, Campbell K, Hradecky P, Huang Y, Kaminker J, Millburn G, Prochnik S, et al: Annotation of the Drosophila melanogaster euchromatic genome: a systematic review. Genome Biol. 2002, 3 (12): research0083.0081-0083.0022. 10.1186/gb-2002-3-12-research0083.

    Article  Google Scholar 

  9. Chen JJ, Sun M, Kent WJ, Huang XQ, Xie HQ, Wang WQ, Zhou GL, Shi RZ, Rowley JD: Over 20% of human transcripts might form sense-antisense pairs. Nucleic Acids Research. 2004, 32 (16): 4812-4820. 10.1093/nar/gkh818.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Osato N, Yamada H, Satoh K, Ooka H, Yamamoto M, Suzuki K, Kawai J, Carninci P, Ohtomo Y, Murakami K, et al: Antisense transcripts with rice full-length cDNAs. Genome Biol. 2004, 5 (1): R5-10.1186/gb-2003-5-1-r5.

    Article  PubMed Central  Google Scholar 

  11. Gibbings JG, Cook BP, Dufault MR, Madden SL, Khuri S, Turnbull CJ, Dunwell JM: Global transcript analysis of rice leaf and seed using SAGE technology. Plant Biotechnology Journal. 2003, 1 (4): 271-285. 10.1046/j.1467-7652.2003.00026.x.

    Article  CAS  PubMed  Google Scholar 

  12. Meyers BC, Vu TH, Tej SS, Ghazal H, Matvienko M, Agrawal V, Ning JC, Haudenschild CD: Analysis of the transcriptional complexity of Arabidopsis thaliana by massively parallel signature sequencing. Nat Biotechnol. 2004, 22 (8): 1006-1011. 10.1038/nbt992.

    Article  CAS  PubMed  Google Scholar 

  13. Yamada K, Lim J, Dale JM, Chen HM, Shinn P, Palm CJ, Southwick AM, Wu HC, Kim C, Nguyen M, et al: Empirical analysis of transcriptional activity in the Arabidopsis genome. Science. 2003, 302 (5646): 842-846. 10.1126/science.1088305.

    Article  CAS  PubMed  Google Scholar 

  14. Poole R, Barker G, Werner K, Biggi G, Coghill J, Gibbings J, Berry S, Dunwell J, Edwards K: Analysis of wheat SAGE tags reveals evidence for widespread antisense transcription. BMC Genomics. 2008, 9: 475-10.1186/1471-2164-9-475.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Yelin R, Dahary D, Sorek R, Levanon EY, Goldstein O, Shoshan A, Diber A, Biton S, Tamir Y, Khosravi R, et al: Widespread occurrence of antisense transcription in the human genome. Nat Biotechnol. 2003, 21 (4): 379-386. 10.1038/nbt808.

    Article  CAS  PubMed  Google Scholar 

  16. Kiyosawa H, Mise N, Iwase S, Hayashizaki Y, Abe K: Disclosing hidden transcripts: Mouse natural sense-antisense transcripts tend to be poly(A) negative and nuclear localized. Genome Research. 2005, 15 (4): 463-474. 10.1101/gr.3155905.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Werner A, Schmutzler G, Carlile M, Miles CG, Peters H: Expression profiling of antisense transcripts on DNA arrays. Physiological Genomics. 2007, 28: 294-300. 10.1152/physiolgenomics.00127.2006.

    Article  CAS  PubMed  Google Scholar 

  18. Ge XJ, Rubinstein WS, Jung YC, Wu QF: Genome-wide analysis of antisense transcription with Affymetrix exon array. BMC Genomics. 2008, 9: 27-10.1186/1471-2164-9-27.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Coram TE, Settles ML, Chen XM: Transcriptome analysis of high-temperature adult-plant resistance conditioned by Yr39 during the wheat-Puccinia striiformis f. sp tritici interaction. Molecular Plant Pathology. 2008, 9 (4): 479-493. 10.1111/j.1364-3703.2008.00476.x.

    Article  CAS  PubMed  Google Scholar 

  20. Coram TE, Wang MN, Chen XM: Transcriptome analysis of the wheat – Puccinia striiformis f. sp tritici interaction. Molecular Plant Pathology. 2008, 9 (2): 157-169. 10.1111/j.1364-3703.2007.00453.x.

    Article  CAS  PubMed  Google Scholar 

  21. Coram T, Settles M, Wang M, Chen X: Surveying expression level polymorphism and single-feature polymorphism in near-isogenic wheat lines differing for the Yr5 stripe rust resistance locus. Theoretical and Applied Genetics. 2008, 117: 401-411. 10.1007/s00122-008-0784-5.

    Article  CAS  PubMed  Google Scholar 

  22. Parman C, Halling C, Gentleman R: affyQCReport:QC Report Generation for affyBatch objects. R package version 1.16.0. [http://www.bioconductor.org/packages/2.4/bioc/manuals/affyQCReport/man/affyQCReport.pdf]

  23. Gentleman R, Carey V, Bates D, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: Open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Kramer C, Loros J, Dunlap J, Crosthwaite S: Role for antisense RNA in regulating circadian clock function in Neurospora crassa. Nature. 2004, 421: 948-952. 10.1038/nature01427.

    Article  Google Scholar 

  25. Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10: 57-63. 10.1038/nrg2484.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Gautier L, Cope L, Bolstad B, Irizarry R: affy-analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315. 10.1093/bioinformatics/btg405.

    Article  CAS  PubMed  Google Scholar 

  27. Bolstad B, Irizarry R, Astrand M, Speed T: A Comparison of Normalization Methods for High Density Oligonucleotide Array Data Based on Bias and Variance. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.

    Article  CAS  PubMed  Google Scholar 

  28. Irizarry R, Bolstad B, Collin F, Cope L, Hobbs B, Speed T: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995, 57: 289-300.

    Google Scholar 

  30. Wise RP, Caldo RA, Hong L, Shen L, Cannon E, Dickerson JA: BarleyBase/PLEXdb – A unified expression profiling database for plants and plant pathogens. Methods in Molecular Biology. 2007, 347-363.

    Google Scholar 

  31. Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Haddad F, Qin AQX, Giger JM, Guo HY, Baldwin KM: Potential pitfalls in the accuracy of analysis of natural sense-antisense RNA pairs by reverse transcription-PCR. BMC Biotechnology. 2007, 7: 14-10.1186/1472-6750-7-21.

    Article  Google Scholar 

Download references

Acknowledgements

This research was supported in part by the US Department of Agriculture (USDA), Agricultural Research Service (ARS) (Project No. 5348-22000-014-00D), USDA-ARS Postdoctoral Program, and Washington Wheat Commission (Project No. 13C-3061-3923). PPNS No. 0510, Department of Plant Pathology, College of Agricultural, Human, and Natural Resources Research Center, Project Numbers: WNP00823. The authors acknowledge Derek Pouchnik (Washington State University, School of Molecular Biosciences) for carrying out GeneChip hybridizations and scanning, and thank Dr. Scot Hulbert and Dr. Patricia Okubara for critical reviews of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Tristan E Coram.

Additional information

Authors' contributions

TEC carried out the tissue collection, RNA extraction, RT-PCR, data analysis, and drafted the manuscript. MLS participated in the design of the study and performed the statistical analysis. TEC and XC conceived the study. XC participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2008_2137_MOESM1_ESM.xls

Additional file 1: List of the 110 sense-antisense transcript pairs with Present calls (P< 0.01), >8.76 3'IVT assay log2 expression value and >9.86 WT assay log2 expression value. The data provided represents annotation and statistical values for the candidate sense-antisense transcript pairs. (XLS 38 KB)

12864_2008_2137_MOESM2_ESM.xls

Additional file 2: List of the 80 potential antisense-specific transcripts with Present calls (P< 0.01) and >9.86 WT assay log2 expression value. The data provided represents annotation and statistical values for the potential antisense-specific transcripts. (XLS 29 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Coram, T.E., Settles, M.L. & Chen, X. Large-scale analysis of antisense transcription in wheat using the Affymetrix GeneChip Wheat Genome Array. BMC Genomics 10, 253 (2009). https://doi.org/10.1186/1471-2164-10-253

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-253

Keywords