A Gene-Based Positive Selection Detection Approach to Identify Vaccine Candidates Using Toxoplasma gondii as a Test Case Protozoan Pathogen

Genetic Diversity of Merozoite Surface Antigens in Global Babesia bovis Populations

Cattle can be severely infected with the tick-borne protozoa Babesia bovis, giving rise to serious economic losses. Invasion of the host's RBCs by the parasite merozoite/sporozoites depends largely on the MSA (merozoite surface antigens) gene family, which comprises various fragments, e.g., MSA-1, MSA-2a1, MSA-2a2, MSA-2b and MSA-2c, highlighting the importance of these antigens as vaccine candidates. However, experimental trials documented the failure of some developed MSA-based vaccines to fully protect animals from B. bovis infection. One reason for this failure may be related to the genetic structure of the parasite. In the present study, all MSA-sequenced B. bovis isolates on the GenBank were collected and subjected to various analyses to evaluate their genetic diversity and population structure. The analyses were conducted on 199 MSA-1, 24 MSA-2a1, 193 MSA-2b and 148 MSA-2c isolates from geographically diverse regions. All these fragments displayed high nucleotide and haplotype diversities, but the MSA-1 was the most hypervariable and had the lowest inter- and intra-population gene flow values. This fragment also displayed a strong positive selection when testing its isolates for the natural selection, which suggests the potential occurrence of more genetic variations. On the contrary, the MSA-2c was the most conserved in comparison to the other fragments, and displayed the highest inter- and intra-population gene flow values, which was evidenced by a significantly negative selection and negative neutrality indices (Fu's Fs and Tajima's D). The majority of the MSA-2c tested isolates had two conserved amino acid repeats, and earlier reports have found these repeats to be highly immunogenic, which underlines the importance of this fragment in developing vaccines against B. bovis. Results of the MSA-2a1 analyses were also promising, but many more MSA-2a1 sequenced isolates are required to validating this assumption. The genetic analyses conducted for the MSA-2b fragment displayed borderline values when compared to the other fragments.

A state-of-the-art methodology for high-throughput in silico vaccine discovery against protozoan parasites and exemplified with discovered candidates for Toxoplasma gondii

Vaccine discovery against eukaryotic parasites is not trivial as highlighted by the limited number of known vaccines compared to the number of protozoal diseases that need one. Only three of 17 priority diseases have commercial vaccines. Live and attenuated vaccines have proved to be more effective than subunit vaccines but adversely pose more unacceptable risks. One promising approach for subunit vaccines is in silico vaccine discovery, which predicts protein vaccine candidates given thousands of target organism protein sequences. This approach, nonetheless, is an overarching concept with no standardised guidebook on implementation. No known subunit vaccines against protozoan parasites exist as a result of this approach, and consequently none to emulate. The study goal was to combine current in silico discovery knowledge specific to protozoan parasites and develop a workflow representing a state-of-the-art approach. This approach reflectively integrates a parasite's biology, a host's immune system defences, and importantly, bioinformatics programs needed to predict vaccine candidates. To demonstrate the workflow effectiveness, every Toxoplasma gondii protein was ranked in its capacity to provide long-term protective immunity. Although testing in animal models is required to validate these predictions, most of the top ranked candidates are supported by publications reinforcing our confidence in the approach.

Genomic coevolution of papillomavirus and immune system in placental mammals indicates the role of IFN-γ in the emergence of new variants

