Macroalgal deep genomics illuminate multiple paths to aquatic, photosynthetic multicellularity

Macroalgae are multicellular, aquatic autotrophs that play vital roles in global climate maintenance and have diverse applications in biotechnology and eco-engineering, which are directly linked to their multicellularity phenotypes. However, their genomic diversity and the evolutionary mechanisms underlying multicellularity in these organisms remain uncharacterized. In this study, we sequenced 110 macroalgal genomes from diverse climates and phyla, and identified key genomic features that distinguish them from their microalgal relatives. Genes for cell adhesion, extracellular matrix formation, cell polarity, transport, and cell differentiation distinguish macroalgae from microalgae across all three major phyla, constituting conserved and unique gene sets supporting multicellular processes. Adhesome genes show phylum- and climate-specific expansions that may facilitate niche adaptation. Collectively, our study reveals genetic determinants of convergent and divergent evolutionary trajectories that have shaped morphological diversity in macroalgae and provides genome-wide frameworks to understand photosynthetic multicellular evolution in aquatic environments.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays or AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold-Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

Basal interaction of the orphan receptor GPR101 with arrestins leads to constitutive internalization

GPR101 is an orphan G protein-coupled receptor that promotes growth hormone secretion in the pituitary. The microduplication of the GPR101 gene has been linked with the X-linked acrogigantism, or X-LAG, syndrome. This disease is characterized by excessive growth hormone secretion and abnormal rapid growth beginning early in life. Mechanistically, GPR101 induces growth hormone secretion through constitutive activation of multiple heterotrimeric G proteins. However, the full scope of GPR101 signaling remains largely elusive. Herein, we investigated the association of GPR101 to multiple transducers and uncovered an important basal interaction with Arrestin 2 (β-arrestin 1) and Arrestin 3 (β-arrestin 2). By using a GPR101 mutant lacking the C-terminus and cell lines with an Arrestin 2/3 null background, we show that the arrestin association leads to constitutive clathrin- and dynamin-mediated GPR101 internalization. To further highlight GPR101 intracellular fate, we assessed the colocalization of GPR101 with Rab protein markers. Internalized GPR101 was mainly colocalized with the early endosome markers, Rab5 and EEA-1, and to a lesser degree with the late endosome marker Rab7. However, GPR101 was not colocalized with the recycling endosome marker Rab11. These findings show that the basal arrestin recruitment by GPR101 C-terminal tail drives the receptor constitutive clathrin-mediated internalization. Intracellularly, GPR101 concentrates in the endosomal compartment and is degraded through the lysosomal pathway. In conclusion, we uncovered a constitutive intracellular trafficking of GPR101 that potentially represents an important layer of regulation of its signaling and function.

Identification of small molecule antivirals against HTLV-1 by targeting the hDLG1-Tax-1 protein-protein interaction

Human T-cell leukemia virus type-1 (HTLV-1) is the first pathogenic retrovirus discovered in human. Although HTLV-1-induced diseases are well-characterized and linked to the encoded Tax-1 oncoprotein, there is currently no strategy to target Tax-1 functions with small molecules. Here, we analyzed the binding of Tax-1 to the human homolog of the drosophila discs large tumor suppressor (hDLG1/SAP97), a multi-domain scaffolding protein involved in Tax-1-transformation ability. We have solved the structures of the PDZ binding motif (PBM) of Tax-1 in complex with the PDZ1 and PDZ2 domains of hDLG1 and assessed the binding of 10 million molecules by virtual screening. Among the 19 experimentally confirmed compounds, one systematically inhibited the Tax-1-hDLG1 interaction in different biophysical and cellular assays, as well as HTLV-1 cell-to-cell transmission in a T-cell model. Thus, our work demonstrates that interactions involving Tax-1 PDZ-domains are amenable to small-molecule inhibition, which provides a framework for the design of targeted therapies for HTLV-1-induced diseases.

The helicase-like transcription factor redirects the autophagic flux and restricts human T cell leukemia virus type 1 infection

Retroviruses and their host have coevolved in a delicate balance between viral replication and survival of the infected cell. In this equilibrium, restriction factors expressed by infected cells control different steps of retroviral replication such as entry, uncoating, nuclear import, expression, or budding. Here, we describe a mechanism of restriction against human T cell leukemia virus type 1 (HTLV-1) by the helicase-like transcription factor (HLTF). We show that RNA and protein levels of HLTF are reduced in primary T cells of HTLV-1-infected subjects, suggesting a clinical relevance. We further demonstrate that the viral oncogene Tax represses HLTF transcription via the Enhancer of zeste homolog 2 methyltransferase of the Polycomb repressive complex 2. The Tax protein also directly interacts with HLTF and induces its proteasomal degradation. RNA interference and gene transduction in HTLV-1-infected T cells derived from patients indicate that HLTF is a restriction factor. Restoring the normal levels of HLTF expression induces the dispersal of the Golgi apparatus and overproduction of secretory granules. By synergizing with Tax-mediated NF-κB activation, physiologically relevant levels of HLTF intensify the autophagic flux. Increased vesicular trafficking leads to an enlargement of the lysosomes and the production of large vacuoles containing viral particles. HLTF induction in HTLV-1-infected cells significantly increases the percentage of defective virions. In conclusion, HLTF-mediated activation of the autophagic flux blunts the infectious replication cycle of HTLV-1, revealing an original mode of viral restriction.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays and AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

Homo cerevisiae-Leveraging Yeast for Investigating Protein-Protein Interactions and Their Role in Human Disease

