Search

Abstract Drug type and dosing recommendation have been designed and optimized based on average response in the general population. Yet, there is significant inter‐individual variability in drug response, which results in treatment inefficacy or adverse drug reactions in a subset of patients. This is partly due to genetic factors that typically affect drug metabolism or clearance. To verify the relevance and applicability of international pharmacogenetic guidelines in the Swiss population, we genotyped 1533 patients from a hospital‐based biobank who received at least 30 different drugs, as documented in their electronic health record. We then assessed the prevalence of clinically actionable variants in 13 high‐risk pharmacogenes. We compared the allele frequencies obtained in the hospital‐based cohort with those of a Swiss population‐based cohort of 4791 individuals. The prevalence of clinically actionable variants was comparable between the two cohorts, with most study participants (97.3%) carrying at least one actionable pharmacogenetic variant. We then assessed the frequency of high‐risk prescriptions due to actionable gene–drug interactions and observed that 31% of patients in the hospital‐based cohort were prescribed at least one drug for which they carried a high‐risk variant, and for which international guidelines recommend a change of drug or dosage. Our analysis confirms the high prevalence of actionable pharmacogenetic variants in the Swiss population. It also shows that a substantial minority of patients are exposed to drugs for which they carry potentially problematic variants. Implementing a genetically informed approach to drug prescribing could have a positive impact on the quality of healthcare delivery.

Abstract The Mycobacterium tuberculosis complex (MTBC) comprises nine human-adapted lineages that differ in their geographical distribution. Local adaptation of specific MTBC genotypes to the respective human host population has been invoked in this context. We aimed to assess if bacterial genetics governs MTBC pathogenesis or if local co-adaptation translates into differential susceptibility of human macrophages to infection by different MTBC genotypes. We generated macrophages from cryopreserved blood mononuclear cells of Tanzanian tuberculosis patients, from which the infecting MTBC strains had previously been phylogenetically characterized. We infected these macrophages ex vivo with a phylogenetically similar MTBC strain (“matched infection”) or with strains representative of other MTBC lineages (“mismatched infection”). We found that L1 infections resulted in a significantly lower bacterial burden and that the intra-cellular replication rate of L2 strains was significantly higher compared the other MTBC lineages, irrespective of the MTBC lineage originally infecting the patients. Moreover, L4-infected macrophages released significantly greater amounts of TNF-α, IL-6, IL-10, MIP-1β, and IL-1β compared to macrophages infected by all other strains. While our results revealed no measurable effect of local adaptation, they further highlight the strong impact of MTBC phylogenetic diversity on the variable outcome of the host–pathogen interaction in human tuberculosis.

The pathogenesis of HIV-1 infection is governed by a highly dynamic, time-dependent interaction between the host and the viral genome. In this study, we developed a novel systematic approach to assess the host-virus interaction, using average pairwise viral diversity as a proxy for time since infection, and applied this method to nearly whole viral genome sequences (n = 4,464), human leukocyte antigen (HLA) genotyping data (n = 1,044), and viral RNA load (VL) measurements during the untreated chronic phase (n = 829) of Swiss HIV Cohort Study participants. Our systematic genome-wide screen revealed for 98 HLA/viral-variant pairs a signature of immune-driven selection in the form of an HLA-dependent effect of infection time on the presence of HIV amino acid variants. Of these pairs, 12 were found to have an effect on VL. Furthermore, 28/58 pairs were validated by time-to-event analyses and 48/92 by computational HLA-epitope predictions. Our diversity-based approach allows a powerful and systematic investigation of the interaction between the virus and cellular immunity, revealing a notable subset of such interaction effects. From an evolutionary perspective, these observations underscore the complexity of HLA-mediated selection pressures on the virus that shape viral evolution and pathogenesis.

Summary: Human genetic studies of critical COVID-19 pneumonia have revealed the essential role of type I interferon-dependent innate immunity to SARS-CoV-2 infection. Conversely, an association between the HLA-B∗15:01 allele and asymptomatic SARS-CoV-2 infection in unvaccinated individuals was recently reported, suggesting a contribution of pre-existing T cell-dependent adaptive immunity. We report a lack of association of classical HLA alleles, including HLA-B∗15:01, with pre-omicron asymptomatic SARS-CoV-2 infection in unvaccinated participants in a prospective population-based study in the United States (191 asymptomatic vs. 945 symptomatic COVID-19 cases). Moreover, we found no such association in the international COVID Human Genetic Effort cohort (206 asymptomatic vs. 574 mild or moderate COVID-19 cases and 1,625 severe or critical COVID-19 cases). Finally, in the Human Challenge Characterisation study, the three HLA-B∗15:01 individuals infected with SARS-CoV-2 developed symptoms. As with other acute primary infections studied, no classical HLA alleles favoring an asymptomatic course of SARS-CoV-2 infection were identified.

