Crohn's disease loci are common targets of protozoa-driven selection

Association between Taxonomic Composition of Gut Microbiota and Host Single Nucleotide Polymorphisms in Crohn's Disease Patients from Russia

Crohn's disease (CD) is a chronic relapsing inflammatory bowel disease of unknown etiology. Genetic predisposition and dysbiotic gut microbiota are important factors in the pathogenesis of CD. In this study, we analyzed the taxonomic composition of the gut microbiota and genotypes of 24 single nucleotide polymorphisms (SNP) associated with the risk of CD. The studied cohorts included 96 CD patients and 24 healthy volunteers from Russia. Statistically significant differences were found in the allele frequencies for 8 SNPs and taxonomic composition of the gut microbiota in CD patients compared with controls. In addition, two types of gut microbiota communities were identified in CD patients. The main distinguishing driver of bacterial families for the first community type are Bacteroidaceae and unclassified members of the Clostridiales order, and the second type is characterized by increased abundance of Streptococcaceae and Enterobacteriaceae. Differences in the allele frequencies of the rs9858542 (BSN), rs3816769 (STAT3), and rs1793004 (NELL1) were also found between groups of CD patients with different types of microbiota communities. These findings confirm the complex multifactorial nature of CD.

Prioritizing autoimmunity risk variants for functional analyses by fine-mapping mutations under natural selection

Pathogen-driven selection shaped adaptive mutations in immunity genes, including those contributing to inflammatory disorders. Functional characterization of such adaptive variants can shed light on disease biology and past adaptations. This popular idea, however, was difficult to test due to challenges in pinpointing adaptive mutations in selection footprints. In this study, using a local-tree-based approach, we show that 28% of risk loci (153/535) in 21 inflammatory disorders bear footprints of moderate and weak selection, and part of them are population specific. Weak selection footprints allow partial fine-mapping, and we show that in 19% (29/153) of the risk loci under selection, candidate disease variants are hitchhikers, and only in 39% of cases they are likely selection targets. We predict function for a subset of these selected SNPs and highlight examples of antagonistic pleiotropy. We conclude by offering disease variants under selection that can be tested functionally using infectious agents and other stressors to decipher the poorly understood link between environmental stressors and genetic risk in inflammatory conditions.

Evidence for the loss and recovery of SLAMF9 during human evolution: implications on Dollo's law

Signaling lymphocyte activation molecule family member 9 (SLAMF9) is a cell surface protein of the CD2/SLAM family of leukocyte surface receptors. It is conserved throughout mammals and has roles in the initiation of inflammatory responses and regulation of plasmacytoid dendritic cell function. Through comparison of reference sequences encoding SLAMF9 in human, mouse, and primate sequences, we have determined that the SLAMF9 gene underwent successive mutation events, resulting in the loss of the protein and subsequent recovery of a less stable version. The mutations included a single base pair deletion in the second exon and a change in the splice acceptor site of that same exon. These changes would have had the effect of creating and later repairing a frameshift in the coding sequence. These events took place since the divergence of the human lineage from the chimpanzee-human last common ancestor and represent the first known case of the functional loss and recovery of a gene within the human lineage.

Identification of differentially expressed genes between the colon and ileum of patients with inflammatory bowel disease by gene co-expression analysis

OBJECTIVE

We aimed to identify differentially expressed genes (DEG) in patients with inflammatory bowel disease (IBD).

METHODS

RNA-seq data were obtained from the Array Express database. DEG were identified using the edgeR package. A co-expression network was constructed and key modules with the highest correlation with IBD inflammatory sites were identified for analysis. The Cytoscape MCODE plugin was used to identify key sub-modules of the protein-protein interaction (PPI) network. The genes in the sub-modules were considered hub genes, and functional enrichment analysis was performed. Furthermore, we constructed a drug-gene interaction network. Finally, we visualized the hub gene expression pattern between the colon and ileum of IBD using the ggpubr package and analyzed it using the Wilcoxon test.

RESULTS

