Optimization for Sequencing and Analysis of Degraded FFPE-RNA Samples

Gene expression analysis by RNA sequencing (RNA-seq) enables unique insights into clinical samples that can potentially lead to mechanistic understanding ...

Frederick National Laboratory for Cancer Research

Interleukin-27 (IL-27) is a pleiotropic cytokine which plays important and diverse roles in the immune system. We have previously demonstrated that IL-27 induces potent anti-viral effects against HIV-1, HIV-2, SIV, HSV-2, KSHV and influenza viruses in macrophages. This induction occurred in an interferon (IFN) independent manner and involved down regulation of SPTBN1. MicroRNAs (miRNAs) are critical regulators of mRNA translation and turnover. There have been reports that some miRNAs inhibit viral replication. In this study, we hypothesized that IL-27 could induce the expression of novel miRNAs in macrophages which may have functional relevance in terms of anti-viral activity and primary monocytes were differentiated into macrophages using either M-CSF (M-Mac) or a combination of M-CSF and IL-27 (I-Mac) for seven days. Following this, total RNA was extracted from these cells and deep sequencing was performed, in parallel with gene expression microarrays. Using the novel miRNA discovery software, miRDeep, seven novel miRNAs were discovered in these macrophages. Four of which were preferentially expressed in I-Mac (miR-SX1, -SX2, -SX3 and -SX6) whilst three were detected in both M-Mac and I-Mac (miR-SX4, -SX5 and -SX7). The expression of six of the seven novel miRNAs was highly correlated with qRT-PCR using specific primer/probes designed for the novel miRNAs. Gene expression microarray further demonstrated that a number of genes were potentially targeted by these differentially expressed novel miRNAs. Finally, several of these novel miRNAs (miR-SX1, -SX4, -SX5, -SX6 and -SX7) were shown to target the open reading frames of a number of viruses (including HSV-1, HSV-2 and HHV-8) which may partially explain the anti-viral properties observed.

Interleukin-27 treated human macrophages induce the expression of novel microRNAs which may mediate anti-viral properties.

Although histone H3 lysine 27 trimethylation (H3K27Me3) is associated with gene silencing, whether H3K27Me3 demethylation affects transcription and cell differentiation in vivo has remained elusive. To investigate this, we conditionally inactivated the two H3K27Me3 demethylases, Jmjd3 and Utx, in non-dividing intrathymic CD4(+) T-cell precursors. Here we show that both enzymes redundantly promote H3K27Me3 removal at, and expression of, a specific subset of genes involved in terminal thymocyte differentiation, especially S1pr1, encoding a sphingosine-phosphate receptor required for thymocyte egress. Thymocyte expression of S1pr1 was not rescued in Jmjd3- and Utx-deficient male mice, which carry the catalytically inactive Utx homolog Uty, supporting the conclusion that it requires H3K27Me3 demethylase activity. These findings demonstrate that Jmjd3 and Utx are required for T-cell development, and point to a requirement for their H3K27Me3 demethylase activity in cell differentiation. 

Histone H3 Lysine 27 demethylases Jmjd3 and Utx are required for T-cell differentiation.

Pneumocystis jirovecii is a major cause of life-threatening pneumonia in immunosuppressed patients including transplant recipients and those with HIV/AIDS, yet surprisingly little is known about the biology of this fungal pathogen. Here we report near complete genome assemblies for three Pneumocystis species that infect humans, rats and mice. Pneumocystis genomes are highly compact relative to other fungi, with substantial reductions of ribosomal RNA genes, transporters, transcription factors and many metabolic pathways, but contain expansions of surface proteins, especially a unique and complex surface glycoprotein superfamily, as well as proteases and RNA processing proteins. Unexpectedly, the key fungal cell wall components chitin and outer chain N-mannans are absent, based on genome content and experimental validation. Our findings suggest that Pneumocystis has developed unique mechanisms of adaptation to life exclusively in mammalian hosts, including dependence on the lungs for gas and nutrients and highly efficient strategies to escape both host innate and acquired immune defenses.

Genome analysis of three Pneumocystis species reveals adaptation mechanisms to life exclusively in mammalian hosts.

