Search

Background: Avian influenza virus (AIV) causes both severe outbreaks and endemic disease among poultry and has caused sporadic human infections in Asia, furthermore the routes of transmission in avian species between geographic regions can be numerous and complex. Using nucleotide sequences from the internal protein coding segments of AIV, we performed a Bayesian phylogeographic study to uncover regional routes of transmission and factors predictive of the rate of viral diffusion within China. Results: We found that the Central area and Pan-Pearl River Delta were the two main sources of AIV diffusion, while the East Coast areas especially the Yangtze River delta, were the major targets of viral invasion. Next we investigated the extent to which economic, agricultural, environmental and climatic regional data was predictive of viral diffusion by fitting phylogeographic discrete trait models using generalised linear models. Conclusions: Our results highlighted that the economic-agricultural predictors, especially the poultry population density and the number of farm product markets, are the key determinants of spatial diffusion of AIV in China; high human density and freight transportation are also important predictors of high rates of viral transmission; Climate features (e.g. temperature) were correlated to the viral invasion in the destination to some degree; while little or no impacts were found from natural environment factors (such as surface water coverage). This study uncovers the risk factors and enhances our understanding of the spatial dynamics of AIV in bird populations [ABSTRACT FROM AUTHOR]

Phenotypic differences between species are driven by changes in gene expression and, by extension, by modifications in the regulation of the transcriptome. Investigation of mammalian transcriptome divergence has been restricted to analysis of bulk gene expression levels and gene-internal splicing. Using allele-specific expression analysis in inter-strain hybrids of Mus musculus, we determined the contribution of multiple cellular regulatory systems to transcriptome divergence, including: alternative promoter usage, transcription start site selection, cassette exon usage, alternative last exon usage, and alternative polyadenylation site choice. Between mouse strains, a fifth of genes have variations in isoform usage that contribute to transcriptomic changes, half of which alter encoded amino acid sequence. Virtually all divergence in isoform usage altered the post-transcriptional regulatory instructions in gene UTRs. Furthermore, most genes with isoform differences between strains contain changes originating from multiple regulatory systems. This result indicates widespread cross-talk and coordination exists among different regulatory systems. Overall, isoform usage diverges in parallel with and independently to gene expression evolution, and the cis and trans regulatory contribution to each differs significantly. [ABSTRACT FROM AUTHOR]

Human immunodeficiency virus type 1 (HIV-1) is pandemic, but its contemporary global transmission network has not been characterized. A better understanding of the properties and dynamics of this network is essential for surveillance, prevention, and eventual eradication of HIV. Here, we apply a simple and computationally efficient network-based approach to all publicly available HIV polymerase sequences in the global database, revealing a contemporary picture of the spread of HIV-1 within and between countries. This approach automatically recovered well-characterized transmission clusters and extended other clusters thought to be contained within a single country across international borders. In addition, previously undescribed transmission clusters were discovered. Together, these clusters represent all known modes of HIV transmission. The extent of international linkage revealed by our comprehensive approach demonstrates the need to consider the global diversity of HIV, even when describing local epidemics. Finally, the speed of this method allows for near-real-time surveillance of the pandemic's progression. [ABSTRACT FROM AUTHOR]

Gene expression levels are thought to diverge primarily via regulatory mutations in trans within species, and in cis between species. To test this hypothesis in mammals we used RNA-sequencing to measure gene expression divergence between C57BL/6J and CAST/EiJ mouse strains and allele-specific expression in their F1 progeny. We identified 535 genes with parent-of-origin specific expression patterns, although few of these showed full allelic silencing. This suggests that the number of imprinted genes in a typical mouse somatic tissue is relatively small. In the set of nonimprinted genes, 32% showed evidence of divergent expression between the two strains. Of these, 2% could be attributed purely to variants acting in trans, while 43% were attributable only to variants acting in cis. The genes with expression divergence driven by changes in trans showed significantly higher sequence constraint than genes where the divergence was explained by variants acting in cis. The remaining genes with divergent patterns of expression (55%) were regulated by a combination of variants acting in cis and variants acting in trans. Intriguingly, the changes in expression induced by the cis and trans variants were in opposite directions more frequently than expected by chance, implying that compensatory regulation to stabilize gene expression levels is widespread. We propose that expression levels of genes regulated by this mechanism are fine-tuned by cis variants that arise following regulatory changes in trans, suggesting that many cis variants are not the primary targets of natural selection. [ABSTRACT FROM AUTHOR]

