Hypermutation of the inactive X chromosome is a frequent event in cancer.

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating cells.

Copy Number Analysis in a Large Cohort Suggestive of Inborn Errors of Immunity Indicates a Wide Spectrum of Relevant Chromosomal Losses and Gains.

Inborn errors of immunity (IEI) are genetically driven disorders. With the advancement of sequencing technologies, a rapidly increasing number of gene defects has been identified, thereby mirroring the high heterogeneity in immunological and clinical presentations observed in patients. However, for a large majority of patients, no causative single nucleotide variant (SNV) or small indel can be identified using next-generation sequencing. First studies have shown that also copy number variants (CNVs) can cause IEI. Unfortunately, CNVs are not well examined in many routine diagnostic settings and the aim of this study was to assess the number of clinically relevant chromosomal losses and gains in a large cohort. We identified a total of 20 CNVs using whole exome sequencing data of a cohort of 191 patients with a suspected IEI. A definite molecular diagnosis could be made in five patients (2.6%), including pathogenic deletions affecting ICOS, TNFAIP3, and 22q11.2. CNVs of uncertain significance were observed in fifteen patients (7.9%), including deletions of 11q22.1q22.3 and 16p11.2 but also duplications affecting entire or parts of genes previously associated with IEI. Importantly, five patients carrying a CNV of uncertain significance also carried pathogenic or likely pathogenic SNVs (PIK3R1, NFKB1, NLRC4, DOCK2), or SNVs of unknown significance (NFKB2). This cooccurrence of SNVs and CNVs suggests modifying effects in some patients, and functional follow-up is warranted now in order to better understand phenotypic heterogeneity. In summary, the diagnostic yield of IEI can be increased substantially by evaluating CNVs, which allows an improved therapeutic management in those patients.

Detection of CSF1 rearrangements deleting the 3' UTR in tenosynovial giant cell tumors.

Tenosynovial giant cell tumors (TGCTs) are characterized by rearrangements of CSF1, thought to drive overexpression of macrophage colony-stimulating factor (CSF1), thereby promoting tumor growth and recruitment of non-neoplastic mononuclear and multinucleated inflammatory cells. While fusions to collagen promoters have been described, the mechanism of CSF1 overexpression has been unclear in a majority of cases. Two cohorts of TGCT were investigated for CSF1 rearrangements using fluorescence in situ hybridization (FISH) and either RNA-seq or DNA-seq with Sanger validation. The study comprised 39 patients, including 13 localized TGCT, 21 diffuse TGCT, and five of unspecified type. CSF1 rearrangements were identified by FISH in 30 cases: 13 translocations, 17 3' deletions. Sequencing confirmed CSF1 breakpoints in 28 cases; in all 28 the breakpoint was found to be downstream of exon 5, replacing or deleting a long 3' UTR containing known miRNA and AU-rich element negative regulatory sequences. We also confirmed the presence of CBL exon 8-9 mutations in six of 21 cases. In conclusion, TGCT in our large cohort were characterized by variable alterations, all of which led to truncation of the 3' end of CSF1, instead of the COL6A3-CSF1 fusions previously reported in some TGCTs. The diversity of fusion partners but consistent integrity of CSF1 functional domains encoded by exons 1-5 support a hypothesis that CSF1 overexpression results from transcription of a truncated form of CSF1 lacking 3' negative regulatory sequences. The presence of CBL mutations affecting the linker and RING finger domain suggests an alternative mechanism for increased CSF1/CSF1R signaling in some cases.

Duplicated Enhancer Region Increases Expression of CTSB and Segregates with Keratolytic Winter Erythema in South African and Norwegian Families.

