DGINN, an automated and highly-flexible pipeline for the detection of genetic innovations on protein-coding genes

AOC: Analysis of Orthologous Collections - an application for the characterization of natural selection in protein-coding sequences

Motivation

Modern molecular sequence analysis increasingly relies on automated and robust software tools for interpretation, annotation, and biological insight. The Analysis of Orthologous Collections (AOC) application automates the identification of genomic sites and species/lineages influenced by natural selection in coding sequence analysis. AOC quantifies different types of selection: negative, diversifying or directional positive, or differential selection between groups of branches. We include all steps necessary to go from unaligned homologous sequences to complete results and interactive visualizations that are designed to aid in the useful interpretation and contextualization.

Results

We are motivated by a desire to make evolutionary analyses as simple as possible, and to close the disparity in the literature between genes which draw a significant amount of interest and those that are largely overlooked and underexplored. We believe that such underappreciated and understudied genetic datasets can hold rich biological information and offer substantial insights into the diverse patterns and processes of evolution, especially if domain experts are able to perform the analyses themselves.

Availability and implementation

A Snakemake [Mölder et al., 2021] application implementation is publicly available on GitHub at https://github.com/aglucaci/AnalysisOfOrthologousCollections and is accompanied by software documentation and a tutorial.

Evolutionary immunology to explore original antiviral strategies

Over the past 25 years, the field of evolutionary developmental biology (evo-devo) has used genomics and genetics to gain insight on the developmental mechanisms underlying the evolution of morphological diversity of animals. Evo-devo exploits the key insight that conserved toolkits of development (e.g. Hox genes) are used in animals to produce genetic novelties that provide adaptation to a new environment. Like development, immunity is forged by interactions with the environment, namely the microbial world. Yet, when it comes to the study of immune defence mechanisms in invertebrates, interest primarily focuses on evolutionarily conserved molecules also present in humans. Here, focusing on antiviral immunity, we argue that immune genes not conserved in humans represent an unexplored resource for the discovery of new antiviral strategies. We review recent findings on the cGAS-STING pathway and explain how cyclic dinucleotides produced by cGAS-like receptors may be used to investigate the portfolio of antiviral genes in a broad range of species. This will set the stage for evo-immuno approaches, exploiting the investment in antiviral defences made by metazoans over hundreds of millions of years of evolution. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

Recognition and Cleavage of Human tRNA Methyltransferase TRMT1 by the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) is critical for the production of functional viral proteins during infection and, like many viral proteases, can also target host proteins to subvert their cellular functions. Here, we show that the human tRNA methyltransferase TRMT1 can be recognized and cleaved by SARS-CoV-2 Mpro. TRMT1 installs the N2,N2-dimethylguanosine (m2,2G) modification on mammalian tRNAs, which promotes global protein synthesis and cellular redox homeostasis. We find that Mpro can cleave endogenous TRMT1 in human cell lysate, resulting in removal of the TRMT1 zinc finger domain required for tRNA modification activity in cells. Evolutionary analysis shows that the TRMT1 cleavage site is highly conserved in mammals, except in Muroidea, where TRMT1 may be resistant to cleavage. In primates, regions outside the cleavage site with rapid evolution could indicate adaptation to ancient viral pathogens. We determined the structure of a TRMT1 peptide in complex with Mpro, revealing a substrate binding conformation distinct from the majority of available Mpro-peptide complexes. Kinetic parameters for peptide cleavage showed that the TRMT1(526-536) sequence is cleaved with comparable efficiency to the Mpro-targeted nsp8/9 viral cleavage site. Mutagenesis studies and molecular dynamics simulations together indicate that kinetic discrimination occurs during a later step of Mpro-mediated proteolysis that follows substrate binding. Our results provide new information about the structural basis for Mpro substrate recognition and cleavage that could help inform future therapeutic design and raise the possibility that proteolysis of human TRMT1 during SARS-CoV-2 infection suppresses protein translation and oxidative stress response to impact viral pathogenesis.