DEG were identified between the colon and ileum of IBD patients. Based on the co-expression network, the green module had the highest correlation with IBD inflammatory sites. In total, 379 DEG in the green module were identified for the PPI network. Nineteen hub genes were differentially expressed between the colon and ileum. The drug-gene network identified these hub genes as potential drug targets.

CONCLUSION

Nineteen DEG were identified between the colon and ileum of IBD patients.

Transethnic associations among immune-mediated diseases and single-nucleotide polymorphisms of the aryl hydrocarbon response gene ARNT and the PTPN22 immune regulatory gene

BACKGROUND

Because immune responses are sensitive to environmental changes that drive selection of genetic variants, we hypothesized that polymorphisms of some xenobiotic response and immune response genes may be associated with specific types of immune-mediated diseases (IMD), while others may be associated with IMD as a larger category regardless of specific phenotype or ethnicity.

OBJECTIVE

To examine transethnic gene-IMD associations for single nucleotide polymorphism (SNP) frequencies of prototypic xenobiotic response genes-aryl hydrocarbon receptor (AHR), AHR nuclear translocator (ARNT), AHR repressor (AHRR) - and a prototypic immune response gene, protein tyrosine phosphatase, non-receptor type 22 (PTPN22), in subjects from the Environmental Polymorphisms Registry (EPR).

METHODS

Subjects (n = 3731) were genotyped for 14 SNPs associated with functional variants of the AHR, ARNT, AHRR, and PTPN22 genes, and their frequencies were compared among African Americans (n = 1562), Caucasians (n = 1838), and Hispanics (n = 331) with previously reported data. Of those genotyped, 2015 EPR subjects completed a Health and Exposure survey. SNPs were assessed via PLINK for associations with IMD, which included those with autoimmune diseases, allergic disorders, asthma, or idiopathic pulmonary fibrosis. Transethnic meta-analyses were performed using METAL and MANTRA approaches.

RESULTS

ARNT SNP rs11204735 was significantly associated with autoimmune disease by transethnic meta-analyses using METAL (odds ratio, OR [95% confidence interval] = 1.29 [1.08-1.55]) and MANTRA (ORs ranged from 1.29 to 1.30), whereas ARNT SNP rs1889740 showed a significant association with autoimmune disease by METAL (OR = 1.25 [1.06-1.47]). For Caucasian females, PTPN22 SNP rs2476601 was significantly associated with autoimmune disease by allelic association tests (OR = 1.99, [1.30-3.04]). In Caucasians and Caucasian males, PTPN22 SNP rs3811021 was significantly associated with IMD (OR = 1.39 [1.12-1.72] and 1.50 [1.12-2.02], respectively) and allergic disease (OR = 1.39 [1.12-1.71], and 1.62 [1.19-2.20], respectively). In the transethnic meta-analysis, PTPN22 SNP rs3811021 was significantly implicated in IMD by METAL (OR = 1.31 [1.10-1.56]), and both METAL and MANTRA suggested that rs3811021 was associated with IMD and allergic disease in males across all three ethnic groups (IMD METAL OR = 1.50 [1.15-1.95]; IMD MANTRA ORs ranged from 1.47 to 1.50; allergic disease METAL OR = 1.58 [1.20-2.08]; allergic disease MANTRA ORs ranged from 1.55 to 1.59).

CONCLUSIONS

Some xenobiotic and immune response gene polymorphisms were shown here, for the first time, to have associations across a broad spectrum of IMD and ethnicities. Our findings also suggest a role for ARNT in the development of autoimmune diseases, implicating environmental factors metabolized by this pathway in pathogenesis. Further studies are needed to confirm these data, assess the implications of these findings, define gene-environment interactions, and explore the mechanisms leading to these increasingly prevalent disorders.

Autoimmune diseases - connecting risk alleles with molecular traits of the immune system

Genome-wide strategies have driven the discovery of more than 300 susceptibility loci for autoimmune diseases. However, for almost all loci, understanding of the mechanisms leading to autoimmunity remains limited, and most variants that are likely to be causal are in non-coding regions of the genome. A critical next step will be to identify the in vivo and ex vivo immunophenotypes that are affected by risk variants. To do this, key cell types and cell states that are implicated in autoimmune diseases will need to be defined. Functional genomic annotations from these cell types and states can then be used to resolve candidate genes and causal variants. Together with longitudinal studies, this approach may yield pivotal insights into how autoimmunity is triggered.

