Subgenomic RNA identification in SARS-CoV-2 genomic sequencing data.

We have developed periscope, a tool for the detection and quantification of subgenomic RNA (sgRNA) in SARS-CoV-2 genomic sequence data. The translation of the SARS-CoV-2 RNA genome for most open reading frames (ORFs) occurs via RNA intermediates termed "subgenomic RNAs." sgRNAs are produced through discontinuous transcription, which relies on homology between transcription regulatory sequences (TRS-B) upstream of the ORF start codons and that of the TRS-L, which is located in the 5' UTR. TRS-L is immediately preceded by a leader sequence. This leader sequence is therefore found at the 5' end of all sgRNA. We applied periscope to 1155 SARS-CoV-2 genomes from Sheffield, United Kingdom, and validated our findings using orthogonal data sets and in vitro cell systems. By using a simple local alignment to detect reads that contain the leader sequence, we were able to identify and quantify reads arising from canonical and noncanonical sgRNA. We were able to detect all canonical sgRNAs at the expected abundances, with the exception of ORF10. A number of recurrent noncanonical sgRNAs are detected. We show that the results are reproducible using technical replicates and determine the optimum number of reads for sgRNA analysis. In VeroE6 ACE2+/- cell lines, periscope can detect the changes in the kinetics of sgRNA in orthogonal sequencing data sets. Finally, variants found in genomic RNA are transmitted to sgRNAs with high fidelity in most cases. This tool can be applied to all sequenced COVID-19 samples worldwide to provide comprehensive analysis of SARS-CoV-2 sgRNA.

Real-World Outcomes of Direct-Acting Antiviral Treatment and Retreatment in United Kingdom-Based Patients Infected With Hepatitis C Virus Genotypes/Subtypes Endemic in Africa.

BACKGROUND: Chronic hepatitis C virus (HCV) infection affects 71 million individuals, mostly residing in low- and middle-income countries (LMICs). Direct-acting antivirals (DAAs) give high rates of sustained virological response (SVR) in high-income countries where a restricted range of HCV genotypes/subtypes circulate. METHODS: We studied United Kingdom-resident patients born in Africa to examine DAA effectiveness in LMICs where there is far greater breadth of HCV genotypes/subtypes. Viral genome sequences were determined from 233 patients. RESULTS: Full-length viral genomic sequences for 26 known subtypes and 5 previously unidentified isolates covering 5 HCV genotypes were determined. From 149 patients who received DAA treatment/retreatment, the overall SVR was 93%. Treatment failure was associated primarily with 2 subtypes, gt1l and gt4r, using sofosbuvir/ledipasvir. These subtypes contain natural resistance-associated variants that likely contribute to poor efficacy with this drug combination. Treatment failure was also significantly associated with hepatocellular carcinoma. CONCLUSIONS: DAA combinations give high SVR rates despite the high HCV diversity across the African continent except for subtypes gt1l and gt4r, which respond poorly to sofosbuvir/ledipasvir. These subtypes are widely distributed across Western, Central, and Eastern Africa. Thus, in circumstances where accurate genotyping is absent, ledipasvir and its generic compounds should not be considered as a recommended treatment option.

Real world SOF/VEL/VOX retreatment outcomes and viral resistance analysis for HCV patients with prior failure to DAA therapy.

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.

Amino Acid Substitutions in Genotype 3a Hepatitis C Virus Polymerase Protein Affect Responses to Sofosbuvir.

