Hla-mapper: An application to optimize the mapping of HLA sequences produced by massively parallel sequencing procedures

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.

Optimal population-specific HLA imputation with dimension reduction

Human genomics has quickly evolved, powering genome-wide association studies (GWASs). SNP-based GWASs cannot capture the intense polymorphism of HLA genes, highly associated with disease susceptibility. There are methods to statistically impute HLA genotypes from SNP-genotypes data, but lack of diversity in reference panels hinders their performance. We evaluated the accuracy of the 1000 Genomes data as a reference panel for imputing HLA from admixed individuals of African and European ancestries, focusing on (a) the full dataset, (b) 10 replications from 6 populations, and (c) 19 conditions for the custom reference panels. The full dataset outperformed smaller models, with a good F1-score of 0.66 for HLA-B. However, custom models outperformed the multiethnic or population models of similar size (F1-scores up to 0.53, against up to 0.42). We demonstrated the importance of using genetically specific models for imputing populations, which are currently underrepresented in public datasets, opening the door to HLA imputation for every genetic population.

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-F transcriptional and protein differential expression according to its genetic polymorphisms

Many specificities single out HLA-F: its structure, expression regulation at cell membrane and function. HLA-F mRNA is detected in the most cell types and the protein is localized in the ER and Golgi apparatus. When expressed at cell surface, HLA-F may be associated to β2-microglobulin and peptide or expressed as an open-conformer molecule. HLA-F reaches the membrane upon activation of different primary cell types and cell-lines. HLA-F has its highest affinity for the KIR3DS1-activating NK receptor, but also binds inhibitory immune receptors. Some studies reported that HLA-F expression is associated with its genotype. Higher HLA-F mRNA expression associated with F*01:01:02, and 3 noncoding SNPs, rs1362126, rs2523405, and rs2523393, located in HLA-F-AS1 or upstream the HLA-F sequence were associated with HLA-F mRNA expression. Given the implication of HLA-F in many clinical setting, and the undisclosed process of its expression regulation, we aim to confirm the effect of the aforementioned SNPs with HLA-F transcriptional and protein expression. We analyzed the distribution, frequency and linkage disequilibrium of these SNPs at worldwide scale in the 1000 Genomes Project samples. Influence on the genotype of each SNP on HLA-F expression was explored using RNAseq data from the 1000 Genomes Project, and using Q-PCR and intracellular cytometry in PBMC from healthy individuals. Our results show that the SNPs under studied displayed remarkably different allelic proportion according to geography and confirm that rs1362126, rs2523405, and rs2523393 displayed the most concordant results, with the highest effect size and a double-dose effect.

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.

The + 3010/C single nucleotide polymorphism (rs1710) at the HLA-G 3' untranslated region is associated with a short transcript exhibiting a deletion of 92 nucleotides

The physiological expression of HLA-G is mainly observed in the placenta, playing an essential role in maternal-fetal tolerance. Among the HLA-G mRNA alternative transcripts, the one lacking 92 bases at the HLA-G 3' untranslated region (3'UTR), the 92bDel transcript, is more stable, is associated with increased HLA-G soluble levels, and was observed in individuals presenting a 14 bp insertion (14 bp⁺) at the 3'UTR. We investigated the presence of the 92bDel transcript in placenta samples, correlating its expression levels with the HLA-G polymorphisms at the 3'UTR. The 14 bp⁺ allele correlates with the presence of the 92bDel transcript. However, the polymorphism triggering this alternative splicing is the + 3010/C allele (rs1710, allele C). Most 14 bp⁺ haplotypes (UTR-2/-5/-7) present allele + 3010/C. However, 14 bp^- haplotypes such as UTR-3 are also associated with + 3010/C, and the 92bDel transcript can be detected in homozygous samples for the 14 bp- allele carrying at least one copy of UTR-3. The UTR-3 haplotype is associated with alleles G*01:04 and the HLA-G lineage HG0104, which is a high-expressing lineage. The only HLA-G lineage that is not likely to produce this transcript is HG010101, associated with the + 3010/G allele. This functional difference may be advantageous, considering the high worldwide frequency of the HG010101 lineage. Therefore, HLA-G lineages are functionally distinct regarding the 92bDel transcript expression, and the 3010/C allele triggers the alternative splicing that produces this shorter and more stable transcript.

