Positive Selection on a Human-Specific Transcription Factor Binding Site Regulating IL4 Expression

Interleukin-4 Gene Polymorphisms in Romanian Patients with Inflammatory Bowel Diseases: Association with Disease Risk and Clinical Features

1.

INTRODUCTION

Multiple cytokines have been studied for their role in the propagation of the inflammatory process related to inflammatory bowel diseases (IBD), but the role of interleukin-4 remains controversial. The aim of this study was to evaluate the role of two IL-4 gene single nucleotide polymorphisms (SNPs) in disease susceptibility and phenotypic expression. 2.

MATERIALS AND METHODS

A group of 160 patients with IBD (86CD/74UC) and 160 healthy controls were genotyped for IL-4 rs2243250/-590C/T and rs2070874/-34C/T using real-time polymerase chain reaction with TaqMan assay. 3.

RESULTS

The analysis of IBD patients and controls revealed a significantly reduced frequency of the minor allele T of both SNPs in CD patients (p = 0.03, OR 0.55 and p = 0.02, OR 0.52) and for the entire IBD group (p = 0.01, OR 0.57 and p = 0.01, OR 0.55). Haplotype analysis identified the most frequent haplotype (rs2243250/rs2070874 CC) associated with a high risk for developing IBD (either UC or CD) (p = 0.003). IBD patients with extraintestinal manifestations had significantly increased frequency of the minor alleles T. We also found an association between the presence of allele C of rs2070874 and response to antiTNF treatment. 4.

CONCLUSIONS

This is the first study to investigate the IL-4 gene's relation to IBD susceptibility conducted in Romania. Both SNPs were found to be associated with disease susceptibility and phenotypic features, such as extraintestinal manifestations and response to antiTNF agents.

Frequency distribution of cytokine and associated transcription factor single nucleotide polymorphisms in Zimbabweans: Impact on schistosome infection and cytokine levels

Cytokines mediate T-helper (TH) responses that are crucial for determining the course of infection and disease. The expression of cytokines is regulated by transcription factors (TFs). Here we present the frequencies of single nucleotide polymorphisms (SNPs) in cytokine and TF genes in a Zimbabwean population, and further relate SNPs to susceptibility to schistosomiasis and cytokine levels. Individuals (N = 850) were genotyped for SNPs across the cytokines IL4, IL10, IL13, IL33, and IFNG, and their TFs STAT4, STAT5A/B, STAT6, GATA3, FOXP3, and TBX21 to determine allele frequencies. Circulatory levels of systemic and parasite-specific IL-4, IL-5, IL-10, IL-13, and IFNγ were quantified via enzyme-linked immunosorbent assay. Schistosoma haematobium infection was determined by enumerating parasite eggs excreted in urine by microscopy. SNP allele frequencies were related to infection status by case-control analysis and logistic regression, and egg burdens and systemic and parasite-specific cytokine levels by analysis of variance and linear regression. Novel findings were i) IL4 rs2070874*T’s association with protection from schistosomiasis, as carriage of ≥1 allele gave an odds ratio of infection of 0.597 (95% CIs, 0.421–0.848, p = 0.0021) and IFNG rs2069727*G’s association with susceptibility to schistosomiasis as carriage of ≥1 allele gave an odds ratio of infection of 1.692 (1.229–2.33, p = 0.0013). Neither IL4 rs2070874*T nor IFNG rs2069727*G were significantly associated with cytokine levels. This study found TH2-upregulating SNPs were more frequent among the Zimbabwean sample compared to African and European populations, highlighting the value of immunogenetic studies of African populations in the context of infectious diseases and other conditions, including allergic and atopic disease. In addition, the identification of novel infection-associated alleles in both TH1- and TH2-associated genes highlights the role of both in regulating and controlling responses to Schistosoma.

Genetic Variants of Interleukin-4 in Romanian Patients with Idiopathic Nephrotic Syndrome

Background and objectives: One of the most frequent glomerular diseases in the pediatric population is represented by the idiopathic nephrotic syndrome (INS). The exact mechanisms mediating the disease are still unknown, but several genetic factors have been studied for possible implications. Cytokines are considered to play a pivotal role in mediating INS disease progression, interleukin-4 (IL-4) exhibiting particular interest. The objective of this research project was to investigate the association between two IL-4 gene single-nucleotide polymorphisms (SNPs) and INS susceptibility as well as response to steroid therapy, in a group of Romanian children. Materials and Methods: In total, 75 patients with INS and 160 healthy controls of Romanian origin were genotyped for IL-4 rs2243250/−590C/T and rs2070874/−34C/T using real-time polymerase chain reaction. Association tests were performed using the DeFinetti program and Plink 1.07 software and p-values < 0.05 were considered statistically significant. Results: The analysis of INS patients and controls revealed a similar genotype distribution of the studied SNPs. The minor T alleles were less frequent in the INS group, but not statistically significant (p = 0.1, OR = 0.68 and p = 0.2, OR = 0.74). Regarding the response to steroids, a low frequency of 590*T allele in steroid-resistant patients (7.7%), compared with steroid-sensitive patients (14%) and controls (17.5%), was obtained, but the difference did not reach the statistical significance threshold. The same result was obtained for −34C/T SNP. Conclusions: This is the first study examining the relationship between the IL-4 gene and INS susceptibility conducted in a European population, and particularly in Romania. The investigated SNPs were found to not be associated with disease susceptibility or response to the steroid treatment of pediatric INS.

Bridging the gap between theory and practice in elucidating modular gene regulatory sequence organisation within genomes

Changes to promoter regions probably have been responsible for many morphological evolutionary transitions, especially in animals. This idea is becoming testable, as data from genome projects amass and enable bioinformaticians to conduct comparative sequence analyses and test for correlations between genotypic similarities or differences and phenotypic likeness or disparity. Although such practical pursuits have initiated some theoretical considerations, a conceptual framework for understanding promoter region evolution, potentially effecting morphological evolution, is only starting to emerge, predominantly resulting from computational research. We contribute to this framework by specifying three big problems for promoter region research; reviewing computational research on promoter region evolution; and exemplifying a topic for future promoter region research - module evolution.

Association of cytokine Th2 gene polymorphisms with autoimmune thyroid diseases in Tunisian population

Autoimmune thyroid diseases (AITD) including Graves' disease (GD) and Hashimoto's thyroiditis (HT) are complex genetic diseases. Th2 cytokines act on the development of AITD. This study was conducted on Tunisian patients with AITD to investigate the association of Th2 cytokine gene polymorphisms and haplotype combination with GD or HT risk. A total of 156 controls, 160 patients with HT and 88 patients with GD were genotyped for IL-4 rs2243250, IL-5 rs2069812, IL-6 rs1800796 and IL-13 rs1800925 polymorphisms by PCR-RFLP. The AITD risk was assessed by a logistic regression analysis using the SNP stats statistical program. False-positive report probability (FPRP) was estimated to evaluate significant findings. IL-13 rs1800925 was associated with GD, after adjustment for age and gender, in codominant, dominant and allele genetic models (p = .0072; p = .0018; p = .012, respectively). Significant association of the IL-6 rs1800796C/G genotype with GD was also detected (p = .025). Furthermore, increased risk of HT was still found for IL-13 rs1800925T allele (p = .039, OR = 1.39) and for IL-4 rs2243250T/T genotype both in codominant (p = .033, OR = 2.59) and recessive (p = .011, OR = 2.73) models after adjustment for age and gender. Interestingly, haplotype analysis performed on the IL-4, IL-5 and IL-13 genes revealed a high risk of HT with CTT haplotype (p = .008, OR = 2.12). However, the CCT haplotype is a protective factor (OR = 0.36). Patients carrying the CT haplotype with only one minor allele had a moderate risk of HT (OR = 1.56). The FPRP analysis showed that the association of IL-13 rs1800925 polymorphism with GD and HT and the association of CTT haplotype with HT were noteworthy. In conclusion, the IL-4, IL-5, IL-6 and IL-13 polymorphism may play a role in susceptibility to GD and HT in the Tunisian population. Furthermore, gene-gene interaction between the IL-4, IL-5 and IL-13 significantly increases the risk of AITD. Further studies with larger numbers of individuals are needed to confirm the results.

