Genome-wide association and HLA region fine-mapping studies identify susceptibility loci for multiple common infections

Genome-Wide Determinants of Cellular Immune Responses to Mumps Vaccine

Background

We have previously described genetic polymorphisms in candidate genes that are associated with inter-individual variations in antibody responses to mumps vaccination. To expand upon our previous work, we performed a genome-wide association study (GWAS) to discover host genetic variants associated with mumps vaccine-induced cellular immune responses.

Methods

We performed a GWAS of mumps-specific immune response outcomes (11 secreted cytokines/chemokines) in a cohort of 1,406 subjects.

Results

Among the 11 cytokine/chemokines we studied, four (IFN-γ, IL-2, IL-1β, and TNFα) demonstrated GWAS signals reaching genome-wide significance (p<5 x 10 -8 ). A genomic region (encoding Sialic acid-binding immunoglobulin-type lectins/SIGLEC) located on chromosome 19q13 (p<5×10 -8 ) was associated with both IL-1β and TNFα responses. The SIGLEC5/SIGLEC14 region contained 11 statistically significant single nucleotide polymorphisms (SNPs), including the intronic SIGLEC5 rs872629 (p=1.3E-11) and rs1106476 (p=1.32E-11) whose alternate alleles were significantly associated with decreased levels of mumps-specific IL-1β (rs872629, p=1.77E-09; rs1106476, p=1.78E-09) and TNFα (rs872629, p=1.3E-11; rs1106476, p=1.32E-11) production.

Conclusions

Our results suggest that SNPs in the SIGLEC5/SIGLEC14 genes play a role in cellular and inflammatory immune responses to mumps vaccination. These findings motivate further research into the functional roles of SIGLEC genes in the regulation of mumps vaccine-induced immunity.

Fulminant pulmonary tuberculosis in a previously healthy young woman from the Marshall Islands: Potential risk factors

A 19-year-old woman originally from the Republic of the Marshall Islands presented with diffuse pneumonia and acute hypoxemic respiratory failure. She dies one month into her hospitalization but the diagnosis of pulmonary tuberculosis (TB) was not made until one day before her demise. A contact investigation screened a total of 155 persons with 36 (23%) found to have latent TB infection and seven (4.5%) with active pulmonary TB. This unfortunate case provided the opportunity to analyze the epidemiology of TB in the state of Washington in the context of those who emigrated from the Marshall Islands. The development of fulminant pulmonary TB in this previously healthy young woman also provides a segue to discuss potential risk factors for TB in the index case that include: (i) foreign-born in a TB-endemic country; (ii) race and genetic factors; (iii) age; (iv) body habitus; (v) pregnancy; and (vi) use of glucocorticoids.

Genome-Wide Meta-Analysis Identifies Multiple Novel Rare Variants to Predict Common Human Infectious Diseases Risk

Infectious diseases still threaten global human health, and host genetic factors have been indicated as determining risk factors for observed variations in disease susceptibility, severity, and outcome. We performed a genome-wide meta-analysis on 4624 subjects from the 10,001 Dalmatians cohort, with 14 infection-related traits. Despite a rather small number of cases in some instances, we detected 29 infection-related genetic associations, mostly belonging to rare variants. Notably, the list included the genes CD28, INPP5D, ITPKB, MACROD2, and RSF1, all of which have known roles in the immune response. Expanding our knowledge on rare variants could contribute to the development of genetic panels that could assist in predicting an individual's life-long susceptibility to major infectious diseases. In addition, longitudinal biobanks are an interesting source of information for identifying the host genetic variants involved in infectious disease susceptibility and severity. Since infectious diseases continue to act as a selective pressure on our genomes, there is a constant need for a large consortium of biobanks with access to genetic and environmental data to further elucidate the complex mechanisms behind host-pathogen interactions and infectious disease susceptibility.

Genome‐wide association study of tuberculosis in the western Chinese Han and Tibetan population

Tuberculosis (TB) remains a serious global public health threat. Accumulated evidence has demonstrated that human susceptibility to TB has a strong genetic basis. And different susceptibility single nucleotide polymorphisms (SNP) have been reported in different studies. To gain greater insight into the host susceptibility to TB, we perform a two‐stage genome‐wide association study to identify the susceptible loci of TB. In the discovery stage, 3116 (1532 TB patients and 1584 healthy controls) and 439 (211 TB patients and 228 healthy controls) individuals were genome‐wide genotyped from a western Chinese Han and Tibetan population, respectively. Based on the additive genetic model, we discovered 14 and three independent loci that had potential associations with TB susceptibility in the Chinese Han and Tibetan populations, respectively (p < 1 × 10⁻⁵). Furthermore, we conducted an imputation‐based meta‐analysis on another two East Asia cohorts to replicate our findings. We identified one independent locus harbored by the human leukocyte antigen (HLA) class II genes that was genome‐wide significantly associated with TB (lead SNP rs111875628 with a p‐value of 2.20 × 10⁻⁹). Our findings suggest a novel mechanism of the interaction with the HLA class II genes and reinforce the importance of the HLA class II alleles in response to TB.

Ethical Challenges Associated with Pathogen and Host Genetics in Infectious Disease

The Covid-19 pandemic has demonstrated the potential of genomic technologies for the detection and surveillance of infectious diseases. Pathogen genomics is likely to play a major role in the future of research and clinical implementation of genomic technologies. However, unlike human genetics, the specific ethical and social challenges associated with the implementation of infectious disease genomics has received comparatively little attention. In this paper, we contribute to this literature, focusing on the potential consequences for individuals and communities of the use of these technologies. We concentrate on areas of challenges related to privacy, stigma, discrimination and the return of results in the cases of the surveillance of known pathogens, metagenomics and host genomics.

Altered IL-6 signalling and risk of tuberculosis disease: a meta-analysis and Mendelian randomisation study

IL-6 responses are ubiquitous in Mycobacterium tuberculosis (Mtb) infections, but their role in determining human tuberculosis (TB) disease risk is unknown. We used single nucleotide polymorphisms (SNPs) in and near the IL-6 receptor (IL6R) gene, focusing on the non-synonymous variant, rs2228145, associated with reduced classical IL-6 signalling, to assess the effect of altered IL-6 activity on TB disease risk. We identified 16 genome wide association studies (GWAS) of TB disease collating 17,982 cases of TB disease and 972,389 controls across 4 continents. Meta-analyses and Mendelian randomisation analyses revealed that reduced classical IL-6 signalling was associated with lower odds of TB disease, a finding replicated using multiple, independent SNP instruments and 2 separate exposure variables. Our findings establish a causal relationship between IL-6 signalling and the outcome of Mtb infection, suggesting IL-6 antagonists do not increase the risk of TB disease and should be investigated as adjuncts in treatment.

Relevance of HLA-DP/DQ and INF-λ4 Polymorphisms to COVID-19 Outcomes

Background: Single nucleotide polymorphisms provide information on individuals' potential reactions to environmental factors, infections, diseases, as well as various therapies. A study on SNPs that influence SARS-CoV-2 susceptibility and severity may provide a predictive tool for COVID-19 outcomes and improve the customized coronavirus treatment. Aim: To evaluate the role of human leukocyte antigens DP/DQ and IFNλ4 polymorphisms on COVID-19 outcomes among Egyptian patients. Participants and Methods: The study involved 80 patients with severe COVID-19, 80 patients with mild COVID-19, and 80 non-infected healthy volunteers. Genotyping and allelic discrimination of HLA-DPrs3077 (G/A), HLA-DQrs7453920 (A/G), and IFNλ4 rs73555604 (C/T) SNPs were performed using real-time PCR. Results: Ages were 47.9 ± 8, 44.1 ± 12.1, and 45.8 ± 10 years in severe, mild and non-infected persons. There was a statistically significant association between severe COVID-19 and male gender (p = 0.002). A statistically significant increase in the frequency of HLA-DPrs3077G, HLA-DQrs7453920A, and IFNλ4rs73555604C alleles among severe COVID-19 patients when compared with other groups (p < 0.001). Coexistence of these alleles in the same individual increases the susceptibility to severe COVID-19 by many folds (p < 0.001). Univariate and multivariate logistic regression analysis for the studied parameters showed that old age, male gender, non-vaccination, HLA-DQ rs7453920AG+AA, HLA-DPrs3077GA+GG, and IFNλ4rs73555604CT+CC genotypes are independent risk factors for severe COVID-19 among Egyptian patients. Conclusion: HLA-DQ rs7453920A, HLA-DPrs3077G, and IFNλ4rs73555604C alleles could be used as markers of COVID-19 severity.