Evidence for Epistatic Interaction between HLA-G and LILRB1 in the Pathogenesis of Nonsegmental Vitiligo

Vitiligo is the most frequent cause of depigmentation worldwide. Genetic association studies have discovered about 50 loci associated with disease, many with immunological functions. Among them is HLA-G, which modulates immunity by interacting with specific inhibitory receptors, mainly LILRB1 and LILRB2. Here we investigated the LILRB1 and LILRB2 association with vitiligo risk and evaluated the possible role of interactions between HLA-G and its receptors in this pathogenesis. We tested the association of the polymorphisms of HLA-G, LILRB1, and LILRB2 with vitiligo using logistic regression along with adjustment by ancestry. Further, methods based on the multifactor dimensionality reduction (MDR) approach (MDR v.3.0.2, GMDR v.0.9, and MB-MDR) were used to detect potential epistatic interactions between polymorphisms from the three genes. An interaction involving rs9380142 and rs2114511 polymorphisms was identified by all methods used. The polymorphism rs9380142 is an HLA-G 3'UTR variant (+3187) with a well-established role in mRNA stability. The polymorphism rs2114511 is located in the exonic region of LILRB1. Although no association involving this SNP has been reported, ChIP-Seq experiments have identified this position as an EBF1 binding site. These results highlight the role of an epistatic interaction between HLA-G and LILRB1 in vitiligo pathogenesis.

Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.

The oldest unvaccinated Covid-19 survivors in South America

BACKGROUND

Although older adults are at a high risk of severe or critical Covid-19, there are many cases of unvaccinated centenarians who had a silent infection or recovered from mild or moderate Covid-19. We studied three Brazilian supercentenarians, older than 110 years, who survived Covid-19 in 2020 before being vaccinated.

RESULTS

Despite their advanced age, humoral immune response analysis showed that these individuals displayed robust levels of IgG and neutralizing antibodies (NAbs) against SARS-CoV-2. Enrichment of plasma proteins and metabolites related to innate immune response and host defense was also observed. None presented autoantibodies (auto-Abs) to type I interferon (IFN). Furthermore, these supercentenarians do not carry rare variants in genes underlying the known inborn errors of immunity, including particular inborn errors of type I IFN.

CONCLUSION

These observations suggest that their Covid-19 resilience might be a combination of their genetic background and their innate and adaptive immunity.

MUC22, HLA-A, and HLA-DOB variants and COVID-19 in resilient super-agers from Brazil

BACKGROUND

Although aging correlates with a worse prognosis for Covid-19, super elderly still unvaccinated individuals presenting mild or no symptoms have been reported worldwide. Most of the reported genetic variants responsible for increased disease susceptibility are associated with immune response, involving type I IFN immunity and modulation; HLA cluster genes; inflammasome activation; genes of interleukins; and chemokines receptors. On the other hand, little is known about the resistance mechanisms against SARS-CoV-2 infection. Here, we addressed polymorphisms in the MHC region associated with Covid-19 outcome in super elderly resilient patients as compared to younger patients with a severe outcome.

METHODS

SARS-CoV-2 infection was confirmed by RT-PCR test. Aiming to identify candidate genes associated with host resistance, we investigated 87 individuals older than 90 years who recovered from Covid-19 with mild symptoms or who remained asymptomatic following positive test for SARS-CoV-2 as compared to 55 individuals younger than 60 years who had a severe disease or died due to Covid-19, as well as to the general elderly population from the same city. Whole-exome sequencing and an in-depth analysis of the MHC region was performed. All samples were collected in early 2020 and before the local vaccination programs started.

RESULTS