ABSTRACT Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) remains one of the deadliest infectious diseases globally. Timely diagnosis is a key step in the management of TB patients and in the prevention of further transmission events. Current diagnostic tools are limited in these regards. There is an urgent need for new accurate non-sputum-based diagnostic tools for the detection of symptomatic as well as subclinical TB. In this study, we recruited 52 symptomatic TB patients (sputum Xpert MTB/RIF positive) and 58 household contacts to assess the accuracy of a sequence-specific hybridization assay that detects the presence of Mtb cell-free DNA in urine. Using sputum Xpert MTB/RIF as a reference test, the magnetic bead-capture assay could discriminate active TB from healthy household contacts with an overall sensitivity of 72.1% [confidence interval (CI) 0.59–0.86] and specificity of 95.5% (CI 0.90–1.02) with a positive predictive value of 93.9% and negative predictive value of 78.2%. The detection of Mtb-specific DNA in urine suggested four asymptomatic TB infection cases that were confirmed in all instances either by concomitant Xpert MTB/RIF sputum testing or by follow-up investigation raising the specificity of the index test to 100%. We conclude that sequence-specific hybridization assays on urine specimens hold promise as non-invasive tests for the detection of subclinical TB.IMPORTANCEThere is an urgent need for a non-sputum-based diagnostic tool allowing sensitive and specific detection of all forms of tuberculosis (TB) infections. In that context, we performed a case-control study to assess the accuracy of a molecular detection method enabling the identification of cell-free DNA from Mycobacterium tuberculosis that is shed in the urine of tuberculosis patients. We present accuracy data that would fulfill the target product profile for a non-sputum test. In addition, recent epidemiological data suggested that up to 50% of individuals secreting live bacilli do not present with symptoms at the time of screening. We report, here, that the investigated index test could also detect instances of asymptomatic TB infections among household contacts.

