Signature pattern analysis: a method for assessing viral sequence relatedness

Evolution and genetic characterization of Seoul virus in wild rats Rattus norvegicus from an urban park in Lyon, France 2020-2022

BACKGROUND

Seoul virus (SEOV) is an orthohantavirus primarily carried by rats. In humans, it may cause hemorrhagic fever with renal syndrome (HFRS). Its incidence is likely underestimated and given the expansion of urban areas, a better knowledge of SEOV circulation in rat populations is called for. Beyond the need to improve human case detection, we need to deepen our comprehension of the ecological, epidemiological, and evolutionary processes involved in the transmission of SEOV.

METHODOLOGY / PRINCIPAL FINDINGS

We performed a comprehensive serological and molecular characterization of SEOV in Rattus norvegicus in a popular urban park within a large city (Lyon, France) to provide essential information to design surveillance strategies regarding SEOV. We sampled rats within the urban park of 'La Tête d'Or' in Lyon city from 2020 to 2022. We combined rat population genetics, immunofluorescence assays, SEOV high-throughput sequencing (S, M, and L segments), and phylogenetic analyses. We found low structuring of wild rat populations within Lyon city. Only one sampling site within the park (building created in 2021) showed high genetic differentiation and deserves further attention. We confirmed the circulation of SEOV in rats from the park with high seroprevalence (17.2%) and high genetic similarity with the strain previously described in 2011 in Lyon city.

CONCLUSION/SIGNIFICANCE

This study confirms the continuous circulation of SEOV in a popular urban park where the risk for SEOV transmission to humans is present. Implementing a surveillance of this virus could provide an efficient early warning system and help prepare risk-based interventions. As we reveal high gene flow between rat populations from the park and the rest of the city, we advocate for SEOV surveillance to be conducted at the scale of the entire city.

BK Polyomavirus Diversity After Hematopoietic Stem Cell Transplantation

BK polyomavirus (BKPyV) infection is common after hematopoietic stem cell transplantation (HSCT) and is associated with the development of hemorrhagic cystitis (HC). The role that BKPyV plays in the pathogenesis of HC is not well characterized. We investigated the impact of BKPyV diversity on the development of HC using a previously established cohort of pediatric HSCT patients. There were 147 urine samples with quantifiable BKPyV at month 1 after HSCT; 137 (93.2%) were amplified using our in-house polymerase chain reaction approach and sent for next-generation sequencing. Subtype Ia was most frequent (61.3%), followed by subtype Ib1 (31.4%). The median viral load of subtype Ia samples was higher than for subtype Ib1 at month 1. Across the protein coding regions, APOBEC-induced mutations and signature patterns associated with HC were identified. This is the largest sequencing study of a single cohort of HSCT patients, providing a vast resource of sequence data for future analyses.

HIV-1 subtypes maintain distinctive physicochemical signatures in Nef domains associated with immunoregulation

BACKGROUND

HIV subtype is associated with varied rates of disease progression. The HIV accessory protein, Nef, continues to be present during antiretroviral therapy (ART) where it has numerous immunoregulatory effects. In this study, we analyzed Nef sequences from HIV subtypes A1, B, C, and D using a machine learning approach that integrates functional amino acid information to identify if unique physicochemical features are associated with Nef functional/structural domains in a subtype-specific manner.

METHODS

2253 sequences representing subtypes A1, B, C, and D were aligned and domains with known functional properties were scored based on amino acid physicochemical properties. Following feature generation, we used statistical pruning and evolved neural networks (ENNs) to determine if we could successfully classify subtypes. Next, we used ENNs to identify the top five key Nef physicochemical features applied to specific immunoregulatory domains that differentiated subtypes. A signature pattern analysis was performed to the assess amino acid diversity in sub-domains that differentiated each subtype.

RESULTS

In validation studies, ENNs successfully differentiated each subtype at A1 (87.2%), subtype B (89.5%), subtype C (91.7%), and subtype D (85.1%). Our feature-based domain scoring, followed by t-tests, and a similar ENN identified subtype-specific domain-associated features. Subtype A1 was associated with alterations in Nef CD4 binding domain; subtype B was associated with alterations with the AP-2 Binding domain; subtype C was associated with alterations in a structural Alpha Helix domain; and, subtype D was associated with alterations in a Beta-Sheet domain.

CONCLUSIONS

Recent studies have focused on HIV Nef as a driver of immunoregulatory disease in those HIV infected and on ART. Nef acts through a complex mixture of interactions that are directly linked to the key features of the subtype-specific domains we identified with the ENN. The study supports the hypothesis that varied Nef subtypes contribute to subtype-specific disease progression.

Infectivity of three Mayaro Virus geographic isolates in human cell lines

Mayaro virus (MAYV) is an emergent arthropod-borne virus that causes an acute febrile illness accompanied by arthralgia, similar to chikungunya virus. Increasing urbanization of MAYV outbreaks in the Americas has led to concerns for geographic expansion and spillover. Given the potential importance of this pathogen, we sought to fill critical gaps in knowledge regarding MAYV infectivity and geographic variation. This study describes the cytopathogenicity of MAYV in human dermal fibroblasts, human skeletal muscle satellite cells, human embryonic kidney cells (HEK), peripherally derived human macrophages, and Vero cells. We found that regional differences between these viruses do not affect replication kinetics, with high titers peaking at 37 h post infection. MAYV-U, did however, cause the most cytopathic effect in a time-dependent manner. Compared to the other two prototypic isolates, MAYV-U harbors unique mutations in the E2 protein, D60G and S205F, that are likely to interact with the host cell receptor and could affect infectivity. We further demonstrate that pre-treatment of cells with interferon-β inhibited viral replication in a dose-dependent manner. Together, these findings advance our understanding of MAYV infection of human target cells and provide initial data regarding variation according to geography.

Intrahost evolution of the HIV-2 capsid correlates with progression to AIDS

HIV-2 infection will progress to AIDS in most patients without treatment, albeit at approximately half the rate of HIV-1 infection. HIV-2 capsid (p26) amino acid polymorphisms are associated with lower viral loads and enhanced processing of T cell epitopes, which may lead to protective Gag-specific T cell responses common in slower progressors. Lower virus evolutionary rates, and positive selection on conserved residues in HIV-2 env have been associated with slower progression to AIDS. In this study we analysed 369 heterochronous HIV-2 p26 sequences from 12 participants with a median age of 30 years at enrolment. CD4% change over time was used to stratify participants into relative faster and slower progressor groups. We analysed p26 sequence diversity evolution, measured site-specific selection pressures and evolutionary rates, and determined if these evolutionary parameters were associated with progression status. Faster progressors had lower CD4% and faster CD4% decline rates. Median pairwise sequence diversity was higher in faster progressors (5.7x10-3 versus 1.4x10-3 base substitutions per site, P<0.001). p26 evolved under negative selection in both groups (dN/dS=0.12). Median virus evolutionary rates were higher in faster than slower progressors - synonymous rates: 4.6x10-3 vs. 2.3x10-3; and nonsynonymous rates: 6.9x10-4 vs. 2.7x10-4 substitutions/site/year, respectively. Virus evolutionary rates correlated negatively with CD4% change rates (ρ = -0.8, P=0.02), but not CD4% level. The signature amino acid at p26 positions 6, 12 and 119 differed between faster (6A, 12I, 119A) and slower (6G, 12V, 119P) progressors. These amino acid positions clustered near to the TRIM5α/p26 hexamer interface surface. p26 evolutionary rates were associated with progression to AIDS and were mostly driven by synonymous substitutions. Nonsynonymous evolutionary rates were an order of magnitude lower than synonymous rates, with limited amino acid sequence evolution over time within hosts. These results indicate HIV-2 p26 may be an attractive therapeutic target.

Seronegative human T-cell lymphotropic virus 1 carriers in blood banks: A potential viral source for silent transmission?

BACKGROUND AND OBJECTIVES

Transfusion-transmitted viruses count among the greatest threats to blood safety. In Argentina, current laws oblige testing all donated blood for the presence of antibodies against human T-cell lymphotropic viruses 1 and 2 (HTLV-1/2). In endemic zones of the country, a high rate of seronegative HTLV-1 individuals with clear evidence of infection because of symptoms and/or presence of tax sequences of HTLV-1 and/or IgG anti-Tax antibodies has been recently described. Migration from endemic to nonendemic zones of Argentina is very frequent.

MATERIALS AND METHODS

During a 1-year period, in the blood bank of Córdoba city, we performed molecular screening of all donors who were born in or arose from endemic zones for HTLV-1/2 in Argentina and neighbouring countries.

RESULTS

By screening 219 bp of HTLV-1/2 tax gene, 0.6% (2/317) of the blood donors proved to be positive for HTLV-1 tax sequence. One of the donors presented anti-Tax antibodies, demonstrating the transcriptional activity of the tax gene, and the other donor was also positive for LTR and pol gene sequences. The HTLV-1 genetic analysis of the LTR sequence determined that it belonged to the Cosmopolitan subtype HTLV-1aA.

CONCLUSION

These findings suggest potential limitations of some currently approved screening assays for HTLV-1 detection applied in some donor populations and the possibility of an HTLV-1 seronegative carrier state with the potential for silent transmission by blood.

