Leidos Biomedical Research Inc. (United States)

Characterization of the molecular function of the human genome and its variation across individuals is essential for identifying the cellular mechanisms that underlie human genetic traits and diseases. The Genotype-Tissue Expression (GTEx) project aims to characterize variation in gene expression levels across individuals and diverse tissues of the human body, many of which are not easily accessible. Here we describe genetic effects on gene expression levels across 44 human tissues. We find that local genetic variation affects gene expression levels for the majority of genes, and we further identify inter-chromosomal genetic effects for 93 genes and 112 loci. On the basis of the identified genetic effects, we characterize patterns of tissue specificity, compare local and distal effects, and evaluate the functional properties of the genetic effects. We also demonstrate that multi-tissue, multi-individual data can be used to identify genes and pathways affected by human disease-associated variation, enabling a mechanistic interpretation of gene regulation and the genetic basis of disease.

Abstract Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature 1 . Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium 2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses 3–15 , enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated—but distinct—DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Abstract Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale 1–3 . Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4–5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter 4 ; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation 5,6 ; analyses timings and patterns of tumour evolution 7 ; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity 8,9 ; and evaluates a range of more-specialized features of cancer genomes 8,10–18 .

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans. Deep whole-genome sequencing of 300 individuals from 142 diverse populations provides insights into key population genetic parameters, shows that all modern human ancestry outside of Africa including in Australasians is consistent with descending from a single founding population, and suggests a higher rate of accumulation of mutations in non-Africans compared to Africans since divergence. Three international collaborations reporting in this issue of Nature describe 787 high-quality genomes from individuals from geographically diverse populations. David Reich and colleagues analysed whole-genome sequences of 300 individuals from 142 populations. Their findings include an accelerated estimated rate of accumulation of mutations in non-Africans compared to Africans since divergence, and that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans but from the same source as that of other non-Africans. Eske Willerlsev and colleagues obtained whole-genome data for 83 Aboriginal Australians and 25 Papuans from the New Guinea Highlands. They estimate that Aboriginal Australians and Papuans diverged from Eurasian populations 51,000–72,000 years ago, following a single out-of-Africa dispersal. Luca Pagani et al. report on a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations. Their analyses support the model by which all non-African populations derive most of their genetic ancestry from a single recent migration out of Africa, although a Papuan contribution suggests a trace of an earlier human expansion.

Scalable, integrative methods to understand mechanisms that link genetic variants with phenotypes are needed. Here we derive a mathematical expression to compute PrediXcan (a gene mapping approach) results using summary data (S-PrediXcan) and show its accuracy and general robustness to misspecified reference sets. We apply this framework to 44 GTEx tissues and 100+ phenotypes from GWAS and meta-analysis studies, creating a growing public catalog of associations that seeks to capture the effects of gene expression variation on human phenotypes. Replication in an independent cohort is shown. Most of the associations are tissue specific, suggesting context specificity of the trait etiology. Colocalized significant associations in unexpected tissues underscore the need for an agnostic scanning of multiple contexts to improve our ability to detect causal regulatory mechanisms. Monogenic disease genes are enriched among significant associations for related traits, suggesting that smaller alterations of these genes may cause a spectrum of milder phenotypes.

Abstract Cancer develops through a process of somatic evolution 1,2 . Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes 3 . Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) 4 , we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Abstract A key mutational process in cancer is structural variation, in which rearrangements delete, amplify or reorder genomic segments that range in size from kilobases to whole chromosomes 1–7 . Here we develop methods to group, classify and describe somatic structural variants, using data from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), which aggregated whole-genome sequencing data from 2,658 cancers across 38 tumour types 8 . Sixteen signatures of structural variation emerged. Deletions have a multimodal size distribution, assort unevenly across tumour types and patients, are enriched in late-replicating regions and correlate with inversions. Tandem duplications also have a multimodal size distribution, but are enriched in early-replicating regions—as are unbalanced translocations. Replication-based mechanisms of rearrangement generate varied chromosomal structures with low-level copy-number gains and frequent inverted rearrangements. One prominent structure consists of 2–7 templates copied from distinct regions of the genome strung together within one locus. Such cycles of templated insertions correlate with tandem duplications, and—in liver cancer—frequently activate the telomerase gene TERT . A wide variety of rearrangement processes are active in cancer, which generate complex configurations of the genome upon which selection can act.

The Genotype-Tissue Expression (GTEx) project, sponsored by the NIH Common Fund, was established to study the correlation between human genetic variation and tissue-specific gene expression in non-diseased individuals. A significant challenge was the collection of high-quality biospecimens for extensive genomic analyses. Here we describe how a successful infrastructure for biospecimen procurement was developed and implemented by multiple research partners to support the prospective collection, annotation, and distribution of blood, tissues, and cell lines for the GTEx project. Other research projects can follow this model and form beneficial partnerships with rapid autopsy and organ procurement organizations to collect high quality biospecimens and associated clinical data for genomic studies. Biospecimens, clinical and genomic data, and Standard Operating Procedures guiding biospecimen collection for the GTEx project are available to the research community.

