Archival single-cell genomics reveals persistent subclones during DCIS progression.

Ductal carcinoma in situ (DCIS) is a common precursor of invasive breast cancer. Our understanding of its genomic progression to recurrent disease remains poor, partly due to challenges associated with the genomic profiling of formalin-fixed paraffin-embedded (FFPE) materials. Here, we developed Arc-well, a high-throughput single-cell DNA-sequencing method that is compatible with FFPE materials. We validated our method by profiling 40,330 single cells from cell lines, a frozen tissue, and 27 FFPE samples from breast, lung, and prostate tumors stored for 3-31 years. Analysis of 10 patients with matched DCIS and cancers that recurred 2-16 years later show that many primary DCIS had already undergone whole-genome doubling and clonal diversification and that they shared genomic lineages with persistent subclones in the recurrences. Evolutionary analysis suggests that most DCIS cases in our cohort underwent an evolutionary bottleneck, and further identified chromosome aberrations in the persistent subclones that were associated with recurrence.

Genomic characterization of metastatic patterns from prospective clinical sequencing of 25,000 patients.

Metastatic progression is the main cause of death in cancer patients, whereas the underlying genomic mechanisms driving metastasis remain largely unknown. Here, we assembled MSK-MET, a pan-cancer cohort of over 25,000 patients with metastatic diseases. By analyzing genomic and clinical data from this cohort, we identified associations between genomic alterations and patterns of metastatic dissemination across 50 tumor types. We found that chromosomal instability is strongly correlated with metastatic burden in some tumor types, including prostate adenocarcinoma, lung adenocarcinoma, and HR+/HER2+ breast ductal carcinoma, but not in others, including colorectal cancer and high-grade serous ovarian cancer, where copy-number alteration patterns may be established early in tumor development. We also identified somatic alterations associated with metastatic burden and specific target organs. Our data offer a valuable resource for the investigation of the biological basis for metastatic spread and highlight the complex role of chromosomal instability in cancer progression.

Data Sanitization to Reduce Private Information Leakage from Functional Genomics.

The generation of functional genomics datasets is surging, because they provide insight into gene regulation and organismal phenotypes (e.g., genes upregulated in cancer). The intent behind functional genomics experiments is not necessarily to study genetic variants, yet they pose privacy concerns due to their use of next-generation sequencing. Moreover, there is a great incentive to broadly share raw reads for better statistical power and general research reproducibility. Thus, we need new modes of sharing beyond traditional controlled-access models. Here, we develop a data-sanitization procedure allowing raw functional genomics reads to be shared while minimizing privacy leakage, enabling principled privacy-utility trade-offs. Our protocol works with traditional Illumina-based assays and newer technologies such as 10x single-cell RNA sequencing. It involves quantifying the privacy leakage in reads by statistically linking study participants to known individuals. We carried out these linkages using data from highly accurate reference genomes and more realistic environmental samples.

Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing.

Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

Sexually dimorphic DNA damage responses and mutation avoidance in the mouse germline.

Germ cells specified during fetal development form the foundation of the mammalian germline. These primordial germ cells (PGCs) undergo rapid proliferation, yet the germline is highly refractory to mutation accumulation compared with somatic cells. Importantly, while the presence of endogenous or exogenous DNA damage has the potential to impact PGCs, there is little known about how these cells respond to stressors. To better understand the DNA damage response (DDR) in these cells, we exposed pregnant mice to ionizing radiation (IR) at specific gestational time points and assessed the DDR in PGCs. Our results show that PGCs prior to sex determination lack a G1 cell cycle checkpoint. Additionally, the response to IR-induced DNA damage differs between female and male PGCs post-sex determination. IR of female PGCs caused uncoupling of germ cell differentiation and meiotic initiation, while male PGCs exhibited repression of piRNA metabolism and transposon derepression. We also used whole-genome single-cell DNA sequencing to reveal that genetic rescue of DNA repair-deficient germ cells (Fancm-/- ) leads to increased mutation incidence and biases. Importantly, our work uncovers novel insights into how PGCs exposed to DNA damage can become developmentally defective, leaving only those genetically fit cells to establish the adult germline.

The impact of age and number of mutations on the size of clonal hematopoiesis.

Clonal hematopoiesis (CH) represents the clonal expansion of hematopoietic stem cells and their progeny driven by somatic mutations. Accurate risk assessment of CH is critical for disease prevention and clinical decision-making. The size of CH has been showed to associate with higher disease risk, yet, factors influencing the size of CH are unknown. In addition, the characteristics of CH in long-lived individuals are not well documented. Here, we report an in-depth analysis of CH in longevous (≥90 y old) and common (60~89 y old) elderly groups. Utilizing targeted deep sequencing, we found that the development of CH is closely related to age and the expression of aging biomarkers. The longevous elderly group exhibited a significantly higher incidence of CH and significantly higher frequency of TET2 and ASXL1 mutations, suggesting that certain CH could be beneficial to prolong life. Intriguingly, the size of CH neither correlates significantly to age, in the range of 60 to 110 y old, nor to the expression of aging biomarkers. Instead, we identified a strong correlation between large CH size and the number of mutations per individual. These findings provide a risk assessment biomarker for CH and also suggest that the evolution of the CH is influenced by factor(s) in addition to age.

CHARR efficiently estimates contamination from DNA sequencing data.

