Principles of chromosome organization for meiotic recombination.

In meiotic cells, chromosomes are organized as chromatin loop arrays anchored to a protein axis. This organization is essential to regulate meiotic recombination, from DNA double-strand break (DSB) formation to their repair. In mammals, it is unknown how chromatin loops are organized along the genome and how proteins participating in DSB formation are tethered to the chromosome axes. Here, we identify three categories of axis-associated genomic sites: PRDM9 binding sites, where DSBs form; binding sites of the insulator protein CTCF; and H3K4me3-enriched sites. We demonstrate that PRDM9 promotes the recruitment of MEI4 and IHO1, two proteins essential for DSB formation. In turn, IHO1 anchors DSB sites to the axis components HORMAD1 and SYCP3. We discovered that IHO1, HORMAD1, and SYCP3 are associated at the DSB ends during DSB repair. Our results highlight how interactions of proteins with specific genomic elements shape the meiotic chromosome organization for recombination.

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.

Genetic architecture of natural variation of cardiac performance from flies to humans.

Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP). Genome-wide associations studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict transcription factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart-specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans.

TOPOVIBL-REC114 interaction regulates meiotic DNA double-strand breaks.

Meiosis requires the formation of programmed DNA double strand breaks (DSBs), essential for fertility and for generating genetic diversity. DSBs are induced by the catalytic activity of the TOPOVIL complex formed by SPO11 and TOPOVIBL. To ensure genomic integrity, DNA cleavage activity is tightly regulated, and several accessory factors (REC114, MEI4, IHO1, and MEI1) are needed for DSB formation in mice. How and when these proteins act is not understood. Here, we show that REC114 is a direct partner of TOPOVIBL, and identify their conserved interacting domains by structural analysis. We then analyse the role of this interaction by monitoring meiotic DSBs in female and male mice carrying point mutations in TOPOVIBL that decrease or disrupt its binding to REC114. In these mutants, DSB activity is strongly reduced genome-wide in oocytes, and only in sub-telomeric regions in spermatocytes. In addition, in mutant spermatocytes, DSB activity is delayed in autosomes. These results suggest that REC114 is a key member of the TOPOVIL catalytic complex, and that the REC114/TOPOVIBL interaction ensures the efficiency and timing of DSB activity.

Epi-MEIF: detecting higher order epistatic interactions for complex traits using mixed effect conditional inference forests.

Understanding the relationship between genetic variations and variations in complex and quantitative phenotypes remains an ongoing challenge. While Genome-wide association studies (GWAS) have become a vital tool for identifying single-locus associations, we lack methods for identifying epistatic interactions. In this article, we propose a novel method for higher-order epistasis detection using mixed effect conditional inference forest (epiMEIF). The proposed method is fitted on a group of single nucleotide polymorphisms (SNPs) potentially associated with the phenotype and the tree structure in the forest facilitates the identification of n-way interactions between the SNPs. Additional testing strategies further improve the robustness of the method. We demonstrate its ability to detect true n-way interactions via extensive simulations in both cross-sectional and longitudinal synthetic datasets. This is further illustrated in an application to reveal epistatic interactions from natural variations of cardiac traits in flies (Drosophila). Overall, the method provides a generalized way to identify higher-order interactions from any GWAS data, thereby greatly improving the detection of the genetic architecture underlying complex phenotypes.

In Depth Exploration of the Alternative Proteome of Drosophila melanogaster.

Recent studies have shown that hundreds of small proteins were occulted when protein-coding genes were annotated. These proteins, called alternative proteins, have failed to be annotated notably due to the short length of their open reading frame (less than 100 codons) or the enforced rule establishing that messenger RNAs (mRNAs) are monocistronic. Several alternative proteins were shown to be biologically active molecules and seem to be involved in a wide range of biological functions. However, genome-wide exploration of the alternative proteome is still limited to a few species. In the present article, we describe a deep peptidomics workflow which enabled the identification of 401 alternative proteins in Drosophila melanogaster. Subcellular localization, protein domains, and short linear motifs were predicted for 235 of the alternative proteins identified and point toward specific functions of these small proteins. Several alternative proteins had approximated abundances higher than their canonical counterparts, suggesting that these alternative proteins are actually the main products of their corresponding genes. Finally, we observed 14 alternative proteins with developmentally regulated expression patterns and 10 induced upon the heat-shock treatment of embryos, demonstrating stage or stress-specific production of alternative proteins.