Preliminary analysis of single-nucleotide polymorphisms in IL-10, IL-4, and IL-4Rα genes and profile of circulating cytokines in patients with gastric Cancer

Background

Gastric Cancer is highly prevalent and deadly worldwide. In Colombia, it is the most lethal form of cancer. Some single-nucleotide polymorphisms in IL-10, IL-4, and IL-4Rα genes have been associated with an anti-inflammatory environment and a Th2 profile in detriment of the antitumor Th1 response. This research sought to detect single-nucleotide polymorphisms in promoter sequences, like − 1082 (G/A), − 592 (C/A), and − 819 (C/T), as well as − 590 (C/T) of the IL-10 and IL-4 genes, respectively; in addition to the IL-4Rα mutation variants, Ile50Val and Q576R, together with circulating levels of IL-4, TNF-α, IL-10, and IFN-γ in patients with gastric carcinoma in Cúcuta, Colombia.

Methods

In a cross-sectional study, 17 patients and 30 healthy individuals were genotyped for the six polymorphisms mentioned through PCR-RFLP of DNA obtained from peripheral blood cells and serum samples were analyzed by sandwich ELISA to quantify cytokines. Statistical difference between groups was determined along with the association between the presence of polymorphisms and the risk of gastric cancer, as well as the mortality in patients, using Mann-Whitney U test and logistic regression analysis, respectively.

Results

An association between the − 1082 (G/A) and the risk of gastric cancer was found (OR = 7.58, range 0.77–74.06, P = 0.08). Furthermore, patients had a significant increase in IL-4 serum levels (P < 0.01) compared to healthy individuals, both variables showed a higher estimated risk of mortality in patients, although without statistical association (P > 0.05).

Conclusion

We infer that two possible biomarkers (one immunological and one genetic) could be considered in association with gastric cancer in our population, which should be confirmed by subsequent studies involving a greater number of individuals.

Role of genetic factors and ethnicity on the multiplicity of Plasmodium falciparum infection in children with asymptomatic malaria in Yaoundé, Cameroon

In this cross-sectional study, we investigated host genetic factors and ethnic variation in circulating Plasmodium falciparum merozoite surface protein 2 (msp-2) clones among children with asymptomatic malaria. Isolates from seventy two asymptomatic malaria children were used for genotyping block 3 of msp-2 gene by nested polymerase chain reaction (PCR). Sickle cell trait and glucose-6-phosphate dehydrogenase (G6PD) deficiency were analysed by restriction fragment length polymorphism of DNA products from PCR targeting codons 6 and 68 of the beta-globin (HBB) and G6PD genes respectively. ABO blood group was typed by agglutination method. A total of forty two msp-2 genotypes (20 for 3D7 and 22 for FC27) were detected for an average (standard error of mean) multiplicity of infection (MOI) of 2.45 (0.16). The MOI was statistically the same among the five identified ethnic groups (P = 0.83). The overall prevalence of sickle cell trait and G6PD deficiency were 12.50 % and 22.22 % respectively. MOI was similar between children with Hb AA and Hb AS genotypes (P = 0.42). MOI was significantly high among children with a mutant G6PD genotype (P = 0.017). MOI was significantly higher in blood group O than group A (P = 0.03). Our findings show that although ethnicity and sickle cell trait have no association with MOI, the association was observed with G6PD genotype and ABO group. The results suggest the need for extension and expansion of the current study in order to investigate the mechanisms involved.

Positive selection of IL-33 in adaptive immunity of domestic Chinese goats

The identification of the candidate genes that play key role in phenotypic variation in livestock populations can provide new information about evolution and positive selection. IL‐33 (71954) (Interleukin) gene is associated with the increased nematode resistance in small ruminants; however, the role of IL‐33 for the genetic control of different diseases in Chinese goat breeds is poorly described in scientific literature. Therefore, the current investigation was performed for the better understanding of the molecular evolution and the positive selection of single‐nucleotide polymorphism in IL‐33 gene. Fixation Index (F_ST)‐based method was used for the outlier loci determination and found that IL‐33 was present in outlier area with the provisional combined allocation of mean heterozygosity and F_ST. Positively selected IL‐33 gene was significantly, that is, p(Simul F_ST < sample F_ST = 0.98*) present in corresponding positive selection area. Hence, our study provided novel information about the nucleotide variations in IL‐33 gene and found to be nonsynonymous which may helpful for the genetic control of diseases by enhancing the immune system in local Chinese goat breeds as well as in other analyzed vertebrate species.

Resolving the etiology of atopic disorders by using genetic analysis of racial ancestry

Atopic dermatitis (AD), food allergy, allergic rhinitis, and asthma are common atopic disorders of complex etiology. The frequently observed atopic march from early AD to asthma, allergic rhinitis, or both later in life and the extensive comorbidity of atopic disorders suggest common causal mechanisms in addition to distinct ones. Indeed, both disease-specific and shared genomic regions exist for atopic disorders. Their prevalence also varies among races; for example, AD and asthma have a higher prevalence in African Americans when compared with European Americans. Whether this disparity stems from true genetic or race-specific environmental risk factors or both is unknown. Thus far, the majority of the genetic studies on atopic diseases have used populations of European ancestry, limiting their generalizability. Large-cohort initiatives and new analytic methods, such as admixture mapping, are currently being used to address this knowledge gap. Here we discuss the unique and shared genetic risk factors for atopic disorders in the context of ancestry variations and the promise of high-throughput "-omics"-based systems biology approach in providing greater insight to deconstruct their genetic and nongenetic etiologies. Future research will also focus on deep phenotyping and genotyping of diverse racial ancestry, gene-environment, and gene-gene interactions.

Genetic signature of strong recent positive selection at interleukin-32 gene in goat

OBJECTIVE

Identification of the candidate genes that play key roles in phenotypic variations can provide new information about evolution and positive selection. Interleukin (IL)-32 is involved in many biological processes, however, its role for the immune response against various diseases in mammals is poorly understood. Therefore, the current investigation was performed for the better understanding of the molecular evolution and the positive selection of single nucleotide polymorphisms in IL-32 gene.

METHODS

By using fixation index (F_ST ) based method, IL-32 (9375) gene was found to be outlier and under significant positive selection with the provisional combined allocation of mean heterozygosity and F_ST . Using nucleotide sequences of 11 mammalian species from National Center for Biotechnology Information database, the evolutionary selection of IL-32 gene was determined using Maximum likelihood model method, through four models (M1a, M2a, M7, and M8) in Codeml program of phylogenetic analysis by maximum liklihood.

RESULTS

IL-32 is detected under positive selection using the F_ST simulations method. The phylogenetic tree revealed that goat IL-32 was in close resemblance with sheep IL-32. The coding nucleotide sequences were compared among 11 species and it was found that the goat IL-32 gene shared identity with sheep (96.54%), bison (91.97%), camel (58.39%), cat (56.59%), buffalo (56.50%), human (56.13%), dog (50.97%), horse (54.04%), and rabbit (53.41%) respectively.

CONCLUSION

This study provides evidence for IL-32 gene as under significant positive selection in goat.

Decoding the non-coding genome: elucidating genetic risk outside the coding genome