Papillomaviruses (PVs) are causative agents for warts and cancers in different parts of the body in the mammalian lineage. Therefore, these viruses are proposed as model organisms to study host immune responses to pathogens causing chronic infections. The virus-associated cancer progression depends on two integral processes namely angiogenesis and immune response (AIR). The angiogenesis process aids in tumour progression through vessel formation and maturation but the host immune response, in contrast, makes every attempt to eliminate pathogens and thereby maintain healthy tissues. However, the evolutionary contribution of individual viral genes and host AIR genes in carcinogenesis is yet to be explored. Here, we applied the evolutionary genomics approach to find correlated evolution between six PV genes and 23 host AIR-related genes. We estimated that IFN-γ is the only host gene evolving in a correlated manner with all six PV genes under study. Furthermore, three papillomavirus genes, L2, E6, and E7, are found to interact with two third of host AIR-related genes. Moreover, a combined differential gene expression analysis and network analysis showed that inflammatory cytokine IFN-γ is the key regulator of hub genes in the PPI network of the differentially expressed genes. Functional enrichment of these hub genes is consistent with their established role in different cancers and viral infections. Overall, we conclude that IFN-γ maintains selective pressure on mammalian PV genes and seems to be a potential biomarker for PV-related cancers. This study demonstrates the evolutionary importance of IFN-γ in deciding the fate of carcinogenic PV variants.

Evaluating the Potential Fitness Effects of Chinook Salmon (Oncorhynchus tshawytscha) Aquaculture Using Non-Invasive Population Genomic Analyses of MHC Nucleotide Substitution Spectra

Genetic diversity plays a vital role in the adaptability of salmon to changing environmental conditions that can introduce new selective pressures on populations. Variability among local subpopulations may increase the chance that certain advantageous genes are passed down to future generations to mitigate susceptibility to novel diseases, warming oceans, loss of genetic stocks, and ocean acidification. Class I and II genes of the major histocompatibility complex (MHC) are crucial for the fitness of Chinook salmon due to the role they play in disease and pathogen resistance. The objective of this study was to assess the DNA sequence variability among wild and hatchery populations of Alaskan Chinook salmon at the class I α1 and class II β1 exons of the MHC. We hypothesized that the 96 wild samples taken from the Deshka River would display greater levels of observed heterozygosity (Ho) relative to expected heterozygosity (He) in suggesting that individuals with similar phenotypes mate with one another more frequently than would be expected under random mating patterns. Conversely, since no mate selection occurs in the William Jack Hernandez Sport Fish hatchery, we would not expect to see this discrepancy (He = Ho) in the 96 hatchery fish tested in this study. Alternatively, we hypothesized that post-mating selection is driving higher levels of observed heterozygosity as opposed to mate selection. If this is the case, we will observe higher than expected levels of heterozygosity among hatchery salmon. Both populations displayed higher levels of observed heterozygosity than expected heterozygosity at the Class I and II loci but genetic differentiation between the spatially distinct communities was minimal. Class I sequences showed evidence of balancing selection, despite high rates of non-synonymous substitutions observed, specifically at the peptide binding regions of both MHC genes.

Population genetic analyses unveiled genetic stratification and differential natural selection signatures across the G-gene of viral hemorrhagic septicemia virus

Introduction: Viral hemorrhagic septicemia virus (VHSV) is the most lethal pathogen in aquaculture, infecting more than 140 fish species in marine, estuarine, and freshwater environments. Viral hemorrhagic septicemia virus is an enveloped RNA virus that belongs to the family Rhabdoviridae and the genus Novirhabdovirus. The current study is designed to infer the worldwide Viral hemorrhagic septicemia virus isolates' genetic diversity and evolutionary dynamics based on G-gene sequences. Methods: The complete G-gene sequences of viral hemorrhagic septicemia virus were retrieved from the public repositories with known timing and geography details. Pairwise statistical analysis was performed using Arlequin. The Bayesian model-based approach implemented in STRUCTURE software was used to investigate the population genetic structure, and the phylogenetic tree was constructed using MEGA X and IQ-TREE. The natural selection analysis was assessed using different statistical approaches, including IFEL, MEME, and SLAC. Results and Discussion: The global Viral hemorrhagic septicemia virus samples are stratified into five genetically distinct subpopulations. The STRUCTURE analysis unveiled spatial clustering of genotype Ia into two distinct clusters at K = 3. However, at K = 5, the genotype Ia samples, deposited from Denmark, showed temporal distribution into two groups. The analyses unveiled that the genotype Ia samples stratified into subpopulations possibly based on spatiotemporal distribution. Several viral hemorrhagic septicemia virus samples are characterized as genetically admixed or recombinant. In addition, differential or subpopulation cluster-specific natural selection signatures were identified across the G-gene codon sites among the viral hemorrhagic septicemia virus isolates. Evidence of low recombination events elucidates that genetic mutations and positive selection events have possibly driven the observed genetic stratification of viral hemorrhagic septicemia virus samples.

