Whole-genome doubling confers unique genetic vulnerabilities on tumour cells.

Whole-genome doubling (WGD) is common in human cancers, occurring early in tumorigenesis and generating genetically unstable tetraploid cells that fuel tumour development1,2. Cells that undergo WGD (WGD+ cells) must adapt to accommodate their abnormal tetraploid state; however, the nature of these adaptations, and whether they confer vulnerabilities that can be exploited therapeutically, is unclear. Here, using sequencing data from roughly 10,000 primary human cancer samples and essentiality data from approximately 600 cancer cell lines, we show that WGD gives rise to common genetic traits that are accompanied by unique vulnerabilities. We reveal that WGD+ cells are more dependent than WGD- cells on signalling from the spindle-assembly checkpoint, DNA-replication factors and proteasome function. We also identify KIF18A, which encodes a mitotic kinesin protein, as being specifically required for the viability of WGD+ cells. Although KIF18A is largely dispensable for accurate chromosome segregation during mitosis in WGD- cells, its loss induces notable mitotic errors in WGD+ cells, ultimately impairing cell viability. Collectively, our results suggest new strategies for specifically targeting WGD+ cancer cells while sparing the normal, non-transformed WGD- cells that comprise human tissue.

Characterization and decontamination of background noise in droplet-based single-cell protein expression data with DecontPro.

Assays such as CITE-seq can measure the abundance of cell surface proteins on individual cells using antibody derived tags (ADTs). However, many ADTs have high levels of background noise that can obfuscate down-stream analyses. In an exploratory analysis of PBMC datasets, we find that some droplets that were originally called 'empty' due to low levels of RNA contained high levels of ADTs and likely corresponded to neutrophils. We identified a novel type of artifact in the empty droplets called a 'spongelet' which has medium levels of ADT expression and is distinct from ambient noise. ADT expression levels in the spongelets correlate to ADT expression levels in the background peak of true cells in several datasets suggesting that they can contribute to background noise along with ambient ADTs. We then developed DecontPro, a novel Bayesian hierarchical model that can decontaminate ADT data by estimating and removing contamination from these sources. DecontPro outperforms other decontamination tools in removing aberrantly expressed ADTs while retaining native ADTs and in improving clustering specificity. Overall, these results suggest that identification of empty drops should be performed separately for RNA and ADT data and that DecontPro can be incorporated into CITE-seq workflows to improve the quality of downstream analyses.

Bulk brain tissue cell-type deconvolution with bias correction for single-nuclei RNA sequencing data using DeTREM.

BACKGROUND: Quantifying cell-type abundance in bulk tissue RNA-sequencing enables researchers to better understand complex systems. Newer deconvolution methodologies, such as MuSiC, use cell-type signatures derived from single-cell RNA-sequencing (scRNA-seq) data to make these calculations. Single-nuclei RNA-sequencing (snRNA-seq) reference data can be used instead of scRNA-seq data for tissues such as human brain where single-cell data are difficult to obtain, but accuracy suffers due to sequencing differences between the technologies. RESULTS: We propose a modification to MuSiC entitled 'DeTREM' which compensates for sequencing differences between the cell-type signature and bulk RNA-seq datasets in order to better predict cell-type fractions. We show DeTREM to be more accurate than MuSiC in simulated and real human brain bulk RNA-sequencing datasets with various cell-type abundance estimates. We also compare DeTREM to SCDC and CIBERSORTx, two recent deconvolution methods that use scRNA-seq cell-type signatures. We find that they perform well in simulated data but produce less accurate results than DeTREM when used to deconvolute human brain data. CONCLUSION: DeTREM improves the deconvolution accuracy of MuSiC and outperforms other deconvolution methods when applied to snRNA-seq data. DeTREM enables accurate cell-type deconvolution in situations where scRNA-seq data are not available. This modification improves characterization cell-type specific effects in brain tissue and identification of cell-type abundance differences under various conditions.

ExperimentSubset: an R package to manage subsets of Bioconductor Experiment objects.

