Pathogenesis and outcome of VA1 astrovirus infection in the human brain are defined by disruption of neural functions and imbalanced host immune responses.

Astroviruses (AstVs) can cause of severe infection of the central nervous system (CNS) in immunocompromised individuals. Here, we identified a human AstV of the VA1 genotype, HAstV-NIH, as the cause of fatal encephalitis in an immunocompromised adult. We investigated the cells targeted by AstV, neurophysiological changes, and host responses by analyzing gene expression, protein expression, and cellular morphology in brain tissue from three cases of AstV neurologic disease (AstV-ND). We demonstrate that neurons are the principal cells targeted by AstV in the brain and that the cerebellum and brainstem have the highest burden of infection. Detection of VA1 AstV in interconnected brain structures such as thalamus, deep cerebellar nuclei, Purkinje cells, and pontine nuclei indicates that AstV may spread between connected neurons transsynaptically. We found transcriptional dysregulation of neural functions and disruption of both excitatory and inhibitory synaptic innervation of infected neurons. Importantly, transcriptional dysregulation of neural functions occurred in fatal cases, but not in a patient that survived AstV-ND. We show that the innate, but not adaptive immune response was transcriptionally driving host defense in the brain of immunocompromised patients with AstV-ND. Both transcriptome and molecular pathology studies showed that most of the cellular changes were associated with CNS-intrinsic cells involved in phagocytosis and injury repair (microglia, perivascular/parenchymal border macrophages, and astrocytes), but not CNS-extrinsic cells (T and B cells), suggesting an imbalance of innate and adaptive immune responses to AstV infection in the brain as a result of the underlying immunodeficiencies. These results show that VA1 AstV infection of the brain in immunocompromised humans is associated with imbalanced host defense responses, disruption of neuronal somatodendritic compartments and synapses and increased phagocytic cellular activity. Improved understanding of the response to viral infections of the human CNS may provide clues for how to manipulate these processes to improve outcomes.

Role of humoral immunity against hepatitis B virus core antigen in the pathogenesis of acute liver failure.

Hepatitis B virus (HBV)-associated acute liver failure (ALF) is a dramatic clinical syndrome leading to death or liver transplantation in 80% of cases. Due to the extremely rapid clinical course, the difficulties in obtaining liver specimens, and the lack of an animal model, the pathogenesis of ALF remains largely unknown. Here, we performed a comprehensive genetic and functional characterization of the virus and the host in liver tissue from HBV-associated ALF and compared the results with those of classic acute hepatitis B in chimpanzees. In contrast with acute hepatitis B, HBV strains detected in ALF livers displayed highly mutated HBV core antigen (HBcAg), associated with increased HBcAg expression ex vivo, which was independent of viral replication levels. Combined gene and miRNA expression profiling revealed a dominant B cell disease signature, with extensive intrahepatic production of IgM and IgG in germline configuration exclusively targeting HBcAg with subnanomolar affinities, and complement deposition. Thus, HBV ALF appears to be an anomalous T cell-independent, HBV core-driven B cell disease, which results from the rare and unfortunate encounter between a host with an unusual B cell response and an infecting virus with a highly mutated core antigen.

A retrospective genomic analysis of drug-resistant strains of M. tuberculosis in a high-burden setting, with an emphasis on comparative diagnostics and reactivation and reinfection status.

