A multiplex two-color real-time PCR method for quality-controlled molecular diagnostic testing of FFPE samples

Factors that drive the increasing use of FFPE tissue in basic and translational cancer research

The decision to use 10% neutral buffered formalin fixed, paraffin embedded (FFPE) archival pathology material may be dictated by the cancer research question or analytical technique, or may be governed by national ethical, legal and social implications (ELSI), biobank, and sample availability and access policy. Biobanked samples of common tumors are likely to be available, but not all samples will be annotated with treatment and outcomes data and this may limit their application. Tumors that are rare or very small exist mostly in FFPE pathology archives. Pathology departments worldwide contain millions of FFPE archival samples, but there are challenges to availability. Pathology departments lack resources for retrieving materials for research or for having pathologists select precise areas in paraffin blocks, a critical quality control step. When samples must be sourced from several pathology departments, different fixation and tissue processing approaches create variability in quality. Researchers must decide what sample quality and quality tolerance fit their specific purpose and whether sample enrichment is required. Recent publications report variable success with techniques modified to examine all common species of molecular targets in FFPE samples. Rigorous quality management may be particularly important in sample preparation for next generation sequencing and for optimizing the quality of extracted proteins for proteomics studies. Unpredictable failures, including unpublished ones, likely are related to pre-analytical factors, unstable molecular targets, biological and clinical sampling factors associated with specific tissue types or suboptimal quality management of pathology archives. Reproducible results depend on adherence to pre-analytical phase standards for molecular in vitro diagnostic analyses for DNA, RNA and in particular, extracted proteins. With continuing adaptations of techniques for application to FFPE, the potential to acquire much larger numbers of FFPE samples and the greater convenience of using FFPE in assays for precision medicine, the choice of material in the future will become increasingly biased toward FFPE samples from pathology archives. Recognition that FFPE samples may harbor greater variation in quality than frozen samples for several reasons, including variations in fixation and tissue processing, requires that FFPE results be validated provided a cohort of frozen tissue samples is available.

A high-throughput method to detect RNA profiling by integration of RT-MLPA with next generation sequencing technology

RNA in formalin-fixed and paraffin-embedded (FFPE) tissues provides large amount of information indicating disease stages, histological tumor types and grades, as well as clinical outcomes. However, Detection of RNA expression levels in formalin-fixed and paraffin-embedded samples is extremely difficult due to poor RNA quality. Here we developed a high-throughput method, Reverse Transcription-Multiple Ligation-dependent Probe Sequencing (RT-MLPSeq), to determine expression levels of multiple transcripts in FFPE samples. By combining Reverse Transcription-Multiple Ligation-dependent Amplification method and next generation sequencing technology, RT-MLPSeq overcomes the limit of probe length in multiplex ligation-dependent probe amplification assay and thus could detect expression levels of transcripts without quantitative limitations. We proved that different RT-MLPSeq probes targeting on the same transcripts have highly consistent results and the starting RNA/cDNA input could be as little as 1 ng. RT-MLPSeq also presented consistent relative RNA levels of selected 13 genes with reverse transcription quantitative PCR. Finally, we demonstrated the application of the new RT-MLPSeq method by measuring the mRNA expression levels of 21 genes which can be used for accurate calculation of the breast cancer recurrence score – an index that has been widely used for managing breast cancer patients.

RNAseq analysis of bronchial epithelial cells to identify COPD-associated genes and SNPs

Background

There is a need for more powerful methods to identify low-effect SNPs that contribute to hereditary COPD pathogenesis. We hypothesized that SNPs contributing to COPD risk through cis-regulatory effects are enriched in genes comprised by bronchial epithelial cell (BEC) expression patterns associated with COPD.

Methods

To test this hypothesis, normal BEC specimens were obtained by bronchoscopy from 60 subjects: 30 subjects with COPD defined by spirometry (FEV1/FVC < 0.7, FEV1% < 80%), and 30 non-COPD controls. Targeted next generation sequencing was used to measure total and allele-specific expression of 35 genes in genome maintenance (GM) genes pathways linked to COPD pathogenesis, including seven TP53 and CEBP transcription factor family members. Shrinkage linear discriminant analysis (SLDA) was used to identify COPD-classification models. COPD GWAS were queried for putative cis-regulatory SNPs in the targeted genes.

Results

