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

MiRNA expression deregulation correlates with the Oncotype DX® DCIS score

Background

Current clinical criteria do not discriminate well between women who will or those who will not develop ipsilateral invasive breast cancer (IBC), or a DCIS recurrence after a ductal carcinoma in situ (DCIS) diagnosis. The 12-gene Oncotype DX® DCIS assay (RT qPCR gene-based scoring system) was established and shown to predict the risk of subsequent ipsilateral IBC or DCIS recurrence. Recent studies have shown that microRNA (miRNA) expression deregulation can contribute to the development of IBC, but very few have evaluated miRNA deregulation in DCIS lesions. In this study, we sought to determine whether specific miRNA expression changes may correlate with Oncotype DX® DCIS scores.

Methods

For this study, we used archived formalin-fixed, paraffin-embedded (FFPE) specimens from 41 women diagnosed with DCIS between 2012 and 2018. The DCIS lesions were stratified into low (n = 26), intermediate (n = 10), and high (n = 5) risk score groups using the Oncotype DX® DCIS assay. Total RNA was extracted from DCIS lesions by macro-dissection of unstained FFPE sections, and next-generation small-RNA sequencing was performed. We evaluated the correlation between miRNA expression data and Oncotype score, as well as patient age. RT-qPCR validations were performed to validate the topmost differentially expressed miRNAs identified between the different risk score groups.

Results

MiRNA sequencing of 32 FFPE DCIS specimens from the three different risk group scores identified a correlation between expression deregulation of 17 miRNAs and Oncotype scores. Our analyses also revealed a correlation between the expression deregulation of 9 miRNAs and the patient’s age. Based on these results, a total of 15 miRNAs were selected for RT-qPCR validation. Of these, miR-190b (p = 0.043), miR-135a (p = 0.05), miR-205 (p = 0.00056), miR-30c (p = 0.011), and miR-744 (p = 0.038) showed a decreased expression in the intermediate/high Oncotype group when compared to the low-risk score group. A composite risk score was established using these 5 miRNAs and indicated a significant association between miRNA expression deregulation and the Oncotype DX® DCIS Score (p < 0.0021), between high/intermediate and low risk groups.

Conclusions

Our analyses identified a subset of 5 miRNAs able to discriminate between Oncotype DX® DCIS score subgroups. Together, our data suggest that miRNA expression analysis may add value to the predictive and prognostic evaluation of DCIS lesions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-022-01558-4.

OXTRHigh stroma fibroblasts control the invasion pattern of oral squamous cell carcinoma via ERK5 signaling

The Pattern Of Invasion (POI) of tumor cells into adjacent normal tissues clinically predicts postoperative tumor metastasis/recurrence of early oral squamous cell carcinoma (OSCC), but the mechanisms underlying the development of these subtypes remain unclear. Focusing on the highest score of POIs (Worst POI, WPOI) present within each tumor, we observe a disease progression-driven shift of WPOI towards the high-risk type 4/5, associated with a mesenchymal phenotype in advanced OSCC. WPOI 4-5-derived cancer-associated fibroblasts (CAFs^WPOI4-5), characterized by high oxytocin receptor expression (OXTR^High), contribute to local-regional metastasis. OXTR^High CAFs induce a desmoplastic stroma and CCL26 is required for the invasive phenotype of CCR3⁺ tumors. Mechanistically, OXTR activates nuclear ERK5 transcription signaling via Gαq and CDC37 to maintain high levels of OXTR and CCL26. ERK5 ablation reprograms the pro-invasive phenotype of OXTR^High CAFs. Therefore, targeting ERK5 signaling in OXTR^High CAFs is a potential therapeutic strategy for OSCC patients with WPOI 4-5.

Worst pattern of invasion (WPOI) is a parameter used to quantify tumor invasiveness of oral squamous cell carcinoma (OSCC). Here the authors show that a fibroblast subset characterized by the expression of the oxytocin receptor is enriched in highly invasive WPOI 4-5 OSCC tumors and can be targeted to reduce the desmoplastic stroma and tumor metastasis.

