Extremely high genetic diversity in a single tumor points to prevalence of non-Darwinian cell evolution.

The prevailing view that the evolution of cells in a tumor is driven by Darwinian selection has never been rigorously tested. Because selection greatly affects the level of intratumor genetic diversity, it is important to assess whether intratumor evolution follows the Darwinian or the non-Darwinian mode of evolution. To provide the statistical power, many regions in a single tumor need to be sampled and analyzed much more extensively than has been attempted in previous intratumor studies. Here, from a hepatocellular carcinoma (HCC) tumor, we evaluated multiregional samples from the tumor, using either whole-exome sequencing (WES) (n = 23 samples) or genotyping (n = 286) under both the infinite-site and infinite-allele models of population genetics. In addition to the many single-nucleotide variations (SNVs) present in all samples, there were 35 "polymorphic" SNVs among samples. High genetic diversity was evident as the 23 WES samples defined 20 unique cell clones. With all 286 samples genotyped, clonal diversity agreed well with the non-Darwinian model with no evidence of positive Darwinian selection. Under the non-Darwinian model, MALL (the number of coding region mutations in the entire tumor) was estimated to be greater than 100 million in this tumor. DNA sequences reveal local diversities in small patches of cells and validate the estimation. In contrast, the genetic diversity under a Darwinian model would generally be orders of magnitude smaller. Because the level of genetic diversity will have implications on therapeutic resistance, non-Darwinian evolution should be heeded in cancer treatments even for microscopic tumors.

Rapid growth of a hepatocellular carcinoma and the driving mutations revealed by cell-population genetic analysis of whole-genome data.

We present the analysis of the evolution of tumors in a case of hepatocellular carcinoma. This case is particularly informative about cancer growth dynamics and the underlying driving mutations. We sampled nine different sections from three tumors and seven more sections from the adjacent nontumor tissues. Selected sections were subjected to exon as well as whole-genome sequencing. Putative somatic mutations were then individually validated across all 9 tumor and 7 nontumor sections. Among the mutations validated, 24 were amino acid changes; in addition, 22 large indels/copy number variants (>1 Mb) were detected. These somatic mutations define four evolutionary lineages among tumor cells. Separate evolution and expansion of these lineages were recent and rapid, each apparently having only one lineage-specific protein-coding mutation. Hence, by using a cell-population genetic definition, this approach identified three coding changes (CCNG1, P62, and an indel/fusion gene) as tumor driver mutations. These three mutations, affecting cell cycle control and apoptosis, are functionally distinct from mutations that accumulated earlier, many of which are involved in inflammation/immunity or cell anchoring. These distinct functions of mutations at different stages may reflect the genetic interactions underlying tumor growth.

Quartet DNA reference materials and datasets for comprehensively evaluating germline variant calling performance.

BACKGROUND: Genomic DNA reference materials are widely recognized as essential for ensuring data quality in omics research. However, relying solely on reference datasets to evaluate the accuracy of variant calling results is incomplete, as they are limited to benchmark regions. Therefore, it is important to develop DNA reference materials that enable the assessment of variant detection performance across the entire genome. RESULTS: We established a DNA reference material suite from four immortalized cell lines derived from a family of parents and monozygotic twins. Comprehensive reference datasets of 4.2 million small variants and 15,000 structural variants were integrated and certified for evaluating the reliability of germline variant calls inside the benchmark regions. Importantly, the genetic built-in-truth of the Quartet family design enables estimation of the precision of variant calls outside the benchmark regions. Using the Quartet reference materials along with study samples, batch effects are objectively monitored and alleviated by training a machine learning model with the Quartet reference datasets to remove potential artifact calls. Moreover, the matched RNA and protein reference materials and datasets from the Quartet project enables cross-omics validation of variant calls from multiomics data. CONCLUSIONS: The Quartet DNA reference materials and reference datasets provide a unique resource for objectively assessing the quality of germline variant calls throughout the whole-genome regions and improving the reliability of large-scale genomic profiling.

Integrated genomic analysis identifies deregulated JAK/STAT-MYC-biosynthesis axis in aggressive NK-cell leukemia.

