Host Genetic Architecture Research Articles

Genome-wide association study identifies genetic variants underlying footrot in Portuguese Merino sheep

BackgroundOvine footrot caused by Dichelobacter nodosus (D. nodosus) is a contagious disease with serious economic and welfare impacts in sheep production systems worldwide. A better understanding of the host genetic architecture regarding footrot resistance/susceptibility is crucial to develop disease control strategies that efficiently reduce infection and its severity. A genome-wide association study was performed using a customized SNP array (47,779 SNPs in total) to identify genetic variants associated to footrot resistance/susceptibility in two Portuguese native breeds, i.e. Merino Branco and Merino Preto, and a population of crossbred animals. A cohort of 1375 sheep sampled across 17 flocks, located in the Alentejo region (southern Portugal), was included in the analyses.ResultsPhenotypes were scored from 0 (healthy) to 5 (severe footrot) based on visual inspection of feet lesions, following the Modified Egerton System. Using a linear mixed model approach, three SNPs located on chromosome 24 reached genome-wide significance after a Bonferroni correction (p < 0.05). Additionally, six genome-wide suggestive SNPs were identified each on chromosomes 2, 4, 7, 8, 9 and 15. The annotation and KEGG pathway analyses showed that these SNPs are located within regions of candidate genes such as the nonsense mediated mRNA decay associated PI3K related kinase (SMG1) (chromosome 24) and the RALY RNA binding protein like (RALYL) (chromosome 9), both involved in immunity, and the heparan sulfate proteoglycan 2 (HSPG2) (chromosome 2) and the Thrombospodin 1 (THBS1) (chromosome 7) implicated in tissue repair and wound healing processes.ConclusionThis is the first attempt to identify molecular markers associated with footrot in Portuguese Merino sheep. These findings provide relevant information on a likely genetic association underlying footrot resistance/susceptibility and the potential candidate genes affecting this trait. Genetic selection strategies assisted on the information obtained from this study could enhance Merino sheep-breeding programs, in combination with farm management strategies, for a more effective and sustainable long-term solution for footrot control.

Association of Viral and Host Genetic Architecture with the Status of Neurocognitive Disorder in HIV-Infected Individuals.

The polymorphisms in host genes such as CCR5, CCR2, stromal derived factor (SDF), and MBL (mannose-binding lectin) as well as the viral nef gene have been shown to influence human immunodeficiency virus (HIV) infection, followed by the development of HIV-associated neurocognitive disorder (HAND). In this preliminary study with a limited number of samples, we have tried to associate the genetic polymorphism from the host and viral genetic factors with the neurocognitive status along with immuno-virological parameters. The total RNA was isolated from 10 unlinked plasma samples containing 5 samples from each group with and without HAND based on the International HIV Dementia Scale (IHDS) score <9.5 and >9.5, respectively. The CCR5, CCR2, SDF, MBL, and HIV nef genes were amplified and digested with restriction enzymes, except for the nef gene amplicon. Restrictions fragment length polymorphism (RFLP) was used to determine whether allelic variations were present in the digested host gene products, while sequencing was done for HIV nef amplicons without digestion. CCR5 delta 32 heterozygous variants were present in two samples from the HAND group. Three samples with HAND showed SDF-1 3' heterozygous allelic variant, while the MBL-2 gene presented with a homozygous mutant allele (D/D) in codon 52, heterozygous mutant allele (A/B) in codon 54, and codon 57 (A/C) for all samples except IHDS-2 irrespective of dementia status. Furthermore, amino acid alignment of Nef sequences confirmed the heterogeneity, while prediction of the human leukocyte antigen binding epitopes further explored its effect on functional motifs with variable binding efficiency such as epitopes GAFDLSFFL (aa 83) and LTFGWCFKL (aa 138) binding with HLA molecules at 60% and 80%, respectively. Thus, host genetics evidently influence predisposition to HIV infection and HAND. The genetic variability in the nef gene from both groups resulted in altering the functionality of specific domains and showing its impact on the progression of the disease, which needs to be explored.

The Genetic Architecture of Adaptation to Leaf and Root Bacterial Microbiota inArabidopsis thaliana

