Search

Vaccination and anti-viral therapy with nucleos(t)ide analogues (NAs) are key approaches to reducing the morbidity, mortality and transmission of hepatitis B virus (HBV) infection. However, the efficacy of these interventions may be reduced by the emergence of drug resistance-associated mutations (RAMs) and/or vaccine escape mutations (VEMs). We have assimilated data on the global prevalence and distribution of HBV RAMs/VEMs from publicly available data and explored the evolution of these mutations. We analysed sequences downloaded from the HBV Database and calculated prevalence of 41 RAMs and 38 VEMs catalogued from published studies. We generated maximum likelihood phylogenetic trees and used treeBreaker to investigate the distribution and estimated the age of selected mutations across tree branches. RAM M204I/V had the highest prevalence, occurring in 3.8% (109/2838) of all HBV sequences in our data set, and a significantly higher rate in genotype C at 5.4% (60/1102, p = 0.0007). VEMs had an overall prevalence of 1.3% (37/2837) and had the highest prevalence in genotype C and in Asia at 2.2% (24/1102; p = 0.002) and 1.6% (34/2109; p = 0.009), respectively. Phylogenetic analysis suggested that RAM/VEMs can arise independently of treatment/vaccine exposure. In conclusion, HBV RAMs/VEMs have been found globally and across genotypes, with the highest prevalence observed in genotype C. Screening for genotype and for resistance-associated mutations may help to improve stratified patient treatment. As NAs and HBV vaccines are increasingly being deployed for HBV prevention and treatment, monitoring for resistance and advocating for better treatment regimens for HBV remains essential.

Pakistan harbours a large burden of hepatitis B virus (HBV) and hepatitis C virus (HCV) infection. We utilised repeat sero‐surveys to assess progress achieved towards hepatitis elimination in Pakistan. Multilevel logistic regression evaluated the change in HBV infection (HBV surface antigen (HBsAg)‐positive) prevalence and HCV exposure (HCV antibody (HCV‐Ab)‐positive) prevalence between two sero‐surveys from 2007 and 2019 for Sindh province and associated risk factors. Adjusted odds ratios (aORs) were estimated and population‐attributable fractions (PAF) for modifiable risk factors for HCV exposure. The 2007 and 2019 surveys included 8855 and 6672 individuals. HBsAg prevalence decreased from 2.6% (95% confidence intervals (95% CI): 2.2–2.9) in 2007 to 1.1% (95% CI: 0.8–1.3) in 2019, while HCV‐Ab prevalence increased from 5.1% (95% CI: 4.6%–5.5%) to 6.2% (95% CI: 5.6%–6.8%). The age and gender‐adjusted HBsAg prevalence decreased by 80% (aOR = 0.2, 95% CI: 0.1–0.4) among children and 60% (aOR = 0.4, 95% CI: 0.3–0.6) among adults over 2007–2019, while HCV‐Ab prevalence decreased by 60% (aOR = 0.4, 95%CI:0.2–0.7) in children and increased by 40% (aOR = 1.4, 95% CI: 1.2–1.7) in adults. HCV‐Ab prevalence was lower in adults with secondary (aOR = 0.6, 95% CI: 0.5–0.8) and higher (aOR = 0.5, 95%CI:0.3–0.8) education compared to illiterates and higher among adults reporting blood transfusion (aOR = 1.7, 95% CI: 1.2–2.4), family history of hepatitis (aOR = 2.5, 95% CI: 1.9–3.3), past year medical injection (aOR = 2.1, 95% CI: 1.6–2.7), being tattooed (aOR = 1.4, 95% CI: 1.0–1.9) and shaved by traditional barber (aOR = 1.2, 95% CI: 1.0–1.5). Modifiable risk factors accounted for 45% of HCV exposure, with medical injection(s) accounting for 38% (95%CI,25.7–48.4%). Overall HCV has increased over 2007–2019 in Sindh province, while HBV prevalence has decreased. Medical injections should be an important focus of prevention activities. [ABSTRACT FROM AUTHOR]

