Impact of HLA class I functional divergence on HIV control.

Heterozygosity of Human leukocyte antigen (HLA) class I genes is linked to beneficial outcomes after HIV infection, presumably through greater breadth of HIV epitope presentation and cytotoxic T cell response. Distinct allotype pairs, however, differ in the extent to which they bind shared sets of peptides. We developed a functional divergence metric that measures pairwise complementarity of allotype-associated peptide binding profiles. Greater functional divergence for pairs of HLA-A and/or HLA-B allotypes was associated with slower AIDS progression and independently with enhanced viral load control. The metric predicts immune breadth at the peptide level rather than gene level and redefines HLA heterozygosity as a continuum differentially affecting disease outcome. Functional divergence may affect response to additional infections, vaccination, immunotherapy, and other diseases where HLA heterozygote advantage occurs.

The beneficial effects of commensal E. coli for colon epithelial cell recovery are related with Formyl peptide receptor 2 (Fpr2) in epithelial cells.

BACKGROUND: Formyl peptide receptor 2 (Fpr2) plays a crucial role in colon homeostasis and microbiota balance. Commensal E. coli is known to promote the regeneration of damaged colon epithelial cells. The aim of the study was to investigate the connection between E. coli and Fpr2 in the recovery of colon epithelial cells. RESULTS: The deficiency of Fpr2 was associated with impaired integrity of the colon mucosa and an imbalance of microbiota, characterized by the enrichment of Proteobacteria in the colon. Two serotypes of E. coli, O22:H8 and O91:H21, were identified in the mouse colon through complete genome sequencing. E. coli O22:H8 was found to be prevalent in the gut of mice and exhibited lower virulence compared to O91:H21. Germ-free (GF) mice that were pre-orally inoculated with E. coli O22:H8 showed reduced susceptibility to chemically induced colitis, increased proliferation of epithelial cells, and improved mouse survival. Following infection with E. coli O22:H8, the expression of Fpr2 in colon epithelial cells was upregulated, and the products derived from E. coli O22:H8 induced migration and proliferation of colon epithelial cells through Fpr2. Fpr2 deficiency increased susceptibility to chemically induced colitis, delayed the repair of damaged colon epithelial cells, and heightened inflammatory responses. Additionally, the population of E. coli was observed to increase in the colons of Fpr2-/- mice with colitis. CONCLUSION: Commensal E. coli O22:H8 stimulated the upregulation of Fpr2 expression in colon epithelial cells, and the products from E. coli induced migration and proliferation of colon epithelial cells through Fpr2. Fpr2 deficiency led to an increased E. coli population in the colon and delayed recovery of damaged colon epithelial cells in mice with colitis. Therefore, Fpr2 is essential for the effects of commensal E. coli on colon epithelial cell recovery.

Distinct frequency patterns of LILRB3 and LILRA6 allelic variants in Europeans.

The leukocyte immunoglobulin-like receptor (LILR)B3 and LILRA6 genes encode homologous myeloid inhibitory and activating orphan receptors, respectively. Both genes exhibit a strikingly high level of polymorphism at the amino acid level and LILRA6 (but not LILRB3) displays copy number variation (CNV). Although multiple alleles have been reported for both genes, limited data is available on frequencies of these alleles among humans. We have sequenced LILRB3/A6 exons encoding signal peptides and ectodomains in 91 healthy blood donors of European descent who carry one or two copies of LILRA6 per diploid genome. Analysis of haplotypes among individuals with two LILRA6 copies, representing the majority in this cohort (N = 86), shows that common LILRB3 and LILRA6 alleles encode some distinct amino acid sequences in homologous regions of the receptors, which could potentially impact their respective functions differentially. Comparison of sequences in individuals with one vs. two copies of LILRA6 supports non-allelic homologous recombination between LILRB3 and LILRA6 as a mechanism for generating LILRA6 CNV and LILRB3 diversity. These data characterize LILRB3/LILRA6 genetic variation in more detail than previously described and underscore the need to determine their ligands.

Population-specific diversity of the immunoglobulin constant heavy G chain (IGHG) genes.