Understanding the role of the host genome in modulating microbiota variation is a need to shed light on the holobiont theory and overcome the current limits on the description of host-microbiota interactions at the genomic and molecular levels. However, the host genetic architecture structuring microbiota is only partly described in plants. In addition, most association genetic studies on microbiota are often carried out outside the native habitats where the host evolves and the identification of signatures of local adaptation on the candidate genes has been overlooked. To fill these gaps and dissect the genetic architecture driving adaptive plant-microbiota interactions, we adopted a genome-environment association (GEA) analysis on 141 whole-genome sequenced natural populations of Arabidopsis thaliana characterized in situ for their leaf and root bacterial communities in fall and spring, and a large range of nonmicrobial ecological factors (i.e., climate, soil, and plant communities). A much higher fraction of among-population microbiota variance was explained by the host genetics than by nonmicrobial ecological factors. Importantly, the relative importance of host genetics and nonmicrobial ecological factors in explaining the presence of particular operational taxonomic units (OTUs) differs between bacterial families and genera. In addition, the polygenic architecture of adaptation to bacterial communities was highly flexible between plant compartments and seasons. Relatedly, signatures of local adaptation were stronger on quantitative trait loci (QTLs) of the root microbiota in spring. Finally, plant immunity appears as a major source of adaptive genetic variation structuring bacterial assemblages in A. thaliana.

Integrating single-cell sequencing data with GWAS summary statistics reveals CD16+monocytes and memory CD8+T cells involved in severe COVID-19

BackgroundUnderstanding the host genetic architecture and viral immunity contributes to the development of effective vaccines and therapeutics for controlling the COVID-19 pandemic. Alterations of immune responses in peripheral blood mononuclear cells play a crucial role in the detrimental progression of COVID-19. However, the effects of host genetic factors on immune responses for severe COVID-19 remain largely unknown.MethodsWe constructed a computational framework to characterize the host genetics that influence immune cell subpopulations for severe COVID-19 by integrating GWAS summary statistics (N = 969,689 samples) with four independent scRNA-seq datasets containing healthy controls and patients with mild, moderate, and severe symptom (N = 606,534 cells). We collected 10 predefined gene sets including inflammatory and cytokine genes to calculate cell state score for evaluating the immunological features of individual immune cells.ResultsWe found that 34 risk genes were significantly associated with severe COVID-19, and the number of highly expressed genes increased with the severity of COVID-19. Three cell subtypes that are CD16+monocytes, megakaryocytes, and memory CD8+T cells were significantly enriched by COVID-19-related genetic association signals. Notably, three causal risk genes of CCR1, CXCR6, and ABO were highly expressed in these three cell types, respectively. CCR1+CD16+monocytes and ABO+ megakaryocytes with significantly up-regulated genes, including S100A12, S100A8, S100A9, and IFITM1, confer higher risk to the dysregulated immune response among severe patients. CXCR6+ memory CD8+ T cells exhibit a notable polyfunctionality including elevation of proliferation, migration, and chemotaxis. Moreover, we observed an increase in cell-cell interactions of both CCR1+ CD16+monocytes and CXCR6+ memory CD8+T cells in severe patients compared to normal controls among both PBMCs and lung tissues. The enhanced interactions of CXCR6+ memory CD8+T cells with epithelial cells facilitate the recruitment of this specific population of T cells to airways, promoting CD8+T cell-mediated immunity against COVID-19 infection.ConclusionsWe uncover a major genetics-modulated immunological shift between mild and severe infection, including an elevated expression of genetics-risk genes, increase in inflammatory cytokines, and of functional immune cell subsets aggravating disease severity, which provides novel insights into parsing the host genetic determinants that influence peripheral immune cells in severe COVID-19.

Whole-Genome Duplication and Host Genotype Affect Rhizosphere Microbial Communities.