Inflammatory and infectious upper respiratory diseases associate with 41 genomic loci and type 2 inflammation

Inflammatory and infectious upper respiratory diseases (ICD-10: J30-J39), such as diseases of the sinonasal tract, pharynx and larynx, are growing health problems yet their genomic similarity is not known. We analyze genome-wide association to eight upper respiratory diseases (61,195 cases) among 260,405 FinnGen participants, meta-analyzing diseases in four groups based on an underlying genetic correlation structure. Aiming to understand which genetic loci contribute to susceptibility to upper respiratory diseases in general and its subtypes, we detect 41 independent genome-wide significant loci, distinguishing impact on sinonasal or pharyngeal diseases, or both. Fine-mapping implicated non-synonymous variants in nine genes, including three linked to immune-related diseases. Phenome-wide analysis implicated asthma and atopic dermatitis at sinonasal disease loci, and inflammatory bowel diseases and other immune-mediated disorders at pharyngeal disease loci. Upper respiratory diseases also genetically correlated with autoimmune diseases such as rheumatoid arthritis, autoimmune hypothyroidism, and psoriasis. Finally, we associated separate gene pathways in sinonasal and pharyngeal diseases that both contribute to type 2 immunological reaction. We show shared heritability among upper respiratory diseases that extends to several immune-mediated diseases with diverse mechanisms, such as type 2 high inflammation.

Genetic adaptation to pathogens and increased risk of inflammatory disorders in post-Neolithic Europe

Ancient genomics can directly detect human genetic adaptation to environmental cues. However, it remains unclear how pathogens have exerted selective pressures on human genome diversity across different epochs and affected present-day inflammatory disease risk. Here, we use an ancestry-aware approximate Bayesian computation framework to estimate the nature, strength, and time of onset of selection acting on 2,879 ancient and modern European genomes from the last 10,000 years. We found that the bulk of genetic adaptation occurred after the start of the Bronze Age, <4,500 years ago, and was enriched in genes relating to host-pathogen interactions. Furthermore, we detected directional selection acting on specific leukocytic lineages and experimentally demonstrated that the strongest negatively selected candidate variant in immunity genes, lipopolysaccharide-binding protein (LBP) D283G, is hypomorphic. Finally, our analyses suggest that the risk of inflammatory disorders has increased in post-Neolithic Europeans, possibly because of antagonistic pleiotropy following genetic adaptation to pathogens.

Haplotypic variants of COVID-19 related genes are associated with blood pressure and metabolites levels