Human immunoglobulin G (IgG) molecules, IgG1, IgG2 and IgG3, exhibit substantial inter-individual variation in their constant heavy chain regions, as discovered by serological methods. This polymorphism is encoded by the IGHG1, IGHG2, and IGHG3 genes and may influence antibody function. We sequenced the coding fragments of these genes in 95 European Americans, 94 African Americans, and 94 Black South Africans. Striking differences were observed between the population groups, including extremely low amino acid sequence variation in IGHG1 among South Africans, and higher IGHG2 and IGHG3 diversity in individuals of African descent compared to individuals of European descent. Molecular definition of the loci illustrates a greater level of allelic polymorphism than previously described, including the presence of common IGHG2 and IGHG3 variants that were indistinguishable serologically. Comparison of our data with the 1000 Genome Project sequences indicates overall agreement between the datasets, although some inaccuracies in the 1000 Genomes Project are likely. These data represent the most comprehensive analysis of IGHG polymorphisms across major populations, which can now be applied to deciphering their functional impact.

Distinct contributions of cathelin-related antimicrobial peptide (CRAMP) derived from epithelial cells and macrophages to colon mucosal homeostasis.

The cathelin-related antimicrobial peptide CRAMP protects the mouse colon from inflammation, inflammation-associated carcinogenesis, and disrupted microbiome balance, as shown in systemic Cnlp-/- mice (also known as Camp-/- mice). However, the mechanistic basis for the role and the cellular source of CRAMP in colon pathophysiology are ill defined. This study, using either epithelial or myeloid conditional Cnlp-/- mice, demonstrated that epithelial cell-derived CRAMP played a major role in supporting normal development of colon crypts, mucus production, and repair of injured mucosa. On the other hand, myeloid cell-derived CRAMP potently supported colon epithelial resistance to bacterial invasion during acute inflammation with exacerbated mucosal damage and higher rate of mouse mortality. Therefore, a well concerted cooperation of epithelial- and myeloid-derived CRAMP is essential for colon mucosal homeostasis. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

FAM3D is essential for colon homeostasis and host defense against inflammation associated carcinogenesis.

The physiological homeostasis of gut mucosal barrier is maintained by both genetic and environmental factors and its impairment leads to pathogenesis such as inflammatory bowel disease. A cytokine like molecule, FAM3D (mouse Fam3D), is highly expressed in mouse gastrointestinal tract. Here, we demonstrate that deficiency in Fam3D is associated with impaired integrity of colonic mucosa, increased epithelial hyper-proliferation, reduced anti-microbial peptide production and increased sensitivity to chemically induced colitis associated with high incidence of cancer. Pretreatment of Fam3D-/- mice with antibiotics significantly reduces the severity of chemically induced colitis and wild type (WT) mice co-housed with Fam3D-/- mice phenocopy Fam3D-deficiency showing increased sensitivity to colitis and skewed composition of fecal microbiota. An initial equilibrium of microbiota in cohoused WT and Fam3D-/- mice is followed by an increasing divergence of the bacterial composition after separation. These results demonstrate the essential role of Fam3D in colon homeostasis, protection against inflammation associated cancer and normal microbiota composition.

HLA tapasin independence: broader peptide repertoire and HIV control.

Human leukocyte antigen (HLA) class I allotypes vary in their ability to present peptides in the absence of tapasin, an essential component of the peptide loading complex. We quantified tapasin dependence of all allotypes that are common in European and African Americans (n = 97), which revealed a broad continuum of values. Ex vivo examination of cytotoxic T cell responses to the entire HIV-1 proteome from infected subjects indicates that tapasin-dependent allotypes present a more limited set of distinct peptides than do tapasin-independent allotypes, data supported by computational predictions. This suggests that variation in tapasin dependence may impact the strength of the immune responses by altering peptide repertoire size. In support of this model, we observed that individuals carrying HLA class I genotypes characterized by greater tapasin independence progress more slowly to AIDS and maintain lower viral loads, presumably due to increased breadth of peptide presentation. Thus, tapasin dependence level, like HLA zygosity, may serve as a means to restrict or expand breadth of the HLA-I peptide repertoire across humans, ultimately influencing immune responses to pathogens and vaccines.

A Prime/Boost Vaccine Regimen Alters the Rectal Microbiome and Impacts Immune Responses and Viremia Control Post-Simian Immunodeficiency Virus Infection in Male and Female Rhesus Macaques.

