Phylovar: toward scalable phylogeny-aware inference of single-nucleotide variations from single-cell DNA sequencing data.

MOTIVATION: Single-nucleotide variants (SNVs) are the most common variations in the human genome. Recently developed methods for SNV detection from single-cell DNA sequencing data, such as SCIΦ and scVILP, leverage the evolutionary history of the cells to overcome the technical errors associated with single-cell sequencing protocols. Despite being accurate, these methods are not scalable to the extensive genomic breadth of single-cell whole-genome (scWGS) and whole-exome sequencing (scWES) data. RESULTS: Here, we report on a new scalable method, Phylovar, which extends the phylogeny-guided variant calling approach to sequencing datasets containing millions of loci. Through benchmarking on simulated datasets under different settings, we show that, Phylovar outperforms SCIΦ in terms of running time while being more accurate than Monovar (which is not phylogeny-aware) in terms of SNV detection. Furthermore, we applied Phylovar to two real biological datasets: an scWES triple-negative breast cancer data consisting of 32 cells and 3375 loci as well as an scWGS data of neuron cells from a normal human brain containing 16 cells and approximately 2.5 million loci. For the cancer data, Phylovar detected somatic SNVs with high or moderate functional impact that were also supported by bulk sequencing dataset and for the neuron dataset, Phylovar identified 5745 SNVs with non-synonymous effects some of which were associated with neurodegenerative diseases. AVAILABILITY AND IMPLEMENTATION: Phylovar is implemented in Python and is publicly available at https://github.com/NakhlehLab/Phylovar.

Single-cell mtDNA heteroplasmy in colorectal cancer.

Human mitochondria can be genetically distinct within the same individual, a phenomenon known as heteroplasmy. In cancer, this phenomenon seems exacerbated, and most mitochondrial mutations seem to be heteroplasmic. How this genetic variation is arranged within and among normal and tumor cells is not well understood. To address this question, here we sequenced single-cell mitochondrial genomes from multiple normal and tumoral locations in four colorectal cancer patients. Our results suggest that single cells, both normal and tumoral, can carry various mitochondrial haplotypes. Remarkably, this intra-cell heteroplasmy can arise before tumor development and be maintained afterward in specific tumoral cell subpopulations. At least in the colorectal patients studied here, the somatic mutations in the single-cells do not seem to have a prominent role in tumorigenesis.

Comparative analysis of capture methods for genomic profiling of circulating tumor cells in colorectal cancer.

The genomic profiling of circulating tumor cells (CTCs) in the bloodstream should provide clinically relevant information on therapeutic efficacy and help predict cancer survival. Here, we contrasted the genomic profiles of CTC pools recovered from metastatic colorectal cancer (mCRC) patients using different enrichment strategies (CellSearch, Parsortix, and FACS). Mutations inferred in the CTC pools differed depending on the enrichment strategy and, in all cases, represented a subset of the mutations detected in the matched primary tumor samples. However, the CTC pools from Parsortix, and in part, CellSearch, showed diversity estimates, mutational signatures, and drug-suitability scores remarkably close to those found in matching primary tumor samples. In addition, FACS CTC pools were enriched in apparent sequencing artifacts, leading to much higher genomic diversity estimates. Our results highlight the utility of CTCs to assess the genomic heterogeneity of individual tumors and help clinicians prioritize drugs in mCRC.

BRCA1 represses DNA replication initiation through antagonizing estrogen signaling and maintains genome stability in parallel with WEE1-MCM2 signaling during pregnancy.

The mammary gland undergoes fast cell proliferation during early pregnancy, yet the mechanism to ensure genome integrity during this highly proliferative stage is largely unknown. We show that pregnancy triggers replicative stresses leading to genetic instability in mice carrying a mammary specific disruption of breast cancer associated gene-1 (BRCA1). The fast cell proliferation was correlated with enhanced expression of most genes encoding replisomes, which are positively regulated by estrogen/ERα signaling but negatively regulated by BRCA1. Our further analysis revealed two parallel signaling pathways, which are mediated by ATR-CHK1 and WEE1-MCM2 and are responsible for regulating DNA replication checkpoint. Upon DNA damage, BRCA1 deficiency markedly enhances DNA replication initiation and preferably impairs DNA replication checkpoint mediated by ATR and CHK1. Meanwhile, DNA damage also activates WEE1-MCM2 signaling, which inhibits DNA replication initiation and enables BRCA1-deficient cells to avoid further genomic instability. Finally, we demonstrated that overriding this defense by WEE1 inhibition in combination with cisplatin, which causes DNA damage, serves as a promising therapeutic approach for killing BRCA1-deficient cancer cells.

