Whole-genome gene expression profiling of formalin-fixed, paraffin-embedded tissue samples

Application of microarray in breast cancer: An overview

There are more than 1.15 million cases of breast cancer diagnosed worldwide annually. At present, only small numbers of accurate prognostic and predictive factors are used clinically for managing the patients with breast cancer. DNA microarrays have the potential to assess the expression of thousands of genes simultaneously. Recent preliminary researches indicate that gene expression profiling based on DNA microarray can offer potential and independent prognostic information in patients with newly diagnosed breast cancer. In this paper, an overview upon the applications of microarray techniques in breast cancer is presented.

Discovery of biomarkers of endometrial receptivity through a minimally invasive approach: a validation study with implications for assisted reproduction

OBJECTIVE

To determine whether a minimally invasive approach to sampling endometrial cells that can be applied during an active conception cycle can generate robust biomarker candidates for endometrial receptivity by genomewide gene expression profiling.

DESIGN

Longitudinal study comparing gene expression profiles of cells isolated from uterine aspirates collected during the prereceptive and receptive phases of a natural cycle.

SETTING

University-affiliated hospital.

PATIENT(S)

Healthy volunteers, ≤40 years of age, with regular menstrual cycles and no history of infertility.

INTERVENTION(S)

One menstrual cycle monitored with urinary kits to identify the luteinizing hormone (LH) surge; uterine aspirations collected at LH + 2 days (LH + 2) and at LH + 7; endometrial biopsy obtained on LH + 7; RNA extraction from the cellular material for gene expression profiling, and differential gene expression validated by NanoString assay and cross-validated against a publically available data set.

MAIN OUTCOME MEASURE(S)

Differentially expressed genes between LH + 2 and LH + 7 samples.

RESULT(S)

NanoString assay validated 96% of the 245 genes found differentially expressed at LH + 7. Unsupervised hierarchical clustering of aspiration and biopsy samples demonstrated the concordance of the sampling methods. A predictor gene cassette derived by a shrunken centroid class prediction technique correctly classified the receptive phase within an external data set.

CONCLUSION(S)

Uterine aspiration, which can be performed during an active conception cycle, identified robust candidate biomarkers of endometrial receptivity, and will enable their validation by direct correlation with clinical outcomes.

The transcriptional architecture of early human hematopoiesis identifies multilevel control of lymphoid commitment

Understanding how differentiation programs originate from within the gene expression landscape of hematopoietic stem cells (HSC) is crucial to develop new clinical therapies. We mapped the transcriptional dynamics underlying the first steps of commitment by tracking transcriptome changes in human HSC and eight early progenitor populations. Transcriptional programs are extensively shared, extend across lineage-potential boundaries, and are not strictly lineage-affiliated. Elements of stem, lymphoid and myeloid programs are retained in multi-lymphoid progenitors (MLP), reflecting a hybrid transcriptional state. Based on functional single cell analysis, BCL11A, SOX4 and TEAD1 governed transcriptional networks within MLPs, leading to B cell specification. Overall, we show that integrated transcriptome approaches can identify novel regulators of multipotency and uncover additional complexity in lymphoid commitment.

An optimized workflow for improved gene expression profiling for formalin-fixed, paraffin-embedded tumor samples

BACKGROUND

Whole genome microarray gene expression profiling is the 'gold standard' for the discovery of prognostic and predictive genetic markers for human cancers. However, suitable research material is lacking as most diagnostic samples are preserved as formalin-fixed, paraffin-embedded tissue (FFPET). We tested a new workflow and data analysis method optimized for use with FFPET samples.

METHODS

Sixteen breast tumor samples were split into matched pairs and preserved as FFPET or fresh-frozen (FF). Total RNA was extracted and tested for yield and purity. RNA from FFPET samples was amplified using three different commercially available kits in parallel, and hybridized to Affymetrix GeneChip® Human Genome U133 Plus 2.0 Arrays. The array probe set was optimized in silico to exclude misdesigned and misannotated probes.

RESULTS

FFPET samples processed using the WT-Ovation™ FFPE System V2 (NuGEN) provided 80% specificity and 97% sensitivity compared with FF samples (assuming values of 100%). In addition, in silico probe set redesign improved sequence detection sensitivity and, thus, may rescue potentially significant small-magnitude gene expression changes that could otherwise be diluted by the overall probe set background.

CONCLUSION

In conclusion, our FFPET-optimized workflow enables the detection of more genes than previous, nonoptimized approaches, opening new possibilities for the discovery, validation, and clinical application of mRNA biomarkers in human diseases.

Prognostic gene expression signature for patients with hepatitis C-related early-stage cirrhosis

BACKGROUND & AIMS

Cirrhosis affects 1% to 2% of the world population and is the major risk factor for hepatocellular carcinoma (HCC). Hepatitis C cirrhosis-related HCC is the most rapidly increasing cause of cancer death in the United States. Noninvasive methods have been developed to identify patients with asymptomatic early-stage cirrhosis, increasing the burden of HCC surveillance, but biomarkers are needed to identify patients with cirrhosis who are most in need of surveillance. We investigated whether a liver-derived 186-gene signature previously associated with outcomes of patients with HCC is prognostic for patients with newly diagnosed cirrhosis but without HCC.

METHODS

We performed gene expression profile analysis of formalin-fixed needle biopsy specimens from the livers of 216 patients with hepatitis C-related early-stage (Child-Pugh class A) cirrhosis who were prospectively followed up for a median of 10 years at an Italian center. We evaluated whether the 186-gene signature was associated with death, progression of cirrhosis, and development of HCC.

RESULTS

Fifty-five (25%), 101 (47%), and 60 (28%) patients were classified as having poor-, intermediate-, and good-prognosis signatures, respectively. In multivariable Cox regression modeling, the poor-prognosis signature was significantly associated with death (P = .004), progression to advanced cirrhosis (P < .001), and development of HCC (P = .009). The 10-year rates of survival were 63%, 74%, and 85% and the annual incidence of HCC was 5.8%, 2.2%, and 1.5% for patients with poor-, intermediate-, and good-prognosis signatures, respectively.

CONCLUSIONS

A 186-gene signature used to predict outcomes of patients with HCC is also associated with outcomes of patients with hepatitis C-related early-stage cirrhosis. This signature might be used to identify patients with cirrhosis in most need of surveillance and strategies to prevent the development of HCC.

High-throughput RNA sequencing of a formalin-fixed, paraffin-embedded autopsy lung tissue sample from the 1918 influenza pandemic

Most biopsy and autopsy tissues are formalin-fixed and paraffin-embedded (FFPE), but this process leads to RNA degradation that limits gene expression analysis. The RNA genome of the 1918 pandemic influenza virus was previously determined in a 9-year effort by overlapping RT-PCR from post-mortem samples. Here, the full genome of the 1918 virus at 3000× coverage was determined in one high-throughput sequencing run of a library derived from total RNA of a 1918 FFPE sample after duplex-specific nuclease treatments. Bacterial sequences associated with secondary bacterial pneumonias were also detected. Host transcripts were well represented in the library. Compared to a 2009 pandemic influenza virus FFPE post-mortem library, the 1918 sample showed significant enrichment for host defence and cell death response genes, concordant with prior animal studies. This methodological approach should assist in the analysis of FFPE tissue samples isolated over the past century from a variety of diseases.

Quality assessment metrics for whole genome gene expression profiling of paraffin embedded samples

BACKGROUND

Formalin fixed, paraffin embedded tissues are most commonly used for routine pathology analysis and for long term tissue preservation in the clinical setting. Many institutions have large archives of Formalin fixed, paraffin embedded tissues that provide a unique opportunity for understanding genomic signatures of disease. However, genome-wide expression profiling of Formalin fixed, paraffin embedded samples have been challenging due to RNA degradation. Because of the significant heterogeneity in tissue quality, normalization and analysis of these data presents particular challenges. The distribution of intensity values from archival tissues are inherently noisy and skewed due to differential sample degradation raising two primary concerns; whether a highly skewed array will unduly influence initial normalization of the data and whether outlier arrays can be reliably identified.

FINDINGS

Two simple extensions of common regression diagnostic measures are introduced that measure the stress an array undergoes during normalization and how much a given array deviates from the remaining arrays post-normalization. These metrics are applied to a study involving 1618 formalin-fixed, paraffin-embedded HER2-positive breast cancer samples from the N9831 adjuvant trial processed with Illumina's cDNA-mediated Annealing Selection extension and Ligation assay.