Motivation Target enrichment strategies generate genomic data from multiple pathogens in a single process, greatly improving sensitivity over metagenomic sequencing and enabling cost-effective, high-throughput surveillance and clinical applications. However, uptake by research and clinical laboratories is constrained by an absence of computational tools that are specifically designed for the analysis of multi-pathogen enrichment sequence data. Here we present an analysis pipeline, Castanet, for use with multi-pathogen enrichment sequencing data. Castanet is designed to work with short-read data produced by existing targeted enrichment strategies, but can be readily deployed on any BAM file generated by another methodology. Also included are an optional graphical interface and installer script. Results In addition to genome reconstruction, Castanet reports method-specific metrics that enable quantification of capture efficiency, estimation of pathogen load, differentiation of low-level positives from contamination, and assessment of sequencing quality. Castanet can be used as a traditional end-to-end pipeline for consensus generation, but its strength lies in the ability to process a flexible, pre-defined set of pathogens of interest directly from multi-pathogen enrichment experiments. In our tests, Castanet consensus sequences were accurate reconstructions of reference sequences, including in instances where multiple strains of the same pathogen were present. Castanet performs effectively on standard computers and can process the entire output of a 96-sample enrichment sequencing run (50M reads) using a single batch process command, in $<$2 h. Availability and implementation Source code freely available under GPL-3 license at https://github.com/MultipathogenGenomics/castanet , implemented in Python 3.10 and supported in Ubuntu Linux 22.04. The data underlying this article are available in Europe Nucleotide Archives, at https://www.ebi.ac.uk/ena/browser/view/PRJEB77004. [ABSTRACT FROM AUTHOR]

Host–pathogen coevolution is defined as the reciprocal evolutionary changes in both species due to genotype × genotype (G×G) interactions at the genetic level determining the outcome and severity of infection. While co-analyses of hosts and pathogen genomes (co-genome-wide association studies) allow us to pinpoint the interacting genes, these do not reveal which host genotype(s) is/are resistant to which pathogen genotype(s). The knowledge of this so-called infection matrix is important for agriculture and medicine. Building on established theories of host–pathogen interactions, we here derive four novel indices capturing the characteristics of the infection matrix. These indices can be computed from full genome polymorphism data of randomly sampled uninfected hosts, as well as infected hosts and their pathogen strains. We use these indices in an approximate Bayesian computation method to pinpoint loci with relevant G×G interactions and to infer their underlying interaction matrix. In a combined single nucleotide polymorphism dataset of 451 European humans and their infecting hepatitis C virus (HCV) strains and 503 uninfected individuals, we reveal a new human candidate gene for resistance to HCV and new virus mutations matching human genes. For two groups of significant human–HCV (G×G) associations, we infer a gene-for-gene infection matrix, which is commonly assumed to be typical of plant–pathogen interactions. Our model-based inference framework bridges theoretical models of G×G interactions with host and pathogen genomic data. It, therefore, paves the way for understanding the evolution of key G×G interactions underpinning HCV adaptation to the European human population after a recent expansion. [ABSTRACT FROM AUTHOR]

Broadly neutralising antibodies (bNAbs) targeting HIV show promise for both prevention of infection and treatment. Among these, 10-1074 has shown potential in neutralising a wide range of HIV strains. However, resistant viruses may limit the clinical efficacy of 10-1074. The prevalence of both de novo and emergent 10-1074 resistance will determine its use at a population level both to protect against HIV transmission and as an option for treatment. To help understand this further, we report the prevalence of pre-existing mutations associated with 10-1074 resistance in a bNAb-naive population of 157 individuals presenting to UK HIV centres with primary HIV infection, predominantly B clade, receiving antiretroviral treatment. Single genome analysis of HIV proviral envelope sequences showed that 29% of participants' viruses tested had at least one sequence with 10-1074 resistance-associated mutations. Mutations interfering with the glycan binding site at HIV Env position 332 accounted for 95% of all observed mutations. Subsequent analysis of a larger historic dataset of 2425 Bclade envelope sequences sampled from 1983 to 2019 revealed an increase of these mutations within the population over time. Clinical studies have shown that the presence of pre-existing bNAb mutations may predict diminished therapeutic effectiveness of 10-1074. Therefore, we emphasise the importance of screening for these mutations before initiating 10-1074 therapy, and to consider the implications of pre-existing resistance when designing prevention strategies. [ABSTRACT FROM AUTHOR]

