Search

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]

Respiratory syncytial virus (RSV) is the leading cause of hospitalisation for respiratory infection in young children. RSV disease severity is known to be age-dependent and highest in young infants, but other correlates of severity, particularly the presence of additional respiratory pathogens, are less well understood. In this study, nasopharyngeal swabs were collected from two cohorts of RSV-positive infants <12 months in Spain, the UK, and the Netherlands during 2017–20. We show, using targeted metagenomic sequencing of >100 pathogens, including all common respiratory viruses and bacteria, from samples collected from 433 infants, that burden of additional viruses is common (111/433, 26%) but only modestly correlates with RSV disease severity. In contrast, there is strong evidence in both cohorts and across age groups that presence of Haemophilus bacteria (194/433, 45%) is associated with higher severity, including much higher rates of hospitalisation (odds ratio 4.25, 95% CI 2.03–9.31). There is no evidence for association between higher severity and other detected bacteria, and no difference in severity between RSV genotypes. Our findings reveal the genomic diversity of additional pathogens during RSV infection in infants, and provide an evidence base for future causal investigations of the impact of co-infection on RSV disease severity. The impact of other pathogens on disease outcome was studied in European infants with RSV infection. Additional viruses were commonly co-detected during infection but were weakly linked to severity. However, presence of Haemophilus bacteria strongly associated with severe cases. [ABSTRACT FROM AUTHOR]

Background: Pakistan has one of the highest burdens of Hepatitis C virus (HCV) infection globally. To achieve the World Health Organization's goals for HCV elimination, there is a need for substantial scale-up in testing, treatment, and a reduction in new infections. Data on the population impact of scaling up treatment is not available in Pakistan, nor is there reliable data on the incidence of infection/reinfection. This project will fill this gap by providing important empirical data on the incidence of infection (primary and reinfection) in Pakistan. Then, by using this data in epidemic models, the study will determine whether response rates achieved with affordable therapies (sofosbuvir plus daclatasvir) will be sufficient to eliminate HCV in Pakistan. Methods: This prospective multi-centre cohort study will screen 25,000 individuals for HCV antibody (Ab) and RNA (if Ab-positive) at various centers in Pakistan- Karachi (Sindh) and Punjab, providing estimates of the disease prevalence. HCV positive patients will be treated with sofosbuvir and daclatasvir for 12-weeks, (extended to 24-weeks in those with cirrhosis) and the proportion responding to this first-line treatment estimated. Patients who test HCV Ab negative will be recalled 12 months later to test for new HCV infections, providing estimates of the incidence rate. Patients diagnosed with HCV (~ 4,000) will be treated and tested for Sustained Virological Response (SVR). Questionnaires to assess risk factors, productivity, health care usage and quality of life will be completed at both the initial screening and at 12-month follow-up, allowing mathematical modelling and economic analysis to assess the current treatment strategies. Viral resistance will be analysed and patients who have successfully completed treatment will be retested 12 months later to estimate the rate of re-infection. Conclusion: The HepFREEPak study will provide evidence on the efficacy of available and widely used treatment options in Pakistan. It will also provide data on the incidence rate of primary infections and re-infections. Data on incidence risk factors will allow us to model and incorporate heterogeneity of risk and how that affects screening and treatment strategies. These data will identify any gaps in current test-and-treat programs to achieve HCV elimination in Pakistan. Study registration: This study was registered on clinicaltrials.gov (NCT04943588) on June 29, 2021. [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]

Human respiratory syncytial virus (RSV) is a major cause of lower respiratory tract infection in young children globally, but little is known about within-host RSV diversity. Here, we characterised within-host RSV populations using deep-sequencing data from 319 nasopharyngeal swabs collected during 2017–2020. RSV-B had lower consensus diversity than RSV-A at the population level, while exhibiting greater within-host diversity. Two RSV-B consensus sequences had an amino acid alteration (K68N) in the fusion (F) protein, which has been associated with reduced susceptibility to nirsevimab (MEDI8897), a novel RSV monoclonal antibody under development. In addition, several minor variants were identified in the antigenic sites of the F protein, one of which may confer resistance to palivizumab, the only licensed RSV monoclonal antibody. The differences in within-host virus populations emphasise the importance of monitoring for vaccine efficacy and may help to explain the different prevalences of monoclonal antibody-escape mutants between the two subgroups. Respiratory syncytial virus (RSV) is a common infection in children and older adults but little is known about within-host viral population diversity. Here, the authors perform deep sequencing and find that RSV subgroup B exhibited more diversity than subgroup A, with implications for development of therapeutics and vaccines. [ABSTRACT FROM AUTHOR]

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. [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]