Understanding how genetic variation affects phenotypes represents a major challenge, particularly in the context of human disease. Although numerous disease-associated genes have been identified, the clinical significance of most human variants remains unknown. Despite unparalleled advances in genomics, functional assays often lack sufficient throughput, hindering efficient variant functionalization. There is a critical need for the development of more potent, high-throughput methods for characterizing human genetic variants. Here, we review how yeast helps tackle this challenge, both as a valuable model organism and as an experimental tool for investigating the molecular basis of phenotypic perturbation upon genetic variation. In systems biology, yeast has played a pivotal role as a highly scalable platform which has allowed us to gain extensive genetic and molecular knowledge, including the construction of comprehensive interactome maps at the proteome scale for various organisms. By leveraging interactome networks, one can view biology from a systems perspective, unravel the molecular mechanisms underlying genetic diseases, and identify therapeutic targets. The use of yeast to assess the molecular impacts of genetic variants, including those associated with viral interactions, cancer, and rare and complex diseases, has the potential to bridge the gap between genotype and phenotype, opening the door for precision medicine approaches and therapeutic development.

Assessment of community efforts to advance network-based prediction of protein-protein interactions

Comprehensive understanding of the human protein-protein interaction (PPI) network, aka the human interactome, can provide important insights into the molecular mechanisms of complex biological processes and diseases. Despite the remarkable experimental efforts undertaken to date to determine the structure of the human interactome, many PPIs remain unmapped. Computational approaches, especially network-based methods, can facilitate the identification of previously uncharacterized PPIs. Many such methods have been proposed. Yet, a systematic evaluation of existing network-based methods in predicting PPIs is still lacking. Here, we report community efforts initiated by the International Network Medicine Consortium to benchmark the ability of 26 representative network-based methods to predict PPIs across six different interactomes of four different organisms: A. thaliana, C. elegans, S. cerevisiae, and H. sapiens. Through extensive computational and experimental validations, we found that advanced similarity-based methods, which leverage the underlying network characteristics of PPIs, show superior performance over other general link prediction methods in the interactomes we considered.

A proteome-scale map of the SARS-CoV-2-human contactome

Understanding the mechanisms of coronavirus disease 2019 (COVID-19) disease severity to efficiently design therapies for emerging virus variants remains an urgent challenge of the ongoing pandemic. Infection and immune reactions are mediated by direct contacts between viral molecules and the host proteome, and the vast majority of these virus-host contacts (the 'contactome') have not been identified. Here, we present a systematic contactome map of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with the human host encompassing more than 200 binary virus-host and intraviral protein-protein interactions. We find that host proteins genetically associated with comorbidities of severe illness and long COVID are enriched in SARS-CoV-2 targeted network communities. Evaluating contactome-derived hypotheses, we demonstrate that viral NSP14 activates nuclear factor κB (NF-κB)-dependent transcription, even in the presence of cytokine signaling. Moreover, for several tested host proteins, genetic knock-down substantially reduces viral replication. Additionally, we show for USP25 that this effect is phenocopied by the small-molecule inhibitor AZ1. Our results connect viral proteins to human genetic architecture for COVID-19 severity and offer potential therapeutic targets.

Human papillomavirus E6/E7 oncoproteins promote radiotherapy-mediated tumor suppression by globally hijacking host DNA damage repair

Rationale: Whatever the mucosa primary infected, HPV-positive cancers are traditionally associated with a favorable outcome, attributable to a high sensitivity to radiation therapy. However, the direct impact of viral E6/E7 oncoproteins on the intrinsic cellular radiosensitivity (and, globally, on host DNA repair) remains mostly speculative. Methods: Using several isogenic cell models expressing HPV16 E6 and/or E7, the effect of viral oncoproteins on global DNA damage response was first investigated by in vitro/in vivo approaches. The binary interactome of each individual HPV oncoprotein with factors involved in the various host DNA damage/repair mechanisms was then precisely mapped by Gaussia princeps luciferase complementation assay (and validated by co-immunoprecipitation). The stability/half-life of protein targets for HPV E6 and/or E7 as well as their subcellular localizations were determined. At last, the host genome integrity following E6/E7 expression and the synergy between radiotherapy and compounds targeting DNA repair were analyzed. Results: We first showed that the sole expression of one viral oncoprotein from HPV16 was able to significantly increase the sensitivity to irradiation of cells without affecting their basal viability parameters. In total, 10 novel targets (CHEK2, CLK2, CLK2/3, ERCC3, MNAT1, PER1, RMI1, RPA1, UVSSA and XRCC6) for E6 and 11 (ALKBH2, CHEK2, DNA2, DUT, ENDOV, ERCC3, PARP3, PMS1, PNKP, POLDIP2 and RBBP8) for E7 were identified. Importantly, not degraded following their interaction with E6 or E7, these proteins have been shown to be less linked to host DNA and to colocalize with HPV replication foci, denoting their crucial implication in viral life cycle. Finally, we found that E6/E7 oncoproteins globally jeopardize host genome integrity, increase the cellular sensitivity to DNA repair inhibitors and enhance their synergy with radiotherapy. Conclusion: Taken together, our findings provide a molecular insight into the direct hijacking of host DNA damage/repair responses by HPV oncoproteins, demonstrate the significant impact of this phenomenon on both intrinsic cellular radiosensitivity and host DNA integrity and suggest novel connected therapeutic vulnerabilities.

Closing the Gap between Bio-Based and Petroleum-Based Plastic through Bioengineering