CONCLUSION

Proper assessment of array quality within a research study is crucial for controlling unwanted variability in the data. The metrics proposed in this paper have direct biological interpretations and can be used to identify arrays that should either be removed from analysis all together or down-weighted to reduce their influence in downstream analyses.

MicroRNA paraffin-based studies in osteosarcoma reveal reproducible independent prognostic profiles at 14q32

Background

Although microRNAs (miRNAs) are implicated in osteosarcoma biology and chemoresponse, miRNA prognostic models are still needed, particularly because prognosis is imperfectly correlated with chemoresponse. Formalin-fixed, paraffin-embedded tissue is a necessary resource for biomarker studies in this malignancy with limited frozen tissue availability.

Methods

We performed miRNA and mRNA microarray formalin-fixed, paraffin-embedded assays in 65 osteosarcoma biopsy and 26 paired post-chemotherapy resection specimens and used the only publicly available miRNA dataset, generated independently by another group, to externally validate our strongest findings (n = 29). We used supervised principal components analysis and logistic regression for survival and chemoresponse, and miRNA activity and target gene set analysis to study miRNA regulatory activity.

Results

Several miRNA-based models with as few as five miRNAs were prognostic independently of pathologically assessed chemoresponse (median recurrence-free survival: 59 months versus not-yet-reached; adjusted hazards ratio = 2.90; P = 0.036). The independent dataset supported the reproducibility of recurrence and survival findings. The prognostic value of the profile was independent of confounding by known prognostic variables, including chemoresponse, tumor location and metastasis at diagnosis. Model performance improved when chemoresponse was added as a covariate (median recurrence-free survival: 59 months versus not-yet-reached; hazard ratio = 3.91; P = 0.002). Most prognostic miRNAs were located at 14q32 - a locus already linked to osteosarcoma - and their gene targets display deregulation patterns associated with outcome. We also identified miRNA profiles predictive of chemoresponse (75% to 80% accuracy), which did not overlap with prognostic profiles.

Conclusions

Formalin-fixed, paraffin-embedded tissue-derived miRNA patterns are a powerful prognostic tool for risk-stratified osteosarcoma management strategies. Combined miRNA and mRNA analysis supports a possible role of the 14q32 locus in osteosarcoma progression and outcome. Our study creates a paradigm for formalin-fixed, paraffin-embedded-based miRNA biomarker studies in cancer.

Quality assessment and data handling methods for Affymetrix Gene 1.0 ST arrays with variable RNA integrity

Background

RNA and microarray quality assessment form an integral part of gene expression analysis and, although methods such as the RNA integrity number (RIN) algorithm reliably asses RNA integrity, the relevance of RNA integrity in gene expression analysis as well as analysis methods to accommodate the possible effects of degradation requires further investigation. We investigated the relationship between RNA integrity and array quality on the commonly used Affymetrix Gene 1.0 ST array platform using reliable within-array and between-array quality assessment measures. The possibility of a transcript specific bias in the apparent effect of RNA degradation on the measured gene expression signal was evaluated after either excluding quality-flagged arrays or compensation for RNA degradation at different steps in the analysis.

Results

Using probe-level and inter-array quality metrics to assess 34 Gene 1.0 ST array datasets derived from historical, paired tumour and normal primary colorectal cancer samples, 7 arrays (20.6%), with a mean sample RIN of 3.2 (SD = 0.42), were flagged during array quality assessment while 10 arrays from samples with RINs < 7 passed quality assessment, including one sample with a RIN < 3. We detected a transcript length bias in RNA degradation in only 5.8% of annotated transcript clusters (p-value 0.05, FC ≥ |2|), with longer and shorter than average transcripts under- and overrepresented in quality-flagged samples respectively. Applying compensatory measures for RNA degradation performed at least as well as excluding quality-flagged arrays, as judged by hierarchical clustering, gene expression analysis and Ingenuity Pathway Analysis; importantly, use of these compensatory measures had the significant benefit of enabling lower quality array data from irreplaceable clinical samples to be retained in downstream analyses.

Conclusions

Here, we demonstrate an effective array-quality assessment strategy, which will allow the user to recognize lower quality arrays that can be included in the analysis once appropriate measures are applied to account for known or unknown sources of variation, such as array quality- and batch- effects, by implementing ComBat or Surrogate Variable Analysis. This approach of quality control and analysis will be especially useful for clinical samples with variable and low RNA qualities, with RIN scores ≥ 2.

Immune markers and differential signaling networks in ulcerative colitis and Crohn's disease

BACKGROUND

Cytokine signaling pathways play a central role in the pathogenesis of inflammatory bowel disease (IBD). Ulcerative colitis (UC) and Crohn's disease (CD) have unique as well as overlapping phenotypes, susceptibility genes, and gene expression profiles. This study aimed to delineate patterns within cytokine signaling pathways in colonic mucosa of UC and CD patients, explore molecular diagnostic markers, and identify novel immune mediators in IBD pathogenesis.

METHODS

We quantified 70 selected immune genes that are important in IBD signaling from formalin-fixed, paraffin-embedded (FFPE) colon biopsy samples from normal control subjects and UC and CD patients having either severe colitis or quiescent disease (n = 98 subjects). We utilized and validated a new modified real-time reverse-transcription polymerase chain reaction (RT-PCR) technique for gene quantification.

RESULTS

Expression levels of signaling molecules including IL-6/10/12/13/17/23/33, STAT1/3/6, T-bet, GATA3, Foxp3, SOCS1/3, and downstream inflammatory mediators such as chemokines CCL-2/11/17/20, oxidative stress inducers, proteases, and mucosal genes were differentially regulated between UC and CD and between active and quiescent disease. We also document the possible role of novel genes in IBD, including SHP-1, IRF-1,TARC, Eotaxin, NOX2, arginase I, and ADAM 8.

CONCLUSIONS

This comprehensive approach to quantifying gene expression provides insights into the pathogenesis of IBD by elucidating distinct immune signaling networks in CD and UC. Furthermore, this is the first study demonstrating that gene expression profiling in FFPE colon biopsies might be a practical and effective tool in the diagnosis and prognosis of IBD and may help identify molecular markers that can predict and monitor response to individualized therapeutic treatments.

Expression profiling of archival tumors for long-term health studies

PURPOSE

More than 20 million archival tissue samples are stored annually in the United States as formalin-fixed, paraffin-embedded (FFPE) blocks, but RNA degradation during fixation and storage has prevented their use for transcriptional profiling. New and highly sensitive assays for whole-transcriptome microarray analysis of FFPE tissues are now available, but resulting data include noise and variability for which previous expression array methods are inadequate.

EXPERIMENTAL DESIGN

We present the two largest whole-genome expression studies from FFPE tissues to date, comprising 1,003 colorectal cancer (CRC) and 168 breast cancer samples, combined with a meta-analysis of 14 new and published FFPE microarray datasets. We develop and validate quality control (QC) methods through technical replication, independent samples, comparison to results from fresh-frozen tissue, and recovery of expected associations between gene expression and protein abundance.

RESULTS

Archival tissues from large, multicenter studies showed a much wider range of transcriptional data quality relative to smaller or frozen tissue studies and required stringent QC for subsequent analysis. We developed novel methods for such QC of archival tissue expression profiles based on sample dynamic range and per-study median profile. This enabled validated identification of gene signatures of microsatellite instability and additional features of CRC, and improved recovery of associations between gene expression and protein abundance of MLH1, FASN, CDX2, MGMT, and SIRT1 in CRC tumors.

CONCLUSIONS

These methods for large-scale QC of FFPE expression profiles enable study of the cancer transcriptome in relation to extensive clinicopathological information, tumor molecular biomarkers, and long-term lifestyle and outcome data.

Juxtaglomerular cell tumor: a morphological, immunohistochemical and genetic study of six cases