Background. Many studies of sexual behavior have shown that individuals vary greatly in their number of sexual partners over time, but it has proved difficult to obtain parameter estimates relating to the dynamics of human immunodeficiency virus (HIV) transmission except in small-scale contact tracing studies. Recent developments in molecular phylodynamics have provided new routes to obtain these parameter estimates, and current clinical practice provides suitable data for entire infected populations. Methods. A phylodynamic analysis was performed on partial pol gene sequences obtained for routine clinical care from 14 560 individuals, representing approximately 60% of the HIV-positive men who have sex with men (MSM) under care in the United Kingdom Results. Among individuals linked to others in the data set, 29% are linked to only 1 individual, 41% are linked to 2-10 individuals, and 29% are linked to≤10 individuals. The right-skewed degree distribution can be approximated by a power law, but the data are best fitted by a Waring distribution for all time depths. For time depths of 5-7 years, the distribution parameter q lies within the range that indicates infinite variance. Conclusions. The transmission network among UK MSM is characterized by preferential association such that a randomly distributed intervention would not be expected to stop the epidemic [ABSTRACT FROM AUTHOR]

Background: Economic evaluations alongside clinical trials are becoming increasingly common. Cost data are often collected through the use of postal questionnaires; however, the accuracy of this method is uncertain. We compared postal questionnaires with hospital records for collecting data on physiotherapy service use. Methods: As part of a randomised trial of orthopaedic medicine compared with orthopaedic surgery we collected physiotherapy use data on a group of patients from retrospective postal questionnaires and from hospital records. Results: 315 patients were referred for physiotherapy. Hospital data on attendances was available for 30% (n = 96), compared with 48% (n = 150) of patients completing questionnaire data (95% Cl for difference = 10% to 24%); 19% (n = 59) had data available from both sources. The two methods produced an intraclass correlation coefficient of 0.54 (95% Cl 0.31 to 0.70). However, the two methods produced significantly different estimates of resource use with patient self report recalling a mean of 1.3 extra visits (95% Cl 0.4 to 2.2) compared with hospital records. Conclusions: Using questionnaires in this study produced data on a greater number of patients compared with examination of hospital records. However, the two data sources did differ in the quantity of physiotherapy used and this should be taken into account in any analysis. [ABSTRACT FROM AUTHOR]

Analyzes genetic sequence cluters in HIV-1 reverse transcriptase affected by subsequent virological response to antiretroviral therapy. Description of clustering sequence around a parsimony method; Details of the virological responses for each cluster; Viral load measurements.

HIV incidence in Kazakhstan increased by 73% between 2010 and 2020, with an estimated 35,000 people living with HIV (PLHIV) in 2020. The development of antiretroviral drug resistance is a major threat to effective antiretroviral therapy (ART), yet studies on the prevalence of drug resistance in Kazakhstan are sparse. In this study on the molecular epidemiology of HIV in Kazakhstan, we analyzed 968 partial HIV-1 pol sequences that were collected between 2017 and 2020 from PLHIV across all regions of Kazakhstan, covering almost 3% of PLHIV in 2020. Sequences predominantly represented subtypes A6 (57%) and CRF02_AG (41%), with 32% of sequences exhibiting high-level drug resistance. We further identified distinct drug-resistant mutations (DRMs) in the two subtypes: subtype A6 showed a propensity for DRMs A62V, G190S, K101E, and D67N, while CRF02_AG showed a propensity for K103N and V179E. Codon usage analysis revealed that different mutational pathways for the two subtypes may explain the difference in G190S and V179E frequencies. Phylogenetic analysis highlighted differences in the timing and geographic spread of both subtypes within the country, with A62V-harboring subtype A6 sequences clustering on the phylogeny, indicative of sustained transmission of the mutation. Our findings suggest an HIV epidemic characterized by high levels of drug resistance and differential DRM frequencies between subtypes. This emphasizes the importance of drug resistance monitoring within Kazakhstan, together with DRM and subtype screening at diagnosis, to tailor drug regimens and provide effective, virally suppressive ART. [ABSTRACT FROM AUTHOR]

Background: Economic evaluations alongside clinical trials are becoming increasingly common. Cost data are often collected through the use of postal questionnaires; however, the accuracy of this method is uncertain. We compared postal questionnaires with hospital records for collecting data on physiotherapy service use.Methods: As part of a randomised trial of orthopaedic medicine compared with orthopaedic surgery we collected physiotherapy use data on a group of patients from retrospective postal questionnaires and from hospital records.Results: 315 patients were referred for physiotherapy. Hospital data on attendances was available for 30% (n = 96), compared with 48% (n = 150) of patients completing questionnaire data (95% Cl for difference = 10% to 24%); 19% (n = 59) had data available from both sources. The two methods produced an intraclass correlation coefficient of 0.54 (95% Cl 0.31 to 0.70). However, the two methods produced significantly different estimates of resource use with patient self report recalling a mean of 1.3 extra visits (95% Cl 0.4 to 2.2) compared with hospital records.Conclusions: Using questionnaires in this study produced data on a greater number of patients compared with examination of hospital records. However, the two data sources did differ in the quantity of physiotherapy used and this should be taken into account in any analysis. [ABSTRACT FROM AUTHOR]

