Paternal HLA-C and Maternal Killer-Cell Immunoglobulin-Like Receptor Genotypes in the Development of Autism

Killer-cell immunoglobulin-like receptors (KIRs) are a family of cell surface proteins found on natural killer cells, which are components of the innate immune system. KIRs recognize MHC class I proteins, mainly HLA-C and are further divided into two groups: short-tailed 2/3DS activating receptors and long-tailed 2/3DL inhibitory receptors. Based on the Barker Hypothesis, the origins of illness can be traced back to embryonic development in the uterus, and since KIR:HLA interaction figures prominently in the maternal–fetal interface, we investigated whether specific KIR:HLA combinations may be found in autism spectrum disorders (ASD) children compared with their healthy parents. This study enrolled 49 ASD children from different Israeli families, and their healthy parents. Among the parents, a higher frequency of HLA-C2 allotypes was found in the fathers, while its corresponding ligand 2DS1 was found in higher percentage in the maternal group. However, such skewing in KIR:HLA frequencies did not appear in the ASD children. Additionally, analysis of “overall activation” indicated higher activation in maternal than in paternal cohorts.

FigUre 1 | schematic diagram of activating (green) and inhibitory (red) receptors on natural killer cells, arranged according to centromic/telomeric location. These receptors are known to interact specifically with HLA ligands (HLA CI or C2, HLA A and B BW4). HLA-C group 1 (CI) allotypes serve as ligands for KIRs 2DS2, 2DL2, and 2DL3; HLA-C group 2 (C2) functions as ligands for KIRs 2DL1 and 2DS1. Altered natural killer (NK) cell activity has been reported in the periphery of ASD patients (7). However, these innate immune cells' function in the pathogenesis of autism is yet unclear. NK cells act specifically against infected or transformed cells, with their actions governed by a balance between activating and inhibitory receptors, each recognizing specific molecules on a given target (8). Killer-cell immunoglobulin-like receptors (KIRs) are a group of cell surface molecules involved in regulating activity of NK (and some T) cells. KIRs are classified as either inhibitory or activating. Their nomenclature is based on the number of extracellular Ig-like domains, and whether the cytoplasmic domain is long (L) or short (S). The KIR gene family is clustered within the LRC leukocyte receptor complex on chromosome 19q13.4, encoding a total of 14 genes and 2 pseudo genes. KIR genes exhibit allelic and haplotypic variability, while each mature NK cell expresses a specific set of KIRs.
Killer-cell immunoglobulin-like receptor counterparts, the HLA class I ligands, are also extremely polymorphic. This generates an additional level of functional diversity. Both KIRs and their HLA ligands are encoded on different chromosomes (19 and 6, respectively), and therefore segregate independently of each other. This expands the repertoire of possible profiles within any given population (9). Specific KIR-HLA combinations have been suggested as underlying susceptibility as well as resistance to various diseases and conditions (10); the functional implications are less clear. Genetic associations have been reported in infectious diseases (11), autoimmunity (12), inflammatory conditions, reproductive failure, and cancer. Biological models in such studies need to incorporate these receptor:ligand interactions, beyond mere consideration of KIR-HLA genotyping. Indeed, several methods have been proposed to analyze the overall genotyping detected in a specific subject or whole population (13,14).
The aim of our research was to genotype the KIR receptors known to interact with HLA ligands in autistic children and their non-autistic parents (mothers and fathers, separately), and to compare KIR receptor: HLA ligand frequencies between these groups. Our hypothesis was that a non-random distribution of these genes profiles should be found if interactions between these two classes of molecules are relevant to autism.
According to the Barker Hypothesis (15), the source of pathologies -including mental disorders such as ASD -can be traced back to embryonic existence in the intrauterine environment (16,17). KIRs are expressed by maternal NK cells (decidua, uterine lining), which are in close contact with HLA-expressing fetal tissues (placenta). Such prenatal interaction between maternal KIR and the fetal HLA in the uterus may lead to either enhanced NK activation or NK inhibition, which in turn may exert influence upon the placenta and/or developing fetus (18)(19)(20). Resulting immune activation and effects upon the central nervous system may lead to neurological and psychiatric disorders in predisposed patients. A mouse model for intrauterine maternal immune activation has already been demonstrated as a possible mechanism for autism-and schizophrenia-related behaviors in offspring (21).