Current evidence emerging from genome-wide association studies indicates that the genetic underpinnings of complex traits are likely attributable to genetic variation that changes gene expression, rather than (or in combination with) variation that changes protein-coding sequences. This is particularly compelling with respect to psychiatric disorders, as genetic changes in regulatory regions may result in differential transcriptional responses to developmental cues and environmental/psychosocial stressors. Until recently, however, the link between transcriptional regulation and psychiatric genetic risk has been understudied. Multiple obstacles have contributed to the paucity of research in this area, including challenges in identifying the positions of remote (distal from the promoter) regulatory elements (e.g. enhancers) and their target genes and the underrepresentation of neural cell types and brain tissues in epigenome projects - the availability of high-quality brain tissues for epigenetic and transcriptome profiling, particularly for the adolescent and developing brain, has been limited. Further challenges have arisen in the prediction and testing of the functional impact of DNA variation with respect to multiple aspects of transcriptional control, including regulatory-element interaction (e.g. between enhancers and promoters), transcription factor binding and DNA methylation. Further, the brain has uncommon DNA-methylation marks with unique genomic distributions not found in other tissues - current evidence suggests the involvement of non-CG methylation and 5-hydroxymethylation in neurodevelopmental processes but much remains unknown. We review here knowledge gaps as well as both technological and resource obstacles that will need to be overcome in order to elucidate the involvement of brain-relevant gene-regulatory variants in genetic risk for psychiatric disorders.

Pleiotropic Effects of Immune Responses Explain Variation in the Prevalence of Fibroproliferative Diseases

Many diseases are differentially distributed among human populations. Differential selection on genetic variants in ancestral environments that coincidentally predispose to disease can be an underlying cause of these unequal prevalence patterns. Selected genes may be pleiotropic, affecting multiple phenotypes and resulting in more than one disease or trait. Patterns of pleiotropy may be helpful in understanding the underlying causes of an array of conditions in a population. For example, several fibroproliferative diseases are more prevalent and severe in populations of sub-Saharan ancestry. We propose that this disparity is due to selection for an enhanced Th2 response that confers resistance to helminthic infections, and concurrently increases susceptibility to fibrosis due to the profibrotic action of Th2 cytokines. Many studies on selection of Th2-related genes for host resistance to helminths have been reported, but the pleiotropic impact of this selection on the distribution of fibrotic disorders has not been explicitly investigated. We discuss the disproportionate occurrence of fibroproliferative diseases in individuals of African ancestry and provide evidence that adaptation of the immune system has shaped the genetic structure of these human populations in ways that alter the distribution of multiple fibroproliferative diseases.

Distribution of cytokine gene single nucleotide polymorphisms among a multi-ethnic Iranian population

Background:

Cytokine gene single nucleotide polymorphisms (SNPs) are widely used to study susceptibility to complex diseases and as a tool for anthropological studies.

Materials and Methods:

To investigate cytokine SNPs in an Iranian multi-ethnic population, we have investigated 10 interleukin (IL) SNPs (IL-1β (C-511T, T-31C), IL-2 (G-384T), IL-4 (C-590T), IL-6 (G-174C), IL-8 (T-251A), IL-10 (G-1082A, C-819T, C-592A) and tumor necrosis factor-alpha (TNF-α) (G-308A) in 415 Iranian subjects comprising of 6 different ethnicities. Allelic and genotypic frequencies as well as Hardy-Weinberg equilibrium (HWE) were calculated by PyPop software. Population genetic indices including observed heterozygosity (H_o), expected heterozygosity (H_e), fixation index (F_IS), the effective number of alleles (N_e) and polymorphism information content (PIC) were derived using Popgene 32 software. Multidimensional scaling (MDS) was constructed using Reynold's genetic distance obtained from the frequencies of cytokine gene polymorphism.

Results:

Genotypic distributions were consistent with the HWE assumptions, except for 3 loci (IL-4-590, IL-8-251 and IL-10-819) in Fars and 4 loci (IL-4-590, IL-6-174, IL-10-1082 and TNF-α-308) in Turks. Pairwise assessment of allelic frequencies, detected differences at the IL-4-590 locus in Gilakis versus Kurds (P = 0.028) and Lurs (P = 0.022). Mazanis and Gilakis displayed the highest (H_o= 0.50 ± 0.24) and lowest (H_o= 0.34 ± 0.16) mean observed heterozygosity, respectively.

Conclusions:

MDS analysis of our study population, in comparison with others, revealed that Iranian ethnicities except Kurds and Mazanis were tightly located within a single cluster with closest genetic affinity to Europeans.

Humoral immune responses against the malaria vaccine candidate antigen Plasmodium vivax AMA-1 and IL-4 gene polymorphisms in individuals living in an endemic area of the Brazilian Amazon

BACKGROUND

Several studies have recently demonstrated that the immune responses against malaria is governed by different factors, including the genetic components of the host. The IL-4 gene appears to be a strong candidate factor because of its role in the regulation of the Th2 response. The present study investigated the role of IL-4 polymorphisms in the development of IgG antibodies against PvAMA-1 and the IL-4 levels in individuals infected with Plasmodium vivax in a malaria endemic area in the Brazilian Amazon.

METHODS

The study sample included 83 patients who were diagnosed with P. vivax infection using thick smear and confirmed by nested-PCR. The IL-4 -590C>T and IL-4 -33C>T polymorphisms were genotyped by PCR-RFLP, and the intron 3 VNTR was genotyped by PCR. A standardised ELISA protocol was used to measure the total IgG against PvAMA-1. The cytokine/chemokine levels were measured using a Milliplex multiplex assay (Millipore). All of the subjects were genotyped with 48 ancestry informative markers to determine the proportions of African, European and Amerindian ancestry using STRUCTURE software.

RESULTS

Of the 83 patients, 60 (73%) produced IgG antibodies against PvAMA-1. A significant decrease in the percentage of respondents was observed among the primo-infected individuals. No significant differences were observed in the frequencies of genotypes and haplotypes among individuals who were positive or negative for IgG antibodies against PvAMA-1. Furthermore, no significant correlation was observed between the IL-4 polymorphisms, antibody levels, IL-4 levels, and parasitemia.

CONCLUSIONS

This study indicated that the polymorphisms identified in the IL-4 gene are not likely to play a role in the regulation of the antibody response against PvAMA-1 and IL-4 production in vivax malaria.

Analysis of genetic variants in the IL4 promoter and VNTR loci in Indian patients with Visceral Leishmaniasis

Visceral Leishmaniasis (VL) is the most severest form of Leishmaniasis and resistance to infection is mediated by cellular immune responses. Interleukin 4 (IL-4) orchestrates of Th2 and Th1 immune responses during infections. In this study, we aimed to investigate possible association between three functional IL-4 polymorphisms -590C/T (rs2243250), -34C/T (rs2070874) and 70bp VNTR (rs79071878 in intron3) with VL in an Indian cohort comprising of 197 VL patients and 193 healthy controls. The three investigated IL-4 polymorphisms were in strong linkage disequilibrium. The investigated IL-4 alleles, genotypes and the reconstructed haplotypes were not significantly distributed between the VL patients and healthy controls. Our study signifies no possible association of functional IL-4 polymorphisms with Indian VL and postulate other vital genes involved in the IL-4 pathway may provide genetic clues to elucidate of IL-4 regulation and immune-pathogenesis during VL.

Genetic polymorphisms and asthma: findings from a case-control study in the Madeira island population

BACKGROUND

Asthma is a complex disease influenced by multiple genetic and environmental factors. While Madeira has the highest prevalence of asthma in Portugal (14.6%), the effect of both genetic and environmental factors in this population has never been assessed. We categorized 98 asthma patients according to the Global Initiative for Asthma (GINA) guidelines, established their sensitization profile, and measured their forced expiratory volume in 1 second (FEV1) and forced vital capacity (FVC) indexes. Selected single nucleotide polymorphisms (SNPs) were analysed as potential markers for asthma susceptibility and severity in the interleukin 4 (IL4), interleukin 13 (IL13), beta-2-adrenergic receptor (ADRB2), a disintegrin and metalloprotease 33 (ADAM33), gasdermin-like (GSDML) and the signal transducer and activator of transcription 6 (STAT6) genes comparatively to a population reference set.

RESULTS

Although mites are the major source of allergic sensitization, no significant difference was found amongst asthma severity categories. IL4-590*CT/TT and IL4-RP2*253183/183183 were found to predict the risk (2-fold) and severity (3 to 4-fold) of asthma and were associated with a lower FEV1 index. ADRB2-c.16*AG is a risk factor (3.5-fold), while genotype GSDML-236*TT was protective (4-fold) for moderate-severe asthma. ADAM33-V4*C was associated to asthma and mild asthma by the transmission disequilibrium test (TDT). Finally, ADAM33-V4*CC and STAT6-21*TT were associated with higher sensitization (mean wheal size ≥10 mm) to house dust (1.4-fold) and storage mite (7.8-fold).

CONCLUSION

In Madeira, IL4-590C/T, IL4-RP2 253/183, GSDML-236C/T and ADAM33-V4C/G SNPs are important risk factors for asthma susceptibility and severity, with implications for asthma healthcare management.

Natural variation in abiotic stress responsive gene expression and local adaptation to climate in Arabidopsis thaliana

Gene expression varies widely in natural populations, yet the proximate and ultimate causes of this variation are poorly known. Understanding how variation in gene expression affects abiotic stress tolerance, fitness, and adaptation is central to the field of evolutionary genetics. We tested the hypothesis that genes with natural genetic variation in their expression responses to abiotic stress are likely to be involved in local adaptation to climate in Arabidopsis thaliana. Specifically, we compared genes with consistent expression responses to environmental stress (expression stress responsive, "eSR") to genes with genetically variable responses to abiotic stress (expression genotype-by-environment interaction, "eGEI"). We found that on average genes that exhibited eGEI in response to drought or cold had greater polymorphism in promoter regions and stronger associations with climate than those of eSR genes or genomic controls. We also found that transcription factor binding sites known to respond to environmental stressors, especially abscisic acid responsive elements, showed significantly higher polymorphism in drought eGEI genes in comparison to eSR genes. By contrast, eSR genes tended to exhibit relatively greater pairwise haplotype sharing, lower promoter diversity, and fewer nonsynonymous polymorphisms, suggesting purifying selection or selective sweeps. Our results indicate that cis-regulatory evolution and genetic variation in stress responsive gene expression may be important mechanisms of local adaptation to climatic selective gradients.

Interleukin-4 and interleukin-13 pathway genetics affect disease susceptibility, serum immunoglobulin E levels, and gene expression in asthma

BACKGROUND

Asthma is a common immune disorder characterized by increased IgE levels. The interleukin (IL)-4 and IL-13 pathway is central for IgE regulation, and previous studies have reported many genetic variants of IL-4/IL-13 signaling in relation to asthma, but few have focused on the gene-to-gene interactions that are likely to contribute to disease complexity.

OBJECTIVE

To assess the combined effects of 7 functional single-nucleotide polymorphisms (SNPs) on asthma susceptibility, total serum IgE levels, and gene expression in children.

METHODS

Seven SNPs (rs2243250, rs1800925, rs1805010, rs324011, rs2251746, rs2494262, and rs2427837) were genotyped children with asthma (n = 500) and a control group (n = 523), and total serum IgE levels and gene expressions were measured in children with asthma.

RESULTS

Children with asthma had a likelier possibility of carrying more risk genotypes. Mean IgE levels increased from the minimum of 71.07 KU/L in children with no tested polymorphisms to a maximum of 901.7 KU/L in children carrying 7 risk genotypes. Gene expression analysis showed that patients with 4 SNPs (rs2243250, rs1800925, rs1805010, and rs3224011) had higher expression levels of IL-4, IL-13, and STAT6. Moreover, serum IgE level generally correlated well with IL-4 (r = 0.236, P = .011) and IL-13 (r = 0.211, P = .021) expressions; IL-4 expression correlated positively with IL-13 (r = 0.962, P = .000) and STAT6 (r = 0.190, P = .022) expressions, and STAT6 expression correlated with IL-4RA expression (r = 0.904, P = .000).

CONCLUSION

These data suggest that combinations of multiple SNPs might magnify the impact on disease risk. Only a combined analysis of the variants in the IL-4/IL-13 pathway could show the functional interplay of multiple genes in asthma.

Population genetics of cis-regulatory sequences that operate during embryonic development in the sea urchin Strongylocentrotus purpuratus

Despite the fact that noncoding sequences comprise a substantial fraction of functional sites within all genomes, the evolutionary mechanisms that operate on genetic variation within regulatory elements remain poorly understood. In this study, we examine the population genetics of the core, upstream cis-regulatory regions of eight genes (AN, CyIIa, CyIIIa, Endo16, FoxB, HE, SM30 a, and SM50) that function during the early development of the purple sea urchin, Strongylocentrotus purpuratus. Quantitative and qualitative measures of segregating variation are not conspicuously different between cis-regulatory and closely linked "proxy neutral" noncoding regions containing no known functional sites. Length and compound mutations are common in noncoding sequences; conventional descriptive statistics ignore such mutations, under-representing true genetic variation by approximately 28% for these loci in this population. Patterns of variation in the cis-regulatory regions of six of the genes examined (CyIIa, CyIIIa, Endo16, FoxB, AN, and HE) are consistent with directional selection. Genetic variation within annotated transcription factor binding sites is comparable to, and frequently greater than, that of surrounding sequences. Comparisons of two paralog pairs (CyIIa/CyIIIa and AN/HE) suggest that distinct evolutionary processes have operated on their cis-regulatory regions following gene duplication. Together, these analyses provide a detailed view of the evolutionary mechanisms operating on noncoding sequences within a natural population, and underscore how little is known about how these processes operate on cis-regulatory sequences.

IL-4 haplotype -590T, -34T and intron-3 VNTR R2 is associated with reduced malaria risk among ancestral indian tribal populations

Background

Interleukin 4 (IL-4) is an anti-inflammatory cytokine, which regulates balance between T_H1 and T_H2 immune response, immunoglobulin class switching and humoral immunity. Polymorphisms in this gene have been reported to affect the risk of infectious and autoimmune diseases.

Methods

We have analyzed three regulatory IL-4 polymorphisms; -590C>T, -34C>T and 70 bp intron-3 VNTR, in 4216 individuals; including: (1) 430 ethnically matched case-control groups (173 severe malaria, 101 mild malaria and 156 asymptomatic); (2) 3452 individuals from 76 linguistically and geographically distinct endogamous populations of India, and (3) 334 individuals with different ancestry from outside India (84 Brazilian, 104 Syrian, and 146 Vietnamese).

Results

The -590T, -34T and intron-3 VNTR R2 alleles were found to be associated with reduced malaria risk (P<0.001 for -590C>T and -34C>T, and P = 0.003 for VNTR). These three alleles were in strong LD (r²>0.75) and the TTR2 (-590T, -34T and intron-3 VNTR R2) haplotype appeared to be a susceptibility factor for malaria (P = 0.009, OR = 0.552, 95% CI = 0.356 –0.854). Allele and genotype frequencies differ significantly between caste, nomadic, tribe and ancestral tribal populations (ATP). The distribution of protective haplotype TTR2 was found to be significant (χ²₃ = 182.95, p-value <0.001), which is highest in ATP (40.5%); intermediate in tribes (33%); and lowest in caste (17.8%) and nomadic (21.6%).

Conclusions

Our study suggests that the IL-4 polymorphisms regulate host susceptibility to malaria and disease progression. TTR2 haplotype, which gives protection against malaria, is high among ATPs. Since they inhabited in isolation and mainly practice hunter-gatherer lifestyles and exposed to various parasites, IL-4 TTR2 haplotype might be under positive selection.

Regulation of human formyl peptide receptor 1 synthesis: role of single nucleotide polymorphisms, transcription factors, and inflammatory mediators

The gene encoding the human formyl peptide receptor 1 (FPR1) is heterogeneous, containing numerous single nucleotide polymorphisms (SNPs). Here, we examine the effect of these SNPs on gene transcription and protein translation. We also identify gene promoter sequences and putative FPR1 transcription factors. To test the effect of codon bias and codon pair bias on FPR1 expression, four FPR1 genetic variants were expressed in human myeloid U937 cells fused to a reporter gene encoding firefly luciferase. No significant differences in luciferase activity were detected, suggesting that the translational regulation and protein stability of FPR1 are modulated by factors other than the SNP codon bias and the variant amino acid properties. Deletion and mutagenesis analysis of the FPR1 promoter showed that a CCAAT box is not required for gene transcription. A −88/41 promoter construct resulted in the strongest transcriptional activity, whereas a −72/41 construct showed large reduction in activity. The region between −88 and −72 contains a consensus binding site for the transcription factor PU.1. Mutagenesis of this site caused significant reduction in reporter gene expression. The PU.1 binding was confirmed in vivo by chromatin immunoprecipitation, and the binding to nucleotides −84 to −76 (TTCCTATTT) was confirmed in vitro by an electrophoretic mobility shift assay. Thus, similar to many other myeloid genes, FPR1 promoter activity requires PU.1. Two single nucleotide polymorphisms at −56 and −54 did not significantly affect FPR1 gene expression, despite differences in binding of transcription factor IRF1 in vitro. Inflammatory mediators such as interferon-γ, tumor necrosis factor-α, and lipopolysaccharide did not increase FPR1 promoter activity in myeloid cells, whereas differentiation induced by DMSO and retinoic acid enhanced the activity. This implies that the expression of FPR1 in myeloid cells is developmentally regulated, and that the differentiated cells are equipped for immediate response to microbial infections.

Evidence that purifying selection acts on promoter sequences

We tested whether functionally important sites in bacterial, yeast, and animal promoters are more conserved than their neighbors. We found that substitutions are predominantly seen in less important sites and that those that occurred tended to have less impact on gene expression than possible alternatives. These results suggest that purifying selection operates on promoter sequences.

Associations of IL-4, IL-4R, and IL-13 gene polymorphisms in coal workers' pneumoconiosis in China: a case-control study

BACKGROUND

The IL-4, IL-4 receptor (IL4R), and IL-13 genes are crucial immune factors and may influence the course of various diseases. In the present study, we investigated the association between the potential functional polymorphisms in IL-4, IL-4R, and IL-13 and coal workers' pneumoconiosis (CWP) risk in a Chinese population.

METHODS

Six polymorphisms (C-590T in IL-4, Ile50Val, Ser478Pro, and Gln551Arg in IL-4R, C-1055T and Arg130Gln in IL-13) were genotyped and analyzed in a case-control study of 556 CWP and 541 control subjects.

RESULTS

Our results revealed that the IL-4 CT/CC genotypes were associated with a significantly decreased risk of CWP (odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.58-0.95), compared with the TT genotype, particularly among subgroups of age <65 years (OR = 0.68, 95%CI = 0.46-0.99) and dust exposure years ≥26 years (OR = 0.69, 95%CI = 0.50-0.94). Moreover, the polymorphism was significantly associated with risk of CWP patients with stage I. In addition, a combined effect was observed in a dose-dependent manner with increasing numbers of risk variant alleles (P(trend) = 0.023), and individuals with 11-12 risk alleles had a 47% higher risk of CWP than those with 0-8 risk alleles (OR = 1.47, 95% CI = 1.05-2.05).

CONCLUSIONS

Our results suggest that the IL-4 C-590T polymorphism is involved in the etiology of CWP and susceptibility to this disease. Larger studies are warranted to validate our findings.

Evolution of IL4 and pathogen antagonism

Interleukin-4 (IL4) is a pleiotropic cytokine involved in host protection from gastrointestinal nematodes. Here, we review the structure, function, and evolutionary history of IL4. Cumulative evidence indicates that over 100 million years of eutherian mammalian evolution, IL4 has experienced multiple episodes of positive selection. We argue that IL4 may have evolved in conflict with pathogen-derived antagonists, and therefore diversified to escape antagonism while being constrained to maintain binding to its cellular receptors. Selective pressure driving IL4 diversification may have arisen from ancient episodes of conflict with parasitic worm-derived IL4 antagonists. Descendants of such antagonists may still equip the armamentarium of contemporary gastrointestinal nematodes.

Diversification at transcription factor binding sites within a species and the implications for environmental adaptation

Evolution of new cellular functions can be achieved both by changes in protein coding sequences and by alteration of expression patterns. Variation of expression may lead to changes in cellular function with relatively little change in genomic sequence. We therefore hypothesize that one of the first signals of functional divergence should be evolution of transcription factor-binding sites (TFBSs). This adaptation should be detectable as substantial variation in the TFBSs of alleles. New data sets allow the first analyses of intraspecies variation from large number of whole-genome sequences. Using data from the Saccharomyces Genome Resequencing Project, we have analyzed variation in TFBSs. We find a large degree of variation both between these closely related strains and between pairs of duplicated genes. There is a correlation between changes in promoter regions and changes in coding sequences, indicating a coupling of changes in expression and function. We show that 1) the types genes with diverged promoters vary between strains from different environments and 2) that patterns of divergence in promoters consistent with positive selection are detectable in alleles between strains and on duplicate promoters. This variation is likely to reflect adaptation to each strain's natural environment. We conclude that, even within a species, we detect signs of selection acting on promoter regions that may act to alter expression patterns. These changes may indicate functional innovation in multiple genes and across the whole genome. Change in function could represent adaptation to the environment and be a precursor to speciation.

