Search

Archival tissues represent a rich resource for clinical genomic studies, particularly when coupled with comprehensive medical records. Use of these in next generation sequencing (NGS) is a priority. Nine formalin-fixed paraffin-embedded (FFPE) DNA extraction methods were evaluated using twelve FFPE samples of varying tissue types. Quality assessment included total yield, percent dsDNA, fragment analysis and multiplex PCR. After assessment, three tissue types from four FFPE DNA methods were selected for NGS downstream evaluation, targeted and whole exome sequencing. In addition, two low input library protocols were evaluated for WES. Analysis revealed average coverage across the target regions for WES was ~20-30X for all four FFPE DNA extraction methods. For the targeted panels, the highest molecular tag coverage was obtained with the Kingfisher FFPE extraction method. The genotype concordance was 99% for the commonly called variant positions between all four extraction methods with the targeted PCR NGS panel and 96% with WES. Assessing quality of extracted DNA aids in selecting the optimal NGS approach, and the choice of both DNA extraction and library preparation approaches can impact the performance of archival tissue in NGS. [ABSTRACT FROM AUTHOR]

Epstein-Barr virus (EBV) nuclear antigen 1 (EBNA1) is the only viral protein expressed in all EBV-positive tumors as it is essential for the maintenance, replication and transcription of the virus genome. According to the polymorphism of residue 487 in EBNA1 gene, EBV isolates can be classified into five subtypes: P-ala, P-thr, V-val, V-leu and V-pro. Whether these EBNA1 subtypes contribute to different tissue tropism of EBV and are consequently associated with certain malignancies remain to be determined. To elucidate the relationship, one hundred and ten EBV-positive lymphoma tissues of different types from Northern China, a non-NPC endemic area, were tested for the five subtypes by nested-PCR and DNA sequencing. In addition, EBV type 1 and type 2 classification was typed by using standard PCR assays across type-specific regions of the EBNA3C genes. Four EBNA1 subtypes were identified: V-val (68.2%, 75/110), P-thrV (15.5%, 17/110), V-leuV (3.6%, 4/110) and P-ala (10.9%, 12/110). The distribution of the EBNA1 subtypes in the four lymphoma groups was not significantly different (p = 0.075), neither was that of the EBV type 1/type 2 (p = 0.089). Compared with the previous data of gastric carcinoma (GC), nasopharyngeal carcinoma (NPC) and throat washing (TW) from healthy donors, the distribution of EBNA1 subtypes in lymphoma differed significantly (p = 0.016), with a little higher frequency of P-ala subtype. The EBV type distribution between lymphoma and the other three groups was significantly different (p = 0.000, p = 0.000, p = 0.001, respectively). The proportion of type 1 and type 2 mixed infections was higher in lymphoma than that in GC, NPC and TW. In lymphomas, the distribution of EBNA1 subtypes in the three EBV types was not significantly different (p = 0.546). These data suggested that the variation patterns of EBNA1 gene may be geographic-associated rather than tumor-specific and the role of EBNA1 gene variations in tumorigenesis needs more extensive and deep explorations. [ABSTRACT FROM AUTHOR]

Folate receptor alpha (FOLR1) has been identified as a potential prognostic and therapeutic target in a number of cancers. A correlation has been shown between intense overexpression of FOLR1 in breast tumors and poor prognosis, yet there is limited examination of the distribution of FOLR1 across clinically relevant breast cancer subtypes. To explore this further, we used RNA-seq data from multiple patient cohorts to analyze the distribution of FOLR1 mRNA across breast cancer subtypes comprised of estrogen receptor positive (ER+), human epidermal growth factor receptor positive (HER2+), and triple negative (TNBC) tumors. FOLR1 expression varied within breast tumor subtypes; triple negative/basal tumors were significantly associated with increased expression of FOLR1 mRNA, compared to ER+ and HER2+ tumors. However, subsets of high level FOLR1 expressing tumors were observed in all clinical subtypes. These observations were supported by immunohistochemical analysis of tissue microarrays, with the largest number of 3+ positive tumors and highest H-scores of any subtype represented by triple negatives, and lowest by ER+ tumors. FOLR1 expression did not correlate to common clinicopathological parameters such as tumor stage and nodal status. To delineate the importance of FOLR1 overexpression in triple negative cancers, RNA-interference was used to deplete FOLR1 in overexpressing triple negative cell breast lines. Loss of FOLR1 resulted in growth inhibition, whereas FOLR1 overexpression promoted folate uptake and growth advantage in low folate conditions. Taken together, our data suggests patients with triple negative cancers expressing high FOLR1 expression represent an important population of patients that may benefit from targeted anti-FOLR1 therapy. This may prove particularly helpful for a large number of patients who would typically be classified as triple negative and who to this point have been left without any targeted treatment options. [ABSTRACT FROM AUTHOR]