Participants
This study comprised of 49 ASD-diagnosed children (38 males, 11 females, a male:female ratio of 3.4:1) from Jewish Israeli families, and their parents. Another six sets of parents of ASD-diagnosed children were enrolled, though no DNA samples were available from their children. DNA samples were collected in Israeli special education schools and treatment centers, by the Israeli Society for Autistic Children (ALUT). 1 Many of the probands had already been diagnosed with ASD by independent clinicians; two trained clinicians confirmed the diagnosis for autism or Pervasive Developmental Disorder-Not Otherwise Specified (PDD-NOS). All probands were diagnosed using the Autism Diagnostic Observation Scale-Generic (ADOS-G) and the Autism Diagnostic Interview-Revised (ADI-R). In addition, parents were interviewed using the Vineland Adaptive Behavior Scales-second edition (Vineland™-II) survey interview edition (VABS).
Probands ranged in age between 25 months to 33 years and 8 months. None of the subjects were known, according to parental reports, to have known genetic factors such as chromosomal aberrations, tuberous sclerosis, or other medical complications that could be related to autism. Parents were screened for psychiatric illness (including ASD); all results were found to be negative.
The Ethics Committee of the Israeli Health Ministry approved this study. All participants signed informed consent forms in advance.
genotyping DNA was extracted from either peripheral blood or buccal cells using a commercial DNA purification kit (Epicentre's Master Pure™, Madison, WI, USA).
The presence or absence of nine KIR genes known to interact with specific HLA ligands was determined by SSP-PCR (sequence-specific primer polymerase chain reaction), as previously described (22,23), using nine pairs of primers specific for these KIR genes. This mapping enables verification of presence or absence of a KIR, without discriminating specific heterozygous or homozygote status. A 440/750-bp fragment of the CRP/ CAMP, respectively, was included in each PCR as an internal control. The following KIR genes were amplified under identical conditions: 2DL2, 2DL1, 3DL1, 2DS2, 2DS4, and 3DS1. Their PCR thermo-profile comprised 2 min at 92°C; then 25 s at 91°C, 45 s at 65°C, and 45 s at 72°C for 4 cycles; 25 s at 91°C, 45 s at 60°C, and 45 s at 72°C for 26 cycles; 25 s at 91°C, 60 s at 55°C, and 45 s at 72°C for 5 cycles, and 10 min at 72°C (an extension step, in the case of 3DL1 and 3DS1 was set to 2 min). Other KIR amplifications (for 2DL3, 3DL2, and 2DS1) were initiated with denaturation for 2 min at 92°C; then 25 s at 91°C, 45 s at 57°C, and 45 s at 72°C for 4 cycles; 25 s at 91°C, 45 s at 54°C, and 45 s at 72°C for 26 cycles; 25 s at 91°C, 60 s at 50°C, and 45 s at 72°C for 5 cycles, and 10 min at 72°C. All primer sequences, expected PCR products and allelic recoveries are listed in Table 1.
The distinction of HLA-C KIR ligand groups C1 (N 80 ) and C2 (K 80 ) was also performed by SSP-PCR. The common HLA-C forward primer was either paired with the C1 (N 80 ) or C2 (K 80 ) reverse primers, leading to an amplification product of 139 bp ( Table 2). This mapping enables discrimination between heterozygous/homozygous status, i.e., C1-C1, C1-C2, or C2-C2. In each HLA PCR reaction, a 440-bp fragment of the CRP gene was included as an internal PCR control. The thermo-profile comprised an initial denaturation of 3 min at 95°C; then 10 s at 95°C, 30 s at 65°C, and 45 s at 72°C for 10 cycles; 10 s at 95°C, 30 s at 58°C, 45 s at 72°C for 22 cycles, and 10 min at 72°C.
HLA-A BW4 and HLA-B BW4-I 80 /T 80 were identified under the same conditions as HLA-C. Validation for HLA-A and B was also performed by selected samples previously genotyped by high    Table 2.
All amplified SSP-PCR products were visualized by agarose gel electrophoresis using a non-mutagenic fluorescent reagent ("Novel Juice"; GeneDireX, Las Vegas City, NV, USA). Selected sample from each amplified band was sent for sequencing.

statistical analysis
The Wilcoxon Signed-Ranks test (p < 0.05) was determined using SPSS software. When combination of more than one factor to specific receptor was detected, some different aspects were detected. The first refers to all ligands equally, and compares the sum of all ligands (relevant to specific receptor) found in each person to determine the mean value. This is defined as "Ligand Composite" ( Table 7). Secondly, an alternative designation of "Ligand One+" compares existence of at least one ligand to this specific receptor. This was the case in C1 + 2DL2/2DL3/2DS2, C2 + 2DL1/2DS1 ( Table 8). In contrast, for overall signaling ( Table 9) and in the case of A_BW4 + 3DL1/3DS1 and B_BW4 + 3DL1/3DS1 ( Table 7) only the first option is shown. All HLAB genotyping analysis is described as "total, " meaning having I 80 /T 80 .

resUlTs hla and Kirs Frequency comparison
Comparison of the nine different KIR genes and HLA frequencies revealed several main differences between parents of autistic children (Tables 3-5).
The most surprising finding is higher frequency of HLA-C2 allotypes in fathers of ASD children (mean 0.73 vs. 0.52, p = 0.038), in contrast to C1 frequencies, which were the same in both mothers and fathers (0.88 vs. 0.91) ( Table 5). On the other hand, 2DS1, which is known to interact specifically with HLA-C2, was found in higher percentage (0.55 vs. 0.34, p = 0.022) in the maternal group (Table 4). Additionally 3DS1, with its specificity against BW4 allotypes, was found to be higher in mothers than in their ASD children (0.51 vs. 0.36, p = 0.044) ( Table 4). Other KIR frequencies (2DL2, 2DL1, 2DL3, 3DL1, 3DL2, 2DS2, and 2DS4) were similar in both mothers and fathers, and corresponded with values expected in Asian populations. 2 Most KIR and HLA frequencies in ASD children were similar to those detected in parents.
Positive correlation between parent (mother or father) and child, reflecting the likelihood of transferring the gene (Tables 3-5) were found relative to most of genes, which may reflect a high probability of homozygous genotyping for these genes.

Kir signaling Through hla Molecules
Relevant combinations between KIR receptors and their specific HLA ligands were detected ( Table 6). ASD children did not exhibit  high frequency of the combination HLA-C2 + KIR2DS1, despite elevated paternal HLA-C2 and maternal 2DS1 frequencies.
Correlational analysis between parent (mother or father) and child reflected the likelihood of transferring a gene with its corresponding ligand. Significant values were reported in most cases ( Table 6). It is important to mention, these gene families are located on different autosomal chromosomes, supporting homozygous status for most of the genes.
The next evaluation involves the examination of particular HLA ligands with presence/absence of all relevant activating/ inhibitory KIR genes ( Table 7). "Ligand Composite" analysis (see Materials and Methods) was compared between ASD children and their parents, as well as "Ligand One+" values ( Table 8). Significant differences were found between parental groups in the case of C2 + 2DL1/2DS1 "Ligand One+" (p = 0.040).

Overall Balance of activation and inhibition
Final analysis (Table 9) presents a full view of the genotyping. It may allow insight to the overall expected interaction of KIR/ HLA in Israeli autistic individuals and their non-autistic parents,   based on all genotype mapping. Mothers have significantly higher levels of activating signals than their male counterparts (p = 0.049).

DiscUssiOn
Despite the increased attention garnered by autism research in recent years, the condition is currently diagnosed only after the infant reaches 2 or 3 years of age. Findings in autistic patients indicate multifactorial inheritance, while definitive genetic factors can be diagnosed in only a small percentage of cases. Immunological studies revealed a unique and distinguishing profile in ASD patients compared to controls. In the context of NK cells, altered gene expression, absolute number and function have been reported (7,24,25). Reports of modified NK function led us to investigate the highly polymorphic cell surface proteins on NK cells, called KIRs, to search for genetic background underlying the immune abnormalities. By interacting with MHC class I proteins expressed on all cell types, KIRs regulate the killing function of NK.
The most significant finding in this study is the unequal distribution of HLA-C2 (higher in fathers of ASD children) known to interact with KIR 2DS1 (higher in mothers of ASD children) to activate NK function. Analysis of "overall activation" ( Table 9) signaling pointed to higher activation in maternal than paternal cohorts. Most HLA and KIR frequencies in ASD children ranged between those of their parents, indicating regular independent genetic inheritance from parents to children. This reinforced the importance of the specific HLA-KIR combinations between parents, though the mechanisms involved are still unclear. Our limitation in interpreting these findings lies in the absence of unusual prevalence of these genes or their combination in the ASD children, in comparison to their healthy parents.
Some reports correlated these initial interactions between parental cells as a possible mechanism for subsequent disease development or other complications. preeclampsia) (26,27) has already been reported. Specifically, KIR 2DS1 seemed to be protective because of its significant lack in RM women. We find this plausible, since activating KIR 2DS1 compete with inhibitory KIR 2DL1 in interaction with HLA-C2, such that a predominance of KIR 2DS1 can overcome the strong inhibitory effect mediated by KIR 2DL1. All these findings support the importance of the innate immune system (mediated by maternal uterine NK and fetal trophoblast cells), modulated by KIR-HLA interaction to play a beneficial role in reproductive success. Some papers have already pointed out differential segregation of KIR-HLA genes in autistic patients. Our study expands upon this by focusing within ASD families. Torres et al. (28) examined a large autistic cohort, revealing a highly significant increase in the activating KIR gene 2DS1 and its cognate HLA-C2 ligand. This increase in activating KIR may explain, at least partially, the high frequency of autoimmune diseases in autism patients. Guerini et al. (29) compared ASD children to their non-ASD mothers and showed both that activating KIR/HLA combinations were increased and inhibitory KIR/HLA were reduced in ASD children. This difference was proposed as having resulted from the genetic influence of the mothers, as these molecules' expression and interaction are relevant in the surrounding uterus.
Ethnic populations are known to differ in KIR genotype frequencies and genotype content (30). To the best of our knowledge, this is the first study to research KIR and HLA combinations in the Israeli Jewish population. As shown in the Allele Frequency Net Database, 3 our geographical location has higher frequency of HLA-C2 (around 0.4-0.7 in Africa and the Middle East). However, we showed here a significant skewing between maternal and paternal C2 frequencies. It raises an interesting question, since HLA clusters map in autosomal chromosomes, leading to independent segregation without respect to gender.
Genotype mapping of this system has certain limitations. While theoretically all gene segments identified by our SSP-PCR system may be expressed on the cell surface (meaning the functional KIR repertoire may depend on the KIR genotype), the transcription and translation processes may lead to different phenotype profiles. Another limitation of the SSP-PCR method in this context is the inability to determine the heterozygous/ homozygous status of each of these KIRs. Our future plans are to determine not only genotype, but also RNA and protein levels. Another point to take into account is that NK cell function is a very complicated process determined upon the balance of an array of receptors, of which KIR is only one family. Moreover, the affinity of these interactions may differ (e.g., KIR 2DL3: HLA-C1 interaction is thought to be weak, compared to KIR 2DL2:HLA-C1 or KIR 2DL1:HLA-C2), such that one interaction may wield greater influence than another. It is also becoming increasingly clear that genetic factors, such as promoters and epigenetic mechanisms, are important in controlling NK cell receptor expression and function.
This study comprised a relatively small cohort of ASDaffected families -though statistically verified based on the 3 http://www.allelefrequencies.net/ minor allele frequency (MAF) of C2 (39.5%). In the absence of a control group, our findings were compared with the general population data of our geographical area. Additionally, the clinical patient data have not been thoroughly listed here, places limitations upon the reader in drawing parallels to other patient cohorts. The trends revealed in this preliminary investigation need to be further corroborated in ongoing and expanded research efforts.
Finally, we would like to emphasis another interesting aspect. On the one hand, different pregnancy complications (FGR, spontaneous abortion, placental abruption) are known to stem (at least to some degree) from immunological factors. On the other hand, perinatal factors unrelated to the immunological milieu may also be important in the pathogenesis of ASD. Certain complications during gestation were found to occur with increased frequency in mothers of ASD patients; an immunological basis has been suggested as underlying a portion of these. In some cases, intermediate factors may play a partial role in the association (e.g., reduced birth weight).