MOTIVATION: R Experiment objects such as the SummarizedExperiment or SingleCellExperiment are data containers for storing one or more matrix-like assays along with associated row and column data. These objects have been used to facilitate the storage and analysis of high-throughput genomic data generated from technologies such as single-cell RNA sequencing. One common computational task in many genomics analysis workflows is to perform subsetting of the data matrix before applying down-stream analytical methods. For example, one may need to subset the columns of the assay matrix to exclude poor-quality samples or subset the rows of the matrix to select the most variable features. Traditionally, a second object is created that contains the desired subset of assay from the original object. However, this approach is inefficient as it requires the creation of an additional object containing a copy of the original assay and leads to challenges with data provenance. RESULTS: To overcome these challenges, we developed an R package called ExperimentSubset, which is a data container that implements classes for efficient storage and streamlined retrieval of assays that have been subsetted by rows and/or columns. These classes are able to inherently provide data provenance by maintaining the relationship between the subsetted and parent assays. We demonstrate the utility of this package on a single-cell RNA-seq dataset by storing and retrieving subsets at different stages of the analysis while maintaining a lower memory footprint. Overall, the ExperimentSubset is a flexible container for the efficient management of subsets. AVAILABILITY AND IMPLEMENTATION: ExperimentSubset package is available at Bioconductor: https://bioconductor.org/packages/ExperimentSubset/ and Github: https://github.com/campbio/ExperimentSubset. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

NetSig: network-based discovery from cancer genomes.

Methods that integrate molecular network information and tumor genome data could complement gene-based statistical tests to identify likely new cancer genes; but such approaches are challenging to validate at scale, and their predictive value remains unclear. We developed a robust statistic (NetSig) that integrates protein interaction networks with data from 4,742 tumor exomes. NetSig can accurately classify known driver genes in 60% of tested tumor types and predicts 62 new driver candidates. Using a quantitative experimental framework to determine in vivo tumorigenic potential in mice, we found that NetSig candidates induce tumors at rates that are comparable to those of known oncogenes and are ten-fold higher than those of random genes. By reanalyzing nine tumor-inducing NetSig candidates in 242 patients with oncogene-negative lung adenocarcinomas, we find that two (AKT2 and TFDP2) are significantly amplified. Our study presents a scalable integrated computational and experimental workflow to expand discovery from cancer genomes.

scruff: an R/Bioconductor package for preprocessing single-cell RNA-sequencing data.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) enables the high-throughput quantification of transcriptional profiles in single cells. In contrast to bulk RNA-seq, additional preprocessing steps such as cell barcode identification or unique molecular identifier (UMI) deconvolution are necessary for preprocessing of data from single cell protocols. R packages that can easily preprocess data and rapidly visualize quality metrics and read alignments for individual cells across multiple samples or runs are still lacking. RESULTS: Here we present scruff, an R/Bioconductor package that preprocesses data generated from the CEL-Seq or CEL-Seq2 protocols and reports comprehensive data quality metrics and visualizations. scruff rapidly demultiplexes, aligns, and counts the reads mapped to genome features with deduplication of unique molecular identifier (UMI) tags. scruff also provides novel and extensive functions to visualize both pre- and post-alignment data quality metrics for cells from multiple experiments. Detailed read alignments with corresponding UMI information can be visualized at specific genome coordinates to display differences in isoform usage. The package also supports the visualization of quality metrics for sequence alignment files for multiple experiments generated by Cell Ranger from 10X Genomics. scruff is available as a free and open-source R/Bioconductor package. CONCLUSIONS: scruff streamlines the preprocessing of scRNA-seq data in a few simple R commands. It performs data demultiplexing, alignment, counting, quality report and visualization systematically and comprehensively, ensuring reproducible and reliable analysis of scRNA-seq data.

Genomic approaches to accelerate cancer interception.

Although major advances have been reported in the last decade in the treatment of late-stage cancer with targeted and immune-based therapies, there is a crucial unmet need to develop new approaches to improve the prevention and early detection of cancer. Advances in genomics and computational biology offer unprecedented opportunities to understand the earliest molecular events associated with carcinogenesis, enabling novel strategies to intercept the development of invasive cancers. This Series paper will highlight emerging big data genomic approaches with the potential to accelerate advances in cancer prevention, screening, and early detection across various tumour types, and the challenges inherent in the development of these tools for clinical use. Through coordinated multicentre consortia, these genomic approaches are likely to transform the landscape of cancer interception in the coming years.