We found that the resilient super elderly group displayed a higher frequency of some missense variants in the MUC22 gene (a member of the mucins' family) as one of the strongest signals in the MHC region as compared to the severe Covid-19 group and the general elderly control population. For example, the missense variant rs62399430 at MUC22 is two times more frequent among the resilient super elderly (p = 0.00002, OR = 2.24).

CONCLUSION

Since the pro-inflammatory basal state in the elderly may enhance the susceptibility to severe Covid-19, we hypothesized that MUC22 might play an important protective role against severe Covid-19, by reducing overactive immune responses in the senior population.

Follow-up of young adult monozygotic twins after simultaneous critical coronavirus disease 2019: A case report

Background

The influence of the host genome on coronavirus disease 2019 (COVID-19) susceptibility and severity is supported by reports on monozygotic (MZ) twins where both were infected simultaneously with similar disease outcomes, including several who died due to the SARS-CoV-2 infection within days apart. However, successive exposures to pathogens throughout life along with other environmental factors make the immune response unique for each individual, even among MZ twins.

Case presentation and methods

Here we report a case of a young adult monozygotic twin pair, who caught attention since both presented simultaneously severe COVID-19 with the need for oxygen support despite age and good health conditions. One of the twins, who spent more time hospitalized, reported symptoms of long-COVID even 7 months after infection. Immune cell profile and specific responses to SARS-CoV-2 were evaluated as well as whole exome sequencing.

Conclusion

Although the MZ twin brothers shared the same genetic mutations which may be associated with their increased risk of developing severe COVID-19, their clinical progression was different, reinforcing the role of both immune response and genetics in the COVID-19 presentation and course. Besides, post-COVID syndrome was observed in one of them, corroborating an association between the duration of hospitalization and the occurrence of long-COVID symptoms.

How HLA diversity is apportioned: influence of selection and relevance to transplantation

In his 1972 paper ‘The apportionment of human diversity’, Lewontin showed that, when averaged over loci, genetic diversity is predominantly attributable to differences among individuals within populations. However, selection can alter the apportionment of diversity of specific genes or genomic regions. We examine genetic diversity at the human leucocyte antigen (HLA) loci, located within the major histocompatibility complex (MHC) region. HLA genes code for proteins that are critical to adaptive immunity and are well-documented targets of balancing selection. The single-nucleotide polymorphisms (SNPs) within HLA genes show strong signatures of balancing selection on large timescales and are broadly shared among populations, displaying low F_ST values. However, when we analyse haplotypes defined by these SNPs (which define ‘HLA alleles’), we find marked differences in frequencies between geographic regions. These differences are not reflected in the F_ST values because of the extreme polymorphism at HLA loci, illustrating challenges in interpreting F_ST. Differences in the frequency of HLA alleles among geographic regions are relevant to bone-marrow transplantation, which requires genetic identity at HLA loci between patient and donor. We discuss the case of Brazil's bone marrow registry, where a deficit of enrolled volunteers with African ancestry reduces the chance of finding donors for individuals with an MHC region of African ancestry.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

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.

Recurrence of COVID-19 associated with reduced T-cell responses in a monozygotic twin pair

Recurrence of COVID-19 in recovered patients has been increasingly reported. However, the immune mechanisms behind the recurrence have not been thoroughly investigated. The presence of neutralizing antibodies (nAbs) in recurrence/reinfection cases suggests that other types of immune response are involved in protection against recurrence. Here, we investigated the innate type I/III interferon (IFN) response, binding and nAb assays and T-cell responses to severe acute respiratory distress syndrome coronavirus 2 (SARS-CoV-2) with IFN gamma (IFNγ) enzyme-linked spot assay (ELISPOT) in three pairs of young adult monozygotic (MZ) twins with previous confirmed COVID-19, one of them presenting a severe recurrence four months after the initial infection. Twin studies have been of paramount importance to comprehend the immunogenetics of infectious diseases. Each MZ twin pair was previously exposed to SARS-CoV-2, as seen by clinical reports. The six individuals presented similar overall recovered immune responses except for the recurrence case, who presented a drastically reduced number of recognized SARS-CoV-2 T-cell epitopes on ELISPOT as compared to her twin sister and the other twin pairs. Our results suggest that the lack of a broad T-cell response to initial infection may have led to recurrence, emphasizing that an effective SARS-CoV-2-specific T-cell immune response is key for complete viral control and avoidance of clinical recurrence of COVID-19.

HLA associations, somatic loss of HLA expression, and clinical outcomes in immune aplastic anemia

Immune aplastic anemia (AA) features somatic loss of HLA class I allele expression on bone marrow cells, consistent with a mechanism of escape from T-cell-mediated destruction of hematopoietic stem and progenitor cells. The clinical significance of HLA abnormalities has not been well characterized. We examined the somatic loss of HLA class I alleles and correlated HLA loss and mutation-associated HLA genotypes with clinical presentation and outcomes after immunosuppressive therapy in 544 AA patients. HLA class I allele loss was detected in 92 (22%) of the 412 patients tested, in whom there were 393 somatic HLA gene mutations and 40 instances of loss of heterozygosity. Most frequently affected was HLA-B*14:02, followed by HLA-A*02:01, HLA-B*40:02, HLA-B*08:01, and HLA-B*07:02. HLA-B*14:02, HLA-B*40:02, and HLA-B*07:02 were also overrepresented in AA. High-risk clonal evolution was correlated with HLA loss, HLA-B*14:02 genotype, and older age, which yielded a valid prediction model. In 2 patients, we traced monosomy 7 clonal evolution from preexisting clones harboring somatic mutations in HLA-A*02:01 and HLA-B*40:02. Loss of HLA-B*40:02 correlated with higher blood counts. HLA-B*07:02 and HLA-B*40:01 genotypes and their loss correlated with late-onset of AA. Our results suggest the presence of specific immune mechanisms of molecular pathogenesis with clinical implications. HLA genotyping and screening for HLA loss may be of value in the management of immune AA. This study was registered at clinicaltrials.gov as NCT00001964, NCT00061360, NCT00195624, NCT00260689, NCT00944749, NCT01193283, and NCT01623167.

Approaching Genetics Through the MHC Lens: Tools and Methods for HLA Research

The current SARS-CoV-2 pandemic era launched an immediate and broad response of the research community with studies both about the virus and host genetics. Research in genetics investigated HLA association with COVID-19 based on in silico, population, and individual data. However, they were conducted with variable scale and success; convincing results were mostly obtained with broader whole-genome association studies. Here, we propose a technical review of HLA analysis, including basic HLA knowledge as well as available tools and advice. We notably describe recent algorithms to infer and call HLA genotypes from GWAS SNPs and NGS data, respectively, which opens the possibility to investigate HLA from large datasets without a specific initial focus on this region. We thus hope this overview will empower geneticists who were unfamiliar with HLA to run MHC-focused analyses following the footsteps of the Covid-19|HLA & Immunogenetics Consortium.

HLA-G genetic diversity and evolutive aspects in worldwide populations

HLA-G is a promiscuous immune checkpoint molecule. The HLA-G gene presents substantial nucleotide variability in its regulatory regions. However, it encodes a limited number of proteins compared to classical HLA class I genes. We characterized the HLA-G genetic variability in 4640 individuals from 88 different population samples across the globe by using a state-of-the-art method to characterize polymorphisms and haplotypes from high-coverage next-generation sequencing data. We also provide insights regarding the HLA-G genetic diversity and a resource for future studies evaluating HLA-G polymorphisms in different populations and association studies. Despite the great haplotype variability, we demonstrated that: (1) most of the HLA-G polymorphisms are in introns and regulatory sequences, and these are the sites with evidence of balancing selection, (2) linkage disequilibrium is high throughout the gene, extending up to HLA-A, (3) there are few proteins frequently observed in worldwide populations, with lack of variation in residues associated with major HLA-G biological properties (dimer formation, interaction with leukocyte receptors). These observations corroborate the role of HLA-G as an immune checkpoint molecule rather than as an antigen-presenting molecule. Understanding HLA-G variability across populations is relevant for disease association and functional studies.

MHC Variants Associated With Symptomatic Versus Asymptomatic SARS-CoV-2 Infection in Highly Exposed Individuals

Despite the high number of individuals infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) who develop coronavirus disease 2019 (COVID-19) symptoms worldwide, many exposed individuals remain asymptomatic and/or uninfected and seronegative. This could be explained by a combination of environmental (exposure), immunological (previous infection), epigenetic, and genetic factors. Aiming to identify genetic factors involved in immune response in symptomatic COVID-19 as compared to asymptomatic exposed individuals, we analyzed 83 Brazilian couples where one individual was infected and symptomatic while the partner remained asymptomatic and serum-negative for at least 6 months despite sharing the same bedroom during the infection. We refer to these as "discordant couples". We performed whole-exome sequencing followed by a state-of-the-art method to call genotypes and haplotypes across the highly polymorphic major histocompatibility complex (MHC) region. The discordant partners had comparable ages and genetic ancestry, but women were overrepresented (65%) in the asymptomatic group. In the antigen-presentation pathway, we observed an association between HLA-DRB1 alleles encoding Lys at residue 71 (mostly DRB1*03:01 and DRB1*04:01) and DOB*01:02 with symptomatic infections and HLA-A alleles encoding 144Q/151R with asymptomatic seronegative women. Among the genes related to immune modulation, we detected variants in MICA and MICB associated with symptomatic infections. These variants are related to higher expression of soluble MICA and low expression of MICB. Thus, quantitative differences in these molecules that modulate natural killer (NK) activity could contribute to susceptibility to COVID-19 by downregulating NK cell cytotoxic activity in infected individuals but not in the asymptomatic partners.

