Immunogenetics of HLA-B: SNP, allele, and haplotype diversity in populations from different continents and ancestry backgrounds.

HLA-B is among the most variable gene in the human genome. This gene encodes a key molecule for antigen presentation to CD8+ T lymphocytes and NK cell modulation. Despite the myriad of studies evaluating its coding region (with an emphasis on exons 2 and 3), few studies evaluated introns and regulatory sequences in real population samples. Thus, HLA-B variability is probably underestimated. We applied a bioinformatics pipeline tailored for HLA genes on 5347 samples from 80 different populations, which includes more than 1000 admixed Brazilians, to evaluate the HLA-B variability (SNPs, indels, MNPs, alleles, and haplotypes) in exons, introns, and regulatory regions. We observed 610 variable sites throughout HLA-B; the most frequent variants are shared worldwide. However, the haplotype distribution is geographically structured. We detected 920 full-length haplotypes (exons, introns, and untranslated regions) encoding 239 different protein sequences. HLA-B gene diversity is higher in admixed populations and Europeans while lower in African ancestry individuals. Each HLA-B allele group is associated with specific promoter sequences. This HLA-B variation resource may improve HLA imputation accuracy and disease-association studies and provide evolutionary insights regarding HLA-B genetic diversity in human populations.

Evaluation of MC1R high-throughput nucleotide sequencing data generated by the 1000 Genomes Project

Abstract The advent of next-generation sequencing allows simultaneous processing of several genomic regions/individuals, increasing the availability and accuracy of whole-genome data. However, these new approaches may present some errors and bias due to alignment, genotype calling, and imputation methods. Despite these flaws, data obtained by next-generation sequencing can be valuable for population and evolutionary studies of specific genes, such as genes related to how pigmentation evolved among populations, one of the main topics in human evolutionary biology. Melanocortin-1 receptor (MC1R) is one of the most studied genes involved in pigmentation variation. As MC1R has already been suggested to affect melanogenesis and increase risk of developing melanoma, it constitutes one of the best models to understand how natural selection acts on pigmentation. Here we employed a locally developed pipeline to obtain genotype and haplotype data for MC1R from the raw sequencing data provided by the 1000 Genomes FTP site. We also compared such genotype data to Phase 3 VCF to evaluate its quality and discover any polymorphic sites that may have been overlooked. In conclusion, either the VCF file or one of the presently described pipelines could be used to obtain reliable and accurate genotype calling from the 1000 Genomes Phase 3 data.

Evaluation of MC1R high-throughput nucleotide sequencing data generated by the 1000 Genomes Project.

The advent of next-generation sequencing allows simultaneous processing of several genomic regions/individuals, increasing the availability and accuracy of whole-genome data. However, these new approaches may present some errors and bias due to alignment, genotype calling, and imputation methods. Despite these flaws, data obtained by next-generation sequencing can be valuable for population and evolutionary studies of specific genes, such as genes related to how pigmentation evolved among populations, one of the main topics in human evolutionary biology. Melanocortin-1 receptor (MC1R) is one of the most studied genes involved in pigmentation variation. As MC1R has already been suggested to affect melanogenesis and increase risk of developing melanoma, it constitutes one of the best models to understand how natural selection acts on pigmentation. Here we employed a locally developed pipeline to obtain genotype and haplotype data for MC1R from the raw sequencing data provided by the 1000 Genomes FTP site. We also compared such genotype data to Phase 3 VCF to evaluate its quality and discover any polymorphic sites that may have been overlooked. In conclusion, either the VCF file or one of the presently described pipelines could be used to obtain reliable and accurate genotype calling from the 1000 Genomes Phase 3 data.

HLA-G 3' untranslated region polymorphisms are associated with systemic lupus erythematosus in 2 Brazilian populations.

OBJECTIVE: HLA-G has well recognized tolerogenic properties in physiological and nonphysiological conditions. The 3' untranslated region (3'UTR) of the HLA-G gene has at least 3 polymorphic sites (14-bpINS/DEL, +3142C/G, and +3196C/G) described as associated with posttranscriptional influence on messenger RNA production; however, only the 14-bpINS/DEL and +3142C/G sites have been studied in systemic lupus erythematosus (SLE). METHODS: We investigated the HLA-G 3'UTR polymorphic sites (14-bpINS/DEL, +3003C/T, +3010C/G, +3027A/C, +3035C/T, +3142C/G, +3187A/G, and +3196C/G) in 190 Brazilian patients with SLE and 282 healthy individuals in allele, genotype, and haplotype analyses. A multiple logistic regression model was used to assess the association of the disease features with the HLA-G 3'UTR haplotypes. RESULTS: Increased frequencies were observed of the 14-bpINS (p = 0.053), +3010C (p = 0.008), +3142G (p = 0.006), and +3187A (p = 0.013) alleles, and increased frequencies of the 14-bpINS-INS (p = 0.094), +3010 C-C (p = 0.033), +3142 G-G (p = 0.021), and +3187 A-A (p = 0.035) genotypes. After Bonferroni correction, only the +3142G (p = 0.05) and +3010C (p = 0.06) alleles were overrepresented in SLE patients. The UTR-1 haplotype (14-bpDEL/+3003T/+3010G/+3027C/+3035C/+3142C/+3187G/+3196C) was underrepresented in SLE (pcorr = 0.035). CONCLUSION: These results indicate that HLA-G 3'UTR polymorphic sites, particularly +3142G and +3010C alleles, were associated with SLE susceptibility, whereas UTR-1 was associated with protection against development of SLE.