Assessment of microRNA differential expression and detection in multiplexed small RNA sequencing data.

Small RNA sequencing can be used to gain an unprecedented amount of detail into the microRNA transcriptome. The relatively high cost and low throughput of sequencing bases technologies can potentially be offset by the use of multiplexing. However, multiplexing involves a trade-off between increased number of sequenced samples and reduced number of reads per sample (i.e., lower depth of coverage). To assess the effect of different sequencing depths owing to multiplexing on microRNA differential expression and detection, we sequenced the small RNA of lung tissue samples collected in a clinical setting by multiplexing one, three, six, nine, or 12 samples per lane using the Illumina HiSeq 2000. As expected, the numbers of reads obtained per sample decreased as the number of samples in a multiplex increased. Furthermore, after normalization, replicate samples included in distinct multiplexes were highly correlated (R > 0.97). When detecting differential microRNA expression between groups of samples, microRNAs with average expression >1 reads per million (RPM) had reproducible fold change estimates (signal to noise) independent of the degree of multiplexing. The number of microRNAs detected was strongly correlated with the log2 number of reads aligning to microRNA loci (R = 0.96). However, most additional microRNAs detected in samples with greater sequencing depth were in the range of expression which had lower fold change reproducibility. These findings elucidate the trade-off between increasing the number of samples in a multiplex with decreasing sequencing depth and will aid in the design of large-scale clinical studies exploring microRNA expression and its role in disease.

SubmiRine: assessing variants in microRNA targets using clinical genomic data sets.

MicroRNAs (miRNAs) regulate gene expression by binding to partially complementary sequences on target mRNA transcripts, thereby causing their degradation, deadenylation, or inhibiting their translation. Genomic variants can alter miRNA regulation by modifying miRNA target sites, and multiple human disease phenotypes have been linked to such miRNA target site variants (miR-TSVs). However, systematic genome-wide identification of functional miR-TSVs is difficult due to high false positive rates; functional miRNA recognition sequences can be as short as six nucleotides, with the human genome encoding thousands of miRNAs. Furthermore, while large-scale clinical genomic data sets are becoming increasingly commonplace, existing miR-TSV prediction methods are not designed to analyze these data. Here, we present an open-source tool called SubmiRine that is designed to perform efficient miR-TSV prediction systematically on variants identified in novel clinical genomic data sets. Most importantly, SubmiRine allows for the prioritization of predicted miR-TSVs according to their relative probability of being functional. We present the results of SubmiRine using integrated clinical genomic data from a large-scale cohort study on chronic obstructive pulmonary disease (COPD), making a number of high-scoring, novel miR-TSV predictions. We also demonstrate SubmiRine's ability to predict and prioritize known miR-TSVs that have undergone experimental validation in previous studies.

MicroRNA 4423 is a primate-specific regulator of airway epithelial cell differentiation and lung carcinogenesis.

Smoking is a significant risk factor for lung cancer, the leading cause of cancer-related deaths worldwide. Although microRNAs are regulators of many airway gene-expression changes induced by smoking, their role in modulating changes associated with lung cancer in these cells remains unknown. Here, we use next-generation sequencing of small RNAs in the airway to identify microRNA 4423 (miR-4423) as a primate-specific microRNA associated with lung cancer and expressed primarily in mucociliary epithelium. The endogenous expression of miR-4423 increases as bronchial epithelial cells undergo differentiation into mucociliary epithelium in vitro, and its overexpression during this process causes an increase in the number of ciliated cells. Furthermore, expression of miR-4423 is reduced in most lung tumors and in cytologically normal epithelium of the mainstem bronchus of smokers with lung cancer. In addition, ectopic expression of miR-4423 in a subset of lung cancer cell lines reduces their anchorage-independent growth and significantly decreases the size of the tumors formed in a mouse xenograft model. Consistent with these phenotypes, overexpression of miR-4423 induces a differentiated-like pattern of airway epithelium gene expression and reverses the expression of many genes that are altered in lung cancer. Together, our results indicate that miR-4423 is a regulator of airway epithelium differentiation and that the abrogation of its function contributes to lung carcinogenesis.