BACKGROUND & AIMS: Sofosbuvir is a frequently used pan-genotype inhibitor of hepatitis C virus (HCV) polymerase. This drug eliminates most chronic HCV infections, and resistance-associated substitutions in the polymerase are rare. However, HCV genotype 3 responds slightly less well to sofosbuvir-based therapies than other genotypes. We collected data from England's National Health Service Early Access Program to search for virus factors associated with sofosbuvir treatment failure. METHODS: We collected patient serum samples and used the capture-fusion assay to assess viral sensitivity to sofosbuvir in 14 HCV genotype 3 samples. We identified polymorphisms associated with reduced response and created modified forms of HCV and replicons containing the substitutions of interest and tested their sensitivity to sofosbuvir and ribavirin. We examined the effects of these polymorphisms by performing logistic regression multivariate analysis on their association with sustained virologic response in a separate cohort of 411 patients with chronic HCV genotype 3 infection who had been treated with sofosbuvir and ribavirin, with or without pegylated interferon. RESULTS: We identified a substitution in the HCV genotype 3a NS5b polymerase at amino acid 150 (alanine [A] to valine [V]), V at position 150 was observed in 42% of patients) with a reduced response to sofosbuvir in virus replication assays. In patients treated with sofosbuvir-containing regimens, the A150V variant was associated with a reduced response to treatment with sofosbuvir and ribavirin, with or without pegylated interferon. In 326 patients with V at position 150, 71% achieved an sustained virologic response compared to 88% with A at position 150. In cells, V at position 150 reduced the response to sofosbuvir 7-fold. We found that another rare substitution, glutamic acid (E) at position 206, significantly reduced the response to sofosbuvir (8.34-fold reduction); the combinations of V at position 150 and E at position 206 reduced the virus response to sofosbuvir 35.77-fold. Additionally, in a single patient, we identified 5 rare polymorphisms that reduced sensitivity to sofosbuvir our cell system. CONCLUSIONS: A common polymorphism, V at position 150 in the HCV genotype 3a NS5b polymerase, combined with other variants, reduces the virus response to sofosbuvir. Clinically, infection with HCV genotype 3 containing this variant reduces odds of sustained virologic response. In addition, we identified rare combinations of variants in HCV genotype 3 that reduce response to sofosbuvir.

A polymorphic residue that attenuates the antiviral potential of interferon lambda 4 in hominid lineages.

As antimicrobial signalling molecules, type III or lambda interferons (IFNλs) are critical for defence against infection by diverse pathogens, including bacteria, fungi and viruses. Counter-intuitively, expression of one member of the family, IFNλ4, is associated with decreased clearance of hepatitis C virus (HCV) in the human population; by contrast, a natural frameshift mutation that abrogates IFNλ4 production improves HCV clearance. To further understand how genetic variation between and within species affects IFNλ4 function, we screened a panel of all known extant coding variants of human IFNλ4 for their antiviral potential and identify three that substantially affect activity: P70S, L79F and K154E. The most notable variant was K154E, which was found in African Congo rainforest 'Pygmy' hunter-gatherers. K154E greatly enhanced in vitro activity in a range of antiviral (HCV, Zika virus, influenza virus and encephalomyocarditis virus) and gene expression assays. Remarkably, E154 is the ancestral residue in mammalian IFNλ4s and is extremely well conserved, yet K154 has been fixed throughout evolution of the hominid genus Homo, including Neanderthals. Compared to chimpanzee IFNλ4, the human orthologue had reduced activity due to amino acid K154. Comparison of published gene expression data from humans and chimpanzees showed that this difference in activity between K154 and E154 in IFNλ4 correlates with differences in antiviral gene expression in vivo during HCV infection. Mechanistically, our data show that the human-specific K154 negatively affects IFNλ4 activity through a novel means by reducing its secretion and potency. We thus demonstrate that attenuated activity of IFNλ4 is conserved among humans and postulate that differences in IFNλ4 activity between species contribute to distinct host-specific responses to-and outcomes of-infection, such as HCV infection. The driver of reduced IFNλ4 antiviral activity in humans remains unknown but likely arose between 6 million and 360,000 years ago in Africa.

Viral genotype correlates with distinct liver gene transcription signatures in chronic hepatitis C virus infection.

BACKGROUND: Chronic hepatitis C virus (HCV) infection of the liver with either genotype 1 or genotype 3 gives rise to distinct pathologies, and the two viral genotypes respond differently to antiviral therapy. METHODS: To understand these clinical differences, we compared gene transcription profiles in liver biopsies from patients infected with either gt1 or gt3, and uninfected controls. RESULTS: Gt1-infected biopsies displayed elevated levels of transcripts regulated by type I and type III interferons (IFN), including genes that predict response to IFN-α therapy. In contrast, genes controlled by IFN-γ were induced in gt3-infected biopsies. Moreover, IFN-γ levels were higher in gt3-infected biopsies. Analysis of hepatocyte-derived cell lines confirmed that the genes upregulated in gt3 infection were preferentially induced by IFN-γ. The transcriptional profile of gt3 infection was unaffected by IFNL4 polymorphisms, providing a rationale for the reduced predictive power of IFNL genotyping in gt3-infected patients. CONCLUSIONS: The interactions between HCV genotypes 1 and 3 and hepatocytes are distinct. These unique interactions provide avenues to explore the biological mechanisms that drive viral genotype-specific differences in disease progression and treatment response. A greater understanding of the distinct host-pathogen interactions of the different HCV genotypes is required to facilitate optimal management of HCV infection.