Genetics of autoimmune diseases: insights from population genetics

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Autoimmune diseases (ADs) are a family of complex heterogeneous disorders with similar underlying mechanisms characterized by immune responses against self. Collectively, ADs are common, exhibit gender and ethnic disparities, and increasing incidence. As natural selection is an important influence on human genetic variation, and immune function genes are enriched for signals of positive selection, it is thought that the prevalence of AD risk alleles seen in different population is partially the result of differing selective pressures (for example, due to pathogens). With the advent of high-throughput technologies, new analytical methodologies and large-scale projects, evidence for the role of natural selection in contributing to the heritable component of ADs keeps growing. This review summarizes the genetic regions associated with susceptibility to different ADs and concomitant evidence for selection, including known agents of selection exerting selective pressure in these regions. Examples of specific adaptive variants with phenotypic effects are included as an evidence of natural selection increasing AD susceptibility. Many of the complexities of gene effects in different ADs can be explained by population genetics phenomena. Integrating AD susceptibility studies with population genetics to investigate how natural selection has contributed to genetic variation that influences disease risk will help to identify functional variants and elucidate biological mechanisms. As such, the study of population genetics in human population holds untapped potential for elucidating the genetic causes of human disease and more rapidly focusing to personalized medicine.

The contribution of natural selection to present-day susceptibility to chronic inflammatory and autoimmune disease

Chronic inflammatory and autoimmune diseases have been the focus of many genome-wide association studies (GWAS) because they represent a significant cause of illness and morbidity, and many are heritable. Almost a decade of GWAS studies suggests that the pathological inflammation associated with these diseases is controlled by a limited number of networked immune system genes. Chronic inflammatory and autoimmune diseases are enigmatic from an evolutionary perspective because they exert a negative affect on reproductive fitness. The persistence of these conditions may be partially explained by the important roles the implicated immune genes play in pathogen defense and other functions thought to be under strong natural selection in humans. The evolutionary reasons for chronic inflammatory and autoimmune disease persistence and uneven distribution across populations are the focus of this review.

Autoimmune effector memory T cells: the bad and the good

Immunological memory is a hallmark of adaptive immunity, a defense mechanism endowed to vertebrates during evolution. However, an autoimmune pathogenic role of memory lymphocytes is also emerging with accumulating evidence, despite reasonable skepticism on their existence in a chronic setting of autoimmune damage. It is conceivable that autoimmune memory would be particularly harmful since memory cells would constantly "remember" and attack the body's healthy tissues. It is even more detrimental given the resistance of memory T cells to immunomodulatory therapies. In this review, we focus on self-antigen-reactive CD(+) effector memory T (TEM) cells, surveying the evidence for the role of the T(EM) compartment in autoimmune pathogenesis. We will also discuss the role of T(EM) cells in chronic and acute infectious disease settings and how they compare to their counterparts in autoimmune diseases. With their long-lasting potency, the autoimmune T(EM) cells could also play a critical role in anti-tumor immunity, which may be largely based on their reactivity to self-antigens. Therefore, although autoimmune T(EM) cells are "bad" due to their role in relentless perpetration of tissue damage in autoimmune disease settings, they are unlikely a by-product of industrial development along the modern surge of autoimmune disease prevalence. Rather, they may be a product of evolution for their "good" in clearing damaged host cells in chronic infections and malignant cells in cancer settings.

Human gene copy number variation and infectious disease

