Conditional Immune Escape during Chronic Simian Immunodeficiency Virus Infection.

UNLABELLED: Anti-HIV CD8 T cells included in therapeutic treatments will need to target epitopes that do not accumulate escape mutations. Identifying the epitopes that do not accumulate variants but retain immunogenicity depends on both host major histocompatibility complex (MHC) genetics and the likelihood for an epitope to tolerate variation. We previously found that immune escape during acute SIV infection is conditional; the accumulation of mutations in T cell epitopes is limited, and the rate of accumulation depends on the number of epitopes being targeted. We have now tested the hypothesis that conditional immune escape extends into chronic SIV infection and that epitopes with a preserved wild-type sequence have the potential to elicit epitope-specific CD8 T cells. We deep sequenced simian immunodeficiency virus (SIV) from Mauritian cynomolgus macaques (MCMs) that were homozygous and heterozygous for the M3 MHC haplotype and had been infected with SIV for about 1 year. When interrogating variation within individual epitopes restricted by M3 MHC alleles, we found three categories of epitopes, which we called categories A, B, and C. Category B epitopes readily accumulated variants in M3-homozygous MCMs, but this was less common in M3-heterozygous MCMs. We then determined that chronic CD8 T cells specific for these epitopes were more likely preserved in the M3-heterozygous MCMs than M3-homozygous MCMs. We provide evidence that epitopes known to escape from chronic CD8 T cell responses in animals that are homozygous for a set of MHC alleles are preserved and retain immunogenicity in a host that is heterozygous for the same MHC alleles. IMPORTANCE: Anti-HIV CD8 T cells that are part of therapeutic treatments will need to target epitopes that do not accumulate escape mutations. Defining these epitope sequences is a necessary precursor to designing approaches that enhance the functionality of CD8 T cells with the potential to control virus replication during chronic infection or after reactivation of latent virus. Using MHC-homozygous and -heterozygous Mauritian cynomolgus macaques, we have now obtained evidence that epitopes known to escape from chronic CD8 T cell responses in animals that are MHC homozygous are preserved and retain immunogenicity in a host that is heterozygous for the same MHC alleles. Importantly, our findings support the conditional immune escape hypothesis, such that the potential to present a greater number of CD8 T cell epitopes within a single animal can delay immune escape in targeted epitopes. As a result, certain epitope sequences can retain immunogenicity into chronic infection.

Arteriviruses, Pegiviruses, and Lentiviruses Are Common among Wild African Monkeys.

UNLABELLED: Nonhuman primates (NHPs) are a historically important source of zoonotic viruses and are a gold-standard model for research on many human pathogens. However, with the exception of simian immunodeficiency virus (SIV) (family Retroviridae), the blood-borne viruses harbored by these animals in the wild remain incompletely characterized. Here, we report the discovery and characterization of two novel simian pegiviruses (family Flaviviridae) and two novel simian arteriviruses (family Arteriviridae) in wild African green monkeys from Zambia (malbroucks [Chlorocebus cynosuros]) and South Africa (vervet monkeys [Chlorocebus pygerythrus]). We examine several aspects of infection, including viral load, genetic diversity, evolution, and geographic distribution, as well as host factors such as age, sex, and plasma cytokines. In combination with previous efforts to characterize blood-borne RNA viruses in wild primates across sub-Saharan Africa, these discoveries demonstrate that in addition to SIV, simian pegiviruses and simian arteriviruses are widespread and prevalent among many African cercopithecoid (i.e., Old World) monkeys. IMPORTANCE: Primates are an important source of viruses that infect humans and serve as an important laboratory model of human virus infection. Here, we discover two new viruses in African green monkeys from Zambia and South Africa. In combination with previous virus discovery efforts, this finding suggests that these virus types are widespread among African monkeys. Our analysis suggests that one of these virus types, the simian arteriviruses, may have the potential to jump between different primate species and cause disease. In contrast, the other virus type, the pegiviruses, are thought to reduce the disease caused by human immunodeficiency virus (HIV) in humans. However, we did not observe a similar protective effect in SIV-infected African monkeys coinfected with pegiviruses, possibly because SIV causes little to no disease in these hosts.