Keratolytic winter erythema (KWE) is a rare autosomal-dominant skin disorder characterized by recurrent episodes of palmoplantar erythema and epidermal peeling. KWE was previously mapped to 8p23.1-p22 (KWE critical region) in South African families. Using targeted resequencing of the KWE critical region in five South African families and SNP array and whole-genome sequencing in two Norwegian families, we identified two overlapping tandem duplications of 7.67 kb (South Africans) and 15.93 kb (Norwegians). The duplications segregated with the disease and were located upstream of CTSB, a gene encoding cathepsin B, a cysteine protease involved in keratinocyte homeostasis. Included in the 2.62 kb overlapping region of these duplications is an enhancer element that is active in epidermal keratinocytes. The activity of this enhancer correlated with CTSB expression in normal differentiating keratinocytes and other cell lines, but not with FDFT1 or NEIL2 expression. Gene expression (qPCR) analysis and immunohistochemistry of the palmar epidermis demonstrated significantly increased expression of CTSB, as well as stronger staining of cathepsin B in the stratum granulosum of affected individuals than in that of control individuals. Analysis of higher-order chromatin structure data and RNA polymerase II ChIA-PET data from MCF-7 cells did not suggest remote effects of the enhancer. In conclusion, KWE in South African and Norwegian families is caused by tandem duplications in a non-coding genomic region containing an active enhancer element for CTSB, resulting in upregulation of this gene in affected individuals.

Decoding the regulatory landscape of medulloblastoma using DNA methylation sequencing.

Epigenetic alterations, that is, disruption of DNA methylation and chromatin architecture, are now acknowledged as a universal feature of tumorigenesis. Medulloblastoma, a clinically challenging, malignant childhood brain tumour, is no exception. Despite much progress from recent genomics studies, with recurrent changes identified in each of the four distinct tumour subgroups (WNT-pathway-activated, SHH-pathway-activated, and the less-well-characterized Group 3 and Group 4), many cases still lack an obvious genetic driver. Here we present whole-genome bisulphite-sequencing data from thirty-four human and five murine tumours plus eight human and three murine normal controls, augmented with matched whole-genome, RNA and chromatin immunoprecipitation sequencing data. This comprehensive data set allowed us to decipher several features underlying the interplay between the genome, epigenome and transcriptome, and its effects on medulloblastoma pathophysiology. Most notable were highly prevalent regions of hypomethylation correlating with increased gene expression, extending tens of kilobases downstream of transcription start sites. Focal regions of low methylation linked to transcription-factor-binding sites shed light on differential transcriptional networks between subgroups, whereas increased methylation due to re-normalization of repressed chromatin in DNA methylation valleys was positively correlated with gene expression. Large, partially methylated domains affecting up to one-third of the genome showed increased mutation rates and gene silencing in a subgroup-specific fashion. Epigenetic alterations also affected novel medulloblastoma candidate genes (for example, LIN28B), resulting in alternative promoter usage and/or differential messenger RNA/microRNA expression. Analysis of mouse medulloblastoma and precursor-cell methylation demonstrated a somatic origin for many alterations. Our data provide insights into the epigenetic regulation of transcription and genome organization in medulloblastoma pathogenesis, which are probably also of importance in a wider developmental and disease context.

Antimalarial Activity of KAF156 in Falciparum and Vivax Malaria.

BACKGROUND: KAF156 belongs to a new class of antimalarial agents (imidazolopiperazines), with activity against asexual and sexual blood stages and the preerythrocytic liver stages of malarial parasites. METHODS: We conducted a phase 2, open-label, two-part study at five centers in Thailand and Vietnam to assess the antimalarial efficacy, safety, and pharmacokinetic profile of KAF156 in adults with acute Plasmodium vivax or P. falciparum malaria. Assessment of parasite clearance rates in cohorts of patients with vivax or falciparum malaria who were treated with multiple doses (400 mg once daily for 3 days) was followed by assessment of the cure rate at 28 days in a separate cohort of patients with falciparum malaria who received a single dose (800 mg). RESULTS: Median parasite clearance times were 45 hours (interquartile range, 42 to 48) in 10 patients with falciparum malaria and 24 hours (interquartile range, 20 to 30) in 10 patients with vivax malaria after treatment with the multiple-dose regimen and 49 hours (interquartile range, 42 to 54) in 21 patients with falciparum malaria after treatment with the single dose. Among the 21 patients who received the single dose and were followed for 28 days, 1 had reinfection and 7 had recrudescent infections (cure rate, 67%; 95% credible interval, 46 to 84). The mean (±SD) KAF156 terminal elimination half-life was 44.1±8.9 hours. There were no serious adverse events in this small study. The most common adverse events included sinus bradycardia, thrombocytopenia, hypokalemia, anemia, and hyperbilirubinemia. Vomiting of grade 2 or higher occurred in 2 patients, 1 of whom discontinued treatment because of repeated vomiting after receiving the single 800-mg dose. More adverse events were reported in the single-dose cohort, which had longer follow-up, than in the multiple-dose cohorts. CONCLUSIONS: KAF156 showed antimalarial activity without evident safety concerns in a small number of adults with uncomplicated P. vivax or P. falciparum malaria. (Funded by Novartis and others; ClinicalTrials.gov number, NCT01753323 .).