An efficacious human immunodeficiency virus (HIV) vaccine will likely require induction of both mucosal and systemic immune responses. We compared the immunogenicity and protective efficacy of two mucosal/systemic vaccine regimens and investigated their effects on the rectal microbiome. Rhesus macaques were primed twice mucosally with replication-competent adenovirus type 5 host range mutant (Ad5hr)-simian immunodeficiency virus (SIV) recombinants and boosted twice intramuscularly with ALVAC-SIV recombinant plus SIV gp120 protein or with DNA for SIV genes and rhesus interleukin-12 plus SIV gp120 protein. Controls received empty Ad5hr vector and alum adjuvant only. Both regimens elicited strong, comparable mucosal and systemic cellular and humoral immunity. Prevaccination rectal microbiomes of males and females differed and significantly changed over the course of immunization, most strongly in females after Ad5hr immunizations. Following repeated low-dose intrarectal SIV challenges, both vaccine groups exhibited modestly but significantly reduced acute viremia. Male and female controls exhibited similar acute viral loads; however, vaccinated females, but not males, exhibited lower levels of acute viremia, compared to same-sex controls. Few differences in adaptive immune responses were observed between the sexes. Striking differences in correlations of the rectal microbiome of males and females with acute viremia and immune responses associated with protection were seen and point to effects of the microbiome on vaccine-induced immunity and viremia control. Our study clearly demonstrates direct effects of a mucosal SIV vaccine regimen on the rectal microbiome and validates our previously reported SIV vaccine-induced sex bias. Sex and the microbiome are critical factors that should not be overlooked in vaccine design and evaluation.IMPORTANCE Differences in HIV pathogenesis between males and females, including immunity postinfection, have been well documented, as have steroid hormone effects on the microbiome, which is known to influence mucosal immune responses. Few studies have applied this knowledge to vaccine trials. We investigated two SIV vaccine regimens combining mucosal priming immunizations and systemic protein boosting. We again report a vaccine-induced sex bias, with female rhesus macaques but not males displaying significantly reduced acute viremia. The vaccine regimens, especially the mucosal primes, significantly altered the rectal microbiome. The greatest effects were in females. Striking differences between female and male macaques in correlations of prevalent rectal bacteria with viral loads and potentially protective immune responses were observed. Effects of the microbiome on vaccine-induced immunity and viremia control require further study by microbiome transfer. However, the findings presented highlight the critical importance of considering effects of sex and the microbiome in vaccine design and evaluation.

Genome-wide Association Study Identifies HLA-DPB1 as a Significant Risk Factor for Severe Aplastic Anemia.

Severe aplastic anemia (SAA) is a rare disorder characterized by hypoplastic bone marrow and progressive pancytopenia. The etiology of acquired SAA is not understood but is likely related to abnormal immune responses and environmental exposures. We conducted a genome-wide association study of individuals with SAA genetically matched to healthy controls in discovery (359 cases, 1,396 controls) and validation sets (175 cases, 1,059 controls). Combined analyses identified linked SNPs in distinct blocks within the major histocompatibility complex on 6p21. The top SNP encodes p.Met76Val in the P4 binding pocket of the HLA class II gene HLA-DPB1 (rs1042151A>G, odds ratio [OR] 1.75, 95% confidence interval [CI] 1.50-2.03, p = 1.94 × 10-13) and was associated with HLA-DP cell surface expression in healthy individuals (p = 2.04 × 10-6). Phylogenetic analyses indicate that Val76 is not monophyletic and likely occurs in conjunction with different HLA-DP binding groove conformations. Imputation of HLA-DPB1 alleles revealed increased risk of SAA associated with Val76-encoding alleles DPB1∗03:01, (OR 1.66, p = 1.52 × 10-7), DPB1∗10:01 (OR 2.12, p = 0.0003), and DPB1∗01:01 (OR 1.60, p = 0.0008). A second SNP near HLA-B, rs28367832G>A, reached genome-wide significance (OR 1.49, 95% CI 1.22-1.78, p = 7.27 × 10-9) in combined analyses; the association remained significant after excluding cases with clonal copy-neutral loss-of-heterozygosity affecting class I HLA genes (8.6% of cases and 0% of controls). SNPs in the HLA class II gene HLA-DPB1 and possibly class I (HLA-B) are associated with SAA. The replacement of Met76 to Val76 in certain HLA-DPB1 alleles might influence risk of SAA through mechanisms involving DP peptide binding specificity, expression, and/or other factors affecting DP function.

Role of HLA-B exon 1 in graft-versus-host disease after unrelated haemopoietic cell transplantation: a retrospective cohort study.

