High-Throughput Peptide Arrays Identify Potential Diagnostic Autoantibody Signatures in Early-Stage Lung Adenocarcinoma.

BACKGROUND: Early diagnosis is critical to lung adenocarcinoma patients' survival but faces inadequacies in convenient early detection. METHODS: We applied a comprehensive microarray of 130,000 peptides to detect "autoantibody signature" that is autoantibodies binding to mimotopes for early detection of stage 0-I LUAD. Plasma samples were collected from 147 early-stage lung adenocarcinoma (Early-LUAD), 108 benign lung disease (BLD), and 122 normal healthy controls (NHC). Clinical characteristics, low-dose CT (LDCT), and laboratory tests were incorporated into correlation analysis. RESULTS: We identified 143 and 133 autoantibody signatures, distinguishing Early-LUAD from NHC/BLD in the discovery cohort. Autoantibody signatures significantly correlated with age, stage, tumor size, basophil count, and IgM level (P < 0.05). The random forest models based on differential autoantibody signatures displayed AUC of 0.92 and 0.87 to discern Early-LUAD from NHC/BLD in the validation cohort, respectively. Compared with LDCT, combining autoantibody signature and LDCT improved the positive predictive value from 50% to 78.33% (P = 0.049). In addition, autoantibody signatures displayed higher sensitivity of 72.4% to 81.0% compared with the combinational tumor markers (cyfra21.1, NSE, SCC, ProGRP) with a sensitivity of 22.4% (P = 0.000). Proteins matched by differential peptides were enriched in cancer-related PI3K/Akt, MAPK, and Wnt pathways. Overlaps between matched epitopes and autoantibody signatures illustrated the underlying engagement of autoantibodies in immune recognition. CONCLUSIONS: Collectively, autoantibody signatures identified by a high-throughput peptide microarray have the potential to detect Early-LUAD, which could assist LDCT to better diagnose Early-LUAD. IMPACT: Novel sensitive autoantibody signatures can adjuvant LDCT to better diagnose LUAD at very early stage.

XY-Meta: A High-Efficiency Search Engine for Large-Scale Metabolome Annotation with Accurate FDR Estimation.

FDR control has been a huge challenge for large-scale metabolome annotation. Although recent research indicated that the target-decoy strategy could be implemented to estimate FDR, it is hard to perform FDR control due to the difficulty of getting a reliable decoy database because of the complex fragmentation mechanism of metabolites and ubiquitous isomers. To tackle this problem, we developed a decoy generation method, which generates forged spectra from the reference target database by preserving the original reference signals to simulate the presence of isomers of metabolites. Benchmarks on GNPS data sets in Passatutto showed that the decoy database generated by our method is closer to the actual FDR than other methods, especially in the low FDR range (0-0.05). Large-scale metabolite annotation on 35 data sets showed that strict FDR reduced the number of annotated metabolites but increased the spectral efficiency, indicating the necessity of quality control. We recommended that the FDR threshold should be set to 0.01 in large-scale metabolite annotation. We implemented decoy generation, database search, and FDR control into a search engine called XY-Meta. It facilitates large-scale metabolome annotation applications.

Whole-genome and Transcriptome Sequencing of Prostate Cancer Identify New Genetic Alterations Driving Disease Progression.