Genetic variation in Th1/Th2 pathway genes and risk of non-Hodgkin lymphoma: a pooled analysis of three population-based case-control studies

The balance between T-helper 1 (Th1) and T-helper 2 (Th2) activity is critical in lymphoid cell development and differentiation. Immune dysfunction underlies lymphomagenesis, so an alteration in the regulation of key Th1/Th2 cytokines may lead to the development of non-Hodgkin lymphoma (NHL). To study the impact of polymorphisms in Th1/Th2 cytokines on NHL risk, we analyzed 145 tag single nucleotide polymorphisms (SNPs) in 17 Th1/Th2 cytokine and related genes in three population-based case-control studies (1946 cases and 1808 controls). Logistic regression was used to compute odds ratios (OR) for NHL and four major NHL subtypes in relation to tag SNP genotypes and haplotypes. A gene-based analysis adjusting for the number of tag SNPs genotyped in each gene showed significant associations with risk of NHL combined and one or more NHL subtypes for Th1 (IL12A and IL12RB1) and Th2 (IL4, IL10RB, and IL18) genes. The strongest association was for rs485497 in IL12A, which plays a central role in bridging the cellular and humoral pathways of innate resistance and antigen-specific adaptive immune responses (allele risk OR= 1·17; P(trend)= 0·00099). This SNP was also associated specifically with risk of follicular lymphoma (allele risk OR= 1·26; P(trend)= 0·0012). These findings suggest that genetic variation in Th1/Th2 cytokine genes may contribute to lymphomagenesis.

Genetic regulation of parasite infection: empirical evidence of the functional significance of an IL4 gene SNP on nematode infections in wild primates

BACKGROUND

Susceptibility to parasite infection affects fitness-related processes, such as mate choice and survival, yet its genetic regulation remains poorly understood. Interleukin-4 (IL4) plays a central role in the humoral immune defence against nematode parasite infections, inducing IgE switch and regulation of worm expulsion from the intestines. The evolutionary and functional significance of single nucleotide polymorphisms (SNPs) in IL4-genes is known, yet empirical information on the effect of IL4 SNPs on gastro-intestinal infections is lacking. Using samples from a population of wild red-fronted lemurs (Eulemur fulvus rufus, Primates: Lemuridae), from western Madagascar, we explored the association of IL4-gene promoter polymorphisms with nematode infections and investigated a possible functional role of the IL4 polymorphism on male reproductive success.

