Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing.

While alternative splicing is known to diversify the functional characteristics of some genes, the extent to which protein isoforms globally contribute to functional complexity on a proteomic scale remains unknown. To address this systematically, we cloned full-length open reading frames of alternatively spliced transcripts for a large number of human genes and used protein-protein interaction profiling to functionally compare hundreds of protein isoform pairs. The majority of isoform pairs share less than 50% of their interactions. In the global context of interactome network maps, alternative isoforms tend to behave like distinct proteins rather than minor variants of each other. Interaction partners specific to alternative isoforms tend to be expressed in a highly tissue-specific manner and belong to distinct functional modules. Our strategy, applicable to other functional characteristics, reveals a widespread expansion of protein interaction capabilities through alternative splicing and suggests that many alternative "isoforms" are functionally divergent (i.e., "functional alloforms").

Widespread macromolecular interaction perturbations in human genetic disorders.

How disease-associated mutations impair protein activities in the context of biological networks remains mostly undetermined. Although a few renowned alleles are well characterized, functional information is missing for over 100,000 disease-associated variants. Here we functionally profile several thousand missense mutations across a spectrum of Mendelian disorders using various interaction assays. The majority of disease-associated alleles exhibit wild-type chaperone binding profiles, suggesting they preserve protein folding or stability. While common variants from healthy individuals rarely affect interactions, two-thirds of disease-associated alleles perturb protein-protein interactions, with half corresponding to "edgetic" alleles affecting only a subset of interactions while leaving most other interactions unperturbed. With transcription factors, many alleles that leave protein-protein interactions intact affect DNA binding. Different mutations in the same gene leading to different interaction profiles often result in distinct disease phenotypes. Thus disease-associated alleles that perturb distinct protein activities rather than grossly affecting folding and stability are relatively widespread.

A proteome-scale map of the human interactome network.

Just as reference genome sequences revolutionized human genetics, reference maps of interactome networks will be critical to fully understand genotype-phenotype relationships. Here, we describe a systematic map of ?14,000 high-quality human binary protein-protein interactions. At equal quality, this map is ?30% larger than what is available from small-scale studies published in the literature in the last few decades. While currently available information is highly biased and only covers a relatively small portion of the proteome, our systematic map appears strikingly more homogeneous, revealing a "broader" human interactome network than currently appreciated. The map also uncovers significant interconnectivity between known and candidate cancer gene products, providing unbiased evidence for an expanded functional cancer landscape, while demonstrating how high-quality interactome models will help "connect the dots" of the genomic revolution.

A reference map of the human binary protein interactome.

Global insights into cellular organization and genome function require comprehensive understanding of the interactome networks that mediate genotype-phenotype relationships1,2. Here we present a human 'all-by-all' reference interactome map of human binary protein interactions, or 'HuRI'. With approximately 53,000 protein-protein interactions, HuRI has approximately four times as many such interactions as there are high-quality curated interactions from small-scale studies. The integration of HuRI with genome3, transcriptome4 and proteome5 data enables cellular function to be studied within most physiological or pathological cellular contexts. We demonstrate the utility of HuRI in identifying the specific subcellular roles of protein-protein interactions. Inferred tissue-specific networks reveal general principles for the formation of cellular context-specific functions and elucidate potential molecular mechanisms that might underlie tissue-specific phenotypes of Mendelian diseases. HuRI is a systematic proteome-wide reference that links genomic variation to phenotypic outcomes.

Codon-specific translation reprogramming promotes resistance to targeted therapy.