Background HIV-1 circulating recombinant forms (CRFs) represent viral recombinant lineages that play a significant role in the global epidemic. Two of them dominate the epidemic in Burkina Faso: CRF06_cpx (first described in this country) and CRF02_AG. We reconstructed the phylodynamics of both recombinant viruses in Burkina Faso and throughout West Africa. Methods We analysed CRF06_cpx and CRF02_AG sequences (protease/gp41) from early samples collected in Burkina Faso in 1986 together with other GenBank sequences (1984–2013) in 4 datasets: African CRF06_cpx (210/60); down-sampled CRF06_cpx (146/45); Burkina Faso CRF02_AG (130/39) and West/Central African CRF02_AG (691/298). For each dataset, we analysed both protease and gp41 jointly using the BEAST multilocus analysis and conducted phylogeographic analysis to reconstruct the early migration routes between countries. Results The time to the most recent common ancestor (tMRCA) of CRF06_cpx was 1979 (1973–1983) for protease and 1981 (1978–1983) for gp41. The gp41 analysis inferred the origin of CRF06_cpx (or at least its parental subtype G lineage) in the Democratic Republic of Congo but migrated to Burkina Faso soon after (1982). Both genes showed that CRF06_cpx radiated to the rest of West Africa predominantly after around 1990. These results were robust to the oversampling of Burkina Faso sequences as they were confirmed in the down-sampled dataset. The tMRCA of the Burkina Faso CRF02_AG lineage was 1979 (1977–1983) for protease and 1980 (1978–1981) for gp41. However, we reconstructed its presence in West Africa much earlier (mid-1960s), with an initial origin in Cameroon and/or Nigeria, and its phylogeographic analysis revealed much interconnection within the region with a lack of country-specific phylogenetic patterns, which prevents tracking its exact migration routes. Conclusions Burkina Faso presents a relatively young HIV epidemic, with the diversification of the current in-country CRF02_AG and CRF06_cpx lineages taking place around 1980. This country represents the main source of CRF06_cpx in West Africa. The CRF02_AG epidemic started at least a decade earlier and showed much interchange between West African countries (especially involving coastal countries) suggesting great population mobility and an extensive viral spread in the region. [ABSTRACT FROM AUTHOR]

We present 109 near full-length HIV genomes amplified from blood serum samples obtained during early 1986 from across Uganda, which to our knowledge is the earliest and largest population sample from the initial phase of the HIV epidemic in Africa. Consensus sequences were made from paired-end Illumina reads with a target-capture approach to amplify HIV material following poor success with standard approaches. In comparisons with a smaller 'intermediate' genome dataset from 1998 to 1999 and a 'modern' genome dataset from 2007 to 2016, the proportion of subtype D was significantly higher initially, dropping from 67% (73/109), to 57% (26/46) to 17% (82/465) respectively (p < 0.0001). Subtype D has previously been shown to have a faster rate of disease progression than other subtypes in East African population studies, and to have a higher propensity to use the CXCR4 co-receptor ("X4 tropism"); associated with a decrease in time to AIDS. Here we find significant differences in predicted tropism between A1 and D subtypes in all three sample periods considered, which is particularly striking the 1986 sample: 66% (53/80) of subtype D env sequences were predicted to be X4 tropic compared with none of the 24 subtype A1. We also analysed the frequency of subtype in the envelope region of inter-subtype recombinants, and found that subtype A1 is over-represented in env, suggesting recombination and selection have acted to remove subtype D env from circulation. The reduction of subtype D frequency over three decades therefore appears to be a result of selective pressure against X4 tropism and its higher virulence. Lastly, we find a subtype D specific codon deletion at position 24 of the V3 loop, which may explain the higher propensity for subtype D to utilise X4 tropism. [ABSTRACT FROM AUTHOR]

