Using phylogenetics to infer HIV-1 transmission direction between known transmission pairs.

Inferring the transmission direction between linked individuals living with HIV provides unparalleled power to understand the epidemiology that determines transmission. Phylogenetic ancestral-state reconstruction approaches infer the transmission direction by identifying the individual in whom the most recent common ancestor of the virus populations originated. While these methods vary in accuracy, it is unclear why. To evaluate the performance of phylogenetic ancestral-state reconstruction to determine the transmission direction of HIV-1 infection, we inferred the transmission direction for 112 transmission pairs where transmission direction and detailed additional information were available. We then fit a statistical model to evaluate the extent to which epidemiological, sampling, genetic, and phylogenetic factors influenced the outcome of the inference. Finally, we repeated the analysis under real-life conditions with only routinely available data. We found that whether ancestral-state reconstruction correctly infers the transmission direction depends principally on the phylogeny's topology. For example, under real-life conditions, the probability of identifying the correct transmission direction increases from 32%-when a monophyletic-monophyletic or paraphyletic-polyphyletic tree topology is observed and when the tip closest to the root does not agree with the state at the root-to 93% when a paraphyletic-monophyletic topology is observed and when the tip closest to the root agrees with the root state. Our results suggest that documenting larger differences in relative intrahost diversity increases our confidence in the transmission direction inference of linked pairs for population-level studies of HIV. These findings provide a practical starting point to determine our confidence in transmission direction inference from ancestral-state reconstruction.

Emergence of Novel Norovirus GII.4 Variant.

We detected a novel GII.4 variant with an amino acid insertion at the start of epitope A in viral protein 1 of noroviruses from the United States, Gabon, South Africa, and the United Kingdom collected during 2017-2022. Early identification of GII.4 variants is crucial for assessing pandemic potential and informing vaccine development.

Using machine learning and big data to explore the drug resistance landscape in HIV.

Drug resistance mutations (DRMs) appear in HIV under treatment pressure. DRMs are commonly transmitted to naive patients. The standard approach to reveal new DRMs is to test for significant frequency differences of mutations between treated and naive patients. However, we then consider each mutation individually and cannot hope to study interactions between several mutations. Here, we aim to leverage the ever-growing quantity of high-quality sequence data and machine learning methods to study such interactions (i.e. epistasis), as well as try to find new DRMs. We trained classifiers to discriminate between Reverse Transcriptase Inhibitor (RTI)-experienced and RTI-naive samples on a large HIV-1 reverse transcriptase (RT) sequence dataset from the UK (n ≈ 55, 000), using all observed mutations as binary representation features. To assess the robustness of our findings, our classifiers were evaluated on independent data sets, both from the UK and Africa. Important representation features for each classifier were then extracted as potential DRMs. To find novel DRMs, we repeated this process by removing either features or samples associated to known DRMs. When keeping all known resistance signal, we detected sufficiently prevalent known DRMs, thus validating the approach. When removing features corresponding to known DRMs, our classifiers retained some prediction accuracy, and six new mutations significantly associated with resistance were identified. These six mutations have a low genetic barrier, are correlated to known DRMs, and are spatially close to either the RT active site or the regulatory binding pocket. When removing both known DRM features and sequences containing at least one known DRM, our classifiers lose all prediction accuracy. These results likely indicate that all mutations directly conferring resistance have been found, and that our newly discovered DRMs are accessory or compensatory mutations. Moreover, apart from the accessory nature of the relationships we found, we did not find any significant signal of further, more subtle epistasis combining several mutations which individually do not seem to confer any resistance.

Human Immunodeficiency Virus Type 1 Phylodynamics to Detect and Characterize Active Transmission Clusters in North Carolina.