Pregnancy without progesterone in horses defines a second endogenous biopotent progesterone receptor agonist, 5α-dihydroprogesterone.

One of the most widely accepted axioms of mammalian reproductive biology is that pregnancy requires the (sole) support of progesterone, acting in large measure through nuclear progesterone receptors (PRs) in uterine and cervical tissues, without which pregnancy cannot be established or maintained. However, mares lack detectable progesterone in the latter half of pregnancy. Instead of progesterone, several (mainly 5α-reduced) pregnanes are elevated and have long been speculated to provide progestational support in lieu of progesterone itself. To the authors' knowledge, evidence for the bioactivity of a second potent endogenously synthesized pregnane able to support pregnancy in the absence of progesterone has never before been reported. The 5α-reduced progesterone metabolite dihydroprogesterone (DHP) was shown in vivo to stimulate endometrial growth and progesterone-dependent gene expression in the horse at subphysiological concentrations and to maintain equine pregnancy in the absence of luteal progesterone in the third and fourth weeks postbreeding. Results of in vitro studies indicate that DHP is an equally potent and efficacious endogenous progestin in the horse but that the PR evolved with increased agonistic potency for DHP at the expense of potency toward progesterone based on comparisons with human PR responses. Sequence analysis and available literature indicate that the enzyme responsible for DHP synthesis, 5α-reductase type 1, also adapted primarily to metabolize progesterone and thereby to serve diverse roles in the physiology of pregnancy in mammals. Our confirmation that endogenously synthesized DHP is a biopotent progestin in the horse ends decades of speculation, explaining how equine pregnancies survive without measurable circulating progesterone in the last 4 to 5 mo of gestation.

Selection Maintains Low Genomic GC Content in Marine SAR11 Lineages.

The genomic G+C content of ocean bacteria varies from below 30% to over 60%. This broad range of base composition is likely shaped by distinct mutational processes, recombination, effective population size, and selection driven by environmental factors. A number of studies have hypothesized that depletion of G/C in genomes of marine bacterioplankton cells is an adaptation to the nitrogen-poor pelagic oceans, but they failed to disentangle environmental factors from mutational biases and population history. Here, we reconstructed the evolutionary changes of bases at synonymous sites in genomes of two marine SAR11 populations and a freshwater counterpart with its evolutionary origin rooted in the marine lineage. Although they all have similar genome sizes, DNA repair gene repertoire, and base compositions, there is a stronger bias toward A/T changes, a reduced frequency of nitrogenous amino acids, and an exclusive occurrence of polyamine, opine, and taurine transport systems in the ocean populations, consistent with a greater nitrogen stress in surface oceans compared with freshwater lakes. Furthermore, the ratio of nonsynoymous to synonymous nucleotide diversity is not statistically distinguishable among these populations, suggesting that population history has a limited effect. Taken together, the ecological transition of SAR11 from ocean to freshwater habitats makes nitrogen more available to these organisms, and thus relaxation of purifying selection drove a genome-wide reduction in the frequency of G/C to A/T changes in the freshwater population.

SNPGenie: estimating evolutionary parameters to detect natural selection using pooled next-generation sequencing data.

UNLABELLED: New applications of next-generation sequencing technologies use pools of DNA from multiple individuals to estimate population genetic parameters. However, no publicly available tools exist to analyse single-nucleotide polymorphism (SNP) calling results directly for evolutionary parameters important in detecting natural selection, including nucleotide diversity and gene diversity. We have developed SNPGenie to fill this gap. The user submits a FASTA reference sequence(s), a Gene Transfer Format (.GTF) file with CDS information and a SNP report(s) in an increasing selection of formats. The program estimates nucleotide diversity, distance from the reference and gene diversity. Sites are flagged for multiple overlapping reading frames, and are categorized by polymorphism type: nonsynonymous, synonymous, or ambiguous. The results allow single nucleotide, single codon, sliding window, whole gene and whole genome/population analyses that aid in the detection of positive and purifying natural selection in the source population. AVAILABILITY AND IMPLEMENTATION: SNPGenie version 1.2 is a Perl program with no additional dependencies. It is free, open-source, and available for download at https://github.com/hugheslab/snpgenie. CONTACT: nelsoncw@email.sc.edu or austin@biol.sc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Life-history evolution at the molecular level: adaptive amino acid composition of avian vitellogenins.