Epstein-Barr virus (EBV) reactivation is common in sepsis patients but the extent and nature of this remains unresolved. We sought to determine the incidence and correlates of EBV-positivity in a large sepsis cohort. We also hypothesised that EBV reactivation would be increased in patients in whom relative immunosuppression was the major feature of their sepsis response. To identify such patients we aimed to use knowledge of sepsis response subphenotypes based on transcriptomic studies of circulating leukocytes, specifically patients with a Sepsis Response Signature endotype (SRS1) that we have previously shown to be associated with increased mortality and features of immunosuppression. We assayed EBV from the plasma of intensive care unit (ICU) patients with sepsis due to community-acquired pneumonia. In total 730 patients were evaluated by targeted metagenomics (n = 573 patients), digital droplet PCR (n = 565), or both (n = 408). We had previously analysed gene expression in peripheral blood leukocytes for a subset of individuals (n = 390). We observed a 37% incidence of EBV-positivity. EBV reactivation was associated with longer ICU stay (12.9 vs 9.2 days; p = 0.004) and increased organ failure (day 1 SOFA score 6.9 vs 5.9; p = 0.00011). EBV reactivation was associated with the relatively immunosuppressed SRS1 endotype (p = 0.014) and differential expression of a small number of biologically relevant genes. These findings are consistent with the hypothesis that viral reactivation in sepsis is a consequence of immune compromise and is associated with increasing severity of illness although further mechanistic studies are required to definitively illustrate cause and effect. [ABSTRACT FROM AUTHOR]

Tree structures, showing hierarchical relationships and the latent structures between samples, are ubiquitous in genomic and biomedical sciences. A common question in many studies is whether there is an association between a response variable measured on each sample and the latent group structure represented by some given tree. Currently, this is addressed on an ad hoc basis, usually requiring the user to decide on an appropriate number of clusters to prune out of the tree to be tested against the response variable. Here, we present a statistical method with statistical guarantees that tests for association between the response variable and a fixed tree structure across all levels of the tree hierarchy with high power while accounting for the overall false positive error rate. This enhances the robustness and reproducibility of such findings. [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]

In eukaryotes, detailed surveys of recombination rates have shown variation at multiple genomic scales and the presence of “hotspots” of highly elevated recombination. In bacteria, studies of recombination rate variation are less developed, in part because there are few analysis methods that take into account the clonal context within which bacterial evolution occurs. Here, we focus in particular on identifying “hot regions” of the genome where DNA is transferred frequently between isolates. We present a computationally efficient algorithm based on the recently developed “chromosome painting” algorithm, which characterizes patterns of haplotype sharing across a genome. We compare the average genome wide painting, which principally reflects clonal descent, with the painting for each site which additionally reflects the specific deviations at the site due to recombination. Using simulated data, we show that hot regions have consistently higher deviations from the genome wide average than normal regions. We applied our approach to previously analyzed Escherichia coli genomes and revealed that the new method is highly correlated with the number of recombination events affecting each site inferred by ClonalOrigin, a method that is only applicable to small numbers of genomes. Furthermore, we analyzed recombination hot regions in Campylobacter jejuni by using 200 genomes. We identified three recombination hot regions, which are enriched for genes related to membrane proteins. Our approach and its implementation, which is downloadable from https://github.com/bioprojects/orderedPainting, will help to develop a new phase of population genomic studies of recombination in prokaryotes. [ABSTRACT FROM PUBLISHER]