RESULTS

Using sequence analyses of lemur DNA we detected a new SNP in the IL4 gene promoter area. Carriers of the genotype T/T showed higher nematode infection intensities than individuals of genotypes C/T and C/C. Genetic population analyses using data from more than 10 years, suggested higher reproductive success of T/T males than expected.

CONCLUSIONS

Our results suggest a regulatory effect of an IL4 gene promoter polymorphism on the intensity of parasite infections in a natural population of red-fronted lemurs, with a seemingly disadvantageous genotype represented in low frequencies. Long-term population analyses, however, point in the direction of a negative frequency-dependent association, giving a fitness advantage to the rare genotype. Due to low frequencies of the genotype in question conclusive evidence of a functional role of IL4 polymorphism cannot be drawn here; still, we suggest the use of IL4 polymorphism as a new molecular tool for quick assessment of individual genetic constitution with regard to nematode infection intensities, contributing to a better understanding of the actual components of the immune response that mediate protection against gastro-intestinal parasites.

Genetic variation in Th1/Th2 pathway genes and risk of non-Hodgkin lymphoma: a pooled analysis of three population-based case-control studies

The balance between T-helper 1 (Th1) and T-helper 2 (Th2) activity is critical in lymphoid cell development and differentiation. Immune dysfunction underlies lymphomagenesis, so an alteration in the regulation of key Th1/Th2 cytokines may lead to the development of non-Hodgkin lymphoma (NHL). To study the impact of polymorphisms in Th1/Th2 cytokines on NHL risk, we analyzed 145 tag single nucleotide polymorphisms (SNPs) in 17 Th1/Th2 cytokine and related genes in three population-based case-control studies (1946 cases and 1808 controls). Logistic regression was used to compute odds ratios (OR) for NHL and four major NHL subtypes in relation to tag SNP genotypes and haplotypes. A gene-based analysis adjusting for the number of tag SNPs genotyped in each gene showed significant associations with risk of NHL combined and one or more NHL subtypes for Th1 (IL12A and IL12RB1) and Th2 (IL4, IL10RB, and IL18) genes. The strongest association was for rs485497 in IL12A, which plays a central role in bridging the cellular and humoral pathways of innate resistance and antigen-specific adaptive immune responses (allele risk OR= 1·17; P(trend)= 0·00099). This SNP was also associated specifically with risk of follicular lymphoma (allele risk OR= 1·26; P(trend)= 0·0012). These findings suggest that genetic variation in Th1/Th2 cytokine genes may contribute to lymphomagenesis.

Role of inflammation-related gene polymorphisms in patients with central retinal vein occlusion

OBJECTIVE

Central retinal vein occlusion (CRVO) is a vision-threatening disease, primarily occurring among patients aged more than 60 years. Several risk factors, including arterial hypertension and diabetes mellitus, have been identified. Compression of the central retinal vein by an atherosclerotic retinal artery at the lamina cribrosa also has been implicated in the pathogenesis of the disease. Functional gene polymorphisms of cytokines or chemokines previously shown to affect atherogenesis or hemostasis are potential risk factors for CRVO. The present study investigates a hypothesized association between inflammation-related gene polymorphisms and the presence of CRVO in a relatively large cohort of patients.

DESIGN

Case-control study.

PARTICIPANTS

The study group consisted of 315 patients with CRVO and 335 control subjects.

METHODS

Determination of genotypes was done by 5' exonuclease assay (TaqMan).

MAIN OUTCOME MEASURES

Genotypes of interleukin (IL)1β -511C>T, IL1 receptor antagonist (IL1RN) 1018T>C, IL4 -584C>T, IL6 -174G>C, IL10 -592C>A, IL18 183A>G, tumor necrosis factor (TNF)-α -308G>A, monocyte chemoattractant protein (MCP)-1/CCL2 -2518A>G, IL8 -251A>T, and RANTES (CCL5) -403G>A polymorphisms.

RESULTS

Genotype distributions and allele frequencies of the investigated gene polymorphisms did not significantly differ between both groups (P>0.05). Arterial hypertension, diabetes mellitus, and cigarette smoking were significantly more frequent in patients with CRVO than among control subjects (arterial hypertension: 67.0% vs. 52.2%, P<0.001; diabetes mellitus: 16.8% vs. 6.3%, P<0.001, cigarette smoking: 32.1% vs. 23.6%, P = 0.02). In a logistic regression analysis, the presence of arterial hypertension was associated with an odds ratio (OR) of 1.75 (95% confidence interval [CI], 1.26-2.44) in those with CRVO, whereas an OR of 2.52 (95% CI, 1.46-4.35) was found in those with diabetes mellitus. A history of cigarette smoking was associated with an OR of 1.57 (95% CI, 1.09 - 2.25) for CRVO.

CONCLUSIONS

Our data suggest that the investigated inflammation-related gene polymorphisms are unlikely major risk factors for CRVO.

FINANCIAL DISCLOSURE(S)

The author(s) have no proprietary or commercial interest in any materials discussed in this article.

Diversifying selection and functional analysis of interleukin-4 suggests antagonism-driven evolution at receptor-binding interfaces

Background

Interleukin-4 (IL4) is a secreted immunoregulatory cytokine critically involved in host protection from parasitic helminths [1]. Reasoning that helminths may have evolved mechanisms to antagonize IL4 to maximize their dispersal, we explored mammalian IL4 evolution.

Results

This analysis revealed evidence of diversifying selection at 15 residues, clustered in epitopes responsible for IL4 binding to its Type I and Type II receptors. Such a striking signature of selective pressure suggested either recurrent episodes of pathogen antagonism or ligand/receptor co-evolution. To test the latter possibility, we performed detailed functional analysis of IL4 allotypes expressed by Mus musculus musculus and Mus musculus castaneus, which happen to differ at 5 residues (including three at positively selected sites) in and adjacent to the site 1 epitope that binds the IL4Rα subunit shared by the Type I and Type II IL4 receptors. We show that this intra-species variation affects the ability of IL4 neither to bind IL4 receptor alpha (IL4Rα) nor to signal biological responses through its Type I receptor.

Conclusions

Our results -- reminiscent of clustered positively selected sites revealing functionally important residues at host-virus interaction interfaces -- are consistent with IL4 having evolved to avoid recurrent pathogen antagonism, while maintaining the capacity to bind and signal through its cognate receptor. This work exposes what may be a general feature of evolutionary conflicts fought by pathogen antagonists at host protein-protein interaction interfaces involved in immune signaling: the emergence of receptor-binding ligand epitopes capable of buffering amino acid variation.

Genome-wide scans for footprints of natural selection

Detecting recent selected ‘genomic footprints’ applies directly to the discovery of disease genes and in the imputation of the formative events that molded modern population genetic structure. The imprints of historic selection/adaptation episodes left in human and animal genomes allow one to interpret modern and ancestral gene origins and modifications. Current approaches to reveal selected regions applied in genome-wide selection scans (GWSSs) fall into eight principal categories: (I) phylogenetic footprinting, (II) detecting increased rates of functional mutations, (III) evaluating divergence versus polymorphism, (IV) detecting extended segments of linkage disequilibrium, (V) evaluating local reduction in genetic variation, (VI) detecting changes in the shape of the frequency distribution (spectrum) of genetic variation, (VII) assessing differentiating between populations (F_ST), and (VIII) detecting excess or decrease in admixture contribution from one population. Here, we review and compare these approaches using available human genome-wide datasets to provide independent verification (or not) of regions found by different methods and using different populations. The lessons learned from GWSSs will be applied to identify genome signatures of historic selective pressures on genes and gene regions in other species with emerging genome sequences. This would offer considerable potential for genome annotation in functional, developmental and evolutionary contexts.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.

