Gene-gene interaction and functional impact of polymorphisms on innate immune genes in controlling Plasmodium falciparum blood infection level

Population genetic and biophysical evidences reveal that purifying selection shapes the genetic landscape of Plasmodium falciparum RH ligands in Chhattisgarh and West Bengal, India

Background

Reticulocyte binding protein-like homologs (RHs) are currently being evaluated as anti-erythrocytic stage vaccine targets against Plasmodium falciparum malaria. Present study explores the possible evolutionary drivers shaping the genetic organization of Pfrhs in Indian parasite population. It simultaneously evaluates a putative gain-of-function variant of PfRH5, a keystone member of PfRH family.

Methods

Receptor binding regions of Pfrh1, Pfrh2a/b, Pfrh4 and whole Pfrh5 were amplified using blood samples of P. falciparum malaria patients from Chhattisgarh and West Bengal and sequenced. Assembled sequences were analysed using MEGA7 and DnaSPv6. Binding affinities of recombinant PfRH5 proteins with basigin (BSG) were compared using in silico (CHARMM and AUTODOCK) and in vitro (Circular dichroism, fluorescence spectroscopy and isothermal titration calorimetry) methods.

Results

Pfrh1 (0.5), Pfrh2a/b (0.875), Pfrh4 (0.667) and Pfrh5 (0.778) sequence changes corresponded to low frequency (< 0.05) variants which resulted in an overall negative Tajima’s D. Since mismatch distribution of none of the Pfrh loci corroborated with the model of demographic expansion, a possible role of natural selection formulating Pfrh sequence diversity was investigated. Among the 5 members, Pfrh5 displayed very high dN/dS (5.7) ratio. Nevertheless, the model of selective sweep due to presence of any advantageous substitutions could not be invoked as polymorphic nonsynonymous sites (17/18) for Pfrh5 exceeded significantly over the divergent (62/86) ones (p = 0.0436). The majority of extant PfRH5 sequences (52/83) differed from the reference Pf3D7 allele by a single amino acid mismatch (C203Y). This non-conservative alteration was predicted to lower the total interaction energy of that PfRH5_variant with BSG, compared to PfRH5_3D7. Biophysical evidences validated the proposition that PfRH5_variant formed a more stable complex with BSG. Thermodynamic association constant for interaction of BSG with PfRH5_variant was also found to be higher (Ka_variant = 3.63E6 ± 2.02E6 M⁻¹ and Ka_3D7 = 1.31E6 ± 1.21E6 M⁻¹).

Conclusions

Together, the study indicates that the genetic architecture of Pfrhs is principally shaped by purifying selection. The most abundant and ubiquitous PfRH5 variant harbouring 203Y, exhibits a greater affinity for BSG compared to PfRH5_3D7 possessing 203C allele. The study underscores the importance of selecting the functional allele that best represents circulating strains in natural parasite populations as vaccine targets.

Adaptive genetic diversity and evidence of population genetic structure in the endangered Sierra Madre Sparrow (Xenospiza baileyi)

The magnitude and distribution of genetic diversity through space and time can provide useful information relating to evolutionary potential and conservation status in threatened species. In assessing genetic diversity in species that are of conservation concern, several studies have focused on the use of Toll-like receptors (TLRs). TLRs are innate immune genes related to pathogen resistance, and polymorphisms may reflect not only levels of functional diversity, but may also be used to assess genetic diversity within and among populations. Here, we combined four potentially adaptive markers (TLRs) with one mitochondrial (COI) marker to evaluate genetic variation in the endangered Sierra Madre Sparrow (Xenospiza baileyi). This species offers an ideal model to investigate population and evolutionary genetic processes that may be occurring in a habitat restricted endangered species with disjunct populations (Mexico City and Durango), the census sizes of which differ by an order of magnitude. TLRs diversity in the Sierra Madre Sparrow was relatively high, which was not expected given its two small, geographically isolated populations. Genetic diversity was different (but not significantly so) between the two populations, with less diversity seen in the smaller Durango population. Population genetic structure between populations was due to isolation and different selective forces acting on different TLRs; population structure was also evident in COI. Reduction of genetic diversity in COI was observed over 20 years in the Durango population, a result likely caused by habitat loss, a factor which may be the main cause of diversity decline generally. Our results provide information related to the ways in which adaptive variation can be altered by demographic changes due to human-mediated habitat alterations. Furthermore, our findings may help to guide conservation schemes for both populations and their restricted habitat.

