Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood.

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

Histologic and Quality Assessment of Genotype-Tissue Expression (GTEx) Research Samples: A Large Postmortem Tissue Collection.

CONTEXT.­: The National Institutes of Health Genotype-Tissue Expression (GTEx) project was developed to elucidate how genetic variation influences gene expression in multiple normal tissues procured from postmortem donors. OBJECTIVE.­: To provide critical insight into a biospecimen's suitability for subsequent analysis, each biospecimen underwent quality assessment measures that included evaluation for underlying disease and potential effects introduced by preanalytic factors. DESIGN.­: Electronic images of each tissue collected from nearly 1000 postmortem donors were evaluated by board-certified pathologists for the extent of autolysis, tissue purity, and the type and abundance of any extraneous tissue. Tissue-specific differences in the severity of autolysis and RNA integrity were evaluated, as were potential relationships between these markers and the duration of postmortem interval and rapidity of death. RESULTS.­: Tissue-specific challenges in the procurement and preservation of the nearly 30 000 tissue specimens collected during the GTEx project are summarized. Differences in the degree of autolysis and RNA integrity number were observed among the 40 tissue types evaluated, and tissue-specific susceptibilities to the duration of postmortem interval and rapidity of death were observed. CONCLUSIONS.­: Ninety-five percent of tissues were of sufficient quality to support RNA sequencing analysis. Biospecimens, annotated whole slide images, de-identified clinical data, and genomic data generated for GTEx represent a high-quality and comprehensive resource for the scientific community that has contributed to its use in approximately 1695 articles. Biospecimens and data collected under the GTEx project are available via the GTEx portal and authorized access to the Database of Genotypes and Phenotypes; procedures and whole slide images are available from the National Cancer Institute.

Peripheral Blood Mononuclear Cell Gene Expression Associated with Pulmonary Microvascular Perfusion: The Multi-Ethnic Study of Atherosclerosis Chronic Obstructive Pulmonary Disease.

Rationale: Chronic obstructive pulmonary disease (COPD) and emphysema are associated with endothelial damage and altered pulmonary microvascular perfusion. The molecular mechanisms underlying these changes are poorly understood in patients, in part because of the inaccessibility of the pulmonary vasculature. Peripheral blood mononuclear cells (PBMCs) interact with the pulmonary endothelium. Objectives: To test the association between gene expression in PBMCs and pulmonary microvascular perfusion in COPD. Methods: The Multi-Ethnic Study of Atherosclerosis (MESA) COPD Study recruited two independent samples of COPD cases and controls with ⩾10 pack-years of smoking history. In both samples, pulmonary microvascular blood flow, pulmonary microvascular blood volume, and mean transit time were assessed on contrast-enhanced magnetic resonance imaging, and PBMC gene expression was assessed by microarray. Additional replication was performed in a third sample with pulmonary microvascular blood volume measures on contrast-enhanced dual-energy computed tomography. Differential expression analyses were adjusted for age, gender, race/ethnicity, educational attainment, height, weight, smoking status, and pack-years of smoking. Results: The 79 participants in the discovery sample had a mean age of 69 ± 6 years, 44% were female, 25% were non-White, 34% were current smokers, and 66% had COPD. There were large PBMC gene expression signatures associated with pulmonary microvascular perfusion traits, with several replicated in the replication sets with magnetic resonance imaging (n = 47) or dual-energy contrast-enhanced computed tomography (n = 157) measures. Many of the identified genes are involved in inflammatory processes, including nuclear factor-κB and chemokine signaling pathways. Conclusions: PBMC gene expression in nuclear factor-κB, inflammatory, and chemokine signaling pathways was associated with pulmonary microvascular perfusion in COPD, potentially offering new targetable candidates for novel therapies.