The hepatitis C virus (HCV) epidemic in Western countries is primarily perpetuated by the sub-populations of men who have sex with men (MSM) and people who inject drugs (PWID). Understanding the dynamics of transmission in these communities is crucial for removing the remaining hurdles towards HCV elimination. We sequenced 269 annotated HCV plasma samples using probe enrichment and next-generation sequencing, obtaining 224 open reading frames of HCV (OR497849-OR498072). Maximum likelihood phylogenies were generated on the four most prevalent subtypes in this study (HCV1a, 1b, 3a, 4d) with a subsequent transmission cluster analysis. The highest rate of clustering was observed for HCV4d samples (13/17 (76.47%)). The second highest rate of clustering was observed in HCV1a samples (42/78 (53.85%)) with significant association with HIV-positive MSM. HCV1b and HCV3a had very low rates of clustering (2/83 (2.41%) and (0/29)). The spread of the prevalent subtype HCV1b appears to have been largely curtailed, and we demonstrate the onwards transmission of HCV1a and HCV4d in the HIV-positive MSM population across municipal borders. More systematic data collection and sequencing is needed to allow a better understanding of the HCV transmission among the community of PWID and overcome the remaining barriers for HCV elimination in Belgium. [ABSTRACT FROM AUTHOR]

The host genetic background for hepatocellular carcinoma (HCC) is incompletely understood. We aimed to determine if four germline genetic polymorphisms, rs429358 in apolipoprotein E (APOE), rs2642438 in mitochondrial amidoxime reducing component 1 (MARC1), rs2792751 in glycerol‐3‐phosphate acyltransferase (GPAM), and rs187429064 in transmembrane 6 superfamily member 2 (TM6SF2), previously associated with progressive alcohol‐related and nonalcoholic fatty liver disease, are also associated with HCC. Four HCC case‐control data sets were constructed, including two mixed etiology data sets (UK Biobank and FinnGen); one hepatitis C virus (HCV) cohort (STOP‐HCV), and one alcohol‐related HCC cohort (Dresden HCC). The frequency of each variant was compared between HCC cases and cirrhosis controls (i.e., patients with cirrhosis without HCC). Population controls were also considered. Odds ratios (ORs) associations were calculated using logistic regression, adjusting for age, sex, and principal components of genetic ancestry. Fixed‐effect meta‐analysis was used to determine the pooled effect size across all data sets. Across four case‐control data sets, 2,070 HCC cases, 4,121 cirrhosis controls, and 525,779 population controls were included. The rs429358:C allele (APOE) was significantly less frequent in HCC cases versus cirrhosis controls (OR, 0.71; 95% confidence interval [CI], 0.61‐0.84; P = 2.9 × 10−5). Rs187429064:G (TM6SF2) was significantly more common in HCC cases versus cirrhosis controls and exhibited the strongest effect size (OR, 2.03; 95% CI, 1.45‐2.86; P = 3.1 × 10−6). In contrast, rs2792751:T (GPAM) was not associated with HCC (OR, 1.01; 95% CI, 0.90‐1.13; P = 0.89), whereas rs2642438:A (MARC1) narrowly missed statistical significance (OR, 0.91; 95% CI, 0.84‐1.00; P = 0.043). Conclusion: This study associates carriage of rs429358:C (APOE) with a reduced risk of HCC in patients with cirrhosis. Conversely, carriage of rs187429064:G in TM6SF2 is associated with an increased risk of HCC in patients with cirrhosis. [ABSTRACT FROM AUTHOR]

Studies have shown that hepatitis C virus subtype 3a (HCV-3a) is likely to have been circulating in South Asia before its global spread. However, the time and route of this dissemination remain unclear. For the first time, we generated host and virus genome-wide data for more than 500 patients infected with HCV-3a from the UK, North America, Australia, and New Zealand. We used the host genomic data to infer the ancestry of the patients and used this information to investigate the epidemic history of HCV-3a. We observed that viruses from hosts of South Asian ancestry clustered together near the root of the tree, irrespective of the sampling country, and that they were more diverse than viruses from other host ancestries. We hypothesized that South Asian hosts are more likely to have been infected in South Asia and used the inferred host ancestries to distinguish between the location where the infection was acquired and where the sample was taken. Next, we inferred that three independent transmission events resulted in the spread of the virus from South Asia to the UK, North America, and Oceania. This initial spread happened during or soon after the end of World War II. This was subsequently followed by many independent transmissions between the UK, North America, and Oceania. Using both host and virus genomic information can be highly informative in studying the virus epidemic history, especially in the context of chronic infections where migration histories need to be accounted for. [ABSTRACT FROM AUTHOR]