Juxtaglomerular cell tumors (JGCTs) are rare tumors characterized by renin synthesis, hyperaldosteronism and hypertension. A curious immunohistochemical overlap between JGCT and gastrointestinal stromal tumor (GIST) including the expression of vimentin, CD34, CD117, α-smooth muscle actin was previously reported, prompting us to further investigate JGCT and its phenotypic and molecular genetic characteristics. Virtual karyotyping showed gain of chromosomes 3, 4, 10, 13, 17 and 18 in one JGCT, and fluorescence in situ hybridization (FISH) study confirmed this multiple gain pattern. Additionally, loss of chromosome 9 was observed in four of six cases analyzed with FISH. A whole genome expression analysis revealed 415 up-regulated (including renin, and CD117) and 325 down-regulated genes between the 2 cases. The study confirmed earlier reports on the gain of chromosomes 4 and 10, and provided further evidence of up-regulation of the genes located on these 2 chromosomes. For the first time our study indicated the importance of the loss of chromosome 9 and loss of expression of several tumor suppressor genes located on this chromosome as possible pathogenetic events important in development of JGCT.

Identification and validation of gene expression models that predict clinical outcome in patients with early-stage laryngeal cancer

BACKGROUND

Despite improvement in therapeutic techniques, patients with early-stage laryngeal cancer still recur after treatment. Gene expression prognostic models could suggest which of these patients would be more appropriate for testing adjuvant strategies.

MATERIALS AND METHODS

Expression profiling using whole-genome DASL arrays was carried out on 56 formalin-fixed paraffin-embedded tumor samples of patients with early-stage laryngeal cancer. We split the samples into a training and a validation set. Using the supervised principal components survival analysis in the first cohort, we identified gene expression profiles that predict the risk of recurrence. These profiles were then validated in an independent cohort.

RESULTS

Gene models comprising different number of genes identified a subgroup of patients who were at high risk of recurrence. Of these, the best prognostic model distinguished between a high- and a low-risk group (log-rank P<0.005). The prognostic value of this model was reproduced in the validation cohort (median disease-free survival: 38 versus 161 months, log-rank P=0.018), hazard ratio=5.19 (95% confidence interval 1.14-23.57, P<0.05).

CONCLUSIONS

We have identified gene expression prognostic models that can refine the estimation of a patient's risk of recurrence. These findings, if further validated, should aid in patient stratification for testing adjuvant treatment strategies.

A microRNA activity map of human mesenchymal tumors: connections to oncogenic pathways; an integrative transcriptomic study

Background

MicroRNAs (miRNAs) are nucleic acid regulators of many human mRNAs, and are associated with many tumorigenic processes. miRNA expression levels have been used in profiling studies, but some evidence suggests that expression levels do not fully capture miRNA regulatory activity. In this study we integrate multiple gene expression datasets to determine miRNA activity patterns associated with cancer phenotypes and oncogenic pathways in mesenchymal tumors – a very heterogeneous class of malignancies.

Results

Using a computational method, we identified differentially activated miRNAs between 77 normal tissue specimens and 135 sarcomas and we validated many of these findings with microarray interrogation of an independent, paraffin-based cohort of 18 tumors. We also showed that miRNA activity is imperfectly correlated with miRNA expression levels. Using next-generation miRNA sequencing we identified potential base sequence alterations which may explain differential activity. We then analyzed miRNA activity changes related to the RAS-pathway and found 21 miRNAs that switch from silenced to activated status in parallel with RAS activation. Importantly, nearly half of these 21 miRNAs were predicted to regulate integral parts of the miRNA processing machinery, and our gene expression analysis revealed significant reductions of these transcripts in RAS-active tumors. These results suggest an association between RAS signaling and miRNA processing in which miRNAs may attenuate their own biogenesis.

Conclusions

Our study represents the first gene expression-based investigation of miRNA regulatory activity in human sarcomas, and our findings indicate that miRNA activity patterns derived from integrated transcriptomic data are reproducible and biologically informative in cancer. We identified an association between RAS signaling and miRNA processing, and demonstrated sequence alterations as plausible causes for differential miRNA activity. Finally, our study highlights the value of systems level integrative miRNA/mRNA assessment with high-throughput genomic data, and the applicability of paraffin-tissue-derived RNA for validation of novel findings.

Formalin fixation in the '-omics' era: a primer for the surgeon-scientist

Formalin is the most commonly used tissue fixative worldwide. While it offers excellent morphological preservation for routine histology, it has detrimental effects on nucleic acids. Most studies of nucleic acids have therefore used fresh frozen tissue, the collection and storage of which is resource intensive. The ability to use modern genomic, transcriptomic and epigenomic methods with nucleic acids derived from formalin-fixed, paraffin-embedded (FFPE) tissues would allow enormous archives of routinely stored tissues (usually with well-annotated clinical data) to be used for translational research. This paper outlines the effects of formalin on nucleic acids, describes ways of minimizing nucleic acid degradation and optimizing extraction, and reviews recent studies that have used contemporary techniques to analyse FFPE-derived nucleic acids (with a focus on malignant tissue sources). Simple tips are also offered to ensure the utility of your institution's samples for future studies, and broadly applicable guidelines are listed for those contemplating their own study using FFPE-derived material.

Transcriptional profiling of formalin fixed paraffin embedded tissue: pitfalls and recommendations for identifying biologically relevant changes

Expression profiling techniques have been used to study the biology of many types of cancer but have been limited to some extent by the requirement for collection of fresh tissue. In contrast, formalin fixed paraffin embedded (FFPE) samples are widely available and represent a vast resource of potential material. The techniques used to handle the degraded and modified RNA from these samples are relatively new and all the pitfalls and limitations of this material for whole genome expression profiling are not yet clarified. Here, we analyzed 70 FFPE tongue carcinoma samples and 17 controls using the whole genome DASL array covering nearly 21000 genes. We identified that sample age is related to quality of extracted RNA and that sample quality influences apparent expression levels in a non-random manner related to gene probe sequence, leading to spurious results. However, by removing sub-standard samples and analysing only those 28 cancers and 15 controls that had similar quality we were able to generate a list of 934 genes significantly altered in tongue cancer compared to control samples of tongue. This list contained previously identified changes and was enriched for genes involved in many cancer-related processes such as tissue remodelling, inflammation, differentiation and apoptosis. Four novel genes of potential importance in tongue cancer development and maintenance, SH3BGL2, SLC2A6, SLC16A3 and CXCL10, were independently confirmed, validating our data. Hence, gene expression profiling can be performed usefully on archival material if appropriate quality assurance steps are taken to ensure sample consistency and we present some recommendations for the use of FFPE material based on our findings.

Effective DNA/RNA co-extraction for analysis of microRNAs, mRNAs, and genomic DNA from formalin-fixed paraffin-embedded specimens

Background

Retrospective studies of archived human specimens, with known clinical follow-up, are used to identify predictive and prognostic molecular markers of disease. Due to biochemical differences, however, formalin-fixed paraffin-embedded (FFPE) DNA and RNA have generally been extracted separately from either different tissue sections or from the same section by dividing the digested tissue. The former limits accurate correlation whilst the latter is impractical when utilizing rare or limited archived specimens.

Principal Findings

For effective recovery of genomic DNA and total RNA from a single FFPE specimen, without splitting the proteinase-K digested tissue solution, we optimized a co-extraction method by using TRIzol and purifying DNA from the lower aqueous and RNA from the upper organic phases. Using a series of seven different archived specimens, we evaluated the total amounts of genomic DNA and total RNA recovered by our TRIzol-based co-extraction method and compared our results with those from two commercial kits, the Qiagen AllPrep DNA/RNA FFPE kit, for co-extraction, and the Ambion RecoverAll™ Total Nucleic Acid Isolation kit, for separate extraction of FFPE-DNA and -RNA. Then, to accurately assess the quality of DNA and RNA co-extracted from a single FFPE specimen, we used qRT-PCR, gene expression profiling and methylation assays to analyze microRNAs, mRNAs, and genomic DNA recovered from matched fresh and FFPE MCF10A cells. These experiments show that the TRIzol-based co-extraction method provides larger amounts of FFPE-DNA and –RNA than the two other methods, and particularly provides higher quality microRNAs and genomic DNA for subsequent molecular analyses.

Significance

We determined that co-extraction of genomic DNA and total RNA from a single FFPE specimen is an effective recovery approach to obtain high-quality material for parallel molecular and high-throughput analyses. Our optimized approach provides the option of collecting DNA, which would otherwise be discarded or degraded, for additional or subsequent studies.

Identification of high-confidence somatic mutations in whole genome sequence of formalin-fixed breast cancer specimens