Variability in the susceptibility to infectious disease and its clinical manifestation can be determined by variation in the environment and by genetic variation in the pathogen and the host. Despite several successes based on candidate gene studies, defining the host variation affecting infectious disease has not been as successful as for other multifactorial diseases. Both single nucleotide variation and copy number variation (CNV) of the host contribute to the host's susceptibility to infectious disease. In this review we focus on CNV, particularly on complex multiallelic CNV that is often not well characterised either directly by hybridisation methods or indirectly by analysis of genotypes and flanking single nucleotide variants. We summarise the well-known examples, such as α-globin deletion and susceptibility to severe malaria, as well as more recent controversies, such as the extensive CNV of the chemokine gene CCL3L1 and HIV infection. We discuss the potential biological mechanisms that could underly any genetic association and reflect on the extensive complexity and functional variation generated by a combination of CNV and sequence variation, as illustrated by the Fc gamma receptor genes FCGR3A, FCGR3B and FCGR2C. We also highlight some understudied areas that might prove fruitful areas for further research.

Human genome variability, natural selection and infectious diseases

The recent availability of large-scale sequencing DNA data allowed researchers to investigate how genomic variation is distributed among populations. While demographic factors explain genome-wide population genetic diversity levels, scans for signatures of natural selection pinpointed several regions under non-neutral evolution. Recent studies found an enrichment of immune-related genes subjected to natural selection, suggesting that pathogens and infectious diseases have imposed a strong selective pressure throughout human history. Pathogen-mediated selection often targeted regulatory sites of genes belonging to the same biological pathway. Results from these studies have the potential to identify mutations that modulate infection susceptibility by integrating a population genomic approach with molecular immunology data and large-scale functional annotations.

Update on primary sclerosing cholangitis genetics

PURPOSE OF REVIEW

The pathogenesis of primary sclerosing cholangitis (PSC) involves heritable factors. This review summarizes the recent genetic studies and discusses the implications of identified risk loci.

RECENT FINDINGS

A total of 16 PSC susceptibility loci have been identified in genome-wide association studies and related study designs. At least 33 additional loci are involved in what is increasingly acknowledged to represent a general pool of genetic risk loci for immune-mediated diseases. One important group of genes is part of well characterized immune pathways (e.g. interleukin 2 signaling), whereas for other loci the relationship to PSC pathophysiology is less evident. Importantly, the loci collectively account for only 7.3% of overall PSC liability, thus pointing to a large contribution from environmental factors to PSC development. The individual PSC risk genes cannot be interpreted within a simple cause-effect model used for monogenic traits, but need to be explored for their individual biological correlates, preferably in a disease context. To some extent, as exemplified for the human leukocyte antigen and FUT2 associations, genetic findings may guide the discovery of interacting and co-occuring environmental susceptibility factors.

SUMMARY

Multiple PSC susceptibility loci are now available for exploration in experimental model systems and patient-centered research.

Natural selection and infectious disease in human populations

Infectious pathogens are among the strongest selective forces that shape the human genome. Migrations and cultural changes in the past 100,000 years exposed populations to dangerous new pathogens.

Host genetics influences susceptibility to infectious disease. Evolutionary adaptations for resistance and symbiosis may underlie common immune-mediated diseases.

Signatures of selection and methods to detect them vary with the age, geographical spread and virulence of the pathogen.

A history of selection on a trait adds power to association studies by driving the emergence of common alleles of strong effect. Combining selection and association metrics can further increase power.

Genome-wide association studies (GWASs) of susceptibility to pathogens that are moderately old (1,000–50,000 years ago), geographically limited in history and exerted strong positive selective pressure will have the most power if GWASs can be done in the historically affected population.

An understanding of host–pathogen interactions can inform the development of new therapies for both infectious diseases and common immune-mediated diseases.

The impact of various infectious agents on human survival and reproduction over thousands of years has exerted selective pressure on numerous regions of the human genome. This Review describes how such signatures of selection can be detected and integrated with data from complementary approaches, such as genome-wide association studies, to provide biological insights into host–pathogen interactions.

The ancient biological 'arms race' between microbial pathogens and humans has shaped genetic variation in modern populations, and this has important implications for the growing field of medical genomics. As humans migrated throughout the world, populations encountered distinct pathogens, and natural selection increased the prevalence of alleles that are advantageous in the new ecosystems in both host and pathogens. This ancient history now influences human infectious disease susceptibility and microbiome homeostasis, and contributes to common diseases that show geographical disparities, such as autoimmune and metabolic disorders. Using new high-throughput technologies, analytical methods and expanding public data resources, the investigation of natural selection is leading to new insights into the function and dysfunction of human biology.