BACKGROUND: Recurrence of drug-resistant tuberculosis (DR-TB) after treatment occurs through relapse of the initial infection or reinfection by a new drug-resistant strain. Outbreaks of DR-TB in high burden regions present unique challenges in determining recurrence status for effective disease management and treatment. In the Republic of Moldova the burden of DR-TB is exceptionally high, with many cases presenting as recurrent. METHODS: We performed a retrospective analysis of Mycobacterium tuberculosis from Moldova to better understand the genomic basis of drug resistance and its effect on the determination of recurrence status in a high DR-burden environment. To do this we analyzed genomes from 278 isolates collected from 189 patients, including 87 patients with longitudinal samples. These pathogen genomes were sequenced using Illumina technology, and SNP panels were generated for each sample for use in phylogenetic and network analysis. Discordance between genomic resistance profiles and clinical drug-resistance test results was examined in detail to assess the possibility of mixed infection. RESULTS: There were clusters of multiple patients with 10 or fewer differences among DR-TB samples, which is evidence of person-to-person transmission of DR-TB. Analysis of longitudinally collected isolates revealed that many infections exhibited little change over time, though 35 patients demonstrated reinfection by divergent (number of differences > 10) lineages. Additionally, several same-lineage sample pairs were found to be more divergent than expected for a relapsed infection. Network analysis of the H3/4.2.1 clade found very close relationships among 61 of these samples, making differentiation of reactivation and reinfection difficult. There was discordance between genomic profile and clinical drug sensitivity test results in twelve samples, and four of these had low level (but not statistically significant) variation at DR SNPs suggesting low-level mixed infections. CONCLUSIONS: Whole-genome sequencing provided a detailed view of the genealogical structure of the DR-TB epidemic in Moldova, showing that reinfection may be more prevalent than currently recognized. We also found increased evidence of mixed infection, which could be more robustly characterized with deeper levels of genomic sequencing.

The phylogeny of 48 alleles, experimentally verified at 21 kb, and its application to clinical allele detection.

BACKGROUND: Sequence information generated from next generation sequencing is often computationally phased using haplotype-phasing algorithms. Utilizing experimentally derived allele or haplotype information improves this prediction, as routinely used in HLA typing. We recently established a large dataset of long ERMAP alleles, which code for protein variants in the Scianna blood group system. We propose the phylogeny of this set of 48 alleles and identify evolutionary steps to derive the observed alleles. METHODS: The nucleotide sequence of > 21 kb each was used for all physically confirmed 48 ERMAP alleles that we previously published. Full-length sequences were aligned and variant sites were extracted manually. The Bayesian coalescent algorithm implemented in BEAST v1.8.3 was used to estimate a coalescent phylogeny for these variants and the allelic ancestral states at the internal nodes of the phylogeny. RESULTS: The phylogenetic analysis allowed us to identify the evolutionary relationships among the 48 ERMAP alleles, predict 4243 potential ancestral alleles and calculate a posterior probability for each of these unobserved alleles. Some of them coincide with observed alleles that are extant in the population. CONCLUSIONS: Our proposed strategy places known alleles in a phylogenetic framework, allowing us to describe as-yet-undiscovered alleles. In this new approach, which relies heavily on the accuracy of the alleles used for the phylogenetic analysis, an expanded set of predicted alleles can be used to infer alleles when large genotype data are analyzed, as typically generated by high-throughput sequencing. The alleles identified by studies like ours may be utilized in designing of microarray technologies, imputing of genotypes and mapping of next generation sequencing data.

Diminished viral replication and compartmentalization of hepatitis C virus in hepatocellular carcinoma tissue.

Analysis of hepatitis C virus (HCV) replication and quasispecies distribution within the tumor of patients with HCV-associated hepatocellular carcinoma (HCC) can provide insight into the role of HCV in hepatocarcinogenesis and, conversely, the effect of HCC on the HCV lifecycle. In a comprehensive study of serum and multiple liver specimens from patients with HCC who underwent liver transplantation, we found a sharp and significant decrease in HCV RNA in the tumor compared with surrounding nontumorous tissues, but found no differences in multiple areas of control non-HCC cirrhotic livers. Diminished HCV replication was not associated with changes in miR-122 expression. HCV genetic diversity was significantly higher in livers containing HCC compared with control non-HCC cirrhotic livers. Tracking of individual variants demonstrated changes in the viral population between tumorous and nontumorous areas, the extent of which correlated with the decline in HCV RNA, suggesting HCV compartmentalization within the tumor. In contrast, compartmentalization was not observed between nontumorous areas and serum, or in controls between different areas of the cirrhotic liver or between liver and serum. Our findings indicate that HCV replication within the tumor is restricted and compartmentalized, suggesting segregation of specific viral variants in malignant hepatocytes.