BACKGROUND: Human immunodeficiency virus type 1 (HIV-1) phylodynamics can be used to monitor epidemic trends and help target prevention through identification and characterization of transmission clusters. METHODS: We analyzed HIV-1 pol sequences sampled in North Carolina from 1997 to 2014. Putative clusters were identified using maximum-likelihood trees and dated using Bayesian Markov Chain Monte Carlo inference. Active clusters were defined as clusters including internal nodes from 2009 to 2014. Effective reproductive numbers (Re) were estimated using birth-death models for large clusters that expanded ≥2-fold from 2009 to 2014. RESULTS: Of 14 921 persons, 7508 (50%) sequences were identified in 2264 clusters. Only 288 (13%) clusters were active from 2009 to 2014; 37 were large (10-36 members). Compared to smaller clusters, large clusters were increasingly populated by men and younger persons; however, nearly 60% included ≥1 women. Clusters with ≥3 members demonstrated assortative mixing by sex, age, and sample region. Of 15 large clusters with ≥2-fold expansion, nearly all had Re approximately 1 by 2014. CONCLUSIONS: Phylodynamics revealed transmission cluster expansion in this densely sampled region and allowed estimates of Re to monitor active clusters, showing the propensity for steady, onward propagation. Associations with clustering and cluster characteristics vary by cluster size. Harnessing sequence-derived epidemiologic parameters within routine surveillance could allow refined monitoring of local subepidemics.

Sustained Low-Level Transmission of Zika and Chikungunya Viruses after Emergence in the Fiji Islands.

Zika and chikungunya viruses were first detected in Fiji in 2015. Examining surveillance and phylogenetic and serologic data, we found evidence of low-level transmission of Zika and chikungunya viruses during 2013-2017, in contrast to the major outbreaks caused by closely related virus strains in other Pacific Island countries.

Ross River Virus Antibody Prevalence, Fiji Islands, 2013-2015.

A unique outbreak of Ross River virus (RRV) infection was reported in Fiji in 1979. In 2013, RRV seroprevalence among residents was 46.5% (362/778). Of the residents who were seronegative in 2013 and retested in 2015, 10.9% (21/192) had seroconverted to RRV, suggesting ongoing endemic circulation of RRV in Fiji.

Human MX2 is an interferon-induced post-entry inhibitor of HIV-1 infection.

Animal cells harbour multiple innate effector mechanisms that inhibit virus replication. For the pathogenic retrovirus human immunodeficiency virus type 1 (HIV-1), these include widely expressed restriction factors, such as APOBEC3 proteins, TRIM5-α, BST2 (refs 4, 5) and SAMHD1 (refs 6, 7), as well as additional factors that are stimulated by type 1 interferon (IFN). Here we use both ectopic expression and gene-silencing experiments to define the human dynamin-like, IFN-induced myxovirus resistance 2 (MX2, also known as MXB) protein as a potent inhibitor of HIV-1 infection and as a key effector of IFN-α-mediated resistance to HIV-1 infection. MX2 suppresses infection by all HIV-1 strains tested, has equivalent or reduced effects on divergent simian immunodeficiency viruses, and does not inhibit other retroviruses such as murine leukaemia virus. The Capsid region of the viral Gag protein dictates susceptibility to MX2, and the block to infection occurs at a late post-entry step, with both the nuclear accumulation and chromosomal integration of nascent viral complementary DNA suppressed. Finally, human MX1 (also known as MXA), a closely related protein that has long been recognized as a broadly acting inhibitor of RNA and DNA viruses, including the orthomyxovirus influenza A virus, does not affect HIV-1, whereas MX2 is ineffective against influenza virus. MX2 is therefore a cell-autonomous, anti-HIV-1 resistance factor whose purposeful mobilization may represent a new therapeutic approach for the treatment of HIV/AIDS.

Characterization update of HIV-1 M subtypes diversity and proposal for subtypes A and D sub-subtypes reclassification.