The genetic association of coronavirus disease 2019 (COVID-19) with its complications has not been fully understood. This study aimed to identify variants and haplotypes of candidate genes implicated in COVID-19 related traits by combining the literature review and pathway analysis. To explore such genes, the protein-protein interactions and relevant pathways of COVID-19-associated genes were assessed. A number of variants on candidate genes were identified from Genome-wide association studies (GWASs) which were associated with COVID-19 related traits (p ˂ 10^-6 ). Haplotypic blocks were assessed using haplotypic structures among the 1000 Genomes Project (r²  ≥ 0.8, D' ≥ 0.8). Further functional analyses were performed on the selected variants. The results demonstrated that a group of variants in ACE and AGT genes were significantly correlated with COVID-19 related traits. Three haplotypes were identified to be involved in the blood metabolites levels and the development of blood pressure. Functional analyses revealed that most GWAS index variants were expression quantitative trait loci and had transcription factor binding sites, exonic splicing enhancers or silencer activities. Furthermore, the proxy haplotype variants, rs4316, rs4353, rs4359, and three variants, namely rs2493133, rs2478543, and rs5051, were associated with blood metabolite and systolic blood pressure, respectively. These variants exerted more regulatory effects compared with other GWAS variants. The present study indicates that the genetic variants and candidate haplotypes of COVID-19 related genes are associated with blood pressure and blood metabolites. However, further observational studies are warranted to confirm these results.

Chickenpox and multiple sclerosis: A Mendelian randomization study

Observational studies have suggested a suspected association between varicella-zoster virus (VZV) infection and multiple sclerosis (MS), but the connection has remained unclear. The aim of the present study is to evaluate the causal relationship between chickenpox which is caused by VZV infection and MS. We performed a two-sample Mendelian randomization analysis to investigate the association of chickenpox with MS using summary statistics from genome-wide association studies (GWAS). The GWAS summary statistics data for chickenpox was from the 23andMe cohort including 107 769 cases and 15 982 controls. A large summary of statistical data from the International Multiple Sclerosis Genetics Consortium (IMSGC) was used as the outcome GWAS data set, including 14 802 MS cases and 26 703 controls. We found evidence of a significant association between genetically predicted chickenpox and risk of MS (odds ratio [OR] = 35.27, 95% confidence interval [CI] = 22.97-54.17, p = 1.46E-59). Our findings provided evidence indicating a causal effect of chickenpox on MS. Further elucidations of this association and underlying mechanisms are needed for identifying feasible interventions to promote MS prevention.

Large registry-based analysis of genetic predisposition to tuberculosis identifies genetic risk factors at HLA

Tuberculosis is a significant public health concern resulting in the death of over 1 million individuals each year worldwide. While treatment options and vaccines exist, a substantial number of infections still remain untreated or are caused by treatment resistant strains. Therefore, it is important to identify mechanisms that contribute to risk and prognosis of tuberculosis as this may provide tools to understand disease mechanisms and provide novel treatment options for those with severe infection. Our goal was to identify genetic risk factors that contribute to the risk of tuberculosis and to understand biological mechanisms and causality behind the risk of tuberculosis. A total of 1895 individuals in the FinnGen study had International Classification of Diseases-based tuberculosis diagnosis. Genome-wide association study analysis identified genetic variants with statistically significant association with tuberculosis at the human leukocyte antigen (HLA) region (P < 5e-8). Fine mapping of the HLA association provided evidence for one protective haplotype tagged by HLA DQB1*05:01 (P = 1.82E-06, OR = 0.81 [CI 95% 0.74-0.88]), and predisposing alleles tagged by HLA DRB1*13:02 (P = 0.00011, OR = 1.35 [CI 95% 1.16-1.57]). Furthermore, genetic correlation analysis showed association with earlier reported risk factors including smoking (P < 0.05). Mendelian randomization supported smoking as a risk factor for tuberculosis (inverse-variance weighted P < 0.05, OR = 1.83 [CI 95% 1.15-2.93]) with no significant evidence of pleiotropy. Our findings indicate that specific HLA alleles associate with the risk of tuberculosis. In addition, lifestyle risk factors such as smoking contribute to the risk of developing tuberculosis.

The immunogenetics of COVID-19

The worldwide coronavirus disease 2019 pandemic was sparked by the severe acute respiratory syndrome caused by coronavirus 2 (SARS-CoV-2) that first surfaced in December 2019 (COVID-19). The effects of COVID-19 differ substantially not just between patients individually but also between populations with different ancestries. In humans, the human leukocyte antigen (HLA) system coordinates immune regulation. Since HLA molecules are a major component of antigen-presenting pathway, they play an important role in determining susceptibility to infectious disease. It is likely that differential susceptibility to SARS-CoV-2 infection and/or disease course in COVID-19 in different individuals could be influenced by the variations in the HLA genes which are associated with various immune responses to SARS-CoV-2. A growing number of studies have identified a connection between HLA variation and diverse COVID-19 outcomes. Here, we review research investigating the impact of HLA on individual responses to SARS-CoV-2 infection and/or progression, also discussing the significance of MHC-related immunological patterns and its use in vaccine design.

A multi-phenotype analysis reveals 19 susceptibility loci for basal cell carcinoma and 15 for squamous cell carcinoma

Basal cell carcinoma and squamous cell carcinoma are the most common skin cancers, and have genetic overlap with melanoma, pigmentation traits, autoimmune diseases, and blood biochemistry biomarkers. In this multi-trait genetic analysis of over 300,000 participants from Europe, Australia and the United States, we reveal 78 risk loci for basal cell carcinoma (19 previously unknown and replicated) and 69 for squamous cell carcinoma (15 previously unknown and replicated). The previously unknown risk loci are implicated in cancer development and progression (e.g. CDKL1), pigmentation (e.g. TPCN2), cardiometabolic (e.g. FADS2), and immune-regulatory pathways for innate immunity (e.g. IFIH1), and HIV-1 viral load modulation (e.g. CCR5). We also report an optimised polygenic risk score for effective risk stratification for keratinocyte cancer in the Canadian Longitudinal Study of Aging (794 cases and 18139 controls), which could facilitate skin cancer surveillance e.g. in high risk subpopulations such as transplantees.

The Awesome Power of Human Genetics of Infectious Disease

Since the identification of sickle cell trait as a heritable form of resistance to malaria, candidate gene studies, linkage analysis paired with sequencing, and genome-wide association (GWA) studies have revealed many examples of genetic resistance and susceptibility to infectious diseases. GWA studies enabled the identification of many common variants associated with small shifts in susceptibility to infectious diseases. This is exemplified by multiple loci associated with leprosy, malaria, HIV, tuberculosis, and coronavirus disease 2019 (COVID-19), which illuminate genetic architecture and implicate pathways underlying pathophysiology. Despite these successes, most of the heritability of infectious diseases remains to be explained. As the field advances, current limitations may be overcome by applying methodological innovations such as cellular GWA studies and phenome-wide association (PheWA) studies as well as by improving methodological rigor with more precise case definitions, deeper phenotyping, increased cohort diversity, and functional validation of candidate loci in the laboratory or human challenge studies.

Human leukocyte antigen system associations in Malassezia-related skin diseases

BACKGROUND

The human leukocyte antigen system (HLA) is divided into two classes involved in antigen presentation: class I presenting intracellular antigens and class II presenting extracellular antigens. While susceptibility to infections is correlated with the HLA system, data on associations between HLA genotypes and Malassezia-related skin diseases (MRSD) are lacking. Thus, the objective of this study was to investigate associations between HLA alleles and MRSD.

MATERIALS AND METHODS

Participants in The Danish Blood Donor Study (2010-2018) provided questionnaire data on life style, anthropometric measures, and registry data on filled prescriptions. Genotyping was done using Illumina Infinium Global Screening Array, and HLA alleles were imputed using the HIBAG algorithm. Cases and controls were defined using filled prescriptions on topical ketoconazole 2% as a proxy of MRSD. Logistic regressions assessed associations between HLA alleles and MRSD adjusted for confounders and Bonferroni corrected for multiple tests.

RESULTS

A total of 9455 participants were considered MRSD cases and 24,144 participants as controls. We identified four risk alleles B*57:01, OR 1.19 (95% CI: 1.09-1.31), C*01:02, OR 1.19 (95% CI: 1.08-1.32), C*06:02, OR 1.14 (95% CI: 1.08-1.22), and DRB1*01:01, OR 1.10 (95% CI: 1.04-1.17), and two protective alleles, DQB1*02:01, OR 0.89 (95% CI: 0.85-0.94), and DRB1*03:01, OR 0.89 (95% CI: 0.85-0.94).

CONCLUSION

Five novel associations between HLA alleles and MRSD were identified in our cohort, and one previous association was confirmed. Future studies should assess the correlation between Malassezia antigens and antigen-binding properties of the associated HLA alleles.

A genome-wide association study of survival in patients with sepsis

BACKGROUND

Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction and has a high mortality rate in adult intensive care units. Most genetic studies have identified gene variants associated with development and outcomes of sepsis focusing on biological candidates. We conducted the first genome-wide association study (GWAS) of 28-day survival in adult patients with sepsis.

METHODS

This study was conducted in two stages. The first stage was performed on 687 European sepsis patients from the GEN-SEP network and 7.5 million imputed variants. Association testing was conducted with Cox regression models, adjusting by sex, age, and the main principal components of genetic variation. A second stage focusing on the prioritized genetic variants was performed on 2,063 ICU sepsis patients (1362 European Americans and 701 African-Americans) from the MESSI study. A meta-analysis of results from the two stages was conducted and significance was established at p < 5.0 × 10-8. Whole-blood transcriptomic, functional annotations, and sensitivity analyses were evaluated on the identified genes and variants.

FINDINGS

We identified three independent low-frequency variants associated with reduced 28-day sepsis survival, including a missense variant in SAMD9 (hazard ratio [95% confidence interval] = 1.64 [1.37-6.78], p = 4.92 × 10-8). SAMD9 encodes a possible mediator of the inflammatory response to tissue injury.

INTERPRETATION

We performed the first GWAS of 28-day sepsis survival and identified novel variants associated with reduced survival. Larger sample size studies are needed to better assess the genetic effects in sepsis survival and to validate the findings.

Chromosome-scale haplotype-resolved pangenomics

Complete pangenomics is crucial for understanding genetic diversity and evolution across the tree of life. Chromosome-scale, haplotype-resolved pangenomics allows complex structural variations, long-range interactions, and associated functions to be discerned in species populations. We explore the need for high-resolution pangenomes, discuss computational strategies for their development, and describe applications in biodiversity and human health.

Maternal infection during pregnancy and likelihood of autism and intellectual disability in children in Sweden: a negative control and sibling comparison cohort study

BACKGROUND

Maternal infections during pregnancy are associated with intellectual disability and autism in exposed children. Whether these associations are causal, and therefore should be targets of preventive strategies, remains unknown. We aimed to investigate these associations, to determine whether there is a causal role of maternal infection during pregnancy for children's risk of autism and intellectual disability, by accounting for unmeasured familial factors.

METHODS

We used a register-based cohort study design, and included children living in Stockholm County, Sweden, who were born in 1987-2010. We excluded children not born in Sweden, adopted children, and children with unknown biological mothers or fathers. Maternal infections during pregnancy, defined by ICD-8, ICD-9, and ICD-10 codes, were identified in the National Patient Register and Medical Birth Register. Children were followed up from birth to an outcome or a censoring event (death, migration from Stockholm, age 18 years, or Dec 31, 2016, whichever occurred first). The primary outcomes were diagnosis of autism or diagnosis of intellectual disability. We did a survival analysis to examine the association between inpatient and outpatient specialised care for any infection during pregnancy and likelihood of autism or intellectual disability in the child. To address potential residual confounding, we also estimated the relationship between maternal infection in the year preceding pregnancy as a negative control exposure and conducted a matched sibling analysis of sibling pairs who were discordant for autism or intellectual disability.

FINDINGS

647 947 children living in Stockholm County were identified and, after excluding 97 980 children, we included 549 967 in the study (267 995 [48·7%] were female and 281 972 [51·3%] were male; mean age at censoring 13·5 years [SD 5·0; range <1 to 18]; 142 597 [25·9%] had a mother who was not born in Sweden). 445 (1·3%) of 34 013 children exposed to maternal infection during pregnancy were diagnosed with intellectual disability and 1123 (3·3%) with autism. 5087 (1·0%) of 515 954 unexposed children were diagnosed with intellectual disability and 13 035 (2·5%) with autism. Maternal infection during pregnancy was associated with autism (hazard ratio [HR] 1·16, 95% CI 1·09-1·23) and intellectual disability (1·37, 1·23-1·51) in exposed children compared with unexposed children. Maternal infection in the year before pregnancy (negative control exposure) was also associated with autism (HR 1·25, 95% CI 1·14-1·36), but was not associated with intellectual disability (1·09, 0·94-1·27). In sibling comparisons, the associations with maternal infection during pregnancy were attenuated for autism (HR 0·94, 95% CI 0·82-1·08; n=21 864), but not to the same extent for intellectual disability (1·15, 0·95-1·40; n=9275).

INTERPRETATION

Although infections in pregnant women are associated with both autism and intellectual disability in their children, the association with autism does not appear to reflect a causal relationship, but is more likely to be explained by factors shared between family members such as genetic variation or aspects of the shared environment. Thus, infection prevention is not expected to reduce autism incidence. For intellectual disability, unmeasured familial factors might not fully explain the observed associations, and a causal role of maternal infections cannot be excluded. Causal effects of specific but rare infections or infections not requiring health care contact cannot be excluded in either autism or intellectual disability.

FUNDING

Swedish Research Council, Stanley Medical Research Institute, and Autism Speaks.

TRANSLATION

For the Swedish translation of the abstract see Supplementary Materials section.

Genetic architecture of tuberculosis susceptibility: A comprehensive research synopsis, meta-analyses, and epidemiological evidence

To date, many studies have been conducted to investigate associations between variants and tuberculosis risk; however, the results have been inconclusive. Here, we systematically provide a summary of the understanding of the genetic architecture of tuberculosis susceptibility. We searched PubMed, Embase and Web of Science to identify genetic association studies of tuberculosis published through October 31, 2021. We conducted meta-analyses for the genetic association with tuberculosis risk. We graded levels of cumulative epidemiological evidence of significant associations with risk of tuberculosis and false-positive report probability tests. We performed functional annotations for these variants using data from the Encyclopedia of DNA Elements (ENCODE) Project and other databases. We identified 703 eligible articles comprising 298,074 cases and 879,593 controls through screening a total of 24,398 citations. Meta-analyses were conducted for 614 genetic variants in 469 genes or loci. We found 39 variants that were nominally significantly associated with tuberculosis risk. Cumulative epidemiological evidence for a significant association was graded strong for 9 variants in or near 9 genes. Among them, 5 variants were associated with tuberculosis risk in at least three main ethnicity (African, Asian and White) which together explained approximately 9.59% of the familial relative risk of tuberculosis. Data from ENCODE and other databases suggested that 8 of these 9 genetic variants with strong evidence might fall within putative functional regions. Our study summarizes the current literature on the genetic architecture of tuberculosis susceptibility and provides useful data for designing future studies to investigate the genetic association with tuberculosis risk.

Genome-wide association analyses of common infections in a large practice-based biobank

INTRODUCTION

Infectious diseases are common causes of morbidity and mortality worldwide. Susceptibility to infection is highly heritable; however, little has been done to identify the genetic determinants underlying common infectious diseases. One GWAS was performed using 23andMe information about self-reported infections; we set out to confirm previous loci and identify new ones using medically diagnosed infections.

METHODS

We used the electronic health record (EHR)-based biobank at Vanderbilt and diagnosis codes to identify cases of 12 infectious diseases in white patients: urinary tract infection, pneumonia, chronic sinus infections, otitis media, candidiasis, streptococcal pharyngitis, herpes zoster, herpes labialis, hepatitis B, infectious mononucleosis, tuberculosis (TB) or a positive TB test, and hepatitis C. We selected controls from patients with no diagnosis code for the candidate disease and matched by year of birth, sex, and calendar year at first and last EHR visits. We conducted GWAS using SAIGE and transcriptome-wide analysis (TWAS) using S-PrediXcan. We also conducted phenome-wide association study to understand associations between identified genetic variants and clinical phenotypes.

RESULTS

We replicated three 23andMe loci (p ≤ 0.05): herpes zoster and rs7047299-A (p = 2.6 × 10^-3) and rs2808290-C (p = 9.6 × 10^-3;); otitis media and rs114947103-C (p = 0.04). We also identified 2 novel regions (p ≤ 5 × 10^-8): rs113235453-G for otitis media (p = 3.04 × 10^-8), and rs10422015-T for candidiasis (p = 3.11 × 10^-8). In TWAS, four gene-disease associations were significant: SLC30A9 for otitis media (p = 8.06 × 10^-7); LRP3 and WDR88 for candidiasis (p = 3.91 × 10^-7 and p = 1.95 × 10^-6); and AAMDC for hepatitis B (p = 1.51 × 10^-6).

CONCLUSION

We conducted GWAS and TWAS for 12 infectious diseases and identified novel genetic contributors to the susceptibility of infectious diseases.

NFκB1 Polymorphisms Are Associated with Severe Influenza A (H1N1) Virus Infection in a Canadian Population

Background: We examined associations between NFκB1 polymorphisms and influenza A (H1N1) clinical outcomes in Canadian. Methods: A total of thirty-six Caucasian patients admitted to the intensive care unit (ICU) in hospitals in Canada were recruited during the 2009 H1N1 pandemic. Genomic DNA was extracted from the whole blood samples. The NFkB1 gene was targeted for genotyping using next-generation sequencing technology—Roche 454. Results: A total of 136 single nucleotide polymorphisms (SNPs) were discovered within the NFκB1 gene. Among them, 63 SNPs were significantly enriched in patients admitted in the ICU (p < 0.05) compared with the British Caucasian population in the 1000 Genomes study. These enriched SNPs are mainly intron variants, and only two are exon SNPs from the non-transcribing portion of the NFκB1 gene. Conclusions: Genetic variations in the NFκB1 gene could influence clinical outcomes of pandemic H1N1 infections. Our findings showed that sequence variations of the NFκB1 gene might influence patient response to influenza infection.

The human genetic epidemiology of COVID-19

Human genetics can inform the biology and epidemiology of coronavirus disease 2019 (COVID-19) by pinpointing causal mechanisms that explain why some individuals become more severely affected by the disease upon infection by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Large-scale genetic association studies, encompassing both rare and common genetic variants, have used different study designs and multiple disease phenotype definitions to identify several genomic regions associated with COVID-19. Along with a multitude of follow-up studies, these findings have increased our understanding of disease aetiology and provided routes for management of COVID-19. Important emergent opportunities include the clinical translatability of genetic risk prediction, the repurposing of existing drugs, exploration of variable host effects of different viral strains, study of inter-individual variability in vaccination response and understanding the long-term consequences of SARS-CoV-2 infection. Beyond the current pandemic, these transferrable opportunities are likely to affect the study of many infectious diseases.

Genome-wide association study of pain sensitivity assessed by questionnaire and the cold pressor test

ABSTRACT

We deployed an online pain sensitivity questionnaire (PSQ) and an at-home version of the cold pressor test (CPT) in a large genotyped cohort. We performed genome-wide association studies on the PSQ score (25,321 participants) and CPT duration (6853). We identified one new genome-wide significant locus associated with the PSQ score, which was located in the TSSC1 (also known as EIPR1 ) gene (rs58194899, OR = 0.950 [0.933-0.967], P -value = 1.9 × 10 -8 ). Although high pain sensitivity measured by both PSQ and CPT was associated with individual history of chronic and acute pains, genetic correlation analyses surprisingly suggested an opposite direction: PSQ score was inversely genetically correlated with neck and shoulder pain ( rg = -0.71), rheumatoid arthritis (-0.68), and osteoarthritis (-0.38), and with known risk factors, such as the length of working week (-0.65), smoking (-0.36), or extreme BMI (-0.23). Gene-based analysis followed by pathway analysis showed that genome-wide association studies results were enriched for genes expressed in the brain and involved in neuronal development and glutamatergic synapse signaling pathways. Finally, we confirmed that females with red hair were more sensitive to pain and found that genetic variation in the MC1R gene was associated with an increase in self-perceived pain sensitivity as assessed by the PSQ.

Infection Polygenic Factors Account for a Small Proportion of the Relationship Between Infections and Mental Disorders

BACKGROUND

Several recent studies have suggested a role for infections in the development of mental disorders; however, the genetic contribution to this association is understudied.

METHODS

We use the iPSYCH case-cohort genotyped sample (n = 65,534) and Danish health care registry data to study the genetic association between infections and mental disorders. To test the hypothesis that these associations are due to genetic pleiotropy, we estimated the genetic correlation between infection and mental disorders. Polygenic risk scores (PRSs) were used to assess whether genetic pleiotropy of infections and mental disorders was mediated by actual infection diagnoses.

RESULTS

We observed that schizophrenia, attention-deficit/hyperactivity disorder, major depressive disorder, bipolar disorder, and posttraumatic stress disorder (r_g ranging between 0.18 and 0.83), but not autism spectrum disorder and anorexia nervosa, were significantly genetically correlated with infection diagnoses. PRSs for infections were associated with modest increase in risk of attention-deficit/hyperactivity disorder, major depressive disorder, and schizophrenia in the iPSYCH case-cohort (hazard ratios = 1.04 to 1.10) but was not associated with risk of anorexia, autism, or bipolar disorder. Using mediation analysis, we show that infection diagnoses account for only a small proportion (6%-14%) of the risk for mental disorders conferred by infection PRSs.

CONCLUSIONS

Infections and mental disorders share a modest genetic architecture. Infection PRSs can predict risk of certain mental disorders; however, this effect is moderate. Finally, recorded infections partially explain the relationship between infection PRSs and mental disorders.

Susceptibility and Severity of COVID-19 Are Both Associated With Lower Overall Viral–Peptide Binding Repertoire of HLA Class I Molecules, Especially in Younger People

An important number of studies have been conducted on the potential association between human leukocyte antigen (HLA) genes and COVID-19 susceptibility and severity since the beginning of the pandemic. However, case–control and peptide-binding prediction methods tended to provide inconsistent conclusions on risk and protective HLA alleles, whereas some researchers suggested the importance of considering the overall capacity of an individual’s HLA Class I molecules to present SARS-CoV-2-derived peptides. To close the gap between these approaches, we explored the distributions of HLA-A, -B, -C, and -DRB1 1st-field alleles in 142 Iranian patients with COVID-19 and 143 ethnically matched healthy controls, and applied in silico predictions of bound viral peptides for each individual’s HLA molecules. Frequency comparison revealed the possible predisposing roles of HLA-A*03, B*35, and DRB1*16 alleles and the protective effect of HLA-A*32, B*58, B*55, and DRB1*14 alleles in the viral infection. None of these results remained significant after multiple testing corrections, except HLA-A*03, and no allele was associated with severity, either. Compared to peptide repertoires of individual HLA molecules that are more likely population-specific, the overall coverage of virus-derived peptides by one’s HLA Class I molecules seemed to be a more prominent factor associated with both COVID-19 susceptibility and severity, which was independent of affinity index and threshold chosen, especially for people under 60 years old. Our results highlight the effect of the binding capacity of different HLA Class I molecules as a whole, and the more essential role of HLA-A compared to HLA-B and -C genes in immune responses against SARS-CoV-2 infection.

The genetic architecture of pneumonia susceptibility implicates mucin biology and a relationship with psychiatric illness

Pneumonia remains one of the leading causes of death worldwide. In this study, we use genome-wide meta-analysis of lifetime pneumonia diagnosis (N = 391,044) to identify four association signals outside of the previously implicated major histocompatibility complex region. Integrative analyses and finemapping of these signals support clinically tractable targets, including the mucin MUC5AC and tumour necrosis factor receptor superfamily member TNFRSF1A. Moreover, we demonstrate widespread evidence of genetic overlap with pneumonia susceptibility across the human phenome, including particularly significant correlations with psychiatric phenotypes that remain significant after testing differing phenotype definitions for pneumonia or genetically conditioning on smoking behaviour. Finally, we show how polygenic risk could be utilised for precision treatment formulation or drug repurposing through pneumonia risk scores constructed using variants mapped to pathways with known drug targets. In summary, we provide insights into the genetic architecture of pneumonia susceptibility and genetics informed targets for drug development or repositioning.

The genetic and evolutionary determinants of COVID-19 susceptibility

Devastating pandemics, such as that due to COVID-19, can provide strong testimony to our knowledge of the genetic and evolutionary determinants of infectious disease susceptibility and severity. One of the most remarkable aspects of such outbreaks is the stunning interindividual variability observed in the course of infection. In recent decades, enormous progress has been made in the field of the human genetics of infectious diseases, and an increasing number of human genetic factors have been reported to explain, to a great extent, the observed variability for a large number of infectious agents. However, our understanding of the cellular, molecular, and immunological mechanisms underlying such disparities between individuals and ethnic groups, remains very limited. Here, we discuss recent findings relating to human genetic predisposition to infectious disease, from an immunological or population genetic perspective, and show how these and other innovative approaches have been applied to deciphering the genetic basis of human susceptibility to COVID-19 and the severity of this disease. From an evolutionary perspective, we show how past demographic and selection events characterizing the history of our species, including admixture with archaic humans, such as Neanderthals, facilitated modern human adaptation to the threats imposed by ancient pathogens. In the context of emerging infectious diseases, these past episodes of genetic adaptation may contribute to some of the observed population differences in the outcome of SARS-CoV-2 infection and the severity of COVID-19 illness.

Temporal changes in gastrointestinal fungi and the risk of autoimmunity during early childhood: the TEDDY study

Fungal infections are a major health problem that often begin in the gastrointestinal tract. Gut microbe interactions in early childhood are critical for proper immune responses, yet there is little known about the development of the fungal population from infancy into childhood. Here, as part of the TEDDY (The Environmental Determinants of Diabetes in the Young) study, we examine stool samples of 888 children from 3 to 48 months and find considerable differences between fungi and bacteria. The metagenomic relative abundance of fungi was extremely low but increased while weaning from milk and formula. Overall fungal diversity remained constant over time, in contrast with the increase in bacterial diversity. Fungal profiles had high temporal variation, but there was less variation from month-to-month in an individual than among different children of the same age. Fungal composition varied with geography, diet, and the use of probiotics. Multiple Candida spp. were at higher relative abundance in children than adults, while Malassezia and certain food-associated fungi were lower in children. There were only subtle fungal differences associated with the subset of children that developed islet autoimmunity or type 1 diabetes. Having proper fungal exposures may be crucial for children to establish appropriate responses to fungi and limit the risk of infection: the data here suggests those gastrointestinal exposures are limited and variable.

Fine-mapping studies distinguish genetic risks for childhood- and adult-onset asthma in the HLA region

BACKGROUND

Genome-wide association studies of asthma have revealed robust associations with variation across the human leukocyte antigen (HLA) complex with independent associations in the HLA class I and class II regions for both childhood-onset asthma (COA) and adult-onset asthma (AOA). However, the specific variants and genes contributing to risk are unknown.

METHODS

We used Bayesian approaches to perform genetic fine-mapping for COA and AOA (n=9432 and 21,556, respectively; n=318,167 shared controls) in White British individuals from the UK Biobank and to perform expression quantitative trait locus (eQTL) fine-mapping in immune (lymphoblastoid cell lines, n=398; peripheral blood mononuclear cells, n=132) and airway (nasal epithelial cells, n=188) cells from ethnically diverse individuals. We also examined putatively causal protein coding variation from protein crystal structures and conducted replication studies in independent multi-ethnic cohorts from the UK Biobank (COA n=1686; AOA n=3666; controls n=56,063).

RESULTS

Genetic fine-mapping revealed both shared and distinct causal variation between COA and AOA in the class I region but only distinct causal variation in the class II region. Both gene expression levels and amino acid variation contributed to risk. Our results from eQTL fine-mapping and amino acid visualization suggested that the HLA-DQA1*03:01 allele and variation associated with expression of the nonclassical HLA-DQA2 and HLA-DQB2 genes accounted entirely for the most significant association with AOA in GWAS. Our studies also suggested a potentially prominent role for HLA-C protein coding variation in the class I region in COA. We replicated putatively causal variant associations in a multi-ethnic cohort.

CONCLUSIONS

We highlight roles for both gene expression and protein coding variation in asthma risk and identified putatively causal variation and genes in the HLA region. A convergence of genomic, transcriptional, and protein coding evidence implicates the HLA-DQA2 and HLA-DQB2 genes and HLA-DQA1*03:01 allele in AOA.

Functional genomic analysis delineates regulatory mechanisms of GWAS-identified bipolar disorder risk variants

Background

Genome-wide association studies (GWASs) have identified multiple risk loci for bipolar disorder (BD). However, pinpointing functional (or causal) variants in the reported risk loci and elucidating their regulatory mechanisms remain challenging.

Methods

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) data from human brain tissues (or neuronal cell lines) and position weight matrix (PWM) data to identify functional single-nucleotide polymorphisms (SNPs). Then, we verified the regulatory effects of these transcription factor (TF) binding–disrupting SNPs (hereafter referred to as “functional SNPs”) through a series of experiments, including reporter gene assays, allele-specific expression (ASE) analysis, TF knockdown, CRISPR/Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis. Finally, we overexpressed PACS1 (whose expression was most significantly associated with the identified functional SNPs rs10896081 and rs3862386) in mouse primary cortical neurons to investigate if PACS1 affects dendritic spine density.