Whole-Genome Sequencing of Mycobacterium tuberculosis Provides Insight into the Evolution and Genetic Composition of Drug-Resistant Tuberculosis in Belarus.

The emergence and spread of drug-resistant Mycobacterium tuberculosis (DR-TB) are critical global health issues. Eastern Europe has some of the highest incidences of DR-TB, particularly multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB. To better understand the genetic composition and evolution of MDR- and XDR-TB in the region, we sequenced and analyzed the genomes of 138 M. tuberculosis isolates from 97 patients sampled between 2010 and 2013 in Minsk, Belarus. MDR and XDR-TB isolates were significantly more likely to belong to the Beijing lineage than to the Euro-American lineage, and known resistance-conferring loci accounted for the majority of phenotypic resistance to first- and second-line drugs in MDR and XDR-TB. Using a phylogenomic approach, we estimated that the majority of MDR-TB was due to the recent transmission of already-resistant M. tuberculosis strains rather than repeated de novo evolution of resistance within patients, while XDR-TB was acquired through both routes. Longitudinal sampling of M. tuberculosis from 34 patients with treatment failure showed that most strains persisted genetically unchanged during treatment or acquired resistance to fluoroquinolones. HIV+ patients were significantly more likely to have multiple infections over time than HIV- patients, highlighting a specific need for careful infection control in these patients. These data provide a better understanding of the genomic composition, transmission, and evolution of MDR- and XDR-TB in Belarus and will enable improved diagnostics, treatment protocols, and prognostic decision-making.

The TB Portals: an Open-Access, Web-Based Platform for Global Drug-Resistant-Tuberculosis Data Sharing and Analysis.

The TB Portals program is an international consortium of physicians, radiologists, and microbiologists from countries with a heavy burden of drug-resistant tuberculosis working with data scientists and information technology professionals. Together, we have built the TB Portals, a repository of socioeconomic/geographic, clinical, laboratory, radiological, and genomic data from patient cases of drug-resistant tuberculosis backed by shareable, physical samples. Currently, there are 1,299 total cases from five country sites (Azerbaijan, Belarus, Moldova, Georgia, and Romania), 976 (75.1%) of which are multidrug or extensively drug resistant and 38.2%, 51.9%, and 36.3% of which contain X-ray, computed tomography (CT) scan, and genomic data, respectively. The top Mycobacterium tuberculosis lineages represented among collected samples are Beijing, T1, and H3, and single nucleotide polymorphisms (SNPs) that confer resistance to isoniazid, rifampin, ofloxacin, and moxifloxacin occur the most frequently. These data and samples have promoted drug discovery efforts and research into genomics and quantitative image analysis to improve diagnostics while also serving as a valuable resource for researchers and clinical providers. The TB Portals database and associated projects are continually growing, and we invite new partners and collaborations to our initiative. The TB Portals data and their associated analytical and statistical tools are freely available at https://tbportals.niaid.nih.gov/.

Nanopore Sequencing as a Rapidly Deployable Ebola Outbreak Tool.

Rapid sequencing of RNA/DNA from pathogen samples obtained during disease outbreaks provides critical scientific and public health information. However, challenges exist for exporting samples to laboratories or establishing conventional sequencers in remote outbreak regions. We successfully used a novel, pocket-sized nanopore sequencer at a field diagnostic laboratory in Liberia during the current Ebola virus outbreak.

Spatiotemporal analysis of Guaroa virus diversity, evolution, and spread in South America.

We conducted phylogeographic modeling to determine the introduction and spread of Guaroa virus in South America. The results suggest a recent introduction of this virus into regions of Peru and Bolivia over the past 60-70 years and emphasize the need for increased surveillance in surrounding areas.

Profibrogenic chemokines and viral evolution predict rapid progression of hepatitis C to cirrhosis.