BACKGROUND: The large and constantly evolving HIV-1 pandemic has led to an increasingly complex diversity. Because of some taxonomic difficulties among the most diverse HIV-1 subtypes, and taking advantage of the large amount of sequence data generated in the recent years, we investigated novel lineage patterns among the main HIV-1 subtypes. RESULTS: All HIV full-length genomes available in public databases were analysed (n = 2017). Maximum likelihood phylogenies and pairwise genetic distance were obtained. Clustering patterns and mean distributions of genetic distances were compared within and across the current groups, subtypes and sub-subtypes of HIV-1 to detect and analyse any divergent lineages within previously defined HIV lineages. The level of genetic similarity observed between most HIV clades was deeply consistent with the current classification. However, both subtypes A and D showed evidence of further intra-subtype diversification not fully described by the nomenclature system at the time and could be divided into several distinct sub-subtypes. CONCLUSIONS: With this work, we propose an updated nomenclature of sub-types A and D better reflecting their current genetic diversity and evolutionary patterns. Allowing a more accurate nomenclature and classification system is a necessary step for easier subtyping of HIV strains and a better detection or follow-up of viral epidemiology shifts.

Single-virion sequencing of lamivudine-treated HBV populations reveal population evolution dynamics and demographic history.

BACKGROUND: Viral populations are complex, dynamic, and fast evolving. The evolution of groups of closely related viruses in a competitive environment is termed quasispecies. To fully understand the role that quasispecies play in viral evolution, characterizing the trajectories of viral genotypes in an evolving population is the key. In particular, long-range haplotype information for thousands of individual viruses is critical; yet generating this information is non-trivial. Popular deep sequencing methods generate relatively short reads that do not preserve linkage information, while third generation sequencing methods have higher error rates that make detection of low frequency mutations a bioinformatics challenge. Here we applied BAsE-Seq, an Illumina-based single-virion sequencing technology, to eight samples from four chronic hepatitis B (CHB) patients - once before antiviral treatment and once after viral rebound due to resistance. RESULTS: With single-virion sequencing, we obtained 248-8796 single-virion sequences per sample, which allowed us to find evidence for both hard and soft selective sweeps. We were able to reconstruct population demographic history that was independently verified by clinically collected data. We further verified four of the samples independently through PacBio SMRT and Illumina Pooled deep sequencing. CONCLUSIONS: Overall, we showed that single-virion sequencing yields insight into viral evolution and population dynamics in an efficient and high throughput manner. We believe that single-virion sequencing is widely applicable to the study of viral evolution in the context of drug resistance and host adaptation, allows differentiation between soft or hard selective sweeps, and may be useful in the reconstruction of intra-host viral population demographic history.

Transmission of Non-B HIV Subtypes in the United Kingdom Is Increasingly Driven by Large Non-Heterosexual Transmission Clusters.

BACKGROUND: The United Kingdom human immunodeficiency virus (HIV) epidemic was historically dominated by HIV subtype B transmission among men who have sex with men (MSM). Now 50% of diagnoses and prevalent infections are among heterosexual individuals and mainly involve non-B subtypes. Between 2002 and 2010, the prevalence of non-B diagnoses among MSM increased from 5.4% to 17%, and this study focused on the drivers of this change. METHODS: Growth between 2007 and 2009 in transmission clusters among 14 000 subtype A1, C, D, and G sequences from the United Kingdom HIV Drug Resistance Database was analysed by risk group. RESULTS: Of 1148 clusters containing at least 2 sequences in 2007, >75% were pairs and >90% were heterosexual. Most clusters (71.4%) did not grow during the study period. Growth was significantly lower for small clusters and higher for clusters of ≥7 sequences, with the highest growth observed for clusters comprising sequences from MSM and people who inject drugs (PWID). Risk group (P< .0001), cluster size (P< .0001), and subtype (P< .01) were predictive of growth in a generalized linear model. DISCUSSION: Despite the increase in non-B subtypes associated with heterosexual transmission, MSM and PWID are at risk for non-B infections. Crossover of subtype C from heterosexuals to MSM has led to the expansion of this subtype within the United Kingdom.

G2/M cell cycle arrest correlates with primate lentiviral Vpr interaction with the SLX4 complex.

