Host Genetic Variation Research Articles

Effects of Coffea canephora genotypes on the microbial community of soil and fruit

In recent years, the role of microbial communities in agricultural systems has received increasing attention, particularly concerning their impact on plant health and productivity. However, the influence of host plant genetic factors on the microbial composition of coffee plants remains largely unexplored. This study provides the first comprehensive investigation into how genotype affects the microbial communities present in the rhizosphere and fruits of Coffea canephora. Conducted on a commercial coffee farm in Brazil, we analyzed six genotypes of C. canephora var. Conilon. Soil and fruit samples were collected from which microbial DNA was extracted and sequenced, targeting the V3-V4 region of the 16 S rDNA and the ITS1 region for fungi. A total of 12,239,769 reads were generated from the 16 S rDNA and ITS1 regions. The PCoA revealed distinct patterns of beta diversity, with genotype 153 exhibiting significant isolation in soil bacterial communities. The dominant bacterial orders included Rhizobiales and Rhodobacterales, while the fungal community comprised diverse taxa from Saccharomycetales and Hypocreales. LEfSe analysis identified key metagenomic biomarkers, highlighting genotype Baiano 4 for its richness in fruit-associated taxa, whereas genotype 153 exhibited lower diversity in both soil and fruit samples. This work enhances our understanding of the microbiomes associated with different coffee genotypes, providing evidence of how host genetic variation influences microbial community composition. Our findings indicate that specific microbial taxa are enriched in the fruits and soil of various genotypes. Future research should focus on identifying these microorganisms and elucidating their specific functions within the rhizosphere and coffee fruits.

Interleukin-1 Beta rs16944 and rs1143634 and Interleukin-6 Receptor rs12083537 Single Nucleotide Polymorphisms as Potential Predictors of COVID-19 Severity.

Host genetic variation has been recognized as a key predictor of diverse clinical sequelae among severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected patients. Insights into the link between the Interleukin-6 receptor (IL-6R) and Interleukin-1 beta (IL-1β) genetic variation and severe coronavirus disease 2019 (COVID-19) are crucial for developing new predictors and therapeutic targets. We aimed to investigate the association of IL-6R rs12083537, IL-1β rs16944, and IL-1β rs1143634 SNPs with the severity of COVID-19. Our study was conducted on 300 COVID-19-negative individuals (control group) and 299 COVID-19-positive cases, classified into mild, moderate, and severe subgroups. Analyses of IL-1β (rs16944, rs1143634) and IL-6R (rs12083537) SNPs' genotypes were performed using qPCR genotyping assays. The IL-1β (rs16944) CC genotype and IL-6R (rs12083537) GG genotype were substantially related to COVID-19 severity, which was also associated with comorbidities and some laboratory parameters (p < 0.001). The IL-1β (rs1143634) TT genotype was found to be protective. Likewise, the IL-1β (rs16944) CC genotype was associated with increased mortality. IL-1β rs16944 and IL-6R rs12083537 SNPs are promising potential predictors of SARS-CoV-2 disease severity. Meanwhile, the rs1143634 SNP T allele was protective against severity and mortality risk.

Host genetic diversity and body condition influence parasite resistance and clearance in a wild marine mammal population.

Host genetic variability can modulate infection resistance, although its role in infection clearance remains unclear. Hookworm disease (Uncinaria sp.) is the leading cause of pup mortality in several otariid species, although the parasite can be cleared through immune-mediated processes. We evaluated the association of host genetic diversity, body condition and immune response with hookworm resistance and/or clearance in the South American fur seal (Arctocephalus australis). Uninfected pups had higher heterozygosity than parasitized individuals, indicating a negative relationship between heterozygosity and the chances of infection. Likewise, pups that died of hookworm infection had lower heterozygosity than those that died of non-infectious causes. Interestingly, once infected, pups that survived hookworm infection had heterozygosities similar to pups that died of hookworm disease. However, pups that cleared the infection had a higher body mass and parasite-specific immunoglobulin G levels than those that did not recover or died of hookworm disease. Thus, although heterozygosity predicted resistance to and mortality from hookworm infections, it did not affect parasite clearance, which was facilitated by better body condition and adaptive immune responses. This demonstrates that host genetic variability and host-environment interactions influence disease dynamics, acting at different, well-defined stages of infection.

Impact of Human Immunological Responses and Viral Genetic Diversity on Outbreak of Human Monkeypox Virus. A Comprehensive Literature Review Study