BACKGROUND: Global disparities in prostate cancer (PCa) incidence highlight the urgent need to identify genomic abnormalities in prostate tumors in different ethnic populations including Asian men. OBJECTIVE: To systematically explore the genomic complexity and define disease-driven genetic alterations in PCa. DESIGN, SETTING, AND PARTICIPANTS: The study sequenced whole-genome and transcriptome of tumor-benign paired tissues from 65 treatment-naive Chinese PCa patients. Subsequent targeted deep sequencing of 293 PCa-relevant genes was performed in another cohort of 145 prostate tumors. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: The genomic alteration landscape in PCa was analyzed using an integrated computational pipeline. Relationships with PCa progression and survival were analyzed using nonparametric test, log-rank, and multivariable Cox regression analyses. RESULTS AND LIMITATIONS: We demonstrated an association of high frequency of CHD1 deletion with a low rate of TMPRSS2-ERG fusion and relatively high percentage of mutations in androgen receptor upstream activator genes in Chinese patients. We identified five putative clustered deleted tumor suppressor genes and provided experimental and clinical evidence that PCDH9, deleted/loss in approximately 23% of tumors, functions as a novel tumor suppressor gene with prognostic potential in PCa. Furthermore, axon guidance pathway genes were frequently deregulated, including gain/amplification of PLXNA1 gene in approximately 17% of tumors. Functional and clinical data analyses showed that increased expression of PLXNA1 promoted prostate tumor growth and independently predicted prostate tumor biochemical recurrence, metastasis, and poor survival in multi-institutional cohorts of patients with PCa. A limitation of this study is that other genetic alterations were not experimentally investigated. CONCLUSIONS: There are shared and salient genetic characteristics of PCa in Chinese and Caucasian men. Novel genetic alterations in PCDH9 and PLXNA1 were associated with disease progression. PATIENT SUMMARY: We reported the first large-scale and comprehensive genomic data of prostate cancer from Asian population. Identification of these genetic alterations may help advance prostate cancer diagnosis, prognosis, and treatment.

Comprehensive Characterization of Oncogenic Drivers in Asian Lung Adenocarcinoma.

INTRODUCTION: The incidence rate of lung adenocarcinoma (LUAD), the predominant histological subtype of lung cancer, is elevated in Asians, particularly in female nonsmokers. The mutation patterns in LUAD in Asians might be distinct from those in LUAD in whites. METHODS: We profiled 271 resected LUAD tumors (mainly stage I) to characterize the genomic landscape of LUAD in Asians with a focus on female nonsmokers. RESULTS: Mutations in EGFR, KRAS, erb-b2 receptor tyrosine kinase 2 gene (ERBB2), and BRAF; gene fusions involving anaplastic lymphoma receptor tyrosine kinase gene (ALK), ROS1, and ret proto-oncogene (RET); and Met Proto-Oncogene Tyrosine Kinase (MET) exon 14 skipping were the major drivers in LUAD in Asians, exhibiting mutually exclusive and differing prevalence from those reported in studies of LUAD in non-Asians. In addition, we identified a novel mutational signature of XNX (the mutated base N in the middle flanked by two identical bases at the 5' and 3' positions) that was overrepresented in LUAD tumors in nonsmokers and negatively correlated with the overall mutational frequency. CONCLUSIONS: In this cohort, approximately 85% of individuals have known driver mutations (EGFR 59.4%, KRAS 7.4%, ALK 7.4%, ERBB2 2.6%, ROS1 2.2%, RET 2.2%, MET 1.8%, BRAF 1.1%, and NRAS 0.4%). Seventy percent of smokers and 90% of nonsmokers had defined oncogenic drivers matching the U.S. Food and Drug Administration-approved targeted therapies.

Genome-wide identification of RNA editing in hepatocellular carcinoma.

We did whole-transcriptome sequencing and whole-genome sequencing on nine pairs of Hepatocellular carcinoma (HCC) tumors and matched adjacent tissues to identify RNA editing events. We identified mean 26,982 editing sites with mean 89.5% canonical A→G edits in each sample using an improved bioinformatics pipeline. The editing rate was significantly higher in tumors than adjacent normal tissues. Comparing the difference between tumor and normal tissues of each patient, we found 7 non-synonymous tissue specific editing events including 4 tumor-specific edits and 3 normal-specific edits in the coding region, as well as 292 edits varying in editing degree. The significant expression changes of 150 genes associated with RNA editing were found in tumors, with 3 of the 4 most significant genes being cancer related. Our results show that editing might be related to higher gene expression. These findings indicate that RNA editing modification may play an important role in the development of HCC.

Discovery of biclonal origin and a novel oncogene SLC12A5 in colon cancer by single-cell sequencing.