Reprogramming of mRNA translation has a key role in cancer development and drug resistance 1 . However, the molecular mechanisms that are involved in this process remain poorly understood. Wobble tRNA modifications are required for specific codon decoding during translation2,3. Here we show, in humans, that the enzymes that catalyse modifications of wobble uridine 34 (U34) tRNA (U34 enzymes) are key players of the protein synthesis rewiring that is induced by the transformation driven by the BRAF V600E oncogene and by resistance to targeted therapy in melanoma. We show that BRAF V600E -expressing melanoma cells are dependent on U34 enzymes for survival, and that concurrent inhibition of MAPK signalling and ELP3 or CTU1 and/or CTU2 synergizes to kill melanoma cells. Activation of the PI3K signalling pathway, one of the most common mechanisms of acquired resistance to MAPK therapeutic agents, markedly increases the expression of U34 enzymes. Mechanistically, U34 enzymes promote glycolysis in melanoma cells through the direct, codon-dependent, regulation of the translation of HIF1A mRNA and the maintenance of high levels of HIF1α protein. Therefore, the acquired resistance to anti-BRAF therapy is associated with high levels of U34 enzymes and HIF1α. Together, these results demonstrate that U34 enzymes promote the survival and resistance to therapy of melanoma cells by regulating specific mRNA translation.

The HTLV-1 viral oncoproteins Tax and HBZ reprogram the cellular mRNA splicing landscape.

Viral infections are known to hijack the transcription and translation of the host cell. However, the extent to which viral proteins coordinate these perturbations remains unclear. Here we used a model system, the human T-cell leukemia virus type 1 (HTLV-1), and systematically analyzed the transcriptome and interactome of key effectors oncoviral proteins Tax and HBZ. We showed that Tax and HBZ target distinct but also common transcription factors. Unexpectedly, we also uncovered a large set of interactions with RNA-binding proteins, including the U2 auxiliary factor large subunit (U2AF2), a key cellular regulator of pre-mRNA splicing. We discovered that Tax and HBZ perturb the splicing landscape by altering cassette exons in opposing manners, with Tax inducing exon inclusion while HBZ induces exon exclusion. Among Tax- and HBZ-dependent splicing changes, we identify events that are also altered in Adult T cell leukemia/lymphoma (ATLL) samples from two independent patient cohorts, and in well-known cancer census genes. Our interactome mapping approach, applicable to other viral oncogenes, has identified spliceosome perturbation as a novel mechanism coordinated by Tax and HBZ to reprogram the transcriptome.

Rapid Whole Genome Sequencing Diagnoses and Guides Treatment in Critically Ill Children in Belgium in Less than 40 Hours.

Rapid Whole Genome Sequencing (rWGS) represents a valuable exploration in critically ill pediatric patients. Early diagnosis allows care to be adjusted. We evaluated the feasibility, turnaround time (TAT), yield, and utility of rWGS in Belgium. Twenty-one unrelated critically ill patients were recruited from the neonatal intensive care units, the pediatric intensive care unit, and the neuropediatric unit, and offered rWGS as a first tier test. Libraries were prepared in the laboratory of human genetics of the University of Liège using Illumina DNA PCR-free protocol. Sequencing was performed on a NovaSeq 6000 in trio for 19 and in duo for two probands. The TAT was calculated from the sample reception to the validation of results. Clinical utility data were provided by treating physicians. A definite diagnosis was reached in twelve (57.5%) patients in 39.80 h on average (range: 37.05-43.7). An unsuspected diagnosis was identified in seven patients. rWGS guided care adjustments in diagnosed patients, including a gene therapy, an off-label drug trial and two condition-specific treatments. We successfully implemented the fastest rWGS platform in Europe and obtained one of the highest rWGS yields. This study establishes the path for a nationwide semi-centered rWGS network in Belgium.

Methylglyoxal, a glycolysis metabolite, triggers metastasis through MEK/ERK/SMAD1 pathway activation in breast cancer.