BACKGROUND: The success of unrelated haemopoietic cell transplantation (HCT) is limited by graft-versus-host disease (GVHD), which is the main post-transplantation challenge when HLA-matched donors are unavailable. A sequence dimorphism in exon 1 of HLA-B gives rise to leader peptides containing methionine (Met; M) or threonine (Thr; T), which differentially influence natural killer and T-cell alloresponses. The main aim of the study was to evaluate the role of the leader dimorphism in GVHD after HLA-B-mismatched unrelated HCT. METHODS: We did a retrospective cohort study of 33 982 patients who received an unrelated HCT done in Australia, Europe, Japan, North America, and the UK between Jan 1, 1988, and Dec 31, 2016. Data were contributed by participants of the International Histocompatibility Working Group in Hematopoietic Cell Transplantation. All cases were included and there were no exclusion criteria. Multivariate regression models were used to assess risks associated with HLA-A, HLA-B, HLA-C, HLA-DRB1, and HLA-DQB1 mismatching. Among the 33 982 transplantations, the risks of GVHD associated with HLA-B M and T leaders were established in 17 100 (50·3%) HLA-matched and 1457 (4·3%) single HLA-B-mismatched transplantations using multivariate regression models. Leader frequencies were defined in 2 004 742 BeTheMatch US registry donors. FINDINGS: Between Jan 20, 2017, and March 11, 2019, we assessed 33 982 HCTs using multivariate regression models for the role of HLA mismatching on outcome. Median follow-up was 1841 days (IQR 909-2963). Mortality and GVHD increased with increasing numbers of HLA mismatches. A single HLA-B mismatch increased grade 3-4 acute GVHD (odds ratio [OR] 1·89, 95% CI 1·53-2·33; p<0·0001). Among the single HLA-B-mismatched transplantations, acute GVHD risk was higher with leader mismatching than with leader matching (OR 1·73, 1·02-2·94; p=0·042 for grade 2-4) and with an M leader shared allotype compared with a T leader shared allotype (OR 1·98, 1·39-2·81; p=0·0001 for grade 3-4). The preferred HLA-B-mismatched donor is leader-matched and shares a T leader allotype. The majority (1 836 939 [91·6%]) of the 2 004 742 US registry donors have the TT or MT genotype. INTERPRETATION: The HLA-B leader informs GVHD risk after HLA-B-mismatched unrelated HCT and differentiates high-risk HLA-B mismatches from those with lower risk. The leader of the matched allotype could be considered to be as important as the leader of the mismatched allotype for GVHD. Prospective identification of leader-matched donors is feasible for most patients in need of a HCT, and could lower GVHD and increase availability of HCT therapy. These findings are being independently validated and warrant further research in prospective trials. FUNDING: The National Institutes of Health, USA.

The MHC in the era of next-generation sequencing: Implications for bridging structure with function.

The MHC continues to have the most disease-associations compared to other regions of the human genome, even in the genome-wide association study (GWAS) and single nucleotide polymorphism (SNP) era. Analysis of non-coding variation and their impact on the level of expression of HLA allotypes has shed new light on the potential mechanisms underlying HLA disease associations and alloreactivity in transplantation. Next-generation sequencing (NGS) technology has the capability of delineating the phase of variants in the HLA antigen-recognition site (ARS) with non-coding regulatory polymorphisms. These relationships are critical for understanding the qualitative and quantitative implications of HLA gene diversity. This article summarizes current understanding of non-coding region variation of HLA loci, the consequences of regulatory variation on HLA expression, the role for evolution in shaping lineage-specific expression, and the impact of HLA expression on disease susceptibility and transplantation outcomes. A role for phased sequencing methods for the MHC, and perspectives for future directions in basic and applied immunogenetic studies of the MHC are presented.

The Antimicrobial Peptide CRAMP Is Essential for Colon Homeostasis by Maintaining Microbiota Balance.