Global gene transcriptome analysis in vaccinated cattle revealed a dominant role of IL-22 for protection against bovine tuberculosis.

Bovine tuberculosis (bTB) is a chronic disease of cattle caused by Mycobacterium bovis, a member of the Mycobacterium tuberculosis complex group of bacteria. Vaccination of cattle might offer a long-term solution for controlling the disease and priority has been given to the development of a cattle vaccine against bTB. Identification of biomarkers in tuberculosis research remains elusive and the goal is to identify host correlates of protection. We hypothesized that by studying global gene expression we could identify in vitro predictors of protection that could help to facilitate vaccine development. Calves were vaccinated with BCG or with a heterologous BCG prime adenovirally vectored subunit boosting protocol. Protective efficacy was determined after M. bovis challenge. RNA was prepared from PPD-stimulated PBMC prepared from vaccinated-protected, vaccinated-unprotected and unvaccinated control cattle prior to M. bovis challenge and global gene expression determined by RNA-seq. 668 genes were differentially expressed in vaccinated-protected cattle compared with vaccinated-unprotected and unvaccinated control cattle. Cytokine-cytokine receptor interaction was the most significant pathway related to this dataset with IL-22 expression identified as the dominant surrogate of protection besides INF-γ. Finally, the expression of these candidate genes identified by RNA-seq was evaluated by RT-qPCR in an independent set of PBMC samples from BCG vaccinated and unvaccinated calves. This experiment confirmed the importance of IL-22 as predictor of vaccine efficacy.

The Transmission Route and Selection Pressure in HCV Subtype 3a and 3b Chinese Infections: Evolutionary Kinetics and Selective Force Analysis.

Hepatitis C virus (HCV) genotype 3 (GT-3) represents 22-30% of all infections and is the second most common genotype among all HCV genotypes. It has two main subtypes, GT-3a and GT-3b, that present epidemiological differences in transmission groups. This report generated 56 GT-3a and 64 GT-3b whole-genome sequences to conduct an evolutionary kinetics and selective force analysis with reference sequences from various countries. Evolutionary analysis showed that HCV GT-3a worldwide might have been transmitted from the Indian subcontinent to South Asia, Europe, North America and then become endemic in China. In China, GT-3a may have been transmitted by intravenous drug users (IDUs) and become endemic in the general population, while GT-3b may have originated from IDUs and then underwent mutual transmission between blood donors (BDs) and IDUs, ultimately becoming independently endemic in IDUs. Furthermore, the spread of GT-3a and GT-3b sequences from BD and IDU populations exhibit different selective pressures: the proportion of positively selected sites (PPSs) in E1 and E2 from IDUs was higher than in BDs. The number of positive selection sites was higher in GT-3b and IDUs. These results indicate that different selective constraints act along with the GT-3a and GT-3b genomes from IDUs and BDs. In addition, GT-3a and GT-3b have different transmission routes in China, which allows us to formulate specific HCV prevention and control strategies in China.

A Human and Rhesus Macaque Interferon-Stimulated Gene Screen Shows That Over-Expression of ARHGEF3/XPLN Inhibits Replication of Hepatitis C Virus and Other Flavivirids.

Natural hepatitis C virus (HCV) infection is restricted to humans, whereas other primates such as rhesus macaques are non-permissive for infection. To identify human and rhesus macaque genes that differ or share the ability to inhibit HCV replication, we conducted a medium-throughput screen of lentivirus-expressed host genes that disrupt replication of HCV subgenomic replicon RNA expressing secreted Gaussia luciferase. A combined total of >800 interferon-stimulated genes (ISGs) were screened. Our findings confirmed established anti-HCV ISGs, such as IRF1, PKR and DDX60. Novel species−specific inhibitors were also identified and independently validated. Using a cell-based system that recapitulates productive HCV infection, we identified that over-expression of the 'Rho Guanine Nucleotide Exchange Factor 3' gene (ARHGEF3) from both species inhibits full-length virus replication. Additionally, replication of two mosquito-borne flaviviruses, yellow fever virus (YFV) and Zika virus (ZIKV), were also reduced in cell lines over-expressing ARHGEF3 compared to controls. In conclusion, we ascribe novel antiviral activity to the cellular gene ARHGEF3 that inhibits replication of HCV and other important human viral pathogens belonging to the Flaviviridae, and which is conserved between humans and rhesus macaques.

