Adenosine-to-inosine (A-to-I) RNA editing is a conserved post-transcriptional mechanism mediated by ADAR enzymes that diversifies the transcriptome by altering selected nucleotides in RNA molecules1. Although many editing sites have recently been discovered2–7, the extent to which most sites are edited and how the editing is regulated in different biological contexts are not fully understood8–10. Here we report dynamic spatiotemporal patterns and new regulators of RNA editing, discovered through an extensive profiling of A-to-I RNA editing in 8,551 human samples (representing 53 body sites from 552 individuals) from the Genotype-Tissue Expression (GTEx) project and in hundreds of other primate and mouse samples. We show that editing levels in non-repetitive coding regions vary more between tissues than editing levels in repetitive regions. Globally, ADAR1 is the primary editor of repetitive sites and ADAR2 is the primary editor of non-repetitive coding sites, whereas the catalytically inactive ADAR3 predominantly acts as an inhibitor of editing. Cross-species analysis of RNA editing in several tissues revealed that species, rather than tissue type, is the primary determinant of editing levels, suggesting stronger cis-directed regulation of RNA editing for most sites, although the small set of conserved coding sites is under stronger trans-regulation. In addition, we curated an extensive set of ADAR1 and ADAR2 targets and showed that many editing sites display distinct tissue-specific regulation by the ADAR enzymes in vivo. Further analysis of the GTEx data revealed several potential regulators of editing, such as AIMP2, which reduces editing in muscles by enhancing the degradation of the ADAR proteins. Collectively, our work provides insights into the complex cis- and trans-regulation of A-to-I editing.

CIViC is an expert-crowdsourced knowledgebase for Clinical Interpretation of Variants in Cancer describing the therapeutic, prognostic, diagnostic and predisposing relevance of inherited and somatic variants of all types. CIViC is committed to open-source code, open-access content, public application programming interfaces (APIs) and provenance of supporting evidence to allow for the transparent creation of current and accurate variant interpretations for use in cancer precision medicine.

Summary The feasibility of performing broad and deep tumour genome sequencing has shed new light into tumour heterogeneity and provided important insights into the evolution of metastases arising from different clones1,2. To add an additional layer of complexity, tumour evolution may be influenced by selective pressure provided by therapy, in a similar fashion as it occurs in infectious diseases. Here, we have studied the tumour genomic evolution in a patient with metastatic breast cancer bearing an activating PIK3CA mutation. The patient was treated with the PI3Kα inhibitor BYL719 and achieved a lasting clinical response, although eventually progressed to treatment and died shortly thereafter. A rapid autopsy was performed and a total of 14 metastatic sites were collected and sequenced. All metastatic lesions, when compared to the pre-treatment tumour, had a copy loss of PTEN, and those lesions that became refractory to BYL719 had additional and different PTEN genetic alterations, resulting in the loss of PTEN expression. Acquired bi-allelic loss of PTEN was found in one additional patient treated with BYL719 whereas in two patients PIK3CA mutations present in the primary tumour were no longer detected at the time of progression. To functionally characterize our findings, inducible PTEN knockdown in sensitive cells resulted in resistance to BYL719, while simultaneous PI3Kp110β blockade reverted this resistance phenotype, both in cell lines and in PTEN-null xenografts derived from our patient. We conclude that parallel genetic evolution of separate sites with different PTEN genomic alterations leads to a convergent PTEN- null phenotype resistant to PI3Kα inhibition.

Summary Tumors are typically sequenced to depths of 75–100× (exome) or 30–50× (whole genome). We demonstrate that current sequencing paradigms are inadequate for tumors that are impure, aneuploid or clonally heterogeneous. To reassess optimal sequencing strategies, we performed ultra-deep (up to ~312×) whole genome sequencing (WGS) and exome capture (up to ~433×) of a primary acute myeloid leukemia, its subsequent relapse, and a matched normal skin sample. We tested multiple alignment and variant calling algorithms and validated ~200,000 putative SNVs by sequencing them to depths of ~1,000×. Additional targeted sequencing provided over 10,000× coverage and ddPCR assays provided up to ~250,000× sampling of selected sites. We evaluated the effects of different library generation approaches, depth of sequencing, and analysis strategies on the ability to effectively characterize a complex tumor. This dataset, representing the most comprehensively sequenced tumor described to date, will serve as an invaluable community resource (dbGaP accession id phs000159).

We present the DeNovoGear software for analyzing de novo mutations from familial and somatic tissue sequencing data. DeNovoGear uses likelihood-based error modeling to reduce the false positive rate of mutation discovery in exome analysis, and fragment information to identify the parental origin of germline mutations. We used our program to create a whole-genome de novo indel callset with a 95% validation rate, producing a direct estimate of the human germline indel mutation rate.