Aggressive NK-cell leukemia (ANKL) is a rare form of NK cell neoplasm that is more prevalent among people from Asia and Central and South America. Patients usually die within days to months, even after receiving prompt therapeutic management. Here we performed the first comprehensive study of ANKL by integrating whole genome, transcriptome and targeted sequencing, cytokine array as well as functional assays. Mutations in the JAK-STAT pathway were identified in 48% (14/29) of ANKL patients, while the extracellular STAT3 stimulator IL10 was elevated by an average of 56-fold (P < 0.0001) in the plasma of all patients examined. Additional frequently mutated genes included TP53 (34%), TET2 (28%), CREBBP (21%) and MLL2 (21%). Patient NK leukemia cells showed prominent activation of STAT3 phosphorylation, MYC expression and transcriptional activities in multiple metabolic pathways. Functionally, STAT3 activation and MYC expression were critical for the proliferation and survival of ANKL cells. STAT signaling regulated the MYC transcription program, and both STAT signaling and MYC transcription were required to maintain the activation of nucleotide synthesis and glycolysis. Collectively, the JAK-STAT pathway represents a major target for genomic alterations and IL10 stimulation in ANKL. This newly discovered JAK/STAT-MYC-biosynthesis axis may provide opportunities for the development of novel therapeutic strategies in treating this subtype of leukemia.

Long non-coding RNA C5orf66-AS1 is downregulated in pituitary null cell adenomas and is associated with their invasiveness.

Pituitary null cell adenoma is a challenging clinical condition, and its pathogenesis remains to be elucidated. We performed this study to determine the roles of C5orf66-AS1, NORAD, and TINCR in the pathogenesis and invasion of pituitary null cell adenomas. Expression of the three long non-coding RNAs in pituitary null cell adenoma tissues of 11 patients and normal pituitary tissues from four donors was examined by performing quantitative reverse transcription-polymerase chain reaction. We found that C5orf66-AS1 expression was lower in pituitary null cell adenoma tissues than in normal pituitary tissues. Moreover, C5orf66-AS1 expression level was significantly lower in invasive pituitary null cell adenomas than in non-invasive ones. After transfection of C5orf66-AS1 into pituitary adenoma cells, assessment of cell viability and invasion suggested that overexpressed C5orf66-AS1 inhibited cell viability and cell invasion. In silico algorithms predicted several cis- and trans-acting target genes of C5orf66-AS1, including PITX1 and SCGB3A1. In addition, expression of some of the predicted target genes was determined using microarray data of another cohort with pituitary null cell adenomas. It showed that some of these target genes were differentially expressed between pituitary null cell adenoma tissues and normal pituitary tissues as well as between invasive and non-invasive tumors. Co-expression analysis in RNA sequencing data showed that PAQR7 was the most correlated gene of C5orf66-AS1 and that several predicted trans-acting target genes, including SCGB3A1, were highly correlated with C5orf66-AS1. NORAD and TINCR expression was not statistically significant in the complete cohort; however, a negative correlation was observed between NORAD expression and maximum tumor diameter in some subgroups. These results indicate that C5orf66-AS1 suppresses the development and invasion of pituitary null cell adenomas. However, our results do not provide enough statistical evidence to support the roles of NORAD and TINCR in the development and invasion of pituitary null cell adenomas.

Identification of functional cooperative mutations of SETD2 in human acute leukemia.

Acute leukemia characterized by chromosomal rearrangements requires additional molecular disruptions to develop into full-blown malignancy, yet the cooperative mechanisms remain elusive. Using whole-genome sequencing of a pair of monozygotic twins discordant for MLL (also called KMT2A) gene-rearranged leukemia, we identified a transforming MLL-NRIP3 fusion gene and biallelic mutations in SETD2 (encoding a histone H3K36 methyltransferase). Moreover, loss-of-function point mutations in SETD2 were recurrent (6.2%) in 241 patients with acute leukemia and were associated with multiple major chromosomal aberrations. We observed a global loss of H3K36 trimethylation (H3K36me3) in leukemic blasts with mutations in SETD2. In the presence of a genetic lesion, downregulation of SETD2 contributed to both initiation and progression during leukemia development by promoting the self-renewal potential of leukemia stem cells. Therefore, our study provides compelling evidence for SETD2 as a new tumor suppressor. Disruption of the SETD2-H3K36me3 pathway is a distinct epigenetic mechanism for leukemia development.