Single-cell phylogenies reveal changes in the evolutionary rate within cancer and healthy tissues.

Cell lineages accumulate somatic mutations during organismal development, potentially leading to pathological states. The rate of somatic evolution within a cell population can vary due to multiple factors, including selection, a change in the mutation rate, or differences in the microenvironment. Here, we developed a statistical test called the Poisson Tree (PT) test to detect varying evolutionary rates among cell lineages, leveraging the phylogenetic signal of single-cell DNA sequencing (scDNA-seq) data. We applied the PT test to 24 healthy and cancer samples, rejecting a constant evolutionary rate in 11 out of 15 cancer and five out of nine healthy scDNA-seq datasets. In six cancer datasets, we identified subclonal mutations in known driver genes that could explain the rate accelerations of particular cancer lineages. Our findings demonstrate the efficacy of scDNA-seq for studying somatic evolution and suggest that cell lineages often evolve at different rates within cancer and healthy tissues.

SB Digestor: a tailored driver gene identification tool for dissecting heterogeneous Sleeping Beauty transposon-induced tumors.

Sleeping Beauty (SB) insertional mutagenesis has been widely used for genome-wide functional screening in mouse models of human cancers, however, intertumor heterogeneity can be a major obstacle in identifying common insertion sites (CISs). Although previous algorithms have been successful in defining some CISs, they also miss CISs in certain situations. A major common characteristic of these previous methods is that they do not take tumor heterogeneity into account. However, intertumoral heterogeneity directly influences the sequence read number for different tumor samples and then affects CIS identification. To precisely detect and define cancer driver genes, we developed SB Digestor, a computational algorithm that overcomes biological heterogeneity to identify more potential driver genes. Specifically, we define the relationship between the sequenced read number and putative gene number to deduce the depth cutoff for each tumor, which can reduce tumor complexity and precisely reflect intertumoral heterogeneity. Using this new tool, we re-analyzed our previously published SB-based screening dataset and identified many additional potent drivers involved in Brca1-related tumorigenesis, including Arhgap42, Tcf12, and Fgfr2. SB Digestor not only greatly enhances our ability to identify and prioritize cancer drivers from SB tumors but also substantially deepens our understanding of the intrinsic genetic basis of cancer.

Somatic variant calling from single-cell DNA sequencing data.

Single-cell sequencing has gained popularity in recent years. Despite its numerous applications, single-cell DNA sequencing data is highly error-prone due to technical biases arising from uneven sequencing coverage, allelic dropout, and amplification error. With these artifacts, the identification of somatic genomic variants becomes a challenging task, and over the years, several methods have been developed explicitly for this type of data. Single-cell variant callers implement distinct strategies, make different use of the data, and typically result in many discordant calls when applied to real data. Here, we review current approaches for single-cell variant calling, emphasizing single nucleotide variants. We highlight their potential benefits and shortcomings to help users choose a suitable tool for their data at hand.

SIEVE: joint inference of single-nucleotide variants and cell phylogeny from single-cell DNA sequencing data.

We present SIEVE, a statistical method for the joint inference of somatic variants and cell phylogeny under the finite-sites assumption from single-cell DNA sequencing. SIEVE leverages raw read counts for all nucleotides and corrects the acquisition bias of branch lengths. In our simulations, SIEVE outperforms other methods in phylogenetic reconstruction and variant calling accuracy, especially in the inference of homozygous variants. Applying SIEVE to three datasets, one for triple-negative breast (TNBC), and two for colorectal cancer (CRC), we find that double mutant genotypes are rare in CRC but unexpectedly frequent in the TNBC samples.

Clonality and timing of relapsing colorectal cancer metastasis revealed through whole-genome single-cell sequencing.