The utilization of archived, formalin-fixed paraffin-embedded (FFPE) tumor samples for massive parallel sequencing has been challenging due to DNA damage and contamination with normal stroma. Here, we perform whole genome sequencing of DNA isolated from two triple-negative breast cancer tumors archived for >11 years as 5 µm FFPE sections and matched germline DNA. The tumor samples show differing amounts of FFPE damaged DNA sequencing reads revealed as relatively high alignment mismatch rates enriched for C · G > T · A substitutions compared to germline samples. This increase in mismatch rate is observable with as few as one million reads, allowing for an upfront evaluation of the sample integrity before whole genome sequencing. By applying innovative quality filters incorporating global nucleotide mismatch rates and local mismatch rates, we present a method to identify high-confidence somatic mutations even in the presence of FFPE induced DNA damage. This results in a breast cancer mutational profile consistent with previous studies and revealing potentially important functional mutations. Our study demonstrates the feasibility of performing genome-wide deep sequencing analysis of FFPE archived tumors of limited sample size such as residual cancer after treatment or metastatic biopsies.

Microarray gene expression analysis of fixed archival tissue permits molecular classification and identification of potential therapeutic targets in diffuse large B-cell lymphoma

Refractory/relapsed diffuse large B-cell lymphoma (DLBCL) has a poor prognosis. Novel drugs targeting the constitutively activated NF-κB pathway characteristic of ABC-DLBCL are promising, but evaluation depends on accurate activated B cell-like (ABC)/germinal center B cell-like (GCB) molecular classification. This is traditionally performed on gene microarray expression profiles of fresh biopsies, which are not routinely collected, or by immunohistochemistry on formalin-fixed, paraffin-embedded (FFPE) tissue, which lacks reproducibility and classification accuracy. We explored the possibility of using routine archival FFPE tissue for gene microarray applications. We examined Affymetrix HG U133 Plus 2.0 gene expression profiles from paired archival FFPE and fresh-frozen tissues of 40 ABC/GCB-classified DLBCL cases to compare classification accuracy and test the potential for this approach to aid the discovery of therapeutic targets and disease classifiers in DLBCL. Unsupervised hierarchical clustering of unselected present probe sets distinguished ABC/GCB in FFPE with remarkable accuracy, and a Bayesian classifier correctly assigned 32 of 36 cases with >90% probability. Enrichment for NF-κB genes was appropriately seen in ABC-DLBCL FFPE tissues. The top discriminatory genes expressed in FFPE separated cases with high statistical significance and contained novel biology with potential therapeutic insights, warranting further investigation. These results support a growing understanding that archival FFPE tissues can be used in microarray experiments aimed at molecular classification, prognostic biomarker discovery, and molecular exploration of rare diseases.

Age-dependent brain gene expression and copy number anomalies in autism suggest distinct pathological processes at young versus mature ages

Autism is a highly heritable neurodevelopmental disorder, yet the genetic underpinnings of the disorder are largely unknown. Aberrant brain overgrowth is a well-replicated observation in the autism literature; but association, linkage, and expression studies have not identified genetic factors that explain this trajectory. Few studies have had sufficient statistical power to investigate whole-genome gene expression and genotypic variation in the autistic brain, especially in regions that display the greatest growth abnormality. Previous functional genomic studies have identified possible alterations in transcript levels of genes related to neurodevelopment and immune function. Thus, there is a need for genetic studies involving key brain regions to replicate these findings and solidify the role of particular functional pathways in autism pathogenesis. We therefore sought to identify abnormal brain gene expression patterns via whole-genome analysis of mRNA levels and copy number variations (CNVs) in autistic and control postmortem brain samples. We focused on prefrontal cortex tissue where excess neuron numbers and cortical overgrowth are pronounced in the majority of autism cases. We found evidence for dysregulation in pathways governing cell number, cortical patterning, and differentiation in young autistic prefrontal cortex. In contrast, adult autistic prefrontal cortex showed dysregulation of signaling and repair pathways. Genes regulating cell cycle also exhibited autism-specific CNVs in DNA derived from prefrontal cortex, and these genes were significantly associated with autism in genome-wide association study datasets. Our results suggest that CNVs and age-dependent gene expression changes in autism may reflect distinct pathological processes in the developing versus the mature autistic prefrontal cortex. Our results raise the hypothesis that genetic dysregulation in the developing brain leads to abnormal regional patterning, excess prefrontal neurons, cortical overgrowth, and neural dysfunction in autism.

Autism is a disorder characterized by aberrant social, communication, and restricted and repetitive behaviors. It develops clinically in the first years of life. Toddlers and children with autism often exhibit early brain enlargement and excess neuron numbers in the prefrontal cortex. Adults with autism generally do not display enlargement but instead may have a smaller brain size. Thus, we investigated DNA and mRNA patterns in prefrontal cortex from young versus adult postmortem individuals with autism to identify age-related gene expression differences as well as possible genetic correlates of abnormal brain enlargement, excess neuron numbers, and abnormal functioning in this disorder. We found abnormalities in genetic pathways governing cell number, neurodevelopment, and cortical lateralization in autism. We also found that the key pathways associated with autism are different between younger and older autistic individuals. These findings suggest that dysregulated gene pathways in the early stages of neurodevelopment could lead to later behavioral and cognitive deficits associated with autism.

Clinical relevance of Ki67 gene expression analysis using formalin-fixed paraffin-embedded breast cancer specimens

BACKGROUND

Ki67 is a protein associated with cell cycle activity and shows a good correlation with the growth fraction, which has been proposed as a prognostic or predictive marker in breast cancer. In this study, we aimed to analyze the expression levels of Ki67 (MKI67) messenger RNA (mRNA) derived from formalin-fixed paraffin-embedded (FFPE) tissues for comparison with the immunohistochemical Ki67 labeling index, and investigate the correlation coefficients with clinical outcomes.

METHODS

We analyzed the data of Ki67 mRNA from FFPE and matched fresh-frozen (FF) tissues based on a real-time quantitative reverse-transcription polymerase chain reaction (RT-qPCR) assay system in 203 cases of primary invasive breast cancer.

RESULTS

The correlation between Ki67 mRNA expression of either FFPE or FF specimens and Ki67 labeling index was positive, as was the correlation between the FFPE and FF results (P < 0.0001). Ki67 mRNA expression of FFPE specimens was significantly associated with clinicopathological characteristics: tumor size, lymph node status, nuclear grade, hormone receptors, human epidermal growth factor receptor 2 (Her2) status, and tumor subtype. In prognostic results, Ki67 gene expression in the FFPE specimens revealed almost similar patterns of significance in Kaplan-Meier curves and univariate and multivariate relapse-free survival results as the Ki67 labeling index.

CONCLUSIONS

Gene expression analysis of Ki67 of FFPE specimens could be successfully performed using RT-qPCR, closely resembling the significant clinical characteristics of Ki67 labeling index. These results confirm that Ki67 gene expression of FFPE specimens has potential for evaluation of cell cycle activity of breast cancer specimens.

Preprocessing and Quality Control Strategies for Illumina DASL Assay-Based Brain Gene Expression Studies with Semi-Degraded Samples

Available statistical preprocessing or quality control analysis tools for gene expression microarray datasets are known to greatly affect downstream data analysis, especially when degraded samples, unique tissue samples, or novel expression assays are used. It is therefore important to assess the validity and impact of the assumptions built in to preprocessing schemes for a dataset. We developed and assessed a data preprocessing strategy for use with the Illumina DASL-based gene expression assay with partially degraded postmortem prefrontal cortex samples. The samples were obtained from individuals with autism as part of an investigation of the pathogenic factors contributing to autism. Using statistical analysis methods and metrics such as those associated with multivariate distance matrix regression and mean inter-array correlation, we developed a DASL-based assay gene expression preprocessing pipeline to accommodate and detect problems with microarray-based gene expression values obtained with degraded brain samples. Key steps in the pipeline included outlier exclusion, data transformation and normalization, and batch effect and covariate corrections. Our goal was to produce a clean dataset for subsequent downstream differential expression analysis. We ultimately settled on available transformation and normalization algorithms in the R/Bioconductor package lumi based on an assessment of their use in various combinations. A log2-transformed, quantile-normalized, and batch and seizure-corrected procedure was likely the most appropriate for our data. We empirically tested different components of our proposed preprocessing strategy and believe that our results suggest that a preprocessing strategy that effectively identifies outliers, normalizes the data, and corrects for batch effects can be applied to all studies, even those pursued with degraded samples.