Objectives: The sequencing by the PolyA selection is the most common approach for library preparation. With limited amount or degraded RNA, alternative protocols such as the NuGEN have been developed. However, it is not yet clear how the different library preparations affect the downstream analyses of the broad applications of RNA sequencing. Methods and Materials: Eight human mammary epithelial cell (HMEC) lines with high quality RNA were sequenced by Illumina’s mRNA-Seq PolyA selection and NuGEN ENCORE library preparation. The following analyses and comparisons were conducted: 1) the numbers of genes captured by each protocol; 2) the impact of protocols on differentially expressed gene detection between biological replicates; 3) expressed single nucleotide variant (SNV) detection; 4) non-coding RNAs, particularly lincRNA detection; and 5) intragenic gene expression. Results: Sequences from the NuGEN protocol had lower (75%) alignment rate than the PolyA (over 90%). The NuGEN protocol detected fewer genes (12–20% less) with a significant portion of reads mapped to non-coding regions. A large number of genes were differentially detected between the two protocols. About 17–20% of the differentially expressed genes between biological replicates were commonly detected between the two protocols. Significantly higher numbers of SNVs (5–6 times) were detected in the NuGEN samples, which were largely from intragenic and intergenic regions. The NuGEN captured fewer exons (25% less) and had higher base level coverage variance. While 6.3% of reads were mapped to intragenic regions in the PolyA samples, the percentages were much higher (20–25%) for the NuGEN samples. The NuGEN protocol did not detect more known non-coding RNAs such as lincRNAs, but targeted small and “novel” lincRNAs. Conclusion: Different library preparations can have significant impacts on downstream analysis and interpretation of RNA-seq data. The NuGEN provides an alternative for limited or degraded RNA but it has limitations for some RNA-seq applications. [ABSTRACT FROM AUTHOR]

Advantages of RNA-Seq over array based platforms are quantitative gene expression and discovery of expressed single nucleotide variants (eSNVs) and fusion transcripts from a single platform, but the sensitivity for each of these characteristics is unknown. We measured gene expression in a set of manually degraded RNAs, nine pairs of matched fresh-frozen, and FFPE RNA isolated from breast tumor with the hybridization based, NanoString nCounter (226 gene panel) and with whole transcriptome RNA-Seq using RiboZeroGold ScriptSeq V2 library preparation kits. We performed correlation analyses of gene expression between samples and across platforms. We then specifically assessed whole transcriptome expression of lincRNA and discovery of eSNVs and fusion transcripts in the FFPE RNA-Seq data. For gene expression in the manually degraded samples, we observed Pearson correlations of >0.94 and >0.80 with NanoString and ScriptSeq protocols, respectively. Gene expression data for matched fresh-frozen and FFPE samples yielded mean Pearson correlations of 0.874 and 0.783 for NanoString (226 genes) and ScriptSeq whole transcriptome protocols respectively, p<2x10-16. Specifically for lincRNAs, we observed superb Pearson correlation (0.988) between matched fresh-frozen and FFPE pairs. FFPE samples across NanoString and RNA-Seq platforms gave a mean Pearson correlation of 0.838. In FFPE libraries, we detected 53.4% of high confidence SNVs and 24% of high confidence fusion transcripts. Sensitivity of fusion transcript detection was not overcome by an increase in depth of sequencing up to 3-fold (increase from ~56 to ~159 million reads). Both NanoString and ScriptSeq RNA-Seq technologies yield reliable gene expression data for degraded and FFPE material. The high degree of correlation between NanoString and RNA-Seq platforms suggests discovery based whole transcriptome studies from FFPE material will produce reliable expression data. The RiboZeroGold ScriptSeq protocol performed particularly well for lincRNA expression from FFPE libraries, but detection of eSNV and fusion transcripts was less sensitive. [ABSTRACT FROM AUTHOR]

Our goal in these analyses was to use genomic features from a test set of primary breast tumors to build an integrated transcriptome landscape model that makes relevant hypothetical predictions about the biological and/or clinical behavior of HER2-positive breast cancer. We interrogated RNA-Seq data from benign breast lesions, ER+, triple negative, and HER2-positive tumors to identify 685 differentially expressed genes, 102 alternatively spliced genes, and 303 genes that expressed single nucleotide sequence variants (eSNVs) that were associated with the HER2-positive tumors in our survey panel. These features were integrated into a transcriptome landscape model that identified 12 highly interconnected genomic modules, each of which represents a cellular processes pathway that appears to define the genomic architecture of the HER2-positive tumors in our test set. The generality of the model was confirmed by the observation that several key pathways were enriched in HER2-positive TCGA breast tumors. The ability of this model to make relevant predictions about the biology of breast cancer cells was established by the observation that integrin signaling was linked to lapatinib sensitivity in vitro and strongly associated with risk of relapse in the NCCTG N9831 adjuvant trastuzumab clinical trial dataset. Additional modules from the HER2 transcriptome model, including ubiquitin-mediated proteolysis, TGF-beta signaling, RHO-family GTPase signaling, and M-phase progression, were linked to response to lapatinib and paclitaxel in vitro and/or risk of relapse in the N9831 dataset. These data indicate that an integrated transcriptome landscape model derived from a test set of HER2-positive breast tumors has potential for predicting outcome and for identifying novel potential therapeutic strategies for this breast cancer subtype. [ABSTRACT FROM AUTHOR]