Diversity analysis of MSP1 identifies conserved epitope organization in block 2 amidst high sequence variability in Indian Plasmodium falciparum isolates

Background

Despite its immunogenicity, the polymorphic nature of merozoite surface protein 1, an important vaccine candidate for Plasmodium falciparum malaria, remains a concern. This study analyses the impact of genetic variability and parasite population structure on epitope organization of different MSP1 segments.

Methods

Altogether 98 blood samples collected from P. falciparum infected mild and severe malaria patients of Chhattisgarh and West Bengal were used to sequence regions encoding block 2 and MSP1-19 of msp1. Sequences were analysed using MEGA7, DnaSPv5, Arlequin3.5 and BepiPred.

Results

All three major MSP1 block 2 allele families namely K1, MAD20 and RO33 were detected in the samples and they together resulted in 41 indel variants. Chhattisgarh samples displayed an average MOI of 2.07 ± 1.59 which was higher in mild malaria and in age group < 18 years. Ultra-structure of block 2 alleles revealed that mutation and repeat expansion were two major mechanisms responsible for allelic variability of K1 and MAD20. Regions flanking block 2 were highly variable in Chhattisgarh with average mismatch differences (k) ranging from 1.198 to 5.156 for three families. In contrast, region encompassing MSP1-19 exhibited limited heterogeneity (k_Chhattisgarh = 1.45, k_{West Bengal} = 1.363). Of the 50 possible B cell linear epitopes predicted from block 2 variants, 94.9% (131 of 138) of the parasites could be represented by three conserved antigens.

Conclusions

Present data indicates that natural selection and transmission intensity jointly play a role in controlling allelic diversity of MSP1 in Indian parasite isolates. Despite remarkable genetic variability, a limited number of predominant and conserved epitopes are present in Indian parasite isolates reinstating the importance of MSP1 as a promising malaria vaccine candidate.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2592-y) contains supplementary material, which is available to authorized users.

A functional promoter polymorphism in interleukin 12B gene is associated with an increased risk of ovarian endometriosis

OBJECTIVE

To explore the role of functional genetic variant of IL12B in the pathogenesis of ovarian endometriosis.

STUDY DESIGN

This study included 815 patients with ovarian endometriosis and 788 women in the control group. Two polymorphisms were genotyped by a method of polymerase chain reaction and ligase detection reaction. To assess the biological significance of polymorphisms, we detected the level of IL12p40 protein expression in patients' eutopic endometrial tissues using quantitative real-time polymerase chain reaction (qRT-PCR).

RESULT(S)

Compared with the CTCTAA/CTCTAA genotype, the GC/GC and GC/CTCTAA genotypes of rs17860508 could significantly increase the risk of ovarian endometriosis development (OR = 1.58, 95% CI = 1.19-2.09; OR = 1.37, 95% CI = 1.08-1.74, respectively). Moreover, the IL12B mRNA was expressed at significantly higher levels in the eutopic endometrial tissue of patients with the GC/GC genotype than in patients with the CTCTAA/CTCTAA genotype. However, the rs3212227 polymorphism may not be associated with a risk of ovarian endometriosis.

CONCLUSION(S)

The rs17860508 polymorphism in the IL12B promoter region may influence the risk of developing ovarian endometriosis by altering the endometrial expression of IL12B of in Northern Chinese women.

Cord Blood IL-12 Confers Protection to Clinical Malaria in Early Childhood Life