ABSTRACTThe composition of microbial communities found in association with plants is influenced by host phenotype and genotype. However, the ways in which specific genetic architectures of host plants shape microbiomes are unknown. Genome duplication events are common in the evolutionary history of plants and influence many important plant traits, and thus, they may affect associated microbial communities. Using experimentally induced whole-genome duplication (WGD), we tested the effect of WGD on rhizosphere bacterial communities in Arabidopsis thaliana. We performed 16S rRNA amplicon sequencing to characterize differences between microbiomes associated with specific host genetic backgrounds (Columbia versus Landsberg) and ploidy levels (diploid versus tetraploid). We modeled relative abundances of bacterial taxa using a hierarchical Bayesian approach. We found that host genetic background and ploidy level affected rhizosphere community composition. We then tested to what extent microbiomes derived from a specific genetic background or ploidy level affected plant performance by inoculating sterile seedlings with microbial communities harvested from a prior generation. We found a negative effect of the tetraploid Columbia microbiome on growth of all four plant genetic backgrounds. These findings suggest an interplay between host genetic background and ploidy level and bacterial community assembly with potential ramifications for host fitness. Given the prevalence of ploidy-level variation in both wild and managed plant populations, the effects on microbiomes of this aspect of host genetic architecture could be a widespread driver of differences in plant microbiomes.IMPORTANCE Plants influence the composition of their associated microbial communities, yet the underlying host-associated genetic determinants are typically unknown. Genome duplication events are common in the evolutionary history of plants and affect many plant traits. Using Arabidopsis thaliana, we characterized how whole-genome duplication affected the composition of rhizosphere bacterial communities and how bacterial communities associated with two host plant genetic backgrounds and ploidy levels affected subsequent plant growth. We observed an interaction between ploidy level and genetic background that affected both bacterial community composition and function. This research reveals how genome duplication, a widespread genetic feature of both wild and crop plant species, influences bacterial assemblages and affects plant growth.

Genomic and Ancestral Variation Underlies the Severity of COVID-19 Clinical Manifestation in Individuals of European Descent.

The coronavirus disease (COVID-19) caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is characterized by a wide spectrum of clinical phenotypes ranging from asymptomatic to symptomatic with mild or moderate presentation and severe disease. COVID-19 susceptibility, severity and recovery have demonstrated high variability worldwide. Variances in the host genetic architecture may underlie the inter-individual and population-scale differences in COVID-19 presentation. We performed a genome-wide association analysis employing the genotyping data from AncestryDNA for COVID-19 patients of European descent and used asymptomatic subjects as the control group. We identified 621 genetic variants that were significantly distinct between asymptomatic and acutely symptomatic COVID-19 patients (multiple-testing corrected p-value < 0.001). These variants were found to be associated with pathways governing host immunity, such as interferon, interleukin and cytokine signalling, and known COVID-19 comorbidities, such as obesity and cholesterol metabolism. Further, our ancestry analysis revealed that the asymptomatic COVID-19 patients possess discernibly higher proportions of the Ancestral North Eurasian (ANE) and Eastern Hunter-Gatherer (EHG) ancestry, which was introduced to Europe through Bell Beaker culture (Yamnaya related) and lower fractions of Western Hunter-Gatherer (WHG) ancestry, while severely symptomatic patients have higher fractions of WHG and lower ANE/EHG ancestral components, thereby delineating the likely ancestral differences between the two groups.

Modulating the Microbiota as a Therapeutic Intervention for Type 2 Diabetes.

Mounting evidence suggested that the gut microbiota has a significant role in the metabolism and disease status of the host. In particular, Type 2 Diabetes (T2D), which has a complex etiology that includes obesity and chronic low-grade inflammation, is modulated by the gut microbiota and microbial metabolites. Current literature supports that unbalanced gut microbial composition (dysbiosis) is a risk factor for T2D. In this review, we critically summarize the recent findings regarding the role of gut microbiota in T2D. Beyond these associative studies, we focus on the causal relationship between microbiota and T2D established using fecal microbiota transplantation (FMT) or probiotic supplementation, and the potential underlying mechanisms such as byproducts of microbial metabolism. These microbial metabolites are small molecules that establish communication between microbiota and host cells. We critically summarize the associations between T2D and microbial metabolites such as short-chain fatty acids (SCFAs) and trimethylamine N-Oxide (TMAO). Additionally, we comment on how host genetic architecture and the epigenome influence the microbial composition and thus how the gut microbiota may explain part of the missing heritability of T2D found by GWAS analysis. We also discuss future directions in this field and how approaches such as FMT, prebiotics, and probiotics supplementation are being considered as potential therapeutics for T2D.