Bioplastics, which are plastic materials produced from renewable bio-based feedstocks, have been investigated for their potential as an attractive alternative to petroleum-based plastics. Despite the harmful effects of plastic accumulation in the environment, bioplastic production is still underdeveloped. Recent advances in strain development, genome sequencing, and editing technologies have accelerated research efforts toward bioplastic production and helped to advance its goal of replacing conventional plastics. In this review, we highlight bioengineering approaches, new advancements, and related challenges in the bioproduction and biodegradation of plastics. We cover different types of polymers, including polylactic acid (PLA) and polyhydroxyalkanoates (PHAs and PHBs) produced by bacterial, microalgal, and plant species naturally as well as through genetic engineering. Moreover, we provide detailed information on pathways that produce PHAs and PHBs in bacteria. Lastly, we present the prospect of using large-scale genome engineering to enhance strains and develop microalgae as a sustainable production platform.

Bioscreening and pre-clinical evaluation of the impact of bioactive molecules from Ptychotis verticillata on the multilineage potential of mesenchymal stromal cells towards immune- and inflammation-mediated diseases

OBJECTIVE AND DESIGN

Mesenchymal stromal cells (MSCs) are currently used in cell reparative medicine due to their trophic and ant-inflammatory properties. The modulation of stem cell properties by phytochemicals has been suggested as a tool to empower their tissue repair capacity. In vitro, MSCs are characterized by their tri-lineage potential that holds great interest for tissue regeneration. Ptychotis Verticillata (PV), an aromatic and medicinal plant, may be thus used to modulate the in vitro multilineage potential of MSCs.

MATERIALS AND METHODS

We screened the impact of PV-derived essential oil and their bioactive molecules (thymol and carvacrol) on the in vitro multilineage potential of MSCs. Different concentrations and incubation times of these compounds were assessed during the osteogenesis and adipogenesis of MSCs.

RESULTS

The analysis of 75 conditions indicates that these compounds are biologically active by promoting two major differentiation lineages from MSCs. In a time- and dose-dependent manner, thymol and carvacrol increased the osteogenesis and adipogenesis.

CONCLUSION

According to these preliminary observations, the addition of PV extract may stimulate the tissue regenerative and repair functions of MSCs. Further optimization of compound extraction and characterization from PV as well as cell treatment conditions should increase their therapeutic value in combination with MSCs.

HPV infection alters vaginal microbiome through down-regulating host mucosal innate peptides used by Lactobacilli as amino acid sources

Despite the high prevalence of both cervico-vaginal human papillomavirus (HPV) infection and bacterial vaginosis (BV) worldwide, their causal relationship remains unclear. While BV has been presumed to be a risk factor for HPV acquisition and related carcinogenesis for a long time, here, supported by both a large retrospective follow-up study (n = 6,085) and extensive in vivo data using the K14-HPV16 transgenic mouse model, we report a novel blueprint in which the opposite association also exists. Mechanistically, by interacting with several core members (NEMO, CK1 and β-TrCP) of both NF-κB and Wnt/β-catenin signaling pathways, we show that HPV E7 oncoprotein greatly inhibits host defense peptide expression. Physiologically secreted by the squamous mucosa lining the lower female genital tract, we demonstrate that some of these latter are fundamental factors governing host-microbial interactions. More specifically, several innate molecules down-regulated in case of HPV infection are hydrolyzed, internalized and used by the predominant Lactobacillus species as amino acid source sustaining their growth/survival. Collectively, this study reveals a new viral immune evasion strategy which, by its persistent/negative impact on lactic acid bacteria, ultimately causes the dysbiosis of vaginal microbiota.

Extracellular Vesicle Proteomes Shed Light on the Evolutionary, Interactive, and Functional Divergence of Their Biogenesis Mechanisms

Extracellular vesicles (EVs) are membranous structures containing bioactive molecules, secreted by most cells into the extracellular environment. EVs are classified by their biogenesis mechanisms into two major subtypes: ectosomes (enriched in large EVs; lEVs), budding directly from the plasma membrane, which is common in both prokaryotes and eukaryotes, and exosomes (enriched in small EVs; sEVs) generated through the multivesicular bodies via the endomembrane system, which is unique to eukaryotes. Even though recent proteomic analyses have identified key proteins associated with EV subtypes, there has been no systematic analysis, thus far, to support the general validity and utility of current EV subtype separation methods, still largely dependent on physical properties, such as vesicular size and sedimentation. Here, we classified human EV proteomic datasets into two main categories based on distinct centrifugation protocols commonly used for isolating sEV or lEV fractions. We found characteristic, evolutionarily conserved profiles of sEV and lEV proteins linked to their respective biogenetic origins. This may suggest that the evolutionary trajectory of vesicular proteins may result in a membership bias toward specific EV subtypes. Protein-protein interaction (PPI) network analysis showed that vesicular proteins formed distinct clusters with proteins in the same EV fraction, providing evidence for the existence of EV subtype-specific protein recruiters. Moreover, we identified functional modules enriched in each fraction, including multivesicular body sorting for sEV, and mitochondria cellular respiration for lEV proteins. Our analysis successfully captured novel features of EVs embedded in heterogeneous proteomics studies and suggests specific protein markers and signatures to be used as quality controllers in the isolation procedure for subtype-enriched EV fractions.

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.