Natural Killer (NK) cells contribute to HIV control in adults, but HLA-B-mediated T-cell activity has a more substantial impact on disease outcome. However, the HLA-B molecules influencing immune control in adults have less impact on paediatric infection. To investigate the contribution NK cells make to immune control, we studied >300 children living with HIV followed over two decades in South Africa. In children, HLA-B alleles associated with adult protection or disease-susceptibility did not have significant effects, whereas Bw4 (p = 0.003) and low HLA-A expression (p = 0.002) alleles were strongly associated with immunological and viral control. In a comparator adult cohort, Bw4 and HLA-A expression contributions to HIV disease outcome were dwarfed by those of protective and disease-susceptible HLA-B molecules. We next investigated the immunophenotype and effector functions of NK cells in a subset of these children using flow cytometry. Slow progression and better plasma viraemic control were also associated with high frequencies of less terminally differentiated NKG2A+NKp46+CD56dim NK cells strongly responsive to cytokine stimulation and linked with the immunogenetic signature identified. Future studies are indicated to determine whether this signature associated with immune control in early life directly facilitates functional cure in children. Author summary: In adults, immune control of HIV is strongly influenced by antiviral CD8+ T-cell responses restricted by 'protective' HLA class I molecules, such as HLA-B*57, and by 'disease-susceptible' HLA class I molecules such as HLA-B*58:02. By contrast, Natural Killer (NK) cells responses make a smaller, albeit significant, contribution. In this study, we evaluate in children living with HIV the contribution of NK cell responses to immune control of HIV, in an age group where HIV-specific CD8+ T-cell responses have less impact on disease outcome. A cohort of >300 therapy-naïve children living with HIV shows that a genetic signature favouring a KIR-education on NK cells is associated with slow progression and better viraemic control. Consistent with this, we observed control of HIV viraemia and lower total HIV DNA levels among children was associated with a less differentiated NKG2A+NKp46+ CD56dim NK cell population that functionally was highly responsive to cytokine stimulation. Thus, the study identifies a signature that can impact future therapeutic strategies to achieve remission in children. [ABSTRACT FROM AUTHOR]

Sustained viral response (SVR) rates for direct‐acting antiviral (DAA) therapy for hepatitis C virus (HCV) infection routinely exceed 95%. However, a small number of patients require retreatment. Sofosbuvir, velpatasvir and voxilaprevir (SOF/VEL/VOX) is a potent DAA combination primarily used for the retreatment of patients who failed by DAA therapies. Here we evaluate retreatment outcomes and the effects of resistance‐associated substitutions (RAS) in a real‐world cohort, including a large number of genotype (GT)3 infected patients. 144 patients from the UK were retreated with SOF/VEL/VOX following virologic failure with first‐line DAA treatment regimens. Full‐length HCV genome sequencing was performed prior to retreatment with SOF/VEL/VOX. HCV subtypes were assigned and RAS relevant to each genotype were identified. GT1a and GT3a each made up 38% (GT1a n = 55, GT3a n = 54) of the cohort. 40% (n = 58) of patients had liver cirrhosis of whom 7% (n = 4) were decompensated, 10% (n = 14) had hepatocellular carcinoma (HCC) and 8% (n = 12) had received a liver transplant prior to retreatment. The overall retreatment SVR12 rate was 90% (129/144). On univariate analysis, GT3 infection (50/62; SVR = 81%, p =.009), cirrhosis (47/58; SVR = 81%, p =.01) and prior treatment with SOF/VEL (12/17; SVR = 71%, p =.02) or SOF+DCV (14/19; SVR = 74%, p =.012) were significantly associated with retreatment failure, but existence of pre‐retreatment RAS was not when viral genotype was taken into account. Retreatment with SOF/VEL/VOX is very successful for non‐GT3‐infected patients. However, for GT3‐infected patients, particularly those with cirrhosis and failed by initial SOF/VEL treatment, SVR rates were significantly lower and alternative retreatment regimens should be considered. [ABSTRACT FROM AUTHOR]