Results

We identified 16 functional SNPs (in 9 risk loci); these functional SNPs disrupted the binding of 7 TFs, for example, CTCF and REST binding was frequently disrupted. We then identified the potential target genes whose expression in the human brain was regulated by these functional SNPs through eQTL analysis. Of note, we showed dysregulation of some target genes of the identified TF binding–disrupting SNPs in BD patients compared with controls, and overexpression of PACS1 reduced the density of dendritic spines, revealing the possible biological mechanisms of these functional SNPs in BD.

Conclusions

Our study identifies functional SNPs in some reported risk loci and sheds light on the regulatory mechanisms of BD risk variants. Further functional characterization and mechanistic studies of these functional SNPs and candidate genes will help to elucidate BD pathogenesis and develop new therapeutic approaches and drugs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01057-3.

Human herpesvirus diversity is altered in HLA class I binding peptides

Herpesviruses are ubiquitous, genetically diverse DNA viruses, with long-term presence in humans associated with infrequent but significant pathology. Human leukocyte antigen (HLA) class I presents intracellularly derived peptide fragments from infected tissue cells to CD8+ T and natural killer cells, thereby directing antiviral immunity. Allotypes of highly polymorphic HLA class I are distinguished by their peptide binding repertoires. Because this HLA class I variation is a major determinant of herpesvirus disease, we examined if sequence diversity of virus proteins reflects evasion of HLA presentation. Using population genomic data from Epstein–Barr virus (EBV), human cytomegalovirus (HCMV), and Varicella–Zoster virus, we tested whether diversity differed between the regions of herpesvirus proteins that can be recognized, or not, by HLA class I. Herpesviruses exhibit lytic and latent infection stages, with the latter better enabling immune evasion. Whereas HLA binding peptides of lytic proteins are conserved, we found that EBV and HCMV proteins expressed during latency have increased peptide sequence diversity. Similarly, latent, but not lytic, herpesvirus proteins have greater population structure in HLA binding than nonbinding peptides. Finally, we found patterns consistent with EBV adaption to the local HLA environment, with less efficient recognition of EBV isolates by high-frequency HLA class I allotypes. Here, the frequency of CD8+ T cell epitopes inversely correlated with the frequency of HLA class I recognition. Previous analyses have shown that pathogen-mediated natural selection maintains exceptional polymorphism in HLA residues that determine peptide recognition. Here, we show that HLA class I peptide recognition impacts diversity of globally widespread pathogens.