The evolutionary dynamics and epidemiological history of hepatitis C virus genotype 6, including unique strains from the Li community of Hainan Island, China.

Hepatitis C virus (HCV) is a highly diverse pathogen that frequently establishes a chronic long-term infection, but the origins and drivers of HCV diversity in the human population remain unclear. Previously unidentified strains of HCV genotype 6 (gt6) were recently discovered in chronically infected individuals of the Li ethnic group living in Baisha County, Hainan Island, China. The Li community, who were early settlers on Hainan Island, has a distinct host genetic background and cultural identity compared to other ethnic groups on the island and mainland China. In this report, we generated 33 whole virus genome sequences to conduct a comprehensive molecular epidemiological analysis of these novel gt6 strains in the context of gt6 isolates present in Southeast Asia. With the exception of one gt6a isolate, the Li gt6 sequences formed three novel clades from two lineages which constituted 3 newly assigned gt6 subtypes and 30 unassigned strains. Using Bayesian inference methods, we dated the most recent common ancestor for all available gt6 whole virus genome sequences to approximately 2767 bce (95 per cent highest posterior density (HPD) intervals, 3670-1397 bce), which is far earlier than previous estimates. The substitution rate was 1.20 × 10-4 substitutions/site/year (s/s/y), and this rate varied across the genome regions, from 1.02 × 10-5 s/s/y in the 5'untranslated region (UTR) region to 3.07 × 10-4 s/s/y in E2. Thus, our study on an isolated ethnic minority group within a small geographical area of Hainan Island has substantially increased the known diversity of HCV gt6, already acknowledged as the most diverse HCV genotype. The extant HCV gt6 sequences from this study were probably transmitted to the Li through at least three independent events dating perhaps from around 4,000 years ago. This analysis describes deeper insight into basic aspects of HCV gt6 molecular evolution including the extensive diversity of gt6 sequences in the isolated Li ethnic group.

Comparison of SARS-CoV-2 Evolution in Paediatric Primary Airway Epithelial Cell Cultures Compared with Vero-Derived Cell Lines.

SARS-CoV-2 can efficiently infect both children and adults, albeit with morbidity and mortality positively associated with increasing host age and presence of co-morbidities. SARS-CoV-2 continues to adapt to the human population, resulting in several variants of concern (VOC) with novel properties, such as Alpha and Delta. However, factors driving SARS-CoV-2 fitness and evolution in paediatric cohorts remain poorly explored. Here, we provide evidence that both viral and host factors co-operate to shape SARS-CoV-2 genotypic and phenotypic change in primary airway cell cultures derived from children. Through viral whole-genome sequencing, we explored changes in genetic diversity over time of two pre-VOC clinical isolates of SARS-CoV-2 during passage in paediatric well-differentiated primary nasal epithelial cell (WD-PNEC) cultures and in parallel, in unmodified Vero-derived cell lines. We identified a consistent, rich genetic diversity arising in vitro, variants of which could rapidly rise to near fixation within two passages. Within isolates, SARS-CoV-2 evolution was dependent on host cells, with paediatric WD-PNECs showing a reduced diversity compared to Vero (E6) cells. However, mutations were not shared between strains. Furthermore, comparison of both Vero-grown isolates on WD-PNECs disclosed marked growth attenuation mapping to the loss of the polybasic cleavage site (PBCS) in Spike, while the strain with mutations in Nsp12 (T293I), Spike (P812R) and a truncation of Orf7a remained viable in WD-PNECs. Altogether, our work demonstrates that pre-VOC SARS-CoV-2 efficiently infects paediatric respiratory epithelial cells, and its evolution is restrained compared to Vero (E6) cells, similar to the case of adult cells. We highlight the significant genetic plasticity of SARS-CoV-2 while uncovering an influential role for collaboration between viral and host cell factors in shaping viral evolution and ultimately fitness in human respiratory epithelium.