The persistence of HIV-infected cells in individuals on suppressive combination antiretroviral therapy (cART) presents a major barrier for curing HIV infections. HIV integrates its DNA into many sites in the host genome; we identified 2410 integration sites in peripheral blood lymphocytes of five infected individuals on cART. About 40% of the integrations were in clonally expanded cells. Approximately 50% of the infected cells in one patient were from a single clone, and some clones persisted for many years. There were multiple independent integrations in several genes, including MKL2 and BACH2; many of these integrations were in clonally expanded cells. Our findings show that HIV integration sites can play a critical role in expansion and persistence of HIV-infected cells.

Cellular diversity in tumors is a key factor for therapeutic failures and lethal outcomes of solid malignancies. Here, we determined the single-cell transcriptomic landscape of liver cancer biospecimens from 19 patients. We found varying degrees of heterogeneity in malignant cells within and between tumors and diverse landscapes of tumor microenvironment (TME). Strikingly, tumors with higher transcriptomic diversity were associated with patient's worse overall survival. We found a link between hypoxia-dependent vascular endothelial growth factor expression in tumor diversity and TME polarization. Moreover, T cells from higher heterogeneous tumors showed lower cytolytic activities. Consistent results were found using bulk genomic and transcriptomic profiles of 765 liver tumors. Our results offer insight into the diverse ecosystem of liver cancer and its impact on patient prognosis.

Abstract The discovery of drivers of cancer has traditionally focused on protein-coding genes 1–4 . Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium 5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers 6,7 , raise doubts about others and identify novel candidates, including point mutations in the 5′ region of TP53 , in the 3′ untranslated regions of NFKBIZ and TOB1 , focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Background: Resistance to all first-line antibiotics necessitates the use of less effective or more toxic "reserve" agents. Gram-negative bloodstream infections (GNBSIs) harboring such difficult-to-treat resistance (DTR) may have higher mortality than phenotypes that allow for ≥1 active first-line antibiotic. Methods: The Premier Database was analyzed for inpatients with select GNBSIs. DTR was defined as intermediate/resistant in vitro to all ß-lactam categories, including carbapenems and fluoroquinolones. Prevalence and aminoglycoside resistance of DTR episodes were compared with carbapenem-resistant, extended-spectrum cephalosporin-resistant, and fluoroquinolone-resistant episodes using CDC definitions. Predictors of DTR were identified. The adjusted relative risk (aRR) of mortality was examined for DTR, CDC-defined phenotypes susceptible to ≥1 first-line agent, and graded loss of active categories. Results: Between 2009-2013, 471 (1%) of 45011 GNBSI episodes at 92 (53.2%) of 173 hospitals exhibited DTR, ranging from 0.04% for Escherichia coli to 18.4% for Acinetobacter baumannii. Among patients with DTR, 79% received parenteral aminoglycosides, tigecycline, or colistin/polymyxin-B; resistance to all aminoglycosides occurred in 33%. Predictors of DTR included urban healthcare and higher baseline illness. Crude mortality for GNBSIs with DTR was 43%; aRR was higher for DTR than for carbapenem-resistant (1.2; 95% confidence interval, 1.0-1.4; P = .02), extended-spectrum cephalosporin-resistant (1.2; 1.1-1.4; P = .001), or fluoroquinolone-resistant (1.2; 1.0-1.4; P = .008) infections. The mortality aRR increased 20% per graded loss of active first-line categories, from 3-5 to 1-2 to 0. Conclusion: Nonsusceptibility to first-line antibiotics is associated with decreased survival in GNBSIs. DTR is a simple bedside prognostic measure of treatment-limiting coresistance.

Preventing xenograft rejection is one of the greatest challenges of transplantation medicine. Here, we describe a reproducible, long-term survival of cardiac xenografts from alpha 1-3 galactosyltransferase gene knockout pigs, which express human complement regulatory protein CD46 and human thrombomodulin (GTKO.hCD46.hTBM), that were transplanted into baboons. Our immunomodulatory drug regimen includes induction with anti-thymocyte globulin and αCD20 antibody, followed by maintenance with mycophenolate mofetil and an intensively dosed αCD40 (2C10R4) antibody. Median (298 days) and longest (945 days) graft survival in five consecutive recipients using this regimen is significantly prolonged over our recently established survival benchmarks (180 and 500 days, respectively). Remarkably, the reduction of αCD40 antibody dose on day 100 or after 1 year resulted in recrudescence of anti-pig antibody and graft failure. In conclusion, genetic modifications (GTKO.hCD46.hTBM) combined with the treatment regimen tested here consistently prevent humoral rejection and systemic coagulation pathway dysregulation, sustaining long-term cardiac xenograft survival beyond 900 days.

