Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing.

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

Mutational Strand Asymmetries in Cancer Genomes Reveal Mechanisms of DNA Damage and Repair.

Mutational processes constantly shape the somatic genome, leading to immunity, aging, cancer, and other diseases. When cancer is the outcome, we are afforded a glimpse into these processes by the clonal expansion of the malignant cell. Here, we characterize a less explored layer of the mutational landscape of cancer: mutational asymmetries between the two DNA strands. Analyzing whole-genome sequences of 590 tumors from 14 different cancer types, we reveal widespread asymmetries across mutagenic processes, with transcriptional ("T-class") asymmetry dominating UV-, smoking-, and liver-cancer-associated mutations and replicative ("R-class") asymmetry dominating POLE-, APOBEC-, and MSI-associated mutations. We report a striking phenomenon of transcription-coupled damage (TCD) on the non-transcribed DNA strand and provide evidence that APOBEC mutagenesis occurs on the lagging-strand template during DNA replication. As more genomes are sequenced, studying and classifying their asymmetries will illuminate the underlying biological mechanisms of DNA damage and repair.

Genetic variation in human DNA replication timing.

Genomic DNA replicates in a choreographed temporal order that impacts the distribution of mutations along the genome. We show here that DNA replication timing is shaped by genetic polymorphisms that act in cis upon megabase-scale DNA segments. In genome sequences from proliferating cells, read depth along chromosomes reflected DNA replication activity in those cells. We used this relationship to analyze variation in replication timing among 161 individuals sequenced by the 1000 Genomes Project. Genome-wide association of replication timing with genetic variation identified 16 loci at which inherited alleles associate with replication timing. We call these "replication timing quantitative trait loci" (rtQTLs). rtQTLs involved the differential use of replication origins, exhibited allele-specific effects on replication timing, and associated with gene expression variation at megabase scales. Our results show replication timing to be shaped by genetic polymorphism and identify a means by which inherited polymorphism regulates the mutability of nearby sequences.

Diverse mutational landscapes in human lymphocytes.

The lymphocyte genome is prone to many threats, including programmed mutation during differentiation1, antigen-driven proliferation and residency in diverse microenvironments. Here, after developing protocols for expansion of single-cell lymphocyte cultures, we sequenced whole genomes from 717 normal naive and memory B and T cells and haematopoietic stem cells. All lymphocyte subsets carried more point mutations and structural variants than haematopoietic stem cells, with higher burdens in memory cells than in naive cells, and with T cells accumulating mutations at a higher rate throughout life. Off-target effects of immunological diversification accounted for approximately half of the additional differentiation-associated mutations in lymphocytes. Memory B cells acquired, on average, 18 off-target mutations genome-wide for every on-target IGHV mutation during the germinal centre reaction. Structural variation was 16-fold higher in lymphocytes than in stem cells, with around 15% of deletions being attributable to off-target recombinase-activating gene activity. DNA damage from ultraviolet light exposure and other sporadic mutational processes generated hundreds to thousands of mutations in some memory cells. The mutation burden and signatures of normal B cells were broadly similar to those seen in many B-cell cancers, suggesting that malignant transformation of lymphocytes arises from the same mutational processes that are active across normal ontogeny. The mutational landscape of normal lymphocytes chronicles the off-target effects of programmed genome engineering during immunological diversification and the consequences of differentiation, proliferation and residency in diverse microenvironments.

Cell-free DNA profiling informs all major complications of hematopoietic cell transplantation.

Allogeneic hematopoietic cell transplantation (HCT) provides effective treatment for hematologic malignancies and immune disorders. Monitoring of posttransplant complications is critical, yet current diagnostic options are limited. Here, we show that cell-free DNA (cfDNA) in blood is a versatile analyte for monitoring of the most important complications that occur after HCT: graft-versus-host disease (GVHD), a frequent immune complication of HCT, infection, relapse of underlying disease, and graft failure. We demonstrate that these therapeutic complications are informed from a single assay, low-coverage bisulfite sequencing of cfDNA, followed by disease-specific bioinformatic analyses. To inform GVHD, we profile cfDNA methylation marks to trace the cfDNA tissues-of-origin and to quantify tissue-specific injury. To inform infection, we implement metagenomic cfDNA profiling. To inform cancer relapse, we implement analyses of tumor-specific genomic aberrations. Finally, to detect graft failure, we quantify the proportion of donor- and recipient-specific cfDNA. We applied this assay to 170 plasma samples collected from 27 HCT recipients at predetermined timepoints before and after allogeneic HCT. We found that the abundance of solid-organ-derived cfDNA in the blood at 1 mo after HCT is predictive of acute GVHD (area under the curve, 0.88). Metagenomic profiling of cfDNA revealed the frequent occurrence of viral reactivation in this patient population. The fraction of donor-specific cfDNA was indicative of relapse and remission, and the fraction of tumor-specific cfDNA was informative of cancer relapse. This proof-of-principle study shows that cfDNA has the potential to improve the care of allogeneic HCT recipients by enabling earlier detection and better prediction of the complex array of complications that occur after HCT.

