Levels of interpopulation differentiation among different functional classes of immunologically important genes

Contrasting patterns of selection and drift between two categories of immune genes in prairie-chickens

Immune-receptor genes of the adaptive immune system, such as the major histocompatibility complex (MHC), are involved in recognizing specific pathogens and are known to have high rates of adaptive evolution, presumably as a consequence of rapid co-evolution between hosts and pathogens. In contrast, many 'mediating' genes of the immune system do not interact directly with specific pathogens and are involved in signalling (e.g. cytokines) or controlling immune cell growth. As a consequence, we might expect stronger selection at immune-receptor than mediating genes, but these two types of genes have not been compared directly in wild populations. Here, we tested the hypothesis that selection differs between MHC (class I and II) and mediating genes by comparing levels of population differentiation across the range of greater prairie-chickens (Tympanuchus cupido). As predicted, there was stronger population differentiation and isolation by distance at immune receptor (MHC) than at either mediating genes or neutral microsatellites, suggesting a stronger role of local adaptation at the MHC. In contrast, mediating genes displayed weaker differentiation between populations than neutral microsatellites, consistent with selection favouring similar alleles across populations for mediating genes. In addition to selection, drift also had a stronger effect on immune receptor (MHC) than mediating genes as indicated by the stronger decline of MHC variation in relation to population size. This is the first study in the wild to show that the effects of selection and drift on immune genes vary across populations depending on their functional role.

Patterns of nucleotide and haplotype diversity at ICAM-1 across global human populations with varying levels of malaria exposure

Malaria is one of the strongest selective pressures in recent human evolution. African populations have been and continue to be at risk for malarial infections. However, few studies have re-sequenced malaria susceptibility loci across geographically and genetically diverse groups in Africa. We examined nucleotide diversity at Intercellular adhesion molecule-1 (ICAM-1), a malaria susceptibility candidate locus, in a number of human populations with a specific focus on diverse African ethnic groups. We used tests of neutrality to assess whether natural selection has impacted this locus and tested whether SNP variation at ICAM-1 is correlated with malaria endemicity. We observe differing patterns of nucleotide and haplotype variation in global populations and higher levels of diversity in Africa. Although we do not observe a deviation from neutrality based on the allele frequency distribution, we do observe several alleles at ICAM-1, including the ICAM-1 (Kilifi) allele, that are correlated with malaria endemicity. We show that the ICAM-1 (Kilifi) allele, which is common in Africa and Asia, exists on distinct haplotype backgrounds and is likely to have arisen more recently in Asia. Our results suggest that correlation analyses of allele frequencies and malaria endemicity may be useful for identifying candidate functional variants that play a role in malaria resistance and susceptibility.

Purifying selection modulates the estimates of population differentiation and confounds genome-wide comparisons across single-nucleotide polymorphisms

An improved understanding of the biological and numerical properties of measures of population differentiation across loci is becoming increasingly more important because of their growing use in analyzing genome-wide polymorphism data for detecting population structures, inferring the rates of migration, and identifying local adaptations. In a genome-wide analysis, we discovered that the estimates of population differentiation (e.g., F(ST), θ, and Jost's D) calculated for human single-nucleotide polymorphisms (SNPs) are strongly and positively correlated to the position-specific evolutionary rates measured from multispecies alignments. That is, genomic positions (loci) experiencing higher purifying selection (lower evolutionary rates) produce lower values for the degree of population differentiation than those evolving with faster rates. We show that this pattern is completely mediated by the negative effects of purifying selection on the minor allele frequency (MAF) at individual loci. Our results suggest that inferences and methods relying on the comparison of population differentiation estimates (F(ST), θ, and Jost's D) based on SNPs across genomic positions should be restricted to loci with similar MAFs and/or the rates of evolution in genome scale surveys.

Cytokine gene polymorphisms: methods of detection and biological significance

In recent years, several functional polymorphisms, particularly, SNPs have been identified in cytokines and their receptor genes that regulate levels of cytokine expression. These have been implicated as immune prognostic markers in diseases, including differential response to therapy and as biomarkers of graft outcome following organ and stem cell transplantation. Population distribution of cytokine gene polymorphisms (CGPs) reveals significant variations in allele frequencies in different ethnic groups and this might explain, to some extent, the observed differences in SNP associations with various diseases and immune-pathologies. A number of molecular methods are available for defining CGPs. These include PCR-SSP, AFLP, Taqman probe assays as well as sequencing based typing. Of these, the PCR-based sequence-specific primer based test is the most widely accepted technique. This chapter describes steps involved in this procedure along with sources for procuring essential reagents. An important aspect of CGP analyses is the correct interpretation of results particularly determination of their multilocus haplotypes.

Distribution of events of positive selection and population differentiation in a metabolic pathway: the case of asparagine N-glycosylation

Background

Asparagine N-Glycosylation is one of the most important forms of protein post-translational modification in eukaryotes. This metabolic pathway can be subdivided into two parts: an upstream sub-pathway required for achieving proper folding for most of the proteins synthesized in the secretory pathway, and a downstream sub-pathway required to give variability to trans-membrane proteins, and involved in adaptation to the environment and innate immunity. Here we analyze the nucleotide variability of the genes of this pathway in human populations, identifying which genes show greater population differentiation and which genes show signatures of recent positive selection. We also compare how these signals are distributed between the upstream and the downstream parts of the pathway, with the aim of exploring how forces of population differentiation and positive selection vary among genes involved in the same metabolic pathway but subject to different functional constraints.

Results

Our results show that genes in the downstream part of the pathway are more likely to show a signature of population differentiation, while events of positive selection are equally distributed among the two parts of the pathway. Moreover, events of positive selection are frequent on genes that are known to be at bifurcation points, and that are identified as being in key position by a network-level analysis such as MGAT3 and GCS1.