Monkeypox was caused by Monkeypox Virus (MPXV) and can infect both humans and animals. An understanding of the interplay between host immunity and genetic diversity was necessary to understand the etiology and epidemiology of monkeypox disease. Objective: To clarify how genetic differences and host immune responses interact when a moneypox infection occurs. Furthermore, we also aim to provide insights into individual variability in illness outcomes and possible treatment targets by investigating how distinct genetic profiles affect immune system activation and efficacy. Methods: Recent research on monkeypox, concentrating on the immune response mechanisms of the host and genetic variables linked to virus vulnerability have thoroughly analyzed. For this purpose, the data were searched from various research engines such as google scholar, pubmed, medline etc., by using different key words i.e., monkeypox and host immunity, monkeypox and antibodies interactions, monkeypox outbreak, mnkeypox strains. Conclusions: The way a monkeypox infection progresses and turns out was greatly influenced by the interplay between host genetic differences and immunological responses. Public health initiatives and the creation of tailored treatment plans can both benefit from the identification of genetic markers linked to immunological response profiles and vulnerability.

Infection of In Vivo and In Vitro Pines with the Pinewood Nematode Bursaphelenchus xylophilus and Isolation of Induced Volatiles.

The pinewood nematode (PWN) is a phytoparasite that causes pine wilt disease (PWD) in conifer species. This plant parasitic nematode has heavily contributed to pine deforestation in Asian countries, e.g., Japan, China, and Korea. Over the last two decades, in Europe, Portugal and Spain have been greatly affected. Research on the mechanisms of PWN infection and/or PWD progression in susceptible host species relies on the controlled infection of pine seedlings under greenhouse conditions. This technique is laborious and mobilizes substantial economic and human resources. Additionally, it can be prone to variability that results from the genetic diversity associated with some pine species but also from the interference of external factors. As an alternative, in vitro co-cultures of pine with PWNs offer a more advantageous system for studying biochemical changes since they a) allow controlling single environmental or nutritional variables, b) occupy less space, c) require less time to obtain, and d) arefree from contamination or from host genetic variation. The following protocol details the standard in vivo PWN infection of Pinus pinaster, the maritime pine, and the establishment of the novel in vitro co-cultures of pine shoots with the PWN as an improved methodology to study this phytoparasite influence on pine volatiles. PWN-induced volatiles are extracted from in vivo and in vitro infected pines by hydrodistillation and distillation-extraction, and the emitted volatiles are captured by solid phase microextraction (SPME), using fiber or packed column techniques.

PSVI-6 Genomic determinants of the rumen microbiome and metabolome in Angus steers

Abstract The rumen microbiome provides approximately 70% of energy required for the ruminant host through fermentation byproducts. Due to microbial impacts on host health and productivity, several studies have associated microbiome and metabolome variation with important sustainability phenotypes, such as feed efficiency and methane emissions. Recent studies have indicated the microbial communities present in rumen fluid are under low-to-moderate host genetic control. However, microorganisms present on the rumen wall have been underrepresented in research, negating their unique role in host-microbiome interactions. Current studies investigating heritability of the rumen planktonic microorganisms are plagued by small sample sizes, high management and breed heterogeneity, and low resolution for classifying microbial communities. To determine the heritability of the rumen fluid and epithelium-associated microbiome and metabolome, Angus steers (n = 350) were enrolled in a 70-d feed efficiency trial, where real-time dry matter intake was measured. Body weights were collected throughout the trial to calculate average daily gain and residual feed intake. Mid-trial, a novel surgery method was used to swab the rumen wall to capture epithelium-associated microorganisms. On d 70, 50 mL of rumen fluid was obtained from steers using orogastric tubing for collection of planktonic microbial communities. Metagenomic DNA extracted from both rumen fluid and epimural samples was used for reduced representation sequencing to predict microbial features such as diversity, composition, and putative function. Genomic DNA isolated from whole blood was used to perform low-pass whole genome sequencing for use in downstream heritability estimates and variant association analyses. Many core microbial abundances are under low-to-moderate host genetic control (h2 &amp;gt; 0.15; P &amp;lt; 0.05); however, many microorganisms had low heritability estimates (h2 &amp;lt; 0.10). Preliminary genome-wide associations studies are vastly underpowered to detect host genetic variants associated with microbial features, but data generated using the total number of animals will better refine these polymorphisms. Several key heritable metabolites were identified in amino acid metabolism pathways in rumen fluid (P &amp;lt; 0.05). The rumen microbiome appears to act as a complex polygenic trait and incorporating microbial features into prediction models may improve feed efficiency estimations. Understanding host genetic control of the rumen microbiome and metabolome may empower downstream applications to enhance genomic and phenotypic predictions for feed efficiency and other sustainability traits in beef cattle.