An immunogenetic view of COVID-19

Meeting the challenges brought by the COVID-19 pandemic requires an interdisciplinary approach. In this context, integrating knowledge of immune function with an understanding of how genetic variation influences the nature of immunity is a key challenge. Immunogenetics can help explain the heterogeneity of susceptibility and protection to the viral infection and disease progression. Here, we review the knowledge developed so far, discussing fundamental genes for triggering the innate and adaptive immune responses associated with a viral infection, especially with the SARS-CoV-2 mechanisms. We emphasize the role of the HLA and KIR genes, discussing what has been uncovered about their role in COVID-19 and addressing methodological challenges of studying these genes. Finally, we comment on questions that arise when studying admixed populations, highlighting the case of Brazil. We argue that the interplay between immunology and an understanding of genetic associations can provide an important contribution to our knowledge of COVID-19.

Human leukocyte antigen (HLA)-F and -G gene polymorphisms and haplotypes are associated with malaria susceptibility in the Beninese Toffin children

BACKGROUND

Little attention has been devoted to the role of the immunoregulatory HLA-E/-F/-G genes in malaria. We evaluated the entire HLA-E/-F/-G variability in Beninese children highly exposed to Plasmodium falciparum (P.f.) malaria.

METHODS

154 unrelated children were followed-up for six months and evaluated for the presence and number of malaria episodes. HLA-E/-F/-G genes were genotyped using massively parallel sequencing. Anti P.f. antibodies were evaluated using ELISA.