Transcriptome and genome sequencing uncovers functional variation in humans.

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Dissecting the genomic complexity underlying medulloblastoma.

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour-normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.

A comprehensive assessment of RNA-seq protocols for degraded and low-quantity samples.

BACKGROUND: RNA-sequencing (RNA-seq) has emerged as one of the most sensitive tool for gene expression analysis. Among the library preparation methods available, the standard poly(A) + enrichment provides a comprehensive, detailed, and accurate view of polyadenylated RNAs. However, on samples of suboptimal quality ribosomal RNA depletion and exon capture methods have recently been reported as better alternatives. METHODS: We compared for the first time three commercial Illumina library preparation kits (TruSeq Stranded mRNA, TruSeq Ribo-Zero rRNA Removal, and TruSeq RNA Access) as representatives of these three different approaches using well-established human reference RNA samples from the MAQC/SEQC consortium on a wide range of input amounts (from 100 ng down to 1 ng) and degradation levels (intact, degraded, and highly degraded). RESULTS: We assessed the accuracy of the generated expression values by comparison to gold standard TaqMan qPCR measurements and gained unprecedented insight into the limits of applicability in terms of input quantity and sample quality of each protocol. We found that each protocol generates highly reproducible results (R 2 > 0.92) on intact RNA samples down to input amounts of 10 ng. For degraded RNA samples, Ribo-Zero showed clear performance advantages over the other two protocols as it generated more accurate and better reproducible gene expression results even at very low input amounts such as 1 ng and 2 ng. For highly degraded RNA samples, RNA Access performed best generating reliable data down to 5 ng input. CONCLUSIONS: We found that the ribosomal RNA depletion protocol from Illumina works very well at amounts far below recommendation and over a good range of intact and degraded material. We also infer that the exome-capture protocol (RNA Access, Illumina) performs better than other methods on highly degraded and low amount samples.

Novel mutation in the CHST6 gene causes macular corneal dystrophy in a black South African family.

BACKGROUND: Macular corneal dystrophy (MCD) is a rare autosomal recessive disorder that is characterized by progressive corneal opacity that starts in early childhood and ultimately progresses to blindness in early adulthood. The aim of this study was to identify the cause of MCD in a black South African family with two affected sisters. METHODS: A multigenerational South African Sotho-speaking family with type I MCD was studied using whole exome sequencing. Variant filtering to identify the MCD-causal mutation included the disease inheritance pattern, variant minor allele frequency and potential functional impact. RESULTS: Ophthalmologic evaluation of the cases revealed a typical MCD phenotype and none of the other family members were affected. An average of 127 713 variants per individual was identified following exome sequencing and approximately 1.2 % were not present in any of the investigated public databases. Variant filtering identified a homozygous E71Q mutation in CHST6, a known MCD-causing gene encoding corneal N-acetyl glucosamine-6-O-sulfotransferase. This E71Q mutation results in a non-conservative amino acid change in a highly conserved functional domain of the human CHST6 that is essential for enzyme activity. CONCLUSION: We identified a novel E71Q mutation in CHST6 as the MCD-causal mutation in a black South African family with type I MCD. This is the first description of MCD in a black Sub-Saharan African family and therefore contributes valuable insights into the genetic aetiology of this disease, while improving genetic counselling for this and potentially other MCD families.

Influence of RNA extraction methods and library selection schemes on RNA-seq data.