The Intricacy of the Viral-Human Protein Interaction Networks: Resources, Data, and Analyses.

Viral infections are one of the major causes of human diseases that cause yearly millions of deaths and seriously threaten global health, as we have experienced with the COVID-19 pandemic. Numerous approaches have been adopted to understand viral diseases and develop pharmacological treatments. Among them, the study of virus-host protein-protein interactions is a powerful strategy to comprehend the molecular mechanisms employed by the virus to infect the host cells and to interact with their components. Experimental protein-protein interactions described in the scientific literature have been systematically captured into several molecular interaction databases. These data are organized in structured formats and can be easily downloaded by users to perform further bioinformatic and network studies. Network analysis of available virus-host interactomes allow us to understand how the host interactome is perturbed upon viral infection and what are the key host proteins targeted by the virus and the main cellular pathways that are subverted. In this review, we give an overview of publicly available viral-human protein-protein interactions resources and the community standards, curation rules and adopted ontologies. A description of the main virus-human interactome available is provided, together with the main network analyses that have been performed. We finally discuss the main limitations and future challenges to assess the quality and reliability of protein-protein interaction datasets and resources.

Genetic Predisposition to the Mortality in Septic Shock Patients: From GWAS to the Identification of a Regulatory Variant Modulating the Activity of a CISH Enhancer.

The high mortality rate in septic shock patients is likely due to environmental and genetic factors, which influence the host response to infection. Two genome-wide association studies (GWAS) on 832 septic shock patients were performed. We used integrative bioinformatic approaches to annotate and prioritize the sepsis-associated single nucleotide polymorphisms (SNPs). An association of 139 SNPs with death based on a false discovery rate of 5% was detected. The most significant SNPs were within the CISH gene involved in cytokine regulation. Among the 139 SNPs associated with death and the 1311 SNPs in strong linkage disequilibrium with them, we investigated 1439 SNPs within non-coding regions to identify regulatory variants. The highest integrative weighted score (IW-score) was obtained for rs143356980, indicating that this SNP is a robust regulatory candidate. The rs143356980 region is located in a non-coding region close to the CISH gene. A CRISPR-Cas9-mediated deletion of this region and specific luciferase assays in K562 cells showed that rs143356980 modulates the enhancer activity in K562 cells. These analyses allowed us to identify several genes associated with death in patients with septic shock. They suggest that genetic variations in key genes, such as CISH, perturb relevant pathways, increasing the risk of death in sepsis patients.

MetamORF: a repository of unique short open reading frames identified by both experimental and computational approaches for gene and metagene analyses.

