Systematic analysis of Mendelian disease-associated gene variants reveals new classes of cancer-predisposing genes.

BACKGROUND: Despite the acceleration of somatic driver gene discovery facilitated by recent large-scale tumor sequencing data, the contribution of inherited variants remains largely unexplored, primarily focusing on previously known cancer predisposition genes (CPGs) due to the low statistical power associated with detecting rare pathogenic variant-phenotype associations. METHODS: Here, we introduce a generalized log-regression model to measure the excess of pathogenic variants within genes in cancer patients compared to control samples. It aims to measure gene-level cancer risk enrichment by collapsing rare pathogenic variants after controlling the population differences across samples. RESULTS: In this study, we investigate whether pathogenic variants in Mendelian disease-associated genes (OMIM genes) are enriched in cancer patients compared to controls. Utilizing data from PCAWG and the 1,000 Genomes Project, we identify 103 OMIM genes demonstrating significant enrichment of pathogenic variants in cancer samples (FDR 20%). Through an integrative approach considering three distinct properties, we classify these CPG-like OMIM genes into four clusters, indicating potential diverse mechanisms underlying tumor progression. Further, we explore the function of PAH (a key metabolic enzyme associated with Phenylketonuria), the gene exhibiting the highest prevalence of pathogenic variants in a pan-cancer (1.8%) compared to controls (0.6%). CONCLUSIONS: Our findings suggest a possible cancer progression mechanism through metabolic profile alterations. Overall, our data indicates that pathogenic OMIM gene variants contribute to cancer progression and introduces new CPG classifications potentially underpinning diverse tumorigenesis mechanisms.

Genetic assessment of pathogenic germline alterations in lysosomal genes among Asian patients with pancreatic ductal adenocarcinoma.

BACKGROUND: Lysosomes are closely linked to autophagic activity, which plays a vital role in pancreatic ductal adenocarcinoma (PDAC) biology. The survival of PDAC patients is still poor, and the identification of novel genetic factors for prognosis and treatment is highly required to prevent PDAC-related deaths. This study investigated the germline variants related to lysosomal dysfunction in patients with PDAC and to analyze whether they contribute to the development of PDAC. METHODS: The germline putative pathogenic variants (PPV) in genes involved in lysosomal storage disease (LSD) was compared between patients with PDAC (n = 418) and healthy controls (n = 845) using targeted panel and whole-exome sequencing. Furthermore, pancreatic organoids from wild-type and KrasG12D mice were used to evaluate the effect of lysosomal dysfunction on PDAC development. RNA sequencing (RNA-seq) analysis was performed with established PDAC patient-derived organoids (PDOs) according to the PPV status. RESULTS: The PPV in LSD-related genes was higher in patients with PDAC than in healthy controls (8.13 vs. 4.26%, Log2 OR = 1.65, P = 3.08 × 10-3). The PPV carriers of LSD-related genes with PDAC were significantly younger than the non-carriers (mean age 61.5 vs. 65.3 years, P = 0.031). We further studied a variant of the lysosomal enzyme, galactosylceramidase (GALC), which was the most frequently detected LSD variant in our cohort. Autophagolysosomal activity was hampered when GALC was downregulated, which was accompanied by paradoxically elevated autophagic flux. Furthermore, the number of proliferating Ki-67+ cells increased significantly in pancreatic organoids derived from Galc knockout KrasG12D mice. Moreover, GALC PPV carriers tended to show drug resistance in both PDAC cell line and PDAC PDO, and RNA-seq analysis revealed that various metabolism and gene repair pathways were upregulated in PDAC PDOs harboring a GALC variant. CONCLUSIONS: Genetically defined lysosomal dysfunction is frequently observed in patients with young-onset PDAC. This might contribute to PDAC development by altering metabolism and impairing autophagolysosomal activity, which could be potentially implicated in therapeutic applications for PDAC.

The impact of rare germline variants on human somatic mutation processes.

Somatic mutations are an inevitable component of ageing and the most important cause of cancer. The rates and types of somatic mutation vary across individuals, but relatively few inherited influences on mutation processes are known. We perform a gene-based rare variant association study with diverse mutational processes, using human cancer genomes from over 11,000 individuals of European ancestry. By combining burden and variance tests, we identify 207 associations involving 15 somatic mutational phenotypes and 42 genes that replicated in an independent data set at a false discovery rate of 1%. We associate rare inherited deleterious variants in genes such as MSH3, EXO1, SETD2, and MTOR with two phenotypically different forms of DNA mismatch repair deficiency, and variants in genes such as EXO1, PAXIP1, RIF1, and WRN with deficiency in homologous recombination repair. In addition, we identify associations with other mutational processes, such as APEX1 with APOBEC-signature mutagenesis. Many of the genes interact with each other and with known mutator genes within cellular sub-networks. Considered collectively, damaging variants in the identified genes are prevalent in the population. We suggest that rare germline variation in diverse genes commonly impacts mutational processes in somatic cells.