HIV-1 subtype C Tat exon-1 amino acid residue 24K is a signature for neurocognitive impairment

Variation and differential selection pressures on Tat genes have been shown to alter the biological function of the protein, resulting in pathological consequences in a number of organs including the brain. We evaluated the impact of genetic variation and selection pressure on 147 HIV-1 subtype C Tat exon 1 sequences from monocyte-depleted peripheral lymphocytes on clinical diagnosis of neurocognitive impairment. Genetic analyses identified two signature amino acid residues, lysine at codon 24 (24K) with a frequency of 43.4% and arginine at codon 29 (29R) with a frequency of 34.0% in individuals with HIV-associated neurocognitive impairment. The analyses also revealed two signature residues, asparagine, 24 N (31.9%), and histidine, 29H (21.3%), in individuals without neurocognitive impairment. Both codons, 24 and 29, were associated with high entropy but only codon 29 was under positive selection. The presence of signature K24 increased by 2.08 times the risk of neurocognitive impairment, 3.15 times higher proviral load, and 69% lower absolute CD4 T-cell count compared to those without the signature. The results support a linkage between HIV-1 C Tat N24K polymorphism, proviral load, immunosuppression, and neurocognitive impairment. The signature may induce more neurotoxic effects, which contributes to establishment and severity of HIV-associated neurocognitive impairment.

Origin and evolution of HIV-1 subtype A6

Background

HIV outbreaks in the Former Soviet Union (FSU) countries were characterized by repeated transmission of the HIV variant AFSU, which is now classified as a distinct subtype A sub-subtype called A6. The current study used phylogenetic/phylodynamic and signature mutation analyses to determine likely evolutionary relationship between subtype A6 and other subtype A sub-subtypes.

Methods

For this study, an initial Maximum Likelihood phylogenetic analysis was performed using a total of 553 full-length, publicly available, reverse transcriptase sequences, from A1, A2, A3, A4, A5, and A6 sub-subtypes of subtype A. For phylogenetic clustering and signature mutation analysis, a total of 5961 and 3959 pol and env sequences, respectively, were used.

Results

Phylogenetic and signature mutation analysis showed that HIV-1 sub-subtype A6 likely originated from sub-subtype A1 of African origin. A6 and A1 pol and env genes shared several signature mutations that indicate genetic similarity between the two subtypes. For A6, tMRCA dated to 1975, 15 years later than that of A1.

Conclusion

The current study provides insights into the evolution and diversification of A6 in the backdrop of FSU countries and indicates that A6 in FSU countries evolved from A1 of African origin and is getting bridged outside the FSU region.

Sequence Length of HIV-1 Subtype B Increases over Time: Analysis of a Cohort of Patients with Hemophilia over 30 Years

We aimed to investigate whether the sequence length of HIV-1 increases over time. We performed a longitudinal analysis of full-length coding region sequences (FLs) during an HIV-1 outbreak among patients with hemophilia and local controls infected with the Korean subclade B of HIV-1 (KSB). Genes were amplified by overlapping RT-PCR or nested PCR and subjected to direct sequencing. Overall, 141 FLs were sequentially determined over 30 years in 62 KSB-infected patients. Phylogenetic analysis indicated that within KSB, two FLs from plasma donors O and P comprised two clusters, together with 8 and 12 patients with hemophilia, respectively. Signature pattern analysis of the KSB of HIV-1 revealed 91 signature nucleotide residues (1.1%). In total, 48 and 43 signature nucleotides originated from clusters O and P, respectively. Six positions contained 100% specific nucleotide(s) in clusters O and P. In-depth FL analysis for over 30 years indicated that the KSB FL significantly increased over time before combination antiretroviral therapy (cART) and decreased with cART. This increase occurred due to the significant increase in env and nef genes, originating in the variable regions of both genes. The increase in sequence length of HIV-1 over time suggests an evolutionary direction.

Baseline Hepatitis C Virus NS5A Resistance-Associated Polymorphisms in Patients With and Without Human Immunodeficiency Virus Coinfection in Mexico

Objective: We aimed to evaluate the frequency and associated factors of baseline NS5A resistance-associated substitutions (RASs) in patients coinfected with human immunodeficiency virus (HIV) and hepatitis C virus (HCV) monoinfection with genotype 1b (GT1b) or genotype 1a (GT1a). Moreover, we performed a phylogenetic analysis to evaluate the pattern of clustering among samples of patients with RASs. Results: Fifty-five patients were infected with GT1a, of whom 44 (80%) were HIV-infected patients. RAS prevalence in GT1a was 14% (6/44) and distributed as follows: 5 (11%) harbored M28V and 1 (2%) A92T. Twenty-four patients were infected with HCV GT1b, of whom only 5 (21%) were HIV coinfected; RASs were found in 17/24 (71%) patients, as follows: Y93H+F37L+Q54H (1/24), Y93H+F37L (1/24), P58S (1/24), L31F+F37L (1/24), F37L+H/Q54H (3/24), and F37L (10/24). Only GT1b was significantly associated with RASs (adjusted odds ratio 16.37; 95% confidence interval 2.74-97.48; p = 0.002) in the multivariate analysis. A cluster of sequences from HIV/HCV GT1a patients was found; however, we did not find phylogenetic relationships among sequences with NS5A RASs. Conclusions: In our population of HCV-infected patients, the frequency of NS5A RASs at baseline was somewhat similar to the previously reported worldwide rate. HCV GT1b showed the most significant association with harboring of NS5A RASs. Of note, despite there being clusters among sequences of HIV-coinfected patients, NS5A RASs were not transmitted.

Emergence of porcine delta-coronavirus pathogenic infections among children in Haiti through independent zoonoses and convergent evolution

Coronaviruses have caused three major epidemics since 2003, including the ongoing SARS-CoV-2 pandemic. In each case, coronavirus emergence in our species has been associated with zoonotic transmissions from animal reservoirs^1,2, underscoring how prone such pathogens are to spill over and adapt to new species. Among the four recognized genera of the family Coronaviridae – Alphacoronavirus, Betacoronavirus, Deltacoronavirus, Gammacoronavirus, – human infections reported to date have been limited to alpha- and betacoronaviruses³. We identify, for the first time, porcine deltacoronavirus (PDCoV) strains in plasma samples of three Haitian children with acute undifferentiated febrile illness. Genomic and evolutionary analyses reveal that human infections were the result of at least two independent zoonoses of distinct viral lineages that acquired the same mutational signature in the nsp15 and the spike glycoprotein genes by convergent evolution. In particular, structural analysis predicts that one of the changes in the Spike S1 subunit, which contains the receptor-binding domain, may affect protein’s flexibility and binding to the host cell receptor. Our findings not only underscore the ability of deltacoronaviruses to adapt and potentially lead to human-to-human transmission, but also raise questions about the role of such transmissions in development of pre-existing immunity to other coronaviruses, such as SARS-CoV-2.

Intra-Clade C signature polymorphisms in HIV-1 LTR region: The Indian and African lookout

Polymorphisms occurring in LTR (Long Terminal Repeat) region can profoundly impact pathogenicity, transmission and biology of Human Immunodeficiency Virus Type 1 (HIV-1). We investigated intra-clade polymorphisms, associated with HIV-1 clade-C infections that occur in India and Africa. Plasma samples were obtained from 24 HIV-infected ART-experienced individuals. Next Generation Sequencing was performed on Illumina Hi Seq X system. Sequence analysis was done using MEGA v7. Transcription factor binding sites (TFBS) were investigated to unveil signature sequences. Signature nucleotides in Indian sequences were observed at 19 positions, of which 7 nucleotide signatures occurred in transcription binding sites (TFBS), namely NF-AT-II, NF-AT-III, USF, TCF- 1alpha, Sp1-I and TAR. Intra-clade C variations in HIV-1 LTR that inscribe signature nucleotides in Indian sequences lead to formation monophyletic cluster of Indian sequences. Moreover, occurrence of intra-clade signature nucleotides was observed at the key positions in the transcription factor binding sites in Indian and African clade-C sequences.

Viral protein R polymorphisms in the pathogenesis of HIV-associated acute ischaemic stroke: a case-control study

HIV-1 viral proteins have been implicated in endothelial dysfunction, which is a major determinant of ischaemic stroke risk in HIV-infected individuals. Polymorphisms in HIV-1 viral protein R (Vpr) may alter its potential to promote endothelial dysfunction, by modifying its effects on viral replication, reactivation of latent cells, upregulation of pro-inflammatory cytokines and infection of macrophages. We analysed Vpr polymorphisms and their association with acute ischaemic stroke by comparing Vpr signature amino acids between 54 HIV-infected individuals with acute ischaemic stroke, and 80 age-matched HIV-infected non-stroke controls. Isoleucine at position 22 and serine at position 41 were associated with ischaemic stroke in HIV. Individuals with stroke had lower CD4 counts and CD4 nadirs than controls. These polymorphisms are unique to individuals with stroke compared to South African subtype C and the control group consensus sequences. Signature Vpr polymorphisms are associated with acute ischaemic stroke in HIV. These may increase stroke risk by promoting endothelial dysfunction and susceptibility to opportunistic infections. Therapeutic targeting of HIV-1 viral proteins may present an additional mechanism of decreasing stroke risk in HIV-infected individuals.