Distinct COPD subtypes in former smokers revealed by gene network perturbation analysis.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) varies significantly in symptomatic and physiologic presentation. Identifying disease subtypes from molecular data, collected from easily accessible blood samples, can help stratify patients and guide disease management and treatment. METHODS: Blood gene expression measured by RNA-sequencing in the COPDGene Study was analyzed using a network perturbation analysis method. Each COPD sample was compared against a learned reference gene network to determine the part that is deregulated. Gene deregulation values were used to cluster the disease samples. RESULTS: The discovery set included 617 former smokers from COPDGene. Four distinct gene network subtypes are identified with significant differences in symptoms, exercise capacity and mortality. These clusters do not necessarily correspond with the levels of lung function impairment and are independently validated in two external cohorts: 769 former smokers from COPDGene and 431 former smokers in the Multi-Ethnic Study of Atherosclerosis (MESA). Additionally, we identify several genes that are significantly deregulated across these subtypes, including DSP and GSTM1, which have been previously associated with COPD through genome-wide association study (GWAS). CONCLUSIONS: The identified subtypes differ in mortality and in their clinical and functional characteristics, underlining the need for multi-dimensional assessment potentially supplemented by selected markers of gene expression. The subtypes were consistent across cohorts and could be used for new patient stratification and disease prognosis.

Transcriptome variation in human tissues revealed by long-read sequencing.

Regulation of transcript structure generates transcript diversity and plays an important role in human disease1-7. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure8-16. In this Article, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from Genotype-Tissue Expression (GTEx) tissues and cell lines, complementing the GTEx resource. We identified just over 70,000 novel transcripts for annotated genes, and validated the protein expression of 10% of novel transcripts. We developed a new computational package, LORALS, to analyse the genetic effects of rare and common variants on the transcriptome by allele-specific analysis of long reads. We characterized allele-specific expression and transcript structure events, providing new insights into the specific transcript alterations caused by common and rare genetic variants and highlighting the resolution gained from long-read data. We were able to perturb the transcript structure upon knockdown of PTBP1, an RNA binding protein that mediates splicing, thereby finding genetic regulatory effects that are modified by the cellular environment. Finally, we used this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.

Resources Required for Implementation of SARS-CoV-2 Screening in Massachusetts K-12 Public Schools in Winter/Spring 2021.

BACKGROUND: The financial costs and human resource requirements at the school and district level to implement a SARS-CoV-2 screening program are not well known. METHODS: A consortium of Massachusetts public K-12 schools was formed to implement and evaluate a range of SARS-CoV-2 screening approaches. Participating districts were surveyed weekly about their programs, including: type of assay used, individual vs. pooled screening, approaches to return of results and deconvolution of positive pools, number and type of personnel, and hours spent implementing the screening program, and hours spent on program implementation. RESULTS: In 21 participating districts, over 21 weeks from January to June 2021, the positivity rate was 0.0% to 0.21% among students and 0.0% to 0.13% among educators/staff. The average weekly cost to implement a screening program, including assay and personnel costs, was $17.00 per person tested; this was $46.68 for individual screenings and $15.61 for pooled screenings. The total weekly costs by district ranged from $1,644 to $93,486, and districts screened between 58 and 3675 people per week. CONCLUSIONS: Where screening is recommended for the 2021 to 2022 school year due to high COVID-19 incidence, understanding the human resources and finances required to implement screening will assist district policymakers in planning.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes.

Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown1. The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer2-4 and coronary heart disease5-this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP)6. Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

Transcriptional and Cellular Diversity of the Human Heart.

BACKGROUND: The human heart requires a complex ensemble of specialized cell types to perform its essential function. A greater knowledge of the intricate cellular milieu of the heart is critical to increase our understanding of cardiac homeostasis and pathology. As recent advances in low-input RNA sequencing have allowed definitions of cellular transcriptomes at single-cell resolution at scale, we have applied these approaches to assess the cellular and transcriptional diversity of the nonfailing human heart. METHODS: Microfluidic encapsulation and barcoding was used to perform single nuclear RNA sequencing with samples from 7 human donors, selected for their absence of overt cardiac disease. Individual nuclear transcriptomes were then clustered based on transcriptional profiles of highly variable genes. These clusters were used as the basis for between-chamber and between-sex differential gene expression analyses and intersection with genetic and pharmacologic data. RESULTS: We sequenced the transcriptomes of 287 269 single cardiac nuclei, revealing 9 major cell types and 20 subclusters of cell types within the human heart. Cellular subclasses include 2 distinct groups of resident macrophages, 4 endothelial subtypes, and 2 fibroblast subsets. Comparisons of cellular transcriptomes by cardiac chamber or sex reveal diversity not only in cardiomyocyte transcriptional programs but also in subtypes involved in extracellular matrix remodeling and vascularization. Using genetic association data, we identified strong enrichment for the role of cell subtypes in cardiac traits and diseases. Intersection of our data set with genes on cardiac clinical testing panels and the druggable genome reveals striking patterns of cellular specificity. CONCLUSIONS: Using large-scale single nuclei RNA sequencing, we defined the transcriptional and cellular diversity in the normal human heart. Our identification of discrete cell subtypes and differentially expressed genes within the heart will ultimately facilitate the development of new therapeutics for cardiovascular diseases.

