Multimodal bioinformatic analyses of the neurodegenerative disease-associated TECPR2 gene reveal its diverse roles

BackgroundLoss of tectonin β-propeller repeat-containing 2 (TECPR2) function has been implicated in an array of neurodegenerative disorders, yet its physiological function remains largely unknown. Understanding TECPR2 function is essential for developing much needed precision therapeutics for TECPR2-related diseases.MethodsWe leveraged considerable amounts of functional data to obtain a comprehensive perspective of the role of TECPR2 in health and disease. We integrated expression patterns, population variation, phylogenetic profiling, protein-protein interactions and regulatory network data for a minimally biased multimodal functional analysis. Genes and proteins linked to TECPR2 via multiple lines of evidence were subject to functional enrichment analyses to identify molecular mechanisms involving TECPR2.ResultsTECPR2 was found to be part of a tight neurodevelopmental gene expression programme that includes KIF1A, ATXN1, TOM1L2 and FA2H, all implicated in neurological diseases. Functional enrichment analyses of TECPR2-related genes converged on a role in late autophagy and ribosomal processes. Large-scale population variation data demonstrated that this role is non-redundant.ConclusionsTECPR2 might serve as an indicator for the energy balance between protein synthesis and autophagy, and a marker for diseases associated with their imbalance, such as Alzheimer’s disease and Huntington’s disease. Specifically, we speculate that TECPR2 plays an important role as a proteostasis regulator during synaptogenesis, highlighting its importance in developing neurons. By advancing our understanding of TECPR2 function, this work provides an essential stepping stone towards the development of precision diagnostics and targeted treatment options for TECPR2-related disorders.

GFICLEE: ultrafast tree-based phylogenetic profile method inferring gene function at the genomic-wide level

Background Phylogenetic profiling is widely used to predict novel members of large protein complexes and biological pathways. Although methods combined with phylogenetic trees have significantly improved prediction accuracy, computational efficiency is still an issue that limits its genome-wise application. Results Here we introduce a new tree-based phylogenetic profiling algorithm named GFICLEE, which infers common single and continuous loss (SCL) events in the evolutionary patterns. We validated our algorithm with human pathways from three databases and compared the computational efficiency with current tree-based with 10 different scales genome dataset. Our algorithm has a better predictive performance with high computational efficiency. Conclusions The GFICLEE is a new method to infers genome-wide gene function. The accuracy and computational efficiency of GFICLEE make it possible to explore gene functions at the genome-wide level on a personal computer.

Inverse Potts model improves accuracy of phylogenetic profiling

Motivation: Phylogenetic profiling is a powerful computational method for revealing the functions of function-unknown genes. Although conventional similarity evaluation measures in phylogenetic profiling showed high prediction accuracy, they have two estimation biases: an evolutionary bias and a spurious correlation bias. Existing studies have focused on the evolutionary bias, but the spurious correlation bias has not been analyzed. Results: To eliminate the spurious correlation bias, we applied an evaluation measure based on the inverse Potts model (IPM) to phylogenetic profiling. We also proposed an evaluation measure to remove both the evolutionary and spurious correlation biases using the IPM. In an empirical dataset analysis, we demonstrated that these IPM-based evaluation measures improved the prediction performance of phylogenetic profiling. In addition, we found that the integration of several evaluation measures, including the IPM-based evaluation measures, had superior performance to a single evaluation measure.

Phylogenetic profiling and cellular analyses of ARL16 reveal roles in traffic of IFT140 and INPP5E

The ARF family of regulatory GTPases is ancient, with 16 members predicted to have been present in the last eukaryotic common ancestor. Our phylogenetic profiling of paralogs in diverse species identified four family members whose presence correlates with that of a cilium/flagellum: ARL3, ARL6, ARL13, and ARL16. No prior evidence links ARL16 to cilia or other cell functions, despite its presence throughout eukaryotes. Deletion of ARL16 in MEFs results in decreased ciliogenesis yet increased ciliary length. We also found Arl16 KO in MEFs to alter ciliary protein content, including loss of ARL13B, ARL3, INPP5E, and the IFT-A core component IFT140. Instead, both INPP5E and IFT140 accumulate at the Golgi in Arl16 KO lines, while other IFT proteins do not, suggesting a specific defect in traffic from Golgi to cilia. We propose that ARL16 regulates a Golgi-cilia traffic pathway and is required specifically in the export of IFT140 and INPP5E from the Golgi.

Phylogenetic profiling suggests early origin of the core subunits of Polycomb Repressive Complex 2 (PRC2)