Transcription factor and microRNA interactions in lung cells: an inhibitory link between NK2 homeobox 1, miR-200c and the developmental and oncogenic factors Nfib and Myb.

BACKGROUND: The transcription factor NK2 homeobox 1 (Nkx2-1) plays essential roles in epithelial cell proliferation and differentiation in mouse and human lung development and tumorigenesis. A better understanding of genes and pathways downstream of Nkx2-1 will clarify the multiple roles of this critical lung factor. Nkx2-1 regulates directly or indirectly numerous protein-coding genes; however, there is a paucity of information about Nkx2-1-regulated microRNAs (miRNAs). METHODS AND RESULTS: By miRNA array analyses of mouse epithelial cell lines in which endogenous Nkx2-1 was knocked-down, we revealed that 29 miRNAs were negatively regulated including miR-200c, and 39 miRNAs were positively regulated by Nkx2-1 including miR-1195. Mouse lungs lacking functional phosphorylated Nkx2-1 showed increased expression of miR-200c and alterations in the expression of other top regulated miRNAs. Moreover, chromatin immunoprecipitation assays showed binding of NKX2-1 protein to regulatory regions of these miRNAs. Promoter reporter assays indicated that 1kb of the miR-200c 5' flanking region was transcriptionally active but did not mediate Nkx2-1- repression of miR-200c expression. 3'UTR reporter assays support a direct regulation of the predicted targets Nfib and Myb by miR-200c. CONCLUSIONS: These studies suggest that Nkx2-1 controls the expression of specific miRNAs in lung epithelial cells. In particular, we identified a regulatory link between Nkx2-1, the known tumor suppressor miR-200c, and the developmental and oncogenic transcription factors Nfib and Myb, adding new players to the regulatory mechanisms driven by Nkx2-1 in lung epithelial cells that may have implications in lung development and tumorigenesis.

MODMatcher: multi-omics data matcher for integrative genomic analysis.

Errors in sample annotation or labeling often occur in large-scale genetic or genomic studies and are difficult to avoid completely during data generation and management. For integrative genomic studies, it is critical to identify and correct these errors. Different types of genetic and genomic data are inter-connected by cis-regulations. On that basis, we developed a computational approach, Multi-Omics Data Matcher (MODMatcher), to identify and correct sample labeling errors in multiple types of molecular data, which can be used in further integrative analysis. Our results indicate that inspection of sample annotation and labeling error is an indispensable data quality assurance step. Applied to a large lung genomic study, MODMatcher increased statistically significant genetic associations and genomic correlations by more than two-fold. In a simulation study, MODMatcher provided more robust results by using three types of omics data than two types of omics data. We further demonstrate that MODMatcher can be broadly applied to large genomic data sets containing multiple types of omics data, such as The Cancer Genome Atlas (TCGA) data sets.

Genetic regulation of gene expression in the lung identifies CST3 and CD22 as potential causal genes for airflow obstruction.

BACKGROUND: COPD is a complex chronic disease with poorly understood pathogenesis. Integrative genomic approaches have the potential to elucidate the biological networks underlying COPD and lung function. We recently combined genome-wide genotyping and gene expression in 1111 human lung specimens to map expression quantitative trait loci (eQTL). OBJECTIVE: To determine causal associations between COPD and lung function-associated single nucleotide polymorphisms (SNPs) and lung tissue gene expression changes in our lung eQTL dataset. METHODS: We evaluated causality between SNPs and gene expression for three COPD phenotypes: FEV(1)% predicted, FEV(1)/FVC and COPD as a categorical variable. Different models were assessed in the three cohorts independently and in a meta-analysis. SNPs associated with a COPD phenotype and gene expression were subjected to causal pathway modelling and manual curation. In silico analyses evaluated functional enrichment of biological pathways among newly identified causal genes. Biologically relevant causal genes were validated in two separate gene expression datasets of lung tissues and bronchial airway brushings. RESULTS: High reliability causal relations were found in SNP-mRNA-phenotype triplets for FEV(1)% predicted (n=169) and FEV(1)/FVC (n=80). Several genes of potential biological relevance for COPD were revealed. eQTL-SNPs upregulating cystatin C (CST3) and CD22 were associated with worse lung function. Signalling pathways enriched with causal genes included xenobiotic metabolism, apoptosis, protease-antiprotease and oxidant-antioxidant balance. CONCLUSIONS: By using integrative genomics and analysing the relationships of COPD phenotypes with SNPs and gene expression in lung tissue, we identified CST3 and CD22 as potential causal genes for airflow obstruction. This study also augmented the understanding of previously described COPD pathways.