RNA sequence analysis reveals macroscopic somatic clonal expansion across normal tissues.

How somatic mutations accumulate in normal cells is poorly understood. A comprehensive analysis of RNA sequencing data from ~6700 samples across 29 normal tissues revealed multiple somatic variants, demonstrating that macroscopic clones can be found in many normal tissues. We found that sun-exposed skin, esophagus, and lung have a higher mutation burden than other tested tissues, which suggests that environmental factors can promote somatic mosaicism. Mutation burden was associated with both age and tissue-specific cell proliferation rate, highlighting that mutations accumulate over both time and number of cell divisions. Finally, normal tissues were found to harbor mutations in known cancer genes and hotspots. This study provides a broad view of macroscopic clonal expansion in human tissues, thus serving as a foundation for associating clonal expansion with environmental factors, aging, and risk of disease.

A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles.

We previously piloted the concept of a Connectivity Map (CMap), whereby genes, drugs, and disease states are connected by virtue of common gene-expression signatures. Here, we report more than a 1,000-fold scale-up of the CMap as part of the NIH LINCS Consortium, made possible by a new, low-cost, high-throughput reduced representation expression profiling method that we term L1000. We show that L1000 is highly reproducible, comparable to RNA sequencing, and suitable for computational inference of the expression levels of 81% of non-measured transcripts. We further show that the expanded CMap can be used to discover mechanism of action of small molecules, functionally annotate genetic variants of disease genes, and inform clinical trials. The 1.3 million L1000 profiles described here, as well as tools for their analysis, are available at https://clue.io.

Punctuated evolution of prostate cancer genomes.

The analysis of exonic DNA from prostate cancers has identified recurrently mutated genes, but the spectrum of genome-wide alterations has not been profiled extensively in this disease. We sequenced the genomes of 57 prostate tumors and matched normal tissues to characterize somatic alterations and to study how they accumulate during oncogenesis and progression. By modeling the genesis of genomic rearrangements, we identified abundant DNA translocations and deletions that arise in a highly interdependent manner. This phenomenon, which we term "chromoplexy," frequently accounts for the dysregulation of prostate cancer genes and appears to disrupt multiple cancer genes coordinately. Our modeling suggests that chromoplexy may induce considerable genomic derangement over relatively few events in prostate cancer and other neoplasms, supporting a model of punctuated cancer evolution. By characterizing the clonal hierarchy of genomic lesions in prostate tumors, we charted a path of oncogenic events along which chromoplexy may drive prostate carcinogenesis.

A landscape of driver mutations in melanoma.

Despite recent insights into melanoma genetics, systematic surveys for driver mutations are challenged by an abundance of passenger mutations caused by carcinogenic UV light exposure. We developed a permutation-based framework to address this challenge, employing mutation data from intronic sequences to control for passenger mutational load on a per gene basis. Analysis of large-scale melanoma exome data by this approach discovered six novel melanoma genes (PPP6C, RAC1, SNX31, TACC1, STK19, and ARID2), three of which-RAC1, PPP6C, and STK19-harbored recurrent and potentially targetable mutations. Integration with chromosomal copy number data contextualized the landscape of driver mutations, providing oncogenic insights in BRAF- and NRAS-driven melanoma as well as those without known NRAS/BRAF mutations. The landscape also clarified a mutational basis for RB and p53 pathway deregulation in this malignancy. Finally, the spectrum of driver mutations provided unequivocal genomic evidence for a direct mutagenic role of UV light in melanoma pathogenesis.

Sequence analysis of mutations and translocations across breast cancer subtypes.