Non-synonymous Substitutions in HIV-1 GAG Are Frequent in Epitopes Outside the Functionally Conserved Regions and Associated With Subtype Differences

In 2019, 38 million people lived with HIV-1 infection resulting in 690,000 deaths. Over 50% of this infection and its associated deaths occurred in Sub-Saharan Africa. The West African region is a known hotspot of the HIV-1 epidemic. There is a need to develop an HIV-1 vaccine if the HIV epidemic would be effectively controlled. Few protective cytotoxic T Lymphocytes (CTL) epitopes within the HIV-1 GAG (HIV_gagconsv) have been previously identified to be functionally conserved among the HIV-1 M group. These epitopes are currently the focus of universal HIV-1 T cell-based vaccine studies. However, these epitopes' phenotypic and genetic properties have not been observed in natural settings for HIV-1 strains circulating in the West African region. This information is critical as the usefulness of universal HIV-1 vaccines in the West African region depends on these epitopes' occurrence in strains circulating in the area. This study describes non-synonymous substitutions within and without HIV_gagconsv genes isolated from 10 infected Nigerians at the early stages of HIV-1 infection. Furthermore, we analyzed these substitutions longitudinally in five infected individuals from the early stages of infection till after seroconversion. We identified three non-synonymous substitutions within HIV_gagconsv genes isolated from early HIV infected individuals. Fourteen and nineteen mutations outside the HIV_gagconsv were observed before and after seroconversion, respectively, while we found four mutations within the HIV_gagconsv. These substitutions include previously mapped CTL epitope immune escape mutants. CTL immune pressure likely leaves different footprints on HIV-1 GAG epitopes within and outside the HIV_gagconsv. This information is crucial for universal HIV-1 vaccine designs for use in the West African region.

Independent infections of porcine deltacoronavirus among Haitian children

Coronaviruses have caused three major epidemics since 2003, including the ongoing SARS-CoV-2 pandemic. In each case, the emergence of coronavirus in our species has been associated with zoonotic transmissions from animal reservoirs^1,2, underscoring how prone such pathogens are to spill over and adapt to new species. Among the four recognized genera of the family Coronaviridae, human infections reported so far have been limited to alphacoronaviruses and betacoronaviruses^3-5. Here we identify porcine deltacoronavirus strains in plasma samples of three Haitian children with acute undifferentiated febrile illness. Genomic and evolutionary analyses reveal that human infections were the result of at least two independent zoonoses of distinct viral lineages that acquired the same mutational signature in the genes encoding Nsp15 and the spike glycoprotein. In particular, structural analysis predicts that one of the changes in the spike S1 subunit, which contains the receptor-binding domain, may affect the flexibility of the protein and its binding to the host cell receptor. Our findings highlight the potential for evolutionary change and adaptation leading to human infections by coronaviruses outside of the previously recognized human-associated coronavirus groups, particularly in settings where there may be close human-animal contact.

Implications of viral transmitted/founder (T/F) dynamics on vaccine development

Viral infection typically originates from a limited number of virions known as transmitted/founder (T/F) viruses. Studies of cross-species transmission, and intra-species transmission of antigenically variable viruses, indicates T/F variants may express distinct, transmissibility enhancing phenotypes. However, with evidence that transmissibility is associated with not only intrinsic virological features, such as virion composition, but also extrinsic factors, such as viral population structure, the challenge of resolving T/F signatures that can be targeted by rational vaccine or antiviral design is substantial. Nonetheless, failure to develop vaccines for antigenically variable viruses, such as HIV/HCV, and the ongoing risk of cross-species transmission with pandemic potential, recommends development of T/F targeting vaccines. In this commentary, the T/F phenomena is introduced, explored in both the classical (HIV) and non-canonical (coronaviruses) instances, and discussed in relation to rational and preemptive vaccine design.

Multiple Recombination Events and Strong Purifying Selection at the Origin of SARS-CoV-2 Spike Glycoprotein Increased Correlated Dynamic Movements

Our evolutionary and structural analyses revealed that the severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2) spike gene is a complex mosaic resulting from several recombination events. Additionally, the fixation of variants has mainly been driven by purifying selection, suggesting the presence of conserved structural features. Our dynamic simulations identified two main long-range covariant dynamic movements of the novel glycoprotein, and showed that, as a result of the evolutionary duality, they are preserved. The first movement involves the receptor binding domain with the N-terminal domain and the C-terminal domain 2 and is maintained across human, bat and pangolin coronaviruses. The second is a complex network of long-range dynamics specific to SARS-CoV-2 involving the novel PRRA and the conserved KR*SF cleavage sites, as well as conserved segments in C-terminal domain 3. These movements, essential for host cell binding, are maintained by hinges conserved across human, bat, and pangolin coronaviruses glycoproteins. The hinges, located around Threonine 333 and Proline 527 within the N-terminal domain and C-terminal domain 2, represent candidate targets for the future development of novel pan-coronavirus inhibitors. In summary, we show that while recombination created a new configuration that increased the covariant dynamic movements of the SARS-CoV-2 glycoprotein, negative selection preserved its inter-domain structure throughout evolution in different hosts and inter-species transmissions.

HIV-1 gp41 genetic diversity and enfuvirtide resistance-associated mutations among enfuvirtide-naïve patients in southern China

BACKGROUND

Human immunodeficiency virus (HIV) increasing molecular diversity and emergence of drug resistant mutants remain a major concern in China. Enfuvirtide (ENF/T-20) is the first entry inhibitor used in patients failing highly active antiretroviral therapy (HAART). However, data on HIV-1 gp41genetic diversity and primary ENF resistance-associated mutations among treatment-naïve patients in China is limited. The objective was to identify molecular diversity and ENF resistance patterns of HIV-1 in southern China, using envelope (env) gp41 sequences and bioinformatics tools, which may help optimize antiretroviral therapy.

METHODS

From November 2018 to January 2019, 439 blood plasma samples from ENF-naïve patients were collected from Shenzhen, Wuhan and Chongqing, of which 396 HIV env regions were sequenced and subtyped, and were performed the analysis of drug resistance-associated mutations (DRMs).

RESULTS

The main subtypes were circulating recombinant form (CRF) 01_AE (30.6 %) and CRF07_BC (48.7 %). CRF55_01B had been the fourth subtype in the study, and many rare CRFs were observed. Notably, CRF02_AG and CRF_BF strains typically found in Africa and US respectively were identified amongst Chinese patients. Known DRMs were detected in 27.5 % (109/396) of ENF treatment-naïve patients. One major DRM (L44 M), many secondary DRMs (including N126 K, E137 K, S138A), and lots of polymorphisms were found in the study, which have been proved to elevate resistance to ENF.

CONCLUSIONS

HIV-1 molecular diversity was observed in the study, which indicating that HIV-1 variability is becoming increasingly complex in southern China. A diverse set of primary DRMs discovered in this study described the serious threat to ART, which reminds us the urgent need of timely surveillance of HIV-1 viral diversity and drug resistance in China.

The Signature Amino Acid Residue Serine 31 of HIV-1C Tat Potentiates an Activated Phenotype in Endothelial Cells

The natural cysteine to serine variation at position 31 of Tat in HIV-1C disrupts the dicysteine motif attenuating the chemokine function of Tat. We ask if there exists a trade-off in terms of a gain of function for HIV-1C Tat due to this natural variation. We constructed two Tat-expression vectors encoding Tat proteins discordant for the serine 31 residue (CS-Tat vs. CC-Tat), expressed the proteins in Jurkat cells under doxycycline control, and performed the whole transcriptome analysis to compare the early events of Tat-induced host gene expression. Our analysis delineated a significant enrichment of pathways and gene ontologies associated with the angiogenic signaling events in CS-Tat stable cells. Subsequently, we validated and compared angiogenic signaling events induced by CS- vs. CC-Tat using human umbilical vein endothelial cells (HUVEC) and the human cerebral microvascular endothelial cell line (hCMEC/D3). CS-Tat significantly enhanced the production of CCL2 from HUVEC and induced an activated phenotype in endothelial cells conferring on them enhanced migration, invasion, and in vitro morphogenesis potential. The ability of CS-Tat to induce the activated phenotype in endothelial cells could be of significance, especially in the context of HIV-associated cardiovascular and neuronal disorders. The findings from the present study are likely to help appreciate the functional significance of the SAR (signature amino acid residues) influencing the unique biological properties.

D. de Azevedo S, Delatorre E, Bello G. Few amino acid signatures distinguish HIV-1 subtype B pandemic and non-pandemic strains

The Human Immunodeficiency Virus Type I (HIV-1) subtype B comprises approximately 10% of all HIV infections in the world. The HIV-1 subtype B epidemic comprehends a pandemic variant (named B_PANDEMIC) disseminated worldwide and non-pandemic variants (named B_CAR) that are mostly restricted to the Caribbean. The goal of this work was the identification of amino acid signatures (AAs) characteristic to the B_CAR and B_PANDEMIC variants. To this end, we analyzed HIV-1 subtype B full-length (n = 486) and partial (n = 814) genomic sequences from the Americas classified within the B_CAR and B_PANDEMIC clades and reconstructed the sequences of their most recent common ancestors (MRCA). Analysis of contemporary HIV-1 sequences revealed 13 AAs between B_CAR and B_PANDEMIC variants (four on Gag, three on Pol, three on Rev, and one in Vif, Vpu, and Tat) of which only two (one on Gag and one on Pol) were traced to the MRCA. All AAs correspond to polymorphic sites located outside essential functional proteins domains, except the AAs in Tat. The absence of stringent AAs inherited from their ancestors between modern B_CAR and B_PANDEMIC variants support that ecological factors, rather than viral determinants, were the main driving force behind the successful spread of the B_PANDEMIC strain.