BACKGROUND: Elevated aerobic glycolysis rate is a biochemical alteration associated with malignant transformation and cancer progression. This metabolic shift unavoidably generates methylglyoxal (MG), a potent inducer of dicarbonyl stress through the formation of advanced glycation end products (AGEs). We have previously shown that the silencing of glyoxalase 1 (GLO1), the main MG detoxifying enzyme, generates endogenous dicarbonyl stress resulting in enhanced growth and metastasis in vivo. However, the molecular mechanisms through which MG stress promotes metastasis development remain to be unveiled. METHODS: In this study, we used RNA sequencing analysis to investigate gene-expression profiling of GLO1-depleted breast cancer cells and we validated the regulated expression of selected genes of interest by RT-qPCR. Using in vitro and in vivo assays, we demonstrated the acquisition of a pro-metastatic phenotype related to dicarbonyl stress in MDA-MB-231, MDA-MB-468 and MCF7 breast cancer cellular models. Hyperactivation of MEK/ERK/SMAD1 pathway was evidenced using western blotting upon endogenous MG stress and exogenous MG treatment conditions. MEK and SMAD1 regulation of MG pro-metastatic signature genes in breast cancer cells was demonstrated by RT-qPCR. RESULTS: High-throughput transcriptome profiling of GLO1-depleted breast cancer cells highlighted a pro-metastatic signature that establishes novel connections between MG dicarbonyl stress, extracellular matrix (ECM) remodeling by neoplastic cells and enhanced cell migration. Mechanistically, we showed that these metastasis-related processes are functionally linked to MEK/ERK/SMAD1 cascade activation in breast cancer cells. We showed that sustained MEK/ERK activation in GLO1-depleted cells notably occurred through the down-regulation of the expression of dual specificity phosphatases in MG-stressed breast cancer cells. The use of carnosine and aminoguanidine, two potent MG scavengers, reversed MG stress effects in in vitro and in vivo experimental settings. CONCLUSIONS: These results uncover for the first time the key role of MG dicarbonyl stress in the induction of ECM remodeling and the activation of migratory signaling pathways, both in favor of enhanced metastatic dissemination of breast cancer cells. Importantly, the efficient inhibition of mitogen-activated protein kinase (MAPK) signaling using MG scavengers further emphasizes the need to investigate their therapeutic potential across different malignancies.

Analysis of Genes Associated With Monogenic Primary Immunodeficiency Identifies Rare Variants in XIAP in Patients With Crohn's Disease.

BACKGROUND & AIMS: A few rare monogenic primary immunodeficiencies (PIDs) are characterized by chronic intestinal inflammation that resembles Crohn's disease (CD). We investigated whether 23 genes associated with 10 of these monogenic disorders contain common, low-frequency, or rare variants that increase risk for CD. METHODS: Common and low frequency variants in 1 Mb loci centered on the candidate genes were analyzed using meta-data corresponding to genotypes of approximately 17,000 patients with CD or without CD (controls) in Europe. The contribution of rare variants was assessed by high-throughput sequencing of 4750 individuals, including 660 early-onset and/or familial cases among the 2390 patients with CD. Variants were expressed from vectors in SW480 or HeLa cells and functions of their products were analyzed in immunofluorescence, luciferase, immunoprecipitation, and immunoblot assays. RESULTS: We reproduced the association of the interleukin 10 locus with CD (P = .007), although none of the significantly associated variants modified the coding sequence of interleukin 10. We found XIAP to be significantly enriched for rare coding mutations in patients with CD vs controls (P = .02). We identified 4 previously unreported missense variants associated with CD. Variants in XIAP cause the PID X-linked lymphoproliferative disease type 2, yet none of the carriers of these variants had all the clinical features of X-linked lymphoproliferative disease type 2. Identified XIAP variants S123N, R233Q, and P257A were associated with an impaired activation of NOD2 signaling after muramyl dipeptide stimulation. CONCLUSIONS: In a systematic analysis of variants in 23 PID-associated genes, we confirmed the association of variants in XIAP with CD. Further screenings for CD-associated variants and analyses of their functions could increase our understanding of the relationship between PID-associated genes and CD pathogenesis.

Proto-genes and de novo gene birth.

Novel protein-coding genes can arise either through re-organization of pre-existing genes or de novo. Processes involving re-organization of pre-existing genes, notably after gene duplication, have been extensively described. In contrast, de novo gene birth remains poorly understood, mainly because translation of sequences devoid of genes, or 'non-genic' sequences, is expected to produce insignificant polypeptides rather than proteins with specific biological functions. Here we formalize an evolutionary model according to which functional genes evolve de novo through transitory proto-genes generated by widespread translational activity in non-genic sequences. Testing this model at the genome scale in Saccharomyces cerevisiae, we detect translation of hundreds of short species-specific open reading frames (ORFs) located in non-genic sequences. These translation events seem to provide adaptive potential, as suggested by their differential regulation upon stress and by signatures of retention by natural selection. In line with our model, we establish that S. cerevisiae ORFs can be placed within an evolutionary continuum ranging from non-genic sequences to genes. We identify ~1,900 candidate proto-genes among S. cerevisiae ORFs and find that de novo gene birth from such a reservoir may be more prevalent than sporadic gene duplication. Our work illustrates that evolution exploits seemingly dispensable sequences to generate adaptive functional innovation.