Association of orthodenticle with natural variation for early embryonic patterning in Drosophila melanogaster

Although it is well established that cis-acting regulatory variation contributes to morphological evolution between species, few concrete examples of polymorphism affecting developmental patterning within species have been demonstrated. Early embryogenesis in Drosophila is initiated by a gradient of Bicoid morphogen activity that results in differential expression of multiple target genes. In a screen for genetic variation affecting this process, we surveyed 96 wild-type lines of Drosophila melanogaster for polymorphisms in binding sites within 16 Bicoid cis-regulatory response elements. One common polymorphism in the orthodenticle (otd) early head enhancer is associated with a complex series of indels/substitutions that define two distinct haplotypes. The middle region of this enhancer exhibits an unusual pattern of nucleotide diversity that does not easily fit into standard models of selection and demography. Population Gene Expression Maps, generated by extracting binary expression profiles from normalized embryo images, revealed a ventral reduction of otd transcript abundance in one of the haplotypes that was recapitulated in expression of transgenic constructs containing the two alleles. We thus demonstrate that even a process as robust as early developmental patterning is affected by standing genetic variation, intriguingly involving otd, whose morphogenetic function bicoid is thought to have displaced during dipteran evolution.

IL4 gene polymorphism and previous malaria experiences manipulate anti-Plasmodium falciparum antibody isotype profiles in complicated and uncomplicated malaria

Background

The IL4-590 gene polymorphism has been shown to be associated with elevated levels of anti-Plasmodium falciparum IgG antibodies and parasite intensity in the malaria protected Fulani of West Africa. This study aimed to investigate the possible impact of IL4-590C/T polymorphism on anti-P. falciparum IgG subclasses and IgE antibodies levels and the alteration of malaria severity in complicated and uncomplicated malaria patients with or without previous malaria experiences.

Methods

Anti-P.falciparum IgG subclasses and IgE antibodies in plasma of complicated and uncomplicated malaria patients with or without previous malaria experiences were analysed using ELISA. IL4-590 polymorphisms were genotyped using RFLP-PCR. Statistical analyses of the IgG subclass levels were done by Oneway ANOVA. Genotype differences were tested by Chi-squared test.

Results

The IL4-590T allele was significantly associated with anti-P. falciparum IgG3 antibody levels in patients with complicated (P = 0.031), but not with uncomplicated malaria (P = 0.622). Complicated malaria patients with previous malaria experiences carrying IL4-590TT genotype had significantly lower levels of anti-P. falciparum IgG3 (P = 0.0156), while uncomplicated malaria patients with previous malaria experiences carrying the same genotype had significantly higher levels (P = 0.0206) compared to their IL4-590 counterparts. The different anti-P. falciparum IgG1 and IgG3 levels among IL4 genotypes were observed. Complicated malaria patients with previous malaria experiences tended to have lower IgG3 levels in individuals carrying TT when compared to CT genotypes (P = 0.075). In contrast, complicated malaria patients without previous malaria experiences carrying CC genotype had significantly higher anti-P. falciparum IgG1 than those carrying either CT or TT genotypes (P = 0.004, P = 0.002, respectively).

Conclusion

The results suggest that IL4-590C or T alleles participated differently in the regulation of anti-malarial antibody isotype profiles in primary and secondary malaria infection and, therefore, could play an important role in alteration of malaria severity.

Constructing genomic maps of positive selection in humans: where do we go from here?

Identifying targets of positive selection in humans has, until recently, been frustratingly slow, relying on the analysis of individual candidate genes. Genomics, however, has provided the necessary resources to systematically interrogate the entire genome for signatures of natural selection. To date, 21 genome-wide scans for recent or ongoing positive selection have been performed in humans. A key challenge is to begin synthesizing these newly constructed maps of positive selection into a coherent narrative of human evolutionary history and derive a deeper mechanistic understanding of how natural populations evolve. Here, I chronicle the recent history of the burgeoning field of human population genomics, critically assess genome-wide scans for positive selection in humans, identify important gaps in knowledge, and discuss both short- and long-term strategies for traversing the path from the low-resolution, incomplete, and error-prone maps of selection today to the ultimate goal of a detailed molecular, mechanistic, phenotypic, and population genetics characterization of adaptive alleles.

Evolutionary trajectories of primate genes involved in HIV pathogenesis

The current availability of five complete genomes of different primate species allows the analysis of genetic divergence over the last 40 million years of evolution. We hypothesized that the interspecies differences observed in susceptibility to HIV-1 would be influenced by the long-range selective pressures on host genes associated with HIV-1 pathogenesis. We established a list of human genes (n = 140) proposed to be involved in HIV-1 biology and pathogenesis and a control set of 100 random genes. We retrieved the orthologous genes from the genome of humans and of four nonhuman primates (Pan troglodytes, Pongo pygmaeus abeli, Macaca mulatta, and Callithrix jacchus) and analyzed the nucleotide substitution patterns of this data set using codon-based maximum likelihood procedures. In addition, we evaluated whether the candidate genes have been targets of recent positive selection in humans by analyzing HapMap Phase 2 single-nucleotide polymorphisms genotyped in a region centered on each candidate gene. A total of 1,064 sequences were used for the analyses. Similar median K(A)/K(S) values were estimated for the set of genes involved in HIV-1 pathogenesis and for control genes, 0.19 and 0.15, respectively. However, genes of the innate immunity had median values of 0.37 (P value = 0.0001, compared with control genes), and genes of intrinsic cellular defense had K(A)/K(S) values around or greater than 1.0 (P value = 0.0002). Detailed assessment allowed the identification of residues under positive selection in 13 proteins: AKT1, APOBEC3G, APOBEC3H, CD4, DEFB1, GML, IL4, IL8RA, L-SIGN/CLEC4M, PTPRC/CD45, Tetherin/BST2, TLR7, and TRIM5alpha. A number of those residues are relevant for HIV-1 biology. The set of 140 genes involved in HIV-1 pathogenesis did not show a significant enrichment in signals of recent positive selection in humans (intraspecies selection). However, we identified within or near these genes 24 polymorphisms showing strong signatures of recent positive selection. Interestingly, the DEFB1 gene presented signatures of both interspecies positive selection in primates and intraspecies recent positive selection in humans. The systematic assessment of long-acting selective pressures on primate genomes is a useful tool to extend our understanding of genetic variation influencing contemporary susceptibility to HIV-1.

Host sequence motifs shared by HIV predict response to antiretroviral therapy

BACKGROUND

The HIV viral genome mutates at a high rate and poses a significant long term health risk even in the presence of combination antiretroviral therapy. Current methods for predicting a patient's response to therapy rely on site-directed mutagenesis experiments and in vitro resistance assays. In this bioinformatics study we treat response to antiretroviral therapy as a two-body problem: response to therapy is considered to be a function of both the host and pathogen proteomes. We set out to identify potential responders based on the presence or absence of host protein and DNA motifs on the HIV proteome.

RESULTS

An alignment of thousands of HIV-1 sequences attested to extensive variation in nucleotide sequence but also showed conservation of eukaryotic short linear motifs on the protein coding regions. The reduction in viral load of patients in the Stanford HIV Drug Resistance Database exhibited a bimodal distribution after 24 weeks of antiretroviral therapy, with 2,000 copies/ml cutoff. Similarly, patients allocated into responder/non-responder categories based on consistent viral load reduction during a 24 week period showed clear separation. In both cases of phenotype identification, a set of features composed of short linear motifs in the reverse transcriptase region of HIV sequence accurately predicted a patient's response to therapy. Motifs that overlap resistance sites were highly predictive of responder identification in single drug regimens but these features lost importance in defining responders in multi-drug therapies.