We present 109 near full-length HIV genomes amplified from blood serum samples obtained during early 1986 from across Uganda, which to our knowledge is the earliest and largest population sample from the initial phase of the HIV epidemic in Africa. Consensus sequences were made from paired-end Illumina reads with a target-capture approach to amplify HIV material following poor success with standard approaches. In comparisons with a smaller 'intermediate' genome dataset from 1998 to 1999 and a 'modern' genome dataset from 2007 to 2016, the proportion of subtype D was significantly higher initially, dropping from 67% (73/109), to 57% (26/46) to 17% (82/465) respectively (p < 0.0001). Subtype D has previously been shown to have a faster rate of disease progression than other subtypes in East African population studies, and to have a higher propensity to use the CXCR4 co-receptor ("X4 tropism"); associated with a decrease in time to AIDS. Here we find significant differences in predicted tropism between A1 and D subtypes in all three sample periods considered, which is particularly striking the 1986 sample: 66% (53/80) of subtype D env sequences were predicted to be X4 tropic compared with none of the 24 subtype A1. We also analysed the frequency of subtype in the envelope region of inter-subtype recombinants, and found that subtype A1 is over-represented in env, suggesting recombination and selection have acted to remove subtype D env from circulation. The reduction of subtype D frequency over three decades therefore appears to be a result of selective pressure against X4 tropism and its higher virulence. Lastly, we find a subtype D specific codon deletion at position 24 of the V3 loop, which may explain the higher propensity for subtype D to utilise X4 tropism. [ABSTRACT FROM AUTHOR]

The Sustainable East Africa Research in Community Health (SEARCH) trial was a universal test-and-treat (UTT) trial in rural Uganda and Kenya, aiming to lower regional HIV-1 incidence. Here, we quantify breakthrough HIV-1 transmissions occurring during the trial from population-based, dried blood spot samples. Between 2013 and 2017, we obtained 549 gag and 488 pol HIV-1 consensus sequences from 745 participants: 469 participants infected prior to trial commencement and 276 SEARCH-incident infections. Putative transmission clusters, with a 1.5% pairwise genetic distance threshold, were inferred from maximum likelihood phylogenies; clusters arising after the start of SEARCH were identified with Bayesian time-calibrated phylogenies. Our phylodynamic approach identified nine clusters arising after the SEARCH start date: eight pairs and one triplet, representing mostly opposite-gender linked (6/9), within-community transmissions (7/9). Two clusters contained individuals with non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance, both linked to intervention communities. The identification of SEARCH-incident, within-community transmissions reveals the role of unsuppressed individuals in sustaining the epidemic in both arms of a UTT trial setting. The presence of transmitted NNRTI resistance, implying treatment failure to the efavirenz-based antiretroviral therapy (ART) used during SEARCH, highlights the need to improve delivery and adherence to up-to-date ART recommendations, to halt HIV-1 transmission. [ABSTRACT FROM AUTHOR]

Highly pathogenic (HP) avian influenza virus (AIV) H7N3 outbreaks occurred 3 times in the Americas in the past 10 years and caused severe economic loss in the affected regions. In June/July 2012, new HP H7N3 outbreaks occurred at commercial farms in Jalisco, Mexico. Outbreaks continued to be identified in neighbouring states in Mexico till August 2013. To explore the origin of this outbreak, time resolved phylogenetic trees were generated from the eight segments of full-length AIV sequences in North America using BEAST. Location, subtype, avian host species and pathogenicity were modelled as discrete traits upon the trees using continuous time Markov chains. A further joint analysis among segments was performed using a hierarchical phylogenetic model (HPM) which allowed trait rates (location, subtype, host species) to be jointly inferred across different segments. The complete spatial diffusion process was visualised through virtual globe software. Our result indicated the Mexico HP H7N3 originated from the large North America low pathogenicity AIV pool through complicated reassortment events. Different segments were contributed by wild waterfowl from different N. American flyways. Five of the eight segments (HA, NA, NP, M, NS) were introduced from wild birds migrating along the central North American flyway, and PB2, PB1 and PA were introduced via the western North American flyway. These results highlight a potential role for Mexico as a hotspot of virus reassortment as it is where wild birds from different migration routes mix during the winter. [ABSTRACT FROM AUTHOR]

