Most large structural variants in cancer genomes can be detected without long reads.

Short-read sequencing is the workhorse of cancer genomics yet is thought to miss many structural variants (SVs), particularly large chromosomal alterations. To characterize missing SVs in short-read whole genomes, we analyzed 'loose ends'-local violations of mass balance between adjacent DNA segments. In the landscape of loose ends across 1,330 high-purity cancer whole genomes, most large (>10-kb) clonal SVs were fully resolved by short reads in the 87% of the human genome where copy number could be reliably measured. Some loose ends represent neotelomeres, which we propose as a hallmark of the alternative lengthening of telomeres phenotype. These pan-cancer findings were confirmed by long-molecule profiles of 38 breast cancer and melanoma cases. Our results indicate that aberrant homologous recombination is unlikely to drive the majority of large cancer SVs. Furthermore, analysis of mass balance in short-read whole genome data provides a surprisingly complete picture of cancer chromosomal structure.

Long-molecule scars of backup DNA repair in BRCA1- and BRCA2-deficient cancers.

Homologous recombination (HR) deficiency is associated with DNA rearrangements and cytogenetic aberrations1. Paradoxically, the types of DNA rearrangements that are specifically associated with HR-deficient cancers only minimally affect chromosomal structure2. Here, to address this apparent contradiction, we combined genome-graph analysis of short-read whole-genome sequencing (WGS) profiles across thousands of tumours with deep linked-read WGS of 46 BRCA1- or BRCA2-mutant breast cancers. These data revealed a distinct class of HR-deficiency-enriched rearrangements called reciprocal pairs. Linked-read WGS showed that reciprocal pairs with identical rearrangement orientations gave rise to one of two distinct chromosomal outcomes, distinguishable only with long-molecule data. Whereas one (cis) outcome corresponded to the copying and pasting of a small segment to a distant site, a second (trans) outcome was a quasi-balanced translocation or multi-megabase inversion with substantial (10 kb) duplications at each junction. We propose an HR-independent replication-restart repair mechanism to explain the full spectrum of reciprocal pair outcomes. Linked-read WGS also identified single-strand annealing as a repair pathway that is specific to BRCA2 deficiency in human cancers. Integrating these features in a classifier improved discrimination between BRCA1- and BRCA2-deficient genomes. In conclusion, our data reveal classes of rearrangements that are specific to BRCA1 or BRCA2 deficiency as a source of cytogenetic aberrations in HR-deficient cells.

Genomic signature of Fanconi anaemia DNA repair pathway deficiency in cancer.

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage1-3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes4-7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas8 (SCCs). Molecular studies of SCCs from individuals with FA (FA SCCs) are limited, and it is unclear how FA SCCs relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or infection with human papillomavirus9 (HPV). Here, by sequencing genomes and exomes of FA SCCs, we demonstrate that the primary genomic signature of FA repair deficiency is the presence of high numbers of structural variants. Structural variants are enriched for small deletions, unbalanced translocations and fold-back inversions, and are often connected, thereby forming complex rearrangements. They arise in the context of TP53 loss, but not in the context of HPV infection, and lead to somatic copy-number alterations of HNSCC driver genes. We further show that FA pathway deficiency may lead to epithelial-to-mesenchymal transition and enhanced keratinocyte-intrinsic inflammatory signalling, which would contribute to the aggressive nature of FA SCCs. We propose that the genomic instability in sporadic HPV-negative HNSCC may arise as a result of the FA repair pathway being overwhelmed by DNA interstrand crosslink damage caused by alcohol and tobacco-derived aldehydes, making FA SCC a powerful model to study tumorigenesis resulting from DNA-crosslinking damage.

Somatic whole genome dynamics of precancer in Barrett's esophagus reveals features associated with disease progression.

While the genomes of normal tissues undergo dynamic changes over time, little is understood about the temporal-spatial dynamics of genomes in premalignant tissues that progress to cancer compared to those that remain cancer-free. Here we use whole genome sequencing to contrast genomic alterations in 427 longitudinal samples from 40 patients with stable Barrett's esophagus compared to 40 Barrett's patients who progressed to esophageal adenocarcinoma (ESAD). We show the same somatic mutational processes are active in Barrett's tissue regardless of outcome, with high levels of mutation, ESAD gene and focal chromosomal alterations, and similar mutational signatures. The critical distinction between stable Barrett's versus those who progress to cancer is acquisition and expansion of TP53-/- cell populations having complex structural variants and high-level amplifications, which are detectable up to six years prior to a cancer diagnosis. These findings reveal the timing of common somatic genome dynamics in stable Barrett's esophagus and define key genomic features specific to progression to esophageal adenocarcinoma, both of which are critical for cancer prevention and early detection strategies.

Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Next-generation characterization of the Cancer Cell Line Encyclopedia.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