Development of HIV-1 drug resistance mutations (HDRMs) is one of the major reasons for the clinical failure of antiretroviral therapy. Treatment success rates can be improved by applying personalized anti-HIV regimens based on a patient's HDRM profile. However, the sensitivity and specificity of the HDRM profile is limited by the methods used for detection. Sanger-based sequencing technology has traditionally been used for determining HDRM profiles at the single nucleotide variant (SNV) level, but with a sensitivity of only ≥ 20% in the HIV population of a patient. Next Generation Sequencing (NGS) technologies offer greater detection sensitivity (~ 1%) and larger scope (hundreds of samples per run). However, NGS technologies produce reads that are too short to enable the detection of the physical linkages of individual SNVs across the haplotype of each HIV strain present. In this article, we demonstrate that the single-molecule long reads generated using the Third Generation Sequencer (TGS), PacBio RS II, along with the appropriate bioinformatics analysis method, can resolve the HDRM profile at a more advanced quasispecies level. The case studies on patients' HIV samples showed that the quasispecies view produced using the PacBio method offered greater detection sensitivity and was more comprehensive for understanding HDRM situations, which is complement to both Sanger and NGS technologies. In conclusion, the PacBio method, providing a promising new quasispecies level of HDRM profiling, may effect an important change in the field of HIV drug resistance research.

Towards Better Precision Medicine: PacBio Single-Molecule Long Reads Resolve the Interpretation of HIV Drug Resistant Mutation Profiles at Explicit Quasispecies (Haplotype) Level.

Gene isoforms are commonly found in both prokaryotes and eukaryotes. Since each isoform may perform a specific function in response to changing environmental conditions, studying the dynamics of gene isoforms is important in understanding biological processes and disease conditions. However, genome-wide identification of gene isoforms is technically challenging due to the high degree of sequence identity among isoforms. Traditional targeted sequencing approach, involving Sanger sequencing of plasmid-cloned PCR products, has low throughput and is very tedious and time-consuming. Next-generation sequencing technologies such as Illumina and 454 achieve high throughput but their short read lengths are a critical barrier to accurate assembly of highly similar gene isoforms, and may result in ambiguities and false joining during sequence assembly. More recently, the third generation sequencer represented by the PacBio platform offers sufficient throughput and long reads covering the full length of typical genes, thus providing a potential to reliably profile gene isoforms. However, the PacBio long reads are error-prone and cannot be effectively analyzed by traditional assembly programs.

Distinguishing highly similar gene isoforms with a clustering-based bioinformatics analysis of PacBio single-molecule long reads.

Intratumor molecular heterogeneity of hepatocellular carcinoma is partly attributed to the presence of hepatic cancer stem cells (CSCs). Different CSC populations defined by various cell surface markers may contain different oncogenic drivers, posing a challenge in defining molecularly targeted therapeutics. We combined transcriptomic and functional analyses of hepatocellular carcinoma cells at the single-cell level to assess the degree of CSC heterogeneity. We provide evidence that hepatic CSCs at the single-cell level are phenotypically, functionally, and transcriptomically heterogeneous. We found that different CSC subpopulations contain distinct molecular signatures. Interestingly, distinct genes within different CSC subpopulations are independently associated with hepatocellular carcinoma prognosis, suggesting that a diverse hepatic CSC transcriptome affects intratumor heterogeneity and tumor progression.

Single-cell analysis reveals cancer stem cell heterogeneity in hepatocellular carcinoma.

Diffuse large B-cell lymphomas (DLBCLs) are phenotypically and genetically heterogeneous. Gene-expression profiling has identified subgroups of DLBCL (activated B-cell-like [ABC], germinal-center B-cell-like [GCB], and unclassified) according to cell of origin that are associated with a differential response to chemotherapy and targeted agents. We sought to extend these findings by identifying genetic subtypes of DLBCL based on shared genomic abnormalities and to uncover therapeutic vulnerabilities based on tumor genetics.

Genetics and Pathogenesis of Diffuse Large B-Cell Lymphoma.

Memory CD4 T cells mediate long-term immunity, and their generation is a key objective of vaccination strategies. However, the transcriptional circuitry controlling the emergence of memory cells from early CD4 antigen-responders remains poorly understood. Here, using single-cell RNA-seq to study the transcriptome of virus-specific CD4 T cells, we identified a gene signature that distinguishes potential memory precursors from effector cells. We found that both that signature and the emergence of memory CD4 T cells required the transcription factor Thpok. We further demonstrated that Thpok cell-intrinsically protected memory cells from a dysfunctional, effector-like transcriptional program, similar to but distinct from the exhaustion pattern of cells responding to chronic infection. Mechanistically, Thpok- bound genes encoding the transcription factors Blimp1 and Runx3 and acted by antagonizing their expression. Thus, a Thpok-dependent circuitry promotes both memory CD4 T cells' differentiation and functional fitness, two previously unconnected critical attributes of adaptive immunity.

The Emergence and Functional Fitness of Memory CD4 T Cells Require the Transcription Factor Thpok.

 derived cell lines are commonly used to enable recombinant protein expression via baculovirus infection to generate materials approved for clinical use and in clinical trials. In order to develop systems biology and genome engineering tools to improve protein expression in this host, we performed de novo genome assembly of the -derived cell line Tni-FNL.

Genome Assembly and Annotation of the Trichoplusia ni Tni-FNL Insect Cell Line Enabled by Long-Read Technologies.

