Protocol optimization of a targeted sequencing panel for genomic profiling of bronchoalveolar lavage fluid in lung cancer.

INTRODUCTION: We investigated a commercially available sequencing panel to study the effect of sequencing depth, variant calling strategy, and targeted sequencing region on identifying tumor-derived variants in cell-free bronchoalveolar lavage (cfBAL) DNA compared with plasma cfDNA. METHODS: Sequencing was performed at low or high coverage using two filtering algorithms to identify tumor variants on two panels targeting 77 and 197 genes respectively. RESULTS: One hundred and four sequencing files from 40 matched DNA samples of cfBAL, plasma, germline leukocytes, and archival tumor specimens in 10 patients with early-stage lung cancer were analyzed. By low-coverage sequencing, tumor-derived cfBAL variants were detected in 5/10 patients (50%) compared with 2/10 (20%) for plasma. High-coverage sequencing did not affect the number of tumor-derived variants detected in either biospecimen type. Accounting for germline mutations eliminated false-positive plasma calls regardless of coverage (0/10 patients with tumor-derived variants identified) and increased the number of cfBAL calls (5/10 patients with tumor-derived variants identified). These results were not affected by the number of targeted genes.

Verification of prognostic expression biomarkers is improved by examining enriched leukemic blasts rather than mononuclear cells from acute myeloid leukemia patients.

BACKGROUND: Studies have not systematically compared the ability to verify performance of prognostic transcripts in paired bulk mononuclear cells versus viable CD34-expressing leukemic blasts from patients with acute myeloid leukemia. We hypothesized that examining the homogenous leukemic blasts will yield different biological information and may improve prognostic performance of expression biomarkers. METHODS: To assess the impact of cellular heterogeneity on expression biomarkers in acute myeloid leukemia, we systematically examined paired mononuclear cells and viable CD34-expressing leukemic blasts from SWOG diagnostic specimens. After enrichment, patients were assigned into discovery and validation cohorts based on availability of extracted RNA. Analyses of RNA sequencing data examined how enrichment impacted differentially expressed genes associated with pre-analytic variables, patient characteristics, and clinical outcomes. RESULTS: Blast enrichment yielded significantly different expression profiles and biological pathways associated with clinical characteristics (e.g., cytogenetics). Although numerous differentially expressed genes were associated with clinical outcomes, most lost their prognostic significance in the mononuclear cells and blasts after adjusting for age and ELN risk, with only 11 genes remaining significant for overall survival in both cell populations (CEP70, COMMD7, DNMT3B, ECE1, LNX2, NEGR1, PIK3C2B, SEMA4D, SMAD2, TAF8, ZNF444). To examine the impact of enrichment on biomarker verification, these 11 candidate biomarkers were examined by quantitative RT/PCR in the validation cohort. After adjusting for ELN risk and age, expression of 4 genes (CEP70, DNMT3B, ECE1, and PIK3CB) remained significantly associated with overall survival in the blasts, while none met statistical significance in mononuclear cells. CONCLUSIONS: This study provides insights into biological information gained/lost by examining viable CD34-expressing leukemic blasts versus mononuclear cells from the same patient and shows an improved verification rate for expression biomarkers in blasts.

Somatic mutation but not aneuploidy differentiates lung cancer in never-smokers and smokers.

Lung cancer in never-smokers disproportionately affects older women. To understand the mutational landscape of this cohort, we performed detailed genome characterization of 73 lung adenocarcinomas from participants of the Women’s Health Initiative (WHI). We find enrichment of EGFR mutations in never-/light-smokers and KRAS mutations in heavy smokers as expected, but we also show that the specific variants of these genes differ by smoking status, with important therapeutic implications. Mutational signature analysis revealed signatures of clock, APOBEC, and DNA repair deficiency in never-/light-smokers; however, the mutational load of these signatures did not differ significantly from those found in smokers. Last, tumors from both smokers and never-/light-smokers shared copy number subtypes, with no significant differences in aneuploidy. Thus, the genomic landscape of lung cancer in never-/light-smokers and smokers is predominantly differentiated by somatic mutations and not copy number alterations.