Background. Clostridium difficile is a major cause of nosocomial diarrhea, with 30-day mortality reaching 30%. The cell surface comprises a paracrystalline proteinaceous S-layer encoded by the slpA gene within the cell wall protein (cwp) gene cluster. Our purpose was to understand the diversity and evolution of slpA and nearby genes also encoding immunodominant cell surface antigens.Methods. Whole-genome sequences were determined for 57 C. difficile isolates representative of the population structure and different clinical phenotypes. Phylogenetic analyses were performed on their genomic region (>63 kb) spanning the cwp cluster.Results. Genetic diversity across the cwp cluster peaked within slpA, cwp66 (adhesin), and secA2 (secretory translocase). These genes formed a 10-kb cassette, of which 12 divergent variants were found. Homologous recombination involving this cassette caused it to associate randomly with genotype. One cassette contained a novel insertion (length, approximately 24 kb) that resembled S-layer glycosylation gene clusters.Conclusions. Genetic exchange of S-layer cassettes parallels polysaccharide capsular switching in other species. Both cause major antigenic shifts, while the remainder of the genome is unchanged. C. difficile genotype is therefore not predictive of antigenic type. S-layer switching and immune escape could help explain temporal and geographic variation in C. difficile epidemiology and may inform genotyping and vaccination strategies. [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]

Advancing interventions to tackle the huge global burden of hepatitis B virus (HBV) infection depends on improved insights into virus epidemiology, transmission, within-host diversity, drug resistance and pathogenesis, all of which can be advanced through the large-scale generation of full-length virus genome data. Here we describe advances to a protocol that exploits the circular HBV genome structure, using isothermal rolling-circle amplification to enrich HBV DNA, generating concatemeric amplicons containing multiple successive copies of the same genome. We show that this product is suitable for Nanopore sequencing as single reads, as well as for generating short-read Illumina sequences. Nanopore reads can be used to implement a straightforward method for error correction that reduces the per-read error rate, by comparing multiple genome copies combined into a single concatemer and by analysing reads generated from plus and minus strands. With this approach, we can achieve an improved consensus sequencing accuracy of 99.7% and resolve intra-sample sequence variants to form whole-genome haplotypes. Thus while Illumina sequencing may still be the most accurate way to capture within-sample diversity, Nanopore data can contribute to an understanding of linkage between polymorphisms within individual virions. The combination of isothermal amplification and Nanopore sequencing also offers appealing potential to develop point-of-care tests for HBV, and for other viruses. [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]

Background: Sustainable Development Goals set a challenge for the elimination of hepatitis B virus (HBV) infection as a public health concern by the year 2030. Deployment of a robust prophylactic vaccine and enhanced interventions for prevention of mother to child transmission (PMTCT) are cornerstones of elimination strategy. However, in light of the estimated global burden of 290 million cases, enhanced efforts are required to underpin optimisation of public health strategy. Robust analysis of population epidemiology is particularly crucial for populations in Africa made vulnerable by HIV co-infection, poverty, stigma and poor access to prevention, diagnosis and treatment.Methods: We here set out to evaluate the current and future role of HBV vaccination and PMTCT as tools for elimination. We first investigated the current impact of paediatric vaccination in a cohort of children with and without HIV infection in Kimberley, South Africa. Second, we used these data to inform a new parsimonious model to simulate the ongoing impact of preventive interventions. By applying these two approaches in parallel, we are able to determine both the current impact of interventions, and the future projected outcome of ongoing preventive strategies over time.Results: Existing efforts have been successful in reducing paediatric prevalence of HBV infection in this setting to < 1%, demonstrating the success of the existing vaccine campaign. Our model predicts that, if consistently deployed, combination efforts of vaccination and PMTCT can significantly reduce population prevalence (HBsAg) by 2030, such that a major public health impact is possible even without achieving elimination. However, the prevalence of HBV e-antigen (HBeAg)-positive carriers will decline more slowly, representing a persistent population reservoir. We show that HIV co-infection significantly reduces titres of vaccine-mediated antibody, but has a relatively minor role in influencing the projected time to elimination. Our model can also be applied to other settings in order to predict impact and time to elimination based on specific interventions.Conclusions: Through extensive deployment of preventive strategies for HBV, significant positive public health impact is possible, although time to HBV elimination as a public health concern is likely to be substantially longer than that proposed by current goals. [ABSTRACT FROM AUTHOR]