Even though immunoglobulins are critical for immune responses and human survival, the diversity of the immunoglobulin heavy chain gene ( IGH ) is poorly known and mostly characterized only by serological methods. Moreover, this genomic region is not well-covered in genomic databases and genome-wide association studies due to particularities that impose technical difficulties for its analysis. Therefore, the IGH gene has never been systematically sequenced across populations. Here, we deliver an unprecedented and comprehensive characterization of the diversity of the IGHG1, IGHG2 , and IGHG3 gene segments, which encode the constant region of the most abundant circulating immunoglobulins: IgG1, IgG2, and IgG3, respectively. We used Sanger sequencing to analyze 357 individuals from seven different Brazilian populations, including five Amerindian, one Japanese-descendant and one Euro-descendant population samples. We discovered 28 novel IGHG alleles and provided evidence that some of them may have been originated by gene conversion between common alleles of different gene segments. The rate of synonymous substitutions was significantly higher than the rate of the non-synonymous substitutions for IGHG1 and IGHG2 ( p = 0.01 and 0.03, respectively), consistent with purifying selection. Fay and Wu's test showed significant negative values for most populations ( p < 0.001), which indicates that positive selection in an adjacent position may be shaping IGHG variation by hitchhiking of variants in the vicinity, possibly the regions that encode the Ig variable regions. This study shows that the variation in the IGH gene is largely underestimated. Therefore, exploring its nucleotide diversity in populations may provide valuable information for comprehension of its evolution, its impact on diseases and vaccine research.

Regulation of NK cell activity is mediated through killer-cell immunoglobulin-like receptors (KIR) ability to recognize human leukocyte antigen (HLA) class I molecules as ligands. Interaction of KIR and HLA is implicated in viral infections, autoimmunity, and reproduction and there is growing evidence of the coevolution of these two independently segregating gene families. By leveraging KIR and HLA-C data from 1000 Genomes consortium we observed that the KIR2DL1 variant rs2304224 * T is associated with lower expression of HLA-C in individuals carrying the ligand HLA-C2 (p = 0.0059). Using flow cytometry, we demonstrated that this variant is also associated with higher expression of KIR2DL1 on the NK cell surface (p = 0.0002). Next, we applied next generation sequencing to analyze KIR2DL1 sequence variation in 109 Euro and 75 Japanese descendants. Analyzing the extended haplotype homozygosity, we show signals of positive selection for rs4806553 * G and rs687000 * G, which are in linkage disequilibrium with rs2304224 * T. Our results suggest that lower expression of HLA-C2 ligands might be compensated for higher expression of the receptor KIR2DL1 and bring new insights into the coevolution of KIR and HLA.

No abstract

The human leukocyte antigen (HLA) are the most polymorphic genes in the human genome. Because of their importance for antigen recognition, HLA molecules play a central role in host defense and graft rejection upon transplantation. The aim of this study was to characterize allelic diversity of the classical HLA genes HLA‐A, ‐B, ‐C, ‐DRA, ‐DRB1, ‐DQA1, ‐DQB1, ‐DPA1, ‐DPB1, and the non‐classical class I genes HLA‐E, ‐F and ‐G at high‐resolution for a population of predominantly European ancestry from Curitiba, Brazil. Genotyping of 108 individuals was performed by next‐generation sequencing on the MiSeq platform and also by Sanger sequencing. The genotype distributions of all loci were in accordance with Hardy‐Weinberg equilibrium (P > 0.05) and a total of 202 HLA variants at second field resolution were observed for the 12 loci. The strongest linkage disequilibrium (r2 = 1.0, P < 10−5) was observed for the following pairs of alleles: HLA‐B*42:01:01 ~ HLA‐DRB1*03:02:01; HLA‐B*14:02:01 ~ HLA‐C*08:02:01; B*42:01:01 ~ HLA‐C*17:01:01; HLA‐DRB1*03:01:01 ~ HLA‐DQB1*02:01:01 ~ DRB1*03:01:01 ~ HLA‐DQB1*02:01:01; DRB1*13:01:01~ HLA‐DQB1*06:03:01 and HLA‐DRB1*09:01:02 ~ HLA‐DQA1*03:02. This is the first study to characterize all 12 HLA genes at high resolution in a single population. On the basis of the allelic frequencies of worldwide populations and principal component analysis, we confirmed the similarity of the study population to European and other Euro‐descendant populations.

No abstract