CONCLUSION

HIV sequence mutates in a way that preferentially preserves peptide sequence motifs that are also found in the human proteome. The presence and absence of such motifs at specific regions of the HIV sequence is highly predictive of response to therapy. Some of these predictive motifs overlap with known HIV-1 resistance sites. These motifs are well established in bioinformatics databases and hence do not require identification via in vitro mutation experiments.

Promoter polymorphisms of IL2, IL4, and risk of gastric cancer in a high-risk Chinese population

Interleukin 2 (IL2) is a typical Th1 cytokine, and interleukin 4 (IL4) is an inducible Th2 cytokine. These cytokines are critical mediators of the Th1/Th2 balance and apoptosis potential and involved in the process of inflammation-mediated carcinogenesis in human organs, including the gastrointestinal tract. Therefore, we tested the hypothesis that functional variants in IL2 and IL4 were associated with risk of gastric cancer by genotyping two promoter polymorphisms in IL2 G-330T (rs2069762) and IL4 T-168C (rs2070874) in a case-control study of 1045 patients with incident gastric cancer and 1100 cancer-free controls in a high-risk Han Chinese population. We found that, compared with the IL4 -168TT genotype, heterozygous -168TC and combined -168TC/CC genotypes were associated with a significantly decreased gastric cancer risk [adjusted odds ratio (OR) = 0.81, 95% confidence interval (CI) = 0.67-0.98 for -168TC; OR = 0.83, 95% CI = 0.69-1.00 for -168TC/CC, respectively]. Furthermore, this significant protective effect was more evident for gastric cardia cancer patients (adjusted OR = 0.73, 95% CI = 0.56-0.95 for -168TC/CC vs. -168TT). For IL2 G-330T, subjects carrying GT/TT genotypes also had a significantly reduced risk of gastric cardia cancer (adjusted OR = 0.68, 95% CI = 0.46-0.99), compared with those carrying the GG genotype. Our results indicate that IL4 T-168C and IL2 G-330T promoter polymorphisms may contribute to the etiology of gastric cardia cancer in Chinese populations.

Introgression from modern hybrid varieties into landrace populations of maize (Zea mays ssp. mays L.) in central Italy

Landraces are domesticated local plant varieties that did not experience a deliberate and intensive selection during a formal breeding programme. In Europe, maize landraces are still cultivated, particularly in marginal areas where traditional farming is often practiced. Here, we have studied the evolution of flint maize landraces from central Italy over 50 years of on-farm cultivation, when dent hybrid varieties were introduced and their use was widespread. We have compared an 'old' collection, obtained during the 1950s, before the introduction of hybrids, and a recent collection of maize landraces. For comparison, a sample of maize landraces from north Italy, and of improved germplasm, including hybrids and inbred lines were also used. A total of 296 genotypes were analysed using 21 microsatellites. Our results show that the maize landraces collected in the last 5-10 years have evolved directly from the flint landrace gene pool cultivated in central Italy before the introduction of modern hybrids. The population structure, diversity and linkage disequilibrium analyses indicate a significant amount of introgression from hybrid varieties into the recent landrace populations. No evidence of genetic erosion of the maize landraces was seen, suggesting that in situ conservation of landraces is an efficient strategy for preserving genetic diversity. Finally, the level of introgression detected was very variable among recent landraces, with most of them showing a low level of introgression; this suggests that coexistence between different types of agriculture is possible, with the adoption of correct practices that are aimed at avoiding introgression from undesired genetic sources.

The -589C>T polymorphism in the interleukin-4 gene (IL-4) is associated with a reduced risk of myocardial infarction in young individuals

BACKGROUND

Inflammatory reactions contribute to the development of arterial disease. We investigated the role of interleukin-4 (IL-4) in the development of myocardial infarction (MI) by genotyping patients with MI and control subjects for the -589C>T (rs2243250) single nucleotide polymorphism (SNP), which tags a functional haplotype of IL-4.

METHODS AND RESULTS

Study of Myocardial Infarctions Leiden (SMILE) included 560 men with a first MI and 646 control subjects. The Valencia study included 305 patients with MI at T genotype was found [odds ratio (OR) 0.84; 95% CI 0.37-1.95 for -589TT and 0.82; 95% CI 0.62-1.07 for -589CT compared with -589CC]. In patients younger than 50 years, carriership of one or two -589T alleles was associated with a reduced risk of MI (OR 0.57: 95% CI 0.34-0.95). This result was replicated in the Valencia study, where carriers of one or two -589T alleles had a reduced risk of MI (OR 0.67: 95% CI 0.47-0.95), with a strong protective effect of the -598T allele in homozygous -589T (OR 0.33: 95% CI 0.10-1.05). In the control subjects of the Valencia study, the -589T allele was associated with reduced levels of F1+2.

CONCLUSION

Our data indicate that the IL-4 haplotype tagged by the -589T allele reduces the risk of MI in young individuals.

Worldwide population differentiation at disease-associated SNPs

BACKGROUND

Recent genome-wide association (GWA) studies have provided compelling evidence of association between genetic variants and common complex diseases. These studies have made use of cases and controls almost exclusively from populations of European ancestry and little is known about the frequency of risk alleles in other populations. The present study addresses the transferability of disease associations across human populations by examining levels of population differentiation at disease-associated single nucleotide polymorphisms (SNPs).

METHODS

We genotyped ~1000 individuals from 53 populations worldwide at 25 SNPs which show robust association with 6 complex human diseases (Crohn's disease, type 1 diabetes, type 2 diabetes, rheumatoid arthritis, coronary artery disease and obesity). Allele frequency differences between populations for these SNPs were measured using Fst. The Fst values for the disease-associated SNPs were compared to Fst values from 2750 random SNPs typed in the same set of individuals.

RESULTS

On average, disease SNPs are not significantly more differentiated between populations than random SNPs in the genome. Risk allele frequencies, however, do show substantial variation across human populations and may contribute to differences in disease prevalence between populations. We demonstrate that, in some cases, risk allele frequency differences are unusually high compared to random SNPs and may be due to the action of local (i.e. geographically-restricted) positive natural selection. Moreover, some risk alleles were absent or fixed in a population, which implies that risk alleles identified in one population do not necessarily account for disease prevalence in all human populations.

CONCLUSION

Although differences in risk allele frequencies between human populations are not unusually large and are thus likely not due to positive local selection, there is substantial variation in risk allele frequencies between populations which may account for differences in disease prevalence between human populations.

Identification and analysis of genomic regions with large between-population differentiation in humans

The primary aim of genetic association and linkage studies is to identify genetic variants that contribute to phenotypic variation within human populations. Since the overwhelming majority of human genetic variation is found within populations, these methods are expected to be effective and can likely be extrapolated from one human population to another. However, they may lack power in detecting the genetic variants that contribute to phenotypes that differ greatly between human populations. Phenotypes that show large differences between populations are expected to be associated with genomic regions exhibiting large allele frequency differences between populations. Thus, from genome-wide polymorphism data genomic regions with large allele frequency differences between populations can be identified, and evaluated as candidates for large between-population phenotypic differences. Here we use allele frequency data from approximately 1.5 million SNPs from three human populations, and present an algorithm that identifies genomic regions containing SNPs with extreme Fst. We demonstrate that our candidate regions have reduced heterozygosity in Europeans and Chinese relative to African-Americans, and are likely enriched with genes that have experienced positive natural selection. We identify genes that are likely responsible for phenotypes known to differ dramatically between human populations and present several candidates worthy of future investigation. Our list of high Fst genomic regions is a first step in identifying the genetic variants that contribute to large phenotypic differences between populations, many of which have likely experienced positive natural selection. Our approach based on between population differences can compliment traditional within population linkage and association studies to uncover novel genotype-phenotype relationships.