Association of Prenatal Antibiotics and Mode of Birth With Otolaryngology Surgery in Offspring: A National Data Linkage Study

BACKGROUND

Pediatric otolaryngology surgery is commonly performed after recurrent infections and allergy/atopy. Prenatal antibiotic exposure and cesarean section deliveries increase the risk of severe infection and allergy/atopy in the offspring, but the relationship with common, related surgical outcomes is unknown. This study measures the associations between prenatal antibiotic use and mode of birth with common pediatric otolaryngology surgery.

METHODS

Data linkage analysis of all live-born, singleton children, born between 2008 and 2018 was done using Norwegian national health registry data. Exposures of interest were prenatal antibiotics and mode of birth. The primary outcome was common otolaryngology surgery before 10 years of age. Exposure-outcome associations were estimated through multivariable Cox proportional hazards models adjusting for predefined covariates. Interaction between exposures was explored.

RESULTS

Of 539,390 children, 146,832 (27.2%) had mothers who were prescribed antibiotics during pregnancy, 83,473 (15.5%) were delivered via cesarean section, and 48,565 (9.0%) underwent an otolaryngology surgery during the study period. Prenatal antibiotic exposure [adjusted hazard ratio (aHR), 1.22; 95% CI: 1.20-1.24] and cesarean section (aHR, 1.14; 95% CI: 1.11-1.16) were each associated with otolaryngology surgery after mutual adjustment. There was some evidence of an interaction between the 2 exposures (P = 0.03).

