Preparation of archival formalin-fixed paraffin-embedded mouse liver samples for use with the Agilent gene expression microarray platform

A Streamlined and High-Throughput Error-Corrected Next-Generation Sequencing Method for Low Variant Allele Frequency Quantitation

Quantifying mutant or variable allele frequencies (VAFs) of ≤10-3 using next-generation sequencing (NGS) has utility in both clinical and nonclinical settings. Two common approaches for quantifying VAFs using NGS are tagged single-strand sequencing and duplex sequencing. While duplex sequencing is reported to have sensitivity up to 10-8 VAF, it is not a quick, easy, or inexpensive method. We report a method for quantifying VAFs that are ≥10-4 that is as easy and quick for processing samples as standard sequencing kits, yet less expensive than the kits. The method was developed using PCR fragment-based VAFs of Kras codon 12 in log10 increments from 10-5 to 10-1, then applied and tested on native genomic DNA. For both sources of DNA, there is a proportional increase in the observed VAF to input VAF from 10-4 to 100% mutant samples. Variability of quantitation was evaluated within experimental replicates and shown to be consistent across sample preparations. The error at each successive base read was evaluated to determine if there is a limit of read length for quantitation of ≥10-4, and it was determined that read lengths up to 70 bases are reliable for quantitation. The method described here is adaptable to various oncogene or tumor suppressor gene targets, with the potential to implement multiplexing at the initial tagging step. While easy to perform manually, it is also suited for robotic handling and batch processing of samples, facilitating detection and quantitation of genetic carcinogenic biomarkers before tumor formation or in normal-appearing tissue.

Bronchoalveolar lavage (BAL) cells in idiopathic pulmonary fibrosis express a complex pro-inflammatory, pro-repair, angiogenic activation pattern, likely associated with macrophage iron accumulation

Idiopathic pulmonary fibrosis (IPF) is a chronic lung disease of unknown cause characterized by alveolar epithelial damage, patchy interstitial fibrosis and diffuse microvascular abnormalities. In IPF, alveolar clustering of iron-laden alveolar macrophages—a common sign of microhemorrhage, has been associated with vascular abnormalities and worsening of pulmonary hypertension. As iron-dependent ROS generation has been shown to induce unrestrained macrophage activation in disease models of vascular damage, we explored alveolar macrophage activation phenotype in IPF patients (n = 16) and healthy controls (CTR, n = 7) by RNA sequencing of bronchoalveolar lavage (BAL) cells. The frequencies of macrophages in BAL cells were 86+4% and 83.4+8% in IPF and CTR groups, respectively (p-value = 0.41). In IPF patients, BAL cells showed increased iron-dependent ROS generation (p-value<0.05 vs CTR). Gene expression analysis showed overrepresentation of Gene Ontology processes/functions and KEGG pathways enriched in upregulated M1-type inflammatory (p-value<0.01), M2-type anti-inflammatory/tissue remodeling (p-value<0.0001), and MTPP-type chronic inflammatory/angiogenic (p-value<0.0001) chemokine and cytokine genes. The ex vivo finding was confirmed by the induction of iron-dependent ROS generation and chemokine/cytokine overexpression of Ccl4, Cxcl10 (M1), Il1rn (M2), Cxcl2, and Cxcl7 (MTPP) in MH-S murine immortalized alveolar macrophages exposed to ferric ammonium citrate in culture (p-value<0.05 vs CTR). The data show alveolar macrophage expression of a pro-inflammatory, tissue remodeling and angiogenic complex activation pattern, suggesting that iron accumulation may play a role in macrophage activation.

Toxicogenomics Case Study: Furan

Development of pragmatic methodologies for human health risk assessment is required to address current regulatory challenges. We applied three toxicogenomic approaches—quantitative, predictive, and mechanistic—to a case study in mice exposed for 3 weeks to the hepatocarcinogen furan. We modeled the dose response of a variety of transcriptional endpoints and found that they produced benchmark doses similar to the furan-dependent cancer benchmark doses. Meta-analyses showed strong similarity between furan-dependent gene expression changes and those associated with several hepatic pathologies. Molecular pathways facilitated the development of a molecular mode of action for furan-induced hepatocellular carcinogenicity. Finally, we compared transcriptomic profiles derived from formalin-fixed and paraffin-embedded (FFPE) samples with those from high-quality frozen samples to evaluate whether archival samples are a viable option for toxicogenomic studies. The advantage of using FFPE tissues is that they are very well characterized (phenotypically); the disadvantage is that formalin degrades biomacromolecules, including RNA. We found that FFPE samples can be used for toxicogenomics using a ribo-depletion RNA-seq protocol. Our case study demonstrates the utility of toxicogenomics data to human health risk assessment, the potential of archival FFPE tissue samples, and identifies viable strategies toward the reduction of animal usage in chemical testing.

Mining the Archives: A Cross-Platform Analysis of Gene Expression Profiles in Archival Formalin-Fixed Paraffin-Embedded Tissues

Formalin-fixed paraffin-embedded (FFPE) tissue samples represent a potentially invaluable resource for transcriptomic research. However, use of FFPE samples in genomic studies has been limited by technical challenges resulting from nucleic acid degradation. Here we evaluated gene expression profiles derived from fresh-frozen (FRO) and FFPE mouse liver tissues preserved in formalin for different amounts of time using 2 DNA microarray protocols and 2 whole-transcriptome sequencing (RNA-seq) library preparation methodologies. The ribo-depletion protocol outperformed the other methods by having the highest correlations of differentially expressed genes (DEGs), and best overlap of pathways, between FRO and FFPE groups. The effect of sample time in formalin (18 h or 3 weeks) on gene expression profiles indicated that test article treatment, not preservation method, was the main driver of gene expression profiles. Meta- and pathway analyses indicated that biological responses were generally consistent for 18 h and 3 week FFPE samples compared with FRO samples. However, clear erosion of signal intensity with time in formalin was evident, and DEG numbers differed by platform and preservation method. Lastly, we investigated the effect of time in paraffin on genomic profiles. Ribo-depletion RNA-seq analysis of 8-, 19-, and 26-year-old control blocks resulted in comparable quality metrics, including expected distributions of mapped reads to exonic, untranslated region, intronic, and ribosomal fractions of the transcriptome. Overall, our results indicate that FFPE samples are appropriate for use in genomic studies in which frozen samples are not available, and that ribo-depletion RNA-seq is the preferred method for this type of analysis in archival and long-aged FFPE samples.