Common Genetic Variants Found in HLA and KIR Immune Genes in Autism Spectrum Disorder
- 1Center for Persons with Disabilities, Utah State University, Logan, UT, USA
- 2Biology Department, Utah State University, Brigham City, UT, USA
- 3BSP CCR Genetics Core, SAIC-Frederick, Inc., Frederick National Laboratory for Cancer Research, Frederick, MD, USA
- 4Center for Integrated BioSystems, Utah State University, Logan, UT, USA
- 5Division of Research, Kaiser Permanente of Northern California, Oakland, CA, USA
The “common variant—common disease” hypothesis was proposed to explain diseases with strong inheritance. This model suggests that a genetic disease is the result of the combination of several common genetic variants. Common genetic variants are described as a 5% frequency differential between diseased vs. matched control populations. This theory was recently supported by an epidemiology paper stating that about 50% of genetic risk for autism resides in common variants. However, rare variants, rather than common variants, have been found in numerous genome wide genetic studies and many have concluded that the “common variant—common disease” hypothesis is incorrect. One interpretation is that rare variants are major contributors to genetic diseases and autism involves the interaction of many rare variants, especially in the brain. It is obvious there is much yet to be learned about autism genetics. Evidence has been mounting over the years indicating immune involvement in autism, particularly the HLA genes on chromosome 6 and KIR genes on chromosome 19. These two large multigene complexes have important immune functions and have been shown to interact to eliminate unwanted virally infected and malignant cells. HLA proteins have important functions in antigen presentation in adaptive immunity and specific epitopes on HLA class I proteins act as cognate ligands for KIR receptors in innate immunity. Data suggests that HLA alleles and KIR activating genes/haplotypes are common variants in different autism populations. For example, class I allele (HLA-A2 and HLA-G 14 bp-indel) frequencies are significantly increased by more than 5% over control populations (Table 2). The HLA-DR4 Class II and shared epitope frequencies are significantly above the control populations (Table 2). Three activating KIR genes: 3DS1, 2DS1, and 2DS2 have increased frequencies of 15, 22, and 14% in autism populations, respectively. There is a 6% increase in total activating KIR genes in autism over control subjects. And, more importantly there is a 12% increase in activating KIR genes and their cognate HLA alleles over control populations (Torres et al., 2012a). These data suggest the interaction of HLA ligand/KIR receptor pairs encoded on two different chromosomes is more significant as a ligand/receptor complex than separately in autism.
Understanding genetic contribution in autism spectrum disorder (ASD) is a major challenge due to genetic complexities of a multifaceted disease. ASD is a puzzling neurodevelopmental condition characterized behaviorally by deficits in social-communication and the presences of restrictive stereotyped behaviors with involvement of genetic, environmental, and neurological components. Although, multiple etiologies are purported to exist, several studies have noted a strong familial clustering suggesting that heredity is of major importance (Bolton et al., 1994; Bailey et al., 1995; Constantino et al., 2013).
The study of human genetics has advanced rapidly over the last 30 years and researchers have adopted new and better methods in an attempt to unravel the extraordinary genomic complexity of autism. An early attempt to find DNA regions that associate with ASD involved the examination of several hundred microsatellites (also known as short tandem repeats) across the genome (Lamb et al., 2000). It became apparent early on that this was not answering questions about autism genetics and this approach was largely abandoned. A significant flaw in these early studies involved the use of siblings as control populations (Philippe et al., 1999). It is now thought that siblings of subjects with autism are unsuitable controls as they have many of the same risk polymorphisms as their autistic siblings. The next attempt involved the studies of purported autism candidate genes. This research typically involved the examination of single nucleotide polymorphisms (SNPs) between subjects and controls in a candidate gene thought to be involved in the etiology of autism. Although a few candidate genes were identified with apparent association to autism, overall the studies hypothesizing the involvement of purported candidate genes were often contradictory and proved to be unreliable (State, 2010).
Genome wide association studies (GWAS) that examine over half a million SNPs have been a little more successful in identifying rare genetic variants (rare variants) in a small number of ASD cases; however, common genetic variants (common variants) have not clearly been identified in GWAS studies (Goldstein, 2009; Manolio et al., 2009; Tao et al., 2016). One interpretation is that rare variants are major contributors to genetic diseases (McClellan and King, 2010) and that ASD involves the interaction of many rare variants, especially in the brain (Irimia et al., 2014; Krystal and State, 2014). The newest research tool involves whole genome sequencing (WGS) to determine the DNA sequence of an individual's entire genome. On the surface this sounds all encompassing; however, weaknesses of this approach besides cost include distorted amplification of certain DNA targets and the difficulty in aligning SNPs to the proper chromosome (Laver et al., 2016).
Several studies have noted a strong familial clustering suggesting that heredity of risk-associated genes is about 50–60% and that about 50% of this risk resides in common variants (Gaugler et al., 2014). The “common variant-common disease” hypothesis was proposed nearly a 100 years ago to explain common diseases with a strong genetic association (Fisher, 1918). This model suggests that each common risk variant (greater or less than 5% of the control population) confers a small degree of risk and that a genetic disease is the result of the combination of several common variants. Common variants (low penetrance) are genetically old and found in multiple populations, whereas rare variants (high penetrance) are new and found in specific populations (Schork et al., 2009). The examination of whole-genome sequence data of autism quartet families (two autism affected siblings and parents) suggested that the majority (69%) of autism affected siblings carried different rare variant mutations (Yuen et al., 2015), supporting the idea that these rare variants are not established in the studied population. Over the years, much attention has been given to discovering the genetic underpinnings of ASD; however, immune abnormalities have consistently been found in this disorder, and recently the number and impact of immune studies has reached a crescendo, including in the area of immunogenetics (Westover et al., 2011; Sweeten and McDougle, 2016; Young et al., 2016).
Basic Immune System
Immunology is the study of a system of molecules, tissue and cells that recognize and attack foreign invaders and abnormal cells that endanger the individual. The entire set of immune genes that contribute to immune function is now referred to as the “Immunome” consisting of about 900 genes across various chromosomes that encode a variety of different proteins to accomplish immune surveillance (Ortutay and Vihinen, 2009).
A major component of the immune system involves genes in the major histocompatibility complex (MHC) on chromosome 6 in humans. The MHC encodes human leukocyte antigen (HLA) proteins and the terms HLA and MHC are often used interchangeably. The human immune system has two main components, the innate immune system and the adaptive immune system. Both systems have humoral and cellular components (Figure 1), but innate immunity is generally considered to be non-specific whereas humoral immunity provides a specific response to pathogens and foreign antigens. Innate immunity is germ-line encoded, does not adapt after exposure to antigens, and is constitutively expressed leading to an immediate and rapid response to antigens. This is the first line of defense against pathogens, utilizing complement activation and other mechanisms to destroy most microorganisms, typically within hours of exposure.
Figure 1. The two branches of the immune system: Boxes display humoral and cellular components in innate and adaptive immunity.
The classical class I HLA proteins, expressed on essentially all nucleated cells of the body, serve as transducers by presenting antigens to T-cell and B-cell receptors, thereby activating the adaptive immune system for cellular cytotoxicity and antibody production, respectively. The adaptive HLA class II genes produce proteins expressed by specialized immune cells to present antigens to T-cells. HLA-class II genes can recombine by somatic mutation allowing for broad antigen adaptation. HLA-class I and II alleles are some of the most genetically variable coding loci in mammals.
Innate Immunity Components
The innate system has soluble complement proteins that help or “complement” the clearance of pathogens in various ways, including acting as an opsonin and by forming a membrane attack complex to kill cells. The complement system contains over 25 proteins and protein fragments that are mainly synthesized in the liver, but also exist in the brain. The innate cellular components include a variety of cell surface and cytoplasmic pattern recognition receptors (PRR) expressed by granulocytes, macrophages, and dendritic cells which bind to pathogen-associated molecular patterns (PAMPs) associated with components of microorganisms. An innate cellular system important in removing aberrant cells (virally infected and transformed) involves the interaction of HLA cognate receptors with killer-cell immunoglobulin-like receptors (KIRs) on NK-cells.
Adaptive Immunity Components
HLA proteins are important in presenting antigen to CD4+ (helper T-cells) and CD8+ (cytotoxic T-cells) and the extraordinary diversity in HLA alleles allows for exquisitely sensitive antigen recognition. Adaptive immune responses are slower than those elicited by the innate immune system; however, long term immunological memory is invoked. Memory involves the recognition of antigens from invading organism for extended periods of time, often decades. This is accomplished by the creation of tailor-made memory cells. Antibodies are the main soluble component of the adaptive system; however, cytokines are also important. Memory cells will continue to survive in the body for decades then re-exposure to the antigen activate them to produce the respective antibodies in unlimited numbers. This production of numerous antibodies is accomplished by gene recombination in B-cells followed by a complicated cellular selection process which results in antibodies that are very specific to the foreign antigens. The adaptive immune system is comprised of the cellular components T and B-cells. T-cells have cell-surface receptors that recognize peptides bound to HLA class I and class II proteins, a critical step in cell-mediated immunity.
The cellular components of the adaptive immune system involve T and B-cells. T-cells are important in cell-mediated immunity and are distinguished by the T-cell receptors which recognize short peptides bound to HLA class I and class II proteins. Several subsets of T-cells with distinct functions exist. For example, cytotoxic T-cells seek and destroy virally infected cells and tumor cells while helper T-cells activate macrophages and assist B-cells to mature into antibody producing plasma cells and memory B-cells.
Immunity and Autism
One of the earliest immune associations observed in ASD was by Stubbs et al. (1985) when they noted that pairs of ASD affected children share HLA haplotypes more often than non-affected pairs. The HLA complex is of major interest in medical research as genes/proteins in this region are involved in many immune processes such as organ transplantation, autoimmunity, resistance to specific pathogens, inflammation, ligands for immune cell receptors, and the complement cascade. There is now considerable evidence from association and linkage studies that suggests the involvement of HLA alleles/genes/haplotypes in the etiology of ASD (Torres et al., 2006, 2012a).
HLA Gene Complex
The HLA complex is fascinating as there are about 200 genes that encode proteins with various functions including ligands, receptors, cytokines, signaling factors, heat shock proteins, transcription regulators, and most importantly the recognition of self from non-self. The HLA gene map in Figure 2 (Shiina et al., 2009) illustrates this complexity. The entire 200 HLA gene region is inherited as a single 4.5 million base pair extended haplotype of which there are thousands; most very rare. A genetic haplotype is a contiguous combination of DNA sequences inherited on a single chromosome: eight of the more common extended haplotypes with strong disease associations have been entirely DNA sequenced (Horton et al., 2008). Various extended haplotypes containing millions of base pairs can have identical smaller haplotypes within their borders meaning that they must be carefully examined. Haplotype structure can be important in the understanding of allele-specific events such as protein structure, methylation, outcomes in transplantation, and disease prediction (McQueen et al., 2007). The HLA gene map does not fully illustrate the complexity of the class I and class II regions that have thousands of alleles important in determining self from non-self.
Figure 2. Gene map of the human leukocyte antigen (HLA) region. The major histocompatibility complex (MHC) gene map corresponds to the genomic coordinates of 29,677,984 (GABBR1) to 33,485,635 (KIFC1) in the human genome build 36.3 of the National Center for Biotechnology Information (NCBI) map viewer. The regions separated by arrows show the HLA subregions such as extended class I, classical class I, class III, classical class II, and extended class II regions from telomere (left and top side) to centromere (right and bottom side). White, gray, striped and black boxes show expressed genes, gene candidates, non-coding genes and pseudogenes, respectively. The location of the alpha, beta, and kappa blocks containing the cluster of duplicated HLA class I genes in the class I region are indicated. Reprinted with permission from Shiina et al. (2009).
HLA Class I ASD Associations
There are numerous ASD associations with alleles/genes in the HLA region. For example, three different research groups have demonstrated a significant association between ASD and the HLA-A2 (A2) allele (Ferrante et al., 2003; Torres et al., 2006; Al-Hakbany et al., 2014). Gourraud et al. (2014) noted that there are over 100,000 SNPs in the 4.5 million base pairs of HLA (1 out of every 45 bases) and concluded that whole genome screening methods like GWAS are not suited for interrogating HLA polymorphisms. They referred to this extreme HLA diversity as a “genome within a genome.” Table 1 illustrates the complexity of SNPs in a small region of the HLA-A-locus. Site specific primers (SSP) PCR assays clearly and inexpensively identify these complex HLA polymorphisms (Bunce and Passey, 2013). Other HLA class I and II loci have similar SNP complexity as do KIR genes, but none of these sets of polymorphisms are readily detected by GWAS techniques.
It has recently been reported in the Thai population that certain HLA-B alleles (B13:02, B35:02, B44:03, and B56:01) are significantly associated with ASD while two others (B18:02 and B46:12) are associated in a protective manner (Puangpetch et al., 2015). This is not surprising as certain HLA extended haplotypes which have about 200 genes including specific HLA-A, B, and C alleles have very strong associations with ASD. For example, in the Caucasian population HLA-B44 and B15 alleles as part of extended haplotypes 44.1 and 62.1, respectively, are associated with ASD (Torres et al., 2012a).
HLA-G non-classical I molecules are expressed during pregnancy on trophoblast cells that interact with leukocyte-associated immunoglobulin-like receptor (LAIR) and KIR (2DL4) molecules expressed by maternal NK-cells at the uterine fetal/maternal interface to suppress normal immune responses (Tilburgs et al., 2015). A 14 bp deletion in the HLA-G 3′-UTR is associated with higher levels of HLA-G expression whereas a 14 bp insertion associates with reduced HLA-G levels. Guerini et al. (2015) have shown that ASD subjects have an increase in the 14 bp insertion and lower levels of soluble HLA-G protein. This suggests that ASD subjects and their mothers have less HLA-G-mediated immune tolerance during pregnancy. Human non-classical HLA I molecules could be important in brain development as it has recently been shown that non-classical MHC class I molecules are expressed in the olfactory bulb, hippocampus, cerebellum, and nerve nuclei in the developing embryonic mouse brain (Liu et al., 2015).
A main function of HLA class I molecules is to present antigens to CD8+ T-cells, thus starting the complex immune process that creates cytotoxic T-cells to attack specific targets. Any gene in an extended haplotype may be the candidate ASD gene and other HLA alleles may be passengers, however, the identification of foreign peptides that bind to A2 could suggest microorganisms that may be involved in the etiology of ASD. Our data suggests that the HLA class I proteins cognate ligands for KIR are important in autism and that other NK-cell killing receptors could also be important (Figure 4). This interaction does not involve peptide binding as KIR receptors recognize specific amino acid sequences of certain HLA cognate protein ligands (Parham and Moffett, 2013). HLA class I A, B, C alleles all behave as antigen presenting ligands for self and nonself-peptides; however, only certain HLA-A, B, C alleles are ligands for KIR cell surface proteins. Although the HLA cognate ligand site slightly overlaps the peptide binding site, it has a separate function (Figure 3). Specific HLA ligands that bind to KIR activating receptors have been shown to be increased in ASD (Torres et al., 2012b, 2016; Guerini et al., 2014). It should be mentioned that the HLA-A2 does not behave as a KIR ligand.
Figure 3. Demonstrates the overlapping regions for KIR ligand binding and antigen-presenting portions of the HLA protein.
Figure 4. HLA gene complex encode ligands that react with receptors on the Leukocyte Receptor Complex (LCR) and Natural Killer Cell Gene Complex (NKC) to influence killing. Blue lines indicate areas of completed research. HLA-A, B, C ligands bind to various KIR receptors encoded on chromosome 19. The DRB1*0401 shared epitope (SE*) is strongly associated with ASD (RR = 19.8) and helps increase TH17 cells. The dashed lines suggest future ASD research to get a more complete picture on NK-cell killing and antibody production. HLA-DR alleles bind to CD4+ helper T-cells to influence B-cells. HLA-E and MHC Class I like- sequence A (MICA) proteins are ligands for activating molecules NKG2(3) and NKG2D, respectively, that influence NK-cell killing (chromosome 12).
It was only a few years ago that the brain was considered immune privileged, being protected from immune fluctuations by the blood brain barrier. Today there is increasing evidence that the systemic immune system is involved in CNS functioning. Pro-inflammatory cytokines such as TNF-α can easily cross the blood brain barrier and complement proteins such as C1q and C4 are active in brain function (Johnson et al., 1994; Torres et al., 2001). C1q appears to be especially important in eliminating unwanted synapses (Stevens et al., 2007). Perhaps the most important research from an immune-brain perspective is the involvement of HLA class I proteins in brain development. The expression of HLA class I proteins occurs throughout the brain on neurons and even synaptic membranes (Elmer et al., 2013). It is yet unknown if there are differences in class I alleles in brain development and this is an area clearly deserving of further research.
HLA Class II ASD Associations
A maternal-fetal immune interaction has been suggested by Lee et al. (2006) in which boys with autism and their mothers were shown to have an increased DR4 frequency over the control subjects (odds ratios 4.20 and 5.54, respectively). Johnson et al. (2009) later reported significant transmission disequilibrium for HLA-DR4 (odds ratio 4.67) from maternal grandparents to mothers of children with autism, more evidence that HLA-DR4 has a maternal-fetal interaction. Finally, it has been shown in Han Chinese that the HLA-DRB1 allelic frequencies including DR4 are different in autism subjects vs. control subjects (Chien et al., 2011). It should be noted that it is common for racial groups to have different HLA alleles associated with autoimmune diseases.
Warren et al. (1996) also noted that the third hypervariable region (HVR-3) of certain class II DR4 alleles, referred to as the shared epitope (SE), were strongly associated with ASD (p < 0.01; relative risk = 19.8). In this review, odds ratio and relative risk, two similar mathematical models to examine disease risk, are considered equivalent. The HLA-DRB1 SE is a 5 amino acid motif shared by 5 DRB1 alleles. The SE also strongly associates with the majority of severe rheumatoid arthritis (RA) patients. The basis for the SE association with ASD and RA is unknown; however, it has been proposed that the SE peptide acts as a signal transduction ligand that activates nitric oxide (NO) and reactive oxygen species production (de Almeida et al., 2011). It is important to mention that De Almeida et al. (2010) concluded, in an earlier publication, that the HLA DRB1 shared epitope peptide is a potent immune-stimulatory ligand that polarizes naive helper T-cell toward the potent inflammatory TH17 lineage resulting in higher IL-17 levels detected in autism subjects (Al-Ayadhi and Mostafa, 2012; Onore et al., 2012).
It should be noted that mothers with SE alleles on either chromosome were more likely to give birth to an ASD child, referred to by Warren et al. (1996) as a “maternal attack against fetal tissue.” The SE observation is important as it suggested that small peptide epitopes contributed to the high relative risk in ASD years before the TH17 cell observation was made.
Autoantibodies in Autism
One of the most interesting areas of current ASD research is the observation in at least 8 studies that up to about 10% of mothers with ASD children and only 0–2% of controls have humoral antibodies against fetal brain proteins (Croen et al., 2008; Braunschweig and Van de Water, 2012). Because HLA class II molecules are important in antibody production it would be interesting to know if there is an association between mothers who make these antibodies and certain DRB1 alleles. The entire area of maternal antibodies against fetal proteins, the production of autoantibodies, as well as antibodies to microbial pathogens should be studied more in ASD (Grether et al., 2016).
Class III ASD Associations
Four of the 25 proteins in the complement system that help or “complement” antibodies and phagocytic cells to clear pathogens from the organism are encoded within the HLA complex (C4A, C4B, C2, and Bf). The C4 complement proteins (C4A and C4B) are encoded by two separate genes in the HLA class III region (Torres et al., 2001). It was reported 25 years ago that subjects with autism had a significant increase in the deletion of the HLA class III C4B gene (C4B null allele) compared to control subjects (Warren et al., 1991; Odell et al., 2005). ASD subjects were also noted to have a significant deficiency in the plasma C4B protein as determined by an Elisa assay (Warren et al., 1994). Mostafa and Shehab (2010) have also reported a significant increase in the deletion of the C4B gene in the Egyptian ASD population. However, the gene copy number data from qPCR assays do not entirely agree with the results from the C4B protein assay, a not uncommon occurrence when comparing genomic and protein expression data (Kendrick, 2014), thus the solidity of these results are still uncertain. While the data of Warren et al. (1996) suggest that the C4B null allele is another common variant, as it is part of the 44.1 extended haplotype that also contains the SE, further clarification of the association of C4B null alleles is required, preferably using samples with known C4A/C4B DNA complotypes.
There are four proinflammatory cytokines: tumor necrosis factor-alpha (TNFα), lymphotoxin alpha and beta (LTA and LTB), and leukocyte specific transcript-1 (LST1) in the class III region that are important in inflammation and infection. TNFα is a proinflammatory, multifunctional cytokine that is synthesized in numerous cells including NK-cells, macrophages, granulocytes and T-cells. In the CNS, TNFα is made in NK-cells, microglia (brain macrophages), astrocytes, and neurons and is necessary for neural cell differentiation and neuron maturation and may be critical for proper synaptic pruning (Cacci et al., 2005). Elevated levels of TNFα are present in ASD patients, as well as in numerous neurological disorders including multiple sclerosis, Alzheimer's Disease, Parkinson's Disease, ischemia, and traumatic brain injury (Montgomery and Bowers, 2012). The three adjacent proinflammatory cytokine genes: lymphotoxin alpha and beta (LTA and LTB) and leukocyte specific transcript-1 (LST1) have not been studied as extensively as TNFα (Ovsyannikova et al., 2010) and deserve further investigation.
There is a developing consensus that there are elevated plasma levels of proinflammatory cytokines such as IL-1β, IL-6, IL-8, as well as CCL2 and CCL5 chemokines in ASD (Grigorenko et al., 2008; Ashwood et al., 2011). Decreased levels of anti-inflammatory cytokines, most notably TGF-β and IL-10, have also been reported in ASD (Okada et al., 2007; Abdallah et al., 2013; Jyonouchi et al., 2014). In addition to these proinflammatory cytokines, it is interesting to note that an increase in interferon-gamma (INFγ) and TNFα has been observed in brain tissue and cerebral spinal fluid of in ASD subjects (Li et al., 2009). These two cytokines are of special interest, as TNFα is in the HLA class III region and activates NK-cells to make large amounts of INFγ (Schoenborn and Wilson, 2007). TNFα also may be important in ASD as it appears to be a key mediator in neuroinflammation and blood-brain barrier deterioration (McCoy and Tansey, 2008). TNFα concentrations are also elevated in several brain disorders including multiple sclerosis, Alzheimer's Disease, Parkinson's Disease, ischemia, and traumatic brain injury (Montgomery and Bowers, 2012). It is unclear if TNFα contributes to the disease state or if the higher concentrations limit brain damage. Finally, another avenue for investigation, is the class III gene encoding the heat shock protein 70 (HSP70) which is induced in many CNS disorders such as stroke, epilepsy, trauma and ASD (El-Ansary and Al-Ayadhi, 2012).
INFγ may contribute to ASD pathology by increasing the enzyme nitric oxide synthase in astrocytes and microglia, and thereby increasing brain levels of NO. It is thought that activated immune cells secrete high levels of the free radical NO to damage pathogens. In the brain however, NO is typically produced by neurons and it acts as an intercellular messenger modulating synaptogenesis, dendrite, and axonal growth and neuronal release of various neurotransmitters (Hess et al., 1993; Lizasoain et al., 1996). Problems with synaptic development and plasticity are implicated in the neuropathology of ASD (Ebert and Greenberg, 2013).
It was reported a couple of years ago that serum levels of IL-17 were elevated in ASD subjects over control subjects suggesting that T helper 17 (TH17) cells are involved in the etiology of ASD (Al-Ayadhi and Mostafa, 2012; Onore et al., 2012). TH17 cells are CD4+ helper T-cells characterized by a unique cytokine profile, mainly high levels of IL-17 thought to be critical for the development of autoimmunity (Bedoya et al., 2013).
The Importance of HLA Immunity in Autism
The HLA gene complex has been considered by many to be the most important region for immune function in the human genome. However, the system is not perfect and certain HLA allotypes and haplotypes have been shown to be strongly associated with a wide spectrum of autoimmune diseases, such as Type I diabetes, celiac disease, psoriasis, rheumatoid arthritis (RA), lupus erthymatosis, myasthenia gravid, Sjogren Syndrome, and now ASD. The HVR-3 (or shared epitope) that is increased in ASD (Warren et al., 1996) is also associated with several autoimmune diseases such as RA, psoriatic arthritis, and systemic lupus erythematosus (de Almeida et al., 2011). Various HLA disease associations often have the relative risk factors 15 to 20-fold higher in subjects with autoimmune diseases than control subjects (Feitsma et al., 2008). The HVR-3 (DRB1*0401) association with ASD had a comparable relative risk of 19.8 and 31.5% increase in the ASD population relative to control subjects, fitting the criteria for a “common variant–common disease” association.
HLA/KIR Common Genetic Variants
There is currently an intense interest in identifying genetic variants that increase the risk of developing ASD. The goal of this paper is to summarize data which suggests that certain HLA/KIR immune alleles/genes fit the criteria as common genetic variants. For example, the HLA-A2 allele in independent studies from three different groups has a >10% increased frequency in the ASD population than controls (Table 2A). Interestingly, it was over three decades ago that Stubbs and Magenis (1980) first implicated HLA-A2 in ASD. A 14 base pair insertion in the class I non-classical HLA-G loci is strongly associated with ASD (Guerini et al., 2015) and is clearly a common variant at 18% over controls (Table 2A).
Table 2. (A) Suggests that HLA class I and class II alleles are common genetic variants (>5%). (B) Suggests that KIR genes and KIR gene-content haplotypes are common genetic variants.
The DR4 allele has higher frequencies of 18.94 and 7.64% in two ASD populations over control populations suggesting that DR4 is a common variant (Table 2A). Warren et al. (1996) presented data that strongly suggested that the shared epitope is a common genetic variant with a 31.5% increase of ASD over control subjects (Table 2A).
Examination of KIR gene frequencies suggests that certain activating genes are common variants with ASD frequency differentials well above 5% in two different publications (Torres et al., 2012b; Guerini et al., 2014). The KIR region in the leukocyte receptor complex (LRC) is not typical as KIR haplotypes are based on different combination of activating and inhibiting genes that influence NK-cell killing (Pyo et al., 2013; Torres et al., 2016). Gene frequency analyses suggest that three activating genes (3DS1, 2DS1, and 2DS2) are common variants (Table 2B). Additionally, examinations of all activating KIR genes (aKIR) indicate an increase in total activating genes in ASD over control subjects. KIR gene-content haplotypes have centromeric and telomeric parts that join at the junction between 3DP1 and 2DL4 (Figure 5). Examination of the KIR haplotypes suggests that three partial KIR haplotypes (tA01/cA01/cB01) and a complete KIR gene-content haplotype (cB01/tA01) all meet common variant criteria (Table 2B; Torres et al., 2016).
Figure 5. Genes present in KIR gene-content haplotypes. Only genes present in a haplotype are listed under the gene column: for example, 2DS2 is absent in the first two haplotypes, but present in the last four haplotypes. Certain haplotype can have 2DL2 or 2DL3 and other haplotype can have 2DL1 or 2DS1 and 2DS1 or 2DS4. The complete gene-content KIR haplotype (cB01/tA01) is a proposed common variant.
Immunogenetics commonly refers to the strong association of HLA alleles/genes with autoimmune diseases (Shiina et al., 2004). We wish to propose that much of the immune dysfunction in ASD has a basis in the HLA /KIR gene complexes; thus, we have coined the term Autism Immunogenetics. Although, many of the multiple thousands of HLA alleles are rare there are hundreds of weIl-characterized HLA class I and class II alleles that make up a significant part of the genetic diversity of the 900 gene Immunome.
One must remember that the entire 200 genes in the HLA gene complex are inherited as large 4.5 million base pair extended haplotypes and a couple of these large extended haplotypes have been associated with autism (Warren et al., 1996). The association of the A2, DR4 alleles and C4B null allele with ASD may be due to the fact that they are part of the extended haplotype 44.1 that is ASD associated (Figure 4; Table 2). However, the important observation is that all are common genetic variants. It would be interesting to know which extended haplotypes have the HLA-G 14 bp insertion as it strongly associates with ASD. The shared epitope peptide may be the most important as it directs helper T-cells to the TH17 lineage which appears to be important in the pathogenesis of autoimmune neuroinflammatory diseases including ASD (Figure 4; Onore et al., 2012; Choi et al., 2016). A recent publication suggested that maternal TH17 cells could induce an ASD-like phenotype in newborn mice and that blocking antibodies to IL-17a ameliorated these symptoms (Choi et al., 2016). The high relative risk for the SE in ASD noted by Warren et al. (1996) may be due to an increase in TH17cells.
A significant body of evidence has accumulated over the past 30 plus years to demonstrate that specific variants in immune system genes, involving both innate and adaptive immunity pathways, are associated with the incidence of ASD. However, such associations have not been detected by genome wide screening methods, e.g., GWAS, and therefore have been widely ignored by the majority of investigators in the autism research field. The inability to detect immune gene ASD associations by GWAS techniques is fully understandable, considering the massive genetic complexity of the “Immunome,” with more than one hundred thousand polymorphisms localized within a small segment (<5%) of the total genome. Thus, highly focused methods, usually PCR or sequence-based typing, are required to reveal disease associations with immune genes or haplotypes, and these can usually be detected using relatively small number of patient and control samples. Next-Generation Sequencing methods to rapidly sequence the HLA locus are well developed for such autism immunogenetic studies (Nelson et al., 2015; Carapito et al., in press). Indeed, as summarized in Table 2, several HLA and KIR gene variants that meet the “common variant–common disease” criteria have been identified using relatively small sample sizes: in the low hundreds. Perhaps the most interesting observation concerning common genetic variants is the 12% increase of the HLA ligand /KIR activating gene complex in autism over the control populations (Table 2). We argue strongly that additional efforts be focused on identifying and characterizing the role of immune genes in the etiology of ASD, particularly using phenotypically well-defined populations of autistic subjects.
The data summarized in this review raises numerous questions in addition to those raised earlier in this paper. For example, how do different immune gene variants influence brain development and function? While it is clear that many HLA proteins are expressed in or on neural cells and they can modulate the activity of downstream effector cells or proteins, such as NK-cells or pro-inflammatory cytokines, details on how they actually function are unknown. The fact that a single amino acid change in certain HLA class I and II proteins can alter their functionality makes such studies extremely challenging. For example, single DNA base differences in HLA-C alleles encode different amino acids which are ligands for different KIR receptors (below) (Torres et al., 2016).
An additional question that needs to be addressed is why do up to 10% of mothers with ASD children have humoral antibodies against fetal brain proteins (compared to 0–2% of control mothers) and do these autoantibodies adversely affect fetal brain development? Furthermore, does the 14 bp insertion in the HLA G gene that occurs more frequently in ASD subjects and their mothers impair immune tolerance during pregnancy and increase the risk of in utero developmental problems?
HLA Alleles and KIR Genes As Common Genetic Variants in ASD
The observations that several HLA alleles as well as several KIR genes and haplotypes appear to be common genetic variants associated with ASD (Table 2) raises the question as to whether or not other immune genes (for example HLA-E, HLA-F, and MICA) have similar ASD associations. These issues should be addressed in future studies. In summary, by further understanding how immune gene variants participate in the etiology of ASD, it may be possible to: #1 develop biological markers to predict ASD at an earlier stage or even in utero and #2 develop targeted pharmaceutical molecules such as monoclonal antibodies, decoy peptides, and special nucleic acid molecules against SE, HLA, and KIR molecules.
ART conceived of the study, participated in its design and coordination, and drafted the manuscript. TS and JW participated in data collection and helped draft the manuscript. PB, DW, CD, AJT, DO, and LC helped draft the manuscript. RJ performed the statistical analyses. MB participated in data collection and statistical analyses. All authors read and approved the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We wish to acknowledge funding from the Center for Persons with Disabilities, Utah State University, the Utah Autism Foundation, the National Institutes of Environmental Health Sciences grant (5RO1 ESO16669). We also wish to honor Dr. Reed P. Warren who died in 1998 for encouraging several of us (ART, TS, DO, JW) to pursue the study of immune genes in autism. This project has been funded in whole or in part with federal funds from the National Cancer Institute, National Institutes of Health, under contract HHSN261200800001E. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the US Government. This research was supported [in part] by the Intramural Research Program of NIH, National Cancer Institute, Center for Cancer Research.
Abdallah, M. W., Larsen, N., Grove, J., Bonefeld-Jørgensen, E. C., Nørgaard-Pedersen, B., Hougaard, D. M., et al. (2013). Neonatal chemokine levels and risk of autism spectrum disorders: findings from a Danish historic birth cohort follow-up study. Cytokine 61, 370–376. doi: 10.1016/j.cyto.2012.11.015
Al-Hakbany, M., Awadallah, S., and Al-Ayadhi, L. (2014). The relationship of HLA Class I and II alleles and haplotypes with autism: a case control study. Autism Res. Treat. 2014:242048. doi: 10.1155/2014/242048
Ashwood, P., Krakowiak, P., Hertz-Picciotto, I., Hansen, R., Pessah, I. N., and Van de Water, J. (2011). Associations of impaired behaviors with elevated plasma chemokines in autism spectrum disorders. J. Neuroimmunol. 232, 196–199. doi: 10.1016/j.jneuroim.2010.10.025
Bailey, A., Le Couteur, A., Gottesman, I., Bolton, P., Simonoff, E., Yuzda, E., et al. (1995). Autism as a strongly genetic disorder: evidence from a british twin study. Psychol. Med. 25, 63–77. doi: 10.1017/S0033291700028099
Bolton, P., Macdonald, H., Pickles, A., Rios, P., Goode, S., Crowson, M., et al. (1994). A case-control family history study of autism. J. Child Psychol. Psychiatry 35, 877–900. doi: 10.1111/j.1469-7610.1994.tb02300.x
Cacci, E., Claasen, J.-H., and Kokaia, Z. (2005). Microglia-derived tumor necrosis factor-alpha exaggerates death of newborn hippocampal progenitor cells in vitro. J. Neurosci. Res. 80, 789–797. doi: 10.1002/jnr.20531
Carapito, R., Radosavljevic, M., and Bahram, S. (in press). Next-generation sequencing of the HLA locus: methods and impacts on HLA typing, population genetics and disease association studies. Hum. Immunol. doi: 10.1016/j.humimm.2016.04.002
Chien, Y. L., Wu, Y. Y., Chen, C. H., Gau, S. S., Huang, Y. S., Chien, W. H., et al. (2011). Association of HLA-DRB1 alleles and neuropsychological function in autism. Psychiatr. Genet. 22, 46–49. doi: 10.1097/YPG.0b013e32834915ae
Choi, G. B., Yim, Y. S., Wong, H., Kim, S., Kim, H., Kim, S. V., et al. (2016). The maternal interleukin-17a pathway in mice promotes autism-like phenotypes in offspring. Science 351, 933–939. doi: 10.1126/science.aad0314
Constantino, J. N., Todorov, A., and Hilton, C. (2013). Autism recurrence in half siblings: strong support for genetic mechanisms of transmission in ASD. Mol. Psychiatry 18, 137–138. doi: 10.1038/mp.2012.9
Croen, L. A., Braunschweig, D., Haapanen, L., Yoshida, C. K., Fireman, B., Grether, J. K., et al. (2008). Maternal mid-pregnancy autoantibodies to fetal brain protein: the early markers for autism study. Biol. Psychiatry 64, 583–588. doi: 10.1016/j.biopsych.2008.05.006
De Almeida, D. E., Ling, S., Pi, X., Hartmann-Scruggs, A. M., Pumpens, P., and Holoshitz, J. (2010). Immune dysregulation by the rheumatoid arthritis shared epitope. J. Immunol. 185, 1927–1934. doi: 10.4049/jimmunol.0904002
Elmer, B. M., Estes, M. L., Barrow, S. L., and McAllister, A. K. (2013). MHCI requires MEF2 transcription factors to negatively regulate synapse density during development and in disease. J. Neurosci. 33, 13791–13804. doi: 10.1523/JNEUROSCI.2366-13.2013
Feitsma, A. L., van der Helm-van Mil, A. H., Huizinga, T. W., de Vries, R. R., and Toes, R. E. (2008). Protection against rheumatoid arthritis by HLA: nature and nurture. Ann. Rheum. Dis. 67(Suppl. 3), 61–63. doi: 10.1136/ard.2008.098509
Ferrante, P., Saresella, M., Guerini, F. R., Marzorati, M., Musetti, M. C., and Cazzullo, A. G. (2003). Significant association of HLA A2-DR11 with CD4 naive decrease in autistic children. Biomed. Pharmacother. 57, 372–374. doi: 10.1016/S0753-3322(03)00099-4
Gaugler, T., Klei, L., Sanders, S. J., Bodea, C. A., Goldberg, A. P., Lee, A. B., et al. (2014). Most genetic risk for autism resides with common variation. Nat. Genet. 46, 881–885. doi: 10.1038/ng.3039
Grether, J. K., Ashwood, P., Van de Water, J., Yolken, R. H., Anderson, M. C., Torres, A. R., et al. (2016). Prenatal and newborn immunoglobulin levels from mother-child pairs and risk of autism spectrum disorders. Front. Neurosci. 10:218. doi: 10.3389/fnins.2016.00218
Grigorenko, E. L., Han, S. S., Yrigollen, C. M., Leng, L., Mizue, Y., Anderson, G. M., et al. (2008). Macrophage migration inhibitory factor and autism spectrum disorders. Pediatrics 122, E438–E445. doi: 10.1542/peds.2007-3604
Guerini, F. R., Bolognesi, E., Chiappedi, M., Ghezzo, A., Canevini, M. P., Mensi, M. 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. doi: 10.1016/j.bbi.2014.10.002
Guerini, F. R., Bolognesi, E., Chiappedi, M., Manca, S., Ghezzo, A., Agliardi, C., et al. (2014). Activating KIR molecules and their cognate ligands prevail in children with a diagnosis of ASD and in their mothers. Brain Behav. Immun. 36, 54–60. doi: 10.1016/j.bbi.2013.10.006
Hess, D. T., Patterson, S. I., Smith, D. S., and Skene, J. H. (1993). Neuronal growth cone collapse and inhibition of protein fatty acylation by nitric oxide. Nature 366, 562–565. doi: 10.1038/366562a0
Horton, R., Gibson, R., Coggill, P., Miretti, M., Allcock, R. J., Almeida, J., et al. (2008). Variation analysis and gene annotation of eight MHC haplotypes: the MHC Haplotype Project. Immunogenetics 60, 1–18, doi: 10.1007/s00251-007-0262-2
Irimia, M., Weatheritt, R. J., Ellis, J. D., Parikshak, N. N., Gonatopoulos-Pournatzis, T., Babor, M., et al. (2014). A highly conserved program of neuronal microexons is misregulated in autistic brains. Cell 159, 1511–1523. doi: 10.1016/j.cell.2014.11.035
Johnson, W. G., Buyske, S., Mars, A. E., Sreenath, M., Stenroos, E. S., Williams, T. A., et al. (2009). HLA-DR4 as a risk allele for autism acting in mothers of probands possibly during pregnancy. Arch. Pediatr. Adolesc. Med. 163, 542–546. doi: 10.1001/archpediatrics.2009.74
Jyonouchi, H., Geng, L., and Davidow, A. L. (2014). Cytokine profiles by peripheral blood monocytes are associated with changes in behavioral symptoms following immune insults in a subset of ASD subjects: an inflammatory subtype? J. Neuroinflammation 11, 187. doi: 10.1186/s12974-014-0187-2
Laver, T. W., Caswell, R. C., Moore, K. A., Poschmann, J., Johnson, M. B., Owens, M. M., et al. (2016). Pitfalls of haplotype phasing from amplicon-based long-read sequencing. Sci. Rep. 6:21746. doi: 10.1038/srep21746
Lee, L. C., Zachary, A. A., Leffell, M. S., Newschaffer, C. J., Matteson, K. J., and Tyler, J. D. (2006). HLA-DR4 in families with autism. Pediatr. Neurol. 35, 303–307. doi: 10.1016/j.pediatrneurol.2006.06.006
Li, X., Chauhan, A., Sheikh, A. M., Patil, S., Chauhan, V., Li, X. M., et al. (2009). Elevated immune response in the brain of autistic patients. J. Neuroimmunol. 207, 111–116. doi: 10.1016/j.jneuroim.2008.12.002
Liu, J., Shen, Y., Li, M., Lv, D., Zhang, A., Peng, Y., et al. (2015). Spatial-temporal expression of non-classical MHC Class I molecules in the C57 mouse Brain. Neurochem. Res. 40, 1487–1496. doi: 10.1007/s11064-015-1620-3
Lizasoain, I., Weiner, C. P., Knowles, R. G., and Moncada, S. (1996). The ontogeny of cerebral and cerebellar nitric oxide synthase in the guinea pig and rat. Pediatr. Res. 39, 779–783. doi: 10.1203/00006450-199605000-00006
Manolio, T. A., Collins, F. S., Cox, N. J., Goldstein, D. B., Hindorff, L. A., Hunter, D. J., et al. (2009). Finding the missing heritability of complex diseases. Nature 461, 747–753. doi: 10.1038/nature08494
McCoy, M. K., and Tansey, M. G. (2008). TNF signaling inhibition in the CNS: implications for normal brain function and neurodegenerative disease. J. Neuroinflammation 5:45. doi: 10.1186/1742-2094-5-45
McQueen, K. L., Dorighi, K. M., Guethlein, L. A., Wong, R., Sanjanwala, B., and Parham, P. (2007). Donor-recipient combinations of group A and B KIR haplotypes and HLA class I ligand affect the outcome of HLA-matched, sibling donor hemotopoietic cell transplantation. Hum. Immunol. 68, 309–323. doi: 10.1016/j.humimm.2007.01.019
Montgomery, S. L., and Bowers, W. J. (2012). Tumor necrosis factor-α and the roles it plays in homeostatic and degenerative processes within the central nervous system. J. Neuroimmune Pharmacol. 7, 42–59. doi: 10.1007/s11481-011-9287-2
Mostafa, G. A., and Shehab, A. A. (2010). The link of C4B null allele to autism and to a family history of autoimmunity in Egyptian autistic children. J. Neuroimmunol. 223, 115–119. doi: 10.1016/j.jneuroim.2010.03.025
Nelson, W. C., Pyo, C. W., Vogan, D., Wang, R., Pyon, Y. S., Hennessey, C., et al. (2015). An integrated genotyping approach for HLA and other complex genetic systems. Hum. Immunol. 76, 928–938. doi: 10.1016/j.humimm.2015.05.001
Odell, D., Maciulis, A., Cutler, A., Warren, L., McMahon, W. M., Coon, H., et al. (2005). Confirmation of the association of the C4B null allele in autism. Hum. Immunol. 66, 140–145. doi: 10.1016/j.humimm.2004.11.002
Okada, K., Hashimoto, K., Iwata, Y., Nakamura, K., Tsujii, M., Tsuchiya, K. J., et al. (2007). Decreased serum levels of transforming growth factor-beta1 in patients with autism. Prog. Neuropsychopharmacol. Biol. Psychiatry 31, 187–190. doi: 10.1016/j.pnpbp.2006.08.020
Ovsyannikova, I. G., Vierkant, R. A., Pankratz, V. S., Jacobson, R. M., and Poland, G. A. (2010). Extended LTA, TNF, LST1 and HLA gene haplotypes and their association with rubella vaccine-induced immunity. PLoS ONE 5:e1106. doi: 10.1371/journal.pone.0011806
Philippe, A., Martinez, M., Guilloud-Bataille, M., Gillberg, C., Råstam, M., Sponheim, E., et al. (1999). Genome-wide scan for autism susceptibility genes. Paris autism research international sibpair study. Hum. Mol. Genet. 8, 805–812. doi: 10.1093/hmg/8.5.805
Puangpetch, A., Suwannarat, P., Chamnanphol, M., Koomdee, N., Ngamsamut, N., Limsila, P., et al. (2015). Significant association of HLA-B alleles and genotypes in Thai children with autism spectrum disorders: a case-control study. Dis. Markers 2015:724935. doi: 10.1155/2015/724935
Pyo, C.-W., Wang, R., Vu, Q., Cereb, N., Yang, S. Y., Duh, F.-M., et al. (2013). Recombinant structures expand and contract inter and intragenic diversification at the KIR locus. BMC Genomics 14:89. doi: 10.1186/1471-2164-14-89
Shiina, T., Inoko, H., and Kulski, J. K. (2004). An update of the HLA genomic region, locus information and disease associations: 2004. Tissue Antigens 64, 631–649. doi: 10.1111/j.1399-0039.2004.00327.x
Stevens, B., Allen, N. J., Vazquez, L. E., Howell, G. R., Christopherson, K. S., Nouri, N., et al. (2007). The classical complement cascade mediates CNS synapse elimination. Cell 131, 1164–1178. doi: 10.1016/j.cell.2007.10.036
Tao, Y., Gao, H., Ackerman, B., Guo, W., Saffen, D., and Shugart, Y. Y. (2016). Evidence for contribution of common genetic variants within chromosome 8p21.2-8p21.1 to restricted and repetitive behaviors in autism spectrum disorders. BMC Genomics 17:163. doi: 10.1186/s12864-016-2475-y
Tilburgs, T., Evans, J. H., Crespo, Â. C., and Strominger, J. L. (2015). The HLA-G cycle provides for both NK tolerance and immunity at the maternal–fetal interface. Proc. Natl. Acad. Sci. U.S.A. 112, 13312–13317. doi: 10.1073/pnas.1517724112
Torres, A., Westover, J., Benson, M., Johnson, R., and Dykes, A. (2016). A killer immunoglobulin-like receptor gene-content haplotype and a cognate human leukocyte antigen ligand are associated with autism. Autism Open Access. 6:171. doi: 10.4172/2165-7890.1000171
Torres, A. R., Maciulis, A., Stubbs, E. G., Cutler, A., and Odell, D. (2002). The transmission disequilibrium test suggests that HLA-DR4 and DR13 are linked to autism spectrum disorder. Hum. Immunol. 63, 311–316. doi: 10.1016/S0198-8859(02)00374-9
Torres, A. R., Sweeten, T. L., Cutler, A., Bedke, B. J., Fillmore, M., Stubbs, E. G., et al. (2006). The association and linkage of the HLA-A2 class I allele with autism. Hum Immunol. 67, 346–351. doi: 10.1016/j.humimm.2006.01.001
Torres, A. R., Westover, J. B., Gibbons, C., Johnson, R. C., and Ward, D. C. (2012b). Activating killer-cell immunoglobulin-like receptors (KIR) and their cognate HLA ligands are significantly increased in autism. Brain Behav. Immun. 26, 1122–1127. doi: 10.1016/j.bbi.2012.07.014
Warren, R. P., Burger, R. A., Odell, D., Torres, A. R., and Warren, W. L. (1994). Decreased plasma concentrations of the C4B complement protein in autism. Arch. Pediatr. Adolesc. Med. 148, 180–183. doi: 10.1001/archpedi.1994.02170020066011
Warren, R. P., Odell, J. D., Warren, W. L., Burger, R. A., Maciulis, A., Daniels, W. W., et al. (1996). Strong association of the third hypervariable region of HLA-DRB1 with autism. J. Neuroimmunol. 67, 97–102. doi: 10.1016/0165-5728(96)00052-5
Warren, R. P., Singh, V. K., Cole, P., Odell, J. D., Pingree, C. B., Warren, W. L., et al. (1991). Increased frequency of the null allele at the complement C4B locus in autism. Clin. Exp. Immunol. 83, 438–440. doi: 10.1111/j.1365-2249.1991.tb05657.x
Westover, J. B., Sweeten, T. L., Benson, M., Bray-Ward, P., and Torres, A. R. (2011). “Immune dysfunction in autism spectrum disorder,” in Autism: A Neurodevelopmental Journey from Genes to Behavior, ed V. Eapen (Rijeka: Intech), 343–368.
Young, A. M., Chakrabarti, B., Roberts, D., Lai, M. C., Suckling, J., and Baron-Cohen, S. (2016). From molecules to neural morphology: understanding neuroinflammation in autism spectrum condition. Mol. Autism 7:9. doi: 10.1186/s13229-016-0068-x
Keywords: autism, HLA, KIR, haplotype, common genetic variant
Citation: Torres AR, Sweeten TL, Johnson RC, Odell D, Westover JB, Bray-Ward P, Ward DC, Davies CJ, Thomas AJ, Croen LA and Benson M (2016) Common Genetic Variants Found in HLA and KIR Immune Genes in Autism Spectrum Disorder. Front. Neurosci. 10:463. doi: 10.3389/fnins.2016.00463
Received: 08 June 2016; Accepted: 26 September 2016;
Published: 20 October 2016.
Edited by:Joshua Pesach Rosenzweig, Cell-El, Israel
Reviewed by:Harry W. Schroeder, University of Alabama at Birmingham, USA
William G. Johnson, Rutgers, USA
Copyright © 2016 Torres, Sweeten, Johnson, Odell, Westover, Bray-Ward, Ward, Davies, Thomas, Croen and Benson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Anthony R. Torres, firstname.lastname@example.org