Akt inhibitors MK-2206 and nelfinavir overcome mTOR inhibitor resistance in diffuse large B-cell lymphoma

PURPOSE

The mTOR pathway is constitutively activated in diffuse large B-cell lymphoma (DLBCL). mTOR inhibitors have activity in DLBCL, although response rates remain low. We evaluated DLBCL cell lines with differential resistance to the mTOR inhibitor rapamycin: (i) to identify gene expression profile(s) (GEP) associated with resistance to rapamycin, (ii) to understand mechanisms of rapamycin resistance, and (iii) to identify compounds likely to synergize with mTOR inhibitor.

EXPERIMENTAL DESIGN

We sought to identify a GEP of mTOR inhibitor resistance by stratification of eight DLBCL cell lines with respect to response to rapamycin. Then, using pathway analysis and connectivity mapping, we sought targets likely accounting for this resistance and compounds likely to overcome it. We then evaluated two compounds thus identified for their potential to synergize with rapamycin in DLBCL and confirmed mechanisms of activity with standard immunoassays.

RESULTS

We identified a GEP capable of reliably distinguishing rapamycin-resistant from rapamycin-sensitive DLBCL cell lines. Pathway analysis identified Akt as central to the differentially expressed gene network. Connectivity mapping identified compounds targeting Akt as having a high likelihood of reversing the GEP associated with mTOR inhibitor resistance. Nelfinavir and MK-2206, chosen for their Akt-inhibitory properties, yielded synergistic inhibition of cell viability in combination with rapamycin in DLBCL cell lines, and potently inhibited phosphorylation of Akt and downstream targets of activated mTOR.

CONCLUSIONS

GEP identifies DLBCL subsets resistant to mTOR inhibitor therapy. Combined targeting of mTOR and Akt suppresses activation of key components of the Akt/mTOR pathway and results in synergistic cytotoxicity. These findings are readily adaptable to clinical trials.

Angiogenic mRNA and microRNA gene expression signature predicts a novel subtype of serous ovarian cancer

Ovarian cancer is the fifth leading cause of cancer death for women in the U.S. and the seventh most fatal worldwide. Although ovarian cancer is notable for its initial sensitivity to platinum-based therapies, the vast majority of patients eventually develop recurrent cancer and succumb to increasingly platinum-resistant disease. Modern, targeted cancer drugs intervene in cell signaling, and identifying key disease mechanisms and pathways would greatly advance our treatment abilities. In order to shed light on the molecular diversity of ovarian cancer, we performed comprehensive transcriptional profiling on 129 advanced stage, high grade serous ovarian cancers. We implemented a, re-sampling based version of the ISIS class discovery algorithm (rISIS: robust ISIS) and applied it to the entire set of ovarian cancer transcriptional profiles. rISIS identified a previously undescribed patient stratification, further supported by micro-RNA expression profiles, and gene set enrichment analysis found strong biological support for the stratification by extracellular matrix, cell adhesion, and angiogenesis genes. The corresponding "angiogenesis signature" was validated in ten published independent ovarian cancer gene expression datasets and is significantly associated with overall survival. The subtypes we have defined are of potential translational interest as they may be relevant for identifying patients who may benefit from the addition of anti-angiogenic therapies that are now being tested in clinical trials.

Highly parallel genome-wide expression analysis of single mammalian cells

BACKGROUND

We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs.

METHODOLOGY/PRINCIPAL FINDINGS

The method uses tagged priming and template-switching, resulting in the incorporation of universal PCR priming sites at both ends of the synthesized cDNA for global PCR amplification. Coupled with a whole-genome gene expression microarray platform, we routinely obtain expression correlation values of R(2)~0.76-0.80 between individual cells and R(2)~0.69 between 50 pg total RNA replicates. Expression profiles generated from single cells or 50 pg total RNA correlate well with that generated with higher input (1 ng total RNA) (R(2)~0.80). Also, the assay is sufficiently sensitive to detect, in a single cell, approximately 63% of the number of genes detected with 1 ng input, with approximately 97% of the genes detected in the single-cell input also detected in the higher input.

CONCLUSIONS/SIGNIFICANCE

In summary, our method facilitates whole-genome gene expression profiling in contexts where starting material is extremely limiting, particularly in areas such as the study of progenitor cells in early development and tumor stem cell biology.

Investigation of Archived Formalin-Fixed Paraffin-Embedded Pancreatic Tissue with Whole-Genome Gene Expression Microarray

The use of formalin-fixed, paraffin-embedded (FFPE) tissue overcomes the most prominent issues related to research on relatively rare diseases: limited sample size, availability of control tissue, and time frame. The use of FFPE pancreatic tissue in GEM may be especially challenging due to its very high amounts of ribonucleases compared to other tissues/organs. In choosing pancreatic tissue, we therefore indirectly address the applicability of other FFPE tissues to gene expression microarray (GEM). GEM was performed on archived, routinely fixed, FFPE pancreatic tissue from patients with congenital hyperinsulinism (CHI), insulinoma, and deceased age-appropriate neonates, using whole-genome arrays. Although ribonuclease-rich, we obtained biologically relevant and disease-specific, significant genes; cancer-related genes; genes involved in (a) the regulation of insulin secretion and synthesis, (b) amino acid metabolism, and (c) calcium ion homeostasis. These results should encourage future research and GEM studies on FFPE tissue from the invaluable biobanks available at the departments of pathology worldwide.

Gene-expression signature of vascular invasion in hepatocellular carcinoma

BACKGROUND & AIMS

Vascular invasion is a major predictor of tumor recurrence after surgical treatments for hepatocellular carcinoma (HCC). While macroscopic vascular invasion can be detected by radiological techniques, pre-operative detection of microscopic vascular invasion, which complicates 30-40% of patients with early tumors, remains elusive.

METHODS

A total of 214 patients with hepatocellular carcinoma who underwent resection were included in the study. By using genome-wide gene-expression profiling of 79 hepatitis C-related hepatocellular carcinoma samples (training set), a gene-expression signature associated with vascular invasion was defined. The signature was validated in formalin-fixed paraffin-embedded tissues obtained from an independent set of 135 patients with various etiologies.

RESULTS

A 35-gene signature of vascular invasion was defined in the training set, predicting vascular invasion with an accuracy of 69%. The signature was independently associated with the presence of vascular invasion (OR 3.38, 95% CI 1.48-7.71, p=0.003) along with tumor size (diameter greater than 3 cm, OR 2.66, 95% CI 1.17-6.05, p=0.02). In the validation set, the signature discarded the presence of vascular invasion with a negative predictive value of 0.77, and significantly improved the diagnostic power of tumor size alone (p=0.045).

CONCLUSIONS

The assessment of a gene-expression signature obtained from resected biopsied tumor specimens improved the diagnosis of vascular invasion beyond clinical variable-based prediction. The signature may aid in candidate selection for liver transplantation, and guide the design of clinical trials with experimental adjuvant therapies.

Molecular diagnostics in melanoma: current status and perspectives

CONTEXT

In the current "molecular" era, the advent of technology, such as array-based platforms, systems biology, and genome-wide approaches, has made it possible to examine human cancers, including melanoma, for genetic mutations, deletions, amplification, differentially regulated genes, and epigenetic changes. Advancement in current technologies is such that one can now examine ribonucleic acid (RNA), deoxyribonucleic acid (DNA), and protein directly from the patient's own tumor.

OBJECTIVE

To apply these new technologies in advancing molecular diagnostics in melanoma has historically suffered from a major obstacle, namely, the scarcity of fresh frozen, morphologically defined tumor banks, annotated with clinical information. Recently, some of the new platforms have advanced to permit utilization of formalin-fixed, paraffin-embedded (FFPE) tumor specimens as starting material.

DATA SOURCES

This article reviews the latest technologies applied to FFPE melanoma sections, narrowing its focus on the utility of transcriptional profiling, especially for melastatin; comparative genomic hybridization; BRAF and NRAS mutational analysis; and micro ribonucleic acid profiling.

CONCLUSION

New molecular approaches are emerging and are likely to improve the classification of melanocytic neoplasms.

Genome-wide expression assay comparison across frozen and fixed postmortem brain tissue samples

Background