Recurrence of tumor cells following local and systemic therapy is a significant hurdle in cancer. Most patients with metastatic colorectal cancer (mCRC) will relapse, despite resection of the metastatic lesions. A better understanding of the evolutionary history of recurrent lesions is required to identify the spatial and temporal patterns of metastatic progression and expose the genetic and evolutionary determinants of therapeutic resistance. With this goal in mind, here we leveraged a unique single-cell whole-genome sequencing dataset from recurrent hepatic lesions of an mCRC patient. Our phylogenetic analysis confirms that the treatment induced a severe demographic bottleneck in the liver metastasis but also that a previously diverged lineage survived this surgery, possibly after migration to a different site in the liver. This lineage evolved very slowly for two years under adjuvant drug therapy and diversified again in a very short period. We identified several non-silent mutations specific to this lineage and inferred a substantial contribution of chemotherapy to the overall, genome-wide mutational burden. All in all, our study suggests that mCRC subclones can migrate locally and evade resection, keep evolving despite rounds of chemotherapy, and re-expand explosively.

Molecular landscape and subtype-specific therapeutic response of nasopharyngeal carcinoma revealed by integrative pharmacogenomics.

Nasopharyngeal carcinoma (NPC) is a malignant head and neck cancer type with high morbidity in Southeast Asia, however the pathogenic mechanism of this disease is poorly understood. Using integrative pharmacogenomics, we find that NPC subtypes maintain distinct molecular features, drug responsiveness, and graded radiation sensitivity. The epithelial carcinoma (EC) subtype is characterized by activations of microtubule polymerization and defective mitotic spindle checkpoint related genes, whereas sarcomatoid carcinoma (SC) and mixed sarcomatoid-epithelial carcinoma (MSEC) subtypes exhibit enriched epithelial-mesenchymal transition (EMT) and invasion promoting genes, which are well correlated with their morphological features. Furthermore, patient-derived organoid (PDO)-based drug test identifies potential subtype-specific treatment regimens, in that SC and MSEC subtypes are sensitive to microtubule inhibitors, whereas EC subtype is more responsive to EGFR inhibitors, which is synergistically enhanced by combining with radiotherapy. Through combinational chemoradiotherapy (CRT) screening, effective CRT regimens are also suggested for patients showing less sensitivity to radiation. Altogether, our study provides an example of applying integrative pharmacogenomics to establish a personalized precision oncology for NPC subtype-guided therapies.

NOTCH1 activation compensates BRCA1 deficiency and promotes triple-negative breast cancer formation.

BRCA1 mutation carriers have a higher risk of developing triple-negative breast cancer (TNBC), which is a refractory disease due to its non-responsiveness to current clinical targeted therapies. Using the Sleeping Beauty transposon system in Brca1-deficient mice, we identified 169 putative cancer drivers, among which Notch1 is a top candidate for accelerating TNBC by promoting the epithelial-mesenchymal transition (EMT) and regulating the cell cycle. Activation of NOTCH1 suppresses mitotic catastrophe caused by BRCA1 deficiency by restoring S/G2 and G2/M cell cycle checkpoints, which may through activation of ATR-CHK1 signalling pathway. Consistently, analysis of human breast cancer tissue demonstrates NOTCH1 is highly expressed in TNBCs, and the activated form of NOTCH1 correlates positively with increased phosphorylation of ATR. Additionally, we demonstrate that inhibition of the NOTCH1-ATR-CHK1 cascade together with cisplatin synergistically kills TNBC by targeting the cell cycle checkpoint, DNA damage and EMT, providing a potent clinical option for this fatal disease.

Apert's syndrome: Study by whole exome sequencing.

In the present study we attempted a parent-child trio, whole exome sequencing (WES) approach to study Apert's syndrome. Clinical characteristics of the child were noted down and WES was carried out using Ion Torrent System that revealed the presence of previously reported P253R mutation in FGFR2 gene. Presence of two SNPs rs1047057 and rs554851880 in FGFR2 gene with an allelic frequency of 0.5113 and 0.001176 respectively and 161 complete damaging mutations were found. This study is the first reported case of exome sequencing approach on an Apert's syndrome patient aimed at providing better genetic counselling in a non-consanguineous relationship.