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.

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.

Systematic identification of cancer-related long noncoding RNAs and aberrant alternative splicing of quintuple-negative lung adenocarcinoma through RNA-Seq.

OBJECTIVES: Lung adenocarcinoma (LUAD) is a common subtype of non-small cell lung cancer prevalent in Asia. There is a dearth of understanding regarding the transcriptome landscape of LUAD without primary known driver mutations. In this study, LUAD samples without well-known driver mutations occurring in EGFR, KRAS, ALK, ROS1 or RET (quintuple-negative) were used for transcriptome study with a focus on long noncoding RNAs (lncRNAs), alternative splicing and gene fusions. MATERIALS AND METHODS: 24 pairs of LUAD and adjacent normal samples and 13 tumor-only samples derived from 37 quintuple-negative patients were used. Differentially expressed lncRNA transcripts were detected by paired t-test and were validated by qPCR. Functions of lncRNAs were predicted by co-expressed mRNAs. Aberrant splicing events in LUAD were identified using MISO. In addition, gene fusions were screened by SOAPfuse. RESULTS AND CONCLUSION: In total, 90 and 153 up- or down-regulated lncRNA transcripts were detected in LUAD samples in comparison with the adjacent normal samples. The most significantly differentially expressed lncRNA transcript was ENST00000598996.1 (FENDRR) down-regulated in LUAD. By lncRNA-mRNA co-expression analysis, functions of 14 lncRNAs were predicted. The predicted functions included vasculature development, immune response, cell cycle and respiratory gaseous exchange. Furthermore, six co-expressed pairs of lncRNAs and their nearby protein coding genes were identified as associated with lung development. This study also identified two highly recurrent (22 in 24) differential exon skipping events occurring in MYH14 and ESYT2 with exon including isoforms of both genes up-regulated in isoform percentage in LUAD samples. On the other hand, two out of 24 LUAD samples possessed the driver mutation exon 14 skipping of MET. The transcriptional alterations of LUAD samples without well-known driver mutations identified in the study can be used as references for future research. The translational values of these transcriptional changes are also worthy of further investigation.

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.

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.

An effort to use human-based exome capture methods to analyze chimpanzee and macaque exomes.

Non-human primates have emerged as an important resource for the study of human disease and evolution. The characterization of genomic variation between and within non-human primate species could advance the development of genetically defined non-human primate disease models. However, non-human primate specific reagents that would expedite such research, such as exon-capture tools, are lacking. We evaluated the efficiency of using a human exome capture design for the selective enrichment of exonic regions of non-human primates. We compared the exon sequence recovery in nine chimpanzees, two crab-eating macaques and eight Japanese macaques. Over 91% of the target regions were captured in the non-human primate samples, although the specificity of the capture decreased as evolutionary divergence from humans increased. Both intra-specific and inter-specific DNA variants were identified; Sanger-based resequencing validated 85.4% of 41 randomly selected SNPs. Among the short indels identified, a majority (54.6%-77.3%) of the variants resulted in a change of 3 base pairs, consistent with expectations for a selection against frame shift mutations. Taken together, these findings indicate that use of a human design exon-capture array can provide efficient enrichment of non-human primate gene regions. Accordingly, use of the human exon-capture methods provides an attractive, cost-effective approach for the comparative analysis of non-human primate genomes, including gene-based DNA variant discovery.

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.

The function of cortactin in the clustering of acetylcholine receptors at the vertebrate neuromuscular junction.

BACKGROUND: Postsynaptic enrichment of acetylcholine receptors (AChRs) at the vertebrate neuromuscular junction (NMJ) depends on the activation of the muscle receptor tyrosine MuSK by neural agrin. Agrin-stimulation of MuSK is known to initiate an intracellular signaling cascade that leads to the clustering of AChRs in an actin polymerization-dependent manner, but the molecular steps which link MuSK activation to AChR aggregation remain incompletely defined. METHODOLOGY/PRINCIPAL FINDINGS: In this study we used biochemical, cell biological and molecular assays to investigate a possible role in AChR clustering of cortactin, a protein which is a tyrosine kinase substrate and a regulator of F-actin assembly and which has also been previously localized at AChR clustering sites. We report that cortactin was co-enriched at AChR clusters in situ with its target the Arp2/3 complex, which is a key stimulator of actin polymerization in cells. Cortactin was further preferentially tyrosine phosphorylated at AChR clustering sites and treatment of myotubes with agrin significantly enhanced the tyrosine phosphorylation of cortactin. Importantly, forced expression in myotubes of a tyrosine phosphorylation-defective cortactin mutant (but not wild-type cortactin) suppressed agrin-dependent AChR clustering, as did the reduction of endogenous cortactin levels using RNA interference, and introduction of the mutant cortactin into muscle cells potently inhibited synaptic AChR aggregation in response to innervation. CONCLUSION: Our results suggest a novel function of phosphorylation-dependent cortactin signaling downstream from agrin/MuSK in facilitating AChR clustering at the developing NMJ.