RESULTS

Children carrying the G allele at HLA-F (-1499,rs183540921) showed increased P.f. asymptomatic/symptomatic ratio, suggesting that these children experienced more asymptomatic P.f. episodes than symptomatic one. Children carrying HLA-G-UTR-03 haplotype exhibited increased risk for symptomatic P.f. episodes and showed lower IgG2 response against P.f. GLURP-R2 when compared to the non-carriers. No associations were observed for the HLA-E gene.

CONCLUSION

HLA-F associations may be related to the differential expression profiles of the encoded immunomodulatory molecules, and the regulatory sites at the HLA-G 3'UTR may be associated to posttranscriptional regulation of HLA-G and to host humoral response against P.f.

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.

Manifesting carriers of X-linked myotubular myopathy: Genetic modifiers modulating the phenotype

Objective

To analyze the modulation of the phenotype in manifesting carriers of recessive X-linked myotubular myopathy (XLMTM), searching for possible genetic modifiers.

Methods

Twelve Brazilian families with XLMTM were molecularly and clinically evaluated. In 2 families, 4 of 6 and 2 of 5 manifesting female carriers were identified. These females were studied for X chromosome inactivation. In addition, whole-exome sequencing was performed, looking for possible modifier variants. We also determined the penetrance rate among carriers of the mutations responsible for the condition.

Results

Mutations in the MTM1 gene were identified in all index patients from the 12 families, being 4 of them novel. In the heterozygotes, X chromosome inactivation was random in 3 of 4 informative manifesting carriers. The disease penetrance rate was estimated to be 30%, compatible with incomplete penetrance. Exome comparative analyses identified variants within a segment of 4.2 Mb on chromosome 19, containing the killer cell immunoglobulin-like receptor cluster of genes that were present in all nonmanifesting carriers and absent in all manifesting carriers. We hypothesized that these killer cell immunoglobulin-like receptor variants may modulate the phenotype, acting as a protective factor in the nonmanifesting carriers.

Conclusions

Affected XLMTM female carriers have been described with a surprisingly high frequency for a recessive X-linked disease, raising the question about the pattern of inheritance or the role of modifier factors acting on the disease phenotype. We demonstrated the possible existence of genetic mechanisms and variants accountable for the clinical manifestation in these women, which can become future targets for therapies.

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