Polycomb Repressive Complex 2 (PRC2) is involved in establishing transcriptionally silent chromatin states through its ability to methylate lysine 27 of histone H3 by the catalytic subunit Enhancer of zeste [E(z)]. Polycomb group (PcG) proteins play a crucial role in the maintenance of cell identity and in developmental regulation. Previously, the diversity of PRC2 subunits within some eukaryotic lineages has been reported and its presence in early eukaryotic evolution has been hypothesized. So far however, systematic survey of the presence of PRC2 subunits in species of all eukaryotic lineages is missing. Here, we report the diversity of PRC2 core subunit proteins in different eukaryotic supergroups with emphasis on the early-diverged lineages and explore the molecular evolution of PRC2 subunits by phylogenetics. In detail, we investigate the SET-domain protein sequences and their evolution across the four domains of life and particularly focus on the structural diversity of the SET-domain subfamily containing E(z), the catalytic subunit of PRC2. We show that PRC2 subunits are already present in early eukaryotic lineages, strengthening the support for PRC2 emergence prior to diversification of eukaryotes. We identify a common presence of E(z) and ESC, suggesting that Su(z)12 may have emerged later and/or may be dispensable from the evolutionarily conserved functional core of PRC2. Furthermore, our results broaden our understanding of the E(z) evolution within the SET-domain protein family, suggesting possibilities of function evolution. Through this, we shed light on a possible emerging point of the PRC2 and the evolution of its function in eukaryotes.

Integrating taxonomic, functional, and strain-level profiling of diverse microbial communities with bioBakery 3

Culture-independent analyses of microbial communities have progressed dramatically in the last decade, particularly due to advances in methods for biological profiling via shotgun metagenomics. Opportunities for improvement continue to accelerate, with greater access to multi-omics, microbial reference genomes, and strain-level diversity. To leverage these, we present bioBakery 3, a set of integrated, improved methods for taxonomic, strain-level, functional, and phylogenetic profiling of metagenomes newly developed to build on the largest set of reference sequences now available. Compared to current alternatives, MetaPhlAn 3 increases the accuracy of taxonomic profiling, and HUMAnN 3 improves that of functional potential and activity. These methods detected novel disease-microbiome links in applications to CRC (1262 metagenomes) and IBD (1635 metagenomes and 817 metatranscriptomes). Strain-level profiling of an additional 4077 metagenomes with StrainPhlAn 3 and PanPhlAn 3 unraveled the phylogenetic and functional structure of the common gut microbe Ruminococcus bromii, previously described by only 15 isolate genomes. With open-source implementations and cloud-deployable reproducible workflows, the bioBakery 3 platform can help researchers deepen the resolution, scale, and accuracy of multi-omic profiling for microbial community studies.

CladeOScope: functional interactions through the prism of clade-wise co-evolution

Mapping co-evolved genes via phylogenetic profiling (PP) is a powerful approach to uncover functional interactions between genes and to associate them with pathways. Despite many successful endeavors, the understanding of co-evolutionary signals in eukaryotes remains partial. Our hypothesis is that ‘Clades’, branches of the tree of life (e.g. primates and mammals), encompass signals that cannot be detected by PP using all eukaryotes. As such, integrating information from different clades should reveal local co-evolution signals and improve function prediction. Accordingly, we analyzed 1028 genomes in 66 clades and demonstrated that the co-evolutionary signal was scattered across clades. We showed that functionally related genes are frequently co-evolved in only parts of the eukaryotic tree and that clades are complementary in detecting functional interactions within pathways. We examined the non-homologous end joining pathway and the UFM1 ubiquitin-like protein pathway and showed that both demonstrated distinguished co-evolution patterns in specific clades. Our research offers a different way to look at co-evolution across eukaryotes and points to the importance of modular co-evolution analysis. We developed the ‘CladeOScope’ PP method to integrate information from 16 clades across over 1000 eukaryotic genomes and is accessible via an easy to use web server at http://cladeoscope.cs.huji.ac.il.

The Functional Classification of ORF8 in SARS-CoV-2 Replication, Immune Evasion, and Viral Pathogenesis Inferred through Phylogenetic Profiling

ORF8 is a highly variable genomic region of SARS-CoV-2. Although non-essential and the precise functions are unknown, it has been suggested that this protein assists in SARS-CoV-2 replication in the early secretory pathway and in immune evasion. We utilized the binding partners of SARS-CoV-2 proteins in human HEK293T cells and performed genome-wide phylogenetic profiling and clustering analyses in 446 eukaryotic species to predict and discover ORF8 binding partners that share associated functional mechanisms based on co-evolution. Results classified 47 ORF8 binding partner proteins into 3 clusters (groups 1-3), which were conserved in vertebrates (group 1), metazoan (group 2), and eukaryotes (group 3). Gene ontology analysis indicated that group 1 had no significant associated biological processes, while groups 2 and 3 were associated with glycoprotein biosynthesis process and ubiquitin-dependent endoplasmic reticulum-associated degradation pathways, respectively. Collectively, our results classified potential genes that might be associated with SARS-CoV-2 viral pathogenesis, specifically related to acute respiratory distress syndrome, and the secretory pathway. Here, we discuss the possible role of ORF8 in viral pathogenesis and in assisting viral replication and immune evasion via secretory pathway, as well as the possible factors associated with the rapid evolution of ORF8.