Relationship between genome-wide and MHC class I and II genetic diversity and complementarity in a nonhuman primate

Although mate choice is expected to favor partners with advantageous genetic properties, the relative importance of genome‐wide characteristics, such as overall heterozygosity or kinship, versus specific loci, is unknown. To disentangle genome‐wide and locus‐specific targets of mate choice, we must first understand congruence in global and local variation within the same individual. This study compares genetic diversity, both absolute and relative to other individuals (i.e., complementarity), assessed across the genome to that found at the major histocompatibility complex (MHC), a hyper‐variable gene family integral to immune system function and implicated in mate choice across species. Using DNA from 22 captive olive baboons (Papio anubis), we conducted double digest restriction site‐associated DNA sequencing to estimate genome‐wide heterozygosity and kinship, and sequenced two class I and two class II MHC loci. We found that genome‐wide diversity was not associated with MHC diversity, and that diversity at class I MHC loci was not correlated with diversity at class II loci. Additionally, kinship was a significant predictor of the number of MHC alleles shared between dyads at class II loci. Our results provide further evidence of the strong selective pressures maintaining genetic diversity at the MHC in comparison to other randomly selected sites throughout the genome. Furthermore, our results indicate that class II MHC disassortative mate choice may mediate inbreeding avoidance in this population. Our study suggests that mate choice favoring genome‐wide genetic diversity is not always synonymous with mate choice favoring MHC diversity, and highlights the importance of controlling for kinship when investigating MHC‐associated mate choice.

Evidence of episodic positive selection in Corynebacterium diphtheriae complex of species and its implementations in identification of drug and vaccine targets

BACKGROUND

Within the pathogenic bacterial species Corynebacterium genus, six species that can produce diphtheria toxin (C. belfantii, C. diphtheriae, C. pseudotuberculosis, C. rouxii, C. silvaticum and C. ulcerans) form a clade referred to as the C. diphtheria complex. These species have been found in humans and other animals, causing diphtheria or other diseases. Here we show the results of a genome scale analysis to identify positive selection in protein-coding genes that may have resulted in the adaptations of these species to their ecological niches and suggest drug and vaccine targets.

METHODS

Forty genomes were sampled to represent species, subspecies or biovars of Corynebacterium. Ten phylogenetic groups were tested for positive selection using the PosiGene pipeline, including species and biovars from the C. diphtheria complex. The detected genes were tested for recombination and had their sequences alignments and homology manually examined. The final genes were investigated for their function and a probable role as vaccine or drug targets.

RESULTS

Nineteen genes were detected in the species C. diphtheriae (two), C. pseudotuberculosis (10), C. rouxii (one), and C. ulcerans (six). Those were found to be involved in defense, translation, energy production, and transport and in the metabolism of carbohydrates, amino acids, nucleotides, and coenzymes. Fourteen were identified as essential genes, and six as virulence factors. Thirteen from the 19 genes were identified as potential drug targets and four as potential vaccine candidates. These genes could be important in the prevention and treatment of the diseases caused by these bacteria.

An Overview of Current Uses and Future Opportunities for Computer-Assisted Design of Vaccines for Neglected Tropical Diseases

