Stuart D Pepper

A comparison of massively parallel nucleotide sequencing with oligonucleotide microarrays for global transcription profiling

BackgroundRNA-Seq exploits the rapid generation of gigabases of sequence data by Massively Parallel Nucleotide Sequencing, allowing for the mapping and digital quantification of whole transcriptomes. Whilst previous comparisons between RNA-Seq and microarrays have been performed at the level of gene expression, in this study we adopt a more fine-grained approach. Using RNA samples from a normal human breast epithelial cell line (MCF-10a) and a breast cancer cell line (MCF-7), we present a comprehensive comparison between RNA-Seq data generated on the Applied Biosystems SOLiD platform and data from Affymetrix Exon 1.0ST arrays. The use of Exon arrays makes it possible to assess the performance of RNA-Seq in two key areas: detection of expression at the granularity of individual exons, and discovery of transcription outside annotated loci.ResultsWe found a high degree of correspondence between the two platforms in terms of exon-level fold changes and detection. For example, over 80% of exons detected as expressed in RNA-Seq were also detected on the Exon array, and 91% of exons flagged as changing from Absent to Present on at least one platform had fold-changes in the same direction. The greatest detection correspondence was seen when the read count threshold at which to flag exons Absent in the SOLiD data was set to t<1 suggesting that the background error rate is extremely low in RNA-Seq. We also found RNA-Seq more sensitive to detecting differentially expressed exons than the Exon array, reflecting the wider dynamic range achievable on the SOLiD platform. In addition, we find significant evidence of novel protein coding regions outside known exons, 93% of which map to Exon array probesets, and are able to infer the presence of thousands of novel transcripts through the detection of previously unreported exon-exon junctions.ConclusionsBy focusing on exon-level expression, we present the most fine-grained comparison between RNA-Seq and microarrays to date. Overall, our study demonstrates that data from a SOLiD RNA-Seq experiment are sufficient to generate results comparable to those produced from Affymetrix Exon arrays, even using only a single replicate from each platform, and when presented with a large genome.

The utility of MAS5 expression summary and detection call algorithms

BackgroundUsed alone, the MAS5.0 algorithm for generating expression summaries has been criticized for high False Positive rates resulting from exaggerated variance at low intensities.ResultsHere we show, with replicated cell line data, that, when used alongside detection calls, MAS5 can be both selective and sensitive. A set of differentially expressed transcripts were identified that were found to be changing by MAS5, but unchanging by RMA and GCRMA. Subsequent analysis by real time PCR confirmed these changes. In addition, with the Latin square datasets often used to assess expression summary algorithms, filtered MAS5.0 was found to have performance approaching that of its peers.ConclusionWhen used alongside detection calls, MAS5 is a sensitive and selective algorithm for identifying differentially expressed genes.

Amplification protocols introduce systematic but reproducible errors into gene expression studies

The desire to perform microarray experiments with small starting amounts of RNA has led to the development of a variety of protocols for preparing and amplifying mRNA. This has consequences not only for the standardization of experimental design, but also for reproducibility and comparability between experiments. Here we investigate the differences between the Affymetrix standard and small sample protocols and address the data analysis issues that arise when comparing samples and experiments that have been processed in different ways. We show that data generated on the same platform using different protocols are not directly comparable. Further, protocols introduce systematic biases that can be largely accounted for by using the correct data analysis techniques.

Acquisition of biologically relevant gene expression data by Affymetrix microarray analysis of archival formalin-fixed paraffin-embedded tumours.