patient-derived xenotransplantation models. Mechanistically, via selective reduction of immunosuppressive M2-like TAMs, RP-182 improved adaptive and innate antitumor immune responses, including increased cancer cell phagocytosis by reprogrammed TAMs.

Abstract Transcript alterations often result from somatic changes in cancer genomes 1 . Various forms of RNA alterations have been described in cancer, including overexpression 2 , altered splicing 3 and gene fusions 4 ; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) 5 . Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis , of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed ‘bridged’ fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

The bacterial chromosome and bacterial plasmids can be engineered in vivo by homologous recombination using PCR products and synthetic oligonucleotides as substrates. This is possible because bacteriophage-encoded recombination proteins efficiently recombine sequences with homologies as short as 35 to 50 bases. Recombineering allows DNA sequences to be inserted or deleted without regard to location of restriction sites. This unit first describes preparation of electrocompetent cells expressing the recombineering functions and their transformation with dsDNA or ssDNA. It then presents support protocols that describe several two-step selection/counter-selection methods of making genetic alterations without leaving any unwanted changes in the targeted DNA, and a method for retrieving onto a plasmid a genetic marker (cloning by retrieval) from the Escherichia coli chromosome or a co-electroporated DNA fragment. Additional protocols describe methods to screen for unselected mutations, removal of the defective prophage from recombineering strains, and other useful techniques.

Reservoirs of infectious HIV-1 persist despite years of combination antiretroviral therapy and make curing HIV-1 infections a major challenge. Most of the proviral DNA resides in CD4(+)T cells. Some of these CD4(+)T cells are clonally expanded; most of the proviruses are defective. It is not known if any of the clonally expanded cells carry replication-competent proviruses. We report that a highly expanded CD4(+) T-cell clone contains an intact provirus. The highly expanded clone produced infectious virus that was detected as persistent plasma viremia during cART in an HIV-1-infected patient who had squamous cell cancer. Cells containing the intact provirus were widely distributed and significantly enriched in cancer metastases. These results show that clonally expanded CD4(+)T cells can be a reservoir of infectious HIV-1.

To assess the utility of human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) as an in vitro proarrhythmia model, we evaluated the concentration dependence and sources of variability of electrophysiologic responses to 28 drugs linked to low, intermediate, and high torsades de pointes (TdP) risk categories using two commercial cell lines and standardized protocols in a blinded multisite study using multielectrode array or voltage-sensing optical approaches. Logistical and ordinal linear regression models were constructed using drug responses as predictors and TdP risk categories as outcomes. Three of seven predictors (drug-induced arrhythmia-like events and prolongation of repolarization at either maximum tested or maximal clinical exposures) categorized drugs with reasonable accuracy (area under the curve values of receiver operator curves ∼0.8). hiPSC-CM line, test site, and platform had minimal influence on drug categorization. These results demonstrate the utility of hiPSC-CMs to detect drug-induced proarrhythmic effects as part of the evolving Comprehensive In Vitro Proarrhythmia Assay paradigm.

Intratumor molecular heterogeneity of hepatocellular carcinoma is partly attributed to the presence of hepatic cancer stem cells (CSCs). Different CSC populations defined by various cell surface markers may contain different oncogenic drivers, posing a challenge in defining molecularly targeted therapeutics. We combined transcriptomic and functional analyses of hepatocellular carcinoma cells at the single-cell level to assess the degree of CSC heterogeneity. We provide evidence that hepatic CSCs at the single-cell level are phenotypically, functionally, and transcriptomically heterogeneous. We found that different CSC subpopulations contain distinct molecular signatures. Interestingly, distinct genes within different CSC subpopulations are independently associated with hepatocellular carcinoma prognosis, suggesting that a diverse hepatic CSC transcriptome affects intratumor heterogeneity and tumor progression. CONCLUSION: Our work provides unique perspectives into the biodiversity of CSC subpopulations, whose molecular heterogeneity further highlights their role in tumor heterogeneity, prognosis, and hepatic CSC therapy. (Hepatology 2018;68:127-140).

African genomics and skin color Skin color varies among human populations and is thought to be under selection, with light skin maximizing vitamin D production at higher latitudes and dark skin providing UV protection in equatorial zones. To identify the genes that give rise to the palette of human skin tones, Crawford et al. applied genome-wide analyses across diverse African populations (see the Perspective by Tang and Barsh). Genetic variants were identified with likely function in skin phenotypes. Comparison to model organisms verified a conserved function of MFSD12 in pigmentation. A global genetic panel was used to trace how alleles associated with skin color likely moved across the globe as humans migrated, both within and out of Africa. Science , this issue p. eaan8433 ; see also p. 867