UNLABELLED: The accessory gene vpr, common to all primate lentiviruses, induces potent G2/M arrest in cycling cells. A recent study showed that human immunodeficiency virus type 1 (HIV-1) viral protein R (Vpr) mediates this through activation of the SLX4/MUS81/EME1 exonuclease complex that forms part of the Fanconi anemia DNA repair pathway. To confirm these observations, we have examined the G2/M arrest phenotypes of a panel of simian immunodeficiency virus (SIV) Vpr proteins. We show that SIV Vpr proteins differ in their ability to promote cell cycle arrest in human cells. While this is dependent on the DCAF1/DDB1/CUL4 ubiquitin ligase complex, interaction with human DCAF1 does not predict G2/M arrest activity of SIV Vpr in human cells. In all cases, SIV Vpr-mediated cell cycle arrest in human cells correlated with interaction with human SLX4 (huSLX4) and could be abolished by small interfering RNA (siRNA) depletion of any member of the SLX4 complex. In contrast, all but one of the HIV/SIV Vpr proteins tested, including those that lacked activity in human cells, were competent for G2/M arrest in grivet cells. Correspondingly, here cell cycle arrest correlated with interaction with the grivet orthologues of the SLX4 complex, suggesting a level of host adaptation in these interactions. Phylogenetic analyses strongly suggest that G2/M arrest/SLX4 interactions are ancestral activities of primate lentiviral Vpr proteins and that the ability to dysregulate the Fanconi anemia DNA repair pathway is an essential function of Vpr in vivo. IMPORTANCE: The Vpr protein of HIV-1 and related viruses is essential for the virus in vivo. The ability of Vpr to block the cell cycle at mitotic entry is well known, but the importance of this function for viral replication is unclear. Recent data have shown that HIV-1 Vpr targets the Fanconi anemia DNA repair pathway by interacting with and activating an endonuclease complex, SLX4/MUS81/EME1, that processes interstrand DNA cross-links. Here we show that the ability of a panel of SIV Vpr proteins to mediate cell cycle arrest correlates with species-specific interactions with the SLX4 complex in human and primate cells. The results of these studies suggest that the SLX4 complex is a conserved target of primate lentiviral Vpr proteins and that the ability to dysregulate members of the Fanconi anemia DNA repair pathway is essential for HIV/SIV replication in vivo.

Preservation of tetherin and CD4 counter-activities in circulating Vpu alleles despite extensive sequence variation within HIV-1 infected individuals.

The HIV-1 Vpu protein is expressed from a bi-cistronic message late in the viral life cycle. It functions during viral assembly to maximise infectious virus release by targeting CD4 for proteosomal degradation and counteracting the antiviral protein tetherin (BST2/CD317). Single genome analysis of vpu repertoires throughout infection in 14 individuals infected with HIV-1 clade B revealed extensive amino acid diversity of the Vpu protein. For the most part, this variation in Vpu increases over the course of infection and is associated with predicted epitopes of the individual's MHC class I haplotype, suggesting CD8+ T cell pressure is the major driver of Vpu sequence diversity within the host. Despite this variability, the Vpu functions of targeting CD4 and counteracting both physical virus restriction and NF-κB activation by tetherin are rigorously maintained throughout HIV-1 infection. Only a minority of circulating alleles bear lesions in either of these activities at any given time, suggesting functional Vpu mutants are heavily selected against even at later stages of infection. Comparison of Vpu proteins defective for one or several functions reveals novel determinants of CD4 downregulation, counteraction of tetherin restriction, and inhibition of NF-κB signalling. These data affirm the importance of Vpu functions for in vivo persistence of HIV-1 within infected individuals, not simply for transmission, and highlight its potential as a target for antiviral therapy.

Phenotypic characterization of virological failure following lopinavir/ritonavir monotherapy using full-length Gag-protease genes.