Why Formalin-fixed, Paraffin-embedded Biospecimens Must Be Used in Genomic Medicine: An Evidence-based Review and Conclusion

Fresh-frozen tissue is the “gold standard” biospecimen type for next-generation sequencing (NGS). However, collecting frozen tissue is usually not feasible because clinical workflows deliver formalin-fixed, paraffin-embedded (FFPE) tissue blocks. Some clinicians and researchers are reticent to embrace the use of FFPE tissue for NGS because FFPE tissue can yield low quantities of degraded DNA, containing formalin-induced mutations. We describe the process by which formalin-induced deamination can lead to artifactual cytosine (C) to thymine (T) and guanine (G) to adenine (A) (C:G > T:A) mutation calls and perform a literature review of 17 publications that compare NGS data from patient-matched fresh-frozen and FFPE tissue blocks. We conclude that although it is indeed true that sequencing data from FFPE tissue can be poorer than those from frozen tissue, any differences occur at an inconsequential magnitude, and FFPE biospecimens can be used in genomic medicine with confidence:

Effect of tumor microenvironment on pathogenesis of the head and neck squamous cell carcinoma: a systematic review

The tumor microenvironment (TME) is comprised of many different cell populations, such as cancer-associated fibroblasts and various infiltrating immune cells, and non-cell components of extracellular matrix. These crucial parts of the surrounding stroma can function as both positive and negative regulators of all hallmarks of cancer development, including evasion of apoptosis, induction of angiogenesis, deregulation of the energy metabolism, resistance to the immune detection and destruction, and activation of invasion and metastasis. This review represents a summary of recent studies focusing on describing these effects of microenvironment on initiation and progression of the head and neck squamous cell carcinoma, focusing on oral squamous cell carcinoma, since it is becoming clear that an investigation of differences in stromal composition of the head and neck squamous cell carcinoma microenvironment and their impact on cancer development and progression may help better understand the mechanisms behind different responses to therapy and help define possible targets for clinical intervention.

Vimentin and Non-Muscle Myosin IIA are Members of the Neural Precursor Cell Expressed Developmentally Down-Regulated 9 (NEDD9) Interactome in Head and Neck Squamous Cell Carcinoma Cells

Here we demonstrate an interaction between neural precursor cell expressed, developmentally-downregulated 9 (NEDD9) and the cytoskeletal proteins vimentin and non-muscle myosin IIA (NMIIA), based on co-immunoprecipitation and mass spectrometric sequence identification. Vimentin was constitutively phosphorylated at Ser56 but vimentin associated with NEDD9-was not phosphorylated at Ser56. In contrast, NMIIA bound to NEDD9 was phosphorylated on S1943 consistent with its function in invasion and secretion. Treatment of cells with the vimentin-targeting steroidal lactone withaferin A had no effect on vimentin turnover as previously reported, instead causing NEDD9 cleavage and cell death. The NMIIA-selective inhibitor blebbistatin induced cells to form long extensions and attenuated secretion of matrix metalloproteinases (MMPs) 2 and 9. While the site of vimentin interaction on NEDD9 was not defined, NMIIA was found to interact with NEDD9 at its substrate domain. NEDD9 interactions with vimentin and NMIIA are consistent with these proteins having roles in MMP secretion and cell invasion. These findings suggest that a better understanding of NEDD9 signaling is likely to reveal novel therapeutic targets for the prevention of invasion and metastasis.

NEDD9 stimulated MMP9 secretion is required for invadopodia formation in oral squamous cell carcinoma