Human leukocyte antigen-C (HLA-C) is a classical HLA class I molecule that binds and presents peptides to cytotoxic T lymphocytes in the cell surface. HLA-C has a dual function because it also interacts with Killer-cell immunoglobulin-like receptors (KIR) receptors expressed in natural killer and T cells, modulating their activity. The structure and diversity of the HLA-C regulatory regions, as well as the relationship among variants along the HLA-C locus, are poorly addressed, and few population-based studies explored the HLA-C variability in the entire gene in different population samples. Here we present a molecular and bioinformatics method to evaluate the entire HLA-C diversity, including regulatory sequences. Then, we applied this method to survey the HLA-C diversity in two population samples with different demographic histories, one highly admixed from Brazil with major European contribution, and one from Benin with major African contribution. The HLA-C promoter and 3'UTR were very polymorphic with the presence of few, but highly divergent haplotypes. These segments also present conserved sequences that are shared among different primate species. Nucleotide diversity was higher in other segments rather than exons 2 and 3, particularly around exon 5 and the second half of the 3'UTR region. We detected evidence of balancing selection on the entire HLA-C locus and positive selection in the HLA-C leader peptide, for both populations. HLA-C motifs previously associated with KIR interaction and expression regulation are similar between both populations. Each allele group is associated with specific regulatory sequences, reflecting the high linkage disequilibrium along the entire HLA-C locus in both populations.

High-Throughput Sequencing to Investigate Associations Between HLA Genes and Metamizole-Induced Agranulocytosis

Background and Objective: Agranulocytosis is a rare and potentially life-threatening complication of metamizole (dipyrone) intake that is characterized by a loss of circulating neutrophil granulocytes. While the mechanism underlying this adverse drug reaction is not well understood, involvement of the immune system has been suggested. In addition, associations between genetic variants in the Human Leukocyte Antigen (HLA) region and agranulocytosis induced by other drugs have been reported. The aim of the present study was to assess whether genetic variants in classical HLA genes are associated with the susceptibility to metamizole-induced agranulocytosis (MIA) in a European population by targeted resequencing of eight HLA genes.

Design: A case-control cohort of Swiss patients with a history of neutropenia or agranulocytosis associated with metamizole exposure (n = 53), metamizole-tolerant (n = 39) and unexposed controls (n = 161) was recruited for this study. A high-throughput resequencing (HTS) and high-resolution typing method was used to sequence and analyze eight HLA loci in a discovery subset of this cohort (n = 31 cases, n = 38 controls). Identified candidate alleles were investigated in the full Swiss cohort as well as in two independent cohorts from Germany and Spain using HLA imputation from genome-wide SNP array data. In addition, variant calling based on HTS data was performed in the discovery subset for the class I genes HLA-A, -B, and -C using the HLA-specific mapper hla-mapper.

Results: Eight candidate alleles (p < 0.05) were identified in the discovery subset, of which HLA-C^∗04:01 was associated with MIA in the full Swiss cohort (p < 0.01) restricted to agranulocytosis (ANC < 0.5 × 10⁹/L) cases. However, no candidate allele showed a consistent association in the Swiss, German and Spanish cohorts. Analysis of individual sequence variants in class I genes produced consistent results with HLA typing but did not reveal additional small nucleotide variants associated with MIA.

Conclusion: Our results do not support an HLA-restricted T cell-mediated immune mechanism for MIA. However, we established an efficient high-resolution (three-field) eight-locus HTS HLA resequencing method to interrogate the HLA region and demonstrated the feasibility of its application to pharmacogenetic studies.

Contribution of unfixed transposable element insertions to human regulatory variation

Thousands of unfixed transposable element (TE) insertions segregate in the human population, but little is known about their impact on genome function. Recently, a few studies associated unfixed TE insertions to mRNA levels of adjacent genes, but the biological significance of these associations, their replicability across cell types and the mechanisms by which they may regulate genes remain largely unknown. Here, we performed a TE-expression QTL analysis of 444 lymphoblastoid cell lines (LCL) and 289 induced pluripotent stem cells using a newly developed set of genotypes for 2743 polymorphic TE insertions. We identified 211 and 176 TE-eQTL acting in cis in each respective cell type. Approximately 18% were shared across cell types with strongly correlated effects. Furthermore, analysis of chromatin accessibility QTL in a subset of the LCL suggests that unfixed TEs often modulate the activity of enhancers and other distal regulatory DNA elements, which tend to lose accessibility when a TE inserts within them. We also document a case of an unfixed TE likely influencing gene expression at the post-transcriptional level. Our study points to broad and diverse cis-regulatory effects of unfixed TEs in the human population and underscores their plausible contribution to phenotypic variation. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

