eISSN: 1644-4124
ISSN: 1426-3912
Central European Journal of Immunology
Current issue Archive Manuscripts accepted About the journal Special Issues Editorial board Abstracting and indexing Subscription Contact Instructions for authors Publication charge Ethical standards and procedures
Editorial System
Submit your Manuscript
SCImago Journal & Country Rank
2/2020
vol. 45
 
Share:
Share:
Clinical immunology

Non-classical human leukocyte antigen class I in Tunisian children with autism

Najla Kharrat
1
,
Rania Abdelhedi
1
,
Imen Gtif
1
,
Imen Ayadi
1
,
Roberta Rizzo
1
,
Daria Bortolotti
1
,
Nouha Bouayed Abdelmoula
2
,
Farhat Ghribi
1, 3
,
Ahmed Rebai
1
,
Ines Zidi
4

  1. Laboratory of Molecular and Cellular Screening Process, Centre of Biotechnology of Sfax, Sfax University, Tunisia
  2. Research Unit: UR17ES36 Genomics of Signalopathies in the Service of Medicine, Sfax University, Medical University of Sfax, Tunisia
  3. Department of Child and Adolescent Psychiatry, Hédi Chaker University Hospital, Sfax, Tunisia
  4. Laboratory Microorganismes and Active Biomolecules, Sciences Faculty of Tunis, University of Tunis El Manar, Tunis, Tunisia
Cent Eur J Immunol 2020; 45 (2): 176-183
Online publish date: 2020/07/27
Article file
- Non-classical.pdf  [0.16 MB]
Get citation
 
PlumX metrics:
 

Introduction

Autism spectrum disorders (ASD) are neurodevelopmental conditions characterized by qualitative impairments of social interaction and repetitive behavior, with restricted interests diagnosed before the age of 3 [1]. The etiology of ASD is not fully understood. Although several studies have focused on the etiology of this disease, none of them have provided clear answers about the real cause of ASD [2, 3].

Several immune abnormalities have been detected in children with ASD [1, 4]. Indeed, neurobiological studies in ASD highlighted pathways involved in neural development, synapse plasticity, structural brain abnormalities, neurological processes, cognition, and behavior [3, 4]. In fact, HLA molecules, which are expressed on the surface of all nucleated cells, binds two types of receptors [5]. Firstly, the HLA molecules interact with T cell receptors (TCR) presented on the T lymphocytes, stimulating the activation of adaptive responses [6]. Secondly, HLA molecules link to the natural killer (NK) cells via the killer immunoglobulin-like receptor (KIR), leading to a reduced NK cell activity [5]. Classical HLA-C and non-classical HLA-G molecules, binding maternal KIRs, are highly expressed in the fetal trophoblast region. In particular, the HLA-G is responsible for the compliance of mothers to tolerate genetically different fetal tissues during pregnancy [5]. Recently in 2013, Patterson et al. have shown that maternal and prenatal infections increased the risk of developing autism in offspring. Moreover, a prenatal infection can affect brain development of fetus in uterus [7]. Furthermore, non-classical HLA class Ib genes are formed during embryonic growth and are expressed in both the ovocyte and blastocyst stages [8].

Genetic polymorphisms of the major histocompatibility complex (MHC) genes coding for HLA molecules have been reported to be associated with ASD development [3]. The classical HLA class Ia (HLA-A, HLA-B, HLA-Cw, and HLA-DR) polymorphisms were systematically typed in ASD children and their relatives; both stratification and haplotype-based analysis were performed in order to take into account their potential linkage disequilibrium (LD) and co-segregation [3]. However, there is a limited number of studies investigating the relationship between the non-classical HLA class Ib (HLA-F, HLA-G, and HLA-E) and the etiology of autism. Therefore, this present study aimed to report the allele frequencies of six polymorphisms located in non-classical Ib HLA genes (HLA-F*01:02, HLA-F*01:03, HLA-F*01:04, 14pb Ins/Del, +3142C>G, and HLA-E*01:01/01:03) and investigate their potential association with ASD in autistic probands from Tunisian families. Also, a measurement of soluble levels of HLA-G in serum was performed in order to examine the potential correlation with HLA-G genetic polymorphisms.

Material and methods

Study population

Fifteen families with one child having a diagnosis of ASD were enrolled in the study, between January and December 2004, for a total of 45 individuals. Diagnoses were performed at the Department of Child and Adolescent Psychiatry, Hédi Chaker University Hospital, Sfax, Tunisia, based on the international classification of diseases 10th edition (ICD-10) [9]. In ICD-10, ASD include Asperger’s disorder, childhood disintegrative disorder, and mental retardation.