An inter-species protein-protein interaction network across vast evolutionary distance.

In cellular systems, biophysical interactions between macromolecules underlie a complex web of functional interactions. How biophysical and functional networks are coordinated, whether all biophysical interactions correspond to functional interactions, and how such biophysical-versus-functional network coordination is shaped by evolutionary forces are all largely unanswered questions. Here, we investigate these questions using an "inter-interactome" approach. We systematically probed the yeast and human proteomes for interactions between proteins from these two species and functionally characterized the resulting inter-interactome network. After a billion years of evolutionary divergence, the yeast and human proteomes are still capable of forming a biophysical network with properties that resemble those of intra-species networks. Although substantially reduced relative to intra-species networks, the levels of functional overlap in the yeast-human inter-interactome network uncover significant remnants of co-functionality widely preserved in the two proteomes beyond human-yeast homologs. Our data support evolutionary selection against biophysical interactions between proteins with little or no co-functionality. Such non-functional interactions, however, represent a reservoir from which nascent functional interactions may arise.

An efficient algorithm to perform multiple testing in epistasis screening.

BACKGROUND: Research in epistasis or gene-gene interaction detection for human complex traits has grown over the last few years. It has been marked by promising methodological developments, improved translation efforts of statistical epistasis to biological epistasis and attempts to integrate different omics information sources into the epistasis screening to enhance power. The quest for gene-gene interactions poses severe multiple-testing problems. In this context, the maxT algorithm is one technique to control the false-positive rate. However, the memory needed by this algorithm rises linearly with the amount of hypothesis tests. Gene-gene interaction studies will require a memory proportional to the squared number of SNPs. A genome-wide epistasis search would therefore require terabytes of memory. Hence, cache problems are likely to occur, increasing the computation time. In this work we present a new version of maxT, requiring an amount of memory independent from the number of genetic effects to be investigated. This algorithm was implemented in C++ in our epistasis screening software MBMDR-3.0.3. We evaluate the new implementation in terms of memory efficiency and speed using simulated data. The software is illustrated on real-life data for Crohn's disease. RESULTS: In the case of a binary (affected/unaffected) trait, the parallel workflow of MBMDR-3.0.3 analyzes all gene-gene interactions with a dataset of 100,000 SNPs typed on 1000 individuals within 4 days and 9 hours, using 999 permutations of the trait to assess statistical significance, on a cluster composed of 10 blades, containing each four Quad-Core AMD Opteron(tm) Processor 2352 2.1 GHz. In the case of a continuous trait, a similar run takes 9 days. Our program found 14 SNP-SNP interactions with a multiple-testing corrected p-value of less than 0.05 on real-life Crohn's disease (CD) data. CONCLUSIONS: Our software is the first implementation of the MB-MDR methodology able to solve large-scale SNP-SNP interactions problems within a few days, without using much memory, while adequately controlling the type I error rates. A new implementation to reach genome-wide epistasis screening is under construction. In the context of Crohn's disease, MBMDR-3.0.3 could identify epistasis involving regions that are well known in the field and could be explained from a biological point of view. This demonstrates the power of our software to find relevant phenotype-genotype higher-order associations.

'Edgetic' perturbation of a C. elegans BCL2 ortholog.