The rapid evolution of RNA viruses has been long considered to result from a combination of high copying error frequencies during RNA replication, short generation times and the consequent extensive fixation of neutral or adaptive changes over short periods. While both the identities and sites of mutations are typically modelled as being random, recent investigations of sequence diversity of SARS coronavirus 2 (SARS-CoV-2) have identified a preponderance of C->U transitions, proposed to be driven by an APOBEC-like RNA editing process. The current study investigated whether this phenomenon could be observed in datasets of other RNA viruses. Using a 5% divergence filter to infer directionality, 18 from 36 datasets of aligned coding region sequences from a diverse range of mammalian RNA viruses (including Picornaviridae, Flaviviridae, Matonaviridae, Caliciviridae and Coronaviridae) showed a >2-fold base composition normalised excess of C->U transitions compared to U->C (range 2.1x–7.5x), with a consistently observed favoured 5' U upstream context. The presence of genome scale RNA secondary structure (GORS) was the only other genomic or structural parameter significantly associated with C->U/U->C transition asymmetries by multivariable analysis (ANOVA), potentially reflecting RNA structure dependence of sites targeted for C->U mutations. Using the association index metric, C->U changes were specifically over-represented at phylogenetically uninformative sites, potentially paralleling extensive homoplasy of this transition reported in SARS-CoV-2. Although mechanisms remain to be functionally characterised, excess C->U substitutions accounted for 11–14% of standing sequence variability of structured viruses and may therefore represent a potent driver of their sequence diversification and longer-term evolution. Author summary: The rapid evolution of RNA viruses is thought to arise from high mutation frequencies during replication and the rapid accumulation of genetic changes over time in response to its changing environments. This study describes an additional potent factor that contributes to the evolution of RNA infecting mammals, the occurrence of specific C->U mutations in a subset of RNA viruses that possess intensely structured genomes. This mutational process substantially damages the virus's ability to replicate and is potentially akin to "genome editing" of HIV-1 and other retrovirus genome sequences by APOBEC, one of the principal components of vertebrate antiretroviral defence mechanisms. This study however provides evidence for a wider mutational activity against several human and veterinary RNA viruses, including HCV and foot and mouth disease virus. The C->U mutational process accounted for 15–20% of standing sequence variability of these RNA viruses, representing a potent driver of their sequence diversification and longer-term evolution. [ABSTRACT FROM AUTHOR]

Identification of protective T cell responses against SARS-CoV-2 requires distinguishing people infected with SARS-CoV-2 from those with cross-reactive immunity to other coronaviruses. Here we show a range of T cell assays that differentially capture immune function to characterise SARS-CoV-2 responses. Strong ex vivo ELISpot and proliferation responses to multiple antigens (including M, NP and ORF3) are found in 168 PCR-confirmed SARS-CoV-2 infected volunteers, but are rare in 119 uninfected volunteers. Highly exposed seronegative healthcare workers with recent COVID-19-compatible illness show T cell response patterns characteristic of infection. By contrast, >90% of convalescent or unexposed people show proliferation and cellular lactate responses to spike subunits S1/S2, indicating pre-existing cross-reactive T cell populations. The detection of T cell responses to SARS-CoV-2 is therefore critically dependent on assay and antigen selection. Memory responses to specific non-spike proteins provide a method to distinguish recent infection from pre-existing immunity in exposed populations. Understanding the immune response to SARS-CoV-2 is dependent on being able to distinguish COVID-19 immune responses from cross-reactive immune responses to other coronaviruses. Here the authors show that choice of antigens and whether an ICS, ELISPOT or T cell proliferation assay is used has a major effect on this discriminatory ability. [ABSTRACT FROM AUTHOR]

Targeted metagenomics using strand-specific libraries with target enrichment is a sensitive, generalized approach to pathogen sequencing and transcriptome profiling. Using this method, we recovered 13 (76%) complete human respiratory syncytial virus (RSV) genomes from 17 clinical respiratory samples, reconstructed the phylogeny of the infecting viruses, and detected differential gene expression between 2 RSV subgroups, specifically, a lower expression of the P gene and a higher expression of the M2 gene in RSV-A than in RSV-B. This methodology can help to relate viral genetics to clinical phenotype and facilitate ongoing population-level RSV surveillance and vaccine development. Clinical Trials Registration. NCT03627572 and NCT03756766. [ABSTRACT FROM AUTHOR]

Choice of direct acting antiviral (DAA) therapy for Hepatitis C Virus (HCV) in the United Kingdom and similar settings usually requires knowledge of the genotype and, in some cases, antiviral resistance (AVR) profile of the infecting virus. To determine these, most laboratories currently use Sanger technology, but next-generation sequencing (NGS) offers potential advantages in throughput and accuracy. However, NGS poses unique technical challenges, which require idiosyncratic development and technical validation approaches. This applies particularly to virology, where sequence diversity is high and the amount of starting genetic material is low, making it difficult to distinguish real data from artifacts. We describe the development and technical validation of a sequence capture-based HCV whole genome sequencing (WGS) assay to determine viral genotype and AVR profile. We use clinical samples of known subtypes and viral loads, and simulated FASTQ datasets to validate the analytical performances of both the wet laboratory and bioinformatic pipeline procedures. We show high concordance of the WGS assay compared to current "gold standard" Sanger assays. Specificity was 92.3 and 96.1% for AVR and genotyping, respectively. Discordances were due to the inability of Sanger assays to assign the correct subtype or accurately call mixed drug-resistant variants. We show high repeatability and reproducibility with >99.8% sequence similarity between sequence runs as well as high precision for variant frequency detection at >98.8% in the 95th percentile. Post-sequencing bioinformatics quality control workflows allow the accurate distinction between mixed infections, cross-contaminants and recombinant viruses at a threshold of >5% for the minority population. The sequence capture-based HCV WGS assay is more accurate than legacy AVR and genotyping assays. The assay has now been implemented in the clinical pathway of England's National Health Service HCV treatment programs, representing the first validated HCV WGS pipeline in clinical service. The data generated will additionally provide granular national-level genomic information for public health policy making and support the WHO HCV elimination strategy. [ABSTRACT FROM AUTHOR]