Commensal bacteria are critical for physiological functions in the gut, and dysbiosis in the gut may cause diseases. In this article, we report that mice deficient in cathelin-related antimicrobial peptide (CRAMP) were defective in the development of colon mucosa and highly sensitive to dextran sulfate sodium (DSS)-elicited colitis, as well as azoxymethane-mediated carcinogenesis. Pretreatment of CRAMP-/- mice with antibiotics markedly reduced the severity of DSS-induced colitis, suggesting CRAMP as a limiting factor on dysbiosis in the colon. This was supported by observations that wild-type (WT) mice cohoused with CRAMP-/- mice became highly sensitive to DSS-induced colitis, and the composition of fecal microbiota was skewed by CRAMP deficiency. In particular, several bacterial species that are typically found in oral microbiota, such as Mogibacterium neglectum, Desulfovibrio piger, and Desulfomicrobium orale, were increased in feces of CRAMP-/- mice and were transferred to WT mice during cohousing. When littermates of CRAMP+/- parents were examined, the composition of the fecal microbiota of WT pups and heterozygous parents was similar. In contrast, although the difference in fecal microbiota between CRAMP-/- and WT pups was small early on after weaning and single mouse housing, there was an increasing divergence with prolonged single housing. These results indicate that CRAMP is critical in maintaining colon microbiota balance and supports mucosal homeostasis, anti-inflammatory responses, and protection from carcinogenesis.

Posttranscriptional Regulation of HLA-A Protein Expression by Alternative Polyadenylation Signals Involving the RNA-Binding Protein Syncrip.

Genomic variation in the untranslated region (UTR) has been shown to influence HLA class I expression level and associate with disease outcomes. Sequencing of the 3'UTR of common HLA-A alleles indicated the presence of two polyadenylation signals (PAS). The proximal PAS is conserved, whereas the distal PAS is disrupted within certain alleles by sequence variants. Using 3'RACE, we confirmed expression of two distinct forms of the HLA-A 3'UTR based on use of either the proximal or the distal PAS, which differ in length by 100 bp. Specific HLA-A alleles varied in the usage of the proximal versus distal PAS, with some alleles using only the proximal PAS, and others using both the proximal and distal PAS to differing degrees. We show that the short and the long 3'UTR produced similar mRNA expression levels. However, the long 3'UTR conferred lower luciferase activity as compared with the short form, indicating translation inhibition of the long 3'UTR. RNA affinity pull-down followed by mass spectrometry analysis as well as RNA coimmunoprecipitation indicated differential binding of Syncrip to the long versus short 3'UTR. Depletion of Syncrip by small interfering RNA increased surface expression of an HLA-A allotype that uses primarily the long 3'UTR, whereas an allotype expressing only the short form was unaffected. Furthermore, specific blocking of the proximal 3'UTR reduced surface expression without decreasing mRNA expression. These data demonstrate HLA-A allele-specific variation in PAS usage, which modulates their cell surface expression posttranscriptionally.

Sequence and Phylogenetic Analysis of the Untranslated Promoter Regions for HLA Class I Genes.

Polymorphisms located within the MHC have been linked to many disease outcomes by mechanisms not yet fully understood in most cases. Variants located within untranslated regions of HLA genes are involved in allele-specific expression and may therefore underlie some of these disease associations. We determined sequences extending nearly 2 kb upstream of the transcription start site for 68 alleles from 57 major lineages of classical HLA class I genes. The nucleotide diversity within this promoter segment roughly follows that seen within the coding regions, with HLA-B showing the highest (∼1.9%), followed by HLA-A (∼1.8%), and HLA-C showing the lowest diversity (∼0.9%). Despite its greater diversity, HLA-B mRNA expression levels determined in 178 European Americans do not vary in an allele- or lineage-specific manner, unlike the differential expression levels of HLA-A or HLA-C reported previously. Close proximity of promoter sequences in phylogenetic trees is roughly reflected by similarity of expression pattern for most HLA-A and -C loci. Although promoter sequence divergence might impact promoter activity, we observed no clear link between the phylogenetic structures as represented by pairwise nucleotide differences in the promoter regions with estimated differences in mRNA expression levels for the classical class I loci. Further, no pair of class I loci showed coordinated expression levels, suggesting that distinct mechanisms across loci determine their expression level under nonstimulated conditions. These data serve as a foundation for more in-depth analysis of the functional consequences of promoter region variation within the classical HLA class I loci.

Significance of positive selection and gene duplication in adaptive evolution: in memory of Austin L. Hughes.

In this minireview, we highlight the contributions of the late Austin L. Hughes to two areas of molecular evolution: the role of positive (Darwinian) selection, and the impact of gene duplications during genome evolution.

High HLA-DP Expression and Graft-versus-Host Disease.