HLA-H: Transcriptional Activity and HLA-E Mobilization

Little attention is paid to pseudogenes from the highly polymorphic HLA genetic region. The pseudogene HLA-H is defined as a non-functional gene because it is deleted at different frequencies in humans and because it encodes a potentially non-functional truncated protein. However, different studies have shown HLA-H transcriptional activity. We formerly identified 13 novel HLA-H alleles, including the H*02:07 allele, which reaches 19.6% in East Asian populations and encodes a full-length HLA protein. The aims of this study were to explore the expression and possible function of the HLA-H molecule. HLA-H may act as a transmembrane molecule and/or indirectly via its signal peptide by mobilizing HLA-E to the cell surface. We analyzed HLA-H RNA expression in Peripheral Blood Mononuclear Cells (PBMC), Human Bronchial Epithelial Cells (HBEC), and available RNA sequencing data from lymphoblastoid cell lines, and we looked to see whether HLA-E was mobilized at the cell surface by the HLA-H signal peptide. Our data confirmed that HLA-H is transcribed at similar levels to HLA-G. We characterized a hemizygous effect in HLA-H expression, and expression differed according to HLA-H alleles; most interestingly, the HLA-H*02:07 allele had the highest level of mRNA expression. We showed that HLA-H signal peptide incubation mobilized HLA-E molecules at the cell surface of T-Lymphocytes, monocytes, B-Lymphocytes, and primary epithelial cells. Our results suggest that HLA-H may be functional but raises many biological issues that need to be addressed.

HLA-G Polymorphisms Are Associated with Non-segmental Vitiligo among Brazilians

(1) Background: Vitiligo is characterized by white patches on the skin caused by loss of melanocyte activity or the absence of these cells. The available treatments minimize the symptoms by retarding the process of skin depigmentation or re-pigmenting the affected regions. New studies are required for a better comprehension of the mechanisms that trigger the disease and for the development of more efficient treatments. Studies have suggested an autoimmune feature for vitiligo, based on the occurrence of other autoimmune diseases in vitiligo patients and their relatives, and on the involvement of genes related to the immune response. (2) Methods: We evaluated, by massive parallel sequencing, polymorphisms of the HLA-G gene in vitiligo patients and control samples, to verify if variants of this gene could influence the susceptibility to vitiligo. (3) Results: We detected an association with non-segmental vitiligo regarding the haplotype Distal-010101a/G*01:01:01:01/UTR-1, adjusting for population stratification by using ancestry-informative markers (AIMs). (4) Conclusions: It remains unclear whether the HLA-G variants associated with vitiligo were detected because of the high linkage disequilibrium (LD) with HLA-A*02, or if the HLA-A variants previously reported as associated with vitiligo were detected because of the high LD with HLA-G*01:01:01:01/UTR-1, or if both genes jointly contribute to vitiligo susceptibility.

HLA-A promoter, coding, and 3'UTR sequences in a Brazilian cohort, and their evolutionary aspects

HLA-A is the second most polymorphic locus of the human leucocyte antigen (HLA) complex encoding a key molecule for antigen presentation and NK cell modulation. Many studies have evaluated HLA-A variability in worldwide populations, focusing mainly on exons, but the regulatory segments have been poorly characterized. HLA-A variability is particularly high in the segment encoding the peptide-binding groove (exons 2 and 3), which is related to the antigen presentation function and the balancing selection in these segments. Here we evaluate the genetic diversity of the HLA-A gene considering a continuous segment encompassing the extended promoter (1.5 kb upstream of the first translated ATG), all exons and introns, and the entire 3' untranslated region, by using massively parallel sequencing. To achieve this goal, we used a freely available bioinformatics workflow that optimizes read mapping for HLA genes and defines complete sequences using either the phase among variable sites directly observed in sequencing data and probabilistic models. The HLA-A variability detected in a highly admixed population sample from Brazil shows that the HLA-A regulatory segments present few, but divergent sequences. The regulatory segments are in close association with the coding alleles. Both exons and introns are highly variable. Moreover, patterns of molecular diversity suggest that the promoter, in addition to the coding region, might be under the same selective pressure, but a different scenario arises when it comes to exon 4 and the 3'UTR segment.