Tax and HBZ are two viral regulatory proteins encoded by the human T-cell leukemia virus type 1 (HTLV-1) via sense and antisense transcripts, respectively. Both proteins are known to drive oncogenic processes that culminate in a T-cell neoplasm, known as Adult T cell leukemia/lymphoma (ATLL). We measured the effects of Tax and HBZ on host gene expression pathway by analyzing the interactome with cellular transcriptional and post-transcriptional regulators, and the transcriptome and mRNA splicing of cell lines expressing either Tax or HBZ. We compared our results with data obtained from independent cohorts of Japanese and Afro-Caribbean patients, and identified common splicing changes that might represent clinically useful biomarkers for ATLL. Finally, we provide evidence that the viral protein Tax can reprogram initial steps of the T-cell transcriptome diversification by hijacking the U2AF complex, a key cellular regulator of pre-mRNA splicing.

Alternative glycosylation controls endoplasmic reticulum dynamics and tubular extension in mammalian cells

The endoplasmic reticulum (ER) is a central eukaryotic organelle with a tubular network made of hairpin proteins linked by hydrolysis of guanosine triphosphate nucleotides. Among posttranslational modifications initiated at the ER level, glycosylation is the most common reaction. However, our understanding of the impact of glycosylation on the ER structure remains unclear. Here, we show that exostosin-1 (EXT1) glycosyltransferase, an enzyme involved in N-glycosylation, is a key regulator of ER morphology and dynamics. We have integrated multiomics and superresolution imaging to characterize the broad effect of EXT1 inactivation, including the ER shape-dynamics-function relationships in mammalian cells. We have observed that inactivating EXT1 induces cell enlargement and enhances metabolic switches such as protein secretion. In particular, suppressing EXT1 in mouse thymocytes causes developmental dysfunctions associated with the ER network extension. Last, our data illuminate the physical and functional aspects of the ER proteome-glycome-lipidome structure axis, with implications in biotechnology and medicine.

Host PDZ-containing proteins targeted by SARS-CoV-2

Small linear motifs targeting protein interacting domains called PSD‐95/Dlg/ZO‐1 (PDZ) have been identified at the C terminus of the severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) proteins E, 3a, and N. Using a high‐throughput approach of affinity‐profiling against the full human PDZome, we identified sixteen human PDZ binders of SARS‐CoV‐2 proteins E, 3A, and N showing significant interactions with dissociation constants values ranging from 3 to 82 μm. Six of them (TJP1, PTPN13, HTRA1, PARD3, MLLT4, LNX2) are also recognized by SARS‐CoV while three (NHERF1, MAST2, RADIL) are specific to SARS‐CoV‐2 E protein. Most of these SARS‐CoV‐2 protein partners are involved in cellular junctions/polarity and could be also linked to evasion mechanisms of the immune responses during viral infection. Among the binders of the SARS‐CoV‐2 proteins E, 3a, or N, seven significantly affect viral replication under knock down gene expression in infected cells. This PDZ profiling identifying human proteins potentially targeted by SARS‐CoV‐2 can help to understand the multifactorial severity of COVID19 and to conceive effective anti‐coronaviral agents for therapeutic purposes.

Non-cytotoxic 1,2,3-triazole tethered fused heterocyclic ring derivatives display Tax protein inhibition and impair HTLV-1 infected cells

Human T cell lymphotropic virus type 1 (HTLV-1) is a human retrovirus that infects approximately 10-20 million people worldwide and causes an aggressive neoplasia (adult T-cell leukemia/lymphoma - ATL). Therapeutic approaches for the treatment of ATL have variable effectiveness and poor prognosis, thus requiring strategies to identify novel compounds with activity on infected cells. In this sense, we initially screened a small series of 25 1,2,3-triazole derivatives to discover cell proliferation inhibitors and apoptosis inducers in HTLV-1-infected T-cell line (MT-2) for further assessment of their effect on viral tax activity through inducible-tax reporter cell line (Jurkat LTR-GFP). Eight promising compounds (02, 05, 06, 13, 15, 21, 22 and 25) with activity ≥70% were initially selected, based on a suitable cell-based assay using resazurin reduction method, and evaluated towards cell cycle, apoptosis and Tax/GFP expression analyses through flow cytometry. Compound 02 induced S phase cell cycle arrest and compounds 05, 06, 22 and 25 promoted apoptosis. Remarkably, compounds 22 and 25 also reduced GFP expression in an inducible-tax reporter cell, which suggests an effect on Tax viral protein. More importantly, compounds 02, 22 and 25 were not cytotoxic in human hepatoma cell line (Huh-7). Therefore, the discovery of 3 active and non-cytotoxic compounds against HTLV-1-infected cells can potentially contribute, as an initial promising strategy, to the development process of new drugs against ATL.

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 relationships^1,2. Here, we present a human “all-by-all” reference interactome map of human binary protein interactions, or “HuRI”. With ~53,000 high-quality protein-protein interactions (PPIs), HuRI has approximately four times more such interactions than high-quality curated interactions from small-scale studies. Integrating HuRI with genome³, transcriptome⁴, and proteome⁵ data enables the study of cellular function within most physiological or pathological cellular contexts. We demonstrate the utility of HuRI in identifying specific subcellular roles of PPIs. Inferred tissue-specific networks reveal general principles for the formation of cellular context-specific functions and elucidate potential molecular mechanisms underlying tissue-specific phenotypes of Mendelian diseases. HuRI represents a systematic proteome-wide reference linking genomic variation to phenotypic outcomes.

Maximizing binary interactome mapping with a minimal number of assays