Genes and gene products do not function in isolation but within highly interconnected 'interactome' networks, modeled as graphs of nodes and edges representing macromolecules and interactions between them, respectively. We propose to investigate genotype-phenotype associations by methodical use of alleles that lack single interactions, while retaining all others, in contrast to genetic approaches designed to eliminate gene products completely. We describe an integrated strategy based on the reverse yeast two-hybrid system to isolate and characterize such edge-specific, or 'edgetic', alleles. We established a proof of concept with CED-9, a Caenorhabditis elegans BCL2 ortholog. Using ced-9 edgetic alleles, we uncovered a new potential functional link between apoptosis and a centrosomal protein. This approach is amenable to higher throughput and is particularly applicable to interactome network analysis in organisms for which transgenesis is straightforward.

Single-cell RNA sequencing of human liver reveals hepatic stellate cell heterogeneity.

BACKGROUND & AIMS: The multiple vital functions of the human liver are performed by highly specialised parenchymal and non-parenchymal cells organised in complex collaborative sinusoidal units. Although crucial for homeostasis, the cellular make-up of the human liver remains to be fully elucidated. Here, single-cell RNA-sequencing was used to unravel the heterogeneity of human liver cells, in particular of hepatocytes (HEPs) and hepatic stellate cells (HSCs). METHOD: The transcriptome of ~25,000 freshly isolated human liver cells was profiled using droplet-based RNA-sequencing. Recently published data sets and RNA in situ hybridisation were integrated to validate and locate newly identified cell populations. RESULTS: In total, 22 cell populations were annotated that reflected the heterogeneity of human parenchymal and non-parenchymal liver cells. More than 20,000 HEPs were ordered along the portocentral axis to confirm known, and reveal previously undescribed, zonated liver functions. The existence of 2 subpopulations of human HSCs with unique gene expression signatures and distinct intralobular localisation was revealed (i.e. portal and central vein-concentrated GPC3 + HSCs and perisinusoidally located DBH + HSCs). In particular, these data suggest that, although both subpopulations collaborate in the production and organisation of extracellular matrix, GPC3 + HSCs specifically express genes involved in the metabolism of glycosaminoglycans, whereas DBH + HSCs display a gene signature that is reminiscent of antigen-presenting cells. CONCLUSIONS: This study highlights metabolic zonation as a key determinant of HEP transcriptomic heterogeneity and, for the first time, outlines the existence of heterogeneous HSC subpopulations in the human liver. These findings call for further research on the functional implications of liver cell heterogeneity in health and disease. LAY SUMMARY: This study resolves the cellular landscape of the human liver in an unbiased manner and at high resolution to provide new insights into human liver cell biology. The results highlight the physiological heterogeneity of human hepatic stellate cells.

Edgetic perturbation models of human inherited disorders.

Cellular functions are mediated through complex systems of macromolecules and metabolites linked through biochemical and physical interactions, represented in interactome models as 'nodes' and 'edges', respectively. Better understanding of genotype-to-phenotype relationships in human disease will require modeling of how disease-causing mutations affect systems or interactome properties. Here we investigate how perturbations of interactome networks may differ between complete loss of gene products ('node removal') and interaction-specific or edge-specific ('edgetic') alterations. Global computational analyses of approximately 50,000 known causative mutations in human Mendelian disorders revealed clear separations of mutations probably corresponding to those of node removal versus edgetic perturbations. Experimental characterization of mutant alleles in various disorders identified diverse edgetic interaction profiles of mutant proteins, which correlated with distinct structural properties of disease proteins and disease mechanisms. Edgetic perturbations seem to confer distinct functional consequences from node removal because a large fraction of cases in which a single gene is linked to multiple disorders can be modeled by distinguishing edgetic network perturbations. Edgetic network perturbation models might improve both the understanding of dissemination of disease alleles in human populations and the development of molecular therapeutic strategies.

The X-linked trichothiodystrophy-causing gene RNF113A links the spliceosome to cell survival upon DNA damage.