All individuals were born in Tunisia and were of Tunisian descent. Thirteen males and two females (mean age, 4.5 ±2 years) were ASD children, and their score at the child autistic rating scale (CARS) was ≥ 30 [10]. ASD children underwent examinations of mental status and neurological and neuropsychological evaluations by neurologists and psychiatrists from the “ERAHMA” association and Hédi Chaker University Hospital, Sfax, Tunisia. None of the autistic subjects presented an identifiable cause of their disease. Thirty individuals were parents.

Informed parental consent approved by the committee of medical ethics (Hédi Chaker Hospital Ethics Committee, Sfax, Tunisia) was obtained from all participants and legal guardians prior to inclusion into the study. The work has been carried out in accordance with the Declaration of Helsinki for experiments in humans.

Genotyping

Genomic DNA were extracted from EDTA anticoagulated peripheral blood samples (5 ml), using the standard phenol–chloroform method [11], and frozen at –20°C until use.

For the determination of the four following polymorphisms: HLA-E*01:01/01:03 A>G and HLA-F *(01:02 T>C, 01:03 C>T, and 01:04 A>C) variants alleles, the protocols of polymerase chain reaction with sequence specific primers (PCR-SSP) described by Lauterbach et al. and Pan et al. [12, 13], respectively, were used. PCR amplification was carried out in a 25µl volume, containing 10µl master mix (Qiagen), 1 pmol or 1.2 pmol of each specific primer, 1 pmol of each internal control primer, 100 ng genomic DNA, and 0.25µl of Taq DNA polymerase (Boiron GmbH., Germany). Nine and three sets of amplification primers were designed based on the published sequences (http://hla.alleles.org/data/txt/f_gen.txt), respectively, for HLA-F and HLA-E isoforms. Each PCR included specific primers pair and an internal control primers pair (that amplify a 439-bp fragment in the hormone growth factor gene) [14].

To identify the HLA-G +3142 C>G polymorphism, a PCR approach was used, followed by restriction enzyme analysis with 1µl of BaeGI (500 U/µl) (New England Biolabs, USA) with 10 µL of the amplification product in a final volume of 50 µl, as previously described by Graebin et al. [15]. The digestion was carried out overnight at 37°C and checked on 3% agarose gel.

The 14-bp Insertion/ Deletion (Ins/Del) polymorphism genotyping was performed by PCR for exon 8-specific primers according to the conditions previously described by Shankarkumar et al. [16]. Two PCR products were generated: 224-bp and 210-bp bands, respectively corresponding to ins and del. The PCR product was visualized and interpreted by two different observers. The sizing of PCR products was performed using 50 pb DNA ladder.

All PCR reactions were done by using a GenAmp PCR system 9700 thermocycler (Applied Biosystems, Fisher Scientific, USA).

Evaluation of soluble HLA-G serum levels

The sera were obtained from all participants and stored at -80°C for later analytical analysis. Soluble HLA-G (sHLA-G) levels (both secreted HLA-G5 and shed HLA-G1) were measured by a validated enzyme-linked immunosorbent assay (ELISA) method, using anti-HLA-G MEM-G/9 or 5A6G7 (Exbio, Prague, Czech Republic) as capture antibody, and anti-β2-microglobulin antibody conjugated with peroxidase (Dako) as secondary antibody [17]. The concentration of sHLA-G was estimated by comparing the absorbance at 450 nm after subtraction of the blank control absorbance values, with a standard curve. sHLA-G levels were normalized for flushing volume and total protein content (sHLA-G, ng/ml/total flushing sample, ml/total proteins, ng/ml). The standard curve was made by a serial dilution of the standard supernatants of HLA-G/721.221. The limit of sensitivity was 1 ng/ml. All samples were run in duplicate.

Allele frequencies

The allele frequencies of the studied HLA-E and HLA-F polymorphisms in other populations were extracted from the 1000 Genomes Project database, available from http://www.1000genomes.org/, while for the HLA-G polymorphisms, the allele frequencies of other ethnic populations were extracted from the study of Sabbagh et al. [18].

Statistical analysis

Data were analyzed using PLINK v1.07 software, a free open-source whole genome association analysis toolset (http://zzz.bwh.harvard.edu/plink/). The χ2 test was used to assess differences in allele frequencies between Tunisians and other populations, and to test the deviation from the Hardy-Weinberg equilibrium (HWE), which is a theory stating that the genetic variation in a population will remain constant from one generation to the next in the absence of disturbing factors. The genotypic, allelic, and haplotypic frequencies were calculated by PLINK v1.07 software. Continuous variables were analyzed using the Mann-Whitney U test and expressed as mean ± standard error of the mean (SEM). The transmission disequilibrium test (TDT) [19] was used to assess allelic association and performed using PLINK v1.07 software. TDT simultaneously measures linkage and association by comparing the allele transmitted and the one that is not transmitted from a heterozygous parent to the child.

Haploview software 4.2 was used to visualize the structure of pair-wise linkage disequilibrium (LD) between the polymorphism in terms of Lewontin’s D’ and the squared correlation coefficient between allele frequencies (r2) expressed as a function of D’ [20]. Haplotype blocks were defined using the “solid spine” incorporated in the haploview program. All SNPs with minor allele frequencies < 0.01% were excluded, and the minimum haplotype frequency was set at 1%. The confidence interval for strong LD was set between 0.7 and 0.98.

Results

Alleles and genotypes frequencies

HLA-G allele frequencies

The allelic frequency of the14bp-Del allele was 0.41. This frequency is close to the frequency of the minor allele in the African population (0.41) (Table 1).

For the polymorphism +3142C>G, the frequency of allele C was 0.48 in the Tunisian population. This frequency is near to the percentage found in European (Tuscan, Italy) (TSI), 0.45) and African population (Luhya from Webuye, Kenya (LWK), 0.47), but different from Asiatic population, except for the Han Chinese from South China (CHS) (Table 1).

Table 1

Allelic frequencies of the polymorphisms studied in Tunisian and other populations

ParameterHLA-FHLA-GHLA-E
Nucleotide sequence variationHLA-F*01:02HLA-F*01:03HLA-F*01:0414bp Ins/Del+3142C>GHLA-E*01:01/01:03
Nucleotide changeT/CC/TA/CIns/DelC/GA/G
Ancestral alleleCTCInsGG
Minor allele in Tunisian populationTCADelCG
Minor allele frequency (MAF) in Tunisian population (n = 45)T (0)C (0.25)A (0)Del (0.41)C (0.48)G (0.45)
Africa
GUI (n = 60)NANANA0.630.33NA
SER (n = 239)NANANA0.480.24NA
TOR (n = 30)NANANA0.350.43NA
YRI (n = 88)00.2600.360.350.34
YAN (n = 175)NANANA0.300.40NA
LWK (n = 97)0.0150.110.010.410.470.46
ASW (n = 61)0.0080.1200.450.390.53
Europe
POR (n = 60)NANANA0.520.38
IBS (n = 14)00.1100.430.500.5
TSI (n = 98)00.190.020.400.470.41
CEU (n = 85)00.210.030.350.560.39
GBR (n = 89)00.130.010.430.510.38
FIN (n = 93)00.330.020.280.660.48
Asia
JPT (n = 89)00.0100.250.250.70
CHB (n = 97)00.0500.390.340.58
CHS (n = 100)00.01500.290.450.64
America
SEB (n = 155)NANANA0.410.46NA
NEB (n = 128)NANANA0.380.45NA
CLM (n = 60)0.0080.1600.330.440.49
MXL (n = 66)0.0080.0800.430.410.45
PUR (n = 55)00.1900.310.500.55

[i] ASW – people of African ancestry from the Southwestern United States, CEU – Utah residents with Northern and Western European ancestry, CHB – Han Chinese from Beijing, CHS – Han Chinese from South China, CLM – Colombians from Medellin, Colombia, FIN – Finnish from Finland, GBR – British from England and Scotland, GUI – natives of Guinea-Bissau, IBS – Iberian populations from Spain, Ins – insertion, JPT – Japanese from Tokyo, Japan, LWK – Luhya from Webuye, Kenya, MXL – people of Mexican ancestry from Los Angeles, California, USA, n – number of individuals, NA – not available, NEB – Northeastern Brazilians from Recife, Pernambuco, Brazil, POR – Portuguese, PUR – Puerto Ricans from Puerto Rico, SEB – Southeastern Brazilians from Ribeirao Preto, Sao Paulo, Brazil, SER – Serer from Niakhar, Senegal, TOR – Tori from Tori-Bossito, Benin, TSI – from Tuscany, Italy, YAN – Yansi from Bandundu, Democratic Republic of the Congo, YRI – Yoruba from Ibadan, Nigeria. The ancestral allele state was inferred from dbSNP (Build.142). Minor allele frequency (MAF) in Tunisian population is provided by this current study.

[ii] The allele frequencies of the HLA-E and HLA-F polymorphisms were extracted from the 1000 Genomes Project database, available from (http://www.1000genomes.org/). The allele frequencies of the HLA-G polymorphisms were extracted from the study of Sabbagh et al. [18]

HLA-E allele frequencies

The frequency of allele HLA-E*01:01 was 0.45. This frequency is close to the percentage found in Latino-American population (people of Mexican ancestry from Los Angeles (0.45)), but higher than African population and lower than European (Table 1).

HLA-F allele frequencies

Only the HLA-F*01:03 was polymorphic in Tunisian population, with an allelic frequency equaling to 0.25 for the C allele. This frequency was different from those found in Asian and Latino American populations (Table 1). The HLA-F*01:02 and the HLA-F*01:04 provided one genotype [CC], with a similar frequency in comparison with African population, except for the Luhya from Webuye, Kenya, and different from European population for the HLA-F*01:04 allele (Table 1).

Association analysis with autism

All genotype distributions of the studied non-classical HLA-class Ib polymorphisms were in accordance with the HWE (p > 0.05). As shown in Table 2, no association was found between the six studied markers and autism

Table 2

Association analysis between polymorphisms and autism

CHRNucleotide sequence
variation
Positionχ2OR95% CIp-value
6HLA-F*01:0229723501NANANANA
6HLA-F*01:0429723805NANANANA
6HLA-F*01:03297252340.1111.250.3357-4.6550.7389
614bp Ins/Del298308040.06671.1430.4144-3.1520.7963
6+3142C>G298309720.052631.1110.4515-2.7340.8185
6HLA-E*01:01/01:03304902871.3161.7140.6749-4.3540.2513

[i] CHR – chromosome, OR – TDT odds ratio, 95% CI – confidence interval for TDT odds ratio, χ2 – TDT χ2 statistic, NA – not available

(p > 0.05). Twelve non-redundant haplotypes (representing a combination of HLA-F*01:02, HLA-F*01:03, HLA-F*01:04, 14pb Ins/Del, HLA-G +3142C>G, and HLA-E*01:01/01:03) inferred in the Tunisian population, as presented in Table 3. In fact, the two most frequent haplotypes were CGCICG and CGTDGA, with a frequency of 0.17 and 0.16, respectively. Only the haplotype CGTIGA was reported as associated with autism (p = 0.042).

Table 3

Haplotypes and association with autism

HaplotypeFrequencyTUCHISQp-value
CGCDCG0.03170.50.501
CGTDCG0.06171.030.0360.93160.334
CGCICG0.05090.7220.6020.437
CGTICG0.17586.863.601.0130.314
CGTDGG0.08061.880.850.38230.536
CGTICA0.049100NANA
CGCICA0.04011.021.250.023750.877
CGTICA0.12292.853.60.085630.769
CGTDGA0.0791.751.2440.087540.767
CGTDGA0.16363.84.360.03510.851
CGCIGA0.04812020.157
CGTIGA0.09620.535.534.1160.042

[i] T – transmitted minor allele count, U – un-transmitted, CHISQ – χ2 test, TDT – χ2 statistic

Linkage disequilibrium

We derived |D’| values to assess linkage disequilibrium (LD) between alleles of the 4 polymorphic variants (14bp-Ins/Del HLA-G, +3142C>G HLA-G, HLA-E*01:01/01:03, and HLA-F*01:03). The LD plot showed a correlation between 14bp-Ins/Del and +3142C>G polymorphisms of HLA-G gene. Based on pairwise analysis of these neighboring SNPs, we reported a weak LD (r2 = 0.22). However, no blocks of LD could be detected with any of the methods described in haploview 4.2 (Fig. 1).

Fig. 1

Linkage disequilibrium plot for HLA-G, HLA-E, and HLA-F

CEJI-45-41512-g001.tif

Soluble HLA-G dosage

The results showed high prevalence of HLA-G5 levels in offspring compared to the parents (positivity percentage = 27.22% vs. 17.39%, respectively), although this did not reach statistical significance (mean ±SEM: 0.68 ±0.46 vs. 0.80 ±0.6 ng/ml, respectively;p = 0.50). Concerning the levels of sHLA-G and sHLA-G1, no difference was reported between parents and offspring (p = 0.62 and p = 0.85, respectively).

Various differences in total sHLA-G levels among parents and children were observed when HLA-G genetics were considered. Interestingly, children with Ins allele had more enhanced sHLA-G levels and those with Del allele presented low sHLA-G without significance (Table 4). Indeed, parents with two ancestral alleles 14bp Ins (Ins/Ins) had a lower sHLA-G levels than those with the genotype Del/Del and Ins/Del. Interestingly, the offspring homozygous for 14bp Ins allele exhibited enhanced sHLA-G concentrations (Table 4).

Table 4

Total sHLA-G levels according to 14bp-Ins/Del and +3142C>G alleles and genotypes

ParametersHLA-G (ng/ml) in parentssHLA-G (ng/ml) in offspring
Mean ±SEM
14bp Ins/Del
Alleles
Ins0.72 ±0.251.90 ±0.64
Del1.81 ±0.910.72 ±0.55
p0.8700.168
Genotypes
Del/Del2.63 ±2.630 ±0
Del/Ins0.98 ±0.561.21 ±0.89
Ins/Ins0.55 ±0.362.59 ±1.44
p*0.450.26
+3142C>G
Alleles
C0.42 ±0.221.82 ±0.59
G1.78 ±0.711.43 ±0.66
p0.1570.878
Genotypes
CC0.6150 ±0.61501.405 ±1.405
CG0.2950 ±0.26401.205 ±0.8933
GG3.053 ±1.7951.657 ±1.657
p**0.0370.8

[i] bp – base pairs, Del – deletion, Ins – insertion, p – p-value of Mann-Whitney test, SEM – standard error of mean, * Del/Ins + Del/Del vs. Ins/Ins, ** CC + CG vs. GG

Furthermore, referring to +3142C>G polymorphism, parents with G allele showed the highest levels of sHLA-G. The difference between sHLA-G in +3142C>G genotype ([GG] vs. [CG] + [CC]) was statistically significant (p = 0.037). On the other hand, offspring had similar sHLA-G levels for all genotypes, without statistically distinct differences (([GG] vs. ([CG] + [CC])); p = 0.80).

The effect of haplotypes on sHLA-G levels in serum was also studied. The concomitant presence of Del and G alleles was enhanced with sHLA-G levels in parents (2.01 ±0 vs. 0.80 ±0.61 ng/ml). However, offspring with Ins/G haplotypes had the highest sHLA-G levels (3.31 ±1.66 vs. 1.27 ±0.57 ng/ml). Despite this fact, no significant p-value was reported (p = 0.566 and p = 0.197, respectively) (data not shown).

Discussion

The autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by impaired social interaction, restricted communications skills, and unusual repetitive behavior [4]. The underlying causes of autism are still unknown; nevertheless, both genetic and environmental causes are believed to contribute to the risk of autism [21]. Various clues suggest that a family history (like autoimmune disease), immunogenetics, epigenetic, and internal and external environmental factors are the potential triggers for this disease [22]. In fact, a possible autoimmune component was suggested as factor in this etiology [23]. Many studies reported an association between autism and autoimmune thyroiditis or hypothyroidism [24], rheumatic fever [25], rheumatoid arthritis, celiac disease ulcerative colitis, and psoriasis [26, 27].

Maternal infection and exposure to toxic compounds could alter the peripheral immune system, which can occur during fetal development and lead to long-term epigenetic changes in child [28]. The high-rate of autoimmune diseases in mothers of autistic children suggests that the autoimmune response, which takes place in the expectant mothers, adversely affects the developing fetus in utero [22]. In fact, the inflammatory process has been associated with autism, and a high level of inflammatory cytokines (IL-6, IFN-γ, TNF-α, and IL-8) was reported [22, 29]. A significant increase in the autoantibodies was found in the sera of autistic mothers [22, 30, 31]. Also, a raised number of B lymphocytes was observed against a decreased rate in T lymphocytes within autistic patients [22, 32]. Although the origins of abnormal immune activity in autism are still unknown, it was hypothesized that the mechanism might involve HLA genes and their products [23]. The importance of HLA in autoimmune conditions has been known for many years. In fact, many studies have shown some associations of classical MHC class I, II, and III alleles and autism [28, 33-36]. In particular, the widely admitted major susceptibility allele for rheumatoid arthritis (an endophenotype described in autistic patients) now is HLA-DRB1*04 loci.

Even though an evidence of association between the HLA system and autism was reported, no studies were carried out on the HLA class Ib. In fact, the non-classical Ib major histocompatibility complex includes three genes: HLA-G, HLA-E, and HLA-F. The generated molecules are characterized by tolerogenic functions and interact not only with NK cells, but also with T lymphocyte subsets and other immune cells, playing a major role in fetus acceptation [37, 38]. These HLA-Ib genes are characterized by low allelic polymorphisms [28, 37]. For the first time, we investigated the implication of the six popular polymorphisms of these genes within 15 Tunisian autistic families.

Concerning HLA-G gene, two polymorphisms +3241C>G and 14bp Ins/Del were examined, both located in the region of 3’UTR of the exon 8 of the gene. These polymorphic sites play a pivotal role in the regulation of HLA-G expression, on the one hand through modifying HLA-G mRNA stability [39], and on the other hand by the binding of specific microRNAs [40]. Controversial results have been reported among the association of HLA-G gene in different pathologies because of missing variables including ethnicity and genetic background information [18, 40]. In the present study, we reported the allelic frequency of +3241C>G and 14bp Ins/Del in Tunisian population. Furthermore, our frequencies are close to African population and the same observation was reported for the +3142C>G polymorphism. However, these frequencies differ from other populations, especially the Asiatic one [18, 40]. No strong linkage disequilibrium was observed between 14bp Ins/Del and +3142C>G. This result was in discordance to those found in Europe, but in accordance within African, American, and Asian populations [18, 40]. Moreover, the two most frequent haplotypes were CGCICG and CGTDGA, with a frequency of 0.17 and 0.16, respectively. These haplotypes differ at the polymorphic sites within HLA-F (only one SNP) HLA-G and HLA-E genes. This phenomenon, known as “Yin-Yang haplotypes” has been reported by Zhang et al. [41] for many other genes in the human genome and seems to be a general characteristic of eukaryotic genomes. An association with autism was reported for the haplotype CGTIGA. Despite this fact, neither an allelic nor a genotypic association was found in a response to the interaction between the studied SNPs.

HLA-G gene is characterized by seven possible isoforms that are generated by alternative splicing of the primary transcript of HLA-G. Four are membrane-bound (HLA-G1, -G2, -G3, and -G4), and three are soluble isoforms (HLA-G5, -G6, and -G7). A supplemental soluble form of HLA-G is the shedding HLA-G1 generated after metalloproteinase’ cleavage of the HLA-G1. We failed to identify any association between sHLA-G levels in parents with one or two 14bp-Del.

Additionally, we reported high sHLA-G in parents with GG genotype. This preliminary result is contrary to several studies highlighting that G allele is related to low sHLA-G expression because the +3142G allele is predicted to bind three miRNAs (miR-148a, -148b, and -152) [39, 42-44].

In autistics, we reported some haplotypes with increased sHLA-G, including firstly GG Ins/Ins, CC Ins/Ins, and CG Ins/Del. Based on the findings of Xu et al. who demonstrated that G/14bp-Ins haplotype is increased in human papillomavirus (HPV18) infected patient, we suggest a possible interference in sHLA-G levels in autistics [45]. The discrepancy between previous studies and ours’ may be due to ethnicity and sample size.

In the present study, for the HLA-F gene, three essential polymorphisms were studied: HLA-F*01:02 (C to T substitution at position 38 in exon 1), HLA-F*01:03 (C to A mutation at position 212 in exon 2), and HLA-F*01:04 alleles (T to C substitution at position 814 in exon 4) [46]. The HLA-F*01:02 and the HLA-F*01:04 alleles were non-polymorphic like in other populations in the world. The HLA-F*01:03 has the same allele frequency as European or African populations. Similar finding was also observed for the HLA-E*01:03.

Despite the absence of notable difference in alleles frequencies in Tunisian population compared to other populations in the world, our work still is the first to study HLA-G, -E, and -F polymorphisms and their eventual association with autism in Tunisian population. Our result reveals an association at the haplotypic level with CGTIGA haplotype. Nevertheless, no correlation between HLA-Ib (HLA-G, HLA-E, and HLA-F) and autism in Tunisian autistic families was found. Therefore, larger samples are still needed to clearly establish the association of HLA-Ib polymorphisms in autism. Our results suggest a potential association between haplotypes and sHLA-G levels in autistics.

Acknowledgments

We thank the “ERAHMA” association for their assistance in patients’ recruitment. This work was supported by the Ministry of Higher Education, Scientific Research, Information, and Communication Technologies of Tunisia.

Notes

[6] Conflicts of interest The authors declare no conflict of interest.

References

1 

Onore C, Careaga M, Ashwood P (2012): The role of immune dysfunction in the pathophysiology of autism. Brain Behav Immun 26: 383-392.

2 

Puangpetch A, Suwannarat P, Chamnanphol M et al. (2015): Significant Association of HLA-B Alleles and Genotypes in Thai Children with Autism Spectrum Disorders: A Case-Control Study. Dis Markers 1-7. doi:10.1155/2015/724935.

3 

Torres AR, Westover JB, Rosenspire AJ (2012): HLA Immune Function Genes in Autism. Autism Res Treat 1-13. doi:10.1155/2012/959073.

4 

Ebert DH, Greenberg ME (2013): Activity-dependent neuronal signalling and autism spectrum disorder. Nature 493: 327-337.

5 

Guerini FR, Bolognesi E, Chiappedi M, et al. (2015): An HLA-G*14bp insertion/deletion polymorphism associates with the development of autistic spectrum disorders. Brain Behav Immun 44: 207-212.

6 

Shiina T, Hosomichi K, Inoko H, et al. (2009): The HLA genomic loci map: expression, interaction, diversity and disease. J Hum Genet 54: 15-39.

7 

Patterson PH. (2011): Maternal infection and immune involvement in autism. Trends Mol Med 17: 389-394.

8 

Kochan G, Escors D, Breckpot K, et al. (2013): Role of non-classical MHC class I molecules in cancer immunosuppression. Oncoimmunology 2(11):e26491. doi:10.4161/onci.26491

9 

Misès R, Garret-Gloanec N, Coinçon Y (2010): Classification de l’autisme et des psychoses précoces, plaidoyer pour des convergences. Inf Psychiatr 86(3):223. doi:10.3917/inpsy.8603.0223

10 

Schopler E, Reichler RJ, DeVellis RF, et al. (1980): Toward objective classification of childhood autism: Childhood Autism Rating Scale (CARS). J Autism Dev Disord 10: 91-103.

11 

Sambrook J, Russell DW (2001): Molecular Cloning: A Laboratory Manual. third ed. New York, USA: Cold Spring Harbor Laboratory Press.

12 

Lauterbach N, Voorter CE, Tilanus MG (2012): Molecular typing of HLA-E. Methods Mol Biol 882: 143-158.

13 

Pan FH, Liu XX, Tian W (2013): Characterization of HLA-F polymorphism in four distinct populations in Mainland China. Int J Immunogenet 40: 369-376.

14 

Boulanger LM (2004): MHC class I in activity-dependent structural and functional plasticity. Neuron Glia Biol 1: 283-289.

15 

Graebin P, Veit TD, Alho CS, et al. (2012): Polymorphic variants in exon 8 at the 3’ UTR of the HLA-G gene are associated with septic shock in critically ill patients. Crit Care 16: R211. doi:10.1186/cc11845.

16 

Shankarkumar U, Shankarkumar A, Chedda Z, et al. (2011): Role of 14-bp deletion/insertion polymorphism in exon 8 of the HLA-G gene in recurrent spontaneous abortion patients. J Hum Reprod Sci 4: 143-146.

17 

Park GM, Lee S, Park B, et al. (2004): Soluble HLA-G generated by proteolytic shedding inhibits NK-mediated cell lysis. Biochem Biophys Res Commun 313: 606-611.

18 

Sabbagh A, Luisi P, Castelli EC, et al. (2014): Worldwide genetic variation at the 3′ untranslated region of the HLA-G gene: balancing selection influencing genetic diversity. Genes Immun 15: 95-106.

19 

Spielman RS, McGinnis RE, Ewens WJ (1993): Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet 52: 506-516.

20 

Zapata C (2000): The D’ mesure of overall gametic disequilibrium between pairs of multiallelic loci. Evolution 54: 1809-1812.

21 

Gadad BS, Hewitson L, Young KA, et al. (2013): Neuropathology and animal models of autism: genetic and environmental factors. Autism Res Treat 731935. doi:10.1155/2013/731935.

22 

Gesundheit B, Rosenzweig JP, Naor D, et al. (2013): Immunological and autoimmune considerations of Autism Spectrum Disorders. J Autoimmun 44 :1-7.

23 

Guerini FR, Manca S, Sotgiu S, et al. (2006): A Family Based Linkage Analysis of HLA and 5-HTTLPR Gene Polymorphisms in Sardinian Children with Autism Spectrum Disorder. Hum Immunol 67: 108-117.

24 

Molloy CA, Morrow AL, Meinzen-Derr J, et al. (2006): Familial Autoimmune Thyroid Disease as a Risk Factor for Regression in Children with Autism Spectrum Disorder: A CPEA Study. J Autism Dev Disord 36: 317-324.

25 

Keil A, Daniels JL, Forssen U, et al. (2010): Parental autoimmune diseases associated with autism spectrum disorders in offspring. Epidemiology 21: 805-808.

26 

Lenroot RK, Yeung PK (2013): Heterogeneity within Autism Spectrum Disorders: What have We Learned from Neuroimaging Studies? Front Hum Neurosci 7:733.

27 

Atladóttir HO, Pedersen MG, Thorsen P, et al. (2009): Association of family history of autoimmune diseases and autism spectrum disorders. Pediatrics 124: 687-694.

28 

Mostafa GA, Shehab AA, Al-Ayadhi LY (2013): The link between some alleles on human leukocyte antigen system and autism in children. J Neuroimmunol 255(1-2):70-74.

29 

Croonenberghs J, Wauters A, Devreese K et al. (2002): Increased serum albumin, γ globulin, immunoglobulin IgG, and IgG2 and IgG4 in autism. Psychol Med 32: 1457-1463.

30 

Ashwood P, Wills S, Van de Water J (2006): The immune response in autism: a new frontier for autism research. J Leukoc Biol 80: 1-15.

31 

Singer HS, Morris CM, Gause CD, et al. (2008): Antibodies against fetal brain in sera of mothers with autistic children. J Neuroimmunol 194: 165-172.

32 

Ratajczak H V (2011): Theoretical aspects of autism: biomarkers—a review. J Immunotoxicol 8: 80-94.

33 

Daniels WW, Warren RP, Odell JD, et al. (1995): Increased frequency of the extended or ancestral haplotype B44-SC30-DR4 in autism. Neuropsychobiology 32: 120-123.

34 

Torres AR, Maciulis A, Stubbs EG, et al. (2002): The transmission disequilibrium test suggests that HLA-DR4 and DR13 are linked to autism spectrum disorder. Hum Immunol 63: 311-316.

35 

Donadi EA, Castelli EC, Arnaiz-Villena A, et al. (2011): Implications of the polymorphism of HLA-G on its function, regulation, evolution and disease association. Cell Mol Life Sci 68: 369-395.

36 

Yucesoy B, Talzhanov Y, Johnson VJ, et al. (2013): Genetic variants within the MHC region are associated with immune responsiveness to childhood vaccinations. Vaccine 31: 5381-5391.

37 

Morandi F, Pistoia V (2014): Interactions between HLA-G and HLA-E in Physiological and Pathological Conditions. Front Immunol 5: 394.

38 

Moscoso J, Serrano-Vela JI, Pacheco R et al. (2006): HLA-G,-E and-F: Allelism, function and evolution. Transpl Immunol 17: 61-64.

39 

Tan Z, Randall G, Fan J, et al. (2007): Allele-Specific Targeting of microRNAs to HLA-G and Risk of Asthma. Am J Hum Genet 81: 829-834.

40 

Castelli EC, Mendes-Junior CT, Veiga-Castelli LC, et al. (2011): A Comprehensive Study of Polymorphic Sites along the HLA-G Gene: Implication for Gene Regulation and Evolution. Mol Biol Evol 28: 3069-3086.

41 

Zhang J, Rowe WL, Clark AG et al. (2003): Genomewide distribution of high-frequency, completely mismatching SNP haplotype pairs observed to be common across human populations. Am J Hum Genet 73: 1073-81.

42 

Rizzo R, Bortolotti D, Fredj N Ben, et al. (2012): Role of HLA-G 14bp deletion/insertion and +3142C&gt;G polymorphisms in the production of sHLA-G molecules in relapsing-remitting multiple sclerosis. Hum Immunol 73: 1140-1146.

43 

Picard C, Di Cristofaro J, Azzouz DF et al. (2013): Analyzing HLA-G polymorphisms in children from women with scleroderma. Hum Immunol 74: 468-472.

44 

Veit TD, Chies JAB (2009): Tolerance versus immune response — MicroRNAs as important elements in the regulation of the HLA-G gene expression. Transpl Immunol 20: 229-231.

45 

Xu H-H, Shi W-W, Lin A et al. (2014): HLA-G 3’ untranslated region polymorphisms influence the susceptibility for human papillomavirus infection. Tissue Antigens 84: 216-222.

46 

Pyo C-W, Williams LM, Moore Y, et al. (2006): HLA-E, HLA-F, and HLA-G polymorphism: genomic sequence defines haplotype structure and variation spanning the nonclassical class I genes. Immunogenetics 58: 241-251.

Copyright: © 2020 Polish Society of Experimental and Clinical Immunology This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0) License (http://creativecommons.org/licenses/by-nc-sa/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material, provided the original work is properly cited and states its license.
 
Quick links
© 2024 Termedia Sp. z o.o.
Developed by Bentus.