Neglected tropical diseases are infectious diseases that impose high morbidity and mortality rates over 1.5 billion people worldwide. Originally restricted to tropical and subtropical regions, changing climate conditions have increased their potential to emerge elsewhere. Control of their impact suffers from shortages like poor epidemiological surveillance or irregular drug distribution, and some NTDs still lack of appropriate diagnostics and/or efficient therapeutics. For these, availability of vaccines to prevent new infections, or the worsening of those already established, would mean a major breakthrough. However, only dengue and rabies count with approved vaccines at present. Herein, we review the state-of-the-art of vaccination strategies for NTDs, setting the focus on third generation vaccines and the concept of reverse vaccinology. Its capability to address pathogens´ biological complexity, likely contributing to save developmental costs is discussed. The use of computational tools is a fundamental aid to analyze increasingly large datasets aimed at designing vaccine candidates with the highest, possibly, opportunities to succeed. Ultimately, we identify and analyze those studies that took an in silico approach to find vaccine candidates, and experimentally assessed their immunogenicity and/or protection capabilities.

GWideCodeML: A Python Package for Testing Evolutionary Hypotheses at the Genome-Wide Level

One of the most widely used programs for detecting positive selection, at the molecular level, is the program codeml, which is implemented in the Phylogenetic Analysis by Maximum Likelihood (PAML) package. However, it has a limitation when it comes to genome-wide studies, as it runs on a gene-by-gene basis. Furthermore, the size of such studies will depend on the number of orthologous genes the genomes have income and these are often restricted to only account for instances where a one-to-one relationship is observed between the genomes. In this work, we present GWideCodeML, a Python package, which runs a genome-wide codeml with the option of parallelization. To maximize the number of analyzed genes, the package allows for a variable number of taxa in the alignments and will automatically prune the topology to fit each of them, before running codeml.

Detecting sequence variants in clinically important protozoan parasites

Second and third generation sequencing methods are crucial for population genetic studies, and variant detection is a popular approach for exploiting this sequence data. While mini- and microsatellites are historically useful markers for studying important Protozoa such as Toxoplasma and Plasmodium spp., detecting non-repetitive variants such as those found in genes can be fundamental to investigating a pathogen's biology. These variants, namely single nucleotide polymorphisms and insertions and deletions, can help elucidate the genetic basis of an organism's pathogenicity, identify selective pressures, and resolve phylogenetic relationships. They also have the added benefit of possessing a comparatively low mutation rate, which contributes to their stability. However, there is a plethora of variant analysis tools with nuanced pipelines and conflicting recommendations for best practise, which can be confounding. This lack of standardisation means that variant analysis requires careful parameter optimisation, an understanding of its limitations, and the availability of high quality data. This review explores the value of variant detection when applied to non-model organisms such as clinically important protozoan pathogens. The limitations of current methods are discussed, including special considerations that require the end-users' attention to ensure that the results generated are reproducible, and the biological conclusions drawn are valid.

Neospora caninum antigens displaying virus-like particles as a bivalent vaccine candidate against neosporosis

Neospora caninum is a causative and transmissible agent of dog and bovine neosporosis. The resulting reproductive failures in infected cattle lead to significant economic losses worldwide. However, there is no satisfactory treatment or vaccine currently available to combat this pathogen. Thus, the development of appropriate vaccines to manage its infection and transmission is urgently needed. In this study, we expressed Rous sarcoma virus-like particles (RSV-LP) that displayed dual N. caninum antigens in silkworms. The antigen candidates are modified by adding a transmembrane domain of GP64 protein from Bombyx mori nucleopolyhedrovirus (BmNPV) to the C-terminus of surface antigen 1 (NcSAG1) and SAG1-related sequence 2 (NcSRS2). The NcSRS2 alone or the NcSAG1/NcSRS2 bivalent form displaying RSV-LPs were purified using sucrose density gradient centrifugation. These purified VLPs were then used for immunizations in gerbils, Meriones unguiculatus, to evaluate the anti-N. caninum effects in vivo. The results demonstrated that antigens displaying RSV-LPs in immunized gerbils produced the antigen-specific antibody, leading to a relatively lower parasite load after infections of N. caninum. To the best of our knowledge, this is the first study to present an RSV-LP vaccine displaying bivalent antigens from neosporosis. Taken together, our strategy suggests that silkworm-expressed virus-like particles (VLPs) are promising bivalent vaccine candidates against N. caninum infections.