In recent years, genotype I (GI) of Japanese encephalitis virus (JEV) has displaced genotype III (GIII) as the dominant virus genotype throughout Asia. In this study, the largest collection of GIII and GI envelope gene-derived viral sequences assembled to date was used to reconstruct the spatiotemporal chronology of genotype displacement throughout Asia and to determine the evolutionary and epidemiological dynamics underlying this significant event. GI consists of two clades, GI-a and GI-b, with the latter being associated with displacement of GIII as the dominant JEV genotype throughout Asia in the 1990s. Phylogeographic analysis indicated that GI-a diverged in Thailand or Cambodia and has remained confined to tropical Asia, whereas GI-b diverged in Vietnam and then dispersed northwards to China, where it was subsequently dispersed to Japan, Korea, and Taiwan. Molecular adaptation was detected by more than one method at one site (residue 15), and coevolution was detected at two pairs of sites (residues 89 to 360 and 129 to 141) within the GI E gene protein alignment. Viral multiplication and temperature sensitivity analyses in avian and mosquito cells revealed that the GI-b isolate JE-91 had significantly higher infectivity titers in mosquito cells from 24 to 48 h postinfection than did the GI-a and GIII isolates. If the JE-91 isolate is indeed representative of GI-b, an increased multiplicative ability of GI-b viruses compared to that of GIII viruses early in mosquito infection may have resulted in a shortened extrinsic incubation period that led to an increased number of GI enzootic transmission cycles and the subsequent displacement of GIII. [ABSTRACT FROM AUTHOR]

The circulation of vector-borne zoonotic viruses is largely determined by the overlap in the geographical distributions of virus-competent vectors and reservoir hosts. What is less clear are the factors influencing the distribution of virus-specific lineages. Japanese encephalitis virus (JEV) is the most important etiologic agent of epidemic encephalitis worldwide, and is primarily maintained between vertebrate reservoir hosts (avian and swine) and culicine mosquitoes. There are five genotypes of JEV: GI-V. In recent years, GI has displaced GIII as the dominant JEV genotype and GV has re-emerged after almost 60 years of undetected virus circulation. JEV is found throughout most of Asia, extending from maritime Siberia in the north to Australia in the south, and as far as Pakistan to the west and Saipan to the east. Transmission of JEV in temperate zones is epidemic with the majority of cases occurring in summer months, while transmission in tropical zones is endemic and occurs year-round at lower rates. To test the hypothesis that viruses circulating in these two geographical zones are genetically distinct, we applied Bayesian phylogeographic, categorical data analysis and phylogeny-trait association test techniques to the largest JEV dataset compiled to date, representing the envelope (E) gene of 487 isolates collected from 12 countries over 75 years. We demonstrated that GIII and the recently emerged GI-b are temperate genotypes likely maintained year-round in northern latitudes, while GI-a and GII are tropical genotypes likely maintained primarily through mosquito-avian and mosquito-swine transmission cycles. This study represents a new paradigm directly linking viral molecular evolution and climate. Author Summary: Although Japanese encephalitis virus (JEV) is a major cause of death and disability throughout tropical and temperate Asia, little is known about the evolution, geographical distribution and epidemiology of the five JEV genotypes (genetically distinct groups). To address this gap in our knowledge, we performed a genetic-based geographical analysis using the largest JEV sequence dataset assembled to date, including 487 viral sequences sampled from 12 countries over 75 years. We showed that both the newly and previously dominant genotypes of JEV are associated with temperate climates and are maintained throughout the cold winter months in northern Asia, likely by hibernating mosquitoes (survive throughout the winter), vertical transmission in mosquitoes (female to offspring), cold-blooded vertebrates and/or bats. [ABSTRACT FROM AUTHOR]

The circulation of vector-borne zoonotic viruses is largely determined by the overlap in the geographical distributions of virus-competent vectors and reservoir hosts. What is less clear are the factors influencing the distribution of virus-specific lineages. Japanese encephalitis virus (JEV) is the most important etiologic agent of epidemic encephalitis worldwide, and is primarily maintained between vertebrate reservoir hosts (avian and swine) and culicine mosquitoes. There are five genotypes of JEV: GI-V. In recent years, GI has displaced GIII as the dominant JEV genotype and GV has re-emerged after almost 60 years of undetected virus circulation. JEV is found throughout most of Asia, extending from maritime Siberia in the north to Australia in the south, and as far as Pakistan to the west and Saipan to the east. Transmission of JEV in temperate zones is epidemic with the majority of cases occurring in summer months, while transmission in tropical zones is endemic and occurs year-round at lower rates. To test the hypothesis that viruses circulating in these two geographical zones are genetically distinct, we applied Bayesian phylogeographic, categorical data analysis and phylogeny-trait association test techniques to the largest JEV dataset compiled to date, representing the envelope (E) gene of 487 isolates collected from 12 countries over 75 years. We demonstrated that GIII and the recently emerged GI-b are temperate genotypes likely maintained year-round in northern latitudes, while GI-a and GII are tropical genotypes likely maintained primarily through mosquito-avian and mosquito-swine transmission cycles. This study represents a new paradigm directly linking viral molecular evolution and climate. [ABSTRACT FROM AUTHOR]