Detection and Genetic Characterization of Puumala Orthohantavirus S-Segment in Areas of France Non-Endemic for Nephropathia Epidemica

Puumala virus (PUUV) in Europe causes nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome (HFRS). The incidence of NE is highly heterogeneous spatially, whereas the geographic distribution of the wild reservoir of PUUV, the bank vole, is essentially homogeneous. Our understanding of the processes driving this heterogeneity remains incomplete due to gaps in knowledge. Little is known about the current distribution and genetic variation of PUUV in the areas outside the well-identified zones of NE endemicity. We trapped bank voles in four forests in French regions in which NE is considered non-endemic, but sporadic NE cases have been reported recently. We tested bank voles for anti-PUUV IgG and characterized the S segment sequences of PUUV from seropositive animals. Phylogenetic analyses revealed specific amino-acid signatures and genetic differences between PUUV circulating in non-endemic and nearby NE-endemic areas. We also showed, in temporal surveys, that the amino-acid sequences of PUUV had undergone fewer recent changes in areas non-endemic for NE than in endemic areas. The evolutionary history of the current French PUUV clusters was investigated by phylogeographic approaches, and the results were considered in the context of the history of French forests. Our findings highlight the need to monitor the circulation and genetics of PUUV in a larger array of bank vole populations, to improve our understanding of the risk of NE.

Diversity of the hepatitis C virus NS5B gene during HIV co-infection

Viral diversity is an important feature of hepatitis C virus (HCV) infection and an important predictor of disease progression and treatment response. HIV/HCV co-infection is associated with enhanced HCV replication, increased fibrosis, and the development of liver disease. HIV also increases quasispecies diversity of HCV structural genes, although limited data are available regarding the impact of HIV on non-structural genes of HCV, particularly in the absence of direct-acting therapies. The genetic diversity and presence of drug resistance mutations within the RNA-dependent RNA polymerase (NS5B) gene were examined in 3 groups of women with HCV genotype 1a infection, including those with HCV mono-infection, antiretroviral (ART)-naïve women with HIV/HCV co-infection and CD4 cell count <350 cells/mm3, and ART-naïve women with HIV/HCV co-infection and CD4 cell count ≥350 cells/mm3. None had ever been treated for HCV infection. There was evidence of significant diversity across the entire NS5B gene in all women. There were several nucleotides and amino acids with distinct distributions across the three study groups, although no obvious clustering of NS5B sequences was observed based on HIV co-infection or CD4 cell count. Polymorphisms at amino acid positions associated with resistance to dasabuvir and sofosbuvir were limited, although the Q309R variant associated with ribavirin resistance was present in 12 individuals with HCV mono-infection, 8 HIV/HCV co-infected individuals with CD4 <350 cells/mm3, and 12 HIV/HCV co-infected individuals with CD4 ≥350 cells/mm3. Previously reported fitness altering mutations were rare. CD8+ T cell responses against the human leukocyte antigen (HLA) B57-restricted epitopes NS5B2629-2637 and NS5B2936-2944 are critical for HCV control and were completely conserved in 44 (51.8%) and 70 (82.4%) study participants. These data demonstrate extensive variation across the NS5B gene. Genotypic variation may have a profound impact on HCV replication and pathogenesis and deserves careful evaluation.

Lack of association between amino acid sequences of the bovine leukemia virus envelope and varying stages of infection in dairy cattle

We collected 724 blood samples from dairy cattle from six Mexican states, and tested them for the presence of antibodies against BLV using a commercial ELISA test. Our study groups consisted of 32 samples: 12 asymptomatic cows, 12 cows with lymphocytosis and 8 samples of tumor tissue of the abomasum and heart of cattle with lymphoma. We designed three pairs of primers to amplify the complete BLV env gene, and obtained a fragment of 1548 nucleotides in length with the sequenced products. According to the phylogenetic tree we constructed to identify the viral genotype, 96.87 % of the sequences grouped into genotype 1, while a single sample from a cow with lymphocytosis (3.13 %) was associated with genotype 3 sequences. The similarity between the Mexican BLV sequences ranged from 0.985-1.00. In addition, the proportion of non-synonymous and synonymous mutations indicated negative selection. We did not identify any conserved residues in the viral protein sequences that could be related to BLV infection stage in cattle. Proviral quantification was performed using quantitative polymerase chain reaction, and we used Mood´s median test as statistical analysis. We found no significant association between proviral load and phase of infection. The sequences showed high similarity without any association between BLV surface glycoprotein and the different infection stages, nor differences in the proviral load. BLV genotype 1 was identified as prevalent in the studied samples, and for the first time in Mexico, we identified BLV genotype 3 in cattle.

Bovine coronavirus in Uruguay: genetic diversity, risk factors and transboundary introductions from neighboring countries

Bovine coronavirus (BCoV) is a recognized cause of severe neonatal calf diarrhea, with a negative impact on animal welfare, leading to economic losses to the livestock industry. Cattle production is one of the most important economic sectors in Uruguay. The aim of this study was to determine the frequency of BCoV infections and their genetic diversity in Uruguayan calves and to describe the evolutionary history of the virus in South America. The overall detection rate of BCoV in Uruguay was 7.8% (64/824): 7.7% (60/782) in dairy cattle and 9.5% (4/42) in beef cattle. The detection rate of BCoV in samples from deceased and live calves was 10.0% (6/60) and 7.6% (58/763), respectively. Interestingly, there was a lower frequency of BCoV detection in calves born to vaccinated dams (3.3%, 8/240) than in calves born to unvaccinated dams (12.2%, 32/263) (OR: 4.02, 95%CI: 1.81–8.90; p = 0.00026). The frequency of BCoV detection was higher in colder months (11.8%, 44/373) than in warmer months (1.5%, 3/206) (OR: 9.05, 95%CI: 2.77–29.53, p = 0.000013). Uruguayan strains grouped together in two different lineages: one with Argentinean strains and the other with Brazilian strains. Both BCoV lineages were estimated to have entered Uruguay in 2013: one of them from Brazil (95%HPD interval: 2011–2014) and the other from Argentina (95%HPD interval: 2010–2014). The lineages differed by four amino acid changes, and both were divergent from the Mebus reference strain. Surveillance should be maintained to detect possible emerging strains that can clearly diverge at the antigenic level from vaccine strains.

Effects of the Q80K Polymorphism on the Physicochemical Properties of Hepatitis C Virus Subtype 1a NS3 Protease

Hepatitis C virus genotype 1a (HCV-1a) comprises clades I and II. The Q80K polymorphism is found predominantly in clade I but rarely in clade II. Here, we investigated whether natural polymorphisms in HCV-1a clade II entailed structural protein changes when occurrence of the Q80K variant was simulated. Based on HCV-1a clade I and II protein sequences, the structure of the HCV-1a Q80K mutant NS3-4A was obtained by comparative modeling. Its physicochemical properties were studied by molecular dynamics simulations and network analysis. Results demonstrate that, in the presence of the K80 variant, clade II protease polymorphisms A91 and S/G174 led to variations in hydrogen bond occupancies. Structural analyses revealed differences in (i) flexibility of the H57 catalytic residue on the NS3 protease and (ii) correlations between amino acids on the NS3 protease and the NS4A cofactor. The latter indicated possible destabilization of interactions, resulting in increased separation of these proteins. The present findings describe how the relationships between different HCV-1a NS3 protease amino acid residues could affect the appearance of viral variants and the existence of distinct genetic barriers to HCV-1a isolates.

Comparisons of Human Immunodeficiency Virus Type 1 Envelope Variants in Blood and Genital Fluids near the Time of Male-to-Female Transmission

To better understand the transmission of human immunodeficiency virus type 1 (HIV-1), the genetic characteristics of blood and genital viruses from males were compared to those of the imputed founding virus population in their female partners. Initially serodiscordant heterosexual African couples with sequence-confirmed male-to-female HIV-1 transmission and blood and genital specimens collected near the time of transmission were studied. Single viral templates from blood plasma and genital tract RNA and DNA were sequenced across HIV-1 env gp160. Eight of 29 couples examined yielded viral sequences from both tissues. Analysis of these couples' sequences demonstrated, with one exception, that the women's founding viral populations arose from a single viral variant and were CCR5 tropic, even though CXCR4 variants were detected within four males. The median genetic distance of the imputed most recent common ancestor of the women's founder viruses showed that they were closer to the semen viruses than to the blood viruses of their transmitting male partner, but this finding was biased by detection of a greater number of viral clades in the blood. Using multiple assays, the blood and genital viruses were consistently found to be compartmentalized in only two of eight men. No distinct amino acid signatures in the men's viruses were found to link to the women's founders, nor did the women's env sequences have shorter variable loops or fewer N-linked glycosylation sites. The lack of selective factors, except for coreceptor tropism, is consistent with others' findings in male-to-female and high-risk transmissions. The infrequent compartmentalization between the transmitters' blood and semen viruses suggests that cell-free blood virus likely includes HIV-1 sequences representative of those of viruses in semen.IMPORTANCE Mucosal transmissions account for the majority of HIV-1 infections. Identification of the viral characteristics associated with transmission would facilitate vaccine design. This study of HIV strains from transmitting males and their seroconverting female partners found that the males' genital tract viruses were rarely distinct from the blood variants. The imputed founder viruses in women were genetically similar to both the blood and genital tract variants of their male partners, indicating a lack of evidence for genital tract-specific lineages. These findings suggest that targeting vaccine responses to variants found in blood are likely to also protect from genital tract variants.

