Hla-C genetic diversity and evolutionary insights in two samples from Brazil and Benin

A multi-ethnic reference panel to impute HLA classical and non-classical class I alleles in admixed samples: Testing imputation accuracy in an admixed sample from Brazil

The MHC class I region contains crucial genes for the innate and adaptive immune response, playing a key role in susceptibility to many autoimmune and infectious diseases. Genome-wide association studies have identified numerous disease-associated SNPs within this region. However, these associations do not fully capture the immune-biological relevance of specific HLA alleles. HLA imputation techniques may leverage available SNP arrays by predicting allele genotypes based on the linkage disequilibrium between SNPs and specific HLA alleles. Successful imputation requires diverse and large reference panels, especially for admixed populations. This study employed a bioinformatics approach to call SNPs and HLA alleles in multi-ethnic samples from the 1000 genomes (1KG) dataset and admixed individuals from Brazil (SABE), utilising 30X whole-genome sequencing data. Using HIBAG, we created three reference panels: 1KG (n = 2504), SABE (n = 1171), and the full model (n = 3675) encompassing all samples. In extensive cross-validation of these reference panels, the multi-ethnic 1KG reference exhibited overall superior performance than the reference with only Brazilian samples. However, the best results were achieved with the full model. Additionally, we expanded the scope of imputation by developing reference panels for non-classical, MICA, MICB and HLA-H genes, previously unavailable for multi-ethnic populations. Validation in an independent Brazilian dataset showcased the superiority of our reference panels over the Michigan Imputation Server, particularly in predicting HLA-B alleles among Brazilians. Our investigations underscored the need to enhance or adapt reference panels to encompass the target population's genetic diversity, emphasising the significance of multiethnic references for accurate imputation across different populations.

Human Leukocyte Antigen Markers for Distinguishing Pustular Psoriasis and Adult-Onset Immunodeficiency with Pustular Reaction

Pustular skin diseases, with pustular psoriasis (PP) being the prototype, are immune-mediated diseases characterized by the presence of multiple pustules, resulting from neutrophil accumulation in the layer of epidermis. Sterile skin pustular eruption, like PP, is also observed in 20-30% of patients with adult-onset immunodeficiency syndrome (AOID) and anti-interferon γ autoantibodies (IFN-γ), leading to challenges in classification and diagnosis. While the mechanism underlying this similar phenotype remains unknown, genetic factors in relation to the immune system are suspected of playing an important role. Here, the association between human leukocyte antigen (HLA) genes, which play essential roles in antigen presentation, contributing to immune response, and the presence of skin pustules in AOID and PP was revealed. HLA genotyping of 41 patients from multiple centers in Thailand who presented with multiple sterile skin pustules (17 AOID patients and 24 PP patients) was conducted using a next-generation-sequencing-based approach. In comparison to healthy controls, HLA-B*13:01 (OR = 3.825, 95%CI: 2.08-7.035), C*03:04 (OR = 3.665, 95%CI: 2.102-6.39), and DQB1*05:02 (OR = 2.134, 95%CI: 1.326-3.434) were significantly associated with the group of aforementioned conditions having sterile cutaneous pustules, suggesting a common genetic-related mechanism. We found that DPB1*05:01 (OR = 3.851, p = 0.008) and DRB1*15:02 (OR = 3.195, p = 0.033) have a significant association with pustular reaction in AOID patients, with PP patients used as a control. A variant in the DRB1 gene, rs17885482 (OR = 9.073, p = 0.005), was observed to be a risk factor for PP when using AOID patients who had pustular reactions as a control group. DPB1*05:01 and DRB1*15:02 alleles, as well as the rs17885482 variant in the DRB1 gene, were proposed as novel biomarkers to differentiate PP and AOID patients who first present with multiple sterile skin pustules without known documented underlying conditions.

18th International HLA and Immunogenetics Workshop: Report on the SNP-HLA Reference Consortium (SHLARC) component

The SNP-HLA Reference Consortium (SHLARC), a component of the 18th International HLA and Immunogenetics Workshop, is aimed at collecting diverse and extensive human leukocyte antigen (HLA) data to create custom reference panels and enhance HLA imputation techniques. Genome-wide association studies (GWAS) have significantly contributed to identifying genetic associations with various diseases. The HLA genomic region has emerged as the top locus in GWAS, particularly in immune-related disorders. However, the limited information provided by single nucleotide polymorphisms (SNPs), the hallmark of GWAS, poses challenges, especially in the HLA region, where strong linkage disequilibrium (LD) spans several megabases. HLA imputation techniques have been developed using statistical inference in response to these challenges. These techniques enable the prediction of HLA alleles from genotyped GWAS SNPs. Here we present the SHLARC activities, a collaborative effort to create extensive, and multi-ethnic reference panels to enhance HLA imputation accuracy.

"HLA-C: evolution, epigenetics, and pathological implications in the major histocompatibility complex"