Host genetics and gut microbiota synergistically regulate feed utilization in egg-type chickens

BackgroundFeed efficiency is a crucial economic trait in poultry industry. Both host genetics and gut microbiota influence feed efficiency. However, the associations between gut microbiota and host genetics, as well as their combined contributions to feed efficiency in laying hens during the late laying period, remain largely unclear.MethodsIn total, 686 laying hens were used for whole-genome resequencing and liver transcriptome sequencing. 16S rRNA gene sequencing was conducted on gut chyme (duodenum, jejunum, ileum, and cecum) and fecal samples from 705 individuals. Bioinformatic analysis was performed by integrating the genome, transcriptome, and microbiome to screen for key genetic variations, genes, and gut microbiota associated with feed efficiency.ResultsThe heritability of feed conversion ratio (FCR) and residual feed intake (RFI) was determined to be 0.28 and 0.48, respectively. The ileal and fecal microbiota accounted for 15% and 10% of the FCR variance, while the jejunal, cecal, and fecal microbiota accounted for 20%, 11%, and 10% of the RFI variance. Through SMR analysis based on summary data from liver eQTL mapping and GWAS, we further identified four protein-coding genes, SUCLA2, TNFSF13B, SERTM1, and MARVELD3, that influence feed efficiency in laying hens. The SUCLA2 and TNFSF13B genes were significantly associated with SNP 1:25664581 and SNP rs312433097, respectively. SERTM1 showed significant associations with rs730958360 and 1:33542680 and is a potential causal gene associated with the abundance of Corynebacteriaceae in feces. MARVELD3 was significantly associated with the 1:135348198 and was significantly correlated with the abundance of Enterococcus in ileum. Specifically, a lower abundance of Enterococcus in ileum and a higher abundance of Corynebacteriaceae in feces were associated with better feed efficiency.ConclusionsThis study confirms that both host genetics and gut microbiota can drive variations in feed efficiency. A small portion of the gut microbiota often interacts with host genes, collectively enhancing feed efficiency. Therefore, targeting both the gut microbiota and host genetic variation by supporting more efficient taxa and selective breeding could improve feed efficiency in laying hens during the late laying period.

Torque Teno Sus Virus 1: A Potential Surrogate Pathogen to Study Pig-Transmitted Transboundary Animal Diseases.

Understanding the epidemiology and transmission dynamics of transboundary animal diseases (TADs) among wild pigs (Sus scrofa) will aid in preventing the introduction or containment of TADs among wild populations. Given the challenges associated with studying TADs in free-ranging populations, a surrogate pathogen system may predict how pathogens may circulate and be maintained within wild free-ranging swine populations, how they may spill over into domestic populations, and how management actions may impact transmission. We assessed the suitability of Torque teno sus virus 1 (TTSuV1) to serve as a surrogate pathogen for molecular epidemiological studies in wild pigs by investigating the prevalence, persistence, correlation with host health status and genetic variability at two study areas: Archbold's Buck Island Ranch in Florida and Savannah River Site in South Carolina. We then conducted a molecular epidemiological case study within Archbold's Buck Island Ranch site to determine how analysis of this pathogen could inform transmission dynamics of a directly transmitted virus. Prevalence was high in both study areas (40%, n = 190), and phylogenetic analyses revealed high levels of genetic variability within and between study areas. Our case study showed that pairwise host relatedness and geographic distance were highly correlated to pairwise viral genetic similarity. Molecular epidemiological analyses revealed a distinct pattern of direct transmission from pig to pig occurring within and between family groups. Our results suggest that TTSuV1 is highly suitable for molecular epidemiological analyses and will be useful for future studies of transmission dynamics in wild free-ranging pigs.

Analysis of hepatic lentiviral vector transduction; implications for preclinical studies and clinical gene therapy protocols.

