Intrinsically disordered regions are abundant in simplexvirus proteomes and display signatures of positive selection

Population Genetic Structure and Human Adaptation of Kaposi Sarcoma-Associated Herpesvirus

Background

Kaposi sarcoma-associated herpesvirus (KSHV), the etiologic agent of Kaposi sarcoma, is human-specific and is thought to have emerged from primate-infecting gammaherpesviruses. KSHV seroprevalence shows geographic variation, being highest in sub-Saharan Africa, intermediate in the Mediterranean area, and low in most other locations. However, KSHV prevalence is also particularly high in specific regions such as the Miyako Islands (Japan).

Methods

We retrieved KSHV genomes from public repositories and analyzed geographic patterns using principal component analysis and STRUCTURE. Adaptation to the human host was investigated by likelihood ratio tests for positive selection. Protein structures were derived from the HerpesFolds database.

Results

Most non-African genomes are genetically separated by the African genomes, and the latter are divided into 2 main lineages. The African genomes received most of their ancestry from 2 populations showing limited drift, suggesting an African origin for circulating KSHV strains. Several non-African genomes instead have most of their ancestry covered by a highly drifted ancestral population. However, some non-African genomes show similar ancestry proportions to the African ones, including those from Miyako Islands and the variant F subtype sampled in France. Molecular analysis of adaptation to the human host identified core genes as the major selection targets, including 2 viral enzymes that counteract human immune defenses.

Conclusions

We suggest that the genetic diversity of extant strains reflects relatively recent demographic events associated with viral lineage extinctions, which may have influenced KSHV epidemiology. Adaptation to the human host involved changes in core genes, possibly a strategy to optimize protein-protein interactions.

Evolutionary analysis of ZAP and its cofactors identifies intrinsically disordered regions as central elements in host-pathogen interactions

The zinc-finger antiviral protein (ZAP) is an innate immunity sensor of non-self nucleic acids. Its antiviral activity is exerted through the physical interaction with different cofactors, including TRIM25, Riplet and KHNYN. Cellular proteins that interact with infectious agents are expected to be engaged in genetic conflicts that often result in their rapid evolution. To test this possibility and to identify the regions most strongly targeted by natural selection, we applied in silico molecular evolution tools to analyze the evolutionary history of ZAP and cofactors in four mammalian groups. We report evidence of positive selection in all genes and in most mammalian groups. On average, the intrinsically disordered regions (IDRs) embedded in the four proteins evolve significantly faster than folded domains and most positively selected sites fall within IDRs. In ZAP, the PARP domain also shows abundant signals of selection, and independent evolution in different mammalian groups suggests modulation of its ADP-ribose binding ability. Detailed analyses of the biophysical properties of IDRs revealed that chain compaction and conformational entropy are conserved across mammals. The IDRs in ZAP and KHNYN are particularly compact, indicating that they may promote phase separation (PS). In line with this hypothesis, we predicted several PS-promoting regions in ZAP and KHNYN, as well as in TRIM25. Positively selected sites are abundant in these regions, suggesting that PS may be important for the antiviral functions of these proteins and the evolutionary arms race with viruses. Our data shed light into the evolution of ZAP and cofactors and indicate that IDRs represent central elements in host-pathogen interactions.

Positive Selection Drives the Evolution of the Structural Maintenance of Chromosomes (SMC) Complexes

Structural Maintenance of Chromosomes (SMC) complexes are an evolutionary conserved protein family. In most eukaryotes, three SMC complexes have been characterized, as follows: cohesin, condensin, and SMC5/6 complexes. These complexes are involved in a plethora of functions, and defects in SMC genes can lead to an increased risk of chromosomal abnormalities, infertility, and cancer. To investigate the evolution of SMC complex genes in mammals, we analyzed their selective patterns in an extended phylogeny. Signals of positive selection were identified for condensin NCAPG, for two SMC5/6 complex genes (SMC5 and NSMCE4A), and for all cohesin genes with almost exclusive meiotic expression (RAD21L1, REC8, SMC1B, and STAG3). For the latter, evolutionary rates correlate with expression during female meiosis, and most positively selected sites fall in intrinsically disordered regions (IDRs). Our results support growing evidence that IDRs are fast evolving, and that they most likely contribute to adaptation through modulation of phase separation. We suggest that the natural selection signals identified in SMC complexes may be the result of different selective pressures: a host-pathogen arms race in the condensin and SMC5/6 complexes, and an intragenomic conflict for meiotic cohesin genes that is similar to that described for centromeres and telomeres.

Coronavirus accessory protein ORF3 biology and its contribution to viral behavior and pathogenesis