The distribution of a phenotype on a phylogenetic tree is often a quantity of interest. Many phenotypes have imperfect heritability, so that a measurement of the phenotype for an individual can be thought of as a single realization from the phenotype distribution of that individual. If all individuals in a phylogeny had the same phenotype distribution, measured phenotypes would be randomly distributed on the tree leaves. This is, however, often not the case, implying that the phenotype distribution evolves over time. Here we propose a new model based on this principle of evolving phenotype distribution on the branches of a phylogeny, which is different from ancestral state reconstruction where the phenotype itself is assumed to evolve. We develop an efficient Bayesian inference method to estimate the parameters of our model and to test the evidence for changes in the phenotype distribution. We use multiple simulated data sets to show that our algorithm has good sensitivity and specificity properties. Since our method identifies branches on the tree on which the phenotype distribution has changed, it is able to break down a tree into components for which this distribution is unique and constant. We present two applications of our method, one investigating the association between HIV genetic variation and human leukocyte antigen and the other studying host range distribution in a lineage of Salmonella enterica, and we discuss many other potential applications. [ABSTRACT FROM AUTHOR]

Next-generation sequencing has critical applications in virus discovery, diagnostics, and environmental surveillance. We used metagenomic sequence libraries for retrospective screening of plasma samples for the recently discovered human hepegivirus 1 (HHpgV-1). From a cohort of 150 hepatitis C virus (HCV)-positive case-patients, we identified 2 persons with HHpgV-1 viremia and a high frequency of human pegivirus (HPgV) viremia (14%). Detection of HHpgV-1 and HPgV was concordant with parallel PCR-based screening using conserved primers matching groups 1 (HPgV) and 2 (HHPgV-1) nonstructural 3 region sequences. PCR identified 1 HHPgV-1-positive person with viremia from a group of 195 persons with hemophilia who had been exposed to nonvirally inactivated factor VII/IX; 18 (9%) were HPgV-positive. Relative to HCV and HPgV, active infections with HHpgV-1 were infrequently detected in blood, even in groups that had substantial parenteral exposure. Our findings are consistent with lower transmissibility or higher rates of virus clearance for HHpgV-1 than for other bloodborne human flaviviruses. [ABSTRACT FROM AUTHOR]

Recent studies have suggested greater HIV cure potential among infected children than adults. A major obstacle to HIV eradication in adults is that the viral reservoir is largely comprised of HIV-specific cytotoxic T lymphocyte (CTL) escape variants. We here evaluate the potential for CTL in HIV-infected slow-progressor children to play an effective role in “shock-and-kill” cure strategies. Two distinct subgroups of children were identified on the basis of viral load. Unexpectedly, in both groups, as in adults, HIV-specific CTL drove the selection of escape variants across a range of epitopes within the first weeks of infection. However, in HIV-infected children, but not adults, de novo autologous variant-specific CTL responses were generated, enabling the pediatric immune system to “corner” the virus. Thus, even when escape variants are selected in early infection, the capacity in children to generate variant-specific anti-HIV CTL responses maintains the potential for CTL to contribute to effective shock-and-kill cure strategies in pediatric HIV infection.

Patterns of linkage disequilibrium, homoplasy, and incompatibility are difficult to interpret because they depend on several factors, including the recombination process and the population structure. Here we introduce a novel model-based framework to infer recombination properties from such summary statistics in bacterial genomes. The underlying model is sequentially Markovian so that data can be simulated very efficiently, and we use approximate Bayesian computation techniques to infer parameters. As this does not require us to calculate the likelihood function, the model can be easily extended to investigate less probed aspects of recombination. In particular, we extend our model to account for the bias in the recombination process whereby closely related bacteria recombine more often with one another. We show that this model provides a good fit to a data set of Bacillus cereus genomes and estimate several recombination properties, including the rate of bias in recombination. All the methods described in this article are implemented in a software package that is freely available for download at http://code.google.com/p/clonalorigin/. [ABSTRACT FROM AUTHOR]