BACKGROUND: Transplantation of hematopoietic cells from unrelated donors can cure blood disorders but carries a significant risk of acute graft-versus-host disease (GVHD). The risk is higher when the recipient and donor are HLA-DPB1-mismatched, but the mechanisms leading to GVHD are unknown. The HLA-DPB1 regulatory region variant rs9277534 is associated with HLA-DPB1 expression. We tested the hypothesis that the GVHD risk correlates with the rs9277534 allele linked to the mismatched HLA-DPB1 in the recipient. METHODS: We genotyped rs9277534 in 3505 persons to define rs9277534-DPB1 haplotypes. Among 1441 recipients of transplants from HLA-A,B,C,DRB1,DQB1-matched unrelated donors with only one HLA-DPB1 mismatch, linkage of the rs9277534 A and G alleles to the mismatched HLA-DPB1 was determined. HLA-DPB1 expression was assessed by means of a quantitative polymerase-chain-reaction assay. The risk of acute GVHD among recipients whose mismatched HLA-DPB1 allele was linked to rs9277534G (high expression) was compared with the risk among recipients whose mismatched HLA-DPB1 allele was linked to rs9277534A (low expression). RESULTS: The mean HLA-DPB1 expression was lower with rs9277534A than with rs9277534G. Among recipients of transplants from donors with rs9277534A-linked HLA-DPB1, the risk of acute GVHD was higher for recipients with rs9277534G-linked HLA-DPB1 mismatches than for recipients with rs9277534A-linked HLA-DPB1 mismatches (hazard ratio, 1.54; 95% confidence interval [CI], 1.25 to 1.89; P<0.001), as was the risk of death due to causes other than disease recurrence (hazard ratio, 1.25; 95% CI, 1.00 to 1.57; P=0.05). CONCLUSIONS: The risk of GVHD associated with HLA-DPB1 mismatching was influenced by the HLA-DPB1 rs9277534 expression marker. Among recipients of HLA-DPB1-mismatched transplants from donors with the low-expression allele, recipients with the high-expression allele had a high risk of GVHD. (Funded by the National Institutes of Health and others.).

Sequencing rare and common APOL1 coding variants to determine kidney disease risk.

A third of African Americans with sporadic focal segmental glomerulosclerosis (FSGS) or HIV-associated nephropathy (HIVAN) do not carry APOL1 renal risk genotypes. This raises the possibility that other APOL1 variants may contribute to kidney disease. To address this question, we sequenced all APOL1 exons in 1437 Americans of African and European descent, including 464 patients with biopsy-proven FSGS/HIVAN. Testing for association with 33 common and rare variants with FSGS/HIVAN revealed no association independent of strong recessive G1 and G2 effects. Seeking additional variants that might have been under selection by pathogens and could represent candidates for kidney disease risk, we also sequenced an additional 1112 individuals representing 53 global populations. Except for G1 and G2, none of the 7 common codon-altering variants showed evidence of selection or could restore lysis against trypanosomes causing human African trypanosomiasis. Thus, only APOL1 G1 and G2 confer renal risk, and other common and rare APOL1 missense variants, including the archaic G3 haplotype, do not contribute to sporadic FSGS and HIVAN in the US population. Hence, in most potential clinical or screening applications, our study suggests that sequencing APOL1 exons is unlikely to bring additional information compared to genotyping only APOL1 G1 and G2 risk alleles.

Epigenetic regulation of differential HLA-A allelic expression levels.

MHC class I expression levels influence the strength of immune responses and represent another variable in determining outcome to disease beyond peptide binding alone. Identification of the HLA loci that vary in allelic expression levels and delineating the mechanism responsible for expression variation may provide the opportunity to modify their expression therapeutically. We have examined the expression levels of allelic lineages at the HLA-A locus in a sample of 216 European Americans using a real-time polymerase chain reaction assay, which amplifies all HLA-A lineages specifically with equal efficiency, and observed a gradient of expression that associates with HLA-A allelic lineage (R = 0.6, P = 5 × 10(-25)). DNA methylation of the HLA-A gene appears to contribute to the variation in HLA-A mRNA expression levels, as a significant inverse correlation was observed between HLA-A mRNA expression levels in untreated cells and the degree to which expression is increased after treatment of the cells with a DNA methyltransferase inhibitor (R = 0.6, P = 2.8 × 10(-6)). Further, deep-sequencing and immunoprecipitation assays revealed allelic lineage-specific methylation patterns within the HLA-A promoter region where increased DNA methylation levels correlated significantly with reduced HLA-A expression levels (R = 0.89, P = 3.7 × 10(-9)). These data demonstrate HLA-A allelic lineage-specific variation in expression levels, and DNA methylation as a likely factor in contributing to this variation.