Endothelial Erg Regulates Expression of Pulmonary Lymphatic Junctional and Inflammation Genes in Mouse Lungs Impacting Lymphatic Transport.

The ETS transcription factor ERG is a master regulator of endothelial gene specificity and highly enriched in the capillary, vein, and arterial endothelial cells. ERG expression is critical for endothelial barrier function, permeability, and vascular inflammation. A dysfunctional vascular endothelial ERG has been shown to impair lung capillary homeostasis, contributing to pulmonary fibrosis as previously observed in IPF lungs. Our preliminary observations indicate that lymphatic endothelial cells (LEC) in the human IPF lung also lack ERG. To understand the role of ERG in pulmonary LECs, we developed LEC-specific inducible Erg-CKO and Erg-GFP-CKO conditional knockout (CKO) mice under Prox1 promoter. Whole lung microarray analysis, flow cytometry, and qPCR confirmed an inflammatory and pro-lymphvasculogenic predisposition in Erg-CKO lung. FITC-Dextran tracing analysis showed an increased pulmonary interstitial lymphatic fluid transport from the lung to the axial lymph node. Single-cell transcriptomics confirmed that genes associated with cell junction integrity were downregulated in Erg-CKO pre-collector and collector LECs. Integrating Single-cell transcriptomics and CellChatDB helped identify LEC specific communication pathways contributing to pulmonary inflammation, trans-endothelial migration, inflammation, and Endo-MT in Erg-CKO lung. Our findings suggest that downregulation of lymphatic Erg crucially affects LEC function, LEC permeability, pulmonary LEC communication pathways and lymphatic transcriptomics.

Repetitive head impacts induce neuronal loss and neuroinflammation in young athletes.

Repetitive head impacts (RHI) sustained from contact sports are the largest risk factor for chronic traumatic encephalopathy (CTE). Currently, CTE can only be diagnosed after death and the multicellular cascade of events that trigger initial hyperphosphorylated tau (p-tau) deposition remain unclear. Further, the symptoms endorsed by young individuals with early disease are not fully explained by the extent of p-tau deposition, severely hampering development of therapeutic interventions. Here, we show that RHI exposure associates with a multicellular response in young individuals (<51 years old) prior to the onset of CTE p-tau pathology that correlates with number of years of RHI exposure. Leveraging single nucleus RNA sequencing of tissue from 8 control, 9 RHI-exposed, and 11 low stage CTE individuals, we identify SPP1+ inflammatory microglia, angiogenic and inflamed endothelial cell profiles, reactive astrocytes, and altered synaptic gene expression in excitatory and inhibitory neurons in all individuals with exposure to RHI. Surprisingly, we also observe a significant loss of cortical sulcus layer 2/3 neurons in contact sport athletes compared to controls independent of p-tau pathology. These results provide robust evidence that multiple years of RHI exposure is sufficient to induce lasting cellular alterations that may underlie p-tau deposition and help explain the early clinical symptoms observed in young former contact sport athletes. Furthermore, these data identify specific cellular responses to repetitive head impacts that may direct future identification of diagnostic and therapeutic strategies for CTE.

Myc Dysregulation in Activated Macrophages Initiates Iron-Mediated Lipid Peroxidation that Fuels Type I Interferon and Compromises TB Resistance.