Gene expression assays have been shown to yield high quality genome-wide data from partially degraded RNA samples. However, these methods have not yet been applied to postmortem human brain tissue, despite their potential to overcome poor RNA quality and other technical limitations inherent in many assays. We compared cDNA-mediated annealing, selection, and ligation (DASL)- and in vitro transcription (IVT)-based genome-wide expression profiling assays on RNA samples from artificially degraded reference pools, frozen brain tissue, and formalin-fixed brain tissue.

Results

The DASL-based platform produced expression results of greater reliability than the IVT-based platform in artificially degraded reference brain RNA and RNA from frozen tissue-based samples. Although data associated with a small sample of formalin-fixed RNA samples were poor when obtained from both assays, the DASL-based platform exhibited greater reliability in a subset of probes and samples.

Conclusions

Our results suggest that the DASL-based gene expression-profiling platform may confer some advantages on mRNA assays of the brain over traditional IVT-based methods. We ultimately consider the implications of these results on investigations of neuropsychiatric disorders.

Gene expression profiling of human whole blood samples with the Illumina WG-DASL assay

Background

Microarray-based gene expression analysis of peripheral whole blood is a common strategy in the development of clinically relevant biomarker panels for a variety of human diseases. However, the results of such an analysis are often plagued by decreased sensitivity and reliability due to the effects of relatively high levels of globin mRNA in whole blood. Globin reduction assays have been shown to overcome such effects, but they require large amounts of total RNA and may induce distinct gene expression profiles. The Illumina whole genome DASL assay can detect gene expression levels using partially degraded RNA samples and has the potential to detect rare transcripts present in highly heterogeneous whole blood samples without the need for globin reduction. We assessed the utility of the whole genome DASL assay in an analysis of peripheral whole blood gene expression profiles.

Results

We find that gene expression detection is significantly increased with the use of whole genome DASL compared to the standard IVT-based direct hybridization. Additionally, globin-probe negative whole genome DASL did not exhibit significant improvements over globin-probe positive whole genome DASL. Globin reduction further increases the detection sensitivity and reliability of both whole genome DASL and IVT-based direct hybridization with little effect on raw intensity correlations. Raw intensity correlations between total RNA and globin reduced RNA were 0.955 for IVT-based direct hybridization and 0.979 for whole genome DASL.

Conclusions

Overall, the detection sensitivity of the whole genome DASL assay is higher than the IVT-based direct hybridization assay, with or without globin reduction, and should be considered in conjunction with globin reduction methods for future blood-based gene expression studies.

Recovery of primary sporocysts in vivo in the Schistosoma mansoni/Biomphalaria glabrata model using a simple fixation method suitable for extraction of genomic DNA and RNA

Detailed studies of host/parasite interactions are currently limited because in situ gene sequencing or monitoring of parasite gene expression is so far limited to genes presenting a high loci copy number in the Schistosome genome or a high level of expression. Indeed, how to investigate the host parasite molecular interplay when parasites are not directly accessible in vivo? Here we describe a method to circumvent this problem and to analyze DNA and RNA of Schistosoma mansoni during the interaction with its intermediate snail host Biomphalaria glabrata. We propose a technique for improved DNA and RNA extraction from the intra-molluscan stage of the parasite recovered after fixation of infected snails in Raillet-Henry solution. The extractions can be used for genetic analysis, transcription studies and microsatellite genotyping.

Illumina whole-genome complementary DNA-mediated annealing, selection, extension and ligation platform: assessing its performance in formalin-fixed, paraffin-embedded samples and identifying invasion pattern-related genes in oral squamous cell carcinoma

High-throughput gene expression profiling from formalin-fixed, paraffin-embedded tissues has become a reality, and several methods are now commercially available. The Illumina whole-genome complementary DNA-mediated annealing, selection, extension and ligation assay (Illumina, Inc) is a full-transcriptome version of the original 512-gene complementary DNA-mediated annealing, selection, extension and ligation assay, allowing high-throughput profiling of 24,526 annotated genes from degraded and formalin-fixed, paraffin-embedded RNA. This assay has the potential to allow identification of novel gene signatures associated with clinical outcome using banked archival pathology specimen resources. We tested the reproducibility of the whole-genome complementary DNA-mediated annealing, selection, extension and ligation assay and its sensitivity for detecting differentially expressed genes in RNA extracted from matched fresh and formalin-fixed, paraffin-embedded cells, after 1 and 13 months of storage, using the human breast cell lines MCF7 and MCF10A. Then, using tumor worst pattern of invasion as a classifier, 1 component of the "risk model," we selected 12 formalin-fixed, paraffin-embedded oral squamous cell carcinomas for whole-genome complementary DNA-mediated annealing, selection, extension and ligation assay analysis. We profiled 5 tumors with nonaggressive, nondispersed pattern of invasion, and 7 tumors with aggressive dispersed pattern of invasion and satellites scattered at least 1 mm apart. To minimize variability, the formalin-fixed, paraffin-embedded specimens were prepared from snap-frozen tissues, and RNA was obtained within 24 hours of fixation. One hundred four down-regulated genes and 72 up-regulated genes in tumors with aggressive dispersed pattern of invasion were identified. We performed quantitative reverse transcriptase polymerase chain reaction validation of 4 genes using Taqman assays and in situ protein detection of 1 gene by immunohistochemistry. Functional cluster analysis of genes up-regulated in tumors with aggressive pattern of invasion suggests presence of genes involved in cellular cytoarchitecture, some of which already associated with tumor invasion. Identification of these genes provides biologic rationale for our histologic classification, with regard to tumor invasion, and demonstrates that the whole-genome complementary DNA-mediated annealing, selection, extension and ligation assay is a powerful assay for profiling degraded RNA from archived specimens when combined with quantitative reverse transcriptase polymerase chain reaction validation.

Formalin fixation at low temperature better preserves nucleic acid integrity

Fixation with formalin, a widely adopted procedure to preserve tissue samples, leads to extensive degradation of nucleic acids and thereby compromises procedures like microarray-based gene expression profiling. We hypothesized that RNA fragmentation is caused by activation of RNAses during the interval between formalin penetration and tissue fixation. To prevent RNAse activation, a series of tissue samples were kept under-vacuum at 4°C until fixation and then fixed at 4°C, for 24 hours, in formalin followed by 4 hours in ethanol 95%. This cold-fixation (CF) procedure preserved DNA and RNA, so that RNA segments up to 660 bp were efficiently amplified. Histological and immunohistochemical features were fully comparable with those of standard fixation. Microarray-based gene expression profiles were comparable with those obtained on matched frozen samples for probes hybridizing within 700 bases from the reverse transcription start site. In conclusion, CF preserves tissues and nucleic acids, enabling reliable gene expression profiling of fixed tissues.

Potential prognostic biomarkers for bone metastasis from hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) most commonly develops in patients who have a viral infection, especially in the case of hepatitis B virus (HBV), and in patients with a chronic liver disease. HCC patients with bone metastasis (BM) suffer from pain and other symptoms that significantly reduce their quality of life. Identification of patients who are at high risk for BM after undergoing potentially curative treatment for HCC remains challenging. Here, we aimed to identify HCC BM-related genes and proteins to establish prediction biomarkers.

METHODS

RNA was extracted from 48 pairs of intratumoral and peritumoral formalin-fixed, paraffin-embedded tissue from HCC patients with and without BM. A cDNA-mediated annealing, selection, extension and ligation assay containing 502 cancer-related genes was used to identify novel BM-associated genes. An additional independent study with 350 HCC patients who had undergone hepatectomy was conducted to evaluate the expression of candidate genes at the protein level using immunohistochemistry on tissue microarrays (TMAs). Of the 350 patients, 273 (78.0%) were infected with HBV.

RESULTS

Seven intratumoral genes and 17 peritumoral genes were overexpressed in patients with BM, whereas 15 intratumoral genes and 28 peritumoral genes were underexpressed in patients with BM. We selected the following four genes for further analysis because they were differentially expressed in the cancer gene-specific microarray and were previously reported to be associated with BM: connective tissue growth factor (CTGF), matrix metalloproteinase-1 (MMP-1), transforming growth factor β1 (TGF-β1), and interleukin-11 (IL-11). We assessed the protein expression of these selected genes using immunohistochemistry on TMAs including 350 HCC patient specimens. We determined that expression of intratumoral CTGF, intratumoral IL-11, and peritumoral MMP-1 were independent prognostic factors for developing BM in HCC patients. Combining intratumoral CTGF and IL-11 expression was also an independent risk factor for BM development.