CONCLUSIONS

Antibiotic exposure in pregnancy and cesarean section may adversely affect early immune development and increase the risk of recurrent upper airway infections and allergy/atopy that may require otolaryngology surgery. Mechanistic studies are warranted to explore genetic and/or molecular pathways that explain these findings. This may identify potential therapeutic targets to reduce the burden of otolaryngology surgery.

Risk of lower respiratory tract infections: a genome-wide association study with Mendelian randomization analysis in three independent European populations

OBJECTIVE

Lower respiratory tract infections (LRTIs) are a leading cause of morbidity and mortality worldwide. Few studies have previously investigated genetic susceptibility and potential risk factors for LRTI.

METHODS

We used data from the UK Biobank, Trøndelag Health Study (HUNT), and FinnGen to conduct a genome-wide association study (GWAS). Cases were subjects hospitalized with LRTI, and controls were subjects with no such hospitalization. We conducted stratification and interaction analyses to evaluate whether the genetic effect of LRTI differed by sex or smoking. Mendelian randomization (MR) analyses were conducted to identify the unconfounded relationship between cardiometabolic risk factors and LRTI.

RESULTS

A total of 25 320 cases and 575 294 controls were included. The 15q25.1 locus reached genome-wide significance in the meta-analysis (rs10519203: OR 0.94, p 3.87e-11). The protective effect of effect allele of rs10519203 was present among smokers (OR 0.90, 95%CI 0.87-0.92, p 1.38e-15) but not among never-smokers (OR 1.01, 95%CI 0.97-1.06, p 5.20e-01). In MR analyses, we found that increasing body mass index (OR 1.31, 95%CI 1.24-1.40, p 3.78e-18), lifetime smoking (OR 2.83, 95%CI 2.34-3.42, p 6.56e-27), and systolic blood pressure robustly increased the risk of LRTIs (OR 1.11, 95%CI 1.02-1.22, p 1.48e-02).

CONCLUSION

A region in 15q25.1 was strongly associated with LRTI susceptibility. Reduction in the prevalence of smoking, overweight, obesity, and hypertension may reduce the disease burden of LRTIs.

HLA-DR polymorphism in SARS-CoV-2 infection and susceptibility to symptomatic COVID-19

SARS-CoV-2 infection results in different outcomes ranging from asymptomatic infection to mild or severe disease and death. Reasons for this diversity of outcome include differences in challenge dose, age, gender, comorbidity and host genomic variation. Human leukocyte antigen (HLA) polymorphisms may influence immune response and disease outcome. We investigated the association of HLAII alleles with case definition symptomatic COVID-19, virus-specific antibody and T-cell immunity. A total of 1364 UK healthcare workers (HCWs) were recruited during the first UK SARS-CoV-2 wave and analysed longitudinally, encompassing regular PCR screening for infection, symptom reporting, imputation of HLAII genotype and analysis for antibody and T-cell responses to nucleoprotein (N) and spike (S). Of 272 (20%) HCW who seroconverted, the presence of HLA-DRB1*13:02 was associated with a 6·7-fold increased risk of case definition symptomatic COVID-19. In terms of immune responsiveness, HLA-DRB1*15:02 was associated with lower nucleocapsid T-cell responses. There was no association between DRB1 alleles and anti-spike antibody titres after two COVID vaccine doses. However, HLA DRB1*15:01 was associated with increased spike T-cell responses following both first and second dose vaccination. Trial registration: NCT04318314 and ISRCTN15677965.

Understanding the genetic basis of immune responses to fungal infection

INTRODUCTION