Whole-exome sequencing of Nigerian benign prostatic hyperplasia reveals increased alterations in apoptotic pathways.

BACKGROUND: Through whole-exome sequencing of 60 formalin-fixed paraffin-embedded Nigerian (NGRn) benign prostatic hyperplasia (BPH) samples, we identified germline and somatic alterations in apoptotic pathways impacting BPH development and progression. Prostate enlargement is a common occurrence in male aging; however, this enlargement can lead to lower urinary tract symptoms that negatively impact quality of life. This impact is disproportionately present in men of African ancestry. BPH pathophysiology is poorly understood and studies examining non-European populations are lacking. METHODS: In this study, NGRn BPH, normal prostate, and prostate cancer (PCa) tumor samples were sequenced and compared to characterize genetic alterations in NGRn BPH. RESULTS: Two hundred and two nonbenign, ClinVar-annotated germline variants were present in NGRn BPH samples. Six genes [BRCA1 (92%), HSD3B1 (85%), TP53 (37%), PMS2 (23%), BARD1 (20%), and BRCA2 (17%)] were altered in at least 10% of samples; however, compared to NGRn normal and tumor, the frequency of alterations in BPH samples showed no significant differences at the gene or variant level. BRCA2_rs11571831 and TP53_rs1042522 germline alterations had a statistically significant co-occurrence interaction in BPH samples. In at least two BPH samples, 173 genes harbored somatic variants known to be clinically actionable. Three genes (COL18A1, KIF16B, and LRP1) showed a statistically significant (p < 0.05) higher frequency in BPH. NGRn BPH also had five gene pairs (PKD1/KIAA0100, PKHD1/PKD1, DNAH9/LRP1B, NWD1/DCHS2, and TCERG1/LMTK2) with statistically significant co-occurring interactions. Two hundred and seventy-nine genes contained novel somatic variants in NGRn BPH. Three genes (CABP1, FKBP1C, and RP11-595B24.2) had a statistically significant (p < 0.05) higher alteration frequency in NGRn BPH and three were significantly higher in NGRn tumor (CACNA1A, DMKN, and CACNA2D2). Pairwise Fisher's exact tests showed 14 gene pairs with statistically significant (p < 0.05) interactions and four interactions approaching significance (p < 0.10). Mutational patterns in NGRn BPH were similar to COSMIC (Catalog of Somatic Mutations in Cancer) signatures associated with aging and dysfunctional DNA damage repair. CONCLUSIONS: NGRn BPH contained significant germline alteration interactions (BRCA2_rs11571831 and TP53_rs1042522) and increased somatic alteration frequencies (LMTK2, LRP1, COL18A1, CABP1, and FKBP1C) that impact apoptosis. Normal prostate development is maintained by balancing apoptotic and proliferative activity. Dysfunction in either mechanism can lead to abnormal prostate growth. This work is the first to examine genomic sequencing in NGRn BPH and provides data that fill known gaps in the understanding BPH and how it impacts men of African ancestry.

Hybrid Stomach-Intestinal Chromatin States Underlie Human Barrett's Metaplasia.

BACKGROUND & AIMS: Tissue metaplasia is uncommon in adults because established cis-element programs resist rewiring. In Barrett's esophagus, the distal esophageal mucosa acquires a predominantly intestinal character, with notable gastric features, and is predisposed to developing invasive cancers. We sought to understand the chromatin underpinnings of Barrett's metaplasia and why it commonly displays simultaneous gastric and intestinal properties. METHODS: We profiled cis-regulatory elements with active histone modifications in primary human biopsy materials using chromatin immunoprecipitation followed by DNA sequencing. Mutations in Barrett's esophagus were examined in relation to tissue-specific enhancer landscapes using a random forest machine-learning algorithm. We also profiled open chromatin at single-cell resolution in primary Barrett's biopsy specimens using the assay for transposase-accessible chromatin. We used 1- and 2-color immunohistochemistry to examine protein expression of tissue-restricted genes. RESULTS: Barrett's esophagus bears epigenome fingerprints of human stomach and intestinal columnar, but not esophageal squamous, epithelia. Mutational patterns were best explained as arising on the epigenome background of active gastric cis-elements, supporting the view that adjoining stomach epithelium is a likely tissue source. Individual cells in Barrett's metaplasia coexpress gastric and intestinal genes, reflecting concomitant chromatin access at enhancers ordinarily restricted to one or the other epithelium. Protein expression of stomach-specific mucins; CLDN18; and a novel gastric marker, ANXA10, showed extensive tissue and subclonal heterogeneity of dual stomach-intestinal cell states. CONCLUSIONS: These findings reveal mixed and dynamic tissue-restricted chromatin states and phenotypic heterogeneity in Barrett's esophagus. Pervasive intragland variation argues against stem-cell governance of this phenotype.