The landscape of host genetic factors involved in immune response to common viral infections

BackgroundHumans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities.MethodsWe conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort.ResultsSignals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRβ1 at positions 11, 13, 71, and 74 for Epstein-Barr virus (EBV), Varicella zoster virus (VZV), human herpesvirus 7, (HHV7), and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P < 5.0 × 10−8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, and CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P = 5.0 × 10−15 (MCV), NTN5: P = 1.1 × 10−9 (BKV), and P2RY13: P = 1.1 × 10−8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases, from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions.ConclusionsOur study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.

Genome-wide pleiotropy and shared biological pathways for resistance to bovine pathogens

Host genetic architecture is a major factor in resistance to pathogens and parasites. The collection and analysis of sufficient data on both disease resistance and host genetics has, however, been a major obstacle to dissection the genetics of resistance to single or multiple pathogens. A severe challenge in the estimation of heritabilities and genetic correlations from pedigree-based studies has been the confounding effects of the common environment shared among relatives which are difficult to model in pedigree analyses, especially for health traits with low incidence rates. To circumvent this problem we used genome-wide single-nucleotide polymorphism data and implemented the Genomic-Restricted Maximum Likelihood (G-REML) method to estimate the heritabilities and genetic correlations for resistance to 23 different infectious pathogens in calves and cows in populations undergoing natural pathogen challenge. Furthermore, we conducted gene-based analysis and generalized gene-set analysis to understand the biological background of resistance to infectious diseases. The results showed relatively higher heritabilities of resistance in calves than in cows and significant pleiotropy (both positive and negative) among some calf and cow resistance traits. We also found significant pleiotropy between resistance and performance in both calves and cows. Finally, we confirmed the role of the B-lymphocyte pathway as one of the most important biological pathways associated with resistance to all pathogens. These results both illustrate the potential power of these approaches to illuminate the genetics of pathogen resistance in cattle and provide foundational information for future genomic selection aimed at improving the overall production fitness of cattle.

Bayesian Diallel Analysis Reveals Mx1-Dependent and Mx1-Independent Effects on Response to Influenza A Virus in Mice.

Influenza A virus (IAV) is a respiratory pathogen that causes substantial morbidity and mortality during both seasonal and pandemic outbreaks. Infection outcomes in unexposed populations are affected by host genetics, but the host genetic architecture is not well understood. Here, we obtain a broad view of how heritable factors affect a mouse model of response to IAV infection using an 8 × 8 diallel of the eight inbred founder strains of the Collaborative Cross (CC). Expanding on a prior statistical framework for modeling treatment response in diallels, we explore how a range of heritable effects modify acute host response to IAV through 4 d postinfection. Heritable effects in aggregate explained ∼57% of the variance in IAV-induced weight loss. Much of this was attributable to a pattern of additive effects that became more prominent through day 4 postinfection and was consistent with previous reports of antiinfluenza myxovirus resistance 1 (Mx1) polymorphisms segregating between these strains; these additive effects largely recapitulated haplotype effects observed at the Mx1 locus in a previous study of the incipient CC, and are also replicated here in a CC recombinant intercross population. Genetic dominance of protective Mx1 haplotypes was observed to differ by subspecies of origin: relative to the domesticus null Mx1 allele, musculus acts dominantly whereas castaneus acts additively. After controlling for Mx1, heritable effects, though less distinct, accounted for ∼34% of the phenotypic variance. Implications for future mapping studies are discussed.

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

Infectious diseases have a profound impact on our health and many studies suggest that host genetics play a major role in the pathogenesis of most of them. We perform 23 genome-wide association studies for common infections and infection-associated procedures, including chickenpox, shingles, cold sores, mononucleosis, mumps, hepatitis B, plantar warts, positive tuberculosis test results, strep throat, scarlet fever, pneumonia, bacterial meningitis, yeast infections, urinary tract infections, tonsillectomy, childhood ear infections, myringotomy, measles, hepatitis A, rheumatic fever, common colds, rubella and chronic sinus infection, in over 200,000 individuals of European ancestry. We detect 59 genome-wide significant (P < 5 × 10−8) associations in genes with key roles in immunity and embryonic development. We apply fine-mapping analysis to dissect associations in the human leukocyte antigen region, which suggests important roles of specific amino acid polymorphisms in the antigen-binding clefts. Our findings provide an important step toward dissecting the host genetic architecture of response to common infections.

