Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions.

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5 × 10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (p value <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

Spatially Resolved Single-Cell Omics: Methods, Challenges, and Future Perspectives.

Overlaying omics data onto spatial biological dimensions has been a promising technology to provide high-resolution insights into the interactome and cellular heterogeneity relative to the organization of the molecular microenvironment of tissue samples in normal and disease states. Spatial omics can be categorized into three major modalities: (a) next-generation sequencing-based assays, (b) imaging-based spatially resolved transcriptomics RNA approaches including in situ hybridization/in situ sequencing, and (c) imaging-based proteomics. These modalities allow assessment of transcripts and proteins at a cellular level, generating large and computationally challenging datasets. The lack of standardized computational pipelines to analyze and integrate these nonuniform structured data has made it necessary to apply artificial intelligence and machine learning strategies to best visualize and translate their complexity. In this review, we summarize the currently available techniques and computational strategies, highlight their advantages and limitations, and discuss their future prospects in the scientific field.

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia.

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

Multi-lineage heart-chip models drug cardiotoxicity and enhances maturation of human stem cell-derived cardiovascular cells.

Cardiovascular toxicity causes adverse drug reactions and may lead to drug removal from the pharmaceutical market. Cancer therapies can induce life-threatening cardiovascular side effects such as arrhythmias, muscle cell death, or vascular dysfunction. New technologies have enabled cardiotoxic compounds to be identified earlier in drug development. Human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes (CMs) and vascular endothelial cells (ECs) can screen for drug-induced alterations in cardiovascular cell function and survival. However, most existing hiPSC models for cardiovascular drug toxicity utilize two-dimensional, immature cells grown in static culture. Improved in vitro models to mechanistically interrogate cardiotoxicity would utilize more adult-like, mature hiPSC-derived cells in an integrated system whereby toxic drugs and protective agents can flow between hiPSC-ECs that represent systemic vasculature and hiPSC-CMs that represent heart muscle (myocardium). Such models would be useful for testing the multi-lineage cardiotoxicities of chemotherapeutic drugs such as VEGFR2/PDGFR-inhibiting tyrosine kinase inhibitors (VPTKIs). Here, we develop a multi-lineage, fully-integrated, cardiovascular organ-chip that can enhance hiPSC-EC and hiPSC-CM functional and genetic maturity, model endothelial barrier permeability, and demonstrate long-term functional stability. This microfluidic organ-chip harbors hiPSC-CMs and hiPSC-ECs on separate channels that can be subjected to active fluid flow and rhythmic biomechanical stretch. We demonstrate the utility of this cardiovascular organ-chip as a predictive platform for evaluating multi-lineage VPTKI toxicity. This study may lead to the development of new modalities for the evaluation and prevention of cancer therapy-induced cardiotoxicity.

Cis- and trans-eQTL TWAS of breast and ovarian cancer identify more than 100 risk associated genes in the BCAC and OCAC consortia.

Transcriptome-wide association studies (TWAS) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have only considered regulatory effects of risk associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWAS of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents a first look into the role of trans-eQTLs in the complex molecular mechanisms underlying these diseases.

A machine learning one-class logistic regression model to predict stemness for single cell transcriptomics and spatial omics.

Cell annotation is a crucial methodological component to interpreting single cell and spatial omics data. These approaches were developed for single cell analysis but are often biased, manually curated and yet unproven in spatial omics. Here we apply a stemness model for assessing oncogenic states to single cell and spatial omic cancer datasets. This one-class logistic regression machine learning algorithm is used to extract transcriptomic features from non-transformed stem cells to identify dedifferentiated cell states in tumors. We found this method identifies single cell states in metastatic tumor cell populations without the requirement of cell annotation. This machine learning model identified stem-like cell populations not identified in single cell or spatial transcriptomic analysis using existing methods. For the first time, we demonstrate the application of a ML tool across five emerging spatial transcriptomic and proteomic technologies to identify oncogenic stem-like cell types in the tumor microenvironment.

Incidence and prognostic impact of HER2-positivity loss after dual HER2-directed neoadjuvant therapy for HER2+ breast cancer.

BACKGROUND: Loss of HER2 "positivity" can occur in patients with residual disease after neoadjuvant treatment, but the incidence of HER2-positivity loss after neoadjuvant dual HER2-targeted treatment plus chemotherapy, the current standard-of-care for most early stage HER2-positive breast cancers, is not well described. Previous studies that report the HER2 discordance rate after neoadjuvant treatment also do not include the novel HER2-low category. In this retrospective study, we determine the incidence and prognostic impact of HER2-positivity loss, including the evolution to HER2-low disease, after neoadjuvant dual HER2-targeted therapy with chemotherapy. METHODS: Clinicopathologic data for patients with stage I-III HER2+ breast cancer diagnosed between 2015 and 2019 were reviewed in this single institution retrospective study. Patients who received dual HER2-targeted treatment with chemotherapy were included, and HER2 status before and after neoadjuvant therapy was interrogated. RESULTS: A total of 163 female patients were included in the analysis with a median age of 50 years. A pathologic complete response (pCR as defined by ypT0/is) was achieved in 102 (62.5%) of 163 evaluable patients. Among the 61 patients with residual disease after neoadjuvant therapy, 36 (59.0%) had HER2-positive and 25 (41.0%) had HER2-negative residual disease. Of the 25 patients with HER2-negative residual disease, 22 (88%) of patients were classified as HER2-low. After a median follow-up of 3.3 years, patients who retained HER2-positivity after neoadjuvant treatment had a 3-year IDFS rate of 91% (95% CI, 91%-100%), while patients who lost HER2-positivity had a 3-year IDFS rate of 82% (95% CI, 67%-100%). CONCLUSION: Almost half of patients with residual disease following neoadjuvant dual HER2-targeted therapy plus chemotherapy lost HER2-positivity. The loss of HER2-positivity may not confer negative prognostic impact, although the results were limited by short follow-up time. Further research on the HER2 status after neoadjuvant treatment may help guide treatment decisions in the adjuvant setting.

Copy Number Variants Are Ovarian Cancer Risk Alleles at Known and Novel Risk Loci.

BACKGROUND: Known risk alleles for epithelial ovarian cancer (EOC) account for approximately 40% of the heritability for EOC. Copy number variants (CNVs) have not been investigated as EOC risk alleles in a large population cohort. METHODS: Single nucleotide polymorphism array data from 13 071 EOC cases and 17 306 controls of White European ancestry were used to identify CNVs associated with EOC risk using a rare admixture maximum likelihood test for gene burden and a by-probe ratio test. We performed enrichment analysis of CNVs at known EOC risk loci and functional biofeatures in ovarian cancer-related cell types. RESULTS: We identified statistically significant risk associations with CNVs at known EOC risk genes; BRCA1 (PEOC = 1.60E-21; OREOC = 8.24), RAD51C (Phigh-grade serous ovarian cancer [HGSOC] = 5.5E-4; odds ratio [OR]HGSOC = 5.74 del), and BRCA2 (PHGSOC = 7.0E-4; ORHGSOC = 3.31 deletion). Four suggestive associations (P < .001) were identified for rare CNVs. Risk-associated CNVs were enriched (P < .05) at known EOC risk loci identified by genome-wide association study. Noncoding CNVs were enriched in active promoters and insulators in EOC-related cell types. CONCLUSIONS: CNVs in BRCA1 have been previously reported in smaller studies, but their observed frequency in this large population-based cohort, along with the CNVs observed at BRCA2 and RAD51C gene loci in EOC cases, suggests that these CNVs are potentially pathogenic and may contribute to the spectrum of disease-causing mutations in these genes. CNVs are likely to occur in a wider set of susceptibility regions, with potential implications for clinical genetic testing and disease prevention.

Dual Stromal Targeting Sensitizes Pancreatic Adenocarcinoma for Anti-Programmed Cell Death Protein 1 Therapy.

BACKGROUND & AIMS: The stroma in pancreatic ductal adenocarcinoma (PDAC) contributes to its immunosuppressive nature and therapeutic resistance. Herein we sought to modify signaling and enhance immunotherapy efficacy by targeting multiple stromal components through both intracellular and extracellular mechanisms. METHODS: A murine liver metastasis syngeneic model of PDAC was treated with focal adhesion kinase inhibitor (FAKi), anti-programmed cell death protein 1 (PD-1) antibody, and stromal hyaluronan (HA) degradation by PEGylated recombinant human hyaluronidase (PEGPH20) to assess immune and stromal modulating effects of these agents and their combinations. RESULTS: The results showed that HA degradation by PEGPH20 and reduction in phosphorylated FAK expression by FAKi leads to improved survival in PDAC-bearing mice treated with anti-PD-1 antibody. HA degradation in combination with FAKi and anti-PD-1 antibody increases T-cell infiltration and alters T-cell phenotype toward effector memory T cells. FAKi alters the expression of T-cell modulating cytokines and leads to changes in T-cell metabolism and increases in effector T-cell signatures. HA degradation in combination with anti-PD-1 antibody and FAKi treatments reduces granulocytes, including granulocytic- myeloid-derived suppressor cells and decreases C-X-C chemokine receptor type 4 (CXCR4)-expressing myeloid cells, particularly the CXCR4-expressing granulocytes. Anti-CXCR4 antibody combined with FAKi and anti-PD-1 antibody significantly decreases metastatic rates in the PDAC liver metastasis model. CONCLUSIONS: This represents the first preclinical study to identify synergistic effects of targeting both intracellular and extracellular components within the PDAC stroma and supports testing anti-CXCR4 antibody in combination with FAKi as a PDAC treatment strategy.

US Severe Acute Respiratory Syndrome Coronavirus 2 Epsilon Variant: Highly Transmissible but With an Adjusted Muted Host T-Cell Response.

BACKGROUND: The multiple mutations comprising the epsilon variant demonstrate the independent convergent evolution of severe acute respiratory syndrome coronavirus (SARS-CoV-2), with its spike protein mutation L452R present in the delta (L452R), kappa (L452R), and lambda (L452Q) variants. METHODS: Coronavirus disease 2019 (COVID-19) variants were detected in 1017 patients using whole-genome sequencing and were assessed for outcome and severity. The mechanistic effects of the epsilon versus non-epsilon variants were investigated using a multiomic approach including cellular response assays and paired cell and host transcriptomic and proteomic profiling. RESULTS: We found that patients carrying the epsilon variant had increased mortality risk but not increased hospitalizations (P < .02). Cells infected with live epsilon compared with non-epsilon virus displayed increased sensitivity to neutralization antibodies in all patients but a slightly protective response in vaccinated individuals (P < .001). That the epsilon SARS-CoV-2 variant is more infectious but less virulent is supported mechanistically in the down-regulation of viral processing pathways seen by multiomic analyses. Importantly, this paired transcriptomics and proteomic profiling of host cellular response to live virus revealed an altered leukocyte response and metabolic messenger RNA processing with the epsilon variant. To ascertain host response to SARS-CoV-2 infection, primary COVID-19-positive nasopharyngeal samples were transcriptomically profiled and revealed a differential innate immune response (P < .001) and an adjusted T-cell response in patients carrying the epsilon variant (P < .002). In fact, patients infected with SARS-CoV-2 and those vaccinated with the BNT162b2 vaccine have comparable CD4+/CD8+ T-cell immune responses to the epsilon variant (P < .05). CONCLUSIONS: While the epsilon variant is more infectious, by altering viral processing, we showed that patients with COVID-19 have adapted their innate immune response to this fitter variant. A protective T-cell response molecular signature is generated by this more transmissible variant in both vaccinated and unvaccinated patients.

Predicting master transcription factors from pan-cancer expression data.

Critical developmental "master transcription factors" (MTFs) can be subverted during tumorigenesis to control oncogenic transcriptional programs. Current approaches to identifying MTFs rely on ChIP-seq data, which is unavailable for many cancers. We developed the CaCTS (Cancer Core Transcription factor Specificity) algorithm to prioritize candidate MTFs using pan-cancer RNA sequencing data. CaCTS identified candidate MTFs across 34 tumor types and 140 subtypes including predictions for cancer types/subtypes for which MTFs are unknown, including e.g. PAX8, SOX17, and MECOM as candidates in ovarian cancer (OvCa). In OvCa cells, consistent with known MTF properties, these factors are required for viability, lie proximal to superenhancers, co-occupy regulatory elements globally, co-bind loci encoding OvCa biomarkers, and are sensitive to pharmacologic inhibition of transcription. Our predictions of MTFs, especially for tumor types with limited understanding of transcriptional drivers, pave the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Pleiotropy-guided transcriptome imputation from normal and tumor tissues identifies candidate susceptibility genes for breast and ovarian cancer.

Familial, sequencing, and genome-wide association studies (GWASs) and genetic correlation analyses have progressively unraveled the shared or pleiotropic germline genetics of breast and ovarian cancer. In this study, we aimed to leverage this shared germline genetics to improve the power of transcriptome-wide association studies (TWASs) to identify candidate breast cancer and ovarian cancer susceptibility genes. We built gene expression prediction models using the PrediXcan method in 681 breast and 295 ovarian tumors from The Cancer Genome Atlas and 211 breast and 99 ovarian normal tissue samples from the Genotype-Tissue Expression project and integrated these with GWAS meta-analysis data from the Breast Cancer Association Consortium (122,977 cases/105,974 controls) and the Ovarian Cancer Association Consortium (22,406 cases/40,941 controls). The integration was achieved through application of a pleiotropy-guided conditional/conjunction false discovery rate (FDR) approach in the setting of a TWASs. This identified 14 candidate breast cancer susceptibility genes spanning 11 genomic regions and 8 candidate ovarian cancer susceptibility genes spanning 5 genomic regions at conjunction FDR < 0.05 that were >1 Mb away from known breast and/or ovarian cancer susceptibility loci. We also identified 38 candidate breast cancer susceptibility genes and 17 candidate ovarian cancer susceptibility genes at conjunction FDR < 0.05 at known breast and/or ovarian susceptibility loci. The 22 genes identified by our cross-cancer analysis represent promising candidates that further elucidate the role of the transcriptome in mediating germline breast and ovarian cancer risk.

DNA methylation-based signatures classify sporadic pituitary tumors according to clinicopathological features.

BACKGROUND: Distinct genome-wide methylation patterns cluster pituitary neuroendocrine tumors (PitNETs) into molecular groups associated with specific clinicopathological features. Here we aim to identify, characterize, and validate methylation signatures that objectively classify PitNET into clinicopathological groups. METHODS: Combining in-house and publicly available data, we conducted an analysis of the methylome profile of a comprehensive cohort of 177 tumors (Panpit cohort) and 20 nontumor specimens from the pituitary gland. We also retrieved methylome data from an independent PitNET cohort (N = 86) to validate our findings. RESULTS: We identified three methylation clusters associated with adenohypophyseal cell lineages and functional status using an unsupervised approach. Differentially methylated probes (DMP) significantly distinguished the Panpit clusters and accurately assigned the samples of the validation cohort to their corresponding lineage and functional subtypes memberships. The DMPs were annotated in regulatory regions enriched with enhancer elements, associated with pathways and genes involved in pituitary cell identity, function, tumorigenesis, and invasiveness. Some DMPs correlated with genes with prognostic and therapeutic values in other intra- or extracranial tumors. CONCLUSIONS: We identified and validated methylation signatures, mainly annotated in enhancer regions that distinguished PitNETs by distinct adenohypophyseal cell lineages and functional status. These signatures provide the groundwork to develop an unbiased approach to classifying PitNETs according to the most recent classification recommended by the 2017 WHO and to explore their biological and clinical relevance in these tumors.

Ovarian Cancer Risk Variants Are Enriched in Histotype-Specific Enhancers and Disrupt Transcription Factor Binding Sites.

Quantifying the functional effects of complex disease risk variants can provide insights into mechanisms underlying disease biology. Genome-wide association studies have identified 39 regions associated with risk of epithelial ovarian cancer (EOC). The vast majority of these variants lie in the non-coding genome, where they likely function through interaction with gene regulatory elements. In this study we first estimated the heritability explained by known common low penetrance risk alleles for EOC. The narrow sense heritability (hg2) of EOC overall and high-grade serous ovarian cancer (HGSOCs) were estimated to be 5%-6%. Partitioned SNP heritability across broad functional categories indicated a significant contribution of regulatory elements to EOC heritability. We collated epigenomic profiling data for 77 cell and tissue types from Roadmap Epigenomics and ENCODE, and from H3K27Ac ChIP-seq data generated in 26 ovarian cancer and precursor-related cell and tissue types. We identified significant enrichment of risk single-nucleotide polymorphisms (SNPs) in active regulatory elements marked by H3K27Ac in HGSOCs. To further investigate how risk SNPs in active regulatory elements influence predisposition to ovarian cancer, we used motifbreakR to predict the disruption of transcription factor binding sites. We identified 469 candidate causal risk variants in H3K27Ac peaks that are predicted to significantly break transcription factor (TF) motifs. The most frequently broken motif was REST (p value = 0.0028), which has been reported as both a tumor suppressor and an oncogene. Overall, these systematic functional annotations with epigenomic data improve interpretation of EOC risk variants and shed light on likely cells of origin.

A Study of High-Grade Serous Ovarian Cancer Origins Implicates the SOX18 Transcription Factor in Tumor Development.

Fallopian tube secretory epithelial cells (FTSECs) are likely the main precursor cell type of high-grade serous ovarian cancers (HGSOCs), but these tumors may also arise from ovarian surface epithelial cells (OSECs). We profiled global landscapes of gene expression and active chromatin to characterize molecular similarities between OSECs (n = 114), FTSECs (n = 74), and HGSOCs (n = 394). A one-class machine learning algorithm predicts that most HGSOCs derive from FTSECs, with particularly high FTSEC scores in mesenchymal-type HGSOCs (padj < 8 × 10-4). However, a subset of HGSOCs likely derive from OSECs, particularly HGSOCs of the proliferative type (padj < 2 × 10-4), suggesting a dualistic model for HGSOC origins. Super-enhancer (SE) landscapes were also more similar between FTSECs and HGSOCs than between OSECs and HGSOCs (p < 2.2 × 10-16). The SOX18 transcription factor (TF) coincided with a HGSOC-specific SE, and ectopic overexpression of SOX18 in FTSECs caused epithelial-to-mesenchymal transition, indicating that SOX18 plays a role in establishing the mesenchymal signature of fallopian-derived HGSOCs.

Super-Enhancer-Associated LncRNA UCA1 Interacts Directly with AMOT to Activate YAP Target Genes in Epithelial Ovarian Cancer.

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of tumorigenesis, and yet their mechanistic roles remain challenging to characterize. Here, we integrate functional proteomics with lncRNA-interactome profiling to characterize Urothelial Cancer Associated 1 (UCA1), a candidate driver of ovarian cancer development. Reverse phase protein array (RPPA) analysis indicates that UCA1 activates transcription coactivator YAP and its target genes. In vivo RNA antisense purification (iRAP) of UCA1 interacting proteins identified angiomotin (AMOT), a known YAP regulator, as a direct binding partner. Loss-of-function experiments show that AMOT mediates YAP activation by UCA1, as UCA1 enhances the AMOT-YAP interaction to promote YAP dephosphorylation and nuclear translocation. Together, we characterize UCA1 as a lncRNA regulator of Hippo-YAP signaling and highlight the UCA1-AMOT-YAP signaling axis in ovarian cancer development.

Optimizing exosomal RNA isolation for RNA-Seq analyses of archival sera specimens.

Exosomes are endosome-derived membrane vesicles that contain proteins, lipids, and nucleic acids. The exosomal transcriptome mediates intercellular communication, and represents an understudied reservoir of novel biomarkers for human diseases. Next-generation sequencing enables complex quantitative characterization of exosomal RNAs from diverse sources. However, detailed protocols describing exosome purification for preparation of exosomal RNA-sequence (RNA-Seq) libraries are lacking. Here we compared methods for isolation of exosomes and extraction of exosomal RNA from human cell-free serum, as well as strategies for attaining equal representation of samples within pooled RNA-Seq libraries. We compared commercial precipitation with ultracentrifugation for exosome purification and confirmed the presence of exosomes via both transmission electron microscopy and immunoblotting. Exosomal RNA extraction was compared using four different RNA purification methods. We determined the minimal starting volume of serum required for exosome preparation and showed that high quality exosomal RNA can be isolated from sera stored for over a decade. Finally, RNA-Seq libraries were successfully prepared with exosomal RNAs extracted from human cell-free serum, cataloguing both coding and non-coding exosomal transcripts. This method provides researchers with strategic options to prepare RNA-Seq libraries and compare RNA-Seq data quantitatively from minimal volumes of fresh and archival human cell-free serum for disease biomarker discovery.