Gene Polymorphisms of Novel Immunotolerant Molecule BTLA: Distribution of Alleles, Genotypes and Haplotypes in Polish Caucasian Population

B and T lymphocyte attenuator (BTLA) is one of the members of immunoglobulin superfamily which, like CTLA-4 and PD-1, is involved in down regulation of immune response. Despite the important role of BTLA in maintaining immune homeostasis, relatively little studies were devoted to the relationship of polymorphisms in the gene encoding BTLA with susceptibility to autoimmune disease and cancer. Moreover, all published works were done in Asian populations. BTLA gene is located on chromosome 3 in q13.2 and consists of five exons. The aim of this study was to investigate the alleles, genotypes and haplotypes frequency of selected BTLA gene polymorphisms in Caucasian population originating from Poland. For this study, the single-nucleotide polymorphisms (SNPs) were chosen on the basis of literature data. Additionally, the tag dSNP under linkage equilibrium r 2 > 0.8 and available at the National Center for Biotechnology Information (NCBI) for Caucasian population of rare alleles at a frequency greater than 5 % have been chosen using the NCBI database. The ten BTLA SNPs investigated were: rs1844089, rs2705535, rs9288952, rs9288953, rs1982809, rs2633580, rs2705511, rs2705565, rs76844316, rs16859633. For all SNPs selected on the basis of literature data the significantly different distributions of genotypes between Asian and Caucasian populations were observed.

Abstract B and T lymphocyte attenuator (BTLA) is one of the members of immunoglobulin superfamily which, like CTLA-4 and PD-1, is involved in down regulation of immune response. Despite the important role of BTLA in maintaining immune homeostasis, relatively little studies were devoted to the relationship of polymorphisms in the gene encoding BTLA with susceptibility to autoimmune disease and cancer. Moreover, all published works were done in Asian populations. BTLA gene is located on chromosome 3 in q13.2 and consists of five exons. The aim of this study was to investigate the alleles, genotypes and haplotypes frequency of selected BTLA gene polymorphisms in Caucasian population originating from Poland. For this study, the single-nucleotide polymorphisms (SNPs) were chosen on the basis of literature data. Additionally, the tag dSNP under linkage equilibrium r 2 [ 0.8 and available at the National Center for Biotechnology Information (NCBI) for Caucasian population of rare alleles at a frequency greater than 5 % have been chosen using the NCBI database. The ten BTLA SNPs investigated were: rs1844089, rs2705535, rs9288952, rs9288953, rs1982809, rs2633580, rs2705511, rs2705565, rs76844316, rs16859633. For all SNPs selected on the basis of literature data the significantly different distributions of genotypes between Asian and Caucasian populations were observed. Introduction B and T lymphocyte attenuator (BTLA) is a suppressor molecule belonging to the immunoglobulin superfamily which, like cytotoxic leukocyte antigen-4 (CTLA-4) and programmed death-1 (PD-1), is involved in the inhibition of the immune response. BTLA shares only 9-13 % amino acid identity with CTLA-4 and PD-1, but is structurally similar to them and like other co-inhibitory molecules has two immunoreceptor tyrosine-based inhibitory motifs (ITIM) in its cytoplasmic region (Watanabe et al. 2003). BTLA interaction with the ligand herpes virus entry mediator (HVEM) causes phosphorylation of tyrosine within the ITIM and their connection to protein tyrosine phosphatase SHP-1 and SHP-2, and consequently inhibition of T-cell activation and production of cytokines interferon c, interleukin (IL)-2, IL-4, and IL-10 (Murphy et al. 2006;Nelson et al. 2008;Sedy et al. 2005). BTLAmediated inhibition of T-cell activation occurred during both the primary CD4 ? and secondary CD4 ? and CD8 ? T-cell response, suggesting that BTLA ligation provides a constitutive inhibition signal to T cells (Otsuki et al. 2006). BTLA expression is equal on Th1 and Th2 cells (Otsuki et al. 2006). BTLA expression in murine naive B cells is high, but its role in the activation and differentiation of B cells has not been studied. BTLA expression in naive murine T cells is negligible, but its expression increases during T-cell activation. BTLA expression is also significant in mice thymocytes during positive selection (Watanabe et al. 2003). Moreover, BTLA is expressed by macrophages, dendritic cells and natural killer cells (Han et al. 2004). In recent years strong expression of BTLA on malignant B cells in patients with B-cell chronic lymphocytic leukemia has been reported in the literature (M'Hidi et al. 2009).
Because BTLA was relatively recently described in the literature, there are only a few studies that address BTLA gene polymorphisms, and most have investigated its role in susceptibility to autoimmune diseases, such as rheumatoid arthritis Oki et al. 2011), systemic lupus erythematosus and type 1 diabetes mellitus ) and only one in risk of cancer .
The aim of this study was to determine allele, genotype and haplotype frequencies as well as linkage disequilibrium of selected polymorphisms in the BTLA gene.

Materials and Methods
DNA was isolated from venous blood from 417 unrelated healthy volunteers from Western Poland (172 female/245 male) according to the manual procedure for white blood cells using the QIAamp DNA Blood Mini Kit (Qiagen, Germany).
In Poland migration and nationality diversion is very small, therefore Polish cohort was homogenous. The genetic homogenization was reflected in virtually identical frequencies of H-Y polymorphisms in different regions of Poland described by Ploski et al. (2002).