Host genetic architecture and the landscape of microbiome composition: humans weigh in.

Comparative analyses of the control of mammalian microbiomes by host genetic architecture reveal striking conserved features that have implications for the evolution of host–microbiome interactions.See related Research article: http://www.genomebiology.com/2015/16/1/191

Interferon-λ rs12979860 genotype and liver fibrosis in viral and non-viral chronic liver disease.

Tissue fibrosis is a core pathologic process that contributes to mortality in ~45% of the population and is likely to be influenced by the host genetic architecture. Here we demonstrate, using liver disease as a model, that a single-nucleotide polymorphism (rs12979860) in the intronic region of interferon-λ4 (IFNL4) is a strong predictor of fibrosis in an aetiology-independent manner. In a cohort of 4,172 patients, including 3,129 with chronic hepatitis C (CHC), 555 with chronic hepatitis B (CHB) and 488 with non-alcoholic fatty liver disease (NAFLD), those with rs12979860CC have greater hepatic inflammation and fibrosis. In CHC, those with rs12979860CC also have greater stage-constant and stage-specific fibrosis progression rates (P<0.0001 for all). The impact of rs12979860 genotypes on fibrosis is maximal in young females, especially those with HCV genotype 3. These findings establish rs12979860 genotype as a strong aetiology-independent predictor of tissue inflammation and fibrosis.

Evidence for Varied Aetiologies Regulating the Transmission of Prion Disease: Implications for Understanding the Heritable Basis of Prion Incubation Times

BackgroundTransmissible Spongiform Encephalopathies (TSEs) are a group of progressive fatal neurodegenerative disorders, triggered by abnormal folding of the endogenous prion protein molecule. The encoding gene is a major biological factor influencing the length of the asymptomatic period after infection. It remains unclear the extent to which the variation between quantitative trait loci (QTLs) reported in mouse models is due to methodological differences between approaches or genuine differences between traits. With this in mind, our approach to identifying genetic factors has sought to extend the linkage mapping approach traditionally applied, to a series of additional traits, while minimising methodological variability between them. Our approach allows estimations of heritability to be derived, as well as predictions to be made about possible existence of genetic overlap between the various traits.Methodology/Principal FindingsOur data indicate a surprising degree of heritability (up to 60%). Correlations between traits are also identified. A series of QTLs on chromosomes 1, 2, 3, 4, 6, 11 and 18 accompany our heritability estimates. However, only a locus on chromosome 11 has a general effect across all 4 models explored.Conclusions/SignificanceWe have achieved some success in detecting novel and pre-existing QTLs associated with incubation time. However, aside from the general effects described, the model-specific nature of the broader host genetic architecture has also been brought into clearer focus. This suggests that genetic overlap can only partially account for the general heritability of incubation time when factors, such as the nature of the TSE agent and the route of administration are considered. This point is highly relevant to vCJD (a potential threat to public health) where the route of primary importance is oral, while the QTLs being sought derive exclusively from studies of the ic route. Our results highlight the limitations of a single-model approach to QTL-mapping of TSEs.

Host genetic architecture in single and multiple infections

Hosts are often target to multiple simultaneous infections by genetically diverse parasite strains. The interaction among these strains and the interaction of each strain with the host was shown to have profound effects on the evolution of parasite traits. Host factors like genetic architecture of resistance have so far been largely neglected. To see whether genetic architecture differs between different kinds of infections we used joint scaling analysis to compare the genetic components of resistance in the red flour beetle Tribolium castaneum exposed to single and multiple strains of the microsporidian Nosema whitei. Our results indicate that additive, dominance and epistatic components were more important in single infections whereas maternal components play a decisive role in multiple infections. In detail, parameter estimates of additive, dominance and epistatic components correlated positively between single and multiple infections, whereas maternal components correlated negatively. These findings may suggest that specificity of host-parasite interactions are mediated by genetic and especially epistatic components whereas maternal effects constitute a more general form of resistance.