HLA-C, a gene located within the major histocompatibility complex, has emerged as a prominent target in biomedical research due to its involvement in various diseases, including cancer and autoimmune disorders; even though its recent addition to the MHC, the interaction between HLA-C and KIR is crucial for immune responses, particularly in viral infections. This review provides an overview of the structure, origin, function, and pathological implications of HLA-C in the major histocompatibility complex. In the last decade, we systematically reviewed original publications from Pubmed, ScienceDirect, Scopus, and Google Scholar. Our findings reveal that genetic variations in HLA-C can determine susceptibility or resistance to certain diseases. However, the first four exons of HLA-C are particularly susceptible to epigenetic modifications, which can lead to gene silencing and alterations in immune function. These alterations can manifest in diseases such as alopecia areata and psoriasis and can also impact susceptibility to cancer and the effectiveness of cancer treatments. By comprehending the intricate interplay between genetic and epigenetic factors that regulate HLA-C expression, researchers may develop novel strategies for preventing and treating diseases associated with HLA-C dysregulation.

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.

Whole-genome sequencing of 1,171 elderly admixed individuals from São Paulo, Brazil

As whole-genome sequencing (WGS) becomes the gold standard tool for studying population genomics and medical applications, data on diverse non-European and admixed individuals are still scarce. Here, we present a high-coverage WGS dataset of 1,171 highly admixed elderly Brazilians from a census-based cohort, providing over 76 million variants, of which ~2 million are absent from large public databases. WGS enables identification of ~2,000 previously undescribed mobile element insertions without previous description, nearly 5 Mb of genomic segments absent from the human genome reference, and over 140 alleles from HLA genes absent from public resources. We reclassify and curate pathogenicity assertions for nearly four hundred variants in genes associated with dominantly-inherited Mendelian disorders and calculate the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observe that whole-genome and HLA imputation could be significantly improved compared to available datasets since rare variation represents the largest proportion of input from WGS. These results demonstrate that even smaller sample sizes of underrepresented populations bring relevant data for genomic studies, especially when exploring analyses allowed only by WGS.

Whole genome sequencing (WGS) data on non-European and admixed individuals remains scarce. Here, the authors analyse WGS data from 1,171 admixed elderly Brazilians from a census cohort, characterising population-specific genetic variation and exploring the clinical utility of this expanded dataset.

Like Wings of a Bird: Functional Divergence and Complementarity between HLA-A and HLA-B Molecules

Human leukocyte antigen (HLA) genes are among the most polymorphic of our genome, as a likely consequence of balancing selection related to their central role in adaptive immunity. HLA-A and HLA-B genes were recently suggested to evolve through a model of joint divergent asymmetric selection conferring all human populations, including those with severe loss of diversity, an equivalent immune potential. However, the mechanisms by which these two genes might undergo joint evolution while displaying very distinct allelic profiles in populations are still unknown. To address this issue, we carried out extensive data analyses (among which factorial correspondence analysis and linear modeling) on 2,909 common and rare HLA-A, HLA-B, and HLA-C alleles and 200,000 simulated pathogenic peptides by taking into account sequence variation, predicted peptide-binding affinity and HLA allele frequencies in 123 populations worldwide. Our results show that HLA-A and HLA-B (but not HLA-C) molecules maintain considerable functional divergence in almost all populations, which likely plays an instrumental role in their immune defense. We also provide robust evidence of functional complementarity between HLA-A and HLA-B molecules, which display asymmetric relationships in terms of amino acid diversity at both inter- and intraprotein levels and in terms of promiscuous or fastidious peptide-binding specificities. Like two wings of a flying bird, the functional complementarity of HLA-A and HLA-B is a perfect example, in our genome, of duplicated genes sharing their capacity of assuming common vital functions while being submitted to complex and sometimes distinct environmental pressures.

KIR2DL4 genetic diversity in a Brazilian population sample: implications for transcription regulation and protein diversity in samples with different ancestry backgrounds

KIR2DL4 is an important immune modulator expressed in natural killer cells; HLA-G is its main ligand. We have characterized the KIR2DL4 genetic diversity by considering the promoter, all exons, and all introns in a highly admixed Brazilian population sample and by using massively parallel sequencing. We introduce a molecular method to amplify and to sequence the complete KIR2DL4 gene. To avoid the mapping bias and genotype errors commonly observed in gene families, we have developed and validated a bioinformatic pipeline designed to minimize these errors and applied it to survey the variability of 220 individuals from the State of São Paulo, southeastern Brazil. We have also compared the KIR2DL4 genetic diversity in the Brazilian cohort with the diversity previously reported by the 1000Genomes consortium. KIR2DL4 presents high linkage disequilibrium throughout the gene, with coding sequences associated with specific promoters. There are few but divergent promoter haplotypes. We have also detected many new KIR2DL4 sequences, all bearing nucleotide exchanges in introns and encoding previously described proteins. Exons 3 and 4, which encode the external domains, are the most variable. The ancestry background influences the KIR2DL4 allele frequencies and must be considered for association studies regarding KIR2DL4.