Using a well-designed longitudinal cohort, we aimed to identify cytokines that were protective against malaria and to explore how they were influenced by genetic and immunological factors. 349 Mozambican pregnant women and their newborn babies were recruited and followed up for malaria outcomes until 24 months of age. Six Th1 cytokines in cord blood were screened for correlation with malaria incidence, of which IL-12 was selected for further analyses. We genotyped IL-12 polymorphisms in children/mothers and evaluated the genotype-phenotype associations and genetic effects on IL-12 levels. Maternal IL-12 concentrations were also investigated in relation to Plasmodium infections and cord blood IL-12 levels. Our data showed that high background IL-12 levels were prospectively associated with a low incidence of clinical malaria, while IL-12 production after parasite stimulation had the opposite effect on malaria incidence. IL-12 genotypes (IL-12b rs2288831/rs17860508) and the haplotype CGTTAGAG distribution were related to malaria susceptibility and background IL-12 levels. Maternal genotypes also exhibited an evident impact on host genotype-phenotype associations. Finally, a positive correlation in background IL-12 levels between maternal and cord blood was identified. Thus, cord blood background IL-12 concentrations are important for protecting children from clinical malaria, likely mediated by both genotypes (children&mothers) and maternal immunity.

Genetic Polymorphisms in Cytokine Genes in Colombian Patients with Ocular Toxoplasmosis

Toxoplasmosis is caused by infection with the protozoan parasite Toxoplasma gondii, which has the capacity to infect all warm-blooded animals worldwide. Toxoplasmosis is a major cause of visual defects in the Colombian population; however, the association between genetic polymorphisms in cytokine genes and susceptibility to ocular toxoplasmosis has not been studied in this population. This work evaluates the associations between polymorphisms in genes coding for the cytokines tumor necrosis factor alpha (TNF-α) (rs1799964, rs1800629, rs1799724, rs1800630, and rs361525), interleukin 1β (IL-1β) (rs16944, rs1143634, and rs1143627), IL-1α (rs1800587), gamma interferon (IFN-γ) (rs2430561), and IL-10 (rs1800896 and rs1800871) and the presence of ocular toxoplasmosis (OT) in a sample of a Colombian population (61 patients with OT and 116 healthy controls). Genotyping was performed with the "dideoxynucleotide (ddNTP) primer extension" technique. Functional-effect predictions of single nucleotide polymorphisms (SNPs) were done by using FuncPred. A polymorphism in the IL-10 gene promoter (-1082G/A) was significantly more prevalent in OT patients than in controls (P = 1.93e-08; odds ratio [OR] = 5.27e+03; 95% confidence interval [CI] = 3.18 to 8.739; Bonferroni correction [BONF] = 3.48e-07). In contrast, haplotype "AG" of the IL-10 gene promoter polymorphisms (rs1800896 and rs1800871) was present at a lower frequency in OT patients (P = 7e-04; OR = 0.10; 95% CI = 0.03 to 0.35). The +874A/T polymorphism of IFN-γ was associated with OT (P = 3.37e-05; OR = 4.2; 95% CI = 2.478 to 7.12; BONF = 6.07e-04). Haplotype "GAG" of the IL-1β gene promoter polymorphisms (rs1143634, rs1143627, and rs16944) appeared to be significantly associated with OT (P = 0.0494). The IL-10, IFN-γ, and IL-1β polymorphisms influence the development of OT in the Colombian population.

Variants in two gene members of the TNF ligand superfamily and hepatitis C virus chronic disease

AIM

To assess the possible correlation between single nucleotide polymorphisms (SNPs) of two members of TNF ligand superfamily genes, tumor necrosis factor-α (TNF-α) and lymphotoxin-α (LTA), and HCV chronic disease.

BACKGROUND

The causes of disease progression from hepatitis C virus (HCV) infection to chronic liver disease still remains unclear. Abnormal production of the cytokines alleged to be contributed to progression of the disease or viral persistence. Regulatory mechanisms that control the production of cytokines including genetic polymorphisms, especially at coding/regulatory regions of genes, may affect expression and secretion of the cytokines.

METHODS