Lentiviral vector-transduced T-cells were approved by the FDA as gene therapy anti-cancer medications. Little is known about the host genetic variation effects on the safety and efficacy of the lentiviral vector gene delivery system. To narrow this knowledge-gap, we characterized hepatic gene delivery by lentiviral vectors across the Collaborative Cross (CC) mouse genetic reference population. For 24 weeks, we periodically measured hepatic luciferase expression from lentiviral vectors in 41 CC mouse strains. Hepatic and splenic vector copy numbers were determined. We report that CC mouse strains showed highly diverse outcomes following lentiviral gene delivery. For the first time, moderate correlation between mouse strain-specific sleeping patterns and transduction efficiency was observed. We associated two quantitative trait loci (QTLs) with intra-strain variations in transduction phenotypes, which mechanistically relates to the phenomenon of metastable epialleles. An additional QTL was associated with the kinetics of hepatic transgene expression. Genes comprised in the above QTLs are potential targets to personalize gene therapy protocols. Importantly, we identified two mouse strains that open new directions in characterizing continuous viral vector silencing and HIV latency. Our findings suggest that wide-range patient-specific outcomes of viral vector-based gene therapy should be expected. Thus, novel escalating dose-based clinical protocols should be considered.

Rare host variants in ciliary expressed genes contribute to COVID-19 severity in Bulgarian patients

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a pneumonia with extremely heterogeneous clinical presentation, ranging from asymptomatic to severely ill patients. Previous studies have reported links between the presence of host genetic variants and the outcome of the COVID-19 infection. In our study, we used whole exome sequencing in a cohort of 444 SARS-CoV-2 patients, admitted to hospital in the period October-2020–April-2022, to search for associations between rare pathogenic/potentially pathogenic variants and COVID-19 progression. We used gene prioritization-based analysis in genes that have been reported by host genetic studies. Although we did not identify correlation between the presence of rare pathogenic variants and COVID-19 outcome, in critically ill patients we detected known mutations in a number of genes associated with severe disease related to cardiovascular disease, primary ciliary dyskinesia, cystic fibrosis, DNA damage repair response, coagulation, primary immune disorder, hemoglobin subunit β, and others. Additionally, we report 93 novel pathogenic variants found in severely infected patients who required intubation or died. A network analysis showed main component, consisting of 13 highly interconnected genes related to epithelial cilium. In conclusion, we have detected rare pathogenic host variants that may have influenced the COVID-19 outcome in Bulgarian patients.

Variants in the Late Cornified Envelope Gene Locus Are Associated With Elevated T-helper 17 Responses in Patients With Postinfectious Lyme Arthritis.

Postinfectious Lyme arthritis (LA) is associated with dysregulated immunity and autoreactive T- and B-cell responses in joints. Here we explored the role of host genetic variation in this outcome. The frequency of 253 702 single-nucleotide polymorphisms (SNPs) was determined in 147 patients with LA (87 with postinfectious LA and 60 with antibiotic-responsive LA), and for comparison in 90 patients with erythema migrans or the general population (n = 2504). Functional outcome of candidate SNPs was assessed by evaluating their impact on clinical outcome and on immune responses in blood and synovial fluid in patients with LA. Six SNPs associated with late cornified envelope (LCE3) genes were present at greater frequency in patients with postinfectious LA compared to those with antibiotic-responsive LA (70% vs 30%; odds ratio, 2; P < .01). These SNPs were associated with heightened levels of inflammatory Th17 cytokines in serum but lower levels of interleukin 27, a regulatory cytokine, implying that they may contribute to dysregulated Th17 immunity in blood. Moreover, in patients with postinfectious LA, the levels of these Th17 mediators correlated directly with autoantibody responses in synovial fluid, providing a possible link between LCE3 SNPs, maladaptive systemic Th17 immunity, and autoreactive responses in joints. Variation in the LCE3 locus, a known genetic risk factor in psoriasis and psoriatic arthritis, is associated with dysregulated systemic Th17 immunity and heightened autoantibody responses in joints. These findings underscore the importance of host genetic predisposition and systemic Th17 immunity in the pathogenesis of postinfectious (antibiotic-refractory) Lyme arthritis.

Association between TMPRSS2 rs2070788 polymorphism and COVID-19 severity: a case-control study in multiple cities of Iran.