Fungal infections represent a global public health problem that affect millions of people. Despite remarkable advances achieved over the last decades, available diagnostic and therapeutic tools remain insufficient for the optimal management of these diseases. The clinical course of fungal infection is highly variable, and evidence accumulated from patients with rare mutations and cohort-based studies suggests that the trajectory of disease is largely defined by patient genetics and its impact on immune responses. Therefore, there is an urgent need to elucidate the precise mechanisms by which which genetic variants influence the risk, progression, and outcome of fungal infection.

AREAS COVERED

In this review, we highlight recent advances in our understanding of the genetic factors that influence antifungal immune responses based on candidate gene studies and genome-wide approaches performed in different experimental and clinical models.

EXPERT OPINION

Research on genetics of susceptibility to infection is expected to lead to a detailed knowledge framework for the pathogenesis of human fungal infections and unveil novel targets and pathways amenable to clinical intervention.

One genome, many cell states: epigenetic control of innate immunity

A hallmark of the innate immune system is its ability to rapidly initiate short-lived or sustained transcriptional programs in a cell-specific and pathogen-specific manner that is dependent on dynamic chromatin states. Much of the epigenetic landscape is set during cellular differentiation; however, pathogens and other environmental cues also induce changes in chromatin that can either promote tolerance or 'train' innate immune cells for amplified secondary responses. We review chromatin processes that enable innate immune cell differentiation and functional transcriptional responses in naive or experienced cells, in concert with signal transduction and cellular metabolic shifts. We discuss how immune chromatin mechanisms are maladapted in disease and novel therapeutic approaches for cellular reprogramming.

Diet and Respiratory Infections: Specific or Generalized Associations?

Background: Based on our recently reported associations between specific dietary behaviors and the risk of COVID-19 infection in the UK Biobank (UKB) cohort, we further investigate whether these associations are specific to COVID-19 or extend to other respiratory infections. Methods: Pneumonia and influenza diagnoses were retrieved from hospital and death record data linked to the UKB. Baseline, self-reported (2006−2010) dietary behaviors included being breastfed as a baby and intakes of coffee, tea, oily fish, processed meat, red meat (unprocessed), fruit, and vegetables. Logistic regression estimated the odds of pneumonia/influenza from baseline to 31 December 2019 with each dietary component, adjusting for baseline socio-demographic factors, medical history, and other lifestyle behaviors. We considered effect modification by sex and genetic factors related to pneumonia, COVID-19, and caffeine metabolism. Results: Of 470,853 UKB participants, 4.0% had pneumonia and 0.2% had influenza during follow up. Increased consumption of coffee, tea, oily fish, and fruit at baseline were significantly and independently associated with a lower risk of future pneumonia events. Increased consumption of red meat was associated with a significantly higher risk. After multivariable adjustment, the odds of pneumonia (p ≤ 0.001 for all) were lower by 6−9% when consuming 1−3 cups of coffee/day (vs. <1 cup/day), 8−11% when consuming 1+ cups of tea/day (vs. <1 cup/day), 10−12% when consuming oily fish in higher quartiles (vs. the lowest quartile—Q1), and 9−14% when consuming fruit in higher quartiles (vs. Q1); it was 9% higher when consuming red meat in the fourth quartile (vs. Q1). Similar patterns of associations were observed for influenza but only associations with tea and oily fish met statistical significance. The association between fruit and pneumonia risk was stronger in women than in men (p = 0.001 for interaction). Conclusions: In the UKB, consumption of coffee, tea, oily fish, and fruit were favorably associated with incident pneumonia/influenza and red meat was adversely associated. Findings for coffee parallel those we reported previously for COVID-19 infection, while other findings are specific to these more common respiratory infections.

Groh A, S. Frick J, Vehreschild MJGT, U. Ludwig K. Genetic Predisposition and the Variable Course of Infectious Diseases

BACKGROUND

Contact with a pathogen is followed by variable courses of infectious disease, which are only partly explicable by classical risk factors. The susceptibility to infection is variable, as is the course of disease after infection. In this review, we discuss the extent to which this variation is due to genetic factors of the affected individual (the host).

METHODS

Selective review of the literature on host genetics in infectious disease, with special attention to the pathogens SARSCoV- 2, influenza viruses, Mycobacterium tuberculosis, and human immunodeficiency virus (HIV).

RESULTS

Genetic variants of the host contribute to the pathogenesis of infectious diseases. For example, in HIV infection, a relatively common variant leading to a loss of function of the HIV co-receptor CCR5 affects the course of the disease, as do variants in genes of the major histocompatibility complex (MHC) region. Rare monogenic variants of the interferon immune response system contribute to severe disease courses in COVID-19 and influenza (type I interferon in these two cases) and in tuberculosis (type II interferon). An estimated 1.8% of life-threatening courses of COVID-19 in men under age 60 are caused by a deficiency of toll-like receptor 7. The scientific understanding of host genetic factors has already been beneficial to the development of effective drugs. In a small number of cases, genetic information has also been used for individual therapeutic decision-making and for the identification of persons at elevated risk.

CONCLUSION

A comprehensive understanding of host genetics can improve the care of patients with infectious diseases. Until the present, the clinical utility of host genetics has been limited to rare cases; in the future, polygenic risk scores summarizing the relevant genetic variants in each patient will enable a wider benefit. To make this possible, multicenter studies are needed that will systematically integrate clinical and genetic data.

Analysis of Gene Expression Microarray Data Reveals Androgen-Responsive Genes of Muscles in Polycystic Ovarian Syndrome Patients

Polycystic ovarian syndrome (PCOS) is an endocrine disorder that is characterized by hyperandrogenism. Therefore, information about androgen-induced molecular changes can be obtained using the tissues of patients with PCOS. We analyzed two microarray datasets of normal and PCOS muscle samples (GSE8157 and GSE6798) to identify androgen-responsive genes (ARGs). Differentially expressed genes were determined using the t-test and a meta-analysis of the datasets. The overlap between significant results of the meta-analysis and ARGs predicted from an external database was determined, and differential coexpression analysis was then applied between these genes and the other genes. We found 313 significant genes in the meta-analysis using the Benjamini–Hochberg multiple testing correction. Of these genes, 61 were in the list of predicted ARGs. When the differential coexpression between these 61 genes and 13,545 genes filtered by variance was analyzed, 540 significant gene pairs were obtained using the Benjamini–Hochberg correction. While no significant results were obtained regarding the functional enrichment of the differentially expressed genes, top-level gene ontology terms were significantly enriched in the list of differentially coexpressed genes, which indicates that a broad range of cellular processes is affected by androgen administration. Our findings provide valuable information for the identification of ARGs.

Functional genomics elucidates regulatory mechanisms of Parkinson's disease-associated variants

BACKGROUND

Genome-wide association studies (GWASs) have identified multiple risk loci for Parkinson's disease (PD). However, identifying the functional (or potential causal) variants in the reported risk loci and elucidating their roles in PD pathogenesis remain major challenges. To identify the potential causal (or functional) variants in the reported PD risk loci and to elucidate their regulatory mechanisms, we report a functional genomics study of PD.

METHODS

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) (from neuronal cells and human brain tissues) data and GWAS-identified single-nucleotide polymorphisms (SNPs) in PD risk loci. We then conducted a series of experiments and analyses to validate the regulatory effects of these (i.e., functional) SNPs, including reporter gene assays, allele-specific expression (ASE), transcription factor (TF) knockdown, CRISPR-Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis.

RESULTS

We identified 44 SNPs (from 11 risk loci) affecting the binding of 12 TFs and we validated the regulatory effects of 15 TF binding-disrupting SNPs. In addition, we also identified the potential target genes regulated by these TF binding-disrupting SNPs through eQTL analysis. Finally, we showed that 4 eQTL genes of these TF binding-disrupting SNPs were dysregulated in PD cases compared with controls.

CONCLUSION

Our study systematically reveals the gene regulatory mechanisms of PD risk variants (including widespread disruption of CTCF binding), generates the landscape of potential PD causal variants, and pinpoints promising candidate genes for further functional characterization and drug development.

A genome-wide association study in a large community-based cohort identifies multiple loci associated with susceptibility to bacterial and viral infections