Highly conserved intragenic HSV-2 sequences: Results from next-generation sequencing of HSV-2 UL and US regions from genital swabs collected from 3 continents.

INTRODUCTION: Understanding the variability in circulating herpes simplex virus type 2 (HSV-2) genomic sequences is critical to the development of HSV-2 vaccines. METHODS: Genital lesion swabs containing ≥ 107log10 copies HSV DNA collected from Africa, the USA, and South America underwent next-generation sequencing, followed by K-mer based filtering and de novo genomic assembly. Sites of heterogeneity within coding regions in unique long and unique short (UL_US) regions were identified. Phylogenetic trees were created using maximum likelihood reconstruction. RESULTS: Among 46 samples from 38 persons, 1468 intragenic base-pair substitutions were identified. The maximum nucleotide distance between strains for concatenated UL_US segments was 0.4%. Phylogeny did not reveal geographic clustering. The most variable proteins had non-synonymous mutations in < 3% of amino acids. CONCLUSIONS: Unenriched HSV-2 DNA can undergo next-generation sequencing to identify intragenic variability. The use of clinical swabs for sequencing expands the information that can be gathered directly from these specimens.

Worldwide circulation of HSV-2 × HSV-1 recombinant strains.

Homo sapiens harbor two distinct, medically significant species of simplexviruses, herpes simplex virus (HSV)-1 and HSV-2, with estimated divergence 6-8 million years ago (MYA). Unexpectedly, we found that circulating HSV-2 strains can contain HSV-1 DNA segments in three distinct genes. Using over 150 genital swabs from North and South America and Africa, we detected recombinants worldwide. Common, widely distributed gene UL39 genotypes are parsimoniously explained by an initial >457 basepair (bp) HSV-1 × HSV-2 crossover followed by back-recombination to HSV-2. Blocks of >244 and >539 bp of HSV-1 DNA within genes UL29 and UL30, respectively, have reached near fixation, with a minority of strains retaining sequences we posit as ancestral HSV-2. Our data add to previous in vitro and animal work, implying that in vivo cellular co-infection with HSV-1 and HSV-2 yields viable interspecies recombinants in the natural human host.

Tumor-Infiltrating Merkel Cell Polyomavirus-Specific T Cells Are Diverse and Associated with Improved Patient Survival.

Tumor-infiltrating CD8+ T cells are associated with improved survival of patients with Merkel cell carcinoma (MCC), an aggressive skin cancer causally linked to Merkel cell polyomavirus (MCPyV). However, CD8+ T-cell infiltration is robust in only 4% to 18% of MCC tumors. We characterized the T-cell receptor (TCR) repertoire restricted to one prominent epitope of MCPyV (KLLEIAPNC, "KLL") and assessed whether TCR diversity, tumor infiltration, or T-cell avidity correlated with clinical outcome. HLA-A*02:01/KLL tetramer+ CD8+ T cells from MCC patient peripheral blood mononuclear cells (PBMC) and tumor-infiltrating lymphocytes (TIL) were isolated via flow cytometry. TCRß (TRB) sequencing was performed on tetramer+ cells from PBMCs or TILs (n = 14) and matched tumors (n = 12). Functional avidity of T-cell clones was determined by IFNγ production. We identified KLL tetramer+ T cells in 14% of PBMC and 21% of TIL from MCC patients. TRB repertoires were strikingly diverse (397 unique TRBs were identified from 12 patients) and mostly private (only one TCRb clonotype shared between two patients). An increased fraction of KLL-specific TIL (>1.9%) was associated with significantly increased MCC-specific survival P = 0.0009). T-cell cloning from four patients identified 42 distinct KLL-specific TCRa/b pairs. T-cell clones from patients with improved MCC-specific outcomes were more avid (P < 0.05) and recognized an HLA-appropriate MCC cell line. T cells specific for a single MCPyV epitope display marked TCR diversity within and between patients. Intratumoral infiltration by MCPyV-specific T cells was associated with significantly improved MCC-specific survival, suggesting that augmenting the number or avidity of virus-specific T cells may have therapeutic benefit. Cancer Immunol Res; 5(2); 137-47. ©2017 AACR.