Higher order genetic interactions switch cancer genes from two-hit to one-hit drivers.

The classic two-hit model posits that both alleles of a tumor suppressor gene (TSG) must be inactivated to cause cancer. In contrast, for some oncogenes and haploinsufficient TSGs, a single genetic alteration can suffice to increase tumor fitness. Here, by quantifying the interactions between mutations and copy number alterations (CNAs) across 10,000 tumors, we show that many cancer genes actually switch between acting as one-hit or two-hit drivers. Third order genetic interactions identify the causes of some of these switches in dominance and dosage sensitivity as mutations in other genes in the same biological pathway. The correct genetic model for a gene thus depends on the other mutations in a genome, with a second hit in the same gene or an alteration in a different gene in the same pathway sometimes representing alternative evolutionary paths to cancer.

A comprehensive analysis of prefoldins and their implication in cancer.

Prefoldins (PFDNs) are evolutionary conserved co-chaperones, initially discovered in archaea but universally present in eukaryotes. PFDNs are prevalently organized into hetero-hexameric complexes. Although they have been overlooked since their discovery and their functions remain elusive, several reports indicate they act as co-chaperones escorting misfolded or non-native proteins to group II chaperonins. Unlike the eukaryotic PFDNs which interact with cytoskeletal components, the archaeal PFDNs can bind and stabilize a wide range of substrates, possibly due to their great structural diversity. The discovery of the unconventional RPB5 interactor (URI) PFDN-like complex (UPC) suggests that PFDNs have versatile functions and are required for different cellular processes, including an important role in cancer. Here, we summarize their functions across different species. Moreover, a comprehensive analysis of PFDNs genomic alterations across cancer types by using large-scale cancer genomic data indicates that PFDNs are a new class of non-mutated proteins significantly overexpressed in some cancer types.

Lipid-Oriented Live-Cell Distinction of B and T Lymphocytes.

The identification of each cell type is essential for understanding multicellular communities. Antibodies set as biomarkers have been the main toolbox for cell-type recognition, and chemical probes are emerging surrogates. Herein we report the first small-molecule probe, CDgB, to discriminate B lymphocytes from T lymphocytes, which was previously impossible without the help of antibodies. Through the study of the origin of cell specificity, we discovered an unexpected novel mechanism of membrane-oriented live-cell distinction. B cells maintain higher flexibility in their cell membrane than T cells and accumulate the lipid-like probe CDgB more preferably. Because B and T cells share common ancestors, we tracked the cell membrane changes of the progenitor cells and disclosed the dynamic reorganization of the membrane properties over the lymphocyte differentiation progress. This study casts an orthogonal strategy for the small-molecule cell identifier and enriches the toolbox for live-cell distinction from complex cell communities.

Integrated Analysis of Germline and Tumor DNA Identifies New Candidate Genes Involved in Familial Colorectal Cancer.

Colorectal cancer (CRC) shows aggregation in some families but no alterations in the known hereditary CRC genes. We aimed to identify new candidate genes which are potentially involved in germline predisposition to familial CRC. An integrated analysis of germline and tumor whole-exome sequencing data was performed in 18 unrelated CRC families. Deleterious single nucleotide variants (SNV), short insertions and deletions (indels), copy number variants (CNVs) and loss of heterozygosity (LOH) were assessed as candidates for first germline or second somatic hits. Candidate tumor suppressor genes were selected when alterations were detected in both germline and somatic DNA, fulfilling Knudson's two-hit hypothesis. Somatic mutational profiling and signature analysis were also performed. A series of germline-somatic variant pairs were detected. In all cases, the first hit was presented as a rare SNV/indel, whereas the second hit was either a different SNV (3 genes) or LOH affecting the same gene (141 genes). BRCA2, BLM, ERCC2, RECQL, REV3L and RIF1 were among the most promising candidate genes for germline CRC predisposition. The identification of new candidate genes involved in familial CRC could be achieved by our integrated analysis. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.

Systematic discovery of germline cancer predisposition genes through the identification of somatic second hits.