Robust protocols for microarray gene expression profiling of archival formalin-fixed paraffin-embedded tissue (FFPET) are needed to facilitate research when availability of fresh-frozen tissue is limited. Recent reports attest to the feasibility of this approach, but the clinical value of these data is poorly understood. We employed state-of-the-art RNA extraction and Affymetrix microarray technology to examine 34 archival FFPET primary extremity soft tissue sarcomas. Nineteen arrays met stringent QC criteria and were used to model prognostic signatures for metastatic recurrence. Arrays from two paired frozen and FFPET samples were compared: although FFPET sensitivity was low (∼50%), high specificity (95%) and positive predictive value (92%) suggest that transcript detection is reliable. Good agreement between arrays and real time (RT)–PCR was confirmed, especially for abundant transcripts, and RT–PCR validated the regulation pattern for 19 of 24 candidate genes (overall R2=0.4662). RT–PCR and immunohistochemistry on independent cases validated prognostic significance for several genes including RECQL4, FRRS1, CFH and MET – whose combined expression carried greater prognostic value than tumour grade – and cmet and TRKB proteins. These molecules warrant further evaluation in larger series. Reliable clinically relevant data can be obtained from archival FFPET, but protocol amendments are needed to improve the sensitivity and broad application of this approach.

High correspondence between Affymetrix exon and standard expression arrays.

Exon arrays aim to provide comprehensive gene expression data at the level of individual exons, similar to that provided on a per-gene basis by existing expression arrays. This report describes the performance of Affymetrix GeneChip Human Exon 1.0 ST array by using replicated RNA samples from two human cell lines, MCF7 and MCF10A, hybridized both to Exon 1.0 ST and to HG-U133 Plus2 arrays. Cross-comparison between array types requires an appropriate mapping to be found between individual probe sets. Three possible mappings were considered, reflecting different strategies for dealing with probe sets that target different parts of the same transcript. Irrespective of the mapping used, Exon 1.0 ST and HG-U133 Plus2 arrays show a high degree of correspondence. More than 80% of HG-U133 Plus2 probe sets may be mapped to the Exon chip, and fold changes are found well preserved for over 96% of those probe sets detected present. Since HG-U133 Plus2 arrays have already been extensively validated, these results lend a significant degree of confidence to exon arrays.

Methods comparison for high-resolution transcriptional analysis of archival material on Affymetrix Plus 2.0 and Exon 1.0 microarrays.

Microarray gene expression profiling of formalin-fixed paraffin-embedded (FFPE) tissues is a new and evolving technique. This report compares transcript detection rates on Affymetrix U133 Plus 2.0 and Human Exon 1.0 ST GeneChips across several RNA extraction and target labeling protocols, using routinely collected archival FFPE samples. All RNA extraction protocols tested (Ambion-Optimum, Ambion-RecoverAll, and Qiagen-RNeasy FFPE) provided extracts suitable for microarray hybridization. Compared with Affymetrix One-Cycle labeled extracts, NuGEN system protocols utilizing oligo(dT) and random hexamer primers, and cDNA target preparations instead of cRNA, achieved percent present rates up to 55% on Plus 2.0 arrays. Based on two paired-sample analyses, at 90% specificity this equalled an average 30 percentage-point increase (from 50% to 80%) in FFPE transcript sensitivity relative to fresh frozen tissues, which we have assumed to have 100% sensitivity and specificity. The high content of Exon arrays, with multiple probe sets per exon, improved FFPE sensitivity to 92% at 96% specificity, corresponding to an absolute increase of ~600 genes over Plus 2.0 arrays. While larger series are needed to confirm high correspondence between fresh-frozen and FFPE expression patterns, these data suggest that both Plus 2.0 and Exon arrays are suitable platforms for FFPE microarray expression analyses.

Programmed fluctuations in sense/antisense transcript ratios drive sexual differentiation in S. pombe.

Strand‐specific RNA sequencing of S. pombe revealed a highly structured programme of ncRNA expression at over 600 loci. Waves of antisense transcription accompanied sexual differentiation. A substantial proportion of ncRNA arose from mechanisms previously considered to be largely artefactual, including improper 3′ termination and bidirectional transcription. Constitutive induction of the entire spk1+, spo4+, dis1+ and spo6+ antisense transcripts from an integrated, ectopic, locus disrupted their respective meiotic functions. This ability of antisense transcripts to disrupt gene function when expressed in trans suggests that cis production at native loci during sexual differentiation may also control gene function. Consistently, insertion of a marker gene adjacent to the dis1+ antisense start site mimicked ectopic antisense expression in reducing the levels of this microtubule regulator and abolishing the microtubule‐dependent ‘horsetail’ stage of meiosis. Antisense production had no impact at any of these loci when the RNA interference (RNAi) machinery was removed. Thus, far from being simply ‘genome chatter’, this extensive ncRNA landscape constitutes a fundamental component in the controls that drive the complex programme of sexual differentiation in S. pombe.

