High-resolution characterization of allelic and haplotypic HLA frequency distribution in a Spanish population using high-throughput next-generation sequencing

Next-generation sequencing (NGS) at the HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1 and -DRB3/4/5 loci was performed on 282 healthy unrelated individuals from different major regions of Spain. High-resolution HLA genotypes defined by full sequencing of class I loci and extended coverage of class II loci were obtained to determine allele frequencies and also to estimate extended haplotype frequencies. HLA alleles were typed at the highest resolution level (4-field level, 4FL); with exception of a minor deviation in HLA-DPA1, no statistically significant deviations from expected Hardy Weinberg Equilibrium (HWE) proportions were observed for all other HLA loci. This study provides new 4FL-allele and -haplotype frequencies estimated for the first time in the Spanish population. Furthermore, our results describe extended haplotypes (including the less frequently typed HLA-DPA1 and HLA-DQA1 loci) and show distinctive haplotype associations found at 4FL-allele definition in this Spanish population study. The distinctive allelic and haplotypic diversity found at the 4FL reveals the high level of heterozygosity and specific haplotypic associations displayed that were not apparent at 2-field level (2FL). Overall, these results may contribute as a useful reference source for future population studies, for HLA-disease association studies as a healthy control group dataset and for improving donor recruitment strategies of bone marrow registries. HLA genotyping data of this Spanish population cohort was also included in the 17th International Histocompatibility and Immunogenetics Workshop (IHIW) as part of the study of HLA diversity in unrelated worldwide populations using NGS.     

A reverse-engineering strategy utilizing the integration of single antigen beads and NGS HLA genotypes to detect potential antibody inducing epitopes

Central to the idea of antibody recognition is some degree of foreignness of the target antigen compared to the antibody producer. Epitopes are distinct regions on an antigen to which antibody can be elicited and bound. However, for HLA antigens, there is no consensus definition of what represents the minimal functional immunogenic unit of dissimilarity. To assess this in an unbiased way, we developed a reverse engineering software strategy based on donor specific antibodies defined by single antigen beads and full length genomic high resolution HLA typing by NGS of recipients and donors (332 transplant pairs). Starting with the ATG of Exon 1 and moving stepwise one amino acid at a time for each of the following triplets, the algorithm compared every possible amino acid triplet of the recipient and donor for 11 loci (A, B, C, DRB1, DRB3, DRB4, DRB5, DQA1, DQB1, DPA1, DPB1). Results were agnostic with respect to HLA class, not restricted to just the mature protein, and not influenced by existing maps (e.g., IMGT, or epitope models). We also developed web-based functions in the 17th IHIWS database to collect the unbiased triplets so that we could group the transplant pairs with the same donor specific antibodies and find shared triplets within the groups as potential core or essential epitopes that trigger the antibody formation. Profiling the pairs where the same DSA was identified led to identification of discrete amino acid triplets shared among the pairs irrespective of HLA match. The potential epitopes were mapped onto the 3D protein structure for reference.     

HLA-F coding and regulatory segments variability determined by massively parallel sequencing procedures in a Brazilian population sample

Human Leucocyte Antigen F (HLA-F) is a non-classical HLA class I gene distinguished from its classical counterparts by low allelic polymorphism and distinctive expression patterns. Its exact function remains unknown. It is believed that HLA-F has tolerogenic and immune modulatory properties. Currently, there is little information regarding the HLA-F allelic variation among human populations and the available studies have evaluated only a fraction of the HLA-F gene segment and/or have searched for known alleles only. Here we present a strategy to evaluate the complete HLA-F variability including its 5′ upstream, coding and 3′ downstream segments by using massively parallel sequencing procedures. HLA-F variability was surveyed on 196 individuals from the Brazilian Southeast. The results indicate that the HLA-F gene is indeed conserved at the protein level, where thirty coding haplotypes or coding alleles were detected, encoding only four different HLA-F full-length protein molecules. Moreover, a same protein molecule is encoded by 82.45% of all coding alleles detected in this Brazilian population sample. However, the HLA-F nucleotide and haplotype variability is much higher than our current knowledge both in Brazilians and considering the 1000 Genomes Project data. This protein conservation is probably a consequence of the key role of HLA-F in the immune system physiology.