The development of high-throughput technologies revealed the existence of non-canonical short open reading frames (sORFs) on most eukaryotic ribonucleic acids. They are ubiquitous genetic elements conserved across species and suspected to be involved in numerous cellular processes. MetamORF (https://metamorf.hb.univ-amu.fr/) aims to provide a repository of unique sORFs identified in the human and mouse genomes with both experimental and computational approaches. By gathering publicly available sORF data, normalizing them and summarizing redundant information, we were able to identify a total of 1 162 675 unique sORFs. Despite the usual characterization of ORFs as short, upstream or downstream, there is currently no clear consensus regarding the definition of these categories. Thus, the data have been reprocessed using a normalized nomenclature. MetamORF enables new analyses at locus, gene, transcript and ORF levels, which should offer the possibility to address new questions regarding sORF functions in the future. The repository is available through an user-friendly web interface, allowing easy browsing, visualization, filtering over multiple criteria and export possibilities. sORFs can be searched starting from a gene, a transcript and an ORF ID, looking in a genome area or browsing the whole repository for a species. The database content has also been made available through track hubs at UCSC Genome Browser. Finally, we demonstrated an enrichment of genes harboring upstream ORFs among genes expressed in response to reticular stress. Database URL  https://metamorf.hb.univ-amu.fr/.

The role of 3'UTR-protein complexes in the regulation of protein multifunctionality and subcellular localization.

Multifunctional proteins often perform their different functions when localized in different subcellular compartments. However, the mechanisms leading to their localization are largely unknown. Recently, 3'UTRs were found to regulate the cellular localization of newly synthesized proteins through the formation of 3'UTR-protein complexes. Here, we investigate the formation of 3'UTR-protein complexes involving multifunctional proteins by exploiting large-scale protein-protein and protein-RNA interaction networks. Focusing on 238 human 'extreme multifunctional' (EMF) proteins, we predicted 1411 3'UTR-protein complexes involving 54% of those proteins and evaluated their role in regulating protein cellular localization and multifunctionality. We find that EMF proteins lacking localization addressing signals, yet present at both the nucleus and cell surface, often form 3'UTR-protein complexes, and that the formation of these complexes could provide EMF proteins with the diversity of interaction partners necessary to their multifunctionality. Our findings are reinforced by archetypal moonlighting proteins predicted to form 3'UTR-protein complexes. Finally, the formation of 3'UTR-protein complexes that involves up to 17% of the proteins in the human protein-protein interaction network, may be a common and yet underestimated protein trafficking mechanism, particularly suited to regulate the localization of multifunctional proteins.

Understanding protein multifunctionality: from short linear motifs to cellular functions.

Moonlighting proteins perform multiple unrelated functions without any change in polypeptide sequence. They can coordinate cellular activities, serving as switches between pathways and helping to respond to changes in the cellular environment. Therefore, regulation of the multiple protein activities, in space and time, is likely to be important for the homeostasis of biological systems. Some moonlighting proteins may perform their multiple functions simultaneously while others alternate between functions due to certain triggers. The switch of the moonlighting protein's functions can be regulated by several distinct factors, including the binding of other molecules such as proteins. We here review the approaches used to identify moonlighting proteins and existing repositories. We particularly emphasise the role played by short linear motifs and PTMs as regulatory switches of moonlighting functions.

In utero exposure to acetaminophen and ibuprofen leads to intergenerational accelerated reproductive aging in female mice.

Nonsteroidal anti-inflammatory drugs (NSAIDs) and analgesic drugs, such as acetaminophen (APAP), are frequently taken during pregnancy, even in combination. However, they can favour genital malformations in newborn boys and reproductive disorders in adults. Conversely, the consequences on postnatal ovarian development and female reproductive health after in utero exposure are unknown. Here, we found that in mice, in utero exposure to therapeutic doses of the APAP-ibuprofen combination during sex determination led to delayed meiosis entry and progression in female F1 embryonic germ cells. Consequently, follicular activation was reduced in postnatal ovaries through the AKT/FOXO3 pathway, leading in F2 animals to subfertility, accelerated ovarian aging with abnormal corpus luteum persistence, due to decreased apoptosis and increased AKT-mediated luteal cell survival. Our study suggests that administration of these drugs during the critical period of sex determination could lead in humans to adverse effects that might be passed to the offspring.

Post-transcriptional regulatory patterns revealed by protein-RNA interactions.

The coordination of the synthesis of functionally-related proteins can be achieved at the post-transcriptional level by the action of common regulatory molecules, such as RNA-binding proteins (RBPs). Despite advances in the genome-wide identification of RBPs and their binding transcripts, the protein-RNA interaction space is still largely unexplored, thus hindering a broader understanding of the extent of the post-transcriptional regulation of related coding RNAs. Here, we propose a computational approach that combines protein-mRNA interaction networks and statistical analyses to provide an inferred regulatory landscape for more than 800 human RBPs and identify the cellular processes that can be regulated at the post-transcriptional level. We show that 10% of the tested sets of functionally-related mRNAs can be post-transcriptionally regulated. Moreover, we propose a classification of (i) the RBPs and (ii) the functionally-related mRNAs, based on their distinct behaviors in the functional landscape, hinting towards mechanistic regulatory hypotheses. In addition, we demonstrate the usefulness of the inferred functional landscape to investigate the cellular role of both well-characterized and novel RBPs in the context of human diseases.

MoonDB 2.0: an updated database of extreme multifunctional and moonlighting proteins.

MoonDB 2.0 (http://moondb.hb.univ-amu.fr/) is a database of predicted and manually curated extreme multifunctional (EMF) and moonlighting proteins, i.e. proteins that perform multiple unrelated functions. We have previously shown that such proteins can be predicted through the analysis of their molecular interaction subnetworks, their functional annotations and their association to distinct groups of proteins that are involved in unrelated functions. In MoonDB 2.0, we updated the set of human EMF proteins (238 proteins), using the latest functional annotations and protein-protein interaction networks. Furthermore, for the first time, we applied our method to four additional model organisms - mouse, fly, worm and yeast - and identified 54 novel EMF proteins in these species. In addition to novel predictions, this update contains 63 human and yeast proteins that were manually curated from literature, including descriptions of moonlighting functions and associated references. Importantly, MoonDB's interface was fully redesigned and improved, and its entries are now cross-referenced in the UniProt Knowledgebase (UniProtKB). MoonDB will be updated once a year with the novel EMF candidates calculated from the latest available protein interactions and functional annotations.

Mouse REC114 is essential for meiotic DNA double-strand break formation and forms a complex with MEI4.

Programmed formation of DNA double-strand breaks (DSBs) initiates the meiotic homologous recombination pathway. This pathway is essential for proper chromosome segregation at the first meiotic division and fertility. Meiotic DSBs are catalyzed by Spo11. Several other proteins are essential for meiotic DSB formation, including three evolutionarily conserved proteins first identified in Saccharomyces cerevisiae (Mer2, Mei4, and Rec114). These three S. cerevisiae proteins and their mouse orthologs (IHO1, MEI4, and REC114) co-localize on the axes of meiotic chromosomes, and mouse IHO1 and MEI4 are essential for meiotic DSB formation. Here, we show that mouse Rec114 is required for meiotic DSB formation. Moreover, MEI4 forms a complex with REC114 and IHO1 in mouse spermatocytes, consistent with cytological observations. We then demonstrated in vitro the formation of a stable complex between REC114 C-terminal domain and MEI4 N-terminal domain. We further determine the structure of the REC114 N-terminal domain that revealed similarity with Pleckstrin homology domains. These analyses provide direct insights into the architecture of these essential components of the meiotic DSB machinery.

Generation of a versatile BiFC ORFeome library for analyzing protein-protein interactions in live Drosophila.

Transcription factors achieve specificity by establishing intricate interaction networks that will change depending on the cell context. Capturing these interactions in live condition is however a challenging issue that requires sensitive and non-invasive methods.We present a set of fly lines, called 'multicolor BiFC library', which covers most of the Drosophila transcription factors for performing Bimolecular Fluorescence Complementation (BiFC). The multicolor BiFC library can be used to probe two different binary interactions simultaneously and is compatible for large-scale interaction screens. The library can also be coupled with established Drosophila genetic resources to analyze interactions in the developmentally relevant expression domain of each protein partner. We provide proof of principle experiments of these various applications, using Hox proteins in the live Drosophila embryo as a case study. Overall this novel collection of ready-to-use fly lines constitutes an unprecedented genetic toolbox for the identification and analysis of protein-protein interactions in vivo.

Protein complex scaffolding predicted as a prevalent function of long non-coding RNAs.

The human transcriptome contains thousands of long non-coding RNAs (lncRNAs). Characterizing their function is a current challenge. An emerging concept is that lncRNAs serve as protein scaffolds, forming ribonucleoproteins and bringing proteins in proximity. However, only few scaffolding lncRNAs have been characterized and the prevalence of this function is unknown. Here, we propose the first computational approach aimed at predicting scaffolding lncRNAs at large scale. We predicted the largest human lncRNA-protein interaction network to date using the catRAPID omics algorithm. In combination with tissue expression and statistical approaches, we identified 847 lncRNAs (∼5% of the long non-coding transcriptome) predicted to scaffold half of the known protein complexes and network modules. Lastly, we show that the association of certain lncRNAs to disease may involve their scaffolding ability. Overall, our results suggest for the first time that RNA-mediated scaffolding of protein complexes and modules may be a common mechanism in human cells.

Perturbed human sub-networks by Fusobacterium nucleatum candidate virulence proteins.

BACKGROUND: Fusobacterium nucleatum is a gram-negative anaerobic species residing in the oral cavity and implicated in several inflammatory processes in the human body. Although F. nucleatum abundance is increased in inflammatory bowel disease subjects and is prevalent in colorectal cancer patients, the causal role of the bacterium in gastrointestinal disorders and the mechanistic details of host cell functions subversion are not fully understood. RESULTS: We devised a computational strategy to identify putative secreted F. nucleatum proteins (FusoSecretome) and to infer their interactions with human proteins based on the presence of host molecular mimicry elements. FusoSecretome proteins share similar features with known bacterial virulence factors thereby highlighting their pathogenic potential. We show that they interact with human proteins that participate in infection-related cellular processes and localize in established cellular districts of the host-pathogen interface. Our network-based analysis identified 31 functional modules in the human interactome preferentially targeted by 138 FusoSecretome proteins, among which we selected 26 as main candidate virulence proteins, representing both putative and known virulence proteins. Finally, six of the preferentially targeted functional modules are implicated in the onset and progression of inflammatory bowel diseases and colorectal cancer. CONCLUSIONS: Overall, our computational analysis identified candidate virulence proteins potentially involved in the F. nucleatum-human cross-talk in the context of gastrointestinal diseases.

Cancer-associated fibroblast-derived annexin A6+ extracellular vesicles support pancreatic cancer aggressiveness.

The intratumoral microenvironment, or stroma, is of major importance in the pathobiology of pancreatic ductal adenocarcinoma (PDA), and specific conditions in the stroma may promote increased cancer aggressiveness. We hypothesized that this heterogeneous and evolving compartment drastically influences tumor cell abilities, which in turn influences PDA aggressiveness through crosstalk that is mediated by extracellular vesicles (EVs). Here, we have analyzed the PDA proteomic stromal signature and identified a contribution of the annexin A6/LDL receptor-related protein 1/thrombospondin 1 (ANXA6/LRP1/TSP1) complex in tumor cell crosstalk. Formation of the ANXA6/LRP1/TSP1 complex was restricted to cancer-associated fibroblasts (CAFs) and required physiopathologic culture conditions that improved tumor cell survival and migration. Increased PDA aggressiveness was dependent on tumor cell-mediated uptake of CAF-derived ANXA6+ EVs carrying the ANXA6/LRP1/TSP1 complex. Depletion of ANXA6 in CAFs impaired complex formation and subsequently impaired PDA and metastasis occurrence, while injection of CAF-derived ANXA6+ EVs enhanced tumorigenesis. We found that the presence of ANXA6+ EVs in serum was restricted to PDA patients and represents a potential biomarker for PDA grade. These findings suggest that CAF-tumor cell crosstalk supported by ANXA6+ EVs is predictive of PDA aggressiveness, highlighting a therapeutic target and potential biomarker for PDA.

Integration of quantitative proteomics data and interaction networks: Identification of dysregulated cellular functions during cancer progression.

Quantitative proteomics allows the characterization of molecular changes between healthy and disease states. To interpret such datasets, their integration to the protein-protein interaction network provides a more comprehensive understanding of cellular function dysregulation in diseases than just considering lists of dysregulated proteins. Here, we propose a novel computational method, which combines protein interaction network and statistical analyses to establish expression profiles at the network module level rather than at the individual protein level, and to detect and characterize dysregulated network modules through different stages of cancer progression. We applied our approach to two publicly available datasets as case studies.