OBJECTIVES: Major protease mutations are rarely observed following first-line failure with PIs and interpretation of genotyping results in this context may be difficult. We performed extensive phenotyping of viruses from five patients failing lopinavir/ritonavir monotherapy in the MONARK study without major PI mutations by standard genotyping. METHODS: Phenotypic susceptibility testing and viral infectivity assessments were performed using a single-cycle assay and fold changes (FC) relative to a lopinavir-susceptible reference strain were calculated. RESULTS: >10-fold reduced baseline susceptibility to lopinavir occurred in two of five patients and >5-fold in another two. Four of five patients exhibited phylogenetic evidence of a limited viral evolution between baseline and failure, with amino acid changes at drug resistance-associated positions in one: T81A emerged in Gag with M36I in the protease gene, correlating with a reduction in lopinavir susceptibility from FC 7 (95% CI 6-8.35) to FC 13 (95% CI 8.11-17.8). Reductions in darunavir susceptibility (>5 FC) occurred in three individuals. DISCUSSION: This study suggests both baseline reduced susceptibility and evolution of resistance could be contributing factors to PI failure, despite the absence of classical PI resistance mutations by standard testing methods. Use of phenotyping also reveals lower darunavir susceptibility, warranting further study as this agent is commonly used following lopinavir failure.

Automated analysis of phylogenetic clusters.

BACKGROUND: As sequence data sets used for the investigation of pathogen transmission patterns increase in size, automated tools and standardized methods for cluster analysis have become necessary. We have developed an automated Cluster Picker which identifies monophyletic clades meeting user-input criteria for bootstrap support and maximum genetic distance within large phylogenetic trees. A second tool, the Cluster Matcher, automates the process of linking genetic data to epidemiological or clinical data, and matches clusters between runs of the Cluster Picker. RESULTS: We explore the effect of different bootstrap and genetic distance thresholds on clusters identified in a data set of publicly available HIV sequences, and compare these results to those of a previously published tool for cluster identification. To demonstrate their utility, we then use the Cluster Picker and Cluster Matcher together to investigate how clusters in the data set changed over time. We find that clusters containing sequences from more than one UK location at the first time point (multiple origin) were significantly more likely to grow than those representing only a single location. CONCLUSIONS: The Cluster Picker and Cluster Matcher can rapidly process phylogenetic trees containing tens of thousands of sequences. Together these tools will facilitate comparisons of pathogen transmission dynamics between studies and countries.

Restriction of V3 region sequence divergence in the HIV-1 envelope gene during antiretroviral treatment in a cohort of recent seroconverters.

BACKGROUND: Dynamic changes in Human Immunodeficiency Virus 1 (HIV-1) sequence diversity and divergence are associated with immune control during primary infection and progression to AIDS. Consensus sequencing or single genome amplification sequencing of the HIV-1 envelope (env) gene, in particular the variable (V) regions, is used as a marker for HIV-1 genome diversity, but population diversity is only minimally, or semi-quantitatively sampled using these methods. RESULTS: Here we use second generation deep sequencing to determine inter-and intra-patient sequence heterogeneity and to quantify minor variants in a cohort of individuals either receiving or not receiving antiretroviral treatment following seroconversion; the SPARTAC trial. We show, through a cross-sectional study of sequence diversity of the env V3 in 30 antiretroviral-naive patients during primary infection that considerable population structure diversity exists, with some individuals exhibiting highly constrained plasma virus diversity. Diversity was independent of clinical markers (viral load, time from seroconversion, CD4 cell count) of infection. Serial sampling over 60 weeks of non-treated individuals that define three initially different diversity profiles showed that complex patterns of continuing HIV-1 sequence diversification and divergence could be readily detected. Evidence for minor sequence turnover, emergence of new variants and re-emergence of archived variants could be inferred from this analysis. Analysis of viral divergence over the same time period in patients who received short (12 weeks, ART12) or long course antiretroviral therapy (48 weeks, ART48) and a non-treated control group revealed that ART48 successfully suppressed viral divergence while ART12 did not have a significant effect. CONCLUSIONS: Deep sequencing is a sensitive and reliable method for investigating the diversity of the env V3 as an important component of HIV-1 genome diversity. Detailed insights into the complex early intra-patient dynamics of env V3 diversity and divergence were explored in antiretroviral-naïve recent seroconverters. Long course antiretroviral therapy, initiated soon after seroconversion and administered for 48 weeks, restricts HIV-1 divergence significantly. The effect of ART12 and ART48 on clinical markers of HIV infection and progression is currently investigated in the SPARTAC trial.