Gene genealogies for genetic association mapping, with application to Crohn's disease

A gene genealogy describes relationships among haplotypes sampled from a population. Knowledge of the gene genealogy for a set of haplotypes is useful for estimation of population genetic parameters and it also has potential application in finding disease-predisposing genetic variants. As the true gene genealogy is unknown, Markov chain Monte Carlo (MCMC) approaches have been used to sample genealogies conditional on data at multiple genetic markers. We previously implemented an MCMC algorithm to sample from an approximation to the distribution of the gene genealogy conditional on haplotype data. Our approach samples ancestral trees, recombination and mutation rates at a genomic focal point. In this work, we describe how our sampler can be used to find disease-predisposing genetic variants in samples of cases and controls. We use a tree-based association statistic that quantifies the degree to which case haplotypes are more closely related to each other around the focal point than control haplotypes, without relying on a disease model. As the ancestral tree is a latent variable, so is the tree-based association statistic. We show how the sampler can be used to estimate the posterior distribution of the latent test statistic and corresponding latent p-values, which together comprise a fuzzy p-value. We illustrate the approach on a publicly-available dataset from a study of Crohn's disease that consists of genotypes at multiple SNP markers in a small genomic region. We estimate the posterior distribution of the tree-based association statistic and the recombination rate at multiple focal points in the region. Reassuringly, the posterior mean recombination rates estimated at the different focal points are consistent with previously published estimates. The tree-based association approach finds multiple sub-regions where the case haplotypes are more genetically related than the control haplotypes, and that there may be one or multiple disease-predisposing loci.

Disease consequences of human adaptation

Adaptive evolution has provided us with a unique set of characteristics that define us as humans, including morphological, physiological and cellular changes. Yet, natural selection provides no assurances that adaptation is without human health consequences; advantageous mutations will increase in frequency so long as there is a net gain in fitness. As such, the current incidence of human disease can depend on previous adaptations. Here, I review genome-wide and gene-specific studies in which adaptive evolution has played a role in shaping human genetic disease. In addition to the disease consequences of adaptive phenotypes, such as bipedal locomotion and resistance to certain pathogens, I review evidence that adaptive mutations have influenced the frequency of linked disease alleles through genetic hitchhiking. Taken together, the links between human adaptation and disease highlight the importance of their combined influence on functional variation within the human genome and offer opportunities to discover and characterize such variation.

Identifying Darwinian selection acting on different human APOL1 variants among diverse African populations

Disease susceptibility can arise as a consequence of adaptation to infectious disease. Recent findings have suggested that higher rates of chronic kidney disease (CKD) in individuals with recent African ancestry might be attributed to two risk alleles (G1 and G2) at the serum-resistance-associated (SRA)-interacting-domain-encoding region of APOL1. These two alleles appear to have arisen adaptively, possibly as a result of their protective effects against human African trypanosomiasis (HAT), or African sleeping sickness. In order to explore the distribution of potential functional variation at APOL1, we studied nucleotide variation in 187 individuals across ten geographically and genetically diverse African ethnic groups with exposure to two Trypanosoma brucei subspecies that cause HAT. We observed unusually high levels of nonsynonymous polymorphism in the regions encoding the functional domains that are required for lysing parasites. Whereas allele frequencies of G2 were similar across all populations (3%-8%), the G1 allele was only common in the Yoruba (39%). Additionally, we identified a haplotype (termed G3) that contains a nonsynonymous change at the membrane-addressing-domain-encoding region of APOL1 and is present in all populations except for the Yoruba. Analyses of long-range patterns of linkage disequilibrium indicate evidence of recent selection acting on the G3 haplotype in Fulani from Cameroon. Our results indicate that the G1 and G2 variants in APOL1 are geographically restricted and that there might be other functional variants that could play a role in HAT resistance and CKD risk in African populations.