Avian genomes typically encode three distinct vitellogenin (VTG) egg yolk proteins (VTG1, VTG2 and VTG3), which arose by gene duplication prior to the most recent common ancestor of birds. Analysis of VTG sequences from 34 avian species in a phylogenetic framework supported the hypothesis that VTG amino acid composition has co-evolved with embryo incubation time. Embryo incubation time was positively correlated with the proportions of dietary essential amino acids (EAAs) in VTG1 and VTG2, and with the proportion of sulfur-containing amino acids in VTG3. These patterns were seen even when only semi-altricial and/or altricial species were considered, suggesting that the duration of embryo incubation is a major selective factor on the amino acid composition of VTGs, rather than developmental mode alone. The results are consistent with the hypothesis that the level of EAAs provided to the egg represents an adaptation to the loss of amino acids through breakdown over the course of incubation and imply that life-history phenotypes and VTG amino acid composition have co-evolved throughout the evolutionary history of birds.

Stability of social behavior in captive American Flamingos (Phoenicopterus ruber): A quantitative case study.

Over a three-year period a captive American flamingo (Phoenicopterus ruber) group showed recurring seasonal changes in the pattern of agonistic behavior. In spite of seasonal changes in the rates of agonistic behavior, dominance relations were generally stable across seasons. Males had significantly higher dominance status than females, and birds in long-term male-female pair bonds had significantly higher dominance than others. Unresolved agonistic encounters (URI), those with no clear winner or loser, were frequently observed; and their pattern of occurrence changed over the course of the year. URI were most frequent on the breeding "island" during the breeding season, and at the feeder in the post-breeding season. Thus, the frequency of URI reliably tracked seasonal changes regarding the resources most in contention; and, the extent of an individual's involvement in URI was indicative of privileged access to resources over which competition occurred.

Cross-clade simultaneous HIV drug resistance genotyping for reverse transcriptase, protease, and integrase inhibitor mutations by Illumina MiSeq.

BACKGROUND: Viral resistance to antiretroviral therapy threatens our best methods to control and prevent HIV infection. Current drug resistance genotyping methods are costly, optimized for subtype B virus, and primarily detect resistance mutations to protease and reverse transcriptase inhibitors. With the increasing use of integrase inhibitors in first-line therapies, monitoring for integrase inhibitor drug resistance mutations is a priority. We designed a universal primer pair to PCR amplify all major group M HIV-1 viruses for genotyping using Illumina MiSeq to simultaneously detect drug resistance mutations associated with protease, nucleoside reverse transcriptase, non-nucleoside reverse transcriptase, and integrase inhibitors. RESULTS: A universal primer pair targeting the HIV pol gene was used to successfully PCR amplify HIV isolates representing subtypes A, B, C, D, CRF01_AE and CRF02_AG. The universal primers were then tested on 62 samples from a US cohort of injection drug users failing treatment after release from prison. 94% of the samples were successfully genotyped for known drug resistance mutations in the protease, reverse transcriptase and integrase gene products. Control experiments demonstrate that mutations present at ≥ 2% frequency are reliably detected and above the threshold of error for this method. New drug resistance mutations not found in the baseline sample were identified in 54% of the patient samples after treatment failure. 86% of patients with major drug resistance mutations had 1 or more mutations associated with drug resistance to the treatment regimen at the time point of treatment failure. 59% of the emerging mutations were found at frequencies between 2% and 20% of the total sequences generated, below the estimated limit of detection of current FDA-approved genotyping techniques. Primary plasma samples with viral loads as low as 799 copies/ml were successfully genotyped using this method. CONCLUSIONS: Here we present an Illumina MiSeq-based HIV drug resistance genotyping assay. Our data suggests that this universal assay works across all major group M HIV-1 subtypes and identifies all drug resistance mutations in the pol gene known to confer resistance to protease, reverse transcriptase and integrase inhibitors. This high-throughput and sensitive assay could significantly improve access to drug resistance genotyping worldwide.

Evolutionary diversification of the vertebrate transferrin multi-gene family.