Neural precursor cell expressed developmentally downregulated 9 (NEDD9) is a component of the metastatic signatures of melanoma, breast cancer, glioblastoma, lung cancer and head and neck squamous cell carcinoma (HNSCC). Here we tested the efficacy of NEDD9's domains in stimulating matrix metalloproteinase (MMP) secretion and invadopodia formation in cells stably expressing various NEDD9 mutants. Replacement of the 13 YxxP motif substrate domain (SD) tyrosines and the C-terminal Y629 with phenylalanines (F14NEDD9) eliminated tyrosine phosphorylation, MMP9 secretion and loss of invadopodia formation. Mutation of the N-terminal SH3 domain Y12 to glutamic acid (Y12ENEDD9) or phenylalanine (Y12FNEDD9) reduced MMP9 secretion and inhibited invadopodia formation. SH3 domain deletion (∆SH3NEDD9) resulted in the loss of MMP9 secretion and a lack of invadopodia formation. The SH3–SD domain (SSNEDD9) construct exhibited tyrosine phosphorylation and stimulated MMP9 secretion, as did ∆CTNEDD9 which lacked the C-terminus (∆C-terminal; ∆CT). E13NEDD9 expression blocked MMP9 secretion and invadopodia formation. MICAL1 (molecule interacting with Cas-L1) silencing with a short hairpin RNA reduced MMP9 secretion, vimentin and E-cadherin levels while increasing N-cadherin and Rab6 levels, consistent with reduced invasive behavior. These findings indicate that NEDD9 SD phosphorylation and SH3 domain interactions are necessary for increasing MMP9 secretion and invadopodia formation.

Retrospective MicroRNA Sequencing: Complementary DNA Library Preparation Protocol Using Formalin-fixed Paraffin-embedded RNA Specimens

-Archived, clinically classified formalin-fixed paraffin-embedded (FFPE) tissues can provide nucleic acids for retrospective molecular studies of cancer development. By using non-invasive or pre-malignant lesions from patients who later develop invasive disease, gene expression analyses may help identify early molecular alterations that predispose to cancer risk. It has been well described that nucleic acids recovered from FFPE tissues have undergone severe physical damage and chemical modifications, which make their analysis difficult and generally requires adapted assays. MicroRNAs (miRNAs), however, which represent a small class of RNA molecules spanning only up to ~18-24 nucleotides, have been shown to withstand long-term storage and have been successfully analyzed in FFPE samples. Here we present a 3' barcoded complementary DNA (cDNA) library preparation protocol specifically optimized for the analysis of small RNAs extracted from archived tissues, which was recently demonstrated to be robust and highly reproducible when using archived clinical specimens stored for up to 35 years. This library preparation is well adapted to the multiplex analysis of compromised/degraded material where RNA samples (up to 18) are ligated with individual 3' barcoded adapters and then pooled together for subsequent enzymatic and biochemical preparations prior to analysis. All purifications are performed by polyacrylamide gel electrophoresis (PAGE), which allows size-specific selections and enrichments of barcoded small RNA species. This cDNA library preparation is well adapted to minute RNA inputs, as a pilot polymerase chain reaction (PCR) allows determination of a specific amplification cycle to produce optimal amounts of material for next-generation sequencing (NGS). This approach was optimized for the use of degraded FFPE RNA from specimens archived for up to 35 years and provides highly reproducible NGS data.

Apolipoprotein E Promotes Invasion in Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) patients generally have a poor prognosis, because of the invasive nature of these tumors. In comparing transcription profiles between OSCC tumors with a more invasive (worst pattern of tumor invasion 5) versus a less invasive (worst pattern of tumor invasion 3) pattern of invasion, we identified a total of 97 genes that were overexpressed at least 1.5-fold in the more invasive tumor subtype. The most functionally relevant genes were assessed using in vitro invasion assays with an OSCC cell line (UM-SCC-1). Individual siRNA knockdown of 15 of these 45 genes resulted in significant reductions in tumor cell invasion compared to a nontargeting siRNA control. One gene whose knockdown had a strong effect on invasion corresponded to apolipoprotein E (APOE). Both matrix degradation and the number of mature invadopodia were significantly decreased with APOE knockdown. APOE knockdown also resulted in increased cellular cholesterol, consistent with APOE's role in regulating cholesterol efflux. APOE knockdown resulted in decreased levels of phospho-extracellular signal-regulated kinase 1/2, phospho-c-Jun N-terminal kinase, and phospho-cJun, as well as decreased activator protein 1 (AP-1) activity. Expression of matrix metalloproteinase 7 (MMP7), an AP-1 target, was also significantly decreased. Our findings suggest that APOE protein plays a significant role in OSCC tumor invasion because of its effects on cellular cholesterol and subsequent effects on cell signaling and AP-1 activity, leading to changes in the expression of invasion-related proteins, including MMP7.