Genetic Analysis of the Full-Length gag Gene from the Earliest Korean Subclade B of HIV-1: An Outbreak among Korean Hemophiliacs

We determined the earliest full-length HIV-1 gag gene sequences in 110 patients with HIV-1, including 20 hemophiliacs (HPs) and 90 local controls (LCs). The gag gene from stored sera was amplified using RT-PCR, and was subjected to direct sequencing. Phylogenetic analysis indicated that 94 and 16 sequences belonged to the Korean subclade of HIV-1 subtype B (KSB) and subtype B, respectively. A total of 12 signature pattern amino acids were found within the KSB, distinct from the worldwide consensus of subtype B. Within the KSB, the gag gene sequences from donors O and P and those from the 20 HPs comprised two subclusters. In particular, sequences from donor O strongly clustered with those of eight HPs. Moreover, signature pattern analysis indicated that 14 signature nucleotides were shared between the HPs and LCs within KSB (p < 0.01). Among the 14 nucleotides, positions 9 and 5 belonged to clusters O and P, respectively. In conclusion, signature pattern analysis for the gag gene revealed 12 signature pattern residues within the KSB and also confirmed the previous conclusion that the 20 HPs were infected with viruses due to incompletely inactivated clotting factor IX. This study is the first genetic analysis of the HIV-1 gag gene in Korea.

Molecular epidemiology of hepatitis B virus mutants associated with vaccine escape, drug resistance and diagnosis failure

The massive implementation of the vaccine and antiviral agents against hepatitis B virus (HBV), targeting the envelope and viral polymerase genes, induces a selection pressure that might lead to the emergence of variants that impair the effectiveness of the vaccine, diagnostic methods and antiviral therapy. The aim of this study was to evaluate the prevalence of HBV vaccine escape mutants (VEMs), diagnostic failure mutants (DFMs) and treatment resistance mutants (ARMs) among individuals from Buenos Aires, Argentina. HBV surface antigen and polymerase sequences obtained from serum samples of 530 HBV-infected individuals were analysed. Samples belonged to genotypes A (28.1%), D (13.6%) and F (58.3%). VEMs, DMFs and ARMs were present in 40 (7.5%), 57 (10.7%) and 27 (5.1%) samples within the studied population. Additionally, eight nonpreviously reported VEMs and nine DFMs were identified. VEMs and DFMs were biased by genotype, being higher in genotype D (33.3% and 33.3%) compared to genotype A (6% and 17.4%) and genotype F (2.3% and 2.3%) (P > 0.001). On the contrary, there was no association between the presence of ARMs and HBV genotype (P = 0.324). VEMs, DFMs and ARMs create public health concerns. The current study provided valuable information about mutants in surface antigen and polymerase in HBV-infected patients from Argentina where HBV-F is the most prevalent genotype. Consequently, it constitutes an important reference for Latin American clinicians in order to optimize the management of HBV-infected patients.

Partial pol Sequences from Drug Naive HIV-2 Infected Individuals from Maharashtra, India

HIV-2 is important due to its unique challenges in diagnosis, treatment, and drug resistance. The data on Indian HIV-2 pol gene as well as resistance to antiretroviral drugs are limited. Here we report sequence data of protease (PR) and reverse transcriptase (RT) genes from HIV-2 infected treatment naive individuals (N = 32) from Maharashtra, India. These sequences were found to be closely related to HIV-2 subtype A sequences from Guinea Bissau. We observed two unique residues at positions 14 and 70 in the PR region specific to Indian HIV-2. Mutations associated with resistance to RT and protease inhibitors were observed in 3 of 32 (9.37%) samples. To our knowledge, this is the first study from India to report drug resistance among treatment naive HIV-2 infected individuals. The results emphasize need for larger nationwide surveillance for HIV-2 drug resistance to better understand the primary drug resistance among HIV-2 infected individuals.

Hepatitis C virus intrinsic molecular determinants may contribute to the development of cholestatic hepatitis after liver transplantation

Cholestatic hepatitis C (CHC) is a severe form of hepatitis C virus (HCV) infection recurrence that leads to high graft loss rates early after liver transplantation (LT). To investigate the pathogenic mechanisms of CHC, we analysed HCV quasispecies in CHC patients compared to a control group (mild hepatitis C recurrence) by deep pyrosequencing. At the time of LT, NS5B quasispecies complexity was similar between the two groups but, after LT, it decreased more sharply in CHC patients than in the control group. Interestingly, the major variant before LT propagated efficiently and remained as the dominant sequence after LT in 62 % of CHC patients versus 11 % of controls (P=0.031). Sequence analysis of the complete non-structural region in a limited number of patients revealed a potential 12 aa signature specific to the CHC group. These data suggest that intrinsic molecular determinants in the circulating HCV quasispecies may provide a fitness advantage, contributing to the development of CHC.

Molecular characterization of hepatitis B virus in blood donors in Botswana

Hepatitis B virus (HBV) poses a significant threat to blood transfusion safety in sub-Saharan Africa (SSA) where allogeneic blood donations are screened serologically, and more sensitive nucleic acid tests (NATs) are utilized infrequently. HBV strains circulating among blood donors in Botswana are not yet characterized. We designed a cross-sectional study to determine the HBV sub-genotypes and prevalence of hepatitis B surface antigen (HBsAg) among blood donors between November 2014 and October 2015. A total of 12,575 blood donations were screened for HBsAg and 50 consecutive plasma samples were selected for genotyping from confirmed HBsAg⁺ donations. Overlapping Pol and complete S (Pol/S) open reading frames (ORFs) were sequenced from extracted HBV DNA. To identify any signature amino acids, mutations were compared to sequences from a cohort of chronic HBV patients co-infected with HIV and were treatment naïve. The prevalence of HBsAg⁺ blood donors was 1.02% (95% CI 0.9-1.2%), and the circulating sub-genotypes were A1 serotype adw2 (36.1%), D2 serotype ayw2 (2.9%), and D3 serotypes ayw 1/2 (58.3%). Prevalence of escape mutations was 14% from HBV isolates of blood donors and 15% from isolates of HBV/HIV co-infected patients (p = 0.6926). The escape mutations sP120L, sG130R, sY134H, and sD144A were identified predominantly among HBV isolates from blood donors. These escape mutations have been associated with accelerated HBV sequelae [e.g., liver cirrhosis (LC) and hepatocellular carcinoma (HCC)], failure to detect HBsAg, inability to respond to immunoglobulin (Ig) therapy, and HBV vaccine escape. Characterizing the HBV burden, circulating sub-genotypes, and clinically relevant mutations among blood donors in Botswana is important to elucidate the efficacy of currently available vaccines, predicting HBV-transmission patterns, understanding the cohort's risk to HBV-related complications, and to developing prevention strategies and effective genotype-based antiretroviral therapies.

Genome-scale analysis of evolutionary rate and selection in a fast-expanding Spanish cluster of HIV-1 subtype F1

This work is aimed at assessing the presence of positive selection and/or shifts of the evolutionary rate in a fast-expanding HIV-1 subtype F1 transmission cluster affecting men who have sex with men in Spain. We applied Bayesian coalescent phylogenetics and selection analyses to 23 full-coding region sequences from patients belonging to that cluster, along with other 19 F1 epidemiologically-unrelated sequences. A shift in the overall evolutionary rate of the virus, explained by positively selected sites in the cluster, was detected. We also found one substitution in Nef (H89F) that was specific to the cluster and experienced positive selection. These results suggest that fast transmission could have been facilitated by some inherent genetic properties of this HIV-1 variant.

HIV-1 non-group M phenotypic susceptibility to integrase strand transfer inhibitors

Objectives

To determine natural phenotypic susceptibility of non-group M HIV-1 to integrase strand transfer inhibitors (INSTIs) in a large panel of 39 clinical strains from groups O, N and P and to identify genotypic polymorphisms according to susceptibility levels.

Methods

Susceptibility to raltegravir, elvitegravir and dolutegravir was evaluated in 36 HIV-1/O, 2 HIV-1/N and 1 HIV-1/P strains plus an HIV-1/M reference strain. IC50 values were determined after 3 days, and fold changes (FCs) were calculated relative to the HIV-1/M strain. Genotypic polymorphism was determined by amplification of codons 19-263 of the integrase; the natural occurrence of resistance-associated mutations was analysed using the main resistance algorithms and the IAS-USA list. VESPA analysis of the strain sequences was used to determine a signature pattern associated with higher FC.