cOnclUsiOn
The aim of our research was to genotype the KIR receptors and their relevant HLA ligands in autistic children and their non-autistic parents, in order to screen for possible unique combinations. The most interesting finding in this study is higher frequency of HLA-C2 allotypes in the paternal group besides higher percentage of KIR 2DS1 in the maternal group. Such interaction may lead to NK activation. An absence of unusual prevalence of these genes or their combination in the ASD children, limited our ability to interpret these findings. Moreover, analysis of "overall activation" signaling pointed to higher activation in maternal than paternal cohorts.

aUThOr cOnTriBUTiOns
Mrs. MG is a Ph.D. student. She is the major contributor to this manuscript doing the thinking, bench work, and manuscript writing. Mrs. KA contributed to the bench work as well as writing the manuscript. Prof. NY is an initiator of the original study and has performed the clinical diagnosis of all the participants of the study. Prof. RE is the initiator of original study recruited the families and the original head of the lab. Dr. DM is the head of the lab the study was conducted in, supported the study, and contributed to the manuscript writing and editing.
acKnOWleDgMenTs This work was supported by a research grant (#232-12-13B) from the National Institute for Psychobiology in Israel, founded by the Charles E. Smith Family, and the Harris Foundation. The authors wish to thank Mrs. Shoshana Israel, Ph.D., and Mrs. Amal Halabi, M.Sc., for their help with validating HLA primers, and Dr. Lubov Nemanov for her assistance in processing patient samples. MG was supported by the Hoffman Leadership and Responsibility Fund, at the Hebrew University.