In this case-control investigation, 258 individuals with serologically proven chronic HCV infection and 277 healthy controls were studied. Genotyping of rs1799964 variant of TNF-α and rs909253 intronic variant in LTA gene were performed. To confirm the results of genotyping, 10% of the specimens analyzed again by sequencing approach.

RESULTS

In this investigation, a significant association was observed between the TNF-α TC genotype and chronic HCV infection (P = 0.035). Moreover, the frequency of C allele was significantly different between control subjects in comparison with chronic HCV patients (P=0.02). On the other hand, no association was found between LTA gene polymorphism and susceptibility to chronic HCV infection.

CONCLUSION

These findings indicate that genetic variants like single nucleotide polymorphism in TNF-α rs1799964, could be a host factor associated with susceptibility to HCV chronic infection. However, further large scale investigations are needed to confirm this finding.

Toll-like receptor variation in the bottlenecked population of the Seychelles warbler: computer simulations see the 'ghost of selection past' and quantify the 'drift debt'

Balancing selection can maintain immunogenetic variation within host populations, but detecting its signal in a postbottlenecked population is challenging due to the potentially overriding effects of drift. Toll-like receptor genes (TLRs) play a fundamental role in vertebrate immune defence and are predicted to be under balancing selection. We previously characterized variation at TLR loci in the Seychelles warbler (Acrocephalus sechellensis), an endemic passerine that has undergone a historical bottleneck. Five of seven TLR loci were polymorphic, which is in sharp contrast to the low genomewide variation observed. However, standard population genetic statistical methods failed to detect a contemporary signature of selection at any TLR locus. We examined whether the observed TLR polymorphism could be explained by neutral evolution, simulating the population's demography in the software DIYABC. This showed that the posterior distributions of mutation rates had to be unrealistically high to explain the observed genetic variation. We then conducted simulations with an agent-based model using typical values for the mutation rate, which indicated that weak balancing selection has acted on the three TLR genes. The model was able to detect evidence of past selection elevating TLR polymorphism in the prebottleneck populations, but was unable to discern any effects of balancing selection in the contemporary population. Our results show drift is the overriding evolutionary force that has shaped TLR variation in the contemporary Seychelles warbler population, and the observed TLR polymorphisms might be merely the 'ghost of selection past'. Forecast models predict immunogenetic variation in this species will continue to be eroded in the absence of contemporary balancing selection. Such 'drift debt' occurs when a gene pool has not yet reached its new equilibrium level of polymorphism, and this loss could be an important threat to many recently bottlenecked populations.

A roadmap to multifactor dimensionality reduction methods

Complex diseases are defined to be determined by multiple genetic and environmental factors alone as well as in interactions. To analyze interactions in genetic data, many statistical methods have been suggested, with most of them relying on statistical regression models. Given the known limitations of classical methods, approaches from the machine-learning community have also become attractive. From this latter family, a fast-growing collection of methods emerged that are based on the Multifactor Dimensionality Reduction (MDR) approach. Since its first introduction, MDR has enjoyed great popularity in applications and has been extended and modified multiple times. Based on a literature search, we here provide a systematic and comprehensive overview of these suggested methods. The methods are described in detail, and the availability of implementations is listed. Most recent approaches offer to deal with large-scale data sets and rare variants, which is why we expect these methods to even gain in popularity.

TNF, IL12B, and IFNG Gene Polymorphisms in Serbian Patients with Psoriasis

Background

Psoriasis is a common chronic inflammatory skin disease with a strong genetic basis. Cytokines such as tumor necrosis factor alpha (TNF-α), interleukins (ILs) such are IL-12 and IL-23, and interferon gamma (IFN-γ) are released from various inflammatory and resident cells, and have been implicated in the initiation/maintenance of inflammation. Certain alleles of the aforementioned cytokines may be associated with disease susceptibility/severity.

Objective

To investigate the association of three common functional gene polymorphisms, namely TNF -308 G/A (rs1800629), IL12B (encoding the p40 subunit of IL-12/23) +1188 A/C (rs3212227), and IFNG +874 T/A (rs2430561) with psoriasis development and severity in Serbian patients.