On a network basis, TP53 and CEBP transcription factor pathway gene pair network connections, including key DNA repair gene ERCC5, were significantly different in COPD subjects (e.g., Wilcoxon rank sum test for closeness, p-value = 5.0E-11). ERCC5 SNP rs4150275 association with chronic bronchitis was identified in a set of Lung Health Study (LHS) COPD GWAS SNPs restricted to those in putative regulatory regions within the targeted genes, and this association was validated in the COPDgene non-hispanic white (NHW) GWAS. ERCC5 SNP rs4150275 is linked (D’ = 1) to ERCC5 SNP rs17655 which displayed differential allelic expression (DAE) in BEC and is an expression quantitative trait locus (eQTL) in lung tissue (p = 3.2E-7). SNPs in linkage (D’ = 1) with rs17655 were predicted to alter miRNA binding (rs873601). A classifier model that comprised gene features CAT, CEBPG, GPX1, KEAP1, TP73, and XPA had pooled 10-fold cross-validation receiver operator characteristic area under the curve of 75.4% (95% CI: 66.3%–89.3%). The prevalence of DAE was higher than expected (p = 0.0023) in the classifier genes.

Conclusions

GM genes comprised by COPD-associated BEC expression patterns were enriched for SNPs with cis-regulatory function, including a putative cis-rSNP in ERCC5 that was associated with COPD risk. These findings support additional total and allele-specific expression analysis of gene pathways with high prior likelihood for involvement in COPD pathogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12890-018-0603-y) contains supplementary material, which is available to authorized users.

Multianalyte quantitative competitive PCR on optically encoded microspheres for an eight-gene panel related to prostate cancer

Nucleic acid-based tests have a profound impact in every medical discipline. Because multigene tests offer higher diagnostic accuracy and lower overall cost than single assays, they are especially useful for diseases, like prostate cancer, that present variability at the molecular level and diversity of available therapeutic interventions. We have developed a quantitative competitive PCR for an eight-gene panel, related to prostate cancer, that includes five genes of the human tissue kallikrein family (KLKs), prostate-specific membrane antigen (PSMA), prostate cancer antigen 3 (PCA3), and HPRT1 as a reference gene. Using PCR as a synthetic tool, a competitor was prepared for each target sequence containing the same primer binding sites as the target but differing in a short segment to enable discrimination by hybridization. The assay involves multiplex amplification of targets and competitors followed by a multiplex hybridization assay for the 16 amplification products. The assay was performed on optically encoded microspheres with oligonucleotide probes attached to their surface. The microspheres were analyzed rapidly (1 min) by flow cytometry. The signal ratio of the target and cognate competitor is a function of the target copy number in the sample prior to amplification. The multiplexing potential of the proposed method is much higher than real-time PCR and other end-point methods since there are 100 sets of commercially available microspheres.

Straightforward and sensitive RT-qPCR based gene expression analysis of FFPE samples

Fragmented RNA from formalin-fixed paraffin-embedded (FFPE) tissue is a known obstacle to gene expression analysis. In this study, the impact of RNA integrity, gene-specific reverse transcription and targeted cDNA preamplification was quantified in terms of reverse transcription polymerase chain reaction (RT-qPCR) sensitivity by measuring 48 protein coding genes on eight duplicate cultured cancer cell pellet FFPE samples and twenty cancer tissue FFPE samples. More intact RNA modestly increased gene detection sensitivity by 1.6 fold (earlier detection by 0.7 PCR cycles, 95% CI = 0.593-0.850). Application of gene-specific priming instead of whole transcriptome priming during reverse transcription further improved RT-qPCR sensitivity by a considerable 4.0 fold increase (earlier detection by 2.0 PCR cycles, 95% CI = 1.73-2.32). Targeted cDNA preamplification resulted in the strongest increase of RT-qPCR sensitivity and enabled earlier detection by an average of 172.4 fold (7.43 PCR cycles, 95% CI = 6.83-7.05). We conclude that gene-specific reverse transcription and targeted cDNA preamplification are adequate methods for accurate and sensitive RT-qPCR based gene expression analysis of FFPE material. The presented methods do not involve expensive or complex procedures and can be easily implemented in any routine RT-qPCR practice.

Improving the Prediction of Prostate Cancer Overall Survival by Supplementing Readily Available Clinical Data with Gene Expression Levels of IGFBP3 and F3 in Formalin-Fixed Paraffin Embedded Core Needle Biopsy Material

Background