DNA sample contamination is a major issue in clinical and research applications of whole-genome and -exome sequencing. Even modest levels of contamination can substantially affect the overall quality of variant calls and lead to widespread genotyping errors. Currently, popular tools for estimating the contamination level use short-read data (BAM/CRAM files), which are expensive to store and manipulate and often not retained or shared widely. We propose a metric to estimate DNA sample contamination from variant-level whole-genome and -exome sequence data called CHARR, contamination from homozygous alternate reference reads, which leverages the infiltration of reference reads within homozygous alternate variant calls. CHARR uses a small proportion of variant-level genotype information and thus can be computed from single-sample gVCFs or callsets in VCF or BCF formats, as well as efficiently stored variant calls in Hail VariantDataset format. Our results demonstrate that CHARR accurately recapitulates results from existing tools with substantially reduced costs, improving the accuracy and efficiency of downstream analyses of ultra-large whole-genome and exome sequencing datasets.

Somatic Mutagenesis in Mammals and Its Implications for Human Disease and Aging.

DNA mutations as a consequence of errors during DNA damage repair, replication, or mitosis are the substrate for evolution. In multicellular organisms, mutations can occur in the germline and also in somatic tissues, where they are associated with cancer and other chronic diseases and possibly with aging. Recent advances in high-throughput sequencing have made it relatively easy to study germline de novo mutations, but in somatic cells, the vast majority of mutations are low-abundant and can be detected only in clonal lineages, such as tumors, or single cells. Here we review recent results on somatic mutations in normal human and animal tissues with a focus on their possible functional consequences.

Systems Biology of the Vasopressin V2 Receptor: New Tools for Discovery of Molecular Actions of a GPCR.

Systems biology can be defined as the study of a biological process in which all of the relevant components are investigated together in parallel to discover the mechanism. Although the approach is not new, it has come to the forefront as a result of genome sequencing projects completed in the first few years of the current century. It has elements of large-scale data acquisition (chiefly next-generation sequencing-based methods and protein mass spectrometry) and large-scale data analysis (big data integration and Bayesian modeling). Here we discuss these methodologies and show how they can be applied to understand the downstream effects of GPCR signaling, specifically looking at how the neurohypophyseal peptide hormone vasopressin, working through the V2 receptor and PKA activation, regulates the water channel aquaporin-2. The emerging picture provides a detailedframework for understanding the molecular mechanisms involved in water balance disorders, pointing the way to improved treatment of both polyuric disorders and water-retention disorders causing dilutional hyponatremia.

Earth Biogenome Project: present status and future plans.

Earth Biogenome Project (EBP) is an ambitious project targeted to provide high-quality reference genome sequences for all 1.8 million named extant (living) eukaryote species. The project was launched on 1 November 2018 with an initial 2 years' pilot phase (2018-2020) followed by Phase I (2020-2023), during which genomes of 9400 species will be sequenced. The genomes of the remaining ~1.7 million species will be sequenced in a planned manner during Phase II (2024-2027) and Phase III (2028-2030). In view of the excitement generated and the progress already made, the subject was covered in a Special Feature of a recent issue of PNAS (25 January 2022). The present status and future plans of EBP along with challenges faced are briefly discussed in this article.

What Is Metagenomics Teaching Us, and What Is Missed?

Shotgun metagenomic sequencing has revolutionized our ability to detect and characterize the diversity and function of complex microbial communities. In this review, we highlight the benefits of using metagenomics as well as the breadth of conclusions that can be made using currently available analytical tools, such as greater resolution of species and strains across phyla and functional content, while highlighting challenges of metagenomic data analysis. Major challenges remain in annotating function, given the dearth of functional databases for environmental bacteria compared to model organisms, and the technical difficulties of metagenome assembly and phasing in heterogeneous environmental samples. In the future, improvements and innovation in technology and methodology will lead to lowered costs. Data integration using multiple technological platforms will lead to a better understanding of how to harness metagenomes. Subsequently, we will be able not only to characterize complex microbiomes but also to manipulate communities to achieve prosperous outcomes for health, agriculture, and environmental sustainability.

Huntington's disease affects mitochondrial network dynamics predisposing to pathogenic mitochondrial DNA mutations.

Huntington's disease (HD) predominantly affects the brain, causing a mixed movement disorder, cognitive decline and behavioural abnormalities. It also causes a peripheral phenotype involving skeletal muscle. Mitochondrial dysfunction has been reported in tissues of HD models, including skeletal muscle, and lymphoblast and fibroblast cultures from patients with HD. Mutant huntingtin protein (mutHTT) expression can impair mitochondrial quality control and accelerate mitochondrial ageing. Here, we obtained fresh human skeletal muscle, a post-mitotic tissue expressing the mutated HTT allele at physiological levels since birth, and primary cell lines from HTT CAG repeat expansion mutation carriers and matched healthy volunteers to examine whether such a mitochondrial phenotype exists in human HD. Using ultra-deep mitochondrial DNA (mtDNA) sequencing, we showed an accumulation of mtDNA mutations affecting oxidative phosphorylation. Tissue proteomics indicated impairments in mtDNA maintenance with increased mitochondrial biogenesis of less efficient oxidative phosphorylation (lower complex I and IV activity). In full-length mutHTT expressing primary human cell lines, fission-inducing mitochondrial stress resulted in normal mitophagy. In contrast, expression of high levels of N-terminal mutHTT fragments promoted mitochondrial fission and resulted in slower, less dynamic mitophagy. Expression of high levels of mutHTT fragments due to somatic nuclear HTT CAG instability can thus affect mitochondrial network dynamics and mitophagy, leading to pathogenic mtDNA mutations. We show that life-long expression of mutant HTT causes a mitochondrial phenotype indicative of mtDNA instability in fresh post-mitotic human skeletal muscle. Thus, genomic instability may not be limited to nuclear DNA, where it results in somatic expansion of the HTT CAG repeat length in particularly vulnerable cells such as striatal neurons. In addition to efforts targeting the causative mutation, promoting mitochondrial health may be a complementary strategy in treating diseases with DNA instability such as HD.