Complementary assays are required to comprehensively map complex biological entities such as genomes, proteomes and interactome networks. However, how various assays can be optimally combined to approach completeness while maintaining high precision often remains unclear. Here, we propose a framework for binary protein-protein interaction (PPI) mapping based on optimally combining assays and/or assay versions to maximize detection of true positive interactions, while avoiding detection of random protein pairs. We have engineered a novel NanoLuc two-hybrid (N2H) system that integrates 12 different versions, differing by protein expression systems and tagging configurations. The resulting union of N2H versions recovers as many PPIs as 10 distinct assays combined. Thus, to further improve PPI mapping, developing alternative versions of existing assays might be as productive as designing completely new assays. Our findings should be applicable to systematic mapping of other biological landscapes.

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.

Activation of the Orphan G Protein-Coupled Receptor GPR27 by Surrogate Ligands Promotes β-Arrestin 2 Recruitment

G protein-coupled receptors are the most important drug targets for human diseases. An important number of them remain devoid of confirmed ligands. GPR27 is one of these orphan receptors, characterized by a high level of conservation among vertebrates and a predominant expression in the central nervous system. In addition, it has recently been linked to insulin secretion. However, the absence of endogenous or surrogate ligands for GPR27 complicates the examination of its biologic function. Our aim was to validate GPR27 signaling pathways, and therefore we sought to screen a diversity-oriented synthesis library to identify GPR27-specific surrogate agonists. To select an optimal screening assay, we investigated GPR27 ligand-independent activity. Both in G protein-mediated pathways and in β-arrestin 2 recruitment, no ligand-independent activity could be measured. However, we observed a recruitment of β-arrestin 2 to a GPR27V₂ chimera in the presence of membrane-anchored G protein-coupled receptor kinase-2. Therefore, we optimized a firefly luciferase complementation assay to screen against this chimeric receptor. We identified two compounds [N-[4-(anilinocarbonyl)phenyl]-2,4-dichlorobenzamide (ChemBridge, San Diego, CA; ID5128535) and 2,4-dichloro-N-{4-[(1,3-thiazol-2-ylamino)sulfonyl]phenyl}benzamide (ChemBridge ID5217941)] sharing a N-phenyl-2,4-dichlorobenzamide scaffold, which were selective for GPR27 over its closely related family members GPR85 and GPR173. The specificity of the activity was confirmed with a NanoLuc Binary Technology β-arrestin 2 assay, imaging of green fluorescent protein-tagged β-arrestin 2, and PathHunter β-arrestin 2 assay. Interestingly, no G protein activation was detected upon activation of GPR27 by these compounds. Our study provides the first selective surrogate agonists for the orphan GPR27.

The transcription factor ERG recruits CCR4-NOT to control mRNA decay and mitotic progression

Control of mRNA levels, a fundamental aspect in the regulation of gene expression, is achieved through a balance between mRNA synthesis and decay. E26-related gene (Erg) proteins are canonical transcription factors whose previously described functions are confined to the control of mRNA synthesis. Here, we report that ERG also regulates gene expression by affecting mRNA stability and identify the molecular mechanisms underlying this function in human cells. ERG is recruited to mRNAs via interaction with the RNA-binding protein RBPMS, and it promotes mRNA decay by binding CNOT2, a component of the CCR4-NOT deadenylation complex. Transcriptome-wide mRNA stability analysis revealed that ERG controls the degradation of a subset of mRNAs highly connected to Aurora signaling, whose decay during S phase is necessary for mitotic progression. Our data indicate that control of gene expression by mammalian transcription factors may follow a more complex scheme than previously anticipated, integrating mRNA synthesis and degradation.

Systematic interactome mapping of acute lymphoblastic leukemia cancer gene products reveals EXT-1 tumor suppressor as a Notch1 and FBWX7 common interactor

BACKGROUND

Perturbed genotypes in cancer can now be identified by whole genome sequencing of large number of diverse tumor samples, and observed gene mutations can be used for prognosis and classification of cancer subtypes. Although mutations in a few causative genes are directly linked to key signaling pathways perturbation, a global understanding of how known cancer genes drive oncogenesis in human is difficult to assess.

METHODS

We collected available information about mutated genes in Acute Lymphoblastic Leukemia (ALL). Validated human protein interactions (PPI) were collected from IntAct, HPRD and BioGRID interactomics databases, or obtained using yeast two-hybrid screening assay.

RESULTS

We have mapped interconnections between 116 cancer census gene products associated with ALL. Combining protein-protein interactions data and cancer-specific gene mutations information, we observed that 63 ALL-gene products are interconnected and identified 37 human proteins interacting with at least 2 ALL-gene products. We highlighted exclusive and coexistence genetic alterations in key signaling pathways including the PI3K/AKT and the NOTCH pathways. We then used different cell lines and reporter assay systems to validate the involvement of EXT1 in the Notch pathway.

CONCLUSION

We propose that novel ALL-gene candidates can be identified based on their functional association with well-known cancer genes. We identified EXT1, a gene not previously linked to ALL via mutations, as a common interactor of NOTCH1 and FBXW7 regulating the NOTCH pathway in an FBXW7-dependend manner.

A fungal biofilm reactor based on metal structured packing improves the quality of a Gla::GFP fusion protein produced by Aspergillus oryzae