In a phylogenetic analysis of vertebrate transferrins (TFs), six major clades (subfamilies) were identified: (a) S, the mammalian serotransferrins; (b) ICA, the mammalian inhibitor of carbonic anhydrase (ICA) homologs; (c) L, the mammalian lactoferrins; (d) O, the ovotransferrins of birds and reptiles; (e) M, the melanotransferrins of bony fishes, amphibians, reptiles, birds, and mammals; and (f) M-like, a newly identified TF subfamily found in bony fishes, amphibians, reptiles, and birds. A phylogenetic tree based on the joint alignment of N-lobes and C-lobes supported the hypothesis that three separate events of internal duplication occurred in vertebrate TFs: (a) in the common ancestor of the M subfamily, (b) in the common ancestor of the M-like subfamily, and (c) in the common ancestor of other vertebrate TFs. The S, ICA, and L subfamilies were found only in placental mammals, and the phylogenetic analysis supported the hypothesis that these three subfamilies arose by gene duplication after the divergence of placental mammals from marsupials. The M-like subfamily was unusual in several respects, including the presence of a uniquely high proportion of clade-specific conserved residues, including distinctive but conserved residues in the sites homologous to those functioning in carbonate binding of human serotransferrin. The M-like family also showed an unusually high proportion of cationic residues in the positively charged region corresponding to human lactoferrampin, suggesting a distinctive role of this region in the M-like subfamily, perhaps in antimicrobial defense.

A novel hepacivirus with an unusually long and intrinsically disordered NS5A protein in a wild Old World primate.

GB virus B (GBV-B; family Flaviviridae, genus Hepacivirus) has been studied in New World primates as a model for human hepatitis C virus infection, but the distribution of GBV-B and its relatives in nature has remained obscure. Here, we report the discovery of a novel and highly divergent GBV-B-like virus in an Old World monkey, the black-and-white colobus (Colobus guereza), in Uganda. The new virus, guereza hepacivirus (GHV), clusters phylogenetically with GBV-B and recently described hepaciviruses infecting African bats and North American rodents, and it shows evidence of ancient recombination with these other hepaciviruses. Direct sequencing of reverse-transcribed RNA from blood plasma from three of nine colobus monkeys yielded near-complete GHV genomes, comprising two distinct viral variants. The viruses contain an exceptionally long nonstructural 5A (NS5A) gene, approximately half of which codes for a protein with no discernible homology to known proteins. Computational structure-based analyses indicate that the amino terminus of the GHV NS5A protein may serve a zinc-binding function, similar to the NS5A of other viruses within the family Flaviviridae. However, the 521-amino-acid carboxy terminus is intrinsically disordered, reflecting an unusual degree of structural plasticity and polyfunctionality. These findings shed new light on the natural history and evolution of the hepaciviruses and on the extent of structural variation within the Flaviviridae.

Evolutionary diversification of aminopeptidase N in Lepidoptera by conserved clade-specific amino acid residues.

Members of the aminopepidase N (APN) gene family of the insect order Lepidoptera (moths and butterflies) bind the naturally insecticidal Cry toxins produced by the bacterium Bacillus thuringiensis. Phylogenetic analysis of amino acid sequences of seven lepidopteran APN classes provided strong support for the hypothesis that lepidopteran APN2 class arose by gene duplication prior to the most recent common ancestor of Lepidoptera and Diptera. The Cry toxin-binding region (BR) of lepidopteran and dipteran APNs was subject to stronger purifying selection within APN classes than was the remainder of the molecule, reflecting conservation of catalytic site and adjoining residues within the BR. Of lepidopteran APN classes, APN2, APN6, and APN8 showed the strongest evidence of functional specialization, both in expression patterns and in the occurrence of conserved derived amino acid residues. The latter three APN classes also shared a convergently evolved conserved residue close to the catalytic site. APN8 showed a particularly strong tendency towards class-specific conserved residues, including one of the catalytic site residues in the BR and ten others in close vicinity to the catalytic site residues. The occurrence of class-specific sequences along with the conservation of enzymatic function is consistent with the hypothesis that the presence of Cry toxins in the environment has been a factor shaping the evolution of this multi-gene family.

Escape from CD8(+) T cell responses in Mamu-B*00801(+) macaques differentiates progressors from elite controllers.