FREEDA: An automated computational pipeline guides experimental testing of protein innovation

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that lead to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA that provides a simple graphical user interface requiring only a gene name; integrates widely used molecular evolution tools to detect positive selection in rodents, primates, carnivores, birds, and flies; and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 centromere proteins, we find statistical evidence of positive selection within loops and turns of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of mouse CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

FREEDA: an automated computational pipeline guides experimental testing of protein innovation by detecting positive selection

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that leads to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA (Finder of Rapidly Evolving Exons in De novo Assemblies) that provides a simple graphical user interface requiring only a gene name, integrates widely used molecular evolution tools to detect positive selection, and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 mouse centromere proteins, we find evidence of positive selection in intrinsically disordered regions of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

Comparative study of GBP recruitment on two cytosol-dwelling pathogens, Francisella novicida and Shigella flexneri highlights differences in GBP repertoire and in GBP1 motif requirements

Guanylate-Binding Proteins are interferon-inducible GTPases that play a key role in cell autonomous responses against intracellular pathogens. Despite sharing high sequence similarity, subtle differences among GBPs translate into functional divergences that are still largely not understood. A key GBP feature is the formation of supramolecular GBP complexes on the bacterial surface. Such complexes are observed when GBP1 binds lipopolysaccharide (LPS) from Shigella and Salmonella and further recruits GBP2-4. Here, we compared GBP recruitment on two cytosol-dwelling pathogens, Francisella novicida and S. flexneri. Francisella novicida was coated by GBP1 and GBP2 and to a lower extent by GBP4 in human macrophages. Contrary to S. flexneri, F. novicida was not targeted by GBP3, a feature independent of T6SS effectors. Multiple GBP1 features were required to promote targeting to F. novicida while GBP1 targeting to S. flexneri was much more permissive to GBP1 mutagenesis suggesting that GBP1 has multiple domains that cooperate to recognize F. novicida atypical LPS. Altogether our results indicate that the repertoire of GBPs recruited onto specific bacteria is dictated by GBP-specific features and by specific bacterial factors that remain to be identified.

Trim69 is a microtubule regulator that acts as a pantropic viral inhibitor

Through a screen that combines functional and evolutionary analyses, we identified tripartite motif protein (Trim69), a poorly studied member of the Trim family, as a negative regulator of HIV-1 infection in interferon (IFN)-stimulated myeloid cells. Trim69 inhibits the early phases of infection of HIV-1, but also of HIV-2 and SIV_MAC in addition to the negative and positive-strand RNA viruses vesicular stomatitis virus and severe acute respiratory syndrome coronavirus 2, with magnitudes that depend on the combination between cell type and virus. Mechanistically, Trim69 associates directly to microtubules and its antiviral activity is linked to its ability to promote the accumulation of stable microtubules, a program that we uncover to be an integral part of antiviral IFN-I responses in myeloid cells. Overall, our study identifies Trim69 as the antiviral innate defense factor that regulates the properties of microtubules to limit viral spread and highlights the cytoskeleton as an unappreciated battleground in the host-pathogen interactions that underlie viral infections.

Distinct evolutionary trajectories of SARS-CoV-2-interacting proteins in bats and primates identify important host determinants of COVID-19

The coronavirus disease 19 (COVID-19) pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a coronavirus that spilled over from the bat reservoir. However, the host genetic determinants that drive SARS-CoV-2 susceptibility and COVID-19 severity are largely unknown. Understanding how cellular proteins interacting with SARS-CoV-2 have evolved in primates and bats is of primary importance to decipher differences in the infection outcome between humans and the viral reservoir in bats. Here, we performed comparative functional genetic analyses of hundreds of SARS-CoV-2-interacting proteins to study virus–host interface adaptation over millions of years, pointing to genes similarly—or differentially—engaged in evolutionary arms races and that may be at the basis of in vivo pathogenic differences.