This has been updated so that all affiliations are accurate in the HTML version of the Article. The original HTML version of this Article contained errors in the author affiliations, which were correct in the PDF version. These authors jointly supervised this work: Tanya Golubchik, Andrew J. Pollard. [Extracted from the article]

Chronic hepatitis B virus (HBV) infection affects 257 million people globally. Current therapies suppress HBV but viral rebound occurs on cessation of therapy; novel therapeutic strategies are urgently required. To develop a therapeutic HBV vaccine that can induce high magnitude T cells to all major HBV antigens, we have developed a novel HBV vaccine using chimpanzee adenovirus (ChAd) and modified vaccinia Ankara (MVA) viral vectors encoding multiple HBV antigens. ChAd vaccine alone generated very high magnitude HBV specific T cell responses to all HBV major antigens. The inclusion of a shark Invariant (SIi) chain genetic adjuvant significantly enhanced the magnitude of T-cells against HBV antigens. Compared to ChAd alone vaccination, ChAd-prime followed by MVA-boost vaccination further enhanced the magnitude and breadth of the vaccine induced T cell response. Intra-cellular cytokine staining study showed that HBV specific CD8+ and CD4+ T cells were polyfunctional, producing combinations of IFNγ, TNF-α, and IL-2. In summary, we have generated genetically adjuvanted ChAd and MVA vectored HBV vaccines with the potential to induce high-magnitude T cell responses through a prime-boost therapeutic vaccination approach. These pre-clinical studies pave the way for new studies of HBV therapeutic vaccination in humans with chronic hepatitis B infection. [ABSTRACT FROM AUTHOR]

Using deep sequencing technologies such as Illumina's platform, it is possible to obtain reads from the viral RNA population revealing the viral genome diversity within a single host. A range of software tools and pipelines can transform raw deep sequencing reads into Sequence Alignment Mapping (SAM) files. We propose that interpretation tools should process these SAM files, directly translating individual reads to amino acids in order to extract statistics of interest such as the proportion of different amino acid residues at specific sites. This preserves per-read linkage between nucleotide variants at different positions within a codon location. The samReporter is a subsystem of the GLUE software toolkit which follows this direct read translation approach in its processing of SAM files. We test samReporter on a deep sequencing dataset obtained from a cohort of 241 UK HCV patients for whom prior treatment with direct-acting antivirals has failed; deep sequencing and resistance testing have been suggested to be of clinical use in this context. We compared the polymorphism interpretation results of the samReporter against an approach that does not preserve per-read linkage. We found that the samReporter was able to properly interpret the sequence data at resistance-associated locations in nine patients where the alternative approach was equivocal. In three cases, the samReporter confirmed that resistance or an atypical substitution was present at NS5A position 30. In three further cases, it confirmed that the sofosbuvir-resistant NS5B substitution S282T was absent. This suggests the direct read translation approach implemented is of value for interpreting viral deep sequencing data. [ABSTRACT FROM AUTHOR]

Summary The Zika virus (ZIKV) epidemic in the Americas established ZIKV as a major public health threat and uncovered its association with severe diseases, including microcephaly. However, genetic epidemiology in some at-risk regions, particularly Central America and Mexico, remains limited. We report 61 ZIKV genomes from this region, generated using metagenomic sequencing with ZIKV-specific enrichment, and combine phylogenetic, epidemiological, and environmental data to reconstruct ZIKV transmission. These analyses revealed multiple independent ZIKV introductions to Central America and Mexico. One introduction, likely from Brazil via Honduras, led to most infections and the undetected spread of ZIKV through the region from late 2014. Multiple lines of evidence indicate biannual peaks of ZIKV transmission in the region, likely driven by varying local environmental conditions for mosquito vectors and herd immunity. The spatial and temporal heterogeneity of ZIKV transmission in Central America and Mexico challenges arbovirus surveillance and disease control measures. [ABSTRACT FROM AUTHOR]