Potent autologous and heterologous neutralizing antibody responses occur in HIV-2 infection across a broad range of infection outcomes.

Few studies have explored the role of neutralizing antibody (NAb) responses in controlling HIV-2 viremia and disease progression. Using a TZM-bl neutralization assay, we assessed heterologous and autologous NAb responses from a community cohort of HIV-2-infected individuals with a broad range of disease outcomes in rural Guinea-Bissau. All subjects (n = 40) displayed exceptionally high heterologous NAb titers (50% inhibitory plasma dilution or 50% inhibitory concentration [IC(50)], 1:7,000 to 1:1,000,000) against 5 novel primary HIV-2 envelopes and HIV-2 7312A, whereas ROD A and 3 primary envelopes were relatively resistant to neutralization. Most individuals also showed high autologous NAb against contemporaneous envelopes (78% of plasma-envelope combinations in 69 envelopes from 21 subjects), with IC(50)s above 1:10,000. No association between heterologous or autologous NAb titer and greater control of HIV-2 was found. A subset of envelopes was found to be more resistant to neutralization (by plasma and HIV-2 monoclonal antibodies). These envelopes were isolated from individuals with greater intrapatient sequence diversity and were associated with changes in potential N-linked glycosylation sites but not CD4 independence or CXCR4 use. Plasma collected from up to 15 years previously was able to potently neutralize recent autologous envelopes, suggesting a lack of escape from NAb and the persistence of neutralization-sensitive variants over time, despite significant NAb pressure. We conclude that despite the presence of broad and potent NAb responses in HIV-2-infected individuals, these are not the primary forces behind the dichotomous outcomes observed but reveal a limited capacity for adaptive selection and escape from host immunity in HIV-2 infection.

Evolution of a species-specific determinant within human CRM1 that regulates the post-transcriptional phases of HIV-1 replication.

The human immunodeficiency virus type-1 (HIV-1) Rev protein regulates the nuclear export of intron-containing viral RNAs by recruiting the CRM1 nuclear export receptor. Here, we employed a combination of functional and phylogenetic analyses to identify and characterize a species-specific determinant within human CRM1 (hCRM1) that largely overcomes established defects in murine cells to the post-transcriptional stages of the HIV-1 life cycle. hCRM1 expression in murine cells promotes the cytoplasmic accumulation of intron-containing viral RNAs, resulting in a substantial stimulation of the net production of infectious HIV-1 particles. These stimulatory effects require a novel surface-exposed element within HEAT repeats 9A and 10A, discrete from the binding cleft previously shown to engage Rev's leucine-rich nuclear export signal. Moreover, we show that this element is a unique feature of higher primate CRM1 proteins, and discuss how this sequence has evolved from a non-functional, ancestral sequence.

HIV-1 capsid-cyclophilin interactions determine nuclear import pathway, integration targeting and replication efficiency.