Expression and functional characterization of CD33 transcript variants in human acute myeloid leukemia.

With the demonstration of improved survival of some acute myeloid leukemia (AML) patients with the CD33 antibody-drug conjugate, gemtuzumab ozogamicin (GO), CD33 has been validated as a target for antigen-specific immunotherapy. Since previous studies identified a CD33 splice variant missing exon 2 (CD33∆E2) and, consequently, the immune-dominant membrane-distal V-set domain, we investigated the expression and functional characteristics of CD33 transcript variants in AML. In primary AML specimens, we not only found full-length CD33 (CD33FL) and CD33∆E2 but also corresponding variants containing an alternate exon 7 predicted to encode a CD33 protein lacking most of the intracellular domain (CD33E7a and, not previously described, CD33∆E2,E7a) in almost all cases. In acute leukemia cell sublines engineered to express individual CD33 splice variants, all splice variants had endocytic properties. CD33FL and CD33E7a mediated similar degrees of GO cytotoxicity, whereas CD33∆E2 and CD33∆E2,E7a could not serve as target for GO. Co-expression of CD33∆E2 did not interfere with CD33FL endocytosis and did not impact CD33FL-mediated GO cytotoxicity. Together, our findings document a greater-than-previously thought complexity of CD33 expression in human AML. They identify CD33 variants that lack exon 2 and are not recognized by current CD33-directed therapeutics as potential target for future unconjugated or conjugated antibodies.

Complete Unique Genome Sequence, Expression Profile, and Salivary Gland Tissue Tropism of the Herpesvirus 7 Homolog in Pigtailed Macaques.

UNLABELLED: Human herpesvirus 6A (HHV-6A), HHV-6B, and HHV-7 are classified as roseoloviruses and are highly prevalent in the human population. Roseolovirus reactivation in an immunocompromised host can cause severe pathologies. While the pathogenic potential of HHV-7 is unclear, it can reactivate HHV-6 from latency and thus contributes to severe pathological conditions associated with HHV-6. Because of the ubiquitous nature of roseoloviruses, their roles in such interactions and the resulting pathological consequences have been difficult to study. Furthermore, the lack of a relevant animal model for HHV-7 infection has hindered a better understanding of its contribution to roseolovirus-associated diseases. Using next-generation sequencing analysis, we characterized the unique genome of an uncultured novel pigtailed macaque roseolovirus. Detailed genomic analysis revealed the presence of gene homologs to all 84 known HHV-7 open reading frames. Phylogenetic analysis confirmed that the virus is a macaque homolog of HHV-7, which we have provisionally named Macaca nemestrina herpesvirus 7 (MneHV7). Using high-throughput RNA sequencing, we observed that the salivary gland tissue samples from nine different macaques had distinct MneHV7 gene expression patterns and that the overall number of viral transcripts correlated with viral loads in parotid gland tissue and saliva. Immunohistochemistry staining confirmed that, like HHV-7, MneHV7 exhibits a natural tropism for salivary gland ductal cells. We also observed staining for MneHV7 in peripheral nerve ganglia present in salivary gland tissues, suggesting that HHV-7 may also have a tropism for the peripheral nervous system. Our data demonstrate that MneHV7-infected macaques represent a relevant animal model that may help clarify the causality between roseolovirus reactivation and diseases. IMPORTANCE: Human herpesvirus 6A (HHV-6A), HHV-6B, and HHV-7 are classified as roseoloviruses. We have recently discovered that pigtailed macaques are naturally infected with viral homologs of HHV-6 and HHV-7, which we provisionally named MneHV6 and MneHV7, respectively. In this study, we confirm that MneHV7 is genetically and biologically similar to its human counterpart, HHV-7. We determined the complete unique MneHV7 genome sequence and provide a comprehensive annotation of all genes. We also characterized viral transcription profiles in salivary glands from naturally infected macaques. We show that broad transcriptional activity across most of the viral genome is associated with high viral loads in infected parotid glands and that late viral protein expression is detected in salivary duct cells and peripheral nerve ganglia. Our study provides new insights into the natural behavior of an extremely prevalent virus and establishes a basis for subsequent investigations of the mechanisms that cause HHV-7 reactivation and associated disease.