Chronic hepatitis C may follow a mild and stable disease course or progress rapidly to cirrhosis and liver-related death. The mechanisms underlying the different rates of disease progression are unknown. Using serial, prospectively collected samples from cases of transfusion-associated hepatitis C, we identified outcome-specific features that predict long-term disease severity. Slowly progressing disease correlated with an early alanine aminotransferase peak and antibody seroconversion, transient control of viremia, and significant induction of IFN-γ and MIP-1ß, all indicative of an effective, albeit insufficient, adaptive immune response. By contrast, rapidly progressive disease correlated with persistent and significant elevations of alanine aminotransferase and the profibrogenic chemokine MCP-1 (CCL-2), greater viral diversity and divergence, and a higher rate of synonymous substitution. This study suggests that the long-term course of chronic hepatitis C is determined early in infection and that disease severity is predicted by the evolutionary dynamics of hepatitis C virus and the level of MCP-1, a chemokine that appears critical to the induction of progressive fibrogenesis and, ultimately, the ominous complications of cirrhosis.

An epigenetic antimalarial resistance mechanism involving parasite genes linked to nutrient uptake.

Acquired antimalarial drug resistance produces treatment failures and has led to periods of global disease resurgence. In Plasmodium falciparum, resistance is known to arise through genome-level changes such as mutations and gene duplications. We now report an epigenetic resistance mechanism involving genes responsible for the plasmodial surface anion channel, a nutrient channel that also transports ions and antimalarial compounds at the host erythrocyte membrane. Two blasticidin S-resistant lines exhibited markedly reduced expression of clag genes linked to channel activity, but had no genome-level changes. Silencing aborted production of the channel protein and was directly responsible for reduced uptake. Silencing affected clag paralogs on two chromosomes and was mediated by specific histone modifications, allowing a rapidly reversible drug resistance phenotype advantageous to the parasite. These findings implicate a novel epigenetic resistance mechanism that involves reduced host cell uptake and is a worrisome liability for water-soluble antimalarial drugs.

The avian XPR1 gammaretrovirus receptor is under positive selection and is disabled in bird species in contact with virus-infected wild mice.

Xenotropic mouse leukemia viruses (X-MLVs) are broadly infectious for mammals except most of the classical strains of laboratory mice. These gammaretroviruses rely on the XPR1 receptor for entry, and the unique resistance of laboratory mice is due to two mutations in different putative XPR1 extracellular loops. Cells from avian species differ in susceptibility to X-MLVs, and 2 replacement mutations in the virus-resistant chicken XPR1 (K496Q and Q579E) distinguish it from the more permissive duck and quail receptors. These substitutions align with the two mutations that disable the laboratory mouse XPR1. Mutagenesis of the chicken and duck genes confirms that residues at both sites are critical for virus entry. Among 32 avian species, the 2 disabling XPR1 mutations are found together only in the chicken, an omnivorous, ground-dwelling fowl that was domesticated in India and/or Southeast Asia, which is also where X-MLV-infected house mice evolved. The receptor-disabling mutations are also present separately in 5 additional fowl and raptor species, all of which are native to areas of Asia populated by the virus-infected subspecies Mus musculus castaneus. Phylogenetic analysis showed that the avian XPR1 gene is under positive selection at sites implicated in receptor function, suggesting a defensive role for XPR1 in the avian lineage. Contact between bird species and virus-infected mice may thus have favored selection of mouse virus-resistant receptor orthologs in the birds, and our data suggest that similar receptor-disabling mutations were fixed in mammalian and avian species exposed to similar virus challenges.

B cell immunofocusing and repriming in two HIV-1 Env immunization regimens.

Diverse and rapidly mutating viruses pose challenges to immunogen and vaccine design. In this study, we evaluated the ability of memory B-cells obtained from two independent NHP trials to cross-react with individual HIV-1 vaccine components of two different multivalent immunization strategies. We demonstrated that while an HIV-1 Env multiclade, multivalent immunization regimen resulted in a dominant memory B-cell response that converged toward shared epitopes, in a sequential immunization with clonally-related non-stabilized gp140 HIV-1 Envs followed by SOSIP-stabilized gp140 trimers, the change in immunogen format resulted in repriming of the B-cell response.