The genetic causes of cancer include both somatic mutations and inherited germline variants. Large-scale tumor sequencing has revolutionized the identification of somatic driver alterations but has had limited impact on the identification of cancer predisposition genes (CPGs). Here we present a statistical method, ALFRED, that tests Knudson's two-hit hypothesis to systematically identify CPGs from cancer genome data. Applied to ~10,000 tumor exomes the approach identifies known and putative CPGs - including the chromatin modifier NSD1 - that contribute to cancer through a combination of rare germline variants and somatic loss-of-heterozygosity (LOH). Rare germline variants in these genes contribute substantially to cancer risk, including to ~14% of ovarian carcinomas, ~7% of breast tumors, ~4% of uterine corpus endometrial carcinomas, and to a median of 2% of tumors across 17 cancer types.

Cancer type-dependent genetic interactions between cancer driver alterations indicate plasticity of epistasis across cell types.

Cancers, like many diseases, are normally caused by combinations of genetic alterations rather than by changes affecting single genes. It is well established that the genetic alterations that drive cancer often interact epistatically, having greater or weaker consequences in combination than expected from their individual effects. In a stringent statistical analysis of data from > 3,000 tumors, we find that the co-occurrence and mutual exclusivity relationships between cancer driver alterations change quite extensively in different types of cancer. This cannot be accounted for by variation in tumor heterogeneity or unrecognized cancer subtypes. Rather, it suggests that how genomic alterations interact cooperatively or partially redundantly to driver cancer changes in different types of cancers. This re-wiring of epistasis across cell types is likely to be a basic feature of genetic architecture, with important implications for understanding the evolution of multicellularity and human genetic diseases. In addition, if this plasticity of epistasis across cell types is also true for synthetic lethal interactions, a synthetic lethal strategy to kill cancer cells may frequently work in one type of cancer but prove ineffective in another.

Discovery of cellular proteins required for the early steps of HCV infection using integrative genomics.

Successful viral infection requires intimate communication between virus and host cell, a process that absolutely requires various host proteins. However, current efforts to discover novel host proteins as therapeutic targets for viral infection are difficult. Here, we developed an integrative-genomics approach to predict human genes involved in the early steps of hepatitis C virus (HCV) infection. By integrating HCV and human protein associations, co-expression data, and tight junction-tetraspanin web specific networks, we identified host proteins required for the early steps in HCV infection. Moreover, we validated the roles of newly identified proteins in HCV infection by knocking down their expression using small interfering RNAs. Specifically, a novel host factor CD63 was shown to directly interact with HCV E2 protein. We further demonstrated that an antibody against CD63 blocked HCV infection, indicating that CD63 may serve as a new therapeutic target for HCV-related diseases. The candidate gene list provides a source for identification of new therapeutic targets.

Spatial and functional organization of mitochondrial protein network.

Characterizing the spatial organization of the human mitochondrial proteome will enhance our understanding of mitochondrial functions at the molecular level and provide key insight into protein-disease associations. However, the sub-organellar location and possible association with mitochondrial diseases are not annotated for most mitochondrial proteins. Here, we characterized the functional and spatial organization of mitochondrial proteins by assessing their position in the Mitochondrial Protein Functional (MPF) network. Network position was assigned to the MPF network and facilitated the determination of sub-organellar location and functional organization of mitochondrial proteins. Moreover, network position successfully identified candidate disease genes of several mitochondrial disorders. Thus, our data support the use of network position as a novel method to explore the molecular function and pathogenesis of mitochondrial proteins.

Epigenetic epistatic interactions constrain the evolution of gene expression.

Reduced activity of two genes in combination often has a more detrimental effect than expected. Such epistatic interactions not only occur when genes are mutated but also due to variation in gene expression, including among isogenic individuals in a controlled environment. We hypothesized that these 'epigenetic' epistatic interactions could place important constraints on the evolution of gene expression. Consistent with this, we show here that yeast genes with many epistatic interaction partners typically show low expression variation among isogenic individuals and low variation across different conditions. In addition, their expression tends to remain stable in response to the accumulation of mutations and only diverges slowly between strains and species. Yeast promoter architectures, the retention of gene duplicates, and the divergence of expression between humans and chimps are also consistent with selective pressure to reduce the likelihood of harmful epigenetic epistatic interactions. Based on these and previous analyses, we propose that the tight regulation of epistatic interaction network hubs makes an important contribution to the maintenance of a robust, 'canalized' phenotype. Moreover, that epigenetic epistatic interactions may contribute substantially to fitness defects when single genes are deleted.

Evolutionary history of human disease genes reveals phenotypic connections and comorbidity among genetic diseases.

The extent to which evolutionary changes have impacted the phenotypic relationships among human diseases remains unclear. In this work, we report that phenotypically similar diseases are connected by the evolutionary constraints on human disease genes. Human disease groups can be classified into slowly or rapidly evolving classes, where the diseases in the slowly evolving class are enriched with morphological phenotypes and those in the rapidly evolving class are enriched with physiological phenotypes. Our findings establish a clear evolutionary connection between disease classes and disease phenotypes for the first time. Furthermore, the high comorbidity found between diseases connected by similar evolutionary constraints enables us to improve the predictability of the relative risk of human diseases. We find the evolutionary constraints on disease genes are a new layer of molecular connection in the network-based exploration of human diseases.

Rewiring of PDZ domain-ligand interaction network contributed to eukaryotic evolution.

PDZ domain-mediated interactions have greatly expanded during metazoan evolution, becoming important for controlling signal flow via the assembly of multiple signaling components. The evolutionary history of PDZ domain-mediated interactions has never been explored at the molecular level. It is of great interest to understand how PDZ domain-ligand interactions emerged and how they become rewired during evolution. Here, we constructed the first human PDZ domain-ligand interaction network (PDZNet) together with binding motif sequences and interaction strengths of ligands. PDZNet includes 1,213 interactions between 97 human PDZ proteins and 591 ligands that connect most PDZ protein-mediated interactions (98%) in a large single network via shared ligands. We examined the rewiring of PDZ domain-ligand interactions throughout eukaryotic evolution by tracing changes in the C-terminal binding motif sequences of the PDZ ligands. We found that interaction rewiring by sequence mutation frequently occurred throughout evolution, largely contributing to the growth of PDZNet. The rewiring of PDZ domain-ligand interactions provided an effective means of functional innovations in nervous system development. Our findings provide empirical evidence for a network evolution model that highlights the rewiring of interactions as a mechanism for the development of new protein functions. PDZNet will be a valuable resource to further characterize the organization of the PDZ domain-mediated signaling proteome.

The protein interaction network of extracellular vesicles derived from human colorectal cancer cells.

Various mammalian cells including tumor cells secrete extracellular vesicles (EVs), otherwise known as exosomes and microvesicles. EVs are nanosized bilayered proteolipids and play multiple roles in intercellular communication. Although many vesicular proteins have been identified, their functional interrelationships and the mechanisms of EV biogenesis remain unknown. By interrogating proteomic data using systems approaches, we have created a protein interaction network of human colorectal cancer cell-derived EVs which comprises 1491 interactions between 957 vesicular proteins. We discovered that EVs have well-connected clusters with several hub proteins similar to other subcellular networks. We also experimentally validated that direct protein interactions between cellular proteins may be involved in protein sorting during EV formation. Moreover, physically and functionally interconnected protein complexes form functional modules involved in EV biogenesis and functions. Specifically, we discovered that SRC signaling plays a major role in EV biogenesis, and confirmed that inhibition of SRC kinase decreased the intracellular biogenesis and cell surface release of EVs. Our study provides global insights into the cargo-sorting, biogenesis, and pathophysiological roles of these complex extracellular organelles.

Protein localization as a principal feature of the etiology and comorbidity of genetic diseases.

Proteins targeting the same subcellular localization tend to participate in mutual protein-protein interactions (PPIs) and are often functionally associated. Here, we investigated the relationship between disease-associated proteins and their subcellular localizations, based on the assumption that protein pairs associated with phenotypically similar diseases are more likely to be connected via subcellular localization. The spatial constraints from subcellular localization significantly strengthened the disease associations of the proteins connected by subcellular localizations. In particular, certain disease types were more prevalent in specific subcellular localizations. We analyzed the enrichment of disease phenotypes within subcellular localizations, and found that there exists a significant correlation between disease classes and subcellular localizations. Furthermore, we found that two diseases displayed high comorbidity when disease-associated proteins were connected via subcellular localization. We newly explained 7584 disease pairs by using the context of protein subcellular localization, which had not been identified using shared genes or PPIs only. Our result establishes a direct correlation between protein subcellular localization and disease association, and helps to understand the mechanism of human disease progression.

Construction of functional interaction networks through consensus localization predictions of the human proteome.

Characterizing the subcellular localization of a protein provides a key clue for understanding protein function. However, different protein localization prediction programs often deliver conflicting results regarding the localization of the same protein. As the number of available localization prediction programs continues to grow, there is a need for a consensus prediction approach. To address this need, we developed a consensus localization prediction method called ConLoc based on a large-scale, systematic integration of 13 available programs that make predictions for five major subcellular localizations (cytosol, extracellular, mitochondria, nucleus, and plasma membrane). The ability of ConLoc to accurately predict protein localization was substantially better than existing programs. Using ConLoc prediction, we built a localization-guided functional interaction network of the human proteome and mapped known disease associations within this network. We found a high degree of shared disease associations among functionally interacting proteins that are localized to the same cellular compartment. Thus, the use of consensus localization prediction, such as ConLoc, is a new approach for the identification of novel disease associated genes.