Cell-of-origin chromatin organization shapes the mutational landscape of cancer.

Cancer is a disease potentiated by mutations in somatic cells. Cancer mutations are not distributed uniformly along the human genome. Instead, different human genomic regions vary by up to fivefold in the local density of cancer somatic mutations, posing a fundamental problem for statistical methods used in cancer genomics. Epigenomic organization has been proposed as a major determinant of the cancer mutational landscape. However, both somatic mutagenesis and epigenomic features are highly cell-type-specific. We investigated the distribution of mutations in multiple independent samples of diverse cancer types and compared them to cell-type-specific epigenomic features. Here we show that chromatin accessibility and modification, together with replication timing, explain up to 86% of the variance in mutation rates along cancer genomes. The best predictors of local somatic mutation density are epigenomic features derived from the most likely cell type of origin of the corresponding malignancy. Moreover, we find that cell-of-origin chromatin features are much stronger determinants of cancer mutation profiles than chromatin features of matched cancer cell lines. Furthermore, we show that the cell type of origin of a cancer can be accurately determined based on the distribution of mutations along its genome. Thus, the DNA sequence of a cancer genome encompasses a wealth of information about the identity and epigenomic features of its cell of origin.

Coronavirus 2019 and People Living With Human Immunodeficiency Virus: Outcomes for Hospitalized Patients in New York City.

BACKGROUND: There are limited data regarding the clinical impact of coronavirus disease 2019 (COVID-19) on people living with human immunodeficiency virus (PLWH). In this study, we compared outcomes for PLWH with COVID-19 to a matched comparison group. METHODS: We identified 88 PLWH hospitalized with laboratory-confirmed COVID-19 in our hospital system in New York City between 12 March and 23 April 2020. We collected data on baseline clinical characteristics, laboratory values, HIV status, treatment, and outcomes from this group and matched comparators (1 PLWH to up to 5 patients by age, sex, race/ethnicity, and calendar week of infection). We compared clinical characteristics and outcomes (death, mechanical ventilation, hospital discharge) for these groups, as well as cumulative incidence of death by HIV status. RESULTS: Patients did not differ significantly by HIV status by age, sex, or race/ethnicity due to the matching algorithm. PLWH hospitalized with COVID-19 had high proportions of HIV virologic control on antiretroviral therapy. PLWH had greater proportions of smoking (P < .001) and comorbid illness than uninfected comparators. There was no difference in COVID-19 severity on admission by HIV status (P = .15). Poor outcomes for hospitalized PLWH were frequent but similar to proportions in comparators; 18% required mechanical ventilation and 21% died during follow-up (compared with 23% and 20%, respectively). There was similar cumulative incidence of death over time by HIV status (P = .94). CONCLUSIONS: We found no differences in adverse outcomes associated with HIV infection for hospitalized COVID-19 patients compared with a demographically similar patient group.

Novel POLE pathogenic germline variant in a family with multiple primary tumors results in distinct mutational signatures.

We describe a family in which four siblings exhibited multiple or classic colonic polyposis with or without colorectal carcinoma (CRC). One female developed three primary tumors, including CRC and carcinomas of the ovary and breast. Whole-exome sequencing of germline DNA from affected and unaffected individuals revealed a novel missense mutation in the exonuclease domain of POLE (c.833C>A; p.Thr278Lys) associated with a highly penetrant, autosomal-dominant inheritance pattern. Functional studies in yeast and demonstration of a high mutational burden in the available tumors confirmed the pathogenicity of the novel variant. Prominent POLE-deficient somatic mutational signatures were seen in the CRCs, but in contrast, a mutational signature typical of concomitant tumoral loss of POLE and mismatch-repair function (POLE-exo* /MSI) was noted in the breast cancer. The breast cancer also showed distinctive pathological characteristics that reflect the presence of both the germline POLE variant and the secondary somatic MMR alterations.

Mutational heterogeneity in cancer and the search for new cancer-associated genes.