CONCLUSIONS

Sixty-seven genes were differentially expressed in HCC patients with and without BM. High intratumoral CTGF, positive IL-11, and high peritumoral MMP-1 expression were associated with BM after hepatectomy. Intratumoral CTGF expression combined with IL-11 expression may serve as a useful predictive biomarker for HCC BM.

Phase II study of aflibercept in recurrent malignant glioma: a North American Brain Tumor Consortium study

PURPOSE

Antivascular endothelial growth factor (anti-VEGF) therapy is a promising treatment approach for patients with recurrent glioblastoma. This single-arm phase II study evaluated the efficacy of aflibercept (VEGF Trap), a recombinantly produced fusion protein that scavenges both VEGF and placental growth factor in patients with recurrent malignant glioma.

PATIENTS AND METHODS

Forty-two patients with glioblastoma and 16 patients with anaplastic glioma who had received concurrent radiation and temozolomide and adjuvant temozolomide were enrolled at first relapse. Aflibercept 4 mg/kg was administered intravenously on day 1 of every 2-week cycle.

RESULTS

The 6-month progression-free survival rate was 7.7% for the glioblastoma cohort and 25% for patients with anaplastic glioma. Overall radiographic response rate was 24% (18% for glioblastoma and 44% for anaplastic glioma). The median progression-free survival was 24 weeks for patients with anaplastic glioma (95% CI, 5 to 31 weeks) and 12 weeks for patients with glioblastoma (95% CI, 8 to 16 weeks). A total of 14 patients (25%) were removed from the study for toxicity, on average less than 2 months from treatment initiation. The main treatment-related National Cancer Institute Common Terminology Criteria grades 3 and 4 adverse events (38 total) included fatigue, hypertension, and lymphopenia. Two grade 4 CNS ischemias and one grade 4 systemic hemorrhage were reported. Aflibercept rapidly decreases permeability on dynamic contrast enhanced magnetic resonance imaging, and molecular analysis of baseline tumor tissue identified tumor-associated markers of response and resistance.

CONCLUSION

Aflibercept monotherapy has moderate toxicity and minimal evidence of single-agent activity in unselected patients with recurrent malignant glioma.

Dysregulated TRK signalling is a therapeutic target in CYLD defective tumours

Individuals with germline mutations in the tumour suppressor gene CYLD are at high risk of developing disfiguring cutaneous appendageal tumours, the defining tumour being the highly organised cylindroma. Here, we analysed CYLD mutant tumour genomes by array comparative genomic hybridisation (aCGH) and gene expression microarray analysis. CYLD mutant tumours were characterised by an absence of copy number aberrations apart from loss-of-heterozygosity at chromosome 16q, the genomic location of the CYLD gene. Gene expression profiling of CYLD mutant tumours revealed dysregulated tropomyosin kinase (TRK) signalling with overexpression of TRKB and TRKC in tumours when compared to perilesional skin. Immunohistochemical analysis of a tumour microarray demonstrated strong membranous TRKB and TRKC staining in cylindromas, as well as elevated levels of ERK phosphorylation and BCL2 expression. Membranous TRKC overexpression was also observed in 70% of sporadic basal cell carcinomas. RNA interference mediated silencing of TRKB and TRKC, as well as treatment with the small molecule TRK inhibitor lestaurtinib, reduced colony formation and proliferation in three-dimensional primary cell cultures established from CYLD mutant tumours. These results suggest that TRK inhibition could be used as a strategy to treat tumours with loss of functional CYLD.

MicroRNA-based classification of hepatocellular carcinoma and oncogenic role of miR-517a

BACKGROUND & AIMS

Hepatocellular carcinoma (HCC) is a heterogeneous tumor that develops via activation of multiple pathways and molecular alterations. It has been a challenge to identify molecular classes of HCC and design treatment strategies for each specific subtype. MicroRNAs (miRNAs) are involved in HCC pathogenesis, and their expression profiles have been used to classify cancers. We analyzed miRNA expression in human HCC samples to identify molecular subclasses and oncogenic miRNAs.

METHODS

We performed miRNA profiling of 89 HCC samples using a ligation-mediated amplification method. Subclasses were identified by unsupervised clustering analysis. We identified molecular features specific for each subclass using expression pattern (Affymetrix U133 2.0; Affymetrix, Santa Clara, CA), DNA change (Affymetrix STY Mapping Array), mutation (CTNNB1), and immunohistochemical (phosphor[p]-protein kinase B, p-insulin growth factor-IR, p-S6, p-epidermal growth factor receptor, β-catenin) analyses. The roles of selected miRNAs were investigated in cell lines and in an orthotopic model of HCC.

RESULTS

We identified 3 main clusters of HCCs: the wingless-type MMTV integration site (32 of 89; 36%), interferon-related (29 of 89; 33%), and proliferation (28 of 89; 31%) subclasses. A subset of patients with tumors in the proliferation subclass (8 of 89; 9%) overexpressed a family of poorly characterized miRNAs from chr19q13.42. Expression of miR-517a and miR-520c (from ch19q13.42) increased proliferation, migration, and invasion of HCC cells in vitro. MiR-517a promoted tumorigenesis and metastatic dissemination in vivo.

CONCLUSIONS

We propose miRNA-based classification of 3 subclasses of HCC. Among the proliferation class, miR-517a is an oncogenic miRNA that promotes tumor progression. There is rationale for developing therapies that target miR-517a for patients with HCC.

Rutgers EJ, Wessels LFA, Van't Veer LJ. Gene expression profiles from formalin fixed paraffin embedded breast cancer tissue are largely comparable to fresh frozen matched tissue

Background and Methods

Formalin Fixed Paraffin Embedded (FFPE) samples represent a valuable resource for cancer research. However, the discovery and development of new cancer biomarkers often requires fresh frozen (FF) samples. Recently, the Whole Genome (WG) DASL (cDNA-mediated Annealing, Selection, extension and Ligation) assay was specifically developed to profile FFPE tissue. However, a thorough comparison of data generated from FFPE RNA and Fresh Frozen (FF) RNA using this platform is lacking. To this end we profiled, in duplicate, 20 FFPE tissues and 20 matched FF tissues and evaluated the concordance of the DASL results from FFPE and matched FF material.

Methodology and Principal Findings

We show that after proper normalization, all FFPE and FF pairs exhibit a high level of similarity (Pearson correlation >0.7), significantly larger than the similarity between non-paired samples. Interestingly, the probes showing the highest correlation had a higher percentage G/C content and were enriched for cell cycle genes. Predictions of gene expression signatures developed on frozen material (Intrinsic subtype, Genomic Grade Index, 70 gene signature) showed a high level of concordance between FFPE and FF matched pairs. Interestingly, predictions based on a 60 gene DASL list (best match with the 70 gene signature) showed very high concordance with the MammaPrint® results.

Conclusions and Significance

We demonstrate that data generated from FFPE material with the DASL assay, if properly processed, are comparable to data extracted from the FF counterpart. Specifically, gene expression profiles for a known set of prognostic genes for a specific disease are highly comparable between two conditions. This opens up the possibility of using both FFPE and FF material in gene expressions analyses, leading to a vast increase in the potential resources available for cancer research.

Genome-wide gene expression profiling of formalin-fixed paraffin-embedded breast cancer core biopsies using microarrays

The routine workflow for invasive cancer diagnostics includes biopsy processing by formalin fixation and paraffin embedding. It has been shown only recently that this kind of sample can be used for gene expression analysis with microarrays. To support this view, the authors conducted a microarray study using formalin-fixed paraffin-embedded (FFPE) core needle biopsies from breast cancers. Typically, for the 3'-biased chip type that was used, the probe sets interrogate sequences near the poly-A-tail of the transcripts, and this kind of probe turned out to be suitable to measure RNA levels in FFPE biopsies. For ER and HER2, the authors observed strong correlations between RNA levels and protein expression (p = 0.000003 and p = 0.0022). ER and HER2 classification of the biopsies by the RNA levels was feasible with high sensitivity and specificity (AUROC = 0.93 and AUROC = 0.96). Furthermore, a signature of 346 genes was identified that correlated with ER and a signature of 528 genes that correlated with HER2 protein status. Many of these genes (ER: 63%) could be confirmed by analysis of gene expression data from frozen tissues. The findings support the notion that clinically relevant information can be gained from microarray analyses of FFPE cancer biopsies. This opens new opportunities for biomarker detection studies and the integration of microarrays into the workflow of cancer diagnostics.