Host genetic variations have been identified as potential influencers of COVID-19 infection. This study aimed to examine the association between transmembrane serine protease type 2 (TMPRSS2) rs2070788 single nucleotide polymorphism (SNP) and the prognosis of COVID-19 in Iranian populations. This case-control study was performed on 756 COVID-19 patients and 59 healthy individuals across Iran. Clinical data, blood samples, and the presence of the TMPRSS2 rs2070788: G>A SNP were determined using T-ARMS-PCR. Additionally, serum levels of tumor necrosis factor α (TNF-α), C-reactive protein (CRP), interleukin-6 (IL-6), and IL-1β were evaluated in the collected blood samples. No significant association was found between the genotypes and allele frequencies of TMPRSS2 rs2070788 SNP and susceptibility to or mortality from COVID-19 infection. However, we observed a substantial increase in IL-6 and CRP levels associated with the severity of COVID-19, while no such trend was observed for IL-1β and TNF-α. This study showed a considerable rise in TNF-α and IL-1β serum levels exclusively in COVID-19 patients with TT rs2070788 TMPRSS2 SNP genotype compared to healthy controls. In this study conducted across multiple cities in Iran, no significant association was found between the TMPRSS2 rs2070788 SNP genotypes and COVID-19 severity or mortality.

MMS19 and IFIH1 Host Genetic Variants Associate with SARS-CoV-2 Infection in Elderly Residents of Long-Term Care Facilities

MMS19 and IFIH1 Host Genetic Variants Associate with SARS-CoV-2 Infection in Elderly Residents of Long-Term Care Facilities

Host genetic variation and specialized metabolites from wheat leaves enriches for phyllosphere Pseudomonas spp. with enriched antibiotic resistomes.

Antibiotic resistance in plant-associated microbiomes poses significant risks for agricultural ecosystems and human health. Although accumulating evidence suggests a role for plant genotypes in shaping their microbiome, almost nothing is known about how the changes of plant genetic information affect the co-evolved plant microbiome carrying antibiotic resistance genes (ARGs). Here, we selected 16 wheat cultivars and experimentally explored the impact of host genetic variation on phyllosphere microbiome, ARGs, and metabolites. Our results demonstrated that host genetic variation significantly influenced the phyllosphere resistomes. Wheat genotypes exhibiting high phyllosphere ARGs were linked to elevated Pseudomonas populations, along with increased abundances of Pseudomonas aeruginosa biofilm formation genes. Further analysis of 350 Pseudomonas spp. genomes from diverse habitats at a global scale revealed that nearly all strains possess multiple ARGs, virulence factor genes (VFGs), and mobile genetic elements (MGEs) on their genomes, albeit with lower nucleotide diversity compared to other species. These findings suggested that the proliferation of Pseudomonas spp. in the phyllosphere significantly contributed to antibiotic resistance. We further observed direct links between the upregulated leaf metabolite DIMBOA-Glc, Pseudomonas spp., and enrichment of phyllosphere ARGs, which were corroborated by microcosm experiments demonstrating that DIMBOA-Glc significantly enhanced the relative abundance of Pseudomonas species. Overall, alterations in leaf metabolites resulting from genetic variation throughout plant evolution may drive the development of highly specialized microbial communities capable of enriching phyllosphere ARGs. This study enhances our understanding of how plants actively shape microbial communities and clarifies the impact of host genetic variation on the plant resistomes.

On the host specificity and genetic diversity of Iodamoeba bütschlii: Observations from short amplicon-based next-generation sequencing

Iodamoeba is a single-celled intestinal parasite, which is common in humans in certain parts of the world, and also in pigs. For the first time, we provide DNA-based evidence of goat, dromedary, fallow deer, and donkey as hosts of Iodamoeba and show that Iodamoeba-specific nucleotide sequences from these four hosts do not appear to overlap with those of humans, unlike those from pigs. We moreover show that similar strains of Iodamoeba can be found in Madagascar, Western Sahara, and Ecuador and that intra-sample diversity is typically extensive across even small fragments of DNA in both human and non-human hosts.

Genetic architecture of heritable leaf microbes.

