Exploring the interactions between the human and viral genomes

Joint host-pathogen genomic analysis identifies hepatitis B virus mutations associated with human NTCP and HLA class I variation

Evolutionary changes in the hepatitis B virus (HBV) genome could reflect its adaptation to host-induced selective pressure. Leveraging paired human exome and ultra-deep HBV genome-sequencing data from 567 affected individuals with chronic hepatitis B, we comprehensively searched for the signatures of this evolutionary process by conducting "genome-to-genome" association tests between all human genetic variants and viral mutations. We identified significant associations between an East Asian-specific missense variant in the gene encoding the HBV entry receptor NTCP (rs2296651, NTCP S267F) and mutations within the receptor-binding region of HBV preS1. Through in silico modeling and in vitro preS1-NTCP binding assays, we observed that the associated HBV mutations are in proximity to the NTCP variant when bound and together partially increase binding affinity to NTCP S267F. Furthermore, we identified significant associations between HLA-A variation and viral mutations in HLA-A-restricted T cell epitopes. We used in silico binding prediction tools to evaluate the impact of the associated HBV mutations on HLA presentation and observed that mutations that result in weaker binding affinities to their cognate HLA alleles were enriched. Overall, our results suggest the emergence of HBV escape mutations that might alter the interaction between HBV PreS1 and its cellular receptor NTCP during viral entry into hepatocytes and confirm the role of HLA class I restriction in inducing HBV epitope variations.

Life-threatening infections in human newborns: Reconciling age-specific vulnerability and interindividual variability

In humans, the interindividual variability of clinical outcome following exposure to a microorganism is immense, ranging from silent infection to life-threatening disease. Age-specific immune responses partially account for the high incidence of infection during the first 28 days of life and the related high mortality at population level. However, the occurrence of life-threatening disease in individual newborns remains unexplained. By contrast, inborn errors of immunity and their immune phenocopies are increasingly being discovered in children and adults with life-threatening viral, bacterial, mycobacterial and fungal infections. There is a need for convergence between the fields of neonatal immunology, with its in-depth population-wide characterization of newborn-specific immune responses, and clinical immunology, with its investigations of infections in patients at the cellular and molecular levels, to facilitate identification of the mechanisms of susceptibility to infection in individual newborns and the design of novel preventive and therapeutic strategies.

Polymorphisms of CCR5, IL-6, IFN-γ and IL-10 genes in Cuban HIV/AIDS patients

INTRODUCTION

Genetic studies have shown associations of several single nucleotide polymorphisms (SNP) with different rates of progression and variation in susceptibility to HIV infection. This study aimed to estimate the frequency of ccr5Δ32, IL-6-174G/C, IFN-γ+874T/A and IL-10-1082A/G polymorphisms in Cuban HIV-infected patients and a group of sero-discordant couples to assess their influence on risk and disease progression.

METHODS

A cross-sectional study was carried out on 120 subjects registered at the Institute of Tropical Medicine «Pedro Kour» (IPK) and the Ameijeiras Hospital from June 2018 until December 2019. The amplification of fragments of the ccr5, IL-6, IFN-γ and IL-10 genes was performed by polymerase chain reaction followed by identification of polymorphisms using the restriction fragment length polymorphism analysis for IL-6 with the restriction enzymes Nla III. Amplification Refractory Mutation System was used for IFN-γ and IL-10 genes.

RESULTS

The allelic and genotypic distributions of the genes ccr5, IL-6, IFN-γ and IL-10 did not differ significantly between the two groups. Cell counts and plasma viral load values did not differ significantly between genotypes of the ccr5, IL-6, IFN-γ and IL-10 genes. Only the IL-6 GC genotype was associated with higher viral load values. The combination of alleles of the four considered SNPs showed a highly significant increase in the risk of HIV infection for one of them, but with a very low frequency (<1%).

CONCLUSION

This study contributes to evaluating the frequency of these polymorphisms and their influence on biomarkers of the progression of HIV infection in the Cuban HIV-population.