Identification of long-lived proteins retained in cells undergoing repeated asymmetric divisions.

Long-lived proteins have been implicated in age-associated decline in metazoa, but they have only been identified in extracellular matrices or postmitotic cells. However, the aging process also occurs in dividing cells undergoing repeated asymmetric divisions. It was not clear whether long-lived proteins exist in asymmetrically dividing cells or whether they are involved in aging. Here we identify long-lived proteins in dividing cells during aging using the budding yeast, Saccharomyces cerevisiae. Yeast mother cells undergo a limited number of asymmetric divisions that define replicative lifespan. We used stable-isotope pulse-chase and total proteome mass-spectrometry to identify proteins that were both long-lived and retained in aging mother cells after ∼ 18 cells divisions. We identified ∼ 135 proteins that we designate as long-lived asymmetrically retained proteins (LARPS). Surprisingly, the majority of LARPs appeared to be stable fragments of their original full-length protein. However, 15% of LARPs were full-length proteins and we confirmed several candidates to be long-lived and retained in mother cells by time-lapse microscopy. Some LARPs localized to the plasma membrane and remained robustly in the mother cell upon cell division. Other full-length LARPs were assembled into large cytoplasmic structures that had a strong bias to remain in mother cells. We identified age-associated changes to LARPs that include an increase in their levels during aging because of their continued synthesis, which is not balanced by turnover. Additionally, several LARPs were posttranslationally modified during aging. We suggest that LARPs contribute to age-associated phenotypes and likely exist in other organisms.

Phenotypic and transcriptional fidelity of patient-derived colon cancer xenografts in immune-deficient mice.

Xenografts of human colorectal cancer (CRC) in immune-deficient mice have great potential for accelerating the study of tumor biology and therapy. We evaluated xenografts established in NOD/scid/IL2Rγ-null mice from the primary or metastatic tumors of 27 patients with CRC to estimate their capacity for expanding tumor cells for in vitro studies and to assess how faithfully they recapitulated the transcriptional profile of their parental tumors. RNA-seq analysis of parental human CRC tumors and their derivative xenografts demonstrated that reproducible transcriptional changes characterize the human tumor to murine xenograft transition. In most but not all cases, the human stroma, vasculature, and hematopoietic elements were systematically replaced by murine analogues while the carcinoma component persisted. Once established as xenografts, human CRC cells that could be propagated by serial transplantation remained transcriptionally stable. Three histologically atypical xenografts, established from patients with peritoneal metastases, contained abundant human stromal elements and blood vessels in addition to human tumor cells. The transcriptomes of these mixed tumor/stromal xenografts did not closely resemble those of their parental tumors, and attempts to propagate such xenografts by serial transplantation were unsuccessful. Stable expression of numerous genes previously identified as high priority targets for immunotherapy was observed in most xenograft lineages. Aberrant expression in CRC cells of human genes that are normally only expressed in hematopoietic cells was also observed. Our results suggest that human CRC cells expanded in murine xenografts have great utility for studies of tumor immunobiology and targeted therapies such as immunotherapy but also identify potential limitations.

Significance of expression of ITGA5 and its splice variants in acute myeloid leukemia: a report from the Children's Oncology Group.

Acute myeloid leukemia (AML) encompasses a heterogeneous group of diseases, and novel biomarkers for risk refinement and stratification are needed to optimize patient care. To identify novel risk factors, we performed transcriptome sequencing on 68 diagnostic AML samples and identified 2 transcript variants (-E2 and -E2/3) of the α-subunit (ITGA5) of the very late antigen-5 integrin. We then quantified expression of ITGA5 and these splice variants in specimens from participants of the AAML03P1 trial. We found no association between ITGA5 expression and clinical outcome. In contrast, patients with the highest relative expression (Q4) of the -E2/3 ITGA5 splice variant less likely had low-risk disease than Q1-3 patients (21% vs. 38%, P = 0.027). Q4 patients had worse response to chemotherapy with a higher proportion having persistent minimal residual disease (50% vs. 23%, P = 0.003) and inferior overall survival (at 5 years: 48% vs. 67%, P = 0.015); the latter association was limited to low-risk patients (Q4 vs. Q1-3: 56% vs. 85%, P = 0.043) and was not seen in standard-risk (51% vs. 60%, P = 0.340) or high-risk (33% vs. 38%, P = 0.952) patients. Our exploratory studies indicate that transcriptome sequencing is useful for biomarker discovery, as exemplified by the identification of ITGA5 -E2/3 splice variant as potential novel adverse prognostic marker for low-risk AML that, if confirmed, could serve to further risk-stratify this patient subset.