There is limited data on host-specific genetic determinants of susceptibility to bacterial and viral infections. Genome-wide association studies using large population cohorts can be a first step towards identifying patients prone to infectious diseases and targets for new therapies. Genetic variants associated with clinically relevant entities of bacterial and viral infections (e.g., abdominal infections, respiratory infections, and sepsis) in 337,484 participants of the UK Biobank cohort were explored by genome-wide association analyses. Cases (n = 81,179) were identified based on ICD-10 diagnosis codes of hospital inpatient and death registries. Functional annotation was performed using gene expression (eQTL) data. Fifty-seven unique genome-wide significant loci were found, many of which are novel in the context of infectious diseases. Some of the detected genetic variants were previously reported associated with infectious, inflammatory, autoimmune, and malignant diseases or key components of the immune system (e.g., white blood cells, cytokines). Fine mapping of the HLA region revealed significant associations with HLA-DQA1, HLA-DRB1, and HLA-DRB4 locus alleles. PPP1R14A showed strong colocalization with abdominal infections and gene expression in sigmoid and transverse colon, suggesting causality. Shared significant loci across infections and non-infectious phenotypes in the UK Biobank cohort were found, suggesting associations for example between SNPs identified for abdominal infections and CRP, rheumatoid arthritis, and diabetes mellitus. We report multiple loci associated with bacterial and viral infections. A better understanding of the genetic determinants of bacterial and viral infections can be useful to identify patients at risk and in the development of new drugs.

No evidence that herpes zoster is associated with increased risk of dementia diagnosis

OBJECTIVE

To investigate whether herpes zoster (HZ) was associated with subsequent increased risk of dementia diagnosis.

METHODS

We conducted a historical cohort study using primary care electronic health records from the Clinical Practice Research Datalink in the United Kingdom. Individuals with incident HZ aged ≥40 years from 2000 to 2017 were matched with up to four individuals without HZ by age, sex, primary care practise and calendar time. The primary outcome was a new diagnosis of all-cause dementia. We used the Cox proportional hazards regression adjusting for demographic, lifestyle and clinical confounders to assess any association between HZ and dementia. We investigated interactions with sex, frailty index and antiviral treatment and conducted various sensitivity analyses.

RESULTS

The cohort comprised 177,144 individuals with HZ and 706,901 matched unexposed individuals (median age 65 years (IQR 55.1-75.0), 40% male) followed for a median duration of 4.6 years (IQR 2.0-8.1). In total, 26,585 (3%) patients had an incident dementia diagnosis recorded and 113,056 patients died (12.8%). HZ was associated with a small reduction in dementia diagnosis (adjusted HR 0.92 (95% CI 0.89-0.95)), occurring predominantly in frail individuals and females. For patients who were fit (578,115, 65%), no association was seen (adjusted HR 0.97, 95% CI 0.92-1.02). There was no association between HZ and a composite outcome of dementia or death (adjusted HR 1.00, 95% CI 0.99-1.02). Dementia risk did not vary by prescription of antiviral agents. Sensitivity analyses showed consistent results.

INTERPRETATION

HZ was not associated with increased dementia diagnosis in a UK primary care-based cohort.

Trans-ancestral fine-mapping of MHC reveals key amino acids associated with spontaneous clearance of hepatitis C in HLA-DQβ1

Spontaneous clearance of acute hepatitis C virus (HCV) infection is associated with single nucleotide polymorphisms (SNPs) on the MHC class II. We fine-mapped the MHC region in European (n = 1,600; 594 HCV clearance/1,006 HCV persistence) and African (n = 1,869; 340 HCV clearance/1,529 HCV persistence) ancestry individuals and evaluated HCV peptide binding affinity of classical alleles. In both populations, HLA-DQβ1Leu26 (p valueMeta = 1.24 × 10-14) located in pocket 4 was negatively associated with HCV spontaneous clearance and HLA-DQβ1Pro55 (p valueMeta = 8.23 × 10-11) located in the peptide binding region was positively associated, independently of HLA-DQβ1Leu26. These two amino acids are not in linkage disequilibrium (r2 < 0.1) and explain the SNPs and classical allele associations represented by rs2647011, rs9274711, HLA-DQB1∗03:01, and HLA-DRB1∗01:01. Additionally, HCV persistence classical alleles tagged by HLA-DQβ1Leu26 had fewer HCV binding epitopes and lower predicted binding affinities compared to clearance alleles (geometric mean of combined IC50 nM of persistence versus clearance; 2,321 nM versus 761.7 nM, p value = 1.35 × 10-38). In summary, MHC class II fine-mapping revealed key amino acids in HLA-DQβ1 explaining allelic and SNP associations with HCV outcomes. This mechanistic advance in understanding of natural recovery and immunogenetics of HCV might set the stage for much needed enhancement and design of vaccine to promote spontaneous clearance of HCV infection.

The FUT2 Variant c.461G>A (p.Trp154*) Is Associated With Differentially Expressed Genes and Nasopharyngeal Microbiota Shifts in Patients With Otitis Media

Otitis media (OM) is a leading cause of childhood hearing loss. Variants in FUT2, which encodes alpha-(1,2)-fucosyltransferase, were identified to increase susceptibility to OM, potentially through shifts in the middle ear (ME) or nasopharyngeal (NP) microbiotas as mediated by transcriptional changes. Greater knowledge of differences in relative abundance of otopathogens in carriers of pathogenic variants can help determine risk for OM in patients. In order to determine the downstream effects of FUT2 variation, we examined gene expression in relation to carriage of a common pathogenic FUT2 c.461G>A (p.Trp154*) variant using RNA-sequence data from saliva samples from 28 patients with OM. Differential gene expression was also examined in bulk mRNA and single-cell RNA-sequence data from wildtype mouse ME mucosa after inoculation with non-typeable Haemophilus influenzae (NTHi). In addition, microbiotas were profiled from ME and NP samples of 65 OM patients using 16S rRNA gene sequencing. In human carriers of the FUT2 variant, FN1, KMT2D, MUC16 and NBPF20 were downregulated while MTAP was upregulated. Post-infectious expression in the mouse ME recapitulated these transcriptional differences, with the exception of Fn1 upregulation after NTHi-inoculation. In the NP, Candidate Division TM7 was associated with wildtype genotype (FDR-adj-p=0.009). Overall, the FUT2 c.461G>A variant was associated with transcriptional changes in processes related to response to infection and with increased load of potential otopathogens in the ME and decreased commensals in the NP. These findings provide increased understanding of how FUT2 variants influence gene transcription and the mucosal microbiota, and thus contribute to the pathology of OM.

Tonsillectomy in Adults over 40 Years of Age Does Not Increase the Risk of Pneumonia: A Three-Year Longitudinal Follow-Up Study

To evaluate the effect of tonsillectomy on the subsequent risk of pneumonia in an adult population, a longitudinal follow-up case control study was conducted using a national health screening cohort dataset between 2003 and 2012. A total of 1005 tonsillectomy participants were 1:4 matched with 4020 control participants for age, sex, income, and region of residence. The number of pneumonia diagnoses were counted from the index date (ID) to the date after the first-year (post-ID 1y), second-year (post-ID 2y), and third-year (post-ID 3y) periods. Simple linear regression and multiple linear regression were conducted to calculate estimated values (EVs) and 95% confidence intervals for each post-ID pneumonia and compared between the two groups. Subgroup analyses were performed according to age, sex, and the number of pneumonia cases during the year prior to the ID (pre-ID 1y). In the simple linear regression model, post-ID pneumonia did not show a significant correlation with tonsillectomy (post-ID 1y: EV = 0.003; post-ID 2y: EV = 0.007; post-ID 3y: EV = 0.013; all p > 0.05). In the multiple regression model, post-ID pneumonia also did not show a significant correlation with tonsillectomy (post-ID 1y: EV = 0.001; post-ID 2y: EV = 0.006; post-ID 3y: EV = 0.011; all p > 0.05). In the subgroup analyses, tonsillectomy did not show a significant correlation with post-ID pneumonia in either the simple linear regression or multiple linear regression models (all p > 0.05). Tonsillectomy performed in the adult population did not show any effect in increasing the incidence of pneumonia during the first three postoperative years.