IntroductionWe have reanalyzed the genomic data of the International Collaboration for the Genomics of HIV (ICGH), centering on HIV-1 Elite Controllers.MethodsWe performed a genome-wide Association Study comparing 543 HIV Elite Controllers with 3,272 uninfected controls of European descent. Using the latest database for imputation, we analyzed 35,552 Single Nucleotide Polymorphisms (SNPs) within the Major Histocompatibility Complex (MHC) region.ResultsOur analysis identified 2,626 SNPs significantly associated (p

Type 2 diabetes mellitus (T2DM) can be multifactorial where both genetics and environmental factors play a role. We aimed to investigate the use of polygenic risk scores (PRS) in the prediction of pre-transplant T2DM and post-transplant diabetes mellitus (PTDM) among solid organ transplant (SOT) patients. Using non-genetic risk scores alone; and the combination with PRS, separate logistic regression models were built and compared using receiver operator curves. Patients were assessed pre-transplant and in three post-transplant periods: 0–45, 46–365 and >365 days. A higher PRS was significantly associated with increased odds of pre-transplant T2DM. However, no improvement was observed for pre-transplant T2DM prediction when comparing PRS combined with non-genetic risk scores to using non-genetic risk scores alone. This was also true for predictions of PTDM in all three post-transplant periods. This study demonstrated that polygenic risk was only associated with the risk of T2DM among SOT recipients prior to transplant and not for PTDM. Combining PRS with a clinical model of non-genetic risk scores did not significantly improve the predictive ability, indicating its limited clinical utility in identifying patients at high risk for T2DM before transplantation, suggesting that non-genetic or different genetic factors may contribute to PTDM.

Abstract Background Aging in postmitotic tissues is associated with clonal expansion of somatic mitochondrial deletions, the origin of which is not well understood. Such deletions are often flanked by direct nucleotide repeats, but this alone does not fully explain their distribution. Here, we hypothesized that the close proximity of direct repeats on single-stranded mitochondrial DNA (mtDNA) might play a role in the formation of deletions. Results By analyzing human mtDNA deletions in the major arc of mtDNA, which is single-stranded during replication and is characterized by a high number of deletions, we found a non-uniform distribution with a “hot spot” where one deletion breakpoint occurred within the region of 6–9 kb and another within 13–16 kb of the mtDNA. This distribution was not explained by the presence of direct repeats, suggesting that other factors, such as the spatial proximity of these two regions, can be the cause. In silico analyses revealed that the single-stranded major arc may be organized as a large-scale hairpin-like loop with a center close to 11 kb and contacting regions between 6–9 kb and 13–16 kb, which would explain the high deletion activity in this contact zone. The direct repeats located within the contact zone, such as the well-known common repeat with a first arm at 8470–8482 bp (base pair) and a second arm at 13,447–13,459 bp, are three times more likely to cause deletions compared to direct repeats located outside of the contact zone. A comparison of age- and disease-associated deletions demonstrated that the contact zone plays a crucial role in explaining the age-associated deletions, emphasizing its importance in the rate of healthy aging. Conclusions Overall, we provide topological insights into the mechanism of age-associated deletion formation in human mtDNA, which could be used to predict somatic deletion burden and maximum lifespan in different human haplogroups and mammalian species.

Abstract Background We previously reported that impaired type I IFN activity, due to inborn errors of TLR3- and TLR7-dependent type I interferon (IFN) immunity or to autoantibodies against type I IFN, account for 15–20% of cases of life-threatening COVID-19 in unvaccinated patients. Therefore, the determinants of life-threatening COVID-19 remain to be identified in ~ 80% of cases. Methods We report here a genome-wide rare variant burden association analysis in 3269 unvaccinated patients with life-threatening COVID-19, and 1373 unvaccinated SARS-CoV-2-infected individuals without pneumonia. Among the 928 patients tested for autoantibodies against type I IFN, a quarter (234) were positive and were excluded. Results No gene reached genome-wide significance. Under a recessive model, the most significant gene with at-risk variants was TLR7, with an OR of 27.68 (95%CI 1.5–528.7, P = 1.1 × 10−4) for biochemically loss-of-function (bLOF) variants. We replicated the enrichment in rare predicted LOF (pLOF) variants at 13 influenza susceptibility loci involved in TLR3-dependent type I IFN immunity (OR = 3.70[95%CI 1.3–8.2], P = 2.1 × 10−4). This enrichment was further strengthened by (1) adding the recently reported TYK2 and TLR7 COVID-19 loci, particularly under a recessive model (OR = 19.65[95%CI 2.1–2635.4], P = 3.4 × 10−3), and (2) considering as pLOF branchpoint variants with potentially strong impacts on splicing among the 15 loci (OR = 4.40[9%CI 2.3–8.4], P = 7.7 × 10−8). Finally, the patients with pLOF/bLOF variants at these 15 loci were significantly younger (mean age [SD] = 43.3 [20.3] years) than the other patients (56.0 [17.3] years; P = 1.68 × 10−5). Conclusions Rare variants of TLR3- and TLR7-dependent type I IFN immunity genes can underlie life-threatening COVID-19, particularly with recessive inheritance, in patients under 60 years old.

Abstract Schizophrenia is a severe psychiatric disorder determined by a complex mixture of genetic and environmental factors. To better understand the contributions of human genetic variations to schizophrenia, we performed a genome-wide association study (GWAS) of a highly sensitive endophenotype. In this visual masking endophenotype, two vertical bars, slightly shifted in the horizontal direction, are briefly presented (vernier offset). Participants are asked to indicate the offset direction of the bars (either left or right). The bars are followed by a grating mask, which makes the task both spatially and temporally challenging. The inter-stimulus interval (ISI) between the vernier and the mask was determined in 206 patients with schizophrenia, 109 first-order relatives, and 143 controls. Usually, in GWAS studies, patients are compared to controls (i.e., a binary task) without considering the large differences in performance between patients and controls, as it occurs in many paradigms. The masking task allows for a particularly powerful analysis because the differences in ISI within the patient population are large. We genotyped all participants and searched for associations between human polymorphisms and the masking endophenotype using a linear mixed model. We did not identify any genome-wide significant associations (p

Abstract Multiple human pathogens establish chronic, sometimes life-long infections. Even if they are often latent, these infections can trigger some degree of local or systemic immune response, resulting in chronic low-grade inflammation. There remains an incomplete understanding of the potential contribution of both persistent infections and human genetic variation on chronic low-grade inflammation. We searched for potential associations between seropositivity for 13 persistent pathogens and the plasma levels of the inflammatory biomarker C-reactive protein (CRP), using data collected in the context of the UK Biobank and the CoLaus|PsyCoLaus Study, two large population-based cohorts. We performed backward stepwise regression starting with the following potential predictors: serostatus for each pathogen, polygenic risk score for CRP, and demographic and clinical factors known to be associated with CRP. We found evidence for an association between Chlamydia trachomatis (P-value = 5.04e − 3) and Helicobacter pylori (P-value = 8.63e − 4) seropositivity and higher plasma levels of CRP. We also found an association between pathogen burden and CRP levels (P-value = 4.12e − 4). These results improve our understanding of the relationship between persistent infections and chronic inflammation, an important determinant of long-term morbidity in humans.

Summary: Oxylipins are major immunomodulating mediators, yet studies of inflammation focus mainly on cytokines. Here, using a standardized whole-blood stimulation system, we characterized the oxylipin-driven inflammatory responses to various stimuli and their relationships with cytokine responses. We performed a pilot study in 25 healthy individuals using 6 different stimuli: 2 bacterial stimuli (LPS and live BCG), 2 viral stimuli (vaccine-grade poly I:C and live H1N1 attenuated influenza), an enterotoxin superantigen and a Null control. All stimuli induced a strong production of oxylipins but most importantly, bacterial, viral, and T cell immune responses show distinct oxylipin signatures. Integration of the oxylipin and cytokine responses for each condition revealed new immune networks improving our understanding of inflammation regulation. Finally, the oxylipin responses and oxylipin-cytokine networks were compared in patients with active tuberculosis or with latent infection. This revealed different responses to BCG but not LPS stimulation highlighting new regulatory pathways for further investigations.

Summary: Some HIV controllers experience immunologic progression with CD4+ T cell decline. We aimed to identify genetic factors associated with CD4+ T cell lost in HIV controllers. A total of 561 HIV controllers were included, 442 and 119 from the International HIV controllers Study Cohort and the Swiss HIV Cohort Study, respectively. No SNP or gene was associated with the long-term non-progressor HIV spontaneous control phenotype in the individual GWAS or in the meta-analysis. However, SNPs previously associated with natural HIV control linked to HLA-B (rs2395029 [p = 0.005; OR = 1.70], rs59440261 [p = 0.003; OR = 1.78]), MICA (rs112243036 [p = 0.011; OR = 1.45]), and PSORS1C1 loci (rs3815087 [p = 0.017; OR = 1.39]) showed nominal association with this phenotype. Genetic factors associated with the long-term HIV controllers without risk of immunologic progression are those previously related to the overall HIV controller phenotype.

Biological sex and host genetics influence HIV pathogenesis. Females have a higher likelihood of spontaneous viral control and lower set point viral load (spVL). No prior studies have assessed sex-specific genetics of HIV. To address this, we performed a sex-stratified genome-wide association study using data from the ICGH. Although it is the largest collection of genomic data in HIV, this multiethnic sample of 9,705 people is 81.3% male. We sought to identify sex-specific genetic variants and genes associated with HIV spVL and control. We confirmed associations in the HLA and CCR5 regions in males and HLA in females. Gene-based analyses detected associations between HIV spVL and PET100, PCP2, XAB2, and STXBP2 only in males. We detected variants with a significant sex-differential effect on spVL in SDC3 and PUM1 (rs10914268) and PSORS1C2 (rs1265159) and on HIV control in SUB1 (rs687659), AL158151.3, PTPA, and IER5L (rs4387067). Those variants have epigenetic and genetic interactions with relevant genes with both cis and trans effects. In summary, we identified sex-shared associations at the single-variant level, sex-specific associations at the gene-based level, and genetic variants with significant differential effects between the sexes.

In settings with high tuberculosis (TB) endemicity, distinct genotypes of the Mycobacterium tuberculosis complex (MTBC) often differ in prevalence. However, the factors leading to these differences remain poorly understood. Here we studied the MTBC population in Dar es Salaam, Tanzania over a six-year period, using 1,082 unique patient-derived MTBC whole-genome sequences (WGS) and associated clinical data. We show that the TB epidemic in Dar es Salaam is dominated by multiple MTBC genotypes introduced to Tanzania from different parts of the world during the last 300 years. The most common MTBC genotypes deriving from these introductions exhibited differences in transmission rates and in the duration of the infectious period, but little differences in overall fitness, as measured by the effective reproductive number. Moreover, measures of disease severity and bacterial load indicated no differences in virulence between these genotypes during active TB. Instead, the combination of an early introduction and a high transmission rate accounted for the high prevalence of L3.1.1, the most dominant MTBC genotype in this setting. Yet, a longer co-existence with the host population did not always result in a higher transmission rate, suggesting that distinct life-history traits have evolved in the different MTBC genotypes. Taken together, our results point to bacterial factors as important determinants of the TB epidemic in Dar es Salaam.

Abstract People living with human immunodeficiency virus (PLWH) have significantly increased risk for cardiovascular disease in part due to inflammation and immune dysregulation. Clonal hematopoiesis of indeterminate potential (CHIP), the age-related acquisition and expansion of hematopoietic stem cells due to leukemogenic driver mutations, increases risk for both hematologic malignancy and coronary artery disease (CAD). Since increased inflammation is hypothesized to be both a cause and consequence of CHIP, we hypothesized that PLWH have a greater prevalence of CHIP. We searched for CHIP in multi-ethnic cases from the Swiss HIV Cohort Study (SHCS, n = 600) and controls from the Atherosclerosis Risk in the Communities study (ARIC, n = 8111) from blood DNA-derived exome sequences. We observed that HIV is associated with a twofold increase in CHIP prevalence, both in the whole study population and in a subset of 230 cases and 1002 matched controls selected by propensity matching to control for demographic imbalances (SHCS 7%, ARIC 3%, p = 0.005). We also observed that ASXL1 is the most commonly mutated CHIP-associated gene in PLWH. Our results suggest that CHIP may contribute to the excess cardiovascular risk observed in PLWH.

Coronary heart disease (CHD) is one of the most pressing health problems of our time and a major cause of preventable death. CHD results from complex interactions between genetic and environmental factors. Using multiplex serological testing for persistent or frequently recurring infections and genome-wide analysis in a prospective population study, we delineate the respective and combined influences of genetic variation, infections, and low-grade inflammation on the risk of incident CHD. Study participants are enrolled in the CoLaus|PsyCoLaus study, a longitudinal, population-based cohort with baseline assessments from 2003 through 2008 and follow-up visits every 5 years. We analyzed a subgroup of 3459 individuals with available genome-wide genotyping data and immunoglobulin G levels for 22 persistent or frequently recurring pathogens. All reported CHD events were evaluated by a panel of specialists. We identified independent associations with incident CHD using univariable and multivariable stepwise Cox proportional hazards regression analyses. Of the 3459 study participants, 210 (6.07%) had at least one CHD event during the 12 years of follow-up. Multivariable stepwise Cox regression analysis, adjusted for known cardiovascular risk factors, socioeconomic status, and statin intake, revealed that high polygenic risk (hazard ratio [HR] 1.31, 95% CI 1.10–1.56, p=2.64 × 10−3) and infection with Fusobacterium nucleatum (HR 1.63, 95% CI 1.08–2.45, p=1.99 × 10−2) were independently associated with incident CHD. In a prospective, population-based cohort, high polygenic risk and infection with F. nucleatum have a small, yet independent impact on CHD risk.

Existing approaches to sharing of distributed medical data either provide only limited protection of patients’ privacy or sacrifice the accuracy of results. Here, the authors propose a federated analytics system, based on multiparty homomorphic encryption (MHE), to overcome these issues.

Sofosbuvir is a common therapy in hepatitis C virus infection, which targets the NS5B polymerase. Here, Smith et al. analyze the association between whole genome HCV polymorphisms and sofosbuvir treatment failure and identify three common polymorphisms present in non-targeted NS2 and NS3 proteins associated with reduced treatment response.

Summary: Genome-wide association studies (GWASs) have demonstrated that most common diseases have a strong genetic component from many genetic variants each with a small effect size. GWAS summary statistics have allowed the construction of polygenic scores (PGSs) estimating part of the individual risk for common diseases. Here, we propose to improve PGS-based risk estimation by incorporating genetic ancestry derived from genome-wide genotyping data. Our method involves three cohorts: a base (or discovery) for association studies, a target for phenotype/risk prediction, and a map for ancestry mapping; successively, (1) it generates for each individual in the base and target cohorts a set of principal components based on the map cohort—called mapped PCs, (2) it associates in the base cohort the phenotype with the mapped-PCs, and (3) it uses the mapped PCs in the target cohort to generate a phenotypic predictor called the ancestry score. We evaluated the ancestry score by comparing a predictive model using a PGS with one combining a PGS and an ancestry score. First, we performed simulations and found that the ancestry score has a greater impact on traits that correlate with ancestry-specific variants. Second, we showed, using UK Biobank data, that the ancestry score improves genetic prediction for our nine phenotypes to very different degrees. Third, we performed simulations and found that the more heterogeneous the base and target cohorts, the more beneficial the ancestry score is. Finally, we validated our approach under realistic conditions with UK Biobank as the base cohort and Swiss individuals from the CoLaus|PsyCoLaus study as the target cohort.

Abstract Epstein–Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count

Viral infections are the leading cause of childhood acute febrile illnesses motivating consultation in sub-Saharan Africa. The majority of causal viruses are never identified in low-resource clinical settings as such testing is either not part of routine screening or available diagnostic tools have limited ability to detect new/unexpected viral variants. An in-depth exploration of the blood virome is therefore necessary to clarify the potential viral origin of fever in children. Metagenomic next-generation sequencing is a powerful tool for such broad investigations, allowing the detection of RNA and DNA viral genomes. Here, we describe the blood virome of 816 febrile children (

T cell activation markers (TAM) expressed by antigen-specific T cells constitute promising candidates to attest the presence of an active infection by Mycobacterium tuberculosis (Mtb). Reciprocally, their modulation may be used to assess antibiotic treatment efficacy and eventually attest disease resolution. We hypothesized that the phenotype of Mtb-specific T cells may be quantitatively impacted by the load of bacteria present in a patient. We recruited 105 Tanzanian adult tuberculosis (TB) patients and obtained blood before and after 5 months of antibiotic treatment. We studied relationships between patients' clinical characteristics of disease severity and microbiological as well as molecular proxies of bacterial load in sputum at the time of diagnosis. Besides, we measured by flow cytometry the expression of CD38 or CD27 on CD4+ T cells producing interferon gamma (IFN-γ) and/or tumor necrosis factor alpha (TNF-α) in response to a synthetic peptide pool covering the sequences of Mtb antigens ESAT-6, CFP-10, and TB10.4. Reflecting the difficulty to extrapolate bacterial burden from a single end-point read-out, we observed statistically significant but weak correlations between Xpert MTB/RIF, molecular bacterial load assay and time to culture positivity. Unlike CD27, the resolution of CD38 expression by antigen-specific T cells was observed readily following 5 months of antibiotic therapy. However, the intensity of CD38-TAM signals measured at diagnosis did not significantly correlate with Mtb 16S RNA or rpoB DNA detected in patients' sputa. Altogether, our data support CD38-TAM as an accurate marker of infection resolution independently of sputum bacterial load.

The HIV reservoir is a major hurdle for a cure of HIV, but the factors determining its size and dynamics remain unclear. Here the authors show in a large cohort of 610 HIV-1 infected individuals, who are on suppressive ART for a median of 5.4 years, that viral genetic factors contribute substantially to the HIV-1 reservoir size.

Genome-wide association studies rely on the statistical inference of untyped variants, called imputation, to increase the coverage of genotyping arrays. However, the results are often suboptimal in populations underrepresented in existing reference panels and array designs, since the selected single nucleotide polymorphisms (SNPs) may fail to capture population-specific haplotype structures, hence the full extent of common genetic variation. Here, we propose to sequence the full genomes of a small subset of an underrepresented study cohort to inform the selection of population-specific add-on tag SNPs and to generate an internal population-specific imputation reference panel, such that the remaining array-genotyped cohort could be more accurately imputed. Using a Tanzania-based cohort as a proof-of-concept, we demonstrate the validity of our approach by showing improvements in imputation accuracy after the addition of our designed add-on tags to the base H3Africa array.

Background: Considering the remaining threat of drug-resistantmutations (DRMs) to antiretroviral treatment (ART) efficacy, we investigated how the selective pressure of human leukocyte antigen (HLA)-restricted cytotoxic T lymphocytes drives certain DRMs’ emergence and retention. Methods: We systematically screened DRM:HLA class I allele combinations in 3997 ART-naïve Swiss HIV Cohort Study (SHCS) patients. For each pair, a logistic regression model preliminarily tested for an association with the DRM as the outcome. The three HLA:DRM pairs remaining after multiple testing adjustment were analyzed in three ways: cross-sectional logistic regression models to determine any HLA/infection time interaction, survival analyses to examine if HLA type correlated with developing specific DRMs, and via NetMHCpan to find epitope binding evidence of immune escape. Results: Only one pair, RT-E138:HLA-B18, exhibited a significant interaction between infection duration and HLA. The survival analyses predicted two pairs with an increased hazard of developing DRMs: RT-E138:HLA-B18 and RT-V179:HLA-B35. RT-E138:HLA-B18 exhibited the greatest significance in both analyses (interaction term odds ratio [OR] 1.169 [95% confidence interval (CI) 1.075–1.273]; p-value

Abstract Background The gut microbiome is an important determinant of human health. Its composition has been shown to be influenced by multiple environmental factors and likely by host genetic variation. In the framework of the Milieu Intérieur Consortium, a total of 1000 healthy individuals of western European ancestry, with a 1:1 sex ratio and evenly stratified across five decades of life (age 20–69), were recruited. We generated 16S ribosomal RNA profiles from stool samples for 858 participants. We investigated genetic and non-genetic factors that contribute to individual differences in fecal microbiome composition. Results Among 110 demographic, clinical, and environmental factors, 11 were identified as significantly correlated with α-diversity, ß-diversity, or abundance of specific microbial communities in multivariable models. Age and blood alanine aminotransferase levels showed the strongest associations with microbiome diversity. In total, all non-genetic factors explained 16.4% of the variance. We then searched for associations between > 5 million single nucleotide polymorphisms and the same indicators of fecal microbiome diversity, including the significant non-genetic factors as covariates. No genome-wide significant associations were identified after correction for multiple testing. A small fraction of previously reported associations between human genetic variants and specific taxa could be replicated in our cohort, while no replication was observed for any of the diversity metrics. Conclusion In a well-characterized cohort of healthy individuals, we identified several non-genetic variables associated with fecal microbiome diversity. In contrast, host genetics only had a negligible influence. Demographic and environmental factors are thus the main contributors to fecal microbiome composition in healthy individuals. Trial registration ClinicalTrials.gov identifier NCT01699893

The LabEx Milieu Interieur (MI) project is a clinical study centered on the detailed characterization of the baseline and induced immune responses in blood samples from 1,000 healthy donors. Analyses of these samples has lay ground for seminal studies on the genetic and environmental determinants of immunologic variance in a healthy cohort population. In the current study we developed in vitro methods enabling standardized quantification of MI-cohort-derived primary fibroblasts responses. Our results show that in vitro human donor cohort fibroblast responses to stimulation by different MAMPs analogs allows to characterize individual donor immune-phenotype variability. The results provide proof-of-concept foundation to a new experimental framework for such studies. A bio-bank of primary fibroblast lines was generated from 323 out of 1,000 healthy individuals selected from the MI-study cohort. To study inter-donor variability of innate immune response in primary human dermal fibroblasts we chose to measure the TLR3 and TLR4 response pathways, both receptors being expressed and previously studied in fibroblasts. We established high-throughput automation compatible methods for standardized primary fibroblast cell activation, using purified MAMPS analogs, poly I:C and LPS that stimulate TLR3 and TLR4 pathways respectively. These results were in turn compared with a stimulation method using infection by HSV-1 virus. Our “Add-only” protocol minimizes high-throughput automation system variability facilitating whole process automation from cell plating through stimulation to recovery of cell supernatants, and fluorescent labeling. Images were acquired automatically by high-throughput acquisition on an automated high-content imaging microscope. Under these methodological conditions standardized image acquisition provided for quantification of cellular responses allowing biological variability to be measured with low system noise and high biological signal fidelity. Optimal for automated analysis of immuno-phenotype of primary human cell responses our method and experimental framework as reported here is highly compatible to high-throughput screening protocols like those necessary for chemo-genomic screening. In context of primary fibroblasts derived from donors enrolled to the MI-clinical-study our results open the way to assert the utility of studying immune-phenotype characteristics relevant to a human clinical cohort.

Abstract The identification of individuals with null alleles enables studying how the loss of gene function affects infection. We previously described a non‐functional variant in SIGLEC1, which encodes the myeloid‐cell receptor Siglec‐1/CD169 implicated in HIV‐1 cell‐to‐cell transmission. Here we report a significant association between the SIGLEC1 null variant and extrapulmonary dissemination of Mycobacterium tuberculosis (Mtb) in two clinical cohorts comprising 6,256 individuals. Local spread of bacteria within the lung is apparent in Mtb‐infected Siglec‐1 knockout mice which, despite having similar bacterial load, developed more extensive lesions compared to wild type mice. We find that Siglec‐1 is necessary to induce antigen presentation through extracellular vesicle uptake. We postulate that lack of Siglec‐1 delays the onset of protective immunity against Mtb by limiting antigen exchange via extracellular vesicles, allowing for an early local spread of mycobacteria that increases the risk for extrapulmonary dissemination.

Human immunodeficiency virus (HIV) infection is associated with an increased risk of non-Hodgkin lymphoma (NHL). Even in the era of suppressive antiretroviral treatment, HIV-infected individuals remain at higher risk of developing NHL compared to the general population. To identify potential genetic risk loci, we performed case-control genome-wide association studies and a meta-analysis across three cohorts of HIV+ patients of European ancestry, including a total of 278 cases and 1924 matched controls. We observed a significant association with NHL susceptibility in the C-X-C motif chemokine ligand 12 (CXCL12) region on chromosome 10. A fine mapping analysis identified rs7919208 as the most likely causal variant (P = 4.77e-11), with the G>A polymorphism creating a new transcription factor binding site for BATF and JUND. These results suggest a modulatory role of CXCL12 regulation in the increased susceptibility to NHL observed in the HIV-infected population.

Following the rapid dissemination of COVID-19 cases in Switzerland, large-scale non-pharmaceutical interventions (NPIs) were implemented by the cantons and the federal government between 28 February and 20 March 2020. Estimates of the impact of these interventions on SARS-CoV-2 transmission are critical for decision making in this and future outbreaks. We here aim to assess the impact of these NPIs on disease transmission by estimating changes in the basic reproduction number (R0) at national and cantonal levels in relation to the timing of these NPIs. We estimated the time-varying R0 nationally and in eleven cantons by fitting a stochastic transmission model explicitly simulating within-hospital dynamics. We used individual-level data from more than 1000 hospitalised patients in Switzerland and public daily reports of hospitalisations and deaths. We estimated the national R0 to be 2.8 (95% confidence interval 2.1–3.8) at the beginning of the epidemic. Starting from around 7 March, we found a strong reduction in time-varying R0 with a 86% median decrease (95% quantile range [QR] 79–90%) to a value of 0.40 (95% QR 0.3–0.58) in the period of 29 March to 5 April. At the cantonal level, R0 decreased over the course of the epidemic between 53% and 92%. Reductions in time-varying R0 were synchronous with changes in mobility patterns as estimated through smartphone activity, which started before the official implementation of NPIs. We inferred that most of the reduction of transmission is attributable to behavioural changes as opposed to natural immunity, the latter accounting for only about 4% of the total reduction in effective transmission. As Switzerland considers relaxing some of the restrictions of social mixing, current estimates of time-varying R0 well below one are promising. However, as of 24 April 2020, at least 96% (95% QR 95.7–96.4%) of the Swiss population remains susceptible to SARS-CoV-2. These results warrant a cautious relaxation of social distance practices and close monitoring of changes in both the basic and effective reproduction numbers.

Switzerland is among the countries with the highest number of coronavirus disease-2019 (COVID-19) cases per capita in the world. There are likely many people with undetected SARS-CoV-2 infection because testing efforts are currently not detecting all infected people, including some with clinical disease compatible with COVID-19. Testing on its own will not stop the spread of SARS-CoV-2. Testing is part of a strategy. The World Health Organization recommends a combination of measures: rapid diagnosis and immediate isolation of cases, rigorous tracking and precautionary self-isolation of close contacts. In this article, we explain why the testing strategy in Switzerland should be strengthened urgently, as a core component of a combination approach to control COVID-19.

Abstract Background Humoral immune responses to infectious agents or vaccination vary substantially among individuals, and many of the factors responsible for this variability remain to be defined. Current evidence suggests that human genetic variation influences (i) serum immunoglobulin levels, (ii) seroconversion rates, and (iii) intensity of antigen-specific immune responses. Here, we evaluated the impact of intrinsic (age and sex), environmental, and genetic factors on the variability of humoral response to common pathogens and vaccines. Methods We characterized the serological response to 15 antigens from common human pathogens or vaccines, in an age- and sex-stratified cohort of 1000 healthy individuals (Milieu Intérieur cohort). Using clinical-grade serological assays, we measured total IgA, IgE, IgG, and IgM levels, as well as qualitative (serostatus) and quantitative IgG responses to cytomegalovirus, Epstein-Barr virus, herpes simplex virus 1 and 2, varicella zoster virus, Helicobacter pylori, Toxoplasma gondii, influenza A virus, measles, mumps, rubella, and hepatitis B virus. Following genome-wide genotyping of single nucleotide polymorphisms and imputation, we examined associations between ~ 5 million genetic variants and antibody responses using single marker and gene burden tests. Results We identified age and sex as important determinants of humoral immunity, with older individuals and women having higher rates of seropositivity for most antigens. Genome-wide association studies revealed significant associations between variants in the human leukocyte antigen (HLA) class II region on chromosome 6 and anti-EBV and anti-rubella IgG levels. We used HLA imputation to fine map these associations to amino acid variants in the peptide-binding groove of HLA-DRβ1 and HLA-DPβ1, respectively. We also observed significant associations for total IgA levels with two loci on chromosome 2 and with specific KIR-HLA combinations. Conclusions Using extensive serological testing and genome-wide association analyses in a well-characterized cohort of healthy individuals, we demonstrated that age, sex, and specific human genetic variants contribute to inter-individual variability in humoral immunity. By highlighting genes and pathways implicated in the normal antibody response to frequently encountered antigens, these findings provide a basis to better understand disease pathogenesis. Trials registration ClinicalTrials.gov, NCT01699893

Genetic polymorphism in the interferon lambda (IFN-λ) region is associated with spontaneous clearance of hepatitis C virus (HCV) infection and response to interferon-based treatment. Here, we evaluate associations between IFN-λ polymorphism and HCV variation in 8729 patients (Europeans 77%, Asians 13%, Africans 8%) infected with various viral genotypes, predominantly 1a (41%), 1b (22%) and 3a (21%). We searched for associations between rs12979860 genotype and variants in the NS3, NS4A, NS5A and NS5B HCV proteins. We report multiple associations in all tested proteins, including in the interferon-sensitivity determining region of NS5A. We also assessed the combined impact of human and HCV variation on pretreatment viral load and report amino acids associated with both IFN-λ polymorphism and HCV load across multiple viral genotypes. By demonstrating that IFN-λ variation leaves a large footprint on the viral proteome, we provide evidence of pervasive viral adaptation to innate immune pressure during chronic HCV infection.

Fc receptor (FcR) genes collectively have copy number and allelic polymorphisms that have been implicated in multiple inflammatory and autoimmune diseases. This variation might also be involved in etiology of infectious diseases. The protective role of Fc-mediated antibody-function in HIV-1 immunity has led to the investigation of specific polymorphisms in FcR genes on acquisition, disease progression, and vaccine efficacy in natural history cohorts. The purpose of this review is not only to explore these known HIV-1 host genetic associations, but also to re-evaluate them in the context of genome-wide data. In the current era of effective anti-retroviral therapy, the potential impact of such variation on post-treatment cohorts cannot go unheeded and is discussed here in the light of current findings. Specific polymorphisms associating with HIV-1 pathogenesis have previously been genotyped by assays that captured only the single-nucleotide polymorphism (SNP) of interest without relative information of neighboring variants. With recent technological advances, variation within these genes can now be characterized using next-generation sequencing, allowing precise annotation of the whole chromosomal region. We herein also discuss updates in the annotation of common FcR variants that have been previously associated with HIV-1 pathogenesis.

Many viral pathogens target innate sensing cascades and/or cellular transcription factors to suppress antiviral immune responses. Here, we show that the accessory viral protein U (Vpu) of HIV-1 exerts broad immunosuppressive effects by inhibiting activation of the transcription factor NF-κB. Global transcriptional profiling of infected CD4 +T cells revealed that vpu-deficient HIV-1 strains induce substantially stronger immune responses than the respective wild type viruses. Gene set enrichment analyses and cytokine arrays showed that Vpu suppresses the expression of NF-κB targets including interferons and restriction factors. Mutational analyses demonstrated that this immunosuppressive activity of Vpu is independent of its ability to counteract the restriction factor and innate sensor tetherin. However, Vpu-mediated inhibition of immune activation required an arginine residue in the cytoplasmic domain that is critical for blocking NF-κB signaling downstream of tetherin. In summary, our findings demonstrate that HIV-1 Vpu potently suppresses NF-κB-elicited antiviral immune responses at the transcriptional level.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Systems approaches for the study of immune signaling pathways have been traditionally based on purified cells or cultured lines. However, in vivo responses involve the coordinated action of multiple cell types, which interact to establish an inflammatory microenvironment. We employed standardized whole-blood stimulation systems to test the hypothesis that responses to Toll-like receptor ligands or whole microbes can be defined by the transcriptional signatures of key cytokines. We found 44 genes, identified using Support Vector Machine learning, that captured the diversity of complex innate immune responses with improved segregation between distinct stimuli. Furthermore, we used donor variability to identify shared inter-cellular pathways and trace cytokine loops involved in gene expression. This provides strategies for dimension reduction of large datasets and deconvolution of innate immune responses applicable for characterizing immunomodulatory molecules. Moreover, we provide an interactive R-Shiny application with healthy donor reference values for induced inflammatory genes.

Binding of virus, HIV-1, to cellular protein Siglec-1 is important for infection of immune cells. Here the authors show that a natural mutation leading to production of truncated Siglec-1 reduces HIV binding and infectivity transfer in vitro, but does not substantially affect infection or AIDS outcome in patients.

The huge genetic diversity of circulating viruses is a challenge for diagnostic assays for emerging or rare viral diseases. High-throughput technology offers a new opportunity to explore the global virome of patients without preconception about the culpable pathogens. It requires a solid reference dataset to be accurate. Virosaurus has been designed to offer a non-biased, automatized and annotated database for clinical metagenomics studies and diagnosis. Raw viral sequences have been extracted from GenBank, and cleaned up to remove potentially erroneous sequences. Complete sequences have been identified for all genera infecting vertebrates, plants and other eukaryotes (insect, fungus, etc.). To facilitate the analysis of clinically relevant viruses, we have annotated all sequences with official and common virus names, acronym, genotypes, and genomic features (linear, circular, DNA, RNA, etc.). Sequences have been clustered to remove redundancy at 90% or 98% identity. The analysis of clustering results reveals the state of the virus genetic landscape knowledge. Because herpes and poxviruses were under-represented in complete genomes considering their potential diversity in nature, we used genes instead of complete genomes for those in Virosaurus.

Studies of host genetic determinants of pathogen sequence variations can identify sites of genomic conflicts, by highlighting variants that are implicated in immune response on the host side and adaptive escape on the pathogen side. However, systematic genetic differences in host and pathogen populations can lead to inflated type I (false positive) and type II (false negative) error rates in genome-wide association analyses. Here, we demonstrate through a simulation that correcting for both host and pathogen stratification reduces spurious signals and increases power to detect real associations in a variety of tested scenarios. We confirm the validity of the simulations by showing comparable results in an analysis of paired human and HIV genomes.

Two cases of neonatal group B streptococcal (GBS) infection occurring in a large consanguineous kindred are described. The observation suggests that susceptibility to isolated GBS disease in this family may have been the result of an inborn error of immunity inherited as a Mendelian autosomal recessive trait.

[This corrects the article DOI: 10.1371/journal.pbio.2001855.].

HIV-1 set-point viral load-the approximately stable value of viraemia in the first years of chronic infection-is a strong predictor of clinical outcome and is highly variable across infected individuals. To better understand HIV-1 pathogenesis and the evolution of the viral population, we must quantify the heritability of set-point viral load, which is the fraction of variation in this phenotype attributable to viral genetic variation. However, current estimates of heritability vary widely, from 6% to 59%. Here we used a dataset of 2,028 seroconverters infected between 1985 and 2013 from 5 European countries (Belgium, Switzerland, France, the Netherlands and the United Kingdom) and estimated the heritability of set-point viral load at 31% (CI 15%-43%). Specifically, heritability was measured using models of character evolution describing how viral load evolves on the phylogeny of whole-genome viral sequences. In contrast to previous studies, (i) we measured viral loads using standardized assays on a sample collected in a strict time window of 6 to 24 months after infection, from which the viral genome was also sequenced; (ii) we compared 2 models of character evolution, the classical "Brownian motion" model and another model ("Ornstein-Uhlenbeck") that includes stabilising selection on viral load; (iii) we controlled for covariates, including age and sex, which may inflate estimates of heritability; and (iv) we developed a goodness of fit test based on the correlation of viral loads in cherries of the phylogenetic tree, showing that both models of character evolution fit the data well. An overall heritability of 31% (CI 15%-43%) is consistent with other studies based on regression of viral load in donor-recipient pairs. Thus, about a third of variation in HIV-1 virulence is attributable to viral genetic variation.

We evaluated the fraction of variation in HIV-1 set point viral load attributable to viral or human genetic factors by using joint host/pathogen genetic data from 541 HIV infected individuals. We show that viral genetic diversity explains 29% of the variation in viral load while host factors explain 8.4%. Using a joint model including both host and viral effects, we estimate a total of 30% heritability, indicating that most of the host effects are reflected in viral sequence variation.

One out of three pediatric sepsis deaths in high income countries occur in previously healthy children. Primary immunodeficiencies have been postulated to underlie fulminant sepsis, but this concept remains to be confirmed in clinical practice. Pseudomonas aeruginosa (P. aeruginosa) is a common bacterium mostly associated with healthcare-related infections in immunocompromised individuals. However, in rare cases, it can cause sepsis in previously healthy children. We used exome sequencing and bioinformatic analysis to systematically search for genetic factors underpinning severe P. aeruginosa infection in the pediatric population. We collected blood samples from 11 previously healthy children, with no family history of immunodeficiency, who presented with severe sepsis due to community-acquired P. aeruginosa bacteremia. Genomic DNA was extracted from blood or tissue samples obtained intravitam or postmortem. We obtained high-coverage exome sequencing data and searched for rare loss-of-function variants. After rigorous filtrations, 12 potentially causal variants were identified. 2 out of 8 (25%) fatal cases were found to carry novel pathogenic variants in primary immunodeficiency genes, including BTK and DNMT3B. This study demonstrates that exome sequencing allows to identify rare, deleterious human genetic variants responsible for fulminant sepsis in apparently healthy children. Diagnosing primary immunodeficiencies in such patients is of high relevance to survivors and affected families. We propose that unusually severe and fatal sepsis cases in previously healthy children should be considered for exome/genome sequencing to search for underlying primary immunodeficiencies.