Lineage classification and antitubercular drug resistance surveillance of Mycobacterium tuberculosis by whole-genome sequencing in Southern India.

Whole-genome sequencing has created a revolution in tuberculosis management by providing a comprehensive picture of the various genetic polymorphisms with unprecedented accuracy. Studies mapping genomic heterogeneity in clinical isolates of Mycobacterium tuberculosis using a whole-genome sequencing approach from high tuberculosis burden countries are underrepresented. We report whole-genome sequencing results of 242 clinical isolates of culture-confirmed M. tuberculosis isolates from tuberculosis patients referred to a tertiary care hospital in Southern India. Phylogenetic analysis revealed that the isolates in our study belonged to five different lineages, with Indo-Oceanic (lineage 1, n = 122) and East-African Indian (lineage 3, n = 80) being the most prevalent. We report several mutations in genes conferring resistance to first and second line antitubercular drugs including the genes rpoB, katG, ahpC, inhA, fabG1, embB, pncA, rpsL, rrs, and gyrA. The majority of these mutations were identified in relatively high proportions in lineage 1. Our study highlights the utility of whole-genome sequencing as a potential supplemental tool to the existing genotypic and phenotypic methods, in providing expedited comprehensive surveillance of mutations that may be associated with antitubercular drug resistance as well as lineage characterization of M. tuberculosis isolates. Further larger-scale whole-genome datasets with linked minimum inhibition concentration testing are imperative for resolving the discrepancies between whole-genome sequencing and phenotypic drug sensitivity testing results and quantifying the level of the resistance associated with the mutations for optimization of antitubercular drug and precise dose selection in clinics. IMPORTANCE Studies mapping genetic heterogeneity of clinical isolates of M. tuberculosis for determining their strain lineage and drug resistance by whole-genome sequencing are limited in high tuberculosis burden settings. We carried out whole-genome sequencing of 242 M. tuberculosis isolates from drug-sensitive and drug-resistant tuberculosis patients, identified and collected as part of the TB Portals Program, to have a comprehensive insight into the genetic diversity of M. tuberculosis in Southern India. We report several genetic variations in M. tuberculosis that may confer resistance to antitubercular drugs. Further wide-scale efforts are required to fully characterize M. tuberculosis genetic diversity at a population level in high tuberculosis burden settings for providing precise tuberculosis treatment.

A novel specificity protein 1 (SP1)-like gene regulating protein kinase C-1 (Pkc1)-dependent cell wall integrity and virulence factors in Cryptococcus neoformans.

Eukaryotic cells utilize complex signaling systems to detect their environments, responding and adapting as new conditions arise during evolution. The basidiomycete fungus Cryptococcus neoformans is a leading cause of AIDS-related death worldwide and utilizes the calcineurin and protein kinase C-1 (Pkc1) signaling pathways for host adaptation and expression of virulence. In the present studies, a C-terminal zinc finger transcription factor, homologous both to the calcineurin-responsive zinc fingers (Crz1) of ascomycetes and to the Pkc1-dependent specificity protein-1 (Sp1) transcription factors of metazoans, was identified and named SP1 because of its greater similarity to the metazoan factors. Structurally, the Cryptococcus neoformans Sp1 (Cn Sp1) protein was found to have acquired an additional zinc finger motif from that of Crz1 and showed Pkc1-dependent phosphorylation, nuclear localization, and whole genome epistatic associations under starvation conditions. Transcriptional targets of Cn Sp1 shared functional similarities with Crz1 factors, such as cell wall synthesis, but gained the regulation of processes involved in carbohydrate metabolism, including trehalose metabolism, and lost others, such as the induction of autophagy. In addition, overexpression of Cn Sp1 in a pkc1Δ mutant showed restoration of altered phenotypes involved in virulence, including cell wall stability, nitrosative stress, and extracellular capsule production. Cn Sp1 was also found to be important for virulence of the fungus using a mouse model. In summary, these data suggest an evolutionary shift in C-terminal zinc finger proteins during fungal evolution, transforming them from calcineurin-dependent to PKC1-dependent transcription factors, helping to shape the role of fungal pathogenesis of C. neoformans.