G2GSnake: a Snakemake workflow for host-pathogen genomic association studies

Summary

Joint analyses of paired host and pathogen genome sequences have the potential to enhance our understanding of host-pathogen interactions. A systematic approach to conduct such a joint analysis is through a "genome-to-genome" (G2G) association study, which involves testing for associations between all host and pathogen genetic variants. Significant associations reveal host genetic factors that might drive pathogen variation, highlighting biological mechanisms likely to be involved in host control and pathogen escape. Here, we present a Snakemake workflow that allows researchers to conduct G2G studies in a reproducible and scalable manner. In addition, we have developed an intuitive R Shiny application that generates custom summaries of the results, enabling users to derive relevant insights.

Availability and implementation

G2GSnake is freely available at: https://github.com/zmx21/G2GSnake under the MIT license.

Differential genome-wide associated variants and enriched pathways of ECG parameters among people with versus without HIV

OBJECTIVES

People with HIV (PWH) are more likely to develop ECG abnormalities. Substantial evidence exists for genetic contribution to ECG parameters among general population. However, whether and how would host genome associate with ECG parameters among PWH is unclear. Our research aims to analyze and compare genetic variants, mapped genes, and enriched pathways of ECG parameters among PWH and HIV-negative controls.

DESIGN

A cross-sectional study.

METHOD

We performed a large original genome-wide association study (GWAS) of ECG parameters among PWH ( n  = 1730) and HIV-negative controls ( n =  3746). Genome-wide interaction analyses were also conducted.

RESULTS

A total of 18 novel variants were detected among PWH, six for PR interval including rs76345397 at ATL2 , 11 for QRS duration including rs10483994 at KCNK10 and rs2478830 at JCAD , and one for QTc interval (rs9815364). Among HIV-negative controls, we identified variants located at previously reported ECG-related genes ( SCN5A , CNOT1 ). Genetic variants had a significant interaction with HIV infection ( P  < 5 × 10 -8 ), implying that HIV infection and host genome might jointly influence ECG parameters. Mapped genes for PR interval and QRS duration among PWH were enriched in the biological process of viral genome replication and host response to virus, respectively, whereas enriched pathways for PR interval among HIV-negative controls were in the cellular component of voltage-gated sodium channel complex.

CONCLUSION

The present GWAS indicated a distinctive impact of host genome on quantitative ECG parameters among PWH. Different from HIV-negative controls, host genome might influence the cardiac electrical activity by interfering with HIV viral infection, production, and latency among PWH.

Associations of functional human leucocyte antigen class I groups with HIV viral load in a heterogeneous cohort

OBJECTIVE

Human leucocyte antigen (HLA) class I alleles are the main host genetic factors involved in controlling HIV-1 viral load (VL). Nevertheless, HLA diversity has proven a significant challenge in association studies. We assessed how accounting for binding affinities of HLA class I alleles to HIV-1 peptides facilitate association testing of HLA with HIV-1 VL in a heterogeneous cohort.

DESIGN

Cohort from the Strategic Timing of AntiRetroviral Treatment (START) study.

METHODS

We imputed HLA class I alleles from host genetic data (2546 HIV+ participants) and sampled immunopeptidomes from 2079 host-paired viral genomes (targeted amplicon sequencing). We predicted HLA class I binding affinities to HIV-1 and unspecific peptides, grouping alleles into functional clusters through consensus clustering. These functional HLA class I clusters were used to test associations with HIV VL.

RESULTS

We identified four clades totaling 30 HLA alleles accounting for 11.4% variability in VL. We highlight HLA-B∗57:01 and B∗57:03 as functionally similar but yet overrepresented in distinct ethnic groups, showing when combined a protective association with HIV+ VL (log, β -0.25; adj. P-value < 0.05). We further demonstrate only a slight power reduction when using unspecific immunopeptidomes, facilitating the use of the inferred functional HLA groups in other studies.

CONCLUSION