Methods

We genotyped 130 patients with psoriasis (26 of whom also had psoriatic arthritis) and 259 controls; rs1800629 and rs3212227, and rs2430561, by real-time PCR assay.

Results

The TNF GG genotype was detected at a higher frequency in patients with psoriasis compared to control subjects (OR, 1.420; 95% CI, 0.870~2.403) without statistical significance (p=0.191). Lack of the TNF G allele was associated with lower psoriasis severity (p=0.007). The IL12B AC genotype was underrepresented in the patients with psoriatic arthritis compared to healthy subjects (OR, 0.308; 95% CI, 0.090~1.057; p=0.049). The distribution of the rs2430561 allele and genotype frequencies was similar between patients with psoriasis and controls.

Conclusion

Our study demonstrates an effect of the rs1800629 on psoriasis severity, and a marginal impact of the rs3212227 on susceptibility to psoriatic arthritis. Collectively, our results obtained in a Serbian cohort expand current knowledge regarding individual predisposition to psoriatic disease.

Two susceptibility loci to Takayasu arteritis reveal a synergistic role of the IL12B and HLA-B regions in a Japanese population

Takayasu arteritis (TAK) is an autoimmune systemic vasculitis of unknown etiology. Although previous studies have revealed that HLA-B*52:01 has an effect on TAK susceptibility, no other genetic determinants have been established so far. Here, we performed genome scanning of 167 TAK cases and 663 healthy controls via Illumina Infinium Human Exome BeadChip arrays, followed by a replication study consisting of 212 TAK cases and 1,322 controls. As a result, we found that the IL12B region on chromosome 5 (rs6871626, overall p = 1.7 × 10(-13), OR = 1.75, 95% CI 1.42-2.16) and the MLX region on chromosome 17 (rs665268, overall p = 5.2 × 10(-7), OR = 1.50, 95% CI 1.28-1.76) as well as the HLA-B region (rs9263739, a proxy of HLA-B*52:01, overall p = 2.8 × 10(-21), OR = 2.44, 95% CI 2.03-2.93) exhibited significant associations. A significant synergistic effect of rs6871626 and rs9263739 was found with a relative excess risk of 3.45, attributable proportion of 0.58, and synergy index of 3.24 (p ≤ 0.00028) in addition to a suggestive synergistic effect between rs665268 and rs926379 (p ≤ 0.027). We also found that rs6871626 showed a significant association with clinical manifestations of TAK, including increased risk and severity of aortic regurgitation, a representative severe complication of TAK. Detection of these susceptibility loci will provide new insights to the basic mechanisms of TAK pathogenesis. Our findings indicate that IL12B plays a fundamental role on the pathophysiology of TAK in combination with HLA-B(∗)52:01 and that common autoimmune mechanisms underlie the pathology of TAK and other autoimmune disorders such as psoriasis and inflammatory bowel diseases in which IL12B is involved as a genetic predisposing factor.

Structure of tumor necrosis factor-alpha haploblocks in European populations

DNA variants in the tumor necrosis factor-α (TNF) and linked lymphotoxin-α genes, and specific alleles of the highly polymorphic human leukocyte antigen B (HLA-B) gene have been implicated in a plethora of immune and infectious diseases. However, the tight linkage disequilibrium characterizing the central region of the human major histocompatibility complex (MHC) containing these gene loci has made difficult the unequivocal interpretation of genetic association data. To alleviate these difficulties and facilitate the design of more focused follow-up studies, we investigated the structure and distribution of HLA-B-specific MHC haplotypes reconstructed in a European population from unphased genotypes at a set of 25 single nucleotide polymorphism sites spanning a 66-kilobase long region across TNF. Consistent with the published data, we found limited genetic diversity across the so-called TNF block, with the emergence of seven common MHC haplotypes, termed TNF block super-haplotypes. We also found that the ancestral haplotype 8.1 shares a TNF block haplotype with HLA-B*4402. HLA-B*5701, a known protective allele in HIV-1 pathogenesis, occurred in a unique TNF block haplotype.