Comparative Genomic Analysis of Trichinella spiralis Reveals Potential Mechanisms of Adaptive Evolution

Trichinellosis caused by parasitic nematodes of the genus Trichinella may result in human morbidity and mortality worldwide. Deciphering processes that drive species diversity and adaptation are key to understanding parasitism and developing effective control strategies. Our goal was to identify genes that are under positive selection and possible mechanisms of adaptive evolution of Trichinella spiralis genes using a comparative genomic analysis with the genomes of Brugia malayi, Trichuris suis, Ancylostoma ceylanicum, and Caenorhabditis elegans. The CODEML program derived from the PAML package was used to deduce the most probable dN/dS ratio, a measurement to detect genes/proteins undergoing adaptation. For each pair of sequences, those with a dN/dS ratio > 1 were considered positively selected genes (PSGs). Altogether, 986 genes were positively selected (p-value < 0.01). Genes involved in metabolic pathways, signaling pathways, and cytosolic DNA-sensing pathways were significantly enriched among the PSGs. Several PSGs are associated with exploitation of the host: modification of the host's metabolism, creation of new parasite-specific morphological structures between T. spiralis and the host interface, xenobiotic metabolism to combat low oxygen concentrations and host toxicity, muscle cell transformation, cell cycle arrest, DNA repair processes during nurse cell formation, antiapoptotic factors, immunomodulation, and regulation of epigenetic processes. Some of the T. spiralis PSGs have C. elegans orthologs that confer severe or lethal RNAi phenotypes. Fifty-seven PSGs in T. spiralis were analyzed to encode differentially expressed proteins. The present study utilized an overall comparative genomic analysis to discover PSGs within T. spiralis and their relationships with biological function and organism fitness. This analysis adds to our understanding of the possible mechanism that contributes to T. spiralis parasitism and biological adaptation within the host, and thus these identified genes may be potential targets for drug and vaccine development.

Genome-wide analyses reveal genes subject to positive selection in Toxoplasma gondii

Toxoplasma gondii is an important protozoan pathogen that infects many wild and domestic animals and causes infections in immunocompromised humans. However, there has been little investigation of the molecular evolutionary trajectories of this pathogenic protozoa using comparative genomics data. Here, we employed a comparative evolutionary genomics approach to identify genes that are under site- and lineage-specific positive selection in nine strains of T. gondii, including two closely related species, Neospora caninum and Hammondia hammondi. Based on the analyses of five coccidian core genomes, 4.5% of the 5788 core genome genes showed strong signals for positive selection in the site model. In addition, the branch-site model analyses in the nine T. gondii core genomes indicated that 2 to 20 genes underwent significant positive selection along each lineage leading to T. gondii strains. Many of the protein products encoded by the positively selected genes are secretory or surface proteins that have previously been implicated in host pathogenesis. The adaptive changes in these positively selected genes might be related to dynamic interactions between the host immune systems and might play a crucial role in the infection and pathogenic processes of T. gondii.

Advances in the Development of Anti-Toxoplasma gondii Vaccines: Challenges, Opportunities, and Perspectives

Important progress has been made in understanding how immunity is elicited against Toxoplasma gondii - a complex pathogen with multiple mechanisms of immune evasion. Many vaccine candidates have been tested using various strategies in animal models. However, none of these strategies has delivered as yet, and important challenges remain in the development of vaccines that can eliminate the tissue cysts and/or fully block vertical transmission. In this review, we provide an overview of the current understanding of the host immune response to T. gondii infection and summarize the key limitations for the development of an effective, safe, and durable toxoplasmosis vaccine. We also discuss how the successes and failures in developing and testing vaccine candidates have provided a roadmap for future vaccine development.