The outlined computational approach provides a robust and efficient way to incorporate HLA function and peptide diversity, aiding clinical association studies in heterogeneous cohorts. To facilitate access to the proposed methods and results we provide an interactive application for exploring data.

Human immunotypes impose selection on viral genotypes through viral epitope specificity

BACKGROUND

Understanding the genetic interplay between human hosts and infectious pathogens is crucial for how we interpret virulence factors. Here, we tested for associations between HIV and host genetics, and interactive genetic effects on viral load (VL) in HIV+ ART-naive clinical trial participants.

METHODS

HIV genomes were sequenced and the encoded amino acid (AA) variants were associated with VL, human single nucleotide polymorphisms (SNPs) and imputed HLA alleles, using generalized linear models with Bonferroni correction.

RESULTS

Human (388,501 SNPs) and HIV (3,010 variants) genetic data was available for 2,122 persons. Four HIV variants were associated with VL (p-values<1.66×10 -5). Twelve HIV variants were associated with a range of 1-512 human SNPs (p-value<4.28×10 -11). We found 46 associations between HLA alleles and HIV variants (p-values<1.29×10 -7). We found HIV variants and immunotypes when analyzed separately, were associated with lower VL, whereas the opposite was true when analyzed in concert. Epitope binding prediction showed HLA alleles to be weaker binders of associated HIV AA variants relative to alternative variants on the same position.

CONCLUSIONS

Our results show the importance of immunotype specificity on viral antigenic determinants, and the identified genetic interplay puts emphasis that viral and human genetics should be studied in the context of each other.

Immunological lessons from genome-wide association studies of infections

Over the past few years, genome-wide association studies (GWAS) have been increasingly applied to identify host genetic factors influencing clinical and laboratory traits related to immunity and infection, and to understand the interplay between the host and the microbial genomes. By screening large cohorts of individuals suffering from various infectious diseases, GWAS explored resistance against infection, natural history of the disease, development of life-threatening clinical signs, and innate and adaptive immune responses. These efforts provided fundamental insight on the role of major genes in the interindividual variability in the response to infection and on the mechanisms of the immune response against human pathogens both at the individual and population levels.

HIV-1 and human genetic variation

Over the past four decades, research on the natural history of HIV infection has described how HIV wreaks havoc on human immunity and causes AIDS. HIV host genomic research, which aims to understand how human genetic variation affects our response to HIV infection, has progressed from early candidate gene studies to recent multi-omic efforts, benefiting from spectacular advances in sequencing technology and data science. In addition to invading cells and co-opting the host machinery for replication, HIV also stably integrates into our own genome. The study of the complex interactions between the human and retroviral genomes has improved our understanding of pathogenic mechanisms and suggested novel preventive and therapeutic approaches against HIV infection.

The human genetic determinism of life-threatening infectious diseases: genetic heterogeneity and physiological homogeneity?

Multicellular eukaryotes emerged late in evolution from an ocean of viruses, bacteria, archaea, and unicellular eukaryotes. These macroorganisms are exposed to and infected by a tremendous diversity of microorganisms. Those that are large enough can even be infected by multicellular fungi and parasites. Each interaction is unique, if only because it operates between two unique living organisms, in an infinite diversity of circumstances. This is neatly illustrated by the extraordinarily high level of interindividual clinical variability in human infections, even for a given pathogen, ranging from a total absence of clinical manifestations to death. We discuss here the idea that the determinism of human life-threatening infectious diseases can be governed by single-gene inborn errors of immunity, which are rarely Mendelian and frequently display incomplete penetrance. We briefly review the evidence in support of this notion obtained over the last two decades, referring to a number of focused and thorough reviews published by eminent colleagues in this issue of Human Genetics. It seems that almost any life-threatening infectious disease can be driven by at least one, and, perhaps, a great many diverse monogenic inborn errors, which may nonetheless be immunologically related. While the proportions of monogenic cases remain unknown, a picture in which genetic heterogeneity is combined with physiological homogeneity is emerging from these studies. A preliminary sketch of the human genetic architecture of severe infectious diseases is perhaps in sight.