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.

Gonadal failure, along with early pregnancy loss and perinatal death, may be an important filter that limits the propagation of harmful mutations in the human population. We hypothesized that men with spermatogenic impairment, a disease with unknown genetic architecture and a common cause of male infertility, are enriched for rare deleterious mutations compared to men with normal spermatogenesis. After assaying genomewide SNPs and CNVs in 323 Caucasian men with idiopathic spermatogenic impairment and more than 1,100 controls, we estimate that each rare autosomal deletion detected in our study multiplicatively changes a man's risk of disease by 10% (OR 1.10 [1.04–1.16], p<2×10−3), rare X-linked CNVs by 29%, (OR 1.29 [1.11–1.50], p<1×10−3), and rare Y-linked duplications by 88% (OR 1.88 [1.13–3.13], p<0.03). By contrasting the properties of our case-specific CNVs with those of CNV callsets from cases of autism, schizophrenia, bipolar disorder, and intellectual disability, we propose that the CNV burden in spermatogenic impairment is distinct from the burden of large, dominant mutations described for neurodevelopmental disorders. We identified two patients with deletions of DMRT1, a gene on chromosome 9p24.3 orthologous to the putative sex determination locus of the avian ZW chromosome system. In an independent sample of Han Chinese men, we identified 3 more DMRT1 deletions in 979 cases of idiopathic azoospermia and none in 1,734 controls, and found none in an additional 4,519 controls from public databases. The combined results indicate that DMRT1 loss-of-function mutations are a risk factor and potential genetic cause of human spermatogenic failure (frequency of 0.38% in 1306 cases and 0% in 7,754 controls, p = 6.2×10−5). Our study identifies other recurrent CNVs as potential causes of idiopathic azoospermia and generates hypotheses for directing future studies on the genetic basis of male infertility and IVF outcomes.

Nearly all patients with small cell lung cancer (SCLC) eventually relapse with chemoresistant disease. The molecular mechanisms driving chemoresistance in SCLC remain un-characterized. Here, we describe whole-exome sequencing of paired SCLC tumor samples procured at diagnosis and relapse from 12 patients, and unpaired relapse samples from 18 additional patients. Multiple somatic copy number alterations, including gains in ABCC1 and deletions in MYCL, MSH2, and MSH6, are identifiable in relapsed samples. Relapse samples also exhibit recurrent mutations and loss of heterozygosity in regulators of WNT signaling, including CHD8 and APC. Analysis of RNA-sequencing data shows enrichment for an ASCL1-low expression subtype and WNT activation in relapse samples. Activation of WNT signaling in chemosensitive human SCLC cell lines through APC knockdown induces chemoresistance. Additionally, in vitro-derived chemoresistant cell lines demonstrate increased WNT activity. Overall, our results suggest WNT signaling activation as a mechanism of chemoresistance in relapsed SCLC.

In this work, we present the Genome Modeling System (GMS), an analysis information management system capable of executing automated genome analysis pipelines at a massive scale. The GMS framework provides detailed tracking of samples and data coupled with reliable and repeatable analysis pipelines. The GMS also serves as a platform for bioinformatics development, allowing a large team to collaborate on data analysis, or an individual researcher to leverage the work of others effectively within its data management system. Rather than separating ad-hoc analysis from rigorous, reproducible pipelines, the GMS promotes systematic integration between the two. As a demonstration of the GMS, we performed an integrated analysis of whole genome, exome and transcriptome sequencing data from a breast cancer cell line (HCC1395) and matched lymphoblastoid line (HCC1395BL). These data are available for users to test the software, complete tutorials and develop novel GMS pipeline configurations. The GMS is available at https://github.com/genome/gms.