A small number of HIV-infected individuals known as elite controllers experience low levels of chronic phase viral replication and delayed progression to AIDS. Specific HLA class I alleles are associated with elite control, implicating CD8(+) T lymphocytes in the establishment of these low levels of viral replication. Most HIV-infected individuals that express protective HLA class I alleles, however, do not control viral replication. Approximately 50% of Mamu-B*00801(+) Indian rhesus macaques control SIVmac239 replication in the chronic phase in a manner that resembles elite control in humans. We followed both the immune response and viral evolution in SIV-infected Mamu-B*00801(+) animals to better understand the role of CD8(+) T lymphocytes during the acute phase of viral infection, when viral control status is determined. The virus escaped from immunodominant Vif and Nef Mamu-B*00801-restricted CD8(+) T lymphocyte responses during the critical early weeks of acute infection only in progressor animals that did not control viral replication. Thus, early CD8(+) T lymphocyte escape is a hallmark of Mamu-B*00801(+) macaques who do not control viral replication. By contrast, virus in elite controller macaques showed little evidence of variation in epitopes recognized by immunodominant CD8(+) T lymphocytes, implying that these cells play a role in viral control.

Evolutionary diversification of insect innexins.

Phylogenetic analysis of insect innexins supported the hypothesis that six major clades of insect innexins arose by gene duplication prior to the origin of the endopterygote insects. Within one of the six clades (the Zpg Clade), two independent gene duplication events were inferred to have occurred in the lineage of Drosophila, after the most recent common ancestor of the dipteran families Culicidae and Drosophilidae. The relationships among this clades were poorly resolved, except for a sister relationship between ShakB and Ogre. Gene expression data from FlyAtlas supported the hypothesis that the latter gene duplication events gave rise to functional differentiation, with Zpg showing a high level of expression in ovary, and Inx5 and Inx6 showing a high level of expression in testis. Because unduplicated members of this clade in Bombyx mori and Anopheles gambiae showed high levels of expression in both ovary and tests, the expression patterns of the Drosophila members of this clade provide evidence of subdivision of an ancestral gene function after gene duplication.

Analysis of hepatitis C virus intrahost diversity across the coding region by ultradeep pyrosequencing.

Hepatitis C virus (HCV) is the leading cause of liver disease worldwide. In this study, we analyzed four treatment-naïve patients infected with subtype 1a and performed Roche/454 pyrosequencing across the coding region. We report the presence of low-level drug resistance mutations that would most likely have been missed using conventional sequencing methods. The approach described here is broadly applicable to studies of viral diversity and could help to improve the efficacy of direct-acting antiviral agents (DAA) in the treatment of HCV-infected patients.

Conditional CD8+ T cell escape during acute simian immunodeficiency virus infection.

CD8+ T cell responses rapidly select viral variants during acute human immunodeficiency virus (HIV)/simian immunodeficiency virus (SIV) infection. We used pyrosequencing to examine variation within three SIV-derived epitopes (Gag386₋394GW9, Nef103₋111RM9, and Rev59₋68SP10) targeted by immunodominant CD8+ T cell responses in acutely infected Mauritian cynomolgus macaques. In animals recognizing all three epitopes, variation within Rev59₋68SP10 was associated with delayed accumulation of variants in Gag386₋394GW9 but had no effect on variation within Nef103₋111RM9. This demonstrates that the entire T cell repertoire, rather than a single T cell population, influences the timing of immune escape, thereby providing the first example of conditional CD8+ T cell escape in HIV/SIV infection.

Behavioral correlates of post-breeding weight change in a captive flock of American flamingos (Phoenicopterus ruber ruber).

In a captive flock of American flamingos (Phoenicopterus ruber ruber), a pattern of weight loss after the breeding season was observed in the majority of individuals of both sexes, but weights eventually stabilized after several months. Social dominance (as determined by a 6-month behavioral study) was a highly significant predictor of early post-breeding weight loss, with dominant individuals undergoing less severe weight loss over this period; and behavioral observations supported the hypothesis that social dominance enhanced access to food in both males and females. Later in the study period, social dominance was no longer a significant predictor of weight changes, but rather individuals that had undergone previous large weight losses seemed to spend a greater amount of time feeding in an apparent effort to offset weight loss. Over the course of the study, there was a marked decline in both the proportion of agonistic encounters that were resolved and in the proportion of agonistic encounters that occurred at the feeder. These trends occurred even though there was no evidence of a decrease in overall use of the feeder. Thus, dominance-related differences with respect to food access and post-breeding weight loss appeared to be transitory phenomena, which corrected themselves as levels of aggression eventually declined. Intervention to increase the equality of food access in the post-breeding period may not be necessary in captive flamingo flocks, since low social rank in most cases appears not to have harmful long-term consequences.