Statistical Analysis
The evaluation of Hardy-Weinberg equilibrium (HWE) was performed by comparing the observed and expected frequencies of genotypes using v 2 analysis. The haplotype frequencies for pairs of alleles were determined using the SHEsis program (http://202.120.7.14/analysis/myAnalysis. php) (Shi and He 2005). Haplotypes with frequencies lower than 0.03 were not considered. Linkage disequilibrium coefficients r 2 (values for a pair of the most common alleles at each locus) were estimated using http://202.120.7. 14/analysis/myAnalysis.php (Shi and He 2005). The v 2 test was used to compare categorical data between poles and other populations on the basis of literature data. Differences were considered statistically significant if p \ 0.05. For the multiple comparisons, Bonferroni multiple adjustments were employed to the level of significance.
Furthermore, to analyze ethnic relationship between the Polish and other populations [Yoruba in Ibadan, Nigeria; Japanese in Tokyo, Japan; Han Chinese in Beijing, China; CEPH (Utah residents with ancestry from northern and western Europe)], we estimated phylogenetic tree based on nine out of ten single nucleotide polymorphism (SNP) loci (rs76844316 locus was excluded because of lack of data for non-Polish populations) used in present work. The tree was created using SNP allele frequency [for Polish population data were estimated in present study and for others possessed from Hapmap database (Genbank)] with Poptree 2 (Takezaki et al. 2010) using Neigbour Joining algorithm (Saitou and Nei 1987) and D A distance (Nei et al. 1983) as a measure of population differentiation. The tree's nodes significance was evaluated with 1,000 bootstrap analysis.

Results and Discussion
We found that each polymorphism in the BTLA gene was in HWE.
Distributions of alleles and genotypes are presented in Table 1. In five out of ten investigated polymorphisms wild alleles were observed in more than 90 % of the Polish population. In the case of intron polymorphism rs9288953, we observed about 40 % minor allele frequency, while for rs1982809 situated in the 3 0 near gene position and rs2705511 located in the intragenic region the minor allele frequency was about 24 % in both cases. Two other polymorphisms-rs76844316, described in the literature as polymorphic in the Japanese population ; and rs16859633, chosen on the basis of HapMap analysisseem not to be polymorphic in Poles since none of 200 genotyped volunteers were carriers of mutant alleles.
Comparison with data available from HapMap (http:// hapmap.ncbi.nlm.nih.gov/cgi-perl/gbrowse/hapmap24_ B36/#search 2013) for different populations is presented in Table 4. Our results are similar to those reported in Hap-Map for a Caucasian population except for rs2705565. However, distribution of alleles in almost all investigated SNPs varies for Asian populations-Han Chinese and Japanese.
Phylogenetic tree obtained in our study showed clear division into three clades (Fig. 1). First clade was consisted of Polish and Caucasian (CEU) populations, second with Japanese (JPT) and Chinese (CHB) populations and third one with African population (YRI). Obtained results indicated that Polish population is ethnically and genetically closest to Caucasian population (represented in this analysis by CEPH population which was derived from European ancestors). Furthermore, both Caucasian populations were different from Japanese and Chinese populations (which were closest to each other) and all aforementioned populations were distinct from African population.
Arch. Immunol. Ther. Exp. (2015) 63:73-78 75 On the basis of a literature survey Lin et al. 2006;Oki et al. 2011), we compared the genotype frequency for four SNPs (Fig. 2). The SNP rs9288952 was the most commonly investigated and there were data available from Chinese, Taiwanese and Japanese population. This polymorphism is situated in exon 4 and is a missense mutation. The nucleotide exchange causes amino acid residue change Pro-Leu in position 219. The literature data indicate that this SNP is associated with susceptibility to rheumatoid arthritis in the Japanese population and with breast cancer risk in Chinese women. We found statistically significant differences in genotype distribution between Polish and Chinese as well as Polish and Taiwanese populations (Fig. 2). For two other polymorphisms, rs1844089 and rs2705535, connected with breast cancer risk, only data from the Chinese population  were available. We found statistically significant differences in distribution of genotypes between Polish and Chinese populations.
The rs9288953 SNP was investigated in a Japanese population as a candidate SNP for type 1 diabetes mellitus  ). The study showed no associations of that polymorphism with selected autoimmune diseases. In our comparison, we found statistically different distribution of genotypes between Polish and Japanese healthy volunteers. Of note, the genotype distributions reported in this study for controls were not in HWE. Recent papers have indicated the key role of the HVEM-BTLA pathway in inflammation, autoimmunity and infection immunity (del Rio et al. 2010;Shui et al. 2011) as well as immune tolerance, especially in aspects of allotransplantation (Hobo et al. 2012).
On the other hand, it was shown in an animal model that polymorphisms existing in the mouse btla gene were associated with differences in the Ig domain and BTLA expression pattern between different mouse strains and cell lines (Han et al. 2004;Hurchla et al. 2005). Therefore, polymorphisms existing in the human BTLA gene might influence the predisposition to immune-related disease. This work contributes to the overall knowledge of the presence and frequency of occurrence of allelic variants in the Caucasian population precisely in Polish population.