Lentiviruses such as HIV-1 traverse nuclear pore complexes (NPC) and infect terminally differentiated non-dividing cells, but how they do this is unclear. The cytoplasmic NPC protein Nup358/RanBP2 was identified as an HIV-1 co-factor in previous studies. Here we report that HIV-1 capsid (CA) binds directly to the cyclophilin domain of Nup358/RanBP2. Fusion of the Nup358/RanBP2 cyclophilin (Cyp) domain to the tripartite motif of TRIM5 created a novel inhibitor of HIV-1 replication, consistent with an interaction in vivo. In contrast to CypA binding to HIV-1 CA, Nup358 binding is insensitive to inhibition with cyclosporine, allowing contributions from CypA and Nup358 to be distinguished. Inhibition of CypA reduced dependence on Nup358 and the nuclear basket protein Nup153, suggesting that CypA regulates the choice of the nuclear import machinery that is engaged by the virus. HIV-1 cyclophilin-binding mutants CA G89V and P90A favored integration in genomic regions with a higher density of transcription units and associated features than wild type virus. Integration preference of wild type virus in the presence of cyclosporine was similarly altered to regions of higher transcription density. In contrast, HIV-1 CA alterations in another patch on the capsid surface that render the virus less sensitive to Nup358 or TRN-SR2 depletion (CA N74D, N57A) resulted in integration in genomic regions sparse in transcription units. Both groups of CA mutants are impaired in replication in HeLa cells and human monocyte derived macrophages. Our findings link HIV-1 engagement of cyclophilins with both integration targeting and replication efficiency and provide insight into the conservation of viral cyclophilin recruitment.

Hepatitis B virus candidate subgenotype I1 varies in distribution throughout Guangxi, China and may have originated in Long An county, Guangxi.

Sequencing of the complete hepatitis B virus (HBV) genomes from Vietnam, China and Laos led to the identification of a complex recombinant, referred to initially as an aberrant genotype and later proposed to be a new genotype, I. However, epidemiological data regarding this new genotype are lacking. A cross-sectional study was carried out to investigate the epidemiology of HBV candidate genotype I in Guangxi, China using stratified, random cluster sampling. Four thousand five hundred thirteen subjects were recruited from five counties within Guangxi. Three genotypes, B, C, and I, were identified with a prevalence of 32.6% (114/350), 64% (224/350), and 3.4% (12/350), respectively. All the genotype I isolates belong to candidate subgenotype I1 and were found in Bing Yang (15.3%, 9/59) and Na Po (5.0%, 3/60) counties only. The prevalence of this subgenotype is significantly higher in males (5.1%, 10/195) than in females (1.3%, 2/155; X(2) = 3.959, P < 0.05) but does not differ significantly with age. It was found in the Han (4.5%, 9/201) and Zhuang (3.1%, 3/97) ethnic populations only. There is no significant difference from other genotypes in the prevalence of HBV serological markers. Phylogeographic analysis revealed that genotype I1 likely arose in Long An county, then spread later to Bing Yang, Na Po counties and elsewhere in southeast Asia. In conclusion, the distribution of candidate genotype I within Guangxi is not even and it is highly endemic in some counties. Its prevalence is associated with gender and ethnicity. Subgenotype I1 likely originated in Long An county.

Emergence and pandemic spread of small ruminant lentiviruses.

Small ruminant lentiviruses (SRLVs) cause chronic, persistent infections in populations of domestic sheep (Ovis aries) and goats (Capra hircus) worldwide. The vast majority of SRLV infections involve two genotypes (A and B) that spread in association with the emergence of global livestock trade. However, SRLVs have likely been present in Eurasian ruminant populations since at least the early Neolithic period. Here, we use phylogenetic and phylogeographic approaches to reconstruct the origin of pandemic SRLV strains and infer their historical pattern of global spread. We constructed an open computational resource ('Lentivirus-GLUE') via which an up-to-date database of published SRLV sequences, multiple sequence alignments (MSAs), and sequence-associated metadata can be maintained. We used data collated in Lentivirus-GLUE to perform a comprehensive phylogenetic investigation of global SRLV diversity. Phylogenies reconstructed from genome-length alignments reveal that the deep divisions in the SRLV phylogeny are consistent with an ancient split into Eastern (A-like) and Western (B-like) lineages as agricultural systems disseminated out of domestication centres during the Neolithic period. These findings are also consistent with historical and phylogeographic evidence linking the early 20th century emergence of SRLV-A to the international export of Central Asian Karakul sheep. Investigating the global diversity of SRLVs can help reveal how anthropogenic factors have impacted the ecology and evolution of livestock diseases. The open resources generated in our study can expedite these studies and can also serve more broadly to facilitate the use of genomic data in SRLV diagnostics and research.