Gene expression profiling in formalin-fixed, paraffin-embedded tissues using the whole-genome DASL assay

Here, we provide a detailed technical description of a gene expression assay (Whole-Genome DASL (WG-DASL)), which not only enables whole-genome transcriptional profiling of degraded material, such as formalin-fixed, paraffin-embedded tissues, but is also capable of generating robust profiles with low input intact RNA. The WG-DASL assay combines target-specific annealing, extension, and ligation events followed by universal PCR and labeling steps to generate highly multiplexed Cy3-labeled products. These short products, which are single-stranded, are directly hybridized to a whole-genome expression BeadChip (HumanRef-8) containing probe content corresponding to ∼24 K RefSeq transcripts. After washing and imaging, fluorescence emissions are quantitatively recorded for each probe using high-resolution confocal scanners and imaging software. GenomeStudio software allows direct analysis of mRNA expression data and provides results in standard file formats that can be readily exported and analyzed with most standard gene expression analysis software programs. This technology is particularly useful for genome-wide expression profiling in degraded, archived material, including limited quantities of clinical samples, such as microdissected and biopsied materials.

Identification of candidate genes for histiocytoid cardiomyopathy (HC) using whole genome expression analysis: analyzing material from the HC registry

Histiocytoid cardiomyopathy (HC) is a rare but distinctive arrhythmogenic disorder characterized by incessant ventricular tachycardia, cardiomegaly, and often sudden death by age 2 years. The underlying genetic mechanism of HC has eluded researchers for decades. To further identify the potential molecular-genetic bases of HC, molecular analyses of HC hearts and hearts of age-matched controls were performed. Total RNA and genomic DNA were prepared from formalin-fixed, paraffin-embedded cardiac tissue from 12 cases of HC and 12 age-matched controls. To identify genes differentially expressed in HC, whole genome cDNA-mediated annealing, selection, extension, and ligation profiling was performed. TaqMan quantitative polymerase chain reaction confirmed changes in RNA expression. DNA copy number changes were measured by TaqMan copy number variant analysis. Analysis of differential gene expression in HC cases identified 2 significantly downregulated gene sets aligned sequentially along the genome. The 1st gene cluster consisted of genes S100A8 , S100A9 , and S100A12 at 1q21.3c, and the 2nd cluster consisted of genes IL1RL1 ( ST2 ), IL18R1 , and IL18RAP at 2q12.1a. Strong decreases in interleukin 33 expression were also observed. Decreases in copy number of the S100A genes were confirmed by TaqMan copy number variant assays. S100A genes are downstream of the p38-MAPK pathway that can be activated by interleukin 33 signaling. These data suggest a model in which the interleukin 33-IL1RL1/p38-MAPK/ S100A8-S100A9 axis is downregulated in HC cardiac tissue and provide several candidate genes on 1q21.3c and 2q12.1a for inherited mutations that may predispose individuals to HC.

Expression profiling of formalin-fixed paraffin-embedded primary breast tumors using cancer-specific and whole genome gene panels on the DASL® platform

Background

The cDNA-mediated Annealing, extension, Selection and Ligation (DASL) assay has become a suitable gene expression profiling system for degraded RNA from paraffin-embedded tissue. We examined assay characteristics and the performance of the DASL 502-gene Cancer Panel^v1 (1.5K) and 24,526-gene panel (24K) platforms at differentiating nine human epidermal growth factor receptor 2- positive (HER2+) and 11 HER2-negative (HER2-) paraffin-embedded breast tumors.

Methods

Bland-Altman plots and Spearman correlations evaluated intra/inter-panel agreement of normalized expression values. Unequal-variance t-statistics tested for differences in expression levels between HER2 + and HER2 - tumors. Regulatory network analysis was performed using Metacore (GeneGo Inc., St. Joseph, MI).

Results

Technical replicate correlations ranged between 0.815-0.956 and 0.986-0.997 for the 1.5K and 24K panels, respectively. Inter-panel correlations of expression values for the common 498 genes across the two panels ranged between 0.485-0.573. Inter-panel correlations of expression values of 17 probes with base-pair sequence matches between the 1.5K and 24K panels ranged between 0.652-0.899. In both panels, erythroblastic leukemia viral oncogene homolog 2 (ERBB2) was the most differentially expressed gene between the HER2 + and HER2 - tumors and seven additional genes had p-values < 0.05 and log2 -fold changes > |0.5| in expression between HER2 + and HER2 - tumors: topoisomerase II alpha (TOP2A), cyclin a2 (CCNA2), v-fos fbj murine osteosarcoma viral oncogene homolog (FOS), wingless-type mmtv integration site family, member 5a (WNT5A), growth factor receptor-bound protein 7 (GRB7), cell division cycle 2 (CDC2), and baculoviral iap repeat-containing protein 5 (BIRC5). The top 52 discriminating probes from the 24K panel are enriched with genes belonging to the regulatory networks centered around v-myc avian myelocytomatosis viral oncogene homolog (MYC), tumor protein p53 (TP53), and estrogen receptor α (ESR1). Network analysis with a two-step extension also showed that the eight discriminating genes common to the 1.5K and 24K panels are functionally linked together through MYC, TP53, and ESR1.

Conclusions

The relative RNA abundance obtained from two highly differing density gene panels are correlated with eight common genes differentiating HER2 + and HER2 - breast tumors. Network analyses demonstrated biological consistency between the 1.5K and 24K gene panels.

Analyses and interpretation of whole-genome gene expression from formalin-fixed paraffin-embedded tissue: an illustration with breast cancer tissues

Background

We evaluated (a) the feasibility of whole genome cDNA-mediated Annealing, Selection, extension and Ligation (DASL) assay on formalin-fixed paraffin-embedded (FFPE) tissue and (b) whether similar conclusions can be drawn by examining FFPE samples as proxies for fresh frozen (FF) tissues. We used a whole genome DASL assay (addressing 18,391 genes) on a total of 72 samples from paired breast tumor and surrounding healthy tissues from both FF and FFPE samples.

Results

Gene detection was very good with comparable success between the FFPE and FF samples. Reproducibility was also high (r² = 0.98); however, concordance between the two types of samples was low. Only one-third of the differentially expressed genes in tumor tissues (compared to corresponding normal) from FF samples could be detected in FFPE samples and conversely only one-fourth of the differentially expressed genes from FFPE samples could be detected in FF samples. GO-enrichment analysis, gene set enrichment analysis (GSEA) and GO-ANOVA analyses also suggested small overlap between the lead functional groups that were differentially expressed in tumor detectable by examining FFPE and FF samples. In other words, FFPE samples may not be ideal for picking individual target gene(s), but may be used to identify some of the lead functional group(s) of genes that are differentially expressed in tumor. The differentially expressed genes in breast cancer found in our study were biologically meaningful. The "cell cycle" & "cell division" related genes were up-regulated and genes related to "regulation of epithelial cell proliferation" were down-regulated.

Conclusions

Gene expression experiments using the DASL assay can efficiently handle fragmentation issues in the FFPE tissues. However, formalin fixation seems to change RNA and consequently significantly alters gene expression in a number of genes which may not be uniform between tumor and normal tissues. Therefore, considerable caution needs to be taken when interpreting gene expression data from FFPE tissues, especially in relation to specific genes.

Gene expression array testing of FFPE archival breast tumor samples: an optimized protocol for WG-DASL sample preparation

Archived formalin-fixed, paraffin embedded (FFPE) tissues constitute a vast, well-annotated, but underexploited resource for the molecular study of cancer progression, largely because degradation, chemical modification, and cross-linking, render FFPE RNA a suboptimal substrate for conventional analytical methods. We report here a modified protocol for RNA extraction from FFPE tissues which maximized the success rate (with 100% of samples) in the expression profiling of a set of 60 breast cancer samples on the WG-DASL platform; yielding data of sufficient quality such that in hierarchical clustering (a) 12/12 (100%) replicates correctly identified their respective counterparts, with a high self-correlation (r = 0.979), and (b) the overall sample set grouped with high specificity into ER+ (38/40; 95%) and ER- (18/20; 90%) subtypes. These results indicate that a large fraction of decade-old FFPE samples, of diverse institutional origins and processing histories, can yield RNA suitable for gene expression profiling experiments.