Host-associated microbiomes are shaped by both their environment and host genetics, and often impact host performance. The scale of host genetic variation important to microbes is largely unknown yet fundamental to the community assembly of host-associated microbiomes, with implications for the eco-evolutionary dynamics of microbes and hosts. Using Ipomoea hederacea, ivyleaf morning glory, we generated matrilines differing in quantitative genetic variation and leaf shape, which is controlled by a single Mendelian locus. We then investigated the relative roles of Mendelian and quantitative genetic variation in structuring the leaf microbiome and how these two sources of genetic variation contributed to microbe heritability. We found that despite large effects of the environment, both Mendelian and quantitative genetic host variation contribute to microbe heritability and that the cumulative small effect genomic differences due to matriline explained as much or more microbial variation than a single large effect Mendelian locus. Furthermore, our results are the first to suggest that leaf shape itself contributes to variation in the abundances of some phyllosphere microbes.IMPORTANCEWe investigated how host genetic variation affects the assembly of Ipomoea hederacea's natural microbiome. We found that the genetic architecture of leaf-associated microbiomes involves both quantitative genetic variation and Mendelian traits, with similar contributions to microbe heritability. The existence of Mendelian and quantitative genetic variation for host-associated microbes means that plant evolution at the leaf shape locus or other quantitative genetic loci has the potential to shape microbial abundance and community composition.

Stratified analyses refine association between TLR7 rare variants and severe COVID-19

Despite extensive global research into genetic predisposition for severe COVID-19, knowledge on the role of rare host genetic variants and their relation to other risk factors remains limited. Here, 52 genes with prior etiological evidence were sequenced in 1,772 severe COVID-19 cases and 5,347 population-based controls from Spain/Italy. Rare deleterious TLR7 variants were present in 2.4% of young (<60 years) cases with no reported clinical risk factors (n=378), compared to 0.24% of controls (odds ratio (OR)=12.3, p=1.27x10-10). Incorporation of the results of either functional assays or protein modeling led to a pronounced increase in effect size (ORmax=46.5, p=1.74x10-15). Association signals for X-chromosomal TLR7 were also detected in the female-only subgroup, suggesting the existence of additional mechanisms beyond X-linked recessive inheritance in males. Additionally, supporting evidence was generated for a contribution to severe COVID-19 of the previously implicated genes IFNAR2, IFIH1 and TBK1. Our results refine the genetic contribution of rare TLR7 variants to severe COVID-19, and strengthen evidence for the etiological relevance of genes in the interferon signaling pathway.

Genome-wide association study identifies host genetic variants influencing oral microbiota diversity and metabolic health

The microbial communities of the oral cavity are important elements of oral and systemic health. With emerging evidence highlighting the heritability of oral bacterial microbiota, this study aimed to identify host genome variants that influence oral microbial traits. Using data from 16S rRNA gene amplicon sequencing, we performed genome-wide association studies with univariate and multivariate traits of the salivary microbiota from 610 unrelated adults from the Danish ADDITION-PRO cohort. We identified six single nucleotide polymorphisms (SNPs) in human genomes that showed associations with abundance of bacterial taxa at different taxonomical tiers (P < 5 × 10–8). Notably, SNP rs17793860 surpassed our study-wide significance threshold (P < 1.19 × 10–9). Additionally, rs4530093 was linked to bacterial beta diversity (P < 5 × 10–8). Out of these seven SNPs identified, six exerted effects on metabolic traits, including glycated hemoglobin A1c, triglyceride and high-density lipoprotein cholesterol levels, the risk of type 2 diabetes and stroke. Our findings highlight the impact of specific host SNPs on the composition and diversity of the oral bacterial community. Importantly, our results indicate an intricate interplay between host genetics, the oral microbiota, and metabolic health. We emphasize the need for integrative approaches considering genetic, microbial, and metabolic factors.

EQTLs identify regulatory networks and drivers of variation in the individual response to sepsis

Sepsis is a clinical syndrome of life-threatening organ dysfunction caused by a dysregulated response to infection, for which disease heterogeneity is a major obstacle to developing targeted treatments. We have previously identified gene-expression-based patient subgroups (sepsis response signatures [SRS]) informative for outcome and underlying pathophysiology. Here, we aimed to investigate the role of genetic variation in determining the host transcriptomic response and to delineate regulatory networks underlying SRS. Using genotyping and RNA-sequencing data on 638 adult sepsis patients, we report 16,049 independent expression (eQTLs) and 32 co-expression module (modQTLs) quantitative trait loci in this disease context. We identified significant interactions between SRS and genotype for 1,578 SNP-gene pairs and combined transcription factor (TF) binding site information (SNP2TFBS) and predicted regulon activity (DoRothEA) to identify candidate upstream regulators. Overall, these approaches identified putative mechanistic links between host genetic variation, cell subtypes, and the individual transcriptomic response to infection.

History of tuberculosis disease is associated with genetic regulatory variation in Peruvians.

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.