Conclusions

These findings indicate that the upstream part of the Asparagine N-Glycosylation pathway has lower diversity among populations, while the downstream part is freer to tolerate diversity among populations. Moreover, the distribution of signatures of population differentiation and positive selection can change between parts of a pathway, especially between parts that are exposed to different functional constraints. Our results support the hypothesis that genes involved in constitutive processes can be expected to show lower population differentiation, while genes involved in traits related to the environment should show higher variability. Taken together, this work broadens our knowledge on how events of population differentiation and of positive selection are distributed among different parts of a metabolic pathway.

Persistent Inflammatory Pathways Associated with Early Onset Myocardial Infarction in a Medicated Multiethnic Hawaiian Cohort

In spite of current standard therapies to target the major pathomechanisms in myocardial infarction (MI), inflammatory gene expression patterns have been consistently revealed in MI patients. In a multiethnic cohort, we aimed to identify MI-associated pathomechanisms that may be unresponsive to medical treatment to improve diagnosis and therapy. Gene expression profiles in whole blood were analyzed in medicated Asian, African American and Caucasian patients living in Hawaii with a history of early MI and age, ethnicity, risk factor and medication-matched controls. PANTHER ontological and Ingenuity Pathway analysis and functional evaluation of the consistently differentially expressed genes identified coordinated up-regulation of genes for inflammation (LGALS3, PTX3, ZBTB32, BCL2L1), T-cell activation (IL12RB1, VAV3, JAG1, CAMP), immune imbalance (IL-8, IL2RA, CCR7, AHNAK), and active atherosclerosis (NR1H4, BIN1, GSTT1, MARCO) that persist in MI patients in spite of concerted treatment efforts to control vascular pathology. Furthermore, significant ethnic differences appear to exist within the active disease mechanisms that need to be further investigated to identify key targets for effective medical intervention.

Haplotype structure and expression divergence at the Drosophila cellular immune gene eater

The protein Eater plays an important role in microbial recognition and defensive phagocytosis in Drosophila melanogaster. We sequenced multiple alleles of the eater gene from an African and a North American population of D. melanogaster and found signatures of a partial selective sweep in North America that is localized around the second intron. This pattern is consistent with local adaptation to novel selective pressures during range expansion out of Africa. The North American sample is divided into two predominant haplotype groups, and the putatively selected haplotype is associated with a significantly higher gene expression level, suggesting that gene regulation is a possible target of selection. The eater alleles contain from 22 to 40 repeat units that are characterized by the presence of a cysteine-rich NIM motif. NIM repeats in the structural stalk of the protein exhibit concerted evolution as a function of physical location in the repeat array. Several NIM repeats within eater have previously been implicated in binding to microbial ligands, a function which in principle might subject them to special evolutionary pressures. However, we find no evidence of elevated positive selection on these pathogen-interacting units. Our study presents an instance where gene expression rather than protein structure is thought to drive the adaptive evolution of a pathogen recognition molecule in the immune system.

Patterns of population differentiation of candidate genes for cardiovascular disease

Background

The basis for ethnic differences in cardiovascular disease (CVD) susceptibility is not fully understood. We investigated patterns of population differentiation (F_ST) of a set of genes in etiologic pathways of CVD among 3 ethnic groups: Yoruba in Nigeria (YRI), Utah residents with European ancestry (CEU), and Han Chinese (CHB) + Japanese (JPT). We identified 37 pathways implicated in CVD based on the PANTHER classification and 416 genes in these pathways were further studied; these genes belonged to 6 biological processes (apoptosis, blood circulation and gas exchange, blood clotting, homeostasis, immune response, and lipoprotein metabolism). Genotype data were obtained from the HapMap database.

Results

We calculated F_ST for 15,559 common SNPs (minor allele frequency ≥ 0.10 in at least one population) in genes that co-segregated among the populations, as well as an average-weighted F_ST for each gene. SNPs were classified as putatively functional (non-synonymous and untranslated regions) or non-functional (intronic and synonymous sites). Mean F_ST values for common putatively functional variants were significantly higher than F_ST values for nonfunctional variants. A significant variation in F_ST was also seen based on biological processes; the processes of 'apoptosis' and 'lipoprotein metabolism' showed an excess of genes with high F_ST. Thus, putative functional SNPs in genes in etiologic pathways for CVD show greater population differentiation than non-functional SNPs and a significant variance of F_ST values was noted among pairwise population comparisons for different biological processes.

Conclusion

These results suggest a possible basis for varying susceptibility to CVD among ethnic groups.

Frequency distribution of cytokine gene polymorphisms in the healthy North Indian population

The allelic and genotype frequencies corresponding to 22 single nucleotide polymorphisms (SNPs) in 13 cytokine genes interleukin (IL) 1-alpha (T/C -889), IL1-beta (C/T -511, T/C +3962), IL12 (C/A -1188), interferon-gamma (A/T UTR 5644), transforming growth factor-beta (C/T codon 10, G/C codon 25), tumour necrosis factor-alpha (G/A -308, G/A -238), IL2 (T/G -330, G/T +160), IL4 (T/G -1098, T/C -590, T/C -33), IL6 (G/C -174, G/A nt 565), IL10 (G/A -1082, C/T -819, C/A -592), IL1R (C/T pst11970), IL1RA (T/C mspa111100) and IL4RA (G/A +1902) were determined in 130 healthy North Indian subjects. All genomic typings were performed with polymerase chain reaction with sequence-specific assays. An analysis of the allelic and haplotype frequencies in the North Indian population showed a good fit with the Hardy-Weinberg equilibrium for most of the SNPs. The data can be used for anthropological comparisons, as well as for association studies with different diseases and for use in transplant situations involving acute and chronic rejection.