Results

Similar IC50 results were observed for the three drugs. Based on the value for the HIV-1/M reference strain, the data showed FC values <2.5 for raltegravir and dolutegravir, whereas the distribution for elvitegravir was heterogeneous, with FC > 10 for six strains (15%). Analysis of the non-M integrase sequences showed a high level of polymorphism without a major genotypic impact; it also revealed mutations that may be associated with the highest FC values obtained for elvitegravir.

Conclusions

Our phenotypic data showed that non-M strains are globally susceptible to the three currently used INSTIs, but the impact of the high FC values observed for some strains with elvitegravir needs to be explored. Clinical data are now needed to confirm these phenotypic results.

The evolution of subtype B HIV-1 tat in the Netherlands during 1985-2012

For the production of viral genomic RNA, HIV-1 is dependent on an early viral protein, Tat, which is required for high-level transcription. The quantity of viral RNA detectable in blood of HIV-1 infected individuals varies dramatically, and a factor involved could be the efficiency of Tat protein variants to stimulate RNA transcription. HIV-1 virulence, measured by set-point viral load, has been observed to increase over time in the Netherlands and elsewhere. Investigation of tat gene evolution in clinical isolates could discover a role of Tat in this changing virulence. A dataset of 291 Dutch HIV-1 subtype B tat genes, derived from full-length HIV-1 genome sequences from samples obtained between 1985-2012, was used to analyse the evolution of Tat. Twenty-two patient-derived tat genes, and the control Tat_HXB2 were analysed for their capacity to stimulate expression of an LTR-luciferase reporter gene construct in diverse cell lines, as well as for their ability to complement a tat-defective HIV-1_LAI clone. Analysis of 291 historical tat sequences from the Netherlands showed ample amino acid (aa) variation between isolates, although no specific mutations were selected for over time. Of note, however, the encoded protein varied its length over the years through the loss or gain of stop codons in the second exon. In transmission clusters, a selection against the shorter Tat86 ORF was apparent in favour of the more common Tat101 version, likely due to negative selection against Tat86 itself, although random drift, transmission bottlenecks, or linkage to other variants could also explain the observation. There was no correlation between Tat length and set-point viral load; however, the number of non-intermediate variants in our study was small. In addition, variation in the length of Tat did not significantly change its capacity to stimulate transcription. From 1985 till 2012, variation in the length of the HIV-1 subtype B tat gene is increasingly found in the Dutch epidemic. However, as Tat proteins did not differ significantly in their capacity to stimulate transcription elongation in vitro, the increased HIV-1 virulence seen in recent years could not be linked to an evolving viral Tat protein.

New insights into Trypanosoma cruzi evolution, genotyping and molecular diagnostics from satellite DNA sequence analysis

Trypanosoma cruzi has been subdivided into seven Discrete Typing Units (DTUs), TcI-TcVI and Tcbat. Two major evolutionary models have been proposed to explain the origin of hybrid lineages, but while it is widely accepted that TcV and TcVI are the result of genetic exchange between TcII and TcIII strains, the origin of TcIII and TcIV is still a matter of debate. T. cruzi satellite DNA (SatDNA), comprised of 195 bp units organized in tandem repeats, from both TcV and TcVI stocks were found to have SatDNA copies type TcI and TcII; whereas contradictory results were observed for TcIII stocks and no TcIV sequence has been analyzed yet. Herein, we have gone deeper into this matter analyzing 335 distinct SatDNA sequences from 19 T. cruzi stocks representative of DTUs TcI-TcVI for phylogenetic inference. Bayesian phylogenetic tree showed that all sequences were grouped in three major clusters, which corresponded to sequences from DTUs TcI/III, TcII and TcIV; whereas TcV and TcVI stocks had two sets of sequences distributed into TcI/III and TcII clusters. As expected, the lowest genetic distances were found between TcI and TcIII, and between TcV and TcVI sequences; whereas the highest ones were observed between TcII and TcI/III, and among TcIV sequences and those from the remaining DTUs. In addition, signature patterns associated to specific T. cruzi lineages were identified and new primers that improved SatDNA-based qPCR sensitivity were designed. Our findings support the theory that TcIII is not the result of a hybridization event between TcI and TcII, and that TcIV had an independent origin from the other DTUs, contributing to clarifying the evolutionary history of T. cruzi lineages. Moreover, this work opens the possibility of typing samples from Chagas disease patients with low parasitic loads and improving molecular diagnostic methods of T. cruzi infection based on SatDNA sequence amplification.

Probing the compartmentalization of HIV-1 in the central nervous system through its neutralization properties

Compartmentalization of HIV-1 has been observed in the cerebrospinal fluid (CSF) of patients at different clinical stages. Considering the low permeability of the blood-brain barrier, we wondered if a reduced selective pressure by neutralizing antibodies (NAb) in the central nervous system (CNS) could favor the evolution of NAb-sensitive viruses in this compartment. Single genome amplification (SGA) was used to sequence full-length HIV-1 envelope variants (453 sequences) from paired CSF and blood plasma samples in 9 subjects infected by HIV variants of various clades and suffering from diverse neurologic disorders. Dynamics of viral evolution were evaluated with a bayesian coalescent approach for individuals with longitudinal samples. Pseudotyped viruses expressing envelope glycoproteins variants representative of the quasi-species present in each compartment were generated, and their sensitivity to autologous neutralization, broadly neutralizing antibodies (bNAbs) and entry inhibitors was assessed. Significant compartmentalization of HIV populations between blood and CSF were detected in 5 out of 9 subjects. Some of the previously described genetic determinants for compartmentalization in the CNS were observed regardless of the HIV-1 clade. There was no difference of sensitivity to autologous neutralization between blood- and CSF-variants, even for subjects with compartmentalization, suggesting that selective pressure by autologous NAb is not the main driver of HIV evolution in the CNS. However, we observed major differences of sensitivity to sCD4 or to at least one bNAb targeting either the N160-V1V2 site, the N332-V3 site or the CD4bs, between blood- and CSF-variants in all cases. In particular, HIV-1 variants present in the CSF were more resistant to bNAbs than their blood counterpart in some cases. Considering the possible migration from CSF to blood, the CNS could be a reservoir of bNAb resistant viruses, an observation that should be considered for immunotherapeutic approaches.

Signature pattern analysis for the full-length env gene of the earliest Korean subclade B of HIV-1: outbreak among Korean hemophiliacs

The epidemiological link in the hypervariable env gene between viruses infecting HIV-positive hemophiliacs (HPs) and plasma donors was not studied. We determined full-length env gene sequences in 20 HPs, 3 plasma donors whose plasma was used for domestic clotting factor (DCF) production, and 54 local controls (LCs). Env genes from viruses in frozen stored sera obtained 1-3 years after diagnosis and from samples collected several years after infection were amplified via RT-PCR and subjected to direct sequencing. Phylogenetic analysis indicated that all sequences were subtype B, including 133 sequences from 77 cases (20 HPs, 3 plasma donors, and 54 LCs) belonging to the Korean subclade B (KSB) and 6 sequences from 5 cases that did not belong to the KSB. Env gene sequences from donors O and P and those of the 20 HPs comprised 2 subclusters within the KSB, although phylogenetic analysis did not support significant bootstrap values. In contrast, signature pattern analysis indicated signature nucleotides at 43 positions between the HPs and LCs (P < 0.05). In particular, specific signature nucleotides at 4 positions were fully conserved in the HPs, but not in the LCs (P < 0.0001). Furthermore, there were 26 signature residues within the KSB and were distinct from the worldwide consensus for subtype B. In conclusion, signature pattern analysis for the hypervariable env gene revealed an epidemiological link that the 20 HPs in this study had been infected with viruses from the DCF used for treatment, consistent with our previous finding.

Feline Immunodeficiency Virus Evolutionarily Acquires Two Proteins, Vif and Protease, Capable of Antagonizing Feline APOBEC3

The interplay between viral and host proteins has been well studied to elucidate virus-host interactions and their relevance to virulence. Mammalian genes encode apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) proteins, which act as intrinsic restriction factors against lentiviruses. To overcome APOBEC3-mediated antiviral actions, lentiviruses have evolutionarily acquired an accessory protein, viral infectivity factor (Vif), and Vif degrades host APOBEC3 proteins via a ubiquitin/proteasome-dependent pathway. Although the Vif-APOBEC3 interaction and its evolutionary significance, particularly those of primate lentiviruses (including HIV) and primates (including humans), have been well investigated, those of nonprimate lentiviruses and nonprimates are poorly understood. Moreover, the factors that determine lentiviral pathogenicity remain unclear. Here, we focus on feline immunodeficiency virus (FIV), a pathogenic lentivirus in domestic cats, and the interaction between FIV Vif and feline APOBEC3 in terms of viral virulence and evolution. We reveal the significantly reduced diversity of FIV subtype B compared to that of other subtypes, which may associate with the low pathogenicity of this subtype. We also demonstrate that FIV subtype B Vif is less active with regard to feline APOBEC3 degradation. More intriguingly, we further reveal that FIV protease cleaves feline APOBEC3 in released virions. Taken together, our findings provide evidence that a lentivirus encodes two types of anti-APOBEC3 factors, Vif and viral protease.IMPORTANCE During the history of mammalian evolution, mammals coevolved with retroviruses, including lentiviruses. All pathogenic lentiviruses, excluding equine infectious anemia virus, have acquired the vif gene via evolution to combat APOBEC3 proteins, which are intrinsic restriction factors against exogenous lentiviruses. Here we demonstrate that FIV, a pathogenic lentivirus in domestic cats, antagonizes feline APOBEC3 proteins by both Vif and a viral protease. Furthermore, the Vif proteins of an FIV subtype (subtype B) have attenuated their anti-APOBEC3 activity through evolution. Our findings can be a clue to elucidate the complicated evolutionary processes by which lentiviruses adapt to mammals.