Genomic epidemiology reveals multiple introductions of SARS-CoV-2 from mainland Europe into Scotland.

Coronavirus disease 2019 (COVID-19) was first diagnosed in Scotland on 1 March 2020. During the first month of the outbreak, 2,641 cases of COVID-19 led to 1,832 hospital admissions, 207 intensive care admissions and 126 deaths. We aimed to identify the source and number of introductions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into Scotland using a combined phylogenetic and epidemiological approach. Sequencing of 1,314 SARS-CoV-2 viral genomes from available patient samples enabled us to estimate that SARS-CoV-2 was introduced to Scotland on at least 283 occasions during February and March 2020. Epidemiological analysis confirmed that early introductions of SARS-CoV-2 originated from mainland Europe (the majority from Italy and Spain). We identified subsequent early outbreaks in the community, within healthcare facilities and at an international conference. Community transmission occurred after 2 March, 3 weeks before control measures were introduced. Earlier travel restrictions or quarantine measures, both locally and internationally, would have reduced the number of COVID-19 cases in Scotland. The risk of multiple reintroduction events in future waves of infection remains high in the absence of population immunity.

Genetic epidemiology of SARS-CoV-2 transmission in renal dialysis units - A high risk community-hospital interface.

OBJECTIVES: Patients requiring haemodialysis are at increased risk of serious illness with SARS-CoV-2 infection. To improve the understanding of transmission risks in six Scottish renal dialysis units, we utilised the rapid whole-genome sequencing data generated by the COG-UK consortium. METHODS: We combined geographical, temporal and genomic sequence data from the community and hospital to estimate the probability of infection originating from within the dialysis unit, the hospital or the community using Bayesian statistical modelling and compared these results to the details of epidemiological investigations. RESULTS: Of 671 patients, 60 (8.9%) became infected with SARS-CoV-2, of whom 16 (27%) died. Within-unit and community transmission were both evident and an instance of transmission from the wider hospital setting was also demonstrated. CONCLUSIONS: Near-real-time SARS-CoV-2 sequencing data can facilitate tailored infection prevention and control measures, which can be targeted at reducing risk in these settings.

Interpreting Viral Deep Sequencing Data with GLUE.

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.

Interferon lambda 4 impacts the genetic diversity of hepatitis C virus.

Hepatitis C virus (HCV) is a highly variable pathogen that frequently establishes chronic infection. This genetic variability is affected by the adaptive immune response but the contribution of other host factors is unclear. Here, we examined the role played by interferon lambda-4 (IFN-λ4) on HCV diversity; IFN-λ4 plays a crucial role in spontaneous clearance or establishment of chronicity following acute infection. We performed viral genome-wide association studies using human and viral data from 485 patients of white ancestry infected with HCV genotype 3a. We demonstrate that combinations of host genetic variants, which determine IFN-λ4 protein production and activity, influence amino acid variation across the viral polyprotein - not restricted to specific viral proteins or HLA restricted epitopes - and modulate viral load. We also observed an association with viral di-nucleotide proportions. These results support a direct role for IFN-λ4 in exerting selective pressure across the viral genome, possibly by a novel mechanism.

Biomarkers of tuberculosis: a research roadmap.

Tuberculosis (TB) continues to represent a major public health problem worldwide. Prompt and accurate diagnosis and effective treatment are fundamental to reducing morbidity and mortality and curtailing spread of infection. Furthermore, tackling the large reservoir of latent infection is the cornerstone to TB control in many high income low TB incidence countries. However, our existing toolkit for prevention, diagnosis and treatment remains outdated and inadequate. Here, we discuss the key targets for biomarker research and discovery in TB and recent developments in the field. We focus on host biomarkers, in particular: correlates of vaccine efficacy and sterilizing immunity; biomarkers of latent TB infection, including diagnosis, risk of progression to active TB and response to treatment; and markers of active TB, including diagnosis, response to treatment and risk of relapse. Recent scientific and technological advances have contributed to significant recent progression in biomarker discovery. Although there are clear remaining paucities, continued efforts within scientific, translational and clinical studies are likely to yield a number of clinically useful biomarkers of TB in the foreseeable future.