Breast carcinoma is the leading cause of cancer-related mortality in women worldwide, with an estimated 1.38 million new cases and 458,000 deaths in 2008 alone. This malignancy represents a heterogeneous group of tumours with characteristic molecular features, prognosis and responses to available therapy. Recurrent somatic alterations in breast cancer have been described, including mutations and copy number alterations, notably ERBB2 amplifications, the first successful therapy target defined by a genomic aberration. Previous DNA sequencing studies of breast cancer genomes have revealed additional candidate mutations and gene rearrangements. Here we report the whole-exome sequences of DNA from 103 human breast cancers of diverse subtypes from patients in Mexico and Vietnam compared to matched-normal DNA, together with whole-genome sequences of 22 breast cancer/normal pairs. Beyond confirming recurrent somatic mutations in PIK3CA, TP53, AKT1, GATA3 and MAP3K1, we discovered recurrent mutations in the CBFB transcription factor gene and deletions of its partner RUNX1. Furthermore, we have identified a recurrent MAGI3-AKT3 fusion enriched in triple-negative breast cancer lacking oestrogen and progesterone receptors and ERBB2 expression. The MAGI3-AKT3 fusion leads to constitutive activation of AKT kinase, which is abolished by treatment with an ATP-competitive AKT small-molecule inhibitor.

Melanoma genome sequencing reveals frequent PREX2 mutations.

Melanoma is notable for its metastatic propensity, lethality in the advanced setting and association with ultraviolet exposure early in life. To obtain a comprehensive genomic view of melanoma in humans, we sequenced the genomes of 25 metastatic melanomas and matched germline DNA. A wide range of point mutation rates was observed: lowest in melanomas whose primaries arose on non-ultraviolet-exposed hairless skin of the extremities (3 and 14 per megabase (Mb) of genome), intermediate in those originating from hair-bearing skin of the trunk (5-55 per Mb), and highest in a patient with a documented history of chronic sun exposure (111 per Mb). Analysis of whole-genome sequence data identified PREX2 (phosphatidylinositol-3,4,5-trisphosphate-dependent Rac exchange factor 2)--a PTEN-interacting protein and negative regulator of PTEN in breast cancer--as a significantly mutated gene with a mutation frequency of approximately 14% in an independent extension cohort of 107 human melanomas. PREX2 mutations are biologically relevant, as ectopic expression of mutant PREX2 accelerated tumour formation of immortalized human melanocytes in vivo. Thus, whole-genome sequencing of human melanoma tumours revealed genomic evidence of ultraviolet pathogenesis and discovered a new recurrently mutated gene in melanoma.

Exome sequencing identifies recurrent SPOP, FOXA1 and MED12 mutations in prostate cancer.

Prostate cancer is the second most common cancer in men worldwide and causes over 250,000 deaths each year. Overtreatment of indolent disease also results in significant morbidity. Common genetic alterations in prostate cancer include losses of NKX3.1 (8p21) and PTEN (10q23), gains of AR (the androgen receptor gene) and fusion of ETS family transcription factor genes with androgen-responsive promoters. Recurrent somatic base-pair substitutions are believed to be less contributory in prostate tumorigenesis but have not been systematically analyzed in large cohorts. Here, we sequenced the exomes of 112 prostate tumor and normal tissue pairs. New recurrent mutations were identified in multiple genes, including MED12 and FOXA1. SPOP was the most frequently mutated gene, with mutations involving the SPOP substrate-binding cleft in 6-15% of tumors across multiple independent cohorts. Prostate cancers with mutant SPOP lacked ETS family gene rearrangements and showed a distinct pattern of genomic alterations. Thus, SPOP mutations may define a new molecular subtype of prostate cancer.

The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity.

The systematic translation of cancer genomic data into knowledge of tumour biology and therapeutic possibilities remains challenging. Such efforts should be greatly aided by robust preclinical model systems that reflect the genomic diversity of human cancers and for which detailed genetic and pharmacological annotation is available. Here we describe the Cancer Cell Line Encyclopedia (CCLE): a compilation of gene expression, chromosomal copy number and massively parallel sequencing data from 947 human cancer cell lines. When coupled with pharmacological profiles for 24 anticancer drugs across 479 of the cell lines, this collection allowed identification of genetic, lineage, and gene-expression-based predictors of drug sensitivity. In addition to known predictors, we found that plasma cell lineage correlated with sensitivity to IGF1 receptor inhibitors; AHR expression was associated with MEK inhibitor efficacy in NRAS-mutant lines; and SLFN11 expression predicted sensitivity to topoisomerase inhibitors. Together, our results indicate that large, annotated cell-line collections may help to enable preclinical stratification schemata for anticancer agents. The generation of genetic predictions of drug response in the preclinical setting and their incorporation into cancer clinical trial design could speed the emergence of 'personalized' therapeutic regimens.