Impact of Induction Therapy on Circulating T Follicular Helper Cells and Subsequent Donor-Specific Antibody Formation After Kidney Transplant.

INTRODUCTION: The cellular events that contribute to generation of donor-specific anti-HLA antibodies (DSA) post-kidney transplantation (KTx) are not well understood. Characterization of such mechanisms could allow tailoring of immunosuppression to benefit sensitized patients. METHODS: We prospectively monitored circulating T follicular helper (cTFH) cells in KTx recipients who received T-cell depleting (thymoglobulin, n = 54) or T-cell nondepleting (basiliximab, n = 20) induction therapy from pre-KTx to 1 year post-KTx and assessed their phenotypic changes due to induction and DSA occurrence, in addition to healthy controls (n = 13), for a total of 307 blood samples. RESULTS: Before KTx, patients displayed comparable levels of resting, central memory cTFH cells with similar polarization to those of healthy controls. Unlike basiliximab induction, thymoglobulin induction significantly depleted cTFH cells, triggered lymphopenia-induced proliferation that skewed cTFH cells toward increased Th1 polarization, effector memory, and elevated programmed cell death protein 1 (PD-1)int/hi expression, resembling activated phenotypes. Regardless of induction, patients who developed DSA post-KTx, harbored pre-KTx donor-reactive memory interleukin (IL)-21+ cTFH cells and showed higher % cTFH and lower % of T regulatory (TREG) cells post-KTx resulting in elevated cTFH:TREG ratio at DSA occurrence. CONCLUSION: Induction therapy distinctly shapes cTFH cell phenotype post-KTx. Monitoring cTFH cells before and after KTx may help detect those patients prone to DSA generation post-KTx.

Phenotype, function, and differentiation potential of human monocyte subsets.

Human monocytes have been grouped into classical (CD14++CD16-), non-classical (CD14dimCD16++), and intermediate (CD14++CD16+) subsets. Documentation of normal function and variation in this complement of subtypes, particularly their differentiation potential to dendritic cells (DC) or macrophages, remains incomplete. We therefore phenotyped monocytes from peripheral blood of healthy subjects and performed functional studies on high-speed sorted subsets. Subset frequencies were found to be tightly controlled over time and across individuals. Subsets were distinct in their secretion of TNFα, IL-6, and IL-1ß in response to TLR agonists, with classical monocytes being the most producers and non-classical monocytes the least. Monocytes, particularly those of the non-classical subtype, secreted interferon-α (IFN-α) in response to intracellular TLR3 stimulation. After incubation with IL-4 and GM-CSF, classical monocytes acquired monocyte-derived DC (mo-DC) markers and morphology and stimulated allogeneic T cell proliferation in MLR; intermediate and non-classical monocytes did not. After incubation with IL-3 and Flt3 ligand, no subset differentiated to plasmacytoid DC. After incubation with GM-CSF (M1 induction) or macrophage colony-stimulating factor (M-CSF) (M2 induction), all subsets acquired macrophage morphology, secreted macrophage-associated cytokines, and displayed enhanced phagocytosis. From these studies we conclude that classical monocytes are the principal source of mo-DCs, but all subsets can differentiate to macrophages. We also found that monocytes, in particular the non-classical subset, represent an alternate source of type I IFN secretion in response to virus-associated TLR agonists.

Human immunodeficiency virus type 1 Vpr polymorphisms associated with progressor and nonprogressor individuals alter Vpr-associated functions.

Following infection with Human immunodeficiency virus 1 (HIV-1) there is a remarkable variation in virus replication and disease progression. Both host and viral factors have been implicated in the observed differences in disease status. Here, we focus on understanding the contribution of HIV-1 viral protein R (Vpr) by evaluating the disease-associated Vpr polymorphism and its biological functions from HIV-1 positive rapid progressor (RP) and long-term nonprogressor (LTNP) subjects. Results presented here show distinct variation in phenotypes of Vpr alleles from LTNP and RP subjects. Most notably, the polymorphism of Vpr at R36W and L68M associated with RP shows higher levels of oligomerization, and increased virus replication, whereas R77Q exhibits poor replication kinetics. Interestingly, we did not observe correlation with cell cycle arrest function. Together these results indicate that polymorphisms in Vpr in part may contribute to altered virus replication kinetics leading to the observed differences in disease progression in LTNP and RP groups.