Differences in serological IgA responses to recombinant baculovirus-derived human papillomavirus E2 protein in the natural history of cervical neoplasia

Infection with certain types of human papillomavirus (HPV) presents a high risk for the subsequent development of cervical intraepithelial neoplasia (CIN) and cervical carcinoma. Immunological mechanisms are likely to play a role in control of cervical HPV lesions. The HPV E2 protein has roles in virus replication and transcription, and loss of E2 functions may be associated with progression of cervical neoplasia. Accordingly, it is of interest to monitor immune responses to the E2 protein, and previous studies have reported associations between serological reactivity to E2 peptide antigens and cervical neoplasia. In order to investigate serological responses to native, full-length E2 protein, we expressed HPV-16 E2 proteins with and without an N-terminal polyhistidine tag using the baculovirus system. Purified HPV-16 E2 protein was used to develop enzyme-linked immunosorbent assays to detect serological IgG and IgA responses in cervical neoplasia patients and controls. We found that serum IgA levels against the E2 protein were elevated in CIN patients relative to normal control subjects but were not elevated in cervical cancer patients. Moreover, there appeared to be a gradient of response within cervical neoplasia such that the highest antibody levels were seen in lower grades of neoplasia up to CIN 2, whereas lower levels were observed in CIN 3 and still lower levels in cervical carcinoma. These findings suggest that the IgA antibody response to E2 may associate with stage and progression in cervical neoplasia.

Exon-array profiling unlocks clinically and biologically relevant gene signatures from formalin-fixed paraffin-embedded tumour samples.

Background:Degradation and chemical modification of RNA in formalin-fixed paraffin-embedded (FFPE) samples hamper their use in expression profiling studies. This study aimed to show that useful information can be obtained by Exon-array profiling archival FFPE tumour samples.Methods:Nineteen cervical squamous cell carcinoma (SCC) and 9 adenocarcinoma (AC) FFPE samples (10–16-year-old) were profiled using Affymetrix Exon arrays. The gene signature derived was tested on a fresh-frozen non-small cell lung cancer (NSCLC) series. Exploration of biological networks involved gene set enrichment analysis (GSEA). Differential gene expression was confirmed using Quantigene, a multiplex bead-based alternative to qRT–PCR.Results:In all, 1062 genes were higher in SCC vs AC, and 155 genes higher in AC. The 1217-gene signature correctly separated 58 NSCLC into SCC and AC. A gene network centered on hepatic nuclear factor and GATA6 was identified in AC, suggesting a role in glandular cell differentiation of the cervix. Quantigene analysis of the top 26 differentially expressed genes correctly partitioned cervix samples as SCC or AC.Conclusion:FFPE samples can be profiled using Exon arrays to derive gene expression signatures that are sufficiently robust to be applied to independent data sets, identify novel biology and design assays for independent platform validation.

Immunisation of common marmosets with vaccinia virus expressing Epstein-Barr virus (EBV) gp340 and challenge with EBV

Epstein‐Barr virus (EBV) is the cause of infectious mononucleosis and is associated with a variety of life‐threatening diseases in humans. Therefore the development of an effective vaccine is an important objective. Many of the initial studies of vaccine efficacy analyse the ability of vaccine preparations to prevent the induction of lymphomas in cottontop tamarins by the B95‐8 strain of EBV. We used a vaccinia virus recombinant expressing gp340, vMA1, tested previously in the cotton‐top tamarin, to evaluate a common marmoset model in which the challenge virus, M81, resembles more closely the wild‐type strains of EBV in the general population than does the standard B95‐8 strain. We characterised the M81 strain of EBV with respect to the sequence of its gp340/220 gene and in regard to the presence of a region deleted in B95‐8. Replication of the challenge virus in the group vaccinated with vMA1 was decreased when compared to control groups.