Adaptive evolution of Mus Apobec3 includes retroviral insertion and positive selection at two clusters of residues flanking the substrate groove.

Mouse APOBEC3 (mA3) is a cytidine deaminase with antiviral activity. mA3 is linked to the Rfv3 virus resistance factor, a gene responsible for recovery from infection by Friend murine leukemia virus, and mA3 allelic variants differ in their ability to restrict mouse mammary tumor virus. We sequenced mA3 genes from 38 inbred strains and wild mouse species, and compared the mouse sequence and predicted structure with human APOBEC3G (hA3G). An inserted sequence was identified in the virus restrictive C57BL strain allele that disrupts a splice donor site. This insertion represents the long terminal repeat of the xenotropic mouse gammaretrovirus, and was acquired in Eurasian mice that harbor xenotropic retrovirus. This viral regulatory sequence does not alter splicing but is associated with elevated mA3 expression levels in spleens of laboratory and wild-derived mice. Analysis of Mus mA3 coding sequences produced evidence of positive selection and identified 10 codons with very high posterior probabilities of having evolved under positive selection. Six of these codons lie in two clusters in the N-terminal catalytically active cytidine deaminase domain (CDA), and 5 of those 6 codons are polymorphic in Rfv3 virus restrictive and nonrestrictive mice and align with hA3G CDA codons that are critical for deaminase activity. Homology models of mA3 indicate that the two selected codon clusters specify residues that are opposite each other along the predicted CDA active site groove, and that one cluster corresponds to an hAPOBEC substrate recognition loop. Substitutions at these clustered mA3 codons alter antiviral activity. This analysis suggests that mA3 has been under positive selection throughout Mus evolution, and identified an inserted retroviral regulatory sequence associated with enhanced expression in virus resistant mice and specific residues that modulate antiviral activity.

Origin, antiviral function and evidence for positive selection of the gammaretrovirus restriction gene Fv1 in the genus Mus.

The Fv1 virus resistance gene is a coopted endogenous retrovirus (ERV) sequence related to the gag gene of the MuERV-L ERV family. Three major Fv1 resistance alleles have been identified in laboratory mice, and they target virus capsid genes to produce characteristic patterns of resistance to mouse leukemia viruses (MLVs). We identified Fv1 in 3 of the 4 Mus subgenera; its absence from Coelomys and 1 of 3 species of Pyromys indicate Fv1 was acquired shortly after the origin of the Mus genus. We sequenced Fv1 genes from 21 mice representative of the major taxonomic groups of Mus. Two lines of evidence indicate that Fv1 has had antiviral function for 7 million years of evolution. First, 2 species of African pygmy mice (subgenus Nannomys) show an Fv1-like MLV resistance, and transduced cells expressing the Nannomys Fv1 gene reproduce this resistance pattern. Second, sequence comparisons suggest that Fv1 has been involved in genetic conflicts throughout Mus evolution. We found evidence for strong positive selection of Fv1 and identified 6 codons that show evidence of positive selection: 3 codons in the C-terminal region including 2 previously shown to contribute to Fv1 restriction in laboratory mice, and 3 codons in a 10-codon segment overlapping the major homology region of Fv1; this segment is known to be involved in capsid multimerization. This analysis suggests that Fv1 has had an antiviral role throughout Mus evolution predating exposure of mice to the MLVs restricted by laboratory mouse Fv1, and suggests a mechanism for Fv1 restriction.

Patterns of genomic interrelatedness of publicly available samples in the TB portals database.