BACKGROUND: Gene expression analysis by RNA sequencing is now widely used in a number of applications surveying the whole transcriptomes of cells and tissues. The recent introduction of ribosomal RNA depletion protocols, such as RiboZero, has extended the view of the polyadenylated transcriptome to the poly(A)- fraction of the RNA. However, substantial amounts of intronic transcriptional activity has been reported in RiboZero protocols, raising issues regarding their potential nuclear origin and the impact on the actual sequence depth in exonic regions. RESULTS: Using HEK293 human cells as source material, we assessed here the impact of the two commonly used RNA extraction methods and of the library construction protocols (rRNA depletion versus mRNA) on 1) the relative abundance of intronic reads and 2) on the estimation of gene expression values. We benchmarked the rRNA depletion-based sequencing with a specific analysis of the cytoplasmic and nuclear transcriptome fractions, suggesting that the large majority of the intronic reads correspond to unprocessed nuclear transcripts rather than to independent transcriptional units. We show that Qiagen or TRIzol extraction methods retain differentially nuclear RNA species, and that consequently, rRNA depletion-based RNA sequencing protocols are particularly sensitive to the extraction methods. CONCLUSIONS: We could show that the combination of Trizol-based RNA extraction with rRNA depletion sequencing protocols led to the largest fraction of intronic reads, after the sequencing of the nuclear transcriptome. We discuss here the impact of the various strategies on gene expression and alternative splicing estimation measures. Further, we propose guidelines and a double selection strategy for minimizing the expression biases, without loss of information.

A high-resolution anatomical atlas of the transcriptome in the mouse embryo.