Fungal biofilm is known to promote the excretion of secondary metabolites in accordance with solid-state-related physiological mechanisms. This work is based on the comparative analysis of classical submerged fermentation with a fungal biofilm reactor for the production of a Gla::green fluorescent protein (GFP) fusion protein by Aspergillus oryzae. The biofilm reactor comprises a metal structured packing allowing the attachment of the fungal biomass. Since the production of the target protein is under the control of the promoter glaB, specifically induced in solid-state fermentation, the biofilm mode of culture is expected to enhance the global productivity. Although production of the target protein was enhanced by using the biofilm mode of culture, we also found that fusion protein production is also significant when the submerged mode of culture is used. This result is related to high shear stress leading to biomass autolysis and leakage of intracellular fusion protein into the extracellular medium. Moreover, 2-D gel electrophoresis highlights the preservation of fusion protein integrity produced in biofilm conditions. Two fungal biofilm reactor designs were then investigated further, i.e. with full immersion of the packing or with medium recirculation on the packing, and the scale-up potentialities were evaluated. In this context, it has been shown that full immersion of the metal packing in the liquid medium during cultivation allows for a uniform colonization of the packing by the fungal biomass and leads to a better quality of the fusion protein.

Targeting oncogenic interleukin-7 receptor signalling with N-acetylcysteine in T cell acute lymphoblastic leukaemia

Activating mutations of the interleukin-7 receptor (IL7R) occur in approximately 10% of patients with T cell acute lymphoblastic leukaemia (T-ALL). Most mutations generate a cysteine at the transmembrane domain leading to receptor homodimerization through disulfide bond formation and ligand-independent activation of STAT5. We hypothesized that the reducing agent N-acetylcysteine (NAC), a well-tolerated drug used widely in clinical practice to treat acetaminophen overdose, would reduce disulfide bond formation, and inhibit mutant IL7R-mediated oncogenic signalling. We found that treatment with NAC disrupted IL7R homodimerization in IL7R-mutant DND-41 cells as assessed by non-reducing Western blot, as well as in a luciferase complementation assay. NAC led to STAT5 dephosphorylation and cell apoptosis at clinically achievable concentrations in DND-41 cells, and Ba/F3 cells transformed by an IL7R-mutant construct containing a cysteine insertion. The apoptotic effects of NAC could be rescued in part by a constitutively active allele of STAT5. Despite using doses lower than those tolerated in humans, NAC treatment significantly inhibited the progression of human DND-41 cells engrafted in immunodeficient mice. Thus, targeting leukaemogenic IL7R homodimerization with NAC offers a potentially effective and feasible therapeutic strategy that warrants testing in patients with T-ALL.

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.

Abstract 478: Predicting interactome networks perturbations in human cancer: application to gene fusions in acute lymphoblastic leukemia

Genomic variations such point mutations and gene fusions are directly or indirectly associated with human diseases. They are recognized as diagnostic, prognostic markers and as therapeutic targets. However, predicting the functional impact of these genetic alterations beyond affected genes and their products is challenging because diseased phenotypes likely result in complex molecular interaction networks. Using as models three different chromosomal translocations ETV6-RUNX1 (TEL-AML1), BCR-ABL1, and E2A-PBX1 (TCF3-PBX1) frequently found in B-cell acute lymphoblastic leukemia (B-ALL), we develop an approach to extract from gene expression changes, perturbed molecular interactions. We show that the Myc and JunD transcriptional circuits are specifically deregulated following TEL-AML1 and TCF3-PBX1 gene fusions, respectively. We also identified the bulk mRNA NXF1-dependent machinery as a direct target for the TCF3-PBX1 fusion protein. Through a novel approach combining gene expression and interactome data analysis, we provide new insight into TCF3-PBX1 and TEL-AML1 acute lymphoblastic leukemia.

Note: This abstract was not presented at the meeting.

Citation Format: Leon Hajingabo, Sarah Daakour, Maud Martin, Reinhard Grausenburger, Renate Panzer-Grümayer, Franck Dequiedt, Jacques Van Helden, Jean-Claude Twizere, Nicolas Simonis. Predicting interactome networks perturbations in human cancer: application to gene fusions in acute lymphoblastic leukemia. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 478. doi:10.1158/1538-7445.AM2014-478

Predicting interactome network perturbations in human cancer: application to gene fusions in acute lymphoblastic leukemia

Genomic variations such as gene fusions are directly or indirectly associated with human diseases. A method is presented combining gene expression and interactome data analyses to identify specific targets in leukemia. The Myc network and the mRNA export machinery are perturbed in ETV6-RUNX1 and TCF3-PBX1 subtypes of leukemia.

Genomic variations such as point mutations and gene fusions are directly or indirectly associated with human diseases. They are recognized as diagnostic, prognostic markers and therapeutic targets. However, predicting the functional effect of these genetic alterations beyond affected genes and their products is challenging because diseased phenotypes are likely dependent of complex molecular interaction networks. Using as models three different chromosomal translocations—ETV6-RUNX1 (TEL-AML1), BCR-ABL1, and TCF3-PBX1 (E2A-PBX1)—frequently found in precursor-B-cell acute lymphoblastic leukemia (preB-ALL), we develop an approach to extract perturbed molecular interactions from gene expression changes. We show that the MYC and JunD transcriptional circuits are specifically deregulated after ETV6-RUNX1 and TCF3-PBX1 gene fusions, respectively. We also identified the bulk mRNA NXF1-dependent machinery as a direct target for the TCF3-PBX1 fusion protein. Through a novel approach combining gene expression and interactome data analysis, we provide new insight into TCF3-PBX1 and ETV6-RUNX1 acute lymphoblastic leukemia.

Hox protein interactions: screening and network building

Understanding the mode of action of Hox proteins requires the identification of molecular and cellular pathways they take part in. This includes to characterize the networks of protein-protein interactions involving Hox proteins. In this chapter we propose a strategy and methods to map Hox interaction networks, from yeast two-hybrid and high-throughput yeast two-hybrid interaction screening to bioinformatic analyses based on the software platform Cytoscape.