Broadly neutralising antibodies (bNAbs) targeting HIV show promise for both prevention of infection and treatment. Among these, 10-1074 has shown potential in neutralising a wide range of HIV strains. However, resistant viruses may limit the clinical efficacy of 10-1074. The prevalence of both de novo and emergent 10-1074 resistance will determine its use at a population level both to protect against HIV transmission and as an option for treatment. To help understand this further, we report the prevalence of pre-existing mutations associated with 10-1074 resistance in a bNAb-naive population of 157 individuals presenting to UK HIV centres with primary HIV infection, predominantly B clade, receiving antiretroviral treatment. Single genome analysis of HIV proviral envelope sequences showed that 29% of participants' viruses tested had at least one sequence with 10-1074 resistance-associated mutations. Mutations interfering with the glycan binding site at HIV Env position 332 accounted for 95% of all observed mutations. Subsequent analysis of a larger historic dataset of 2425 B-clade envelope sequences sampled from 1983 to 2019 revealed an increase of these mutations within the population over time. Clinical studies have shown that the presence of pre-existing bNAb mutations may predict diminished therapeutic effectiveness of 10-1074. Therefore, we emphasise the importance of screening for these mutations before initiating 10-1074 therapy, and to consider the implications of pre-existing resistance when designing prevention strategies., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Zacharopoulou, Lee, Oliveira, Thornhill, Robinson, Brown, Kinloch, Goulder, Fox, Fidler, Ansari and Frater.)

Hepatitis B virus (HBV) is one of the smallest human DNA viruses and its 3.2 Kb genome encodes multiple overlapping open reading frames, making its viral transcriptome challenging to dissect. Previous studies have combined quantitative PCR and Next Generation Sequencing to identify viral transcripts and splice junctions, however the fragmentation and selective amplification used in short read sequencing precludes the resolution of full length RNAs. Our study coupled an oligonucleotide enrichment protocol with state-of-the-art long read sequencing (PacBio) to identify the repertoire of HBV RNAs. This methodology provides sequencing libraries where up to 25 % of reads are of viral origin and enable the identification of canonical (unspliced), non-canonical (spliced) and chimeric viral-human transcripts. Sequencing RNA isolated from de novo HBV infected cells or those transfected with 1.3 × overlength HBV genomes allowed us to assess the viral transcriptome and to annotate 5' truncations and polyadenylation profiles. The two HBV model systems showed an excellent agreement in the pattern of major viral RNAs, however differences were noted in the abundance of spliced transcripts. Viral-host chimeric transcripts were identified and more commonly found in the transfected cells. Enrichment capture and PacBio sequencing allows the assignment of canonical and non-canonical HBV RNAs using an open-source analysis pipeline that enables the accurate mapping of the HBV transcriptome.

Objective: Hepatocellular carcinoma (HCC) often develops in patients with alcohol-related cirrhosis at an annual risk of up to 2.5%. Some host genetic risk factors have been identified but do not account for the majority of the variance in occurrence. This study aimed to identify novel susceptibility loci for the development of HCC in people with alcohol related cirrhosis., Design: Patients with alcohol-related cirrhosis and HCC (cases: n=1214) and controls without HCC (n=1866), recruited from Germany, Austria, Switzerland, Italy and the UK, were included in a two-stage genome-wide association study using a case-control design. A validation cohort of 1520 people misusing alcohol but with no evidence of liver disease was included to control for possible association effects with alcohol misuse. Genotyping was performed using the InfiniumGlobal Screening Array (V.24v2, Illumina) and the OmniExpress Array (V.24v1-0a, Illumina)., Results: Associations with variants rs738409 in PNPLA3 and rs58542926 in TM6SF2 previously associated with an increased risk of HCC in patients with alcohol-related cirrhosis were confirmed at genome-wide significance. A novel locus rs2242652(A) in TERT (telomerase reverse transcriptase) was also associated with a decreased risk of HCC, in the combined meta-analysis, at genome-wide significance (p=6.41×10 -9 , OR=0.61 (95% CI 0.52 to 0.70). This protective association remained significant after correction for sex, age, body mass index and type 2 diabetes (p=7.94×10 -5 , OR=0.63 (95% CI 0.50 to 0.79). Carriage of rs2242652(A) in TERT was associated with an increased leucocyte telomere length (p=2.12×10 -44 )., Conclusion: This study identifies rs2242652 in TERT as a novel protective factor for HCC in patients with alcohol-related cirrhosis., Competing Interests: Competing interests: JT has received speaking and/or consulting fees from Versantis, Gore, Bayer, Alexion, Norgine, Grifols and CSL Behring., (© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY. Published by BMJ.)