The coronavirus disease 19 (COVID-19) pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a coronavirus that spilled over from the bat reservoir. Despite numerous clinical trials and vaccines, the burden remains immense, and the host determinants of SARS-CoV-2 susceptibility and COVID-19 severity remain largely unknown. Signatures of positive selection detected by comparative functional genetic analyses in primate and bat genomes can uncover important and specific adaptations that occurred at virus–host interfaces. We performed high-throughput evolutionary analyses of 334 SARS-CoV-2-interacting proteins to identify SARS-CoV adaptive loci and uncover functional differences between modern humans, primates, and bats. Using DGINN (Detection of Genetic INNovation), we identified 38 bat and 81 primate proteins with marks of positive selection. Seventeen genes, including the ACE2 receptor, present adaptive marks in both mammalian orders, suggesting common virus–host interfaces and past epidemics of coronaviruses shaping their genomes. Yet, 84 genes presented distinct adaptations in bats and primates. Notably, residues involved in ubiquitination and phosphorylation of the inflammatory RIPK1 have rapidly evolved in bats but not primates, suggesting different inflammation regulation versus humans. Furthermore, we discovered residues with typical virus–host arms race marks in primates, such as in the entry factor TMPRSS2 or the autophagy adaptor FYCO1, pointing to host-specific in vivo interfaces that may be drug targets. Finally, we found that FYCO1 sites under adaptation in primates are those associated with severe COVID-19, supporting their importance in pathogenesis and replication. Overall, we identified adaptations involved in SARS-CoV-2 infection in bats and primates, enlightening modern genetic determinants of virus susceptibility and severity.

OBI: A computational tool for the analysis and systematization of the positive selection in proteins

There are multiple tools for positive selection analysis, including vaccine design and detection of variants of circulating drug-resistant pathogens in population selection. However, applying these tools to analyze a large number of protein families or as part of a comprehensive phylogenomics pipeline could be challenging. Since many standard bioinformatics tools are only available as executables, integrating them into complex Bioinformatics pipelines may not be possible.

We have developed OBI, an open-source tool aimed to facilitate positive selection analysis on a large scale. It can be used as a stand-alone command-line app that can be easily installed and used as a Conda package.

Some advantages of using OBI are:•

It speeds up the analysis by automating the entire process

•

It allows multiple starting points and customization for the analysis

•

It allows the retrieval and linkage of structural and evolutive data for a protein through

We hope to provide with OBI a solution for reliably speeding up large-scale protein evolutionary and structural analysis.

Rapid Evolution of HERC6 and Duplication of a Chimeric HERC5/6 Gene in Rodents and Bats Suggest an Overlooked Role of HERCs in Mammalian Immunity

Studying the evolutionary diversification of mammalian antiviral defenses is of main importance to better understand our innate immune repertoire. The small HERC proteins are part of a multigene family, including HERC5 and HERC6, which have probably diversified through complex evolutionary history in mammals. Here, we performed mammalian-wide phylogenetic and genomic analyses of HERC5 and HERC6, using 83 orthologous sequences from bats, rodents, primates, artiodactyls, and carnivores—the top five representative groups of mammalian evolution. We found that HERC5 has been under weak and differential positive selection in mammals, with only primate HERC5 showing evidences of pathogen-driven selection. In contrast, HERC6 has been under strong and recurrent adaptive evolution in mammals, suggesting past and widespread genetic arms-races with viral pathogens. Importantly, the rapid evolution of mammalian HERC6 spacer domain suggests that it might be a host-pathogen interface, targeting viral proteins and/or being the target of virus antagonists. Finally, we identified a HERC5/6 chimeric gene that arose from independent duplication in rodent and bat lineages and encodes for a conserved HERC5 N-terminal domain and divergent HERC6 spacer and HECT domains. This duplicated chimeric gene highlights adaptations that potentially contribute to rodent and bat immunity. Our findings open new research avenues on the functions of HERC6 and HERC5/6 in mammals, and on their implication in antiviral innate immunity.