RAS mutations are associated with the development of cutaneous squamous cell tumors in patients treated with RAF inhibitors.

PURPOSE: RAF inhibitors are effective against melanomas with BRAF V600E mutations but may induce keratoacanthomas (KAs) and cutaneous squamous cell carcinomas (cSCCs). The potential of these agents to promote secondary malignancies is concerning. We analyzed cSCC and KA lesions for genetic mutations in an attempt to identify an underlying mechanism for their formation. METHODS: Four international centers contributed 237 KA or cSCC tumor samples from patients receiving an RAF inhibitor (either vemurafenib or sorafenib; n = 19) or immunosuppression therapy (n = 53) or tumors that developed spontaneously (n = 165). Each sample was profiled for 396 known somatic mutations across 33 cancer-related genes by using a mass spectrometric-based genotyping platform. RESULTS: Mutations were detected in 16% of tumors (38 of 237), with five tumors harboring two mutations. Mutations in TP53, CDKN2A, HRAS, KRAS, and PIK3CA were previously described in squamous cell tumors. Mutations in MYC, FGFR3, and VHL were identified for the first time. A higher frequency of activating RAS mutations was found in tumors from patients treated with an RAF inhibitor versus populations treated with a non-RAF inhibitor (21.1% v 3.2%; P < .01), although overall mutation rates between treatment groups were similar (RAF inhibitor, 21.1%; immunosuppression, 18.9%; and spontaneous, 17.6%; P = not significant). Tumor histology (KA v cSCC), tumor site (head and neck v other), patient age (≤ 70 v > 70 years), and sex had no significant impact on mutation rate or type. CONCLUSION: Squamous cell tumors from patients treated with an RAF inhibitor have a distinct mutational profile that supports a mechanism of therapy-induced tumorigenesis in RAS-primed cells. Conceivably, cotargeting of MEK together with RAF may reduce or prevent formation of these tumors.

The mutational landscape of head and neck squamous cell carcinoma.

Head and neck squamous cell carcinoma (HNSCC) is a common, morbid, and frequently lethal malignancy. To uncover its mutational spectrum, we analyzed whole-exome sequencing data from 74 tumor-normal pairs. The majority exhibited a mutational profile consistent with tobacco exposure; human papillomavirus was detectable by sequencing DNA from infected tumors. In addition to identifying previously known HNSCC genes (TP53, CDKN2A, PTEN, PIK3CA, and HRAS), our analysis revealed many genes not previously implicated in this malignancy. At least 30% of cases harbored mutations in genes that regulate squamous differentiation (for example, NOTCH1, IRF6, and TP63), implicating its dysregulation as a major driver of HNSCC carcinogenesis. More generally, the results indicate the ability of large-scale sequencing to reveal fundamental tumorigenic mechanisms.

Initial genome sequencing and analysis of multiple myeloma.

Multiple myeloma is an incurable malignancy of plasma cells, and its pathogenesis is poorly understood. Here we report the massively parallel sequencing of 38 tumour genomes and their comparison to matched normal DNAs. Several new and unexpected oncogenic mechanisms were suggested by the pattern of somatic mutation across the data set. These include the mutation of genes involved in protein translation (seen in nearly half of the patients), genes involved in histone methylation, and genes involved in blood coagulation. In addition, a broader than anticipated role of NF-κB signalling was indicated by mutations in 11 members of the NF-κB pathway. Of potential immediate clinical relevance, activating mutations of the kinase BRAF were observed in 4% of patients, suggesting the evaluation of BRAF inhibitors in multiple myeloma clinical trials. These results indicate that cancer genome sequencing of large collections of samples will yield new insights into cancer not anticipated by existing knowledge.