A quarter of human population is infected with Mycobacterium tuberculosis, but less than 10% of those infected develop clinical, mostly pulmonary, TB. To dissect mechanisms of susceptibility in immunocompetent individuals, we developed a genetically defined sst1-susceptible mouse model that uniquely reproduces a defining feature of human TB: development of necrotic lung lesions after infection with virulent Mtb. In this study, we explored the connectivity of the sst1-regulated pathways during prolonged macrophage activation with TNF. We determined that the aberrant response of the sst1-susceptible macrophages to TNF was primarily driven by conflicting Myc and antioxidant response pathways that resulted in a coordinated failure to properly sequester intracellular iron and activate ferroptosis inhibitor enzymes. Consequently, iron-mediated lipid peroxidation fueled IFNß superinduction and sustained the Type I Interferon (IFN-I) pathway hyperactivity that locked the sst1-susceptible macrophages in a state of unresolving stress and compromised their resistance to Mtb. The accumulation of the aberrantly activated, stressed, macrophages within granuloma microenvironment led to the local failure of anti-tuberculosis immunity and tissue necrosis. Our findings suggest a novel link between metabolic dysregulation in macrophages and susceptibility to TB, offering insights into potential therapeutic targets aimed at modulating macrophage function and improving TB control.

A single-sample microarray normalization method to facilitate personalized-medicine workflows.

Gene-expression microarrays allow researchers to characterize biological phenomena in a high-throughput fashion but are subject to technological biases and inevitable variabilities that arise during sample collection and processing. Normalization techniques aim to correct such biases. Most existing methods require multiple samples to be processed in aggregate; consequently, each sample's output is influenced by other samples processed jointly. However, in personalized-medicine workflows, samples may arrive serially, so renormalizing all samples upon each new arrival would be impractical. We have developed Single Channel Array Normalization (SCAN), a single-sample technique that models the effects of probe-nucleotide composition on fluorescence intensity and corrects for such effects, dramatically increasing the signal-to-noise ratio within individual samples while decreasing variation across samples. In various benchmark comparisons, we show that SCAN performs as well as or better than competing methods yet has no dependence on external reference samples and can be applied to any single-channel microarray platform.

Characterization and decontamination of background noise in droplet-based single-cell protein expression data with DecontPro.

Assays such as CITE-seq can measure the abundance of cell surface proteins on individual cells using antibody derived tags (ADTs). However, many ADTs have high levels of background noise that can obfuscate down-stream analyses. Using an exploratory analysis of PBMC datasets, we find that some droplets that were originally called "empty" due to low levels of RNA contained high levels of ADTs and likely corresponded to neutrophils. We identified a novel type of artifact in the empty droplets called a "spongelet" which has medium levels of ADT expression and is distinct from ambient noise. ADT expression levels in the spongelets correlate to ADT expression levels in the background peak of true cells in several datasets suggesting that they can contribute to background noise along with ambient ADTs. We then developed DecontPro, a novel Bayesian hierarchical model that can decontaminate ADT data by estimating and removing contamination from these sources. DecontPro outperforms other decontamination tools in removing aberrantly expressed ADTs while retaining native ADTs and in improving clustering specificity. Overall, these results suggest that identification of empty drops should be performed separately for RNA and ADT data and that DecontPro can be incorporated into CITE-seq workflows to improve the quality of downstream analyses.

Characterization of highly active mutational signatures in tumors from a large Chinese population.

The majority of mutational signatures have been characterized in tumors from Western countries and the degree to which mutational signatures are similar or different in Eastern populations has not been fully explored. We leveraged a large-scale clinical sequencing cohort of tumors from a Chinese population containing 25 tumor types and found that the highly active mutational signatures were similar to those previously characterized1,2. The aristolochic acid signature SBS22 was observed in four soft tissue sarcomas and the POLE-associated signature SBS10 was observed in a gallbladder carcinoma. In lung adenocarcinoma, the polycyclic aromatic hydrocarbon (PAH) signature SBS4 was significantly higher in males compared to females but not associated with smoking status. The UV-associated signature SBS7 was significantly lower in cutaneous melanomas from the Chinese population compared to a similar American cohort. Overall, these results add to our understanding of the mutational processes that contribute to tumors from the Chinese population.