Guidelines for Validation of Next-Generation Sequencing-Based Oncology Panels: A Joint Consensus Recommendation of the Association for Molecular Pathology and College of American Pathologists

Next-generation sequencing (NGS) methods for cancer testing have been rapidly adopted by clinical laboratories. To establish analytical validation best practice guidelines for NGS gene panel testing of somatic variants, a working group was convened by the Association of Molecular Pathology with liaison representation from the College of American Pathologists. These joint consensus recommendations address NGS test development, optimization, and validation, including recommendations on panel content selection and rationale for optimization and familiarization phase conducted before test validation; utilization of reference cell lines and reference materials for evaluation of assay performance; determining of positive percentage agreement and positive predictive value for each variant type; and requirements for minimal depth of coverage and minimum number of samples that should be used to establish test performance characteristics. The recommendations emphasize the role of laboratory director in using an error-based approach that identifies potential sources of errors that may occur throughout the analytical process and addressing these potential errors through test design, method validation, or quality controls so that no harm comes to the patient. The recommendations contained herein are intended to assist clinical laboratories with the validation and ongoing monitoring of NGS testing for detection of somatic variants and to ensure high quality of sequencing results.

Fixing Formalin: A Method to Recover Genomic-Scale DNA Sequence Data from Formalin-Fixed Museum Specimens Using High-Throughput Sequencing

For 150 years or more, specimens were routinely collected and deposited in natural history collections without preserving fresh tissue samples for genetic analysis. In the case of most herpetological specimens (i.e. amphibians and reptiles), attempts to extract and sequence DNA from formalin-fixed, ethanol-preserved specimens—particularly for use in phylogenetic analyses—has been laborious and largely ineffective due to the highly fragmented nature of the DNA. As a result, tens of thousands of specimens in herpetological collections have not been available for sequence-based phylogenetic studies. Massively parallel High-Throughput Sequencing methods and the associated bioinformatics, however, are particularly suited to recovering meaningful genetic markers from severely degraded/fragmented DNA sequences such as DNA damaged by formalin-fixation. In this study, we compared previously published DNA extraction methods on three tissue types subsampled from formalin-fixed specimens of Anolis carolinensis, followed by sequencing. Sufficient quality DNA was recovered from liver tissue, making this technique minimally destructive to museum specimens. Sequencing was only successful for the more recently collected specimen (collected ~30 ybp). We suspect this could be due either to the conditions of preservation and/or the amount of tissue used for extraction purposes. For the successfully sequenced sample, we found a high rate of base misincorporation. After rigorous trimming, we successfully mapped 27.93% of the cleaned reads to the reference genome, were able to reconstruct the complete mitochondrial genome, and recovered an accurate phylogenetic placement for our specimen. We conclude that the amount of DNA available, which can vary depending on specimen age and preservation conditions, will determine if sequencing will be successful. The technique described here will greatly improve the value of museum collections by making many formalin-fixed specimens available for genetic analysis.

Current challenges in bacterial transcriptomics

Over the past decade or so, dramatic developments in our ability to experimentally determine the content and function of genomes have taken place. In particular, next-generation sequencing technologies are now inspiring a new understanding of bacterial transcriptomes on a global scale. In bacterial cells, whole-transcriptome studies have not received attention, owing to the general view that bacterial genomes are simple. However, several recent RNA sequencing results are revealing unexpected levels of complexity in bacterial transcriptomes, indicating that the transcribed regions of genomes are much larger and complex than previously anticipated. In particular, these data show a wide array of small RNAs, antisense RNAs, and alternative transcripts. Here, we review how current transcriptomics are now revolutionizing our understanding of the complexity and regulation of bacterial transcriptomes.