The homeodomain transcription factor Hoxa2 interacts with and promotes the proteasomal degradation of the E3 ubiquitin protein ligase RCHY1

Hox proteins are conserved homeodomain transcription factors known to be crucial regulators of animal development. As transcription factors, the functions and modes of action (co-factors, target genes) of Hox proteins have been very well studied in a multitude of animal models. However, a handful of reports established that Hox proteins may display molecular activities distinct from gene transcription regulation. Here, we reveal that Hoxa2 interacts with 20S proteasome subunits and RCHY1 (also known as PIRH2), an E3 ubiquitin ligase that targets p53 for degradation. We further show that Hoxa2 promotes proteasome-dependent degradation of RCHY1 in an ubiquitin-independent manner. Correlatively, Hoxa2 alters the RCHY1-mediated ubiquitination of p53 and promotes p53 stabilization. Together, our data establish that Hoxa2 can regulate the proteasomal degradation of RCHY1 and stabilization of p53.

PP2A regulatory subunit Bα controls endothelial contractility and vessel lumen integrity via regulation of HDAC7

To supply tissues with nutrients and oxygen, the cardiovascular system forms a seamless, hierarchically branched, network of lumenized tubes. Here, we show that maintenance of patent vessel lumens requires the Bα regulatory subunit of protein phosphatase 2A (PP2A). Deficiency of Bα in zebrafish precludes vascular lumen stabilization resulting in perfusion defects. Similarly, inactivation of PP2A-Bα in cultured ECs induces tubulogenesis failure due to alteration of cytoskeleton dynamics, actomyosin contractility and maturation of cell-extracellular matrix (ECM) contacts. Mechanistically, we show that PP2A-Bα controls the activity of HDAC7, an essential transcriptional regulator of vascular stability. In the absence of PP2A-Bα, transcriptional repression by HDAC7 is abrogated leading to enhanced expression of the cytoskeleton adaptor protein ArgBP2. ArgBP2 hyperactivates RhoA causing inadequate rearrangements of the EC actomyosin cytoskeleton. This study unravels the first specific role for a PP2A holoenzyme in development: the PP2A-Bα/HDAC7/ArgBP2 axis maintains vascular lumens by balancing endothelial cytoskeletal dynamics and cell-matrix adhesion.

Protein interactions of the transcription factor Hoxa1

Background

Hox proteins are transcription factors involved in crucial processes during animal development. Their mode of action remains scantily documented. While other families of transcription factors, like Smad or Stat, are known cell signaling transducers, such a function has never been squarely addressed for Hox proteins.

Results

To investigate the mode of action of mammalian Hoxa1, we characterized its interactome by a systematic yeast two-hybrid screening against ~12,200 ORF-derived polypeptides. Fifty nine interactors were identified of which 45 could be confirmed by affinity co-purification in animal cell lines. Many Hoxa1 interactors are proteins involved in cell-signaling transduction, cell adhesion and vesicular trafficking. Forty-one interactions were detectable in live cells by Bimolecular Fluorescence Complementation which revealed distinctive intracellular patterns for these interactions consistent with the selective recruitment of Hoxa1 by subgroups of partner proteins at vesicular, cytoplasmic or nuclear compartments.

Conclusions

The characterization of the Hoxa1 interactome presented here suggests unexplored roles for Hox proteins in cell-to-cell communication and cell physiology.

Interaction network of antimicrobial peptides of Arabidopsis thaliana, based on high-throughput yeast two-hybrid screening

One mechanism used by plants to respond to infection is the production of antimicrobial peptides (AMPs). In addition to a role in defence, AMPs seem to have other biological functions. Furthermore, the number of cysteine-rich AMP-like peptides appears to have been underpredicted in plant genomes. Such peptides could be involved in plant defence and/or in other biological processes. Here we generated an interaction network between 15 AMPs/AMP-like peptides and ca. 8000 other Arabidopsis thaliana proteins (AtORFeome2.0) and found 53 putative novel interactions. These interactions involve five transcription factors, a subunit of the COP9 signalosome, a heat shock protein, a MAP kinase kinase, a thioredoxin and 4 uncharacterized proteins.

Host-pathogen interactome mapping for HTLV-1 and -2 retroviruses

BACKGROUND

Human T-cell leukemia virus type 1 (HTLV-1) and type 2 both target T lymphocytes, yet induce radically different phenotypic outcomes. HTLV-1 is a causative agent of Adult T-cell leukemia (ATL), whereas HTLV-2, highly similar to HTLV-1, causes no known overt disease. HTLV gene products are engaged in a dynamic struggle of activating and antagonistic interactions with host cells. Investigations focused on one or a few genes have identified several human factors interacting with HTLV viral proteins. Most of the available interaction data concern the highly investigated HTLV-1 Tax protein. Identifying shared and distinct host-pathogen protein interaction profiles for these two viruses would enlighten how they exploit distinctive or common strategies to subvert cellular pathways toward disease progression.

RESULTS

We employ a scalable methodology for the systematic mapping and comparison of pathogen-host protein interactions that includes stringent yeast two-hybrid screening and systematic retest, as well as two independent validations through an additional protein interaction detection method and a functional transactivation assay. The final data set contained 166 interactions between 10 viral proteins and 122 human proteins. Among the 166 interactions identified, 87 and 79 involved HTLV-1 and HTLV-2 -encoded proteins, respectively. Targets for HTLV-1 and HTLV-2 proteins implicate a diverse set of cellular processes including the ubiquitin-proteasome system, the apoptosis, different cancer pathways and the Notch signaling pathway.