Single-cell sequencing is a powerful tool for delineating clonal relationship and identifying key driver genes for personalized cancer management. Here we performed single-cell sequencing analysis of a case of colon cancer. Population genetics analyses identified two independent clones in tumor cell population. The major tumor clone harbored APC and TP53 mutations as early oncogenic events, whereas the minor clone contained preponderant CDC27 and PABPC1 mutations. The absence of APC and TP53 mutations in the minor clone supports that these two clones were derived from two cellular origins. Examination of somatic mutation allele frequency spectra of additional 21 whole-tissue exome-sequenced cases revealed the heterogeneity of clonal origins in colon cancer. Next, we identified a mutated gene SLC12A5 that showed a high frequency of mutation at the single-cell level but exhibited low prevalence at the population level. Functional characterization of mutant SLC12A5 revealed its potential oncogenic effect in colon cancer. Our study provides the first exome-wide evidence at single-cell level supporting that colon cancer could be of a biclonal origin, and suggests that low-prevalence mutations in a cohort may also play important protumorigenic roles at the individual level.

Whole-genome sequencing of matched primary and metastatic hepatocellular carcinomas.

BACKGROUND: To gain biological insights into lung metastases from hepatocellular carcinoma (HCC), we compared the whole-genome sequencing profiles of primary HCC and paired lung metastases. METHODS: We used whole-genome sequencing at 33X-43X coverage to profile somatic mutations in primary HCC (HBV+) and metachronous lung metastases (> 2 years interval). RESULTS: In total, 5,027-13,961 and 5,275-12,624 somatic single-nucleotide variants (SNVs) were detected in primary HCC and lung metastases, respectively. Generally, 38.88-78.49% of SNVs detected in metastases were present in primary tumors. We identified 65-221 structural variations (SVs) in primary tumors and 60-232 SVs in metastases. Comparison of these SVs shows very similar and largely overlapped mutated segments between primary and metastatic tumors. Copy number alterations between primary and metastatic pairs were also found to be closely related. Together, these preservations in genomic profiles from liver primary tumors to metachronous lung metastases indicate that the genomic features during tumorigenesis may be retained during metastasis. CONCLUSIONS: We found very similar genomic alterations between primary and metastatic tumors, with a few mutations found specifically in lung metastases, which may explain the clinical observation that both primary and metastatic tumors are usually sensitive or resistant to the same systemic treatments.

Whole-genome sequencing identifies recurrent mutations in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the most deadly cancers worldwide and has no effective treatment, yet the molecular basis of hepatocarcinogenesis remains largely unknown. Here we report findings from a whole-genome sequencing (WGS) study of 88 matched HCC tumor/normal pairs, 81 of which are Hepatitis B virus (HBV) positive, seeking to identify genetically altered genes and pathways implicated in HBV-associated HCC. We find beta-catenin to be the most frequently mutated oncogene (15.9%) and TP53 the most frequently mutated tumor suppressor (35.2%). The Wnt/beta-catenin and JAK/STAT pathways, altered in 62.5% and 45.5% of cases, respectively, are likely to act as two major oncogenic drivers in HCC. This study also identifies several prevalent and potentially actionable mutations, including activating mutations of Janus kinase 1 (JAK1), in 9.1% of patients and provides a path toward therapeutic intervention of the disease.

Genome-wide survey of recurrent HBV integration in hepatocellular carcinoma.

To survey hepatitis B virus (HBV) integration in liver cancer genomes, we conducted massively parallel sequencing of 81 HBV-positive and 7 HBV-negative hepatocellular carcinomas (HCCs) and adjacent normal tissues. We found that HBV integration is observed more frequently in the tumors (86.4%) than in adjacent liver tissues (30.7%). Copy-number variations (CNVs) were significantly increased at HBV breakpoint locations where chromosomal instability was likely induced. Approximately 40% of HBV breakpoints within the HBV genome were located within a 1,800-bp region where the viral enhancer, X gene and core gene are located. We also identified recurrent HBV integration events (in ≥ 4 HCCs) that were validated by RNA sequencing (RNA-seq) and Sanger sequencing at the known and putative cancer-related TERT, MLL4 and CCNE1 genes, which showed upregulated gene expression in tumor versus normal tissue. We also report evidence that suggests that the number of HBV integrations is associated with patient survival.