Genome reduction by deletion of paralogs in the marine cyanobacterium Prochlorococcus.

Several isolates of the marine cyanobacterial genus Prochlorococcus have smaller genome sizes than those of the closely related genus Synechococcus. In order to test whether loss of protein-coding genes has contributed to genome size reduction in Prochlorococcus, we reconstructed events of gene family evolution over a strongly supported phylogeny of 12 Prochlorococcus genomes and 9 Synechococcus genomes. Significantly, more events both of loss of paralogs within gene families and of loss of entire gene families occurred in Prochlorococcus than in Synechococcus. The number of nonancestral gene families in genomes of both genera was positively correlated with the extent of genomic islands (GIs), consistent with the hypothesis that horizontal gene transfer (HGT) is associated with GIs. However, even when only isolates with comparable extents of GIs were compared, significantly more events of gene family loss and of paralog loss were seen in Prochlorococcus than in Synechococcus, implying that HGT is not the primary reason for the genome size difference between the two genera.

Evolution of the ßGRP/GNBP/ß-1,3-glucanase family of insects.

The ßGRP/GNBP/ß-1,3-glucanase protein family of insects includes several proteins involved in innate immune recognition, such as the ß-glucan recognition proteins of Lepidoptera and the Gram-negative bacteria-binding proteins of Drosophila. A phylogenetic analysis supported the existence of two distinct subfamilies, designated the pattern recognition receptor (PRR) and glucanase subfamilies, which originated by gene duplication prior to the origin of the Holometabola. In the C-terminal region (CTR) shared by both subfamilies, the PRR subfamily has evolved significantly more rapidly at the amino acid sequence level than has the glucanase subfamily, implying a relative lack of constraint on the amino acid sequence of this region in the PRR subfamily. PRR subfamily members also include an N-terminal region (NTR), involved in carbohydrate recognition, which is not shared by glucanase subfamily members. In comparisons between paralogous PRR subfamily members, there were no conserved amino acid residues in the NTR. However, when pairs of putatively orthologous PRR subfamily members were compared, the NTR was most often as conserved as the CTR or more so. This pattern suggests that the NTR may be important in functions specific to the different paralogs, while amino acid sequence changes in the NTR may have been important in functional differentiation among paralogs, specifically with regard to the types of carbohydrates that they recognize.

Transmission of clonal hepatitis C virus genomes reveals the dominant but transitory role of CD8⁺ T cells in early viral evolution.

The RNA genome of the hepatitis C virus (HCV) diversifies rapidly during the acute phase of infection, but the selective forces that drive this process remain poorly defined. Here we examined whether Darwinian selection pressure imposed by CD8(+) T cells is a dominant force driving early amino acid replacement in HCV viral populations. This question was addressed in two chimpanzees followed for 8 to 10 years after infection with a well-defined inoculum composed of a clonal genotype 1a (isolate H77C) HCV genome. Detailed characterization of CD8(+) T cell responses combined with sequencing of recovered virus at frequent intervals revealed that most acute-phase nonsynonymous mutations were clustered in class I epitopes and appeared much earlier than those in the remainder of the HCV genome. Moreover, the ratio of nonsynonymous to synonymous mutations, a measure of positive selection pressure, was increased 50-fold in class I epitopes compared with the rest of the HCV genome. Finally, some mutation of the clonal H77C genome toward a genotype 1a consensus sequence considered most fit for replication was observed during the acute phase of infection, but the majority of these amino acid substitutions occurred slowly over several years of chronic infection. Together these observations indicate that during acute hepatitis C, virus evolution was driven primarily by positive selection pressure exerted by CD8(+) T cells. This influence of immune pressure on viral evolution appears to subside as chronic infection is established and genetic drift becomes the dominant evolutionary force.