Silent dissemination of HTLV-1 in an endemic area of Argentina. Epidemiological and molecular evidence of intrafamilial transmission

Background

Molecular and epidemiological studies of transmission routes and risk factors for infection by HTLV-1 are extremely important in order to implement control measures, especially because of the high prevalence of HTLV-1 in several regions of the world. San Salvador de Jujuy, Northwest Argentina, is a highly endemic area for HTLV-1 and foci of tropical spastic paraparesis/HTLV-1-associated myelopathy.

Objective

To gain further insight into the role of intrafamilial transmission of HTLV-1 in a highly endemic region in Argentina.

Method

Cross-sectional study in Northwest Argentina. Epidemiological data and blood samples were collected from 28 HTLV-1 infected subjects (index cases) and 92 close relatives/cohabitants. HTLV-1 infection was diagnosed by detection of antibodies and proviral DNA. The LTR region was sequenced and analyzed for genetic distances (VESPA software), in addition to determination and identification of polymorphisms to define HTLV-1 family signatures.

Results

Fifty seven of the 120 subjects enrolled had antibodies against HTLV-1 and were typified as HTLV-1 by PCR. The prevalence rate of HTLV-1 infection in family members of infected index cases was 31.52% (29/92). The infection was significantly associated with gender, age and prolonged lactation. Identity of LTR sequences and presence of polymorphisms revealed high prevalence of mother-to-child and interspousal transmission of HTLV-1 among these families.

Conclusion

There is an ongoing and silent transmission of HTLV-1 through vertical and sexual routes within family clusters in Northwest Argentina. This evidence highlights that HTLV-1 infection should be considered as a matter of public health in Argentina, in order to introduce preventive measures as prenatal screening and breastfeeding control.

Transmitted/Founder HIV-1 Subtype C Viruses Show Distinctive Signature Patterns in Vif, Vpr, and Vpu That Are Under Subsequent Immune Pressure During Early Infection

Viral variants that predominate during early infection may exhibit constrained diversity compared with those found during chronic infection and could contain amino acid signature patterns that may enhance transmission, establish productive infection, and influence early events that modulate the infection course. We compared amino acid distributions in 17 patients recently infected with HIV-1C with patients with chronic infection. We found significantly lower entropy in inferred transmitted/founder (t/f) compared with chronic viruses and identified signature patterns in Vif and Vpr from inferred t/f viruses. We investigated sequence evolution longitudinally up to 500 days postseroconversion and compared the impact of selected substitutions on predicted human leukocyte antigen (HLA) binding affinities of published and predicted cytotoxic T-lymphocyte epitopes. Polymorphisms in Vif and Vpr during early infection occurred more frequently at epitope-HLA anchor residues and significantly decreased predicted epitope-HLA binding. Transmission-associated sequence signatures may have implications for novel strategies to prevent HIV-1 transmission.

Naturally occurring NS3 resistance-associated variants in hepatitis C virus genotype 1: Their relevance for developing countries

Hepatitis C virus (HCV) is a major cause of global morbidity and mortality, with an estimated 130-150 million infected individuals worldwide. HCV is a leading cause of chronic liver diseases including cirrhosis and hepatocellular carcinoma. Current treatment options in developing countries involve pegylated interferon-α and ribavirin as dual therapy or in combination with one or more direct-acting antiviral agents (DAA). The emergence of resistance-associated variants (RAVs) after treatment reveals the great variability of this virus leading to a great difficulty in developing effective antiviral strategies. Baseline RAVs detected in DAA treatment-naïve HCV-infected patients could be of great importance for clinical management and outcome prediction. Although the frequency of naturally occurring HCV NS3 protease inhibitor mutations has been addressed in many countries, there are only a few reports on their prevalence in South America. In this study, we investigated the presence of RAVs in the HCV NS3 serine protease region by analysing a cohort of Uruguayan patients with chronic hepatitis C who had not been treated with any DAAs and compare them with the results found for other South American countries. The results of these studies revealed that naturally occurring mutations conferring resistance to NS3 inhibitors exist in a substantial proportion of Uruguayan treatment-naïve patients infected with HCV genotype 1 enrolled in these studies. The identification of these baseline RAVs could be of great importance for patients' management and outcome prediction in developing countries.

Los Alamos Hepatitis C Virus Sequence and Human Immunology Databases: An Expanding Resource for Antiviral Research

The hepatitis C virus (HCV) resource at Los Alamos (hcv.lanl.gov) provides access to multiple databases: one containing annotated sequences and the other a repository of immunogenic epitopes. They are derived from databases originally developed for HIV research (hiv.lanl.gov). HCV and HIV are RNA viruses with relatively compact genomes (around 10 kb) that are extraordinarily variable, both within and between hosts. This diversity requires methods to track and exclude variants from an individual infection or from epidemiologically related infections, and tools to analyse the variation. The HCV immunology database contains a curated inventory of immunogenic epitopes and information about their interaction with the host immune system, with associated retrieval and analysis tools. This interactive resource provides flexible retrieval tools for sequences, epitopes, clinical information, and metadata, as well as utilities for scientific data analysis, to investigators with internet access and a web browser. This paper describes the types of data and the services that these databases offer, the tools they provide, and their configuration and use. Examples of applications to clonal analysis for drug-resistance mutations are shown.

Molecular epidemiology of hepatitis B virus in Misiones, Argentina

Hepatitis B virus (HBV) infection is a major public health problem worldwide. The aims of this study were to describe the molecular epidemiology of HBV in the Province of Misiones, Argentina and estimate the phylodynamic of the main groups in a Bayesian coalescent framework. To this end, partial or complete genome sequences were obtained from 52 blood donor candidates. The phylogenetic analysis based on partial sequences of S/P region showed a predominance of genotype D (65.4%), followed by genotype F (30.8%) and genotype A as a minority (3.8%). At subgenotype level, the circulation of subgenotypes D3 (42.3%), D2 (13.5%), F1b (11.5%) and F4 (9.6%) was mainly identified. The Bayesian coalescent analysis of 29 complete genome sequences for the main groups revealed that the subgenotypes D2 and D3 had several introductions to the region, with ancestors dating back from 1921 to 1969 and diversification events until the late '70s. The genotype F in Misiones has a more recent history; subgenotype F4 isolates were intermixed with sequences from Argentina and neighboring countries and only one significant cluster dated back in 1994 was observed. Subgenotype F1b isolates exhibited low genetic distance and formed a closely related monophyletic cluster, suggesting a very recent introduction. In conclusion, the phylogenetic and coalescent analyses showed that the European genotype D has a higher circulation, a longer history of diversification and may be responsible for the largest proportion of chronic HBV infections in the Province of Misiones. Genotype F, especially subgenotype F1b, had a more recent introduction and its diversification in the last 20years might be related to its involvement in new transmission events.

The Evolving Profile of the Signature Amino Acid Residues in HIV-1 Subtype C Tat

Using several HIV-1 tat exon 1 amino acid sequences available from public databases and additional sequences derived from a southern Indian clinical cohort, we compared the profile of the signature amino acid residues (SAR) between two different time periods, 1986-2004 and 2005-2014. The analysis identified eight positions as signature residues in subtype C Tat and demonstrated a changing pattern at four of these positions between the two periods. At three locations (histidine 29, serine 57, and proline 60), there appears to be a nonuniform negative selection against the SAR. The negative selection appears to be severe, especially against histidine 29 (p < .0001) and moderate against proline 60 (p < .0001). The negative selection against serine 57 is statistically insignificant and appears to have begun recently. At position 63, the frequency of signature residue glutamic acid increased over the past decade, although the difference was not significant. Importantly, at the three locations where the negative selection is in progress, the substitute amino acids are the generic residues present in most of the other HIV-1 subtypes. Our data demonstrate that viral evolution can subject specific amino acid residues to subtle and progressive selection pressures without affecting the prevalence of other amino acid residues.

Comprehensive Screening for Naturally Occurring Hepatitis C Virus Resistance to Direct-Acting Antivirals in the NS3, NS5A, and NS5B Genes in Worldwide Isolates of Viral Genotypes 1 to 6