Coronavirus porcine epidemic diarrhea virus (PEDV) is classified in the genus Alphacoronavirus, family Coronaviridae that encodes the only accessory protein, ORF3 protein. However, how ORF3 contributes to viral pathogenicity, adaptability, and replication is obscure. In this review, we summarize current knowledge and identify gaps in many aspects of ORF3 protein in PEDV, with emphasis on its unique biological features, including membrane topology, Golgi retention mechanism, potential intrinsic disordered property, functional motifs, protein glycosylation, and codon usage phenotypes related to genetic evolution and gene expression. In addition, we propose intriguing questions related to ORF3 protein that we hope to stimulate further studies and encourage collaboration among virologists worldwide to provide constructive knowledge about the unique characteristics and biological functions of ORF3 protein, by which their potential role in clarifying viral behavior and pathogenesis can be possible.

Evolution of the orthopoxvirus core genome

Orthopoxviruses comprise several relevant pathogens, including the causative agent of smallpox and monkeypox virus. Analysis of orthopoxvirus genome evolution mainly focused on gene gains/losses. We instead analyzed core genes, which are conserved in all orthopoxviruses. We show that, despite their strong constraint, some genes involved in viral morphogenesis and transcription/replication were targets of pervasive positive selection, which was relatively uncommon in immunomodulatory genes. However at least three of the positively selected genes, E3L, A24R, and H3L, might have evolved in response to immune selection. Episodic positive selection was particularly common on the internal branches of the orthopox phylogeny and on the monkeypox virus lineage. The latter showed evidence of episodic positive selection at the D14L gene, which encodes a modulator of complement activation (MOPICE). Notably, two genes (B1R and A33R) targeted by episodic selection on more than one branch are involved in forms of intra-genomic conflict. Finally, we found that, in orthopoxvirus proteomes, intrinsically disordered regions (IDRs) tend to be less constrained and are common targets of positive selection. Extension of our analysis to all poxviruses showed no evidence that the IDR fraction differs with host range. Conversely, we found a strong effect of base composition, which was however not sufficient to explain IDR fraction. We thus suggest that, in poxviruses, the IDR fraction is maintained by modulating GC content to accommodate disorder-promoting codons. Overall, our data provide novel insight in orthopoxvirus evolution and provide a list of genes and sites that are expected to modulate viral phenotypes.

Evolution of SLiM-mediated hijack functions in intrinsically disordered viral proteins

Viruses and their hosts are involved in an ‘arms race’ where they continually evolve mechanisms to overcome each other. It has long been proposed that intrinsic disorder provides a substrate for the evolution of viral hijack functions and that short linear motifs (SLiMs) are important players in this process. Here, we review evidence in support of this tenet from two model systems: the papillomavirus E7 protein and the adenovirus E1A protein. Phylogenetic reconstructions reveal that SLiMs appear and disappear multiple times across evolution, providing evidence of convergent evolution within individual viral phylogenies. Multiple functionally related SLiMs show strong coevolution signals that persist across long distances in the primary sequence and occur in unrelated viral proteins. Moreover, changes in SLiMs are associated with changes in phenotypic traits such as host range and tropism. Tracking viral evolutionary events reveals that host switch events are associated with the loss of several SLiMs, suggesting that SLiMs are under functional selection and that changes in SLiMs support viral adaptation. Fine-tuning of viral SLiM sequences can improve affinity, allowing them to outcompete host counterparts. However, viral SLiMs are not always competitive by themselves, and tethering of two suboptimal SLiMs by a disordered linker may instead enable viral hijack. Coevolution between the SLiMs and the linker indicates that the evolution of disordered regions may be more constrained than previously thought. In summary, experimental and computational studies support a role for SLiMs and intrinsic disorder in viral hijack functions and in viral adaptive evolution.

Rapidly evolving viral motifs mostly target biophysically constrained binding pockets of host proteins

Evolutionary changes in host-virus interactions can alter the course of infection, but the biophysical and regulatory constraints that shape interface evolution remain largely unexplored. Here, we focus on viral mimicry of host-like motifs that allow binding to host domains and modulation of cellular pathways. We observe that motifs from unrelated viruses preferentially target conserved, widely expressed, and highly connected host proteins, enriched with regulatory and essential functions. The interface residues within these host domains are more conserved and bind a larger number of cellular proteins than similar motif-binding domains that are not known to interact with viruses. In contrast, rapidly evolving viral-binding human proteins form few interactions with other cellular proteins and display high tissue specificity, and their interfaces have few inter-residue contacts. Our results distinguish between conserved and rapidly evolving host-virus interfaces and show how various factors limit host capacity to evolve, allowing for efficient viral subversion of host machineries.

Simplexviruses successfully adapt to their host by fine-tuning immune responses