CONCLUSIONS

This study constitutes a first pass, with homogeneous data, at comparative analysis of host targets for HTLV-1 and -2 retroviruses, complements currently existing data for formulation of systems biology models of retroviral induced diseases and presents new insights on biological pathways involved in retroviral infection.

Identification of a Brucella spp. secreted effector specifically interacting with human small GTPase Rab2

Bacteria of the Brucella genus are facultative intracellular class III pathogens. These bacteria are able to control the intracellular trafficking of their vacuole, presumably by the use of yet unknown translocated effectors. To identify such effectors, we used a high-throughput yeast two-hybrid screen to identify interactions between putative human phagosomal proteins and predicted Brucella spp. proteins. We identified a specific interaction between the human small GTPase Rab2 and a Brucella spp. protein named RicA. This interaction was confirmed by GST-pull-down with the GDP-bound form of Rab2. A TEM-β-lactamase-RicA fusion was translocated from Brucella abortus to RAW264.7 macrophages during infection. This translocation was not detectable in a strain deleted for the virB operon, coding for the type IV secretion system. However, RicA secretion in a bacteriological culture was still observed in a ΔvirB mutant. In HeLa cells, a ΔricA mutant recruits less GTP-locked myc-Rab2 on its Brucella-containing vacuoles, compared with the wild-type strain. We observed altered kinetics of intracellular trafficking and faster proliferation of the B. abortusΔricA mutant in HeLa cells, compared with the wild-type control. Altogether, the data reported here suggest RicA as the first reported effector with a proposed function for B. abortus.

The HTLV-1 Tax protein inhibits formation of stress granules by interacting with histone deacetylase 6

Human T cell leukemia virus type-1 (HTLV-1) is the causative agent of a fatal adult T-cell leukemia. Through deregulation of multiple cellular signaling pathways the viral Tax protein has a pivotal role in T-cell transformation. In response to stressful stimuli, cells mount a cellular stress response to limit the damage that environmental forces inflict on DNA or proteins. During stress response, cells postpone the translation of most cellular mRNAs, which are gathered into cytoplasmic mRNA-silencing foci called stress granules (SGs) and allocate their available resources towards the production of dedicated stress-management proteins. Here we demonstrate that Tax controls the formation of SGs and interferes with the cellular stress response pathway. In agreement with previous reports, we observed that Tax relocates from the nucleus to the cytoplasm in response to environmental stress. We found that the presence of Tax in the cytoplasm of stressed cells prevents the formation of SGs and counteracts the shutoff of specific host proteins. Unexpectedly, nuclear localization of Tax promotes spontaneous aggregation of SGs, even in the absence of stress. Mutant analysis revealed that the SG inhibitory capacity of Tax is independent of its transcriptional abilities but relies on its interaction with histone deacetylase 6, a critical component of SGs. Importantly, the stress-protective effect of Tax was also observed in the context of HTLV-1 infected cells, which were shown to be less prone to form SGs and undergo apoptosis under arsenite exposure. These observations identify Tax as the first virally encoded inhibitory component of SGs and unravel a new strategy developed by HTLV-1 to deregulate normal cell processes. We postulate that inhibition of the stress response pathway by Tax would favor cell survival under stressful conditions and may have an important role in HTLV-1-induced cellular transformation.

DNA cytosine methylation in the bovine leukemia virus promoter is associated with latency in a lymphoma-derived B-cell line: potential involvement of direct inhibition of cAMP-responsive element (CRE)-binding protein/CRE modulator/activation transcription factor binding

Bovine leukemia virus (BLV) proviral latency represents a viral strategy to escape the host immune system and allow tumor development. Besides the previously demonstrated role of histone deacetylation in the epigenetic repression of BLV expression, we showed here that BLV promoter activity was induced by several DNA methylation inhibitors (such as 5-aza-2'-deoxycytidine) and that overexpressed DNMT1 and DNMT3A, but not DNMT3B, down-regulated BLV promoter activity. Importantly, cytosine hypermethylation in the 5'-long terminal repeat (LTR) U3 and R regions was associated with true latency in the lymphoma-derived B-cell line L267 but not with defective latency in YR2 cells. Moreover, the virus-encoded transactivator Tax(BLV) decreased DNA methyltransferase expression levels, which could explain the lower level of cytosine methylation observed in the L267(LTaxSN) 5'-LTR compared with the L267 5'-LTR. Interestingly, DNA methylation inhibitors and Tax(BLV) synergistically activated BLV promoter transcriptional activity in a cAMP-responsive element (CRE)-dependent manner. Mechanistically, methylation at the -154 or -129 CpG position (relative to the transcription start site) impaired in vitro binding of CRE-binding protein (CREB) transcription factors to their respective CRE sites. Methylation at -129 CpG alone was sufficient to decrease BLV promoter-driven reporter gene expression by 2-fold. We demonstrated in vivo the recruitment of CREB/CRE modulator (CREM) and to a lesser extent activating transcription factor-1 (ATF-1) to the hypomethylated CRE region of the YR2 5'-LTR, whereas we detected no CREB/CREM/ATF recruitment to the hypermethylated corresponding region in the L267 cells. Altogether, these findings suggest that site-specific DNA methylation of the BLV promoter represses viral transcription by directly inhibiting transcription factor binding, thereby contributing to true proviral latency.