There is no comprehensive study available on the natural hepatitis C virus (HCV) polymorphism in sites associated with resistance including all viral genotypes which may present variable susceptibilities to particular direct-acting antivirals (DAAs). This study aimed to analyze the frequencies, genetic barriers, and evolutionary histories of naturally occurring resistance-associated variants (RAVs) in the six main HCV genotypes. A comprehensive analysis of up to 103 RAVs was performed in 2,901, 2,216, and 1,344 HCV isolates for the NS3, NS5A, and NS5B genes, respectively. We report significant intergenotypic differences in the frequencies of natural RAVs for these three HCV genes. In addition, we found a low genetic barrier for the generation of new RAVs, irrespective of the viral genotype. Furthermore, in 1,126 HCV genomes, including sequences spanning the three genes, haplotype analysis revealed a remarkably high frequency of viruses carrying more than one natural RAV to DAAs (53% of HCV-1a, 28.5% of HCV-1b, 67.1% of HCV-6, and 100% of genotype 2, 3, 4, and 5 haplotypes). With the exception of HCV-1a, the most prevalent haplotypes showed RAVs in at least two different viral genes. Finally, evolutionary analyses revealed that, while most natural RAVs appeared recently, others have been efficiently transmitted over time and cluster in well-supported clades. In summary, and despite the observed high efficacy of DAA-based regimens, we show that naturally occurring RAVs are common in all HCV genotypes and that there is an overall low genetic barrier for the selection of resistance mutations. There is a need for natural DAA resistance profiling specific for each HCV genotype.

A novel mutant 10Ala/Arg together with mutant 144Ser/Arg of hepatitis B virus X protein involved in hepatitis B virus-related hepatocarcinogenesis in HepG2 cell lines

Hepatitis B virus (HBV) infection-related hepatocellular carcinoma (HCC) represents a major health problem worldwide. HBV X (HBx) protein is the most common open reading frame that may undergo mutations, resulting in the development of HCC. This study aimed to determine specific HBx mutations that differentiate the central- and para-tumor tissues, and identify their association with HCC development. HBx gene from HCC tumor and para-tumor tissues of 47 HCC patients was amplified, sequenced and statistically analyzed. A novel combination of 2 mutations at residues 10 and 144 was identified which might play a significant role in HCC development. Expression vectors carrying HBx with the specific mutations were constructed and transfected into HepG2 and p53-null HepG2 cells. Compared to wild type (WT) and single mutation of HBx at residue 10 or 144, the 10/144 double mutations strongly up-regulated p21 expression and prolonged G1/S transition in WT- and p53-null HepG2 cells. Apoptosis was also inhibited by HBx harboring 10/44 double-mutation. Binding of 10/144 double-mutant HBx to p53 was lower than WT HBx. Conclusively, the 10/144 double mutation of HBx might play a crucial role in HCC formation.

Genetic variability of hepatitis B virus in Uruguay: D/F, A/F genotype recombinants

Hepatitis B virus (HBV) infection is a serious global health problem. Approximately 2 billion people worldwide have been infected, and approximately 350 million individuals currently suffer from HBV-induced chronic liver infection, which causes 600,000 deaths annually from chronic hepatitis, cirrhosis and hepatocellular carcinoma. HBV is classified in eight genotypes (A-H), and two more have been proposed (I-J). In this paper, complete genome sequences of nine Uruguayan HBV are reported. Five samples belong to genotype F1b and one to genotype A2. Three HBV recombinants were detected: A1/F1b, A2/F1b and D3/F1b. The following mutations were detected: a G1896A substitution, a 33-nucleotide deletion from position 2896 to 2928 in the Pre-S1 region involving Pre-S1 residues 3-13, a 33-nt deletion in the Pre-S1 region involving nt 2913-2945 and Pre-S1 residues 9-19. More F genotypes strains than expected were detected in this study, supporting the hypothesis that there are more people of indigenous origin than declared in our population. Also, one third of the samples analyzed were recombinants. This cannot be explained by the low HBV prevalence in Uruguay, but a high HBV infection rate in drug addicts and dialysis patients could act in favor of multiple-genotype HBV infections that could lead to recombination.

Global Diversity within and between Human Herpesvirus 1 and 2 Glycoproteins

Human herpes simplex viruses 1 and 2 (HSV-1 and HSV-2) are large-genome DNA viruses that establish a persistent infection in sensory neurons and commonly manifest with recurring oral or genital erosions that transmit virus. HSV encodes 12 predicted glycoproteins that serve various functions, including cellular attachment, entry, and egress. Glycoprotein G is currently the target of an antibody test to differentiate HSV-1 from HSV-2; however, this test has shown reduced capacity to differentiate HSV strains in East Africa. Until the recent availability of 26 full-length HSV-1 and 36 full-length HSV-2 sequences, minimal comparative information was available for these viruses. In this study, we use a variety of sequence analysis methods to compare all available sequence data for HSV-1 and HSV-2 glycoproteins, using viruses isolated in Europe, Asia, North America, the Republic of South Africa, and East Africa. We found numerous differences in diversity, nonsynonymous/synonymous substitution rates, and recombination rates between HSV-1 glycoproteins and their HSV-2 counterparts. Phylogenetic analysis revealed that while most global HSV-2 glycoprotein G sequences did not form clusters within or between continents, one clade (supported at 60.5%) contained 37% of the African sequences analyzed. Accordingly, sequences from this African subset contained unique amino acid signatures, not only in glycoprotein G, but also in glycoproteins I and E, which may account for the failure of sensitive antibody tests to distinguish HSV-1 from HSV-2 in some African individuals. Consensus sequences generated in the study can be used to improve diagnostic assays that differentiate HSV-1 from HSV-2 in global populations.

IMPORTANCE

Human herpes simplex viruses 1 and 2 (HSV-1 and HSV-2) are large DNA viruses associated with recurring oral or genital erosions that transmit virus. Up to 12 HSV-1 and HSV-2 glycoproteins are involved in HSV cell entry or are required for viral spread in animals, albeit some are dispensable for replication in vitro. The recent availability of comparable numbers of full-length HSV-1 and HSV-2 sequences enabled comparative analysis of gene diversity of glycoproteins within and between HSV types. Overall, we found less glycoprotein sequence diversity within HSV-2 than within the HSV-1 strains studied, while at the same time, several HSV-2 glycoproteins were evolving under less selective pressure. Because HSV glycoproteins are the focus of antibody tests to detect and differentiate between infections with the two strains and are constituents of vaccines in clinical-stage development, these findings will aid in refining the targets for diagnostic tests and vaccines.

Hepatitis B virus genotype distribution and genotype-specific BCP/preCore substitutions in acute and chronic infections in Argentina

AIM

In order to assess Hepatitis B Virus genotype (g) and subgenotype (sg) implications in the course of infection, 234 HBsAg positive patients in different infection stages were characterized (66 acute infections, 63 HBeAg positive chronic infections and 105 anti-HBe positive chronic infections).

RESULTS

Overall, sgA2 (17.9%), gD (20.9%), sgF1b (34.2%) and sgF4 (19.7%) were the most prevalent. Subgenotype F1b was overrepresented in acute and chronic HBeAg infections (56.1%), whereas gD was the most frequent (40.0%) in anti-HBe positive chronic infections. Among chronic infections, HBeAg positivity rates were 50.0, 12.5, 62.8 and 35.3% for sgA2, gD, sgF1b and sgF4, respectively (p <0.05). A bias toward BCP/preCore mutations was observed among genotypes. In anti-HBe positive chronic infections, sgF1b was more prone to have A1762T/G1764A mutation than sgA2, sgF4 and gD (75.0, 40.0, 33.3 and 31.8%, p<0.005), whereas in the pC region, gD and sgF4 were more likely to have G1896A than sgA2 and sgF1b (81.0, 72.7, 0.0 and 31.3%, p <0.001). The unexpected low frequency of the G1896A mutation in the sgF1b (despite carrying 1858T) prompted us to perform a further analysis in order to identify genotype-specific features that could justify the pattern mutations observed. A region encompassing nucleotides 1720 to 1920 showed the higher dissimilarity between sgF1b and sgF4. Genotypes and subgenotypes carrying the 1727G, 1740C and 1773T polymorphisms were prevented to mutate position 1896.

DISCUSSION

HBeAg seroconversion is a critical event in the natural history of HBV infection. Differences in the HBeAg positivity rate might be relevant since different studies have observed that delayed HBeAg seroconversion is associated with a more severe clinical course of infection, highlighting the critical role that genotypes/subgenotypes might play in the progression of HBV infection. Polymorphisms in the regions 1720 to 1920 could be involved in the molecular mechanisms underlying seroconversion of each genotype/subgenotype.

Development of forensic assay signatures for ebolaviruses

Ebolaviruses are a diverse group of RNA viruses comprising five different species, four of which cause fatal hemorrhagic fever in humans. Because of their high infectivity and lethality, ebolaviruses are considered major biothreat agents. Although detection assays exist, no forensic assays are currently available. Here, we report the development of forensic assays that differentiate ebolaviruses. We performed phylogenetic analyses and identified canonical SNPs for all species, major clades and isolates. TaqMan-MGB allelic discrimination assays based on these SNPs were designed, screened against synthetic RNA templates, and validated against ebolavirus genomic RNAs. A total of 45 assays were validated to provide 100% coverage of the species and variants with additional resolution at the isolate level. These assays enabled accurate forensic analysis on 4 "unknown" ebolaviruses. Unknowns were correctly classified to species and variant. A goal of providing resolution below the isolate level was not successful. These high-resolution forensic assays allow rapid and accurate genotyping of ebolaviruses for forensic investigations.