A previously reported expression signature of three genes (IGFBP3, F3 and VGLL3) was shown to have potential prognostic value in estimating overall and cancer-specific survivals at diagnosis of prostate cancer in a pilot cohort study using freshly frozen Fine Needle Aspiration (FNA) samples.

Methods

We carried out a new cohort study with 241 prostate cancer patients diagnosed from 2004–2007 with a follow-up exceeding 6 years in order to verify the prognostic value of gene expression signature in formalin fixed paraffin embedded (FFPE) prostate core needle biopsy tissue samples. The cohort consisted of four patient groups with different survival times and death causes. A four multiplex one-step RT-qPCR test kit, designed and optimized for measuring the expression signature in FFPE core needle biopsy samples, was used. In archive FFPE biopsy samples the expression differences of two genes (IGFBP3 and F3) were measured. The survival time predictions using the current clinical parameters only, such as age at diagnosis, Gleason score, PSA value and tumor stage, and clinical parameters supplemented with the expression levels of IGFBP3 and F3, were compared.

Results

When combined with currently used clinical parameters, the gene expression levels of IGFBP3 and F3 are improving the prediction of survival time as compared to using clinical parameters alone.

Conclusion

The assessment of IGFBP3 and F3 gene expression levels in FFPE prostate cancer tissue would provide an improved survival prediction for prostate cancer patients at the time of diagnosis.

Control for stochastic sampling variation and qualitative sequencing error in next generation sequencing

BACKGROUND

Clinical implementation of Next-Generation Sequencing (NGS) is challenged by poor control for stochastic sampling, library preparation biases and qualitative sequencing error. To address these challenges we developed and tested two hypotheses.

METHODS

Hypothesis 1: Analytical variation in quantification is predicted by stochastic sampling effects at input of a) amplifiable nucleic acid target molecules into the library preparation, b) amplicons from library into sequencer, or c) both. We derived equations using Monte Carlo simulation to predict assay coefficient of variation (CV) based on these three working models and tested them against NGS data from specimens with well characterized molecule inputs and sequence counts prepared using competitive multiplex-PCR amplicon-based NGS library preparation method comprising synthetic internal standards (IS). Hypothesis 2: Frequencies of technically-derived qualitative sequencing errors (i.e., base substitution, insertion and deletion) observed at each base position in each target native template (NT) are concordant with those observed in respective competitive synthetic IS present in the same reaction. We measured error frequencies at each base position within amplicons from each of 30 target NT, then tested whether they correspond to those within the 30 respective IS.

RESULTS

For hypothesis 1, the Monte Carlo model derived from both sampling events best predicted CV and explained 74% of observed assay variance. For hypothesis 2, observed frequency and type of sequence variation at each base position within each IS was concordant with that observed in respective NTs (R² = 0.93).

CONCLUSION

In targeted NGS, synthetic competitive IS control for stochastic sampling at input of both target into library preparation and of target library product into sequencer, and control for qualitative errors generated during library preparation and sequencing. These controls enable accurate clinical diagnostic reporting of confidence limits and limit of detection for copy number measurement, and of frequency for each actionable mutation.

Operator dependent choice of prostate cancer biopsy has limited impact on a gene signature analysis for the highly expressed genes IGFBP3 and F3 in prostate cancer epithelial cells

Background

Predicting the prognosis of prostate cancer disease through gene expression analysis is receiving increasing interest. In many cases, such analyses are based on formalin-fixed, paraffin embedded (FFPE) core needle biopsy material on which Gleason grading for diagnosis has been conducted. Since each patient typically has multiple biopsy samples, and since Gleason grading is an operator dependent procedure known to be difficult, the impact of the operator's choice of biopsy was evaluated.

Methods

Multiple biopsy samples from 43 patients were evaluated using a previously reported gene signature of IGFBP3, F3 and VGLL3 with potential prognostic value in estimating overall survival at diagnosis of prostate cancer. A four multiplex one-step qRT-PCR test kit, designed and optimized for measuring the signature in FFPE core needle biopsy samples was used. Concordance of gene expression levels between primary and secondary Gleason tumor patterns, as well as benign tissue specimens, was analyzed.

Results

The gene expression levels of IGFBP3 and F3 in prostate cancer epithelial cell-containing tissue representing the primary and secondary Gleason patterns were high and consistent, while the low expressed VGLL3 showed more variation in its expression levels.

Conclusion

The assessment of IGFBP3 and F3 gene expression levels in prostate cancer tissue is independent of Gleason patterns, meaning that the impact of operator's choice of biopsy is low.