Single-cell exome sequencing and monoclonal evolution of a JAK2-negative myeloproliferative neoplasm.

Tumor heterogeneity presents a challenge for inferring clonal evolution and driver gene identification. Here, we describe a method for analyzing the cancer genome at a single-cell nucleotide level. To perform our analyses, we first devised and validated a high-throughput whole-genome single-cell sequencing method using two lymphoblastoid cell line single cells. We then carried out whole-exome single-cell sequencing of 90 cells from a JAK2-negative myeloproliferative neoplasm patient. The sequencing data from 58 cells passed our quality control criteria, and these data indicated that this neoplasm represented a monoclonal evolution. We further identified essential thrombocythemia (ET)-related candidate mutations such as SESN2 and NTRK1, which may be involved in neoplasm progression. This pilot study allowed the initial characterization of the disease-related genetic architecture at the single-cell nucleotide level. Further, we established a single-cell sequencing method that opens the way for detailed analyses of a variety of tumor types, including those with high genetic complex between patients.

Single-cell sequencing analysis characterizes common and cell-lineage-specific mutations in a muscle-invasive bladder cancer.

BACKGROUND: Cancers arise through an evolutionary process in which cell populations are subjected to selection; however, to date, the process of bladder cancer, which is one of the most common cancers in the world, remains unknown at a single-cell level. RESULTS: We carried out single-cell exome sequencing of 66 individual tumor cells from a muscle-invasive bladder transitional cell carcinoma (TCC). Analyses of the somatic mutant allele frequency spectrum and clonal structure revealed that the tumor cells were derived from a single ancestral cell, but that subsequent evolution occurred, leading to two distinct tumor cell subpopulations. By analyzing recurrently mutant genes in an additional cohort of 99 TCC tumors, we identified genes that might play roles in the maintenance of the ancestral clone and in the muscle-invasive capability of subclones of this bladder cancer, respectively. CONCLUSIONS: This work provides a new approach of investigating the genetic details of bladder tumoral changes at the single-cell level and a new method for assessing bladder cancer evolution at a cell-population level.

Frequent mutations of chromatin remodeling genes in transitional cell carcinoma of the bladder.

Transitional cell carcinoma (TCC) is the most common type of bladder cancer. Here we sequenced the exomes of nine individuals with TCC and screened all the somatically mutated genes in a prevalence set of 88 additional individuals with TCC with different tumor stages and grades. In our study, we discovered a variety of genes previously unknown to be mutated in TCC. Notably, we identified genetic aberrations of the chromatin remodeling genes (UTX, MLL-MLL3, CREBBP-EP300, NCOR1, ARID1A and CHD6) in 59% of our 97 subjects with TCC. Of these genes, we showed UTX to be altered substantially more frequently in tumors of low stages and grades, highlighting its potential role in the classification and diagnosis of bladder cancer. Our results provide an overview of the genetic basis of TCC and suggest that aberration of chromatin regulation might be a hallmark of bladder cancer.

Coexpression of delta- and mu-opioid receptors in nociceptive sensory neurons.

Morphine-induced analgesia and antinociceptive tolerance are known to be modulated by interaction between delta-opioid receptors (DORs) and mu-opioid receptors (MORs) in the pain pathway. However, evidence for expression of DORs in nociceptive small-diameter neurons in dorsal root ganglia (DRG) and for coexistence of DORs with MORs and neuropeptides has recently been challenged. We now report, using in situ hybridization, single-cell PCR, and immunostaining, that DORs are widely expressed not only in large DRG neurons but in small ones and coexist with MORs in peptidergic small DRG neurons, with protachykinin-dependent localization in large dense-core vesicles. Importantly, both DOR and MOR agonists reduce depolarization-induced Ca(2+) currents in single small DRG neurons and inhibit afferent C-fiber synaptic transmission in the dorsal spinal cord. Thus, coexistence of DORs and MORs in small DRG neurons is a basis for direct interaction of opioid receptors in modulation of nociceptive afferent transmission and opioid analgesia.