Major international projects are underway that are aimed at creating a comprehensive catalogue of all the genes responsible for the initiation and progression of cancer. These studies involve the sequencing of matched tumour-normal samples followed by mathematical analysis to identify those genes in which mutations occur more frequently than expected by random chance. Here we describe a fundamental problem with cancer genome studies: as the sample size increases, the list of putatively significant genes produced by current analytical methods burgeons into the hundreds. The list includes many implausible genes (such as those encoding olfactory receptors and the muscle protein titin), suggesting extensive false-positive findings that overshadow true driver events. We show that this problem stems largely from mutational heterogeneity and provide a novel analytical methodology, MutSigCV, for resolving the problem. We apply MutSigCV to exome sequences from 3,083 tumour-normal pairs and discover extraordinary variation in mutation frequency and spectrum within cancer types, which sheds light on mutational processes and disease aetiology, and in mutation frequency across the genome, which is strongly correlated with DNA replication timing and also with transcriptional activity. By incorporating mutational heterogeneity into the analyses, MutSigCV is able to eliminate most of the apparent artefactual findings and enable the identification of genes truly associated with cancer.

Genomic characterization of biliary tract cancers identifies driver genes and predisposing mutations.

BACKGROUND & AIMS: Biliary tract cancers (BTCs) are clinically and pathologically heterogeneous and respond poorly to treatment. Genomic profiling can offer a clearer understanding of their carcinogenesis, classification and treatment strategy. We performed large-scale genome sequencing analyses on BTCs to investigate their somatic and germline driver events and characterize their genomic landscape. METHODS: We analyzed 412 BTC samples from Japanese and Italian populations, 107 by whole-exome sequencing (WES), 39 by whole-genome sequencing (WGS), and a further 266 samples by targeted sequencing. The subtypes were 136 intrahepatic cholangiocarcinomas (ICCs), 101 distal cholangiocarcinomas (DCCs), 109 peri-hilar type cholangiocarcinomas (PHCs), and 66 gallbladder or cystic duct cancers (GBCs/CDCs). We identified somatic alterations and searched for driver genes in BTCs, finding pathogenic germline variants of cancer-predisposing genes. We predicted cell-of-origin for BTCs by combining somatic mutation patterns and epigenetic features. RESULTS: We identified 32 significantly and commonly mutated genes including TP53, KRAS, SMAD4, NF1, ARID1A, PBRM1, and ATR, some of which negatively affected patient prognosis. A novel deletion of MUC17 at 7q22.1 affected patient prognosis. Cell-of-origin predictions using WGS and epigenetic features suggest hepatocyte-origin of hepatitis-related ICCs. Deleterious germline mutations of cancer-predisposing genes such as BRCA1, BRCA2, RAD51D, MLH1, or MSH2 were detected in 11% (16/146) of BTC patients. CONCLUSIONS: BTCs have distinct genetic features including somatic events and germline predisposition. These findings could be useful to establish treatment and diagnostic strategies for BTCs based on genetic information. LAY SUMMARY: We here analyzed genomic features of 412 BTC samples from Japanese and Italian populations. A total of 32 significantly and commonly mutated genes were identified, some of which negatively affected patient prognosis, including a novel deletion of MUC17 at 7q22.1. Cell-of-origin predictions using WGS and epigenetic features suggest hepatocyte-origin of hepatitis-related ICCs. Deleterious germline mutations of cancer-predisposing genes were detected in 11% of patients with BTC. BTCs have distinct genetic features including somatic events and germline predisposition.

Patterns and rates of exonic de novo mutations in autism spectrum disorders.

Autism spectrum disorders (ASD) are believed to have genetic and environmental origins, yet in only a modest fraction of individuals can specific causes be identified. To identify further genetic risk factors, here we assess the role of de novo mutations in ASD by sequencing the exomes of ASD cases and their parents (n = 175 trios). Fewer than half of the cases (46.3%) carry a missense or nonsense de novo variant, and the overall rate of mutation is only modestly higher than the expected rate. In contrast, the proteins encoded by genes that harboured de novo missense or nonsense mutations showed a higher degree of connectivity among themselves and to previous ASD genes as indexed by protein-protein interaction screens. The small increase in the rate of de novo events, when taken together with the protein interaction results, are consistent with an important but limited role for de novo point mutations in ASD, similar to that documented for de novo copy number variants. Genetic models incorporating these data indicate that most of the observed de novo events are unconnected to ASD; those that do confer risk are distributed across many genes and are incompletely penetrant (that is, not necessarily sufficient for disease). Our results support polygenic models in which spontaneous coding mutations in any of a large number of genes increases risk by 5- to 20-fold. Despite the challenge posed by such models, results from de novo events and a large parallel case-control study provide strong evidence in favour of CHD8 and KATNAL2 as genuine autism risk factors.