Ascertaining when and where genes are expressed is of crucial importance to understanding or predicting the physiological role of genes and proteins and how they interact to form the complex networks that underlie organ development and function. It is, therefore, crucial to determine on a genome-wide level, the spatio-temporal gene expression profiles at cellular resolution. This information is provided by colorimetric RNA in situ hybridization that can elucidate expression of genes in their native context and does so at cellular resolution. We generated what is to our knowledge the first genome-wide transcriptome atlas by RNA in situ hybridization of an entire mammalian organism, the developing mouse at embryonic day 14.5. This digital transcriptome atlas, the Eurexpress atlas (http://www.eurexpress.org), consists of a searchable database of annotated images that can be interactively viewed. We generated anatomy-based expression profiles for over 18,000 coding genes and over 400 microRNAs. We identified 1,002 tissue-specific genes that are a source of novel tissue-specific markers for 37 different anatomical structures. The quality and the resolution of the data revealed novel molecular domains for several developing structures, such as the telencephalon, a novel organization for the hypothalamus, and insight on the Wnt network involved in renal epithelial differentiation during kidney development. The digital transcriptome atlas is a powerful resource to determine co-expression of genes, to identify cell populations and lineages, and to identify functional associations between genes relevant to development and disease.

Profiling the heterogeneity of colorectal cancer consensus molecular subtypes using spatial transcriptomics.

The consensus molecular subtypes (CMS) of colorectal cancer (CRC) is the most widely-used gene expression-based classification and has contributed to a better understanding of disease heterogeneity and prognosis. Nevertheless, CMS intratumoral heterogeneity restricts its clinical application, stressing the necessity of further characterizing the composition and architecture of CRC. Here, we used Spatial Transcriptomics (ST) in combination with single-cell RNA sequencing (scRNA-seq) to decipher the spatially resolved cellular and molecular composition of CRC. In addition to mapping the intratumoral heterogeneity of CMS and their microenvironment, we identified cell communication events in the tumor-stroma interface of CMS2 carcinomas. This includes tumor growth-inhibiting as well as -activating signals, such as the potential regulation of the ETV4 transcriptional activity by DCN or the PLAU-PLAUR ligand-receptor interaction. Our study illustrates the potential of ST to resolve CRC molecular heterogeneity and thereby help advance personalized therapy.

A Simple, Quick, and Partially Automated Protocol for the Isolation of Single Nuclei from Frozen Mammalian Tissues for Single Nucleus Sequencing.

Single-cell and single-nucleus RNA sequencing have become common laboratory applications due to the wealth of transcriptomic information that they provide. Single nucleus RNA sequencing, particularly, is useful for investigating gene expression in difficult-to-dissociate tissues. Furthermore, this approach is also compatible with frozen (archival) material. Here, we describe a protocol to isolate high-quality single nuclei from frozen mammalian tissues for downstream single nucleus RNA sequencing in a partially-automated manner using commercially available instruments and reagents. Specifically, a robotic dissociator is used to automate and standardize tissue homogenization, followed by an optimized chemical gradient to filter the nuclei. Lastly, we accurately and automatically count the nuclei using an automated fluorescent cell counter. The performance of this protocol is demonstrated on mouse brain, rat kidney, and cynomolgus liver and spleen tissue. This protocol is straightforward, rapid, and readily adaptable to various mammalian tissues without requiring extensive optimization and provides good quality nuclei for downstream single nuclei RNA sequencing.

The BTB and CNC homology 1 (BACH1) target genes are involved in the oxidative stress response and in control of the cell cycle.

The regulation of gene expression in response to environmental signals and metabolic imbalances is a key step in maintaining cellular homeostasis. BTB and CNC homology 1 (BACH1) is a heme-binding transcription factor repressing the transcription from a subset of MAF recognition elements at low intracellular heme levels. Upon heme binding, BACH1 is released from the MAF recognition elements, resulting in increased expression of antioxidant response genes. To systematically address the gene regulatory networks involving BACH1, we combined chromatin immunoprecipitation sequencing analysis of BACH1 target genes in HEK 293 cells with knockdown of BACH1 using three independent types of small interfering RNAs followed by transcriptome profiling using microarrays. The 59 BACH1 target genes identified by chromatin immunoprecipitation sequencing were found highly enriched in genes showing expression changes after BACH1 knockdown, demonstrating the impact of BACH1 repression on transcription. In addition to known and new BACH1 targets involved in heme degradation (HMOX1, FTL, FTH1, ME1, and SLC48A1) and redox regulation (GCLC, GCLM, and SLC7A11), we also discovered BACH1 target genes affecting cell cycle and apoptosis pathways (ITPR2, CALM1, SQSTM1, TFE3, EWSR1, CDK6, BCL2L11, and MAFG) as well as subcellular transport processes (CLSTN1, PSAP, MAPT, and vault RNA). The newly identified impact of BACH1 on genes involved in neurodegenerative processes and proliferation provides an interesting basis for future dissection of BACH1-mediated gene repression in neurodegeneration and virus-induced cancerogenesis.

A simple strand-specific RNA-Seq library preparation protocol combining the Illumina TruSeq RNA and the dUTP methods.

Preserving the original RNA orientation information in RNA-Sequencing (RNA-Seq) experiment is essential to the analysis and understanding of the complexity of mammalian transcriptomes. We describe herein a simple, robust, and time-effective protocol for generating strand-specific RNA-seq libraries suited for the Illumina sequencing platform. We modified the Illumina TruSeq RNA sample preparation by implementing the strand specificity feature using the dUTP method. This protocol uses low amounts of starting material and allows a fast processing within two days. It can be easily implemented and requires only few additional reagents to the original Illumina kit.

Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments.

Alternative splicing, polyadenylation of pre-messenger RNA molecules and differential promoter usage can produce a variety of transcript isoforms whose respective expression levels are regulated in time and space, thus contributing specific biological functions. However, the repertoire of mammalian alternative transcripts and their regulation are still poorly understood. Second-generation sequencing is now opening unprecedented routes to address the analysis of entire transcriptomes. Here, we developed methods that allow the prediction and quantification of alternative isoforms derived solely from exon expression levels in RNA-Seq data. These are based on an explicit statistical model and enable the prediction of alternative isoforms within or between conditions using any known gene annotation, as well as the relative quantification of known transcript structures. Applying these methods to a human RNA-Seq dataset, we validated a significant fraction of the predictions by RT-PCR. Data further showed that these predictions correlated well with information originating from junction reads. A direct comparison with exon arrays indicated improved performances of RNA-Seq over microarrays in the prediction of skipped exons. Altogether, the set of methods presented here comprehensively addresses multiple aspects of alternative isoform analysis. The software is available as an open-source R-package called Solas at http://cmb.molgen.mpg.de/2ndGenerationSequencing/Solas/.

TET2-Driven Clonal Hematopoiesis and Response to Canakinumab: An Exploratory Analysis of the CANTOS Randomized Clinical Trial.

Importance: Clonal hematopoiesis of indeterminate potential (CHIP) is associated with increased risk of atherosclerotic cardiovascular disease, and mouse experiments suggest that CHIP related to Tet2 loss of function in myeloid cells accelerates atherosclerosis via augmented interleukin (IL) 1ß signaling. Objective: To assess whether individuals with CHIP have greater cardiovascular event reduction in response to IL-1ß neutralization in the Canankinumab Anti-inflammatory Thrombosis Outcomes Trial (CANTOS). Design, Setting, and Participants: This randomized clinical trial took place from April 2011 to June 2017 at more than 1000 clinical sites in 39 countries. Targeted deep sequencing of genes previously associated with CHIP in a subset of trial participants using genomic DNA prepared from baseline peripheral blood samples were analyzed. All participants had prior myocardial infarction and elevated high-sensitivity C-reactive protein level above 0.20 mg/dL. Analysis took place between June 2017 and December 2021. Interventions: Canakinumab, an anti-IL-1ß antibody, given at doses of 50, 150, and 300 mg once every 3 months. Main Outcomes and Measures: Major adverse cardiovascular events (MACE). Results: A total of 338 patients (8.6%) were identified in this subset with evidence for clonal hematopoiesis. As expected, the incidence of CHIP increased with age; the mean (SD) age of patients with CHIP was 66.3 (9.2) years and 61.5 (9.6) years in patients without CHIP. Unlike other populations that were not preselected for elevated C-reactive protein, in the CANTOS population variants in TET2 were more common than DNMT3A (119 variants in 103 patients vs 86 variants in 85 patients). Placebo-treated patients with CHIP showed a nonsignificant increase in the rate of MACE compared with patients without CHIP using a Cox proportional hazard model (hazard ratio, 1.32 [95% CI, 0.86-2.04]; P = .21). Exploratory analyses of placebo-treated patients with a somatic variant in either TET2 or DNMT3A (n = 58) showed an equivocal risk for MACE (hazard ratio, 1.65 [95% CI, 0.97-2.80]; P = .06). Patients with CHIP due to somatic variants in TET2 also had reduced risk for MACE while taking canakinumab (hazard ratio, 0.38 [95% CI, 0.15-0.96]) with equivocal difference compared with others (P for interaction = .14). Conclusions and Relevance: These results are consistent with observations of increased risk for cardiovascular events in patients with CHIP and raise the possibility that those with TET2 variants may respond better to canakinumab than those without CHIP. Future studies are required to further substantiate this hypothesis. Trial Registration: ClinicalTrials.gov Identifier: NCT01327846.

Whole genome and exome sequencing reference datasets from a multi-center and cross-platform benchmark study.

With the rapid advancement of sequencing technologies, next generation sequencing (NGS) analysis has been widely applied in cancer genomics research. More recently, NGS has been adopted in clinical oncology to advance personalized medicine. Clinical applications of precision oncology require accurate tests that can distinguish tumor-specific mutations from artifacts introduced during NGS processes or data analysis. Therefore, there is an urgent need to develop best practices in cancer mutation detection using NGS and the need for standard reference data sets for systematically measuring accuracy and reproducibility across platforms and methods. Within the SEQC2 consortium context, we established paired tumor-normal reference samples and generated whole-genome (WGS) and whole-exome sequencing (WES) data using sixteen library protocols, seven sequencing platforms at six different centers. We systematically interrogated somatic mutations in the reference samples to identify factors affecting detection reproducibility and accuracy in cancer genomes. These large cross-platform/site WGS and WES datasets using well-characterized reference samples will represent a powerful resource for benchmarking NGS technologies, bioinformatics pipelines, and for the cancer genomics studies.

Toward best practice in cancer mutation detection with whole-genome and whole-exome sequencing.

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor-normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.