Comparison of clinical targeted next-generation sequence data from formalin-fixed and fresh-frozen tissue specimens

Next-generation sequencing (NGS) has emerged as a powerful technique for the detection of genetic variants in the clinical laboratory. NGS can be performed using DNA from FFPE tissue, but it is unknown whether such specimens are truly equivalent to unfixed tissue for NGS applications. To address this question, we performed hybridization-capture enrichment and multiplexed Illumina NGS for 27 cancer-related genes using DNA from 16 paired fresh-frozen and routine FFPE lung adenocarcinoma specimens and conducted extensive comparisons between the sequence data from each sample type. This analysis revealed small but detectable differences between FFPE and frozen samples. Compared with frozen samples, NGS data from FFPE samples had smaller library insert sizes, greater coverage variability, and an increase in C to T transitions that was most pronounced at CpG dinucleotides, suggesting interplay between DNA methylation and formalin-induced changes; however, the error rate, library complexity, enrichment performance, and coverage statistics were not significantly different. Comparison of base calls between paired samples demonstrated concordances of >99.99%, with 96.8% agreement in the single-nucleotide variants detected and >98% accuracy of NGS data when compared with genotypes from an orthogonal single-nucleotide polymorphism array platform. This study demonstrates that routine processing of FFPE samples has a detectable but negligible effect on NGS data and that these samples can be a reliable substrate for clinical NGS testing.

Whole genome expression profiling based on paraffin embedded tissue can be used to classify diffuse large B-cell lymphoma and predict clinical outcome

This study tested the validity of whole-genome expression profiling (GEP) using RNA from formalin-fixed, paraffin-embedded (FFPE) tissue to sub-classify Diffuse Large B-cell Lymphoma (DLBCL), in a population based cohort of 172 patients. GEP was performed using Illumina Whole Genome cDNA-mediated Annealing, Selection, extension & Ligation, and tumours were classified into germinal centre (GCB), activated B-cell (ABC) and Type-III subtypes. The method was highly reproducible and reliably classified cell lines of known phenotype. GCB and ABC subtypes were each characterized by unique gene expression signatures consistent with previously published data. A significant relationship between subtype and survival was observed, with ABC having the worst clinical outcome and in a multivariate survival model only age and GEP class remained significant. This effect was not seen when tumours were classified by immunohistochemistry. There was a significant association between age and subtype (mean ages ABC - 72·8 years, GC - 68·4 years, Type-III - 64·5 years). Older patients with ABC subtype were also over-represented in patients who died soon after diagnosis. The relationship between prognosis and subtype improved when only patients assigned to the three categories with the highest level of confidence were analysed. This study demonstrates that GEP-based classification of DLBCL can be applied to RNA extracted from routine FFPE samples and has potential for use in stratified medicine trials and clinical practice.

Sequencing intractable DNA to close microbial genomes

Advancement in high throughput DNA sequencing technologies has supported a rapid proliferation of microbial genome sequencing projects, providing the genetic blueprint for in-depth studies. Oftentimes, difficult to sequence regions in microbial genomes are ruled “intractable” resulting in a growing number of genomes with sequence gaps deposited in databases. A procedure was developed to sequence such problematic regions in the “non-contiguous finished” Desulfovibrio desulfuricans ND132 genome (6 intractable gaps) and the Desulfovibrio africanus genome (1 intractable gap). The polynucleotides surrounding each gap formed GC rich secondary structures making the regions refractory to amplification and sequencing. Strand-displacing DNA polymerases used in concert with a novel ramped PCR extension cycle supported amplification and closure of all gap regions in both genomes. The developed procedures support accurate gene annotation, and provide a step-wise method that reduces the effort required for genome finishing.

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.