The TB Portals program is an international collaboration for the collection and dissemination of tuberculosis data from patient cases focused on drug resistance. The central database is a patient-oriented resource containing both patient and pathogen clinical and genomic information. Herein we provide a summary of the pathogen genomic data available through the TB Portals and show one potential application by examining patterns of genomic pairwise distances. Distributions of pairwise distances highlight overall patterns of genome variability within and between Mycobacterium tuberculosis phylogenomic lineages. Closely related isolates (based on whole-genome pairwise distances and time between sample collection dates) from different countries were identified as potential evidence of international transmission of drug-resistant tuberculosis. These high-level views of genomic relatedness provide information that can stimulate hypotheses for further and more detailed research.

Evolution of functional and sequence variants of the mammalian XPR1 receptor for mouse xenotropic gammaretroviruses and the human-derived retrovirus XMRV.

Genetic conflicts between retroviruses and their receptors result in the evolution of novel host entry restrictions and novel virus envelopes, and such variants can influence trans-species transmission. We screened rodents and other mammals for sequence variation in the Xpr1 receptor for the mouse xenotropic or polytropic mouse leukemia viruses (X-MLVs or P-MLVs, respectively) of the gammaretrovirus family and for susceptibility to mouse-derived X/P-MLVs and to XMRV (xenotropic murine leukemia virus-related virus), an X-MLV-like virus isolated from humans with prostate cancer and chronic fatigue syndrome. We identified multiple distinct susceptibility phenotypes; these include the four known Xpr1 variants in Mus and a novel fifth Xpr1 gene found in Mus molossinus and Mus musculus. We describe the geographic and species distribution of the Mus Xpr1 variants but failed to find the X-MLV-restrictive laboratory mouse allele in any wild mouse. We used mutagenesis and phylogenetic analysis to evaluate the functional contributions made by constrained, variable, and deleted residues. Rodent Xpr1 is under positive selection, indicating a history of host-pathogen conflicts; several codons under selection have known roles in virus entry. All non-Mus mammals are susceptible to mouse X-MLVs, but some restrict other members of the X/P-MLV family, and the resistance of hamster and gerbil cells to XMRV indicates that XMRV has unique receptor requirements. We show that the hypervariable fourth extracellular XPR1 loop (ECL4) contains three evolutionarily constrained residues that do not contribute to receptor function, we identify two novel residues important for virus entry (I579 and T583), and we describe a unique pattern of ECL4 variation in the three virus-restrictive Xpr1 variants found in MLV-infected house mice; these mice carry different deletions in ECL4, suggesting either that these sites or loop size affects receptor function.

Patterns of Coevolutionary Adaptations across Time and Space in Mouse Gammaretroviruses and Three Restrictive Host Factors.

The classical laboratory mouse strains are genetic mosaics of three Mus musculus subspecies that occupy distinct regions of Eurasia. These strains and subspecies carry infectious and endogenous mouse leukemia viruses (MLVs) that can be pathogenic and mutagenic. MLVs evolved in concert with restrictive host factors with some under positive selection, including the XPR1 receptor for xenotropic/polytropic MLVs (X/P-MLVs) and the post-entry restriction factor Fv1. Since positive selection marks host-pathogen genetic conflicts, we examined MLVs for counter-adaptations at sites that interact with XPR1, Fv1, and the CAT1 receptor for ecotropic MLVs (E-MLVs). Results describe different co-adaptive evolutionary paths within the ranges occupied by these virus-infected subspecies. The interface of CAT1, and the otherwise variable E-MLV envelopes, is highly conserved; antiviral protection is afforded by the Fv4 restriction factor. XPR1 and X/P-MLVs variants show coordinate geographic distributions, with receptor critical sites in envelope, under positive selection but with little variation in envelope and XPR1 in mice carrying P-ERVs. The major Fv1 target in the viral capsid is under positive selection, and the distribution of Fv1 alleles is subspecies-correlated. These data document adaptive, spatial and temporal, co-evolutionary trajectories at the critical interfaces of MLVs and the host factors that restrict their replication.