Primate herpes simplex viruses are species-specific and relatively harmless to their natural hosts. However, cross-species transmission is often associated with severe disease, as exemplified by the virulence of macacine herpesvirus 1 (B virus) in humans. We performed a genome-wide scan for signals of adaptation of simplexviruses to their hominin hosts. Among core genes, we found evidence of episodic positive selection in three glycoproteins, with several selected sites located in antigenic determinants. Positively selected noncore genes were found to be involved in different immune-escape mechanisms. The herpes simplex virus (HSV)-1/HSV-2 encoded product (ICP47) of one of these genes is known to down-modulate major histocompatibility complex class I expression. This feature is not shared with B virus, which instead up-regulates Human Leukocyte Antigen (HLA)-G, an immunomodulatory molecule. By in vitro expression of different ICP47 mutants, we functionally characterized the selection signals. Results indicated that the selected sites do not represent the sole determinants of binding to the transporter associated with antigen processing (TAP). Conversely, the amino acid status at these sites was sufficient to determine HLA-G up-regulation. In fact, both HSV-1 and HSV-2 ICP47 induced HLA-G when mutated to recapitulate residues in B virus, whereas the mutated version of B virus ICP47 failed to determine HLA-G expression. These differences might contribute to the severity of B virus infection in humans. Importantly, they indicate that the evolution of ICP47 in HSV-1/HSV-2 led to the loss of an immunosuppressive effect. Thus, related simplexviruses finely tune the balance between immunosuppressive and immunostimulatory pathways to promote successful co-existence with their primate hosts.

Tracking the pipeline: immunoinformatics and the COVID-19 vaccine design

With the onset of the COVID-19 pandemic, the amount of data on genomic and proteomic sequences of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) stored in various databases has exponentially grown. A large volume of these data has led to the production of equally immense sets of immunological data, which require rigorous computational approaches to sort through and make sense of. Immunoinformatics has emerged in the recent decades as a field capable of offering this approach by bridging experimental and theoretical immunology with state-of-the-art computational tools. Here, we discuss how immunoinformatics can assist in the development of high-performance vaccines and drug discovery needed to curb the spread of SARS-CoV-2. Immunoinformatics can provide a set of computational tools to extract meaningful connections from the large sets of COVID-19 patient data, which can be implemented in the design of effective vaccines. With this in mind, we represent a pipeline to identify the role of immunoinformatics in COVID-19 treatment and vaccine development. In this process, a number of free databases of protein sequences, structures and mutations are introduced, along with docking web servers for assessing the interaction between antibodies and the SARS-CoV-2 spike protein segments as most commonly considered antigens in vaccine design.

A holistic perspective on herpes simplex virus (HSV) ecology and evolution

Herpes simplex viruses (HSV) cause chronic infection in humans that are characterized by periodic episodes of mucosal shedding and ulcerative disease. HSV causes millions of infections world-wide, with lifelong bouts of viral reactivation from latency in neuronal ganglia. Infected individuals experience different levels of disease severity and frequency of reactivation. There are two distantly related HSV species, with HSV-1 infections historically found most often in the oral niche and HSV-2 infections in the genital niche. Over the last two decades, HSV-1 has emerged as the leading cause of first-episode genital herpes in multiple countries. While HSV-1 has the highest level of genetic diversity among human alpha-herpesviruses, it is not yet known how quickly the HSV-1 viral population in a human host adapts over time, or if there are population bottlenecks associated with viral reactivation and/or transmission. It is also unknown how the ecological environments in which HSV infections occur influence their evolutionary trajectory, or that of co-occurring viruses and microbes. In this review, we explore how HSV accrues genetic diversity within each new infection, and yet maintains its ability to successfully infect most of the human population. A holistic examination of the ecological context of natural human infections can expand our awareness of how HSV adapts as it moves within and between human hosts, and reveal the complexity of these lifelong human-virus interactions. These insights may in turn suggest new areas of exploration for other chronic pathogens that successfully evolve and persist among their hosts.

Proline-Rich Hypervariable Region of Hepatitis E Virus: Arranging the Disorder

The hepatitis E virus (HEV) hypervariable region (HVR) presents the highest divergence of the entire HEV genome. It is characteristically rich in proline, and so is also known as the “polyproline region” (PPR). HEV genotype 3 (HEV-3) exhibits different PPR lengths due to insertions, PPR and/or RNA-dependent RNA polymerase (RdRp) duplications and deletions. A total of 723 PPR-HEV sequences were analyzed, of which 137 HEV-3 sequences were obtained from clinical specimens (from acute and chronic infection) by Sanger sequencing. Eight swine stool/liver samples were also analyzed. N- and C-terminal fragments were confirmed as being conserved, but they harbored differences between genotypes and were not proline-plentiful regions. The genuine PPR is the intermediate region between them. HEV-3 PPR contains a higher percentage (30.4%) of prolines than other genotypes. We describe for the first time: (1) the specific placement of HEV-3 PPR rearrangements in sites 1 to 14 of the PPR, noting that duplications are more frequently attached to sites 11 and 12 (AAs 74–79 and 113–118, respectively); (2) the cadence of repetitions follows a circular-like pattern of blocks A to J, with F, G, H, and I being the most frequent; (3) a previously unreported insertion homologous to apolipoprotein C1; and (4) the increase in frequency of potential N-glycosylation sites and differences in AAs composition related to duplications.