Ng Of Total RNA Research Articles

Abstract 745: Development of a comprehensive gene fusion NGS panel using an integrated microfluidic circuit enabling highly efficient, multiplex PCR enrichment

Abstract Introduction: RNA fusion transcripts result from genomic rearrangements where two distinct genes become juxtaposed and fused. Fusion transcripts are an important class of somatic alterations because they have the potential to create chimeric proteins with altered function, contributing to oncogene activation. These aberrant proteins are expressed in a tumor-specific manner and are thus excellent targets for therapeutic intervention. We have designed an amplicon-based library preparation (LP) kit that targets over 350 fusion gene pairs, representing over 1,000 unique fusion breakpoint events from both solid tumor and hematologic cancers. This panel is specifically designed to run on a newly developed microfluidic integrated fluidic circuit (IFC) that supports flexible amplicon library preparation of 1-6 unique LP panels on a single IFC. Samples can be run against the LP panels in groups of 8, allowing for up to 48 samples per IFC. In this poster we present the results of our analytical validation tests, demonstrating the performance of our comprehensive gene fusion panel on synthetic targets, reference FFPE samples, and subject FFPE samples for use in research studies. Methods: Primers for amplicon library generation targeting known fusion breakpoints were designed in collaboration with Q2 Solutions®. The panel is comprised of over 1,000 amplicons with an average insert size of 170 bp. Primers were divided into 8 assay pools using informatic parameters that minimize generation of off-target products. Assay pools and cDNA from an off-IFC reverse transcription reaction, along with sample barcodes, were dispensed into designated inlets on an IFC and placed in a Juno™ targeted DNA sequencing library preparation system for mixing the nanoscale reactions, thermal cycling, and amplicon harvesting. Harvested amplicons were collected from the IFC, pooled, and prepared for sequencing on an Illumina® NextSeq™. Results: The panel was tested against synthetic targets representing all targeted amplicons to demonstrate that each assay reliably detects the intended target of interest. Reference standards and subject FFPE samples were used to assess performance, including positive predictive agreement (PPA), positive predictive value (PPV), and limit of detection. Using these samples at inputs as low as 10 ng of total RNA, PPA and PPV were greater than 99% for the fusion events tested. In addition, fusion events could be reliably detected with as little as 250 copies of target fusion material loaded onto the IFC. Conclusions: A comprehensive fusion gene panel for targeted next-generation sequencing has been developed using nanoliter-scale PCR based enrichment on a newly developed IFC. Generation of high-quality libraries from a minimum of 10 ng total RNA from FFPE samples for use in research studies has been demonstrated. Citation Format: Michael L. Gonzales, Jian Qin, Xiaohui Wang, Sangpen Chamnongpol, Jeff Jasper, Thomas Goralski, Christopher Kubu. Development of a comprehensive gene fusion NGS panel using an integrated microfluidic circuit enabling highly efficient, multiplex PCR enrichment [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 745.

Janus kinase 1/2 inhibition for the treatment of autoinflammation associated with heterozygous TNFAIP3 mutation

Janus kinase 1/2 inhibition for the treatment of autoinflammation associated with heterozygous TNFAIP3 mutation

A comparative analysis of RNA sequencing methods with ribosome RNA depletion for degraded and low-input total RNA from formalin-fixed and paraffin-embedded samples.

e14654 Background: Formalin-fixed and paraffin-embedded (FFPE) blocks held in clinical laboratories are an invaluable resource for clinical research, especially in the era of personalized medicine. With the development of massive parallel sequencing, RNA expression profiles have become an important source of new biomarkers with potential values in cancer metastasis and disease prognosis. Standard practice in tissue fixing and paraffin-embedding has little impact on the expression analysis of RNA samples, which makes archived FFPE samples valuable and retrospective studies feasible. But it is still difficult for some FFPE samples which have low quality or low amount of RNA to be used in the RNA expression study. Methods: We compared four FFPE RNA library preparation kits (KAPA, TaKaRa, QIAGEN and Vazyme) based on two principles of rRNA depletion, with total RNA from FFPE samples and paired FF samples using two bioinformatics analysis methods. Furthermore, we tried to construct RNA libraries with the two kits allowing for degraded low-input RNA using clinical FFPE samples of different quality. HISAT was used to analyze the data of all the libraries. Results: Our results demonstrated that high concordance in transcript quantifications were shown between FF and FFPE samples using the same kit and more transcripts could be obtained by rRNA depletion after reverse transcript reaction. The gene expression profiling using the TaKaRa kit showed a difference with other kits, which may be due to the different principle of rRNA depletion or the amount of input total RNA. For seriously degraded RNA from FFPE samples, libraries could be constructed with as low as 50 ng of total RNA, although there was higher residual rRNA and lower unique mapping ratio in the libraries. Conclusions: These results indicated that most of clinical FFPE samples could be used in the RNA profiling study. The method of RNA library construction with rRNA depletion can be used for clinical FFPE samples. For degraded and low-input RNA samples, it is still possible to obtain repeatable RNA expression profiling but with a low unique mapping ratio and high residual rRNA.

Androgen induces olfactory expression of prostaglandin E2 receptor Ep1 in the burrow-living fish Bostrychus sinensis

It is well documented that androgens modify olfactory processing in vertebrates. In fish, several lines of evidence indicate that androgens increase olfactory sensitivity to prostaglandin pheromone, but the molecular mechanism is still unclear. Our previous studies showed that prostaglandin E2 (PGE2) is a sex pheromone in the burrowing-living fish Chinese black sleeper (Bostrychus sinensis) and that the PGE2 receptor 1 (Ep1) in the olfactory rosette is a candidate receptor for sensing sex pheromone PGE2. In the present study, we found that testosterone (T) and 11-ketotestosterone (11-KT) exhibited stimulatory effects on the expression of ep1 in the olfactory rosette in vivo and ex vivo. Moreover, the androgen receptor (Ar) agonist R1881 had similar effects to 11-KT on the expression of ep1 ex vivo, suggesting the up-regulatory effect is mediated by Ar. The amount of arα transcripts (˜1500 copies/100 ng total RNA) was greater than that of arβ (˜300 copies/100 ng total RNA) in the olfactory rosette, and the expression levels of arα increased with spermatogenesis and peaked at late meiosis stage. Moreover, activated Arα but not Arβ transactivated a 2k bp ep1 promoter in HEK293T cell, and some OSNs exhibited co-localization of arα mRNA and Ep1 protein signals. Taken together, our results suggest that Arα, but not Arβ, plays a crucial role in mediating the androgen-induced up-regulation of ep1 expression in B. sinensis. The present study is the first to shed light on the molecular mechanisms whereby androgens enhance responsiveness to prostaglandin sex pheromones in teleosts.

Developmental validation of the ParaDNA® Body Fluid ID System—A rapid multiplex mRNA-profiling system for the forensic identification of body fluids

Identifying the biological origin of forensic traces can provide crucial evidence to aid criminal investigations. Current forensic practice for the identification of body fluids mostly uses protein-based presumptive tests. Such tests cannot identify all of the forensically relevant fluids and have issues of cross-reactivity. More recently, messenger RNA methods have been developed that have expanded the range of body fluids that can be positively identified. However, these methods are slow and require expert scientists to run the processes and interpret the results. The ParaDNA® Body Fluid ID System has been designed to provide a simple, fast and robust way to identify forensically relevant body fluids in a lab or field-deployable manner. The system can analyse and detect mRNA targets for six different body fluids: vaginal fluid, seminal fluid, sperm cells, saliva, menstrual blood and peripheral blood. Utilising the ParaDNA Sample Collector and instruments, minimal training is required to enable both forensic scientists and non-specialist personnel to analyse biological samples directly, without prior sample processing, in approximately 90 min. The developmental validation studies described here were designed to address requirements of end users, based on the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines, and tested the sensitivity, reproducibility, accuracy, inhibitor tolerance, and performance of the ParaDNA Body Fluid ID System on a range of mock evidence items. The data collected demonstrate that the Body Fluid ID System can automatically determine the presence of specific body fluid mRNA markers in single-source and mixed samples on multiple substrate types and body fluids could be identified with as little as 0.05ng total RNA and 1μl of the relevant fluid. Results can either be used to support confirmation of source from previously obtained STR DNA profiling results or to improve sample success rates by making better informed evidential submissions.

Abstract 1492: High-efficiency, strand-specific, low-bias library preparation method for transcriptome profiling of low-input RNA

Abstract RNA Sequencing (RNA-seq) is a powerful technique for transcriptome analysis. It is widely used for gene expression analysis, detection of mutations, fusion transcripts, alternative splicing, and the study of post transcriptional modifications. RNA-seq is emerging as a powerful tool for molecular diagnostics as well, providing insights into the biological pathways and molecular mechanisms associated with disease progression including cancer. Recent improvements in next generation sequencing technologies (NGS) and sample barcoding strategies have resulted in significant cost reduction. As RNA-seq is increasingly adopted for molecular diagnostics, the quality and reproducibility of library preparation methods become more important. In addition, demand for library preparation methods that produce successful NGS libraries from very low input RNA or precious clinical samples is increasing. To achieve these goals, we have developed an RNA-seq library preparation method which can be used across a wide range of input RNA (5-1000 ng total RNA). GC content analysis, gene body coverage and gene expression correlation show that these important parameters remain consistent across varying inputs, even though input amounts vary by two orders of magnitude. As a result, our method has increased sensitivity and specificity for low-abundance transcripts, and reduced PCR duplicates and sequence bias, delivering high quality strand-specific data. Our method is compatible with poly A-tail enriched and ribosomal RNA depleted samples as well as RNA extracted from FFPE (formaldehyde fixed paraffin embedded) tissue samples. Citation Format: Anetta Nowosielska-Vecchietti. High-efficiency, strand-specific, low-bias library preparation method for transcriptome profiling of low-input RNA [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1492.

Differential Expression of Coding and Long Noncoding RNAs in Keratoconus-Affected Corneas.

PurposeKeratoconus (KC) is the most common corneal ectasia. We aimed to determine the differential expression of coding and long noncoding RNAs (lncRNAs) in human corneas affected with KC.MethodsFrom the corneas of 10 KC patients and 8 non-KC healthy controls, 200 ng total RNA was used to prepare sequencing libraries with the SMARTer Stranded RNA-Seq kit after ribosomal RNA depletion, followed by paired-end 50-bp sequencing with Illumina Sequencer. Differential analysis was done using TopHat/Cufflinks with a gene file from Ensembl and a lncRNA file from NONCODE. Pathway analysis was performed using WebGestalt. Using the expression level of differentially expressed coding and noncoding RNAs in each sample, we correlated their expression levels in KC and controls separately and identified significantly different correlations in KC against controls followed by visualization using Cytoscape.ResultsUsing |fold change| ≥ 2 and a false discovery rate ≤ 0.05, we identified 436 coding RNAs and 584 lncRNAs with differential expression in the KC-affected corneas. Pathway analysis indicated the enrichment of genes involved in extracellular matrix, protein binding, glycosaminoglycan binding, and cell migration. Our correlation analysis identified 296 pairs of significant KC-specific correlations containing 117 coding genes enriched in functions related to cell migration/motility, extracellular space, cytokine response, and cell adhesion. Our study highlighted the potential roles of several genes (CTGF, SFRP1, AQP5, lnc-WNT4-2:1, and lnc-ALDH3A2-2:1) and pathways (TGF-β, WNT signaling, and PI3K/AKT pathways) in KC pathogenesis.ConclusionsOur RNA-Seq–based differential expression and correlation analyses have identified many potential KC contributing coding and noncoding RNAs.

Effect of different reverse transcription aproaches in Pru p 3 transcripts semiquantitative amplification

Reverse transcriptase transcribes the cDNA based on its previous extraction and standardization. Reverse transcription step is considered to be critical in the workflow of quantification of transcribed genes. The aim of the study was to extract total RNA by different methods and to analyse the results of the subsequent reverse transcription reaction when different commercial RT kits were used to process RNA extracted from pulp of matured peach fruit. Mature peach pulp was used in the study. The fruit of variety Vistarich was collected in summer 2017 in the orchard of Dvory nad Žitavou. Two RNA extraction methods, TRIzol® Reagent and GeneJET Plant RNA Purification Kit, were tested in to determine the suitable method for peach fruit RNA extraction. Three different cDNA reagent sets were used to transcribe 115 ng/500 ng total RNA or 11 ng/115 ng, respectively. Both variants of the primers, random hexamers as well as oligo (dT) 18, were used to anneal the target mRNA of Pru p 3 allergen following the manufacturer instructions. No specific effect was obtained in the case of peach fruit when using ethanol-extracted tissue treatment and the effect of the used extraction method was more significant. The A260/230 ratios were similar for three from four tested methods. In the case of these three methods, the A260/A230 ratios for all the extracted samples were higher than 1.9 which indicates high purity without contamination by polyphenols and polysaccharides. The specificity of obtained amplicons was proved by restriction cleavage using Tse I restriction endonuclease. This provided the cleavage of the 179 bp long product in all amplicons. Working with mature fruit meet a specific situation in the field of RNA extraction and subsequently all the downstream applications. That is, why choosing the most fitting methods and kits is a crucial step. Here, the method for the semi-quantitative analysis of the Pru p 3 allergen expressions was set up in the way that will be directly applicable for Pru p 3 expression analyses.

Reverse transcription loop-mediated isothermal amplification for species-specific detection of tomato chlorotic spot orthotospovirus

Tomato chlorotic spot orthotospovirus (TCSV) is an emerging orthotospovirus that can cause severe disease on tomato plants. There are at least four orthotospoviruses infecting tomato, and mixed infection of two or more orthotospoviruses in a single tomato plant is quite common in the field. With similarity in the symptomatology and cross serological reactivity among tomato-infecting orthotospoviruses, especially between TCSV and groundnut ringspot orthotospovirus (GRSV), the current serological tests could not achieve definite and accurate species-specific determination in disease diagnosis. Here, a one-step reverse transcription loop-mediated isothermal amplification (RT-LAMP) assay was developed for TCSV. Under optimum conditions, the virus was detected in as little as 0.2 ng of total RNA or in 1:10,000 dilution of a simple diluted tissue extract, which was ten times more sensitive than a conventional RT-PCR assay. The RT-LAMP assay was highly specific for TCSV, with no cross reaction with the other two orthotospoviruses: GRSV and tomato spotted wilt orthotospovirus (TSWV). These results demonstrate that this simple and sensitive RT-LAMP could be used to achieve species-specific detection for TCSV under field conditions.

Differentiation potential of Pluripotent Stem Cells correlates to the level of CHD7

Embryonic Stem Cells (ESC) possesses two distinct features; self-renewal and the potential to differentiate. Here we show the differentiation potential and growth rate of ESC correlates positively with the expression level of the gene encoding chromodomain helicase DNA binding protein 7 (CHD7). When ESCs are maintained in feeder-free conditions and single cell seeding, ESC KhES-1 having 4520 copies or more of CHD7 in 5 ng total RNA show differentiation potential, but this is lost when the CHD7 copy number is reduced in KhES-1 to less than 696 by alternative culture conditions. Introduction of siCHD7 reduced differentiation potential and growth rate of KhES-1. Interestingly, KhES-1 underwent spontaneous differentiation when mCHD7 was introduced and we could not obtain CHD7-overexpressing ESC in culture. These data suggest that CHD7 drives differentiation, and there is a lower limit for CHD7 to initiate differentiation and an upper limit for CHD7 if maintained in undifferentiated state, and such upper limit varies depending on culture condition. As CHD7 drives cell growth, ESC with the highest permissible CHD7 level in the given culture become dominant in a couple of passages. Thus, we can select differentiation resistance-free cell clones by optimizing the culture system using CHD7 as an index.

A comprehensive molecular and morphological study of the effects of space flight on human capillary endothelial cells: sample quality assessment and preliminary results.

ESA (ESA ILSRA-2009-1026); ASI (contract no. 5681); Regione Toscana (POR FSE 2007-2013-FORTEC); Kayser Italia; Lions Club International, District 108LA, Toscana, Italy.

Trypanosoma cruzi specific mRNA amplification by in vitro transcription improves parasite transcriptomics in host-parasite RNA mixtures

BackgroundTrypanosomatids are a group of protozoan parasites that includes the etiologic agents of important human illnesses as Chagas disease, sleeping sickness and leishmaniasis. These parasites have a significant distinction from other eukaryotes concerning mRNA structure, since all mature mRNAs have an identical species-specific sequence of 39 nucleotides at the 5′ extremity, named spliced leader (SL). Considering this peculiar aspect of trypanosomatid mRNA, the aim of the present work was to develop a Trypanosoma cruzi specific in vitro transcription (IVT) linear mRNA amplification method in order to improve parasite transcriptomics analyses.MethodsWe designed an oligonucleotide complementary to the last 21 bases of T. cruzi SL sequence, bearing an upstream T7 promoter (T7SL primer), which was used to direct the synthesis of second-strand cDNA. Original mRNA was then amplified by IVT using T7 RNA polymerase. T7SL-amplified RNA from two distinct T. cruzi stages (epimastigotes and trypomastigotes) were deep sequenced in SOLiD platform. Usual poly(A) + RNA and and T7-oligo(dT) amplified RNA (Eberwine method) were also sequenced. RNA-Seq reads were aligned to our new and improved T. cruzi Dm28c genome assembly (PacBio technology) and resulting transcriptome pattern from these three RNA preparation methods were compared, mainly concerning the conservation of mRNA transcritional levels and DEGs detection between epimastigotes and trypomastigotes.ResultsT7SL IVT method detected more potential differentially expressed genes in comparison to either poly(A) + RNA or T7dT IVT, and was also able to produce reliable quantifications of the parasite transcriptome down to 3 ng of total RNA. Furthermore, amplification of parasite mRNA in HeLa/epimastigote RNA mixtures showed that T7SL IVT generates transcriptome quantification with similar detection of differentially expressed genes when parasite RNA mass was only 0.1% of the total mixture (R = 0.78 when compared to poly(A) + RNA).ConclusionsThe T7SL IVT amplification method presented here allows the detection of more potential parasite differentially expressed genes (in comparison to poly(A) + RNA) in host-parasite mixtures or samples with low amount of RNA. This method is especially useful for trypanosomatid transcriptomics because it produces less bias than PCR-based mRNA amplification. Additionally, by simply changing the complementary region of the T7SL primer, the present method can be applied to any trypanosomatid species.

Abstract 5409: Assessment of the MammaPrint 70-gene profile using RNA sequencing technology

Abstract Introduction: Improvements in RNA processing have enabled microarray diagnostics for formalin-fixed, paraffin-embedded (FFPE) tissue. Recently, MammaPrint, a prognostic 70-gene profile for early-stage breast cancer, was successfully translated to FFPE tissue showing to be substantially equivalent to fresh tissue. In recent years, RNA-sequencing (RNA-Seq) became the standard method for transcriptome analysis, because of its low background signal and its ability of quantifying a large dynamic range of expression levels. Here we report a preliminary analysis of the FFPE MammaPrint 70-gene profile using RNA-Seq technology and the comparison with the MammaPrint® microarray diagnostic test in a series of FFPE samples. Methods: RNA-Seq was carried out using a strand-specific RNA library preparation followed by target enrichment of the coding region of the human transcriptome without relying on the presence of poly-A tail. RNA sequencing libraries were prepared starting from a minimal amount of 20 ng of total RNA based on the DV200 metric assessment. The library pools were single-end sequenced on the Illumina HiSeq 2500 instrument at the length of 65bp. The resulting sequences were mapped to the human reference genome (build 38) using TopHat v2.1. Tophat was guided by using a transcriptome index from Ensembl (version 77). The HTSeq-count tool was used to generate the total number of uniquely mapped reads for each gene. Gene expressions were normalized with Count Per Million (CPM) normalization and log2 transformed afterwards. Microarray data of the sample were available for analysis comparison. Results: On average, we obtained 22 million reads assigned to gene per sample (min=15M, max=28M). The number of reads assigned to genes vary from 61% to 70% of the total number of reads. Between 80% and 90% of the reads assigned to genes mapped to protein coding genes which is comparable to fresh frozen material. The 70-gene signature was successfully mapped to the RNA-Seq genes. A median raw read-count of 384 was observed for the 70-gene profile among the samples. Importantly, we observed a high concordance (R2 Pearson correlation=0.97) between the MammaPrint index calculated using the RNA-Seq data and the correspondent Microarray MammaPrint index. Additionally, the BluePrint profile, a microarray diagnostic test for breast cancer molecular subtyping, was successfully translated to the RNA-Seq platform. As with the MammaPrint profile, BluePrint showed high concordance between the two technologies with high correlation values for each of the subtypes (Luminal R2 Pearson correlation=0.98, Basal R2 Pearson correlation=0.97, HER2 R2 Pearson correlation=0.77). Conclusions: Next Generation RNA-sequencing is a feasible technology to assess diagnostic signatures, such as the 70 gene MammaPrint and BluePrint profiles. Citation Format: Lorenza Mittempergher, Jacob B. Spangler, Mireille H. Snel, Leonie J. Delahaye, Iris de Rink, Sun Tian, Annuska M. Glas, Rene Bernards. Assessment of the MammaPrint 70-gene profile using RNA sequencing technology [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5409. doi:10.1158/1538-7445.AM2017-5409

Abstract 3994: Using NGS to characterize 40 NSCLC tumor with gene expression and tumor infiltrating T cell repertoire profiling from FFPE and matched Fresh Frozen samples

Abstract The tumor microenvironment (TME) is made up of stromal cells, immune cells, signaling molecules, and blood vessels surrounding the tumor cells. It has emerged as a key factor in multiple stages of cancer progression, immune-escaping, and disease progression. The composition and activity of TME co-evolve with tumor cells and may likely affect how the cancer responses to immunotherapy. A clear understanding of the functional effect and evolution of the TME necessitates a comprehensive approach to identify key immune cells as well as to characterize signaling and inflammatory (immune-editing) activities at the tumor site. Here we are describing the use of a targeted NGS panel to detect aggregated gene expression and a novel AmpliSeqTM approach to profile the relative abundance of different T cell clones at the TME. Specifically, we measured the expression of 395 relevant genes that capture interferon and chemokine signaling, T and B cell activation, checkpoint pathway, tumor proliferation, and antigen presentation. By looking at the expression of markers specific to different effector cell types, this gene panel offers a high-level view of the composition of different lymphocyte infiltrates. Complementarily, the TCR sequencing assay provides an estimate of T cell diversity and therefore offering a different dimension of the immune response. We studied 40 NSCLC tumor samples, of which we have matched FFPE and fresh frozen specimens. With a targeted panel, we could detect expression of transcripts present as few as 50 copies in 10 ng of total RNA as input. Across 40 NSCLC samples, we were able to measure expression of many low expressing cytokines such as IL2, IL21, IL23, IFNG and TNF. We observed strong co-expression pattern among genes involved in type II interferon signaling, indicating that they’re informative of T cell activation. More specifically, we found strong correlation between CD8 expression and other T cell co-stimulatory receptors (CD28, CD80, CD86, CD40), suggesting that expression of these genes can be reliable surrogates for the protein counterparts as markers for CD8+ T cells. TCRβ was sequenced for each of the matched fresh frozen sample in duplicates. Clonotype abundance of replicates was highly correlated with each other, indicating that the assay was reliable, and the samples were sequenced to appropriated depth. We identified 2000-10,000 unique clones in each tumor sample; with the diversity index moderately correlated with the percentage of tumor infiltrating lymphocytes provided by pathology review. Together, these two assays provide a convenient means to characterize T cells and their activities in combating tumor cells, thereby offering insights into how that tumor may respond to a particular immunotherapy. Note: This abstract was not presented at the meeting. Citation Format: Ann Mongan, Geoffrey Lowman, Yuan-Chieh Ku, Timothy Looney, Lifeng Lin, Xinzhan Peng, Fiona Hyland. Using NGS to characterize 40 NSCLC tumor with gene expression and tumor infiltrating T cell repertoire profiling from FFPE and matched Fresh Frozen samples [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3994. doi:10.1158/1538-7445.AM2017-3994

Abstract 5364: A targeted NGS solution to evaluate gene expression signature of the tumor microenvironment from 40 NSCLC FFPE and matched fresh frozen samples

Abstract Cancer cells and their surrounding non-malignant cells, including immune cells, signaling molecules, stromal and extracellular matrix, create the tumor microenvironment (TME). The composition of this TME plays important roles in tumor progression, evading growth suppressors and activating metastasis. However, the regulatory mechanism and function of each constituent remains poorly understood. With several checkpoint blockade therapy studies, the presence of PD-L1 has been reported to be a promising marker to predict positive response. Current IHC methods to measure PD-L1 are subjective and highly variable. A higher-throughput and standardized solution that can systematically measure gene expression of cells present in the TME has emerged to be a more desirable alternative. Here, we applied the OncomineTM Immune Response Research Assay to measure the expression of 395 genes in non-small cell lung cancer (NSCLC) samples from 40 matched FFPE and fresh frozen sample types. This assay leverages NGS technology to sequence and count reads derived from the original transcript. With an input requirement of 10 ng of total RNA, libraries were generated, templated on the Ion ChefTM and sequenced on the Ion S5TM System. Results showed that, despite small input amount, the expression profiles of FFPE and fresh frozen samples are highly correlated with an average correlation greater than 0.9. We selected 22 genes out of the panel to validate expression with qPCR using FFPE samples. These genes were selected to cover a range of low, medium, and high expressors per our NGS data. Again, we observed a strong correlation (R ~ 0.9) between NGS and qPCR data. Approximately 80% of the 40 samples show moderate to high expression of CD8+ T cell cytokines, IFNG and TNFa. We further found that the expression of CD8A and CD8B are highly correlated with CD4, suggesting the co-presence of both cytotoxic and helper T cells. High correlation among CD4, FOXP3, TGFB1, and IL2RA (CD25) also suggests that their expression can be used as markers for the presence of Treg cells. We conducted a differential expression analysis between a group of samples (n=8) with high percentage of surrounding and infiltrating lymphocytes and another group (n=5) with low stromal content but devoid of infiltrating lymphocytes. Interestingly, we found a large number of genes which annotated as markers for infiltrating lymphocytes (CTSS, CXCR4, CD37, SRGN, FCER1G, SAMHD1, and GZMA) are significantly up-regulated in samples with high percentage of surrounding and infiltrating lymphocytes. In summary, this study highlights the robustness of using a targeted panel to understand the composition and regulatory mechanism of the TME and tumor immune response. Citation Format: Yuan-Chieh Ku, Warren Tom, Yongming Sun, Alex Pankov, Tim Looney, Fiona Hyland, Janice Au-Young, Ann Mongan. A targeted NGS solution to evaluate gene expression signature of the tumor microenvironment from 40 NSCLC FFPE and matched fresh frozen samples [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5364. doi:10.1158/1538-7445.AM2017-5364

Abstract 5363: Measuring gene expression at the tumor microenvironment: A comparison between nCounter PanCancer Immune Profiling Panel and Oncomine Immune Response Research Assay

Abstract The tumor microenvironment, especially infiltrating T lymphocytes and inflammatory molecules, is believed to be highly relevant to the tumor’s sensitivity to cancer checkpoint blockade therapy. At the same time, the exact markers that are predictive of response for each therapeutic agent are still the subject of active investigations. To address the need for better understanding of the effect of different T cell subsets, antigen presentation, and tumor killing, gene expression profiling presents an attractive means to simultaneously evaluate the tumor microenvironment and cancer cells. In this study we compare the results and performance of the nCounter PanCancer Immune Profiling Panel and the Oncomine Immune Response Research Assay, both of which are designed to measure the expression of genes indicative of an immune response and potential immune-editing activities by tumor cells. The nCounter panel detects gene expression by counting unique probes that hybridize target mRNA, while the Oncomine panel employs NGS to sequence and count reads derived from the targets. While both panels are designed to work with FFPE samples, The nCounter panel expects 100 ng of unamplified RNA, while the Oncomine panel requires only 10 ng total RNA with its AmpliSeq technology. The two panels share 254 common genes, which constitute the basis for this comparison. Across 12 cancer samples (breast, H&N, melanoma, NSCLC, and RCC), results show that the Oncomine panel offers 20% higher dynamic range, thereby providing more robust readouts about the differences among samples. More importantly, with virtually no background noise, the underlying NGS technology provides an absolute zero readout for non-expressing genes which significantly improves the sensitivity for detecting low expressing genes, as their presence can be confirmed by as few as two reads. The two technologies show moderate correlation (R ~ 0.7), with the Oncomine panel more strongly correlating with qPCR (R ~ 0.9). Finally, when clustered using all genes on the panel, only the Oncomine panel provides clear stratification of cancer types, thus allowing the panel to be used for tissue type confirmation in addition to evaluating the immune response. Citation Format: Ann Mongan, Warren Tom, Janice Au-Young, Aleksandr Pankov, Gauri Ganpule, Fiona Hyland. Measuring gene expression at the tumor microenvironment: A comparison between nCounter PanCancer Immune Profiling Panel and Oncomine Immune Response Research Assay [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5363. doi:10.1158/1538-7445.AM2017-5363

Abstract 5406: Low-input transcript profiling with enhanced sensitivity using a highly efficient, low-bias and strand-specific RNA-Seq library preparation method

Abstract RNA-seq has become the most popular method for transcriptome analysis and is widely used to study gene expression, detect mutations, fusion transcripts, alternative splicing, and post-transcriptional modifications. It is becoming the method of choice to detect alterations in diseases, especially cancer, to provide insights on the various molecular pathways perturbed by changes in the transcriptome and study their implications. As RNA-seq is being adopted in molecular diagnostics and biomarker identifications, the need for good quality, reproducible library preparation methods using very low amounts of RNA input, or precious clinical samples, is increasing. To address these challenges, we have developed a strand-specific RNA-seq library preparation method that retains information about which strand of DNA is transcribed, from as low as 5 ng total RNA input. Strand specificity is important for correct annotation of genes, identification of antisense transcripts with potential regulatory roles, and accurate determination of gene expression levels in the presence of antisense transcripts. Enhanced sensitivity to detect transcripts with even coverage across their length offers a non-biased approach for accurate quantification of transcript levels. Methods: Enriched poly-A mRNA or ribo-depleted RNA (Universal Human Reference RNA) was used to make libraries with our strand specific library preparation method. Library quality and quantity were analyzed on an Agilent Bioanalyzer, pooled at equimolar ratio and sequenced on Illumina’s Nextseq 500. Paired end reads were mapped to a human reference genome (hg19) using Hisat2 and sequencing metrics were calculated using Picard's RNA-seq Metrics and RSeqC tools. Transcript abundance was measured using Salmon and the Ensembl GRCh38 CDS sequences. Results: Libraries prepared with our streamlined method using inputs that range from 5ng to 1ug show greater than 98% directionality at all input levels. GC content analysis, gene body coverage and gene expression correlation are similar for all inputs tested (5ng to 1ug), even though input amounts vary by over three orders of magnitude. These consistent results are recapitulated with the spiked-in ERCC controls at all inputs. Conclusions: Our library preparation method is streamlined and can be used for a wide range of input RNA without any major modifications to the protocol, making it an easy to follow, convenient method for RNA-seq library preparation. In addition, our method has increased sensitivity and specificity, especially for low-abundance transcripts, reduced PCR duplicates and sequence bias, delivering high quality strand-specific data even for low input RNA. Finally, our method is compatible with both poly A-tail enriched and ribosomal RNA depleted samples, and is amenable to large-scale library construction and automation. Citation Format: Keerthana Krishnan, Erbay Yigit, Mehmet Karaca, Deyra Rodriguez, Bradley Langhorst, Timur Shtatland, Daniela Munafo, Pingfang Liu, Lynne Apone, Vaishnavi Panchapakesa, Karen Duggan, Christine Sumner, Christine Rozzi, Fiona Stewart, Laurie Mazzola, Joanna Bybee, Danielle Rivizzigno, Eileen T. Dimalanta, Theodore B. Davis. Low-input transcript profiling with enhanced sensitivity using a highly efficient, low-bias and strand-specific RNA-Seq library preparation method [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5406. doi:10.1158/1538-7445.AM2017-5406

Optimized Method for Robust Transcriptome Profiling of Minute Tissues Using Laser Capture Microdissection and Low-Input RNA-Seq.

Obtaining high quality RNA from complex biological tissues, such as the brain, is needed for establishing high-fidelity cell-type specific transcriptomes. Although combining genetic labeling techniques with laser capture microdissection (LCM) is generally sufficient, concerns over RNA degradation and limited yields call into question results of many sequencing studies. Here we set out to address both of these issues by: (1) developing a fluorescence-assisted LCM protocol that yields high quality RNA from fresh-frozen tissues; and (2) determining a suitable RNA-Seq library generation method for limited amounts of RNA (1–5 ng total RNA). The latter focused on comparing commercially available kits able to produce libraries of sufficient concentration and complexity while limiting PCR amplification biases. We find that high quality RNA (RNA integrity number, RIN, >9) of sufficient concentration can be isolated from laser-captured material from thinly-sectioned tissues when digestion time and temperature are minimized. Furthermore, we found that library generation approaches that retain ribosomal RNA (rRNA) through cDNA library generation required fewer cycles of PCR, minimizing bias in the resulting libraries. Lastly, end stage depletion of rRNA prior to sequencing enriches for target RNAs, thereby increasing read depth and level of gene detection while decreasing sequencing costs. Here we describe our protocol for generating robust RNA-Seq libraries from laser-captured tissue and demonstrate that with this method, we obtain samples with RNA quality superior to the current standard in the LCM field, and show that low-input RNA-Seq kits that minimize PCR bias produce high fidelity sequencing metrics with less variability compared to current practices.

High-Throughput Mapping of 2'-O-Me Residues in RNA Using Next-Generation Sequencing (Illumina RiboMethSeq Protocol).

Detection of RNA modifications in native RNAs is a tedious and laborious task, since the global level of these residues is low and most of the suitable physico-chemical methods require purification of the RNA of interest almost to homogeneity. To overcome these limitations, methods based on RT-driven primer extension have been developed and successfully used, sometimes in combination with a specific chemical treatment. Nowadays, some of these approaches have been coupled to high-throughput sequencing technologies, allowing the access to transcriptome-wide data. RNA 2'-O-methylation is one of the ubiquitous nucleotide modifications found in many RNA types from bacteria, archaea, and eukarya. Here, we describe a reliable and optimized protocol based on alkaline fragmentation of total RNA coupled to a commonly used ligation approach followed by Illumina sequencing. We describe the methodology for detection and relative quantification of 2'-O-methylations with a high sensitivity and reproducibility even with a limited amount of starting material (1 ng of total RNA). Altogether this technique unlocks a technological barrier since it will be applicable for routine parallel treatment of biological and clinical samples to decipher the functions of 2'-O-methylations in pathologies.

Gene Expression Profiling of Reticulocytes in Patients with Hereditary Spherocytosis after Depleting Globin Gene Transcripts

IntroductionGlobal transcriptome analysis of circulating reticulocytes using microarrays is challenging due to the high abundance of globin transcripts. In order to accurately measure the transcriptome in reticulocytes, we planned to study the reticulocytes transcriptome profile before and after globin depletion. Patients with Hereditary Spherocytosis (HS) were recruited because of high reticulocytosis and also there was no global trancriptome profile available for the disease to the best of our knowledge.MethodsReticulocytes were purified by passing EDTA anticoagulated red cell suspensions through a column of microcrystalline cellulose and α-cellulose mixture. The extraction of total RNA was done fresh prior to microarray analysis using TRIzol reagent (Invitrogen) and RNeasy columns as per manufacturer’s instruction. The globin transcripts were depleted with starting quantity of 3 µg of total RNA.using GLOBINclearTMHuman kit (Ambion, Austin, TX). RNA quality were assessed before and after depletion using Agilent Technologies 2100 Bioanalyzer in each sample used for the study.The Agilent Low Input Quick Amp Labelling kit was used to generate cRNA with a sample input of 200 ng total RNA in single-color microarray analysis (SurePrint G3 Human Gene Expression v3 8x60K Microarray, G4858A- 072363). cRNA was hybridized for 17 hours at 65˚C. After thorough washing, the results were extracted with Feature Extraction Project software and analysed using GeneSpring v13.0. For this study, differentially expressed genes were identified according to the criteria: t-test unpaired p (Corr) cut-off = 0.05 and fold-change cut-off of 2.0. These differentially expressed genes were imported in PANTHER (Protein Analysis Through Evolutionary Relationships) classification system to classify the genes.ResultsTranscriptome profiling of reticulocyte RNA from HS patients revealed 5318 annotated transcripts and 2744 lnc-RNA/uncharacterized transcripts to be differentially regulated before globin depletion. After globin depletion increase (8196 annotated transcripts and 4292 lnc-RNA) in the number of differentially regulated genes were observed. An increase of 54% and 56% was seen in annotated transcripts and noncoding transcripts respectively after globin depletion. This increased coverage may be attributed to the removal of the globin transcript. Gene Ontology analysis for molecular function, biological process, cellular component and protein class did not reveal any difference in the percentage of the category, though the number of transcripts was more after depletion. There were approximately 7% more pathways in the globin transcript depleted samples. Namely, 2-arachidonoylglycerol biosynthesis, biotin biosynthesis, carnitine metabolism, cobalamin biosynthesis, coenzyme a linked carnitine metabolism, cysteine biosynthesis, flavin biosynthesis, phenylalanine biosynthesis, sulfate assimilation, tyrosine biosynthesis pathways, etc. were elucidated only after globin depletion.ConclusionsThe globin transcript reduction followed by microarray analysis represents a robust methodology for studying pathophysiology of hematologic diseases which are not related to globin chain disorders. Furthermore we describe the global transcriptome profile of HS for the first time. Globin transcript depletion elucidates the better number of the transcripts and eventually the pathways which may have been earlier masked under the abundance of globin mRNAs. Pathway analysis and validation experiments are required to explain the role of these differentially expressed genes in the pathophysiology of HS. DisclosuresNo relevant conflicts of interest to declare.