Formalin-fixed paraffin-embedded (FFPE) tissues are among the most widely available clinical specimens. Their potential utility as a source of RNA for transcriptome studies would greatly enhance population-based cancer studies. Although preliminary studies suggest FFPE tissue may be used for RNA sequencing, the effect of storage time on these specimens needs to be determined. We conducted this study to determine whether RNA in archived FFPE high-grade ovarian serous adenocarcinomas from Surveillance, Epidemiology and End Results (SEER) registries was present in sufficient quantity and quality for RNA-Seq analysis. FFPE tissues, stored from 7 to 32 years, were obtained from three SEER sites. RNA was extracted, quantified, quality assessed, and subjected to RNA-Seq (a whole transcriptome sequencing technology). FFPE specimens stored for longer periods of time had poorer RNA sample quality as indicated by negative correlations between specimen storage time and fragment distribution values (DV). In addition, sample contamination was a common issue among the RNA, with 41 of 67 samples having 5% to 48% bacterial contamination. However, regardless of specimen storage time and bacterial contamination, 60% of the samples yielded data that enabled gene expression quantification, identifying more than 10,000 genes, with the correlations among most biological replicates above 0.7. This study demonstrates that FFPE high-grade ovarian serous adenocarcinomas specimens stored in repositories for up to 32 years and under varying storage conditions are a promising source of RNA for RNA-Seq. We also describe certain caveats to be considered when designing RNA-Seq studies using archived FFPE tissues.

Robustness of RNA sequencing on older formalin-fixed paraffin-embedded tissue from high-grade ovarian serous adenocarcinomas.

Herpes simplex virus type 1 (HSV-1) strain MacIntyre has a severe defect in the anterograde spread after replication in the nucleus. To better understand and identify the genetic determinants that lead to the unique phenotypes of the MacIntyre strain, we sequenced its genome with PacBio single-molecule real-time sequencing technology and resolved the complete sequence.

Complete Genome Sequence of Herpes Simplex Virus 1 Strain MacIntyre.

Herpes simplex virus 1 (HSV-1) strain McKrae is highly virulent and relatively neuroinvasive in animal models compared with other wild-type HSV-1 strains. To identify the genetic determinants that lead to the unique phenotypes of the McKrae strain, we sequenced its genome with PacBio single-molecule real-time (SMRT) technology and resolved the complete sequence.

Complete Genome Sequence of Herpes Simplex Virus 1 Strain McKrae.

Cellular diversity in tumors is a key factor for therapeutic failures and lethal outcomes of solid malignancies. Here, we determined the single-cell transcriptomic landscape of liver cancer biospecimens from 19 patients. We found varying degrees of heterogeneity in malignant cells within and between tumors and diverse landscapes of tumor microenvironment (TME). Strikingly, tumors with higher transcriptomic diversity were associated with patient's worse overall survival. We found a link between hypoxia-dependent vascular endothelial growth factor expression in tumor diversity and TME polarization. Moreover, T cells from higher heterogeneous tumors showed lower cytolytic activities. Consistent results were found using bulk genomic and transcriptomic profiles of 765 liver tumors. Our results offer insight into the diverse ecosystem of liver cancer and its impact on patient prognosis.

Tumor Cell Biodiversity Drives Microenvironmental Reprogramming in Liver Cancer.

Most current tumor immunotherapy strategies leverage cytotoxic CD8 T cells. Despite evidence for clinical potential of CD4 tumor-infiltrating lymphocytes (TILs), their functional diversity limits our ability to harness their activity. Here, we use single-cell mRNA sequencing to analyze the response of tumor-specific CD4 TILs and draining lymph node (dLN) T cells. Computational approaches to characterize subpopulations identify TIL transcriptomic patterns strikingly distinct from acute and chronic anti-viral responses and dominated by diversity among T-bet-expressing T helper type 1 (Th1)-like cells. In contrast, the dLN response includes T follicular helper (Tfh) cells but lacks Th1 cells. We identify a type I interferon-driven signature in Th1-like TILs and show that it is found in human cancers, in which it is negatively associated with response to checkpoint therapy. Our study provides a proof-of-concept methodology to characterize tumor-specific CD4 T cell effector programs. Targeting these programs should help improve immunotherapy strategies.

Bao Tran

NCI CCR Sequencing Facility

Frederick National Laboratory for Cancer Research

Optimization for Sequencing and Analysis of Degraded FFPE-RNA Samples

Bao Tran

.css-o250i3{font-size:18px;font-family:Open Sans,sans-serif;line-height:21.6px;}NCI CCR Sequencing Facility

Frederick National Laboratory for Cancer Research

Optimization for Sequencing and Analysis of Degraded FFPE-RNA Samples

NCI CCR Sequencing Facility