Integrative proteomic analysis of serum and peritoneal fluids helps identify proteins that are up-regulated in serum of women with ovarian cancer.

BACKGROUND: We used intensive modern proteomics approaches to identify predictive proteins in ovary cancer. We identify up-regulated proteins in both serum and peritoneal fluid. To evaluate the overall performance of the approach we track the behavior of 20 validated markers across these experiments. METHODOLOGY: Mass spectrometry based quantitative proteomics following extensive protein fractionation was used to compare serum of women with serous ovarian cancer to healthy women and women with benign ovarian tumors. Quantitation was achieved by isotopically labeling cysteine amino acids. Label-free mass spectrometry was used to compare peritoneal fluid taken from women with serous ovarian cancer and those with benign tumors. All data were integrated and annotated based on whether the proteins have been previously validated using antibody-based assays. FINDINGS: We selected 54 quantified serum proteins and 358 peritoneal fluid proteins whose case-control differences exceeded a predefined threshold. Seventeen proteins were quantified in both materials and 14 are extracellular. Of 19 validated markers that were identified all were found in cancer peritoneal fluid and a subset of 7 were quantified in serum, with one of these proteins, IGFBP1, newly validated here. CONCLUSION: Proteome profiling applied to symptomatic ovarian cancer cases identifies a large number of up-regulated serum proteins, many of which are or have been confirmed by immunoassays. The number of currently known validated markers is highest in peritoneal fluid, but they make up a higher percentage of the proteins observed in both serum and peritoneal fluid, suggesting that the 10 additional markers in this group may be high quality candidates.

Brain-specific proteins decline in the cerebrospinal fluid of humans with Huntington disease.

We integrated five sets of proteomics data profiling the constituents of cerebrospinal fluid (CSF) derived from Huntington disease (HD)-affected and -unaffected individuals with genomics data profiling various human and mouse tissues, including the human HD brain. Based on an integrated analysis, we found that brain-specific proteins are 1.8 times more likely to be observed in CSF than in plasma, that brain-specific proteins tend to decrease in HD CSF compared with unaffected CSF, and that 81% of brain-specific proteins have quantitative changes concordant with transcriptional changes identified in different regions of HD brain. The proteins found to increase in HD CSF tend to be liver-associated. These protein changes are consistent with neurodegeneration, microgliosis, and astrocytosis known to occur in HD. We also discuss concordance between laboratories and find that ratios of individual proteins can vary greatly, but the overall trends with respect to brain or liver specificity were consistent. Concordance is highest between the two laboratories observing the largest numbers of proteins.

A statistical method for chromatographic alignment of LC-MS data.

Integrated liquid-chromatography mass-spectrometry (LC-MS) is becoming a widely used approach for quantifying the protein composition of complex samples. The output of the LC-MS system measures the intensity of a peptide with a specific mass-charge ratio and retention time. In the last few years, this technology has been used to compare complex biological samples across multiple conditions. One challenge for comparative proteomic profiling with LC-MS is to match corresponding peptide features from different experiments. In this paper, we propose a new method--Peptide Element Alignment (PETAL) that uses raw spectrum data and detected peak to simultaneously align features from multiple LC-MS experiments. PETAL creates spectrum elements, each of which represents the mass spectrum of a single peptide in a single scan. Peptides detected in different LC-MS data are aligned if they can be represented by the same elements. By considering each peptide separately, PETAL enjoys greater flexibility than time warping methods. While most existing methods process multiple data sets by sequentially aligning each data set to an arbitrarily chosen template data set, PETAL treats all experiments symmetrically and can analyze all experiments simultaneously. We illustrate the performance of PETAL on example data sets.