Prolonged cell survival occurs through the expression of specific protein isoforms generated by alternate splicing of mRNA precursors in cancer cells. How alternate splicing regulates tumor development and resistance to targeted therapies in cancer remain poorly understood. Here we show that RNF113A, whose loss-of-function causes the X-linked trichothiodystrophy, is overexpressed in lung cancer and protects from Cisplatin-dependent cell death. RNF113A is a RNA-binding protein which regulates the splicing of multiple candidates involved in cell survival. RNF113A deficiency triggers cell death upon DNA damage through multiple mechanisms, including apoptosis via the destabilization of the prosurvival protein MCL-1, ferroptosis due to enhanced SAT1 expression, and increased production of ROS due to altered Noxa1 expression. RNF113A deficiency circumvents the resistance to Cisplatin and to BCL-2 inhibitors through the destabilization of MCL-1, which thus defines spliceosome inhibitors as a therapeutic approach to treat tumors showing acquired resistance to specific drugs due to MCL-1 stabilization.

Detection and identification of transgenic events by next generation sequencing combined with enrichment technologies.

Next generation sequencing (NGS) is a promising tool for analysing the quality and safety of food and feed products. The detection and identification of genetically modified organisms (GMOs) is complex, as the diversity of transgenic events and types of structural elements introduced in plants continue to increase. In this paper, we show how a strategy that combines enrichment technologies with NGS can be used to detect a large panel of structural elements and partially or completely reconstruct the new sequence inserted into the plant genome in a single analysis, even at low GMO percentages. The strategy of enriching sequences of interest makes the approach applicable even to mixed products, which was not possible before due to insufficient coverage of the different genomes present. This approach is also the first step towards a more complete characterisation of agrifood products in a single analysis.

Characterization of a novel aphid prenyltransferase displaying dual geranyl/farnesyl diphosphate synthase activity.

We report on the cDNA cloning and characterization of a novel short-chain isoprenyl diphosphate synthase from the aphid Myzus persicae. Of the three IPPS cDNAs we cloned, two yielded prenyltransferase activity following expression in Escherichia coli; these cDNAs encode identical proteins except for the presence, in one of them, of an N-terminal mitochondrial targeting peptide. Although the aphid enzyme was predicted to be a farnesyl diphosphate synthase by BLASTP analysis, rMpIPPS, when isopentenyl diphosphate and dimethylallyl diphosphate are supplied as substrates, typically generated geranyl diphosphate (C10) as its main product, along with significant quantities of farnesyl diphosphate (C15). Analysis of an MpIPPS homology model pointed to substitutions that could confer GPP/FPP synthase activity to the aphid enzyme.

Distantly related lipocalins share two conserved clusters of hydrophobic residues: use in homology modeling.

BACKGROUND: Lipocalins are widely distributed in nature and are found in bacteria, plants, arthropoda and vertebra. In hematophagous arthropods, they are implicated in the successful accomplishment of the blood meal, interfering with platelet aggregation, blood coagulation and inflammation and in the transmission of disease parasites such as Trypanosoma cruzi and Borrelia burgdorferi. The pairwise sequence identity is low among this family, often below 30%, despite a well conserved tertiary structure. Under the 30% identity threshold, alignment methods do not correctly assign and align proteins. The only safe way to assign a sequence to that family is by experimental determination. However, these procedures are long and costly and cannot always be applied. A way to circumvent the experimental approach is sequence and structure analyze. To further help in that task, the residues implicated in the stabilisation of the lipocalin fold were determined. This was done by analyzing the conserved interactions for ten lipocalins having a maximum pairwise identity of 28% and various functions. RESULTS: It was determined that two hydrophobic clusters of residues are conserved by analysing the ten lipocalin structures and sequences. One cluster is internal to the barrel, involving all strands and the 310 helix. The other is external, involving four strands and the helix lying parallel to the barrel surface. These clusters are also present in RaHBP2, a unusual "outlier" lipocalin from tick Rhipicephalus appendiculatus. This information was used to assess assignment of LIR2 a protein from Ixodes ricinus and to build a 3D model that helps to predict function. FTIR data support the lipocalin fold for this protein. CONCLUSION: By sequence and structural analyzes, two conserved clusters of hydrophobic residues in interactions have been identified in lipocalins. Since the residues implicated are not conserved for function, they should provide the minimal subset necessary to confer the lipocalin fold. This information has been used to assign LIR2 to lipocalins and to investigate its structure/function relationship. This study could be applied to other protein families with low pairwise similarity, such as the structurally related fatty acid binding proteins or avidins.