The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection.

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

Immunotherapy targeting B cells and long-lived plasma cells effectively eliminates pre-existing donor-specific allo-antibodies.

Pre-existing anti-human leukocyte antigen (HLA) allo-antibodies constitute a major barrier to transplantation. Current desensitization approaches fail due to ineffective depletion of allo-specific memory B cells (Bmems) and long-lived plasma cells (LLPCs). We evaluate the efficacy of chimeric antigen receptor (CAR) T cells targeting CD19 and B cell maturation antigen (BCMA) to eliminate allo-antibodies in a skin pre-sensitized murine model of islet allo-transplantation. We find that treatment of allo-sensitized hosts with CAR T cells targeting Bmems and LLPCs eliminates donor-specific allo-antibodies (DSAs) and mitigates hyperacute rejection of subsequent islet allografts. We then assess the clinical efficacy of the CAR T therapy for desensitization in patients with multiple myeloma (MM) with pre-existing HLA allo-antibodies who were treated with the combination of CART-BCMA and CART-19 (ClinicalTrials.gov: NCT03549442) and observe clinically meaningful allo-antibody reduction. These findings provide logical rationale for clinical evaluation of CAR T-based immunotherapy in highly sensitized candidates to promote successful transplantation.

Genomic characterization of HLA class I and class II genes in ethnically diverse sub-Saharan African populations: A report on novel HLA alleles.

HLA allelic variation has been well studied and documented in many parts of the world. However, African populations have been relatively under-represented in studies of HLA variation. We have characterized HLA variation from 489 individuals belonging to 13 ethnically diverse populations from rural communities from the African countries of Botswana, Cameroon, Ethiopia, and Tanzania, known to practice traditional subsistence lifestyles using next generation sequencing (Illumina) and long-reads from Oxford Nanopore Technologies. We identified 342 distinct alleles among the 11 HLA targeted genes: HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQA1, -DQB1, -DPA1, and -DPB1, with 140 of those alleles containing novel sequences that were submitted to the IPD-IMGT/HLA database. Sixteen of the 140 alleles contained novel content within the exonic regions of the genes, while 110 alleles contained novel intronic variants. Four alleles were found to be recombinants of already described HLA alleles and 10 alleles extended the sequence content of already described alleles. All 140 alleles include complete allelic sequence from the 5' UTR to the 3' UTR that are inclusive of all exons and introns. This report characterizes the HLA allelic variation from these individuals and describes the novel allelic variation present within these specific African populations.

A combination of HLA-DP α and ß chain polymorphisms paired with a SNP in the DPB1 3' UTR region, denoting expression levels, are associated with atopic dermatitis.

Introduction: Components of the immune response have previously been associated with the pathophysiology of atopic dermatitis (AD), specifically the Human Leukocyte Antigen (HLA) Class II region via genome-wide association studies, however the exact elements have not been identified. Methods: This study examines the genetic variation of HLA Class II genes using next generation sequencing (NGS) and evaluates the resultant amino acids, with particular attention on binding site residues, for associations with AD. The Genetics of AD cohort was used to evaluate HLA Class II allelic variation on 464 subjects with AD and 384 controls. Results: Statistically significant associations with HLA-DP α and ß alleles and specific amino acids were found, some conferring susceptibility to AD and others with a protective effect. Evaluation of polymorphic residues in DP binding pockets revealed the critical role of P1 and P6 (P1: α31M + (ß84G or ß84V) [protection]; α31Q + ß84D [susceptibility] and P6: α11A + ß11G [protection]) and were replicated with a national cohort of children consisting of 424 AD subjects. Independently, AD susceptibility-associated residues were associated with the G polymorphism of SNP rs9277534 in the 3' UTR of the HLA-DPB1 gene, denoting higher expression of these HLA-DP alleles, while protection-associated residues were associated with the A polymorphism, denoting lower expression. Discussion: These findings lay the foundation for evaluating non-self-antigens suspected to be associated with AD as they potentially interact with particular HLA Class II subcomponents, forming a complex involved in the pathophysiology of AD. It is possible that a combination of structural HLA-DP components and levels of expression of these components contribute to AD pathophysiology.

Identification of HLA-DPA1*01:03:01:57 and HLA-DPA1*02:01:01:29 from a case-control study of atopic dermatitis.

DPA1*01:03:01:57 and DPA1*02:01:01:29 differ by a single nucleotide from their closest references, DPA1*01:03:01:02 and DPA1*02:01:01:06.

Human abdominal aortic aneurysm (AAA): Evidence for an autoimmune antigen-driven disease.

Abdominal aortic aneurism (AAA) is a complex immunological disease with a strong genetic component, and one of the ten leading causes of death of individuals 55-74 years old worldwide. Strong evidence has been accumulated suggesting that AAA is an autoimmune specific antigen-driven disease. Mononuclear cells infiltrating AAA lesions comprised of T and B lymphocytes and other cells expressing early-, intermediate- and late-activation antigens, and the presence of antigen-presenting cells have been documented, demonstrating an ongoing immune response. The three components of the trimolecular complex, T-cell receptor (TCR)/peptide (antigen)/HLA have been identified in AAA, and specifically: (i) clonal expansions of T-cell clones in AAA lesions; (ii) the association of AAA with particular HLA Class I and Class II; and (iii) self or nonself putative AAA-associated antigens. IgG autoantibodies recognizing proteins present in normal aortic tissue have been reported in patients with AAA. Molecular mimicry, defined as the sharing of antigenic epitopes between microorganisms (bacteria, viruses) and self antigens, maybe is responsible for T-cell responses and antibody production in AAA. Also, the frequency and the suppressor activity of CD4+ CD25+ FOXP3+ Tregs and the expression of FOXP3 transcripts and protein have been reported to be significantly impaired in AAA patients vs normal donors.

Severe delayed hypersensitivity reactions to IL-1 and IL-6 inhibitors link to common HLA-DRB1*15 alleles.

OBJECTIVES: Drug reaction with eosinophilia and systemic symptoms (DRESS) is a severe, delayed hypersensitivity reaction (DHR). We observed DRESS to inhibitors of interleukin 1 (IL-1) or IL-6 in a small group of patients with Still's disease with atypical lung disease. We sought to characterise features of patients with Still's disease with DRESS compared with drug-tolerant Still's controls. We analysed human leucocyte antigen (HLA) alleles for association to inhibitor-related DHR, including in a small Kawasaki disease (KD) cohort. METHODS: In a case/control study, we collected a multicentre series of patients with Still's disease with features of inhibitor-related DRESS (n=66) and drug-tolerant Still's controls (n=65). We retrospectively analysed clinical data from all Still's subjects and typed 94/131 for HLA. European Still's-DRESS cases were ancestry matched to International Childhood Arthritis Genetics Consortium paediatric Still's cases (n=550) and compared for HLA allele frequencies. HLA association also was analysed using Still's-DRESS cases (n=64) compared with drug-tolerant Still's controls (n=30). KD subjects (n=19) were similarly studied. RESULTS: Still's-DRESS features included eosinophilia (89%), AST-ALT elevation (75%) and non-evanescent rash (95%; 88% involving face). Macrophage activation syndrome during treatment was frequent in Still's-DRESS (64%) versus drug-tolerant Still's (3%; p=1.2×10-14). We found striking enrichment for HLA-DRB1*15 haplotypes in Still's-DRESS cases versus INCHARGE Still's controls (p=7.5×10-13) and versus self-identified, ancestry-matched Still's controls (p=6.3×10-10). In the KD cohort, DRB1*15:01 was present only in those with suspected anakinra reactions. CONCLUSIONS: DRESS-type reactions occur among patients treated with IL-1/IL-6 inhibitors and strongly associate with common HLA-DRB1*15 haplotypes. Consideration of preprescription HLA typing and vigilance for serious reactions to these drugs are warranted.

Unrelated donor α/ß T cell- and B cell-depleted HSCT for the treatment of pediatric acute leukemia.

Unrelated donor (URD) hematopoietic stem cell transplant (HSCT) is associated with an increased risk of severe graft-versus-host disease (GVHD). TCRαß/CD19 depletion may reduce this risk, whereas maintaining graft-versus-leukemia. Outcome data with TCRαß/CD19 depletion generally describe haploidentical donors, with relatively few URDs. We hypothesized that TCRαß/CD19-depletion would attenuate the risks of GVHD and relapse for URD HSCT. Sixty pediatric and young adult (YA) patients with hematologic malignancies who lacked a matched-related donor were enrolled at 2 large pediatric transplantation centers between October 2014 and September 2019. All patients with acute leukemia had minimal residual disease testing, and DP typing was available for 77%. All patients received myeloablative total body irradiation- or busulfan-based conditioning with no posttransplant immune suppression. Engraftment occurred in 98%. Four-year overall survival was 69% (95% confidence interval [CI], 52%-81%), and leukemia-free survival was 64% (95% CI, 48%-76%), with no difference between lymphoid and myeloid malignancies (P = .6297 and P = .5441, respectively). One patient (1.7%) experienced primary graft failure. Relapse occurred in 11 patients (3-year cumulative incidence, 21%; 95% CI, 11-34), and 8 patients (cumulative incidence, 15%; 95% CI, 6.7-26) experienced nonrelapse mortality. Grade III to IV acute GVHD was seen in 8 patients (13%), and 14 patients (26%) developed chronic GVHD, of which 6 (11%) had extensive disease. Nonpermissive DP mismatch was associated with higher likelihood of acute GVHD (odds ratio, 16.50; 95% CI, 1.67-163.42; P = .0166) but not with the development of chronic GVHD. URD TCRαß/CD19-depleted peripheral HSCT is a safe and effective approach to transplantation for children/YAs with leukemia. This trial was registered at www.clinicaltrials.gov as #NCT02323867.

Accurate long-read sequencing allows assembly of the duplicated RHD and RHCE genes harboring variants relevant to blood transfusion.

Next-generation sequencing (NGS) technologies have transformed medical genetics. However, short-read lengths pose a limitation on identification of structural variants, sequencing repetitive regions, phasing of distant nucleotide changes, and distinguishing highly homologous genomic regions. Long-read sequencing technologies may offer improvements in the characterization of genes that are currently difficult to assess. We used a combination of targeted DNA capture, long-read sequencing, and a customized bioinformatics pipeline to fully assemble the RH region, which harbors variation relevant to red cell donor-recipient mismatch, particularly among patients with sickle cell disease. RHD and RHCE are a pair of duplicated genes located within an ∼175 kb region on human chromosome 1 that have high sequence similarity and frequent structural variations. To achieve the assembly, we utilized palindrome repeats in PacBio SMRT reads to obtain consensus sequences of 2.1 to 2.9 kb average length with over 99% accuracy. We used these long consensus sequences to identify 771 assembly markers and to phase the RHD-RHCE region with high confidence. The dataset enabled direct linkage between coding and intronic variants, phasing of distant SNPs to determine RHD-RHCE haplotypes, and identification of known and novel structural variations along with the breakpoints. A limiting factor in phasing is the frequency of heterozygous assembly markers and therefore was most successful in samples from African Black individuals with increased heterogeneity at the RH locus. Overall, this approach allows RH genotyping and de novo assembly in an unbiased and comprehensive manner that is necessary to expand application of NGS technology to high-resolution RH typing.

HLA-DPB1 E69 genotype and exposure in beryllium sensitisation and disease.

OBJECTIVES: Human leukocyte antigen-DP beta 1 (HLA-DPB1) with a glutamic acid at the 69th position of the ß chain (E69) genotype and inhalational beryllium exposure individually contribute to risk of chronic beryllium disease (CBD) and beryllium sensitisation (BeS) in exposed individuals. This retrospective nested case-control study assessed the contribution of genetics and exposure in the development of BeS and CBD. METHODS: Workers with BeS (n=444), CBD (n=449) and beryllium-exposed controls (n=890) were enrolled from studies conducted at nuclear weapons and primary beryllium manufacturing facilities. Lifetime-average beryllium exposure estimates were based on workers' job questionnaires and historical and industrial hygienist exposure estimates, blinded to genotype and case status. Genotyping was performed using sequence-specific primer-PCR. Logistic regression models were developed allowing for over-dispersion, adjusting for workforce, race, sex and ethnicity. RESULTS: Having no E69 alleles was associated with lower odds of both CBD and BeS; every additional E69 allele increased odds for CBD and BeS. Increasing exposure was associated with lower odds of BeS. CBD was not associated with exposure as compared to controls, yet the per cent of individuals with CBD versus BeS increased with increasing exposure. No evidence of a gene-by-exposure interaction was found for CBD or BeS. CONCLUSIONS: Risk of CBD increases with E69 allele frequency and increasing exposure, although no gene by environment interaction was found. A decreased risk of BeS with increasing exposure and lack of exposure response in CBD cases may be due to the limitations of reconstructed exposure estimates. Although reducing exposure may not prevent BeS, it may reduce CBD and the associated health effects, especially in those carrying E69 alleles.

Association of KIR Genes and MHC Class I Ligands with Atopic Dermatitis.

Atopic dermatitis (AD) is a chronic illness that is associated with immune dysregulation. NK cell function has previously been associated with AD. NK cells directly interact with polymorphic HLA class I ligand variants using killer cell Ig-like receptors (KIRs). The purpose of this study was to identify potential associations between NK cell function and AD by evaluating variation in the presence of KIR genes as well as KIR gene interactions with the appropriate HLA class I KIR-specific ligands. Human DNA from the genetics of AD case-control study was used to genotype HLA class I KIR-specific ligands and the presence of KIR genes. In the full cohort, an increased risk of AD was noted for KIR2DL5 (1.51 [1.13, 2.01]), KIR2DS5 (1.72 [1.26, 2.34]), and KIR2DS1 (1.41 [1.04, 1.91]). Individuals with KIR2DS5 or KIR2DS1 and the HLA-C*C2 epitope were at an increased risk of AD (1.74 [1.21, 2.51] and 1.48 [1.04, 2.12], respectively). The HLA-B*-21T (TT) leader sequence increased the risk of AD across ethnicity. African Americans with KIR2DL2, KIR2DS1, KIR2DL5, and KIR2DS5 are more likely to have AD, and the risk increased for KIR2DS1 and KIR2DS5 in the presence of appropriate HLA-C C2 epitope. The risk of AD also increased for individuals with the HLA-B*-21T leader sequence. Future studies should focus on KIR gene allelic variation as well as consider cell-based measurements of KIR and the associated HLA class I epitopes.

Analysis of single nucleotide polymorphisms in chronic beryllium disease.

Sarcoidosis and chronic beryllium disease (CBD) are phenocopies, however the latter one has a clear trigger factor that is beryllium exposure. This study analyses single nucleotide polymorphisms (SNPs) in a large cohort for beryllium-exposed persons. SNPs were chosen for their relevance in sarcoidosis. Even though one of largest cohorts of beryllium-exposed persons was analysed, no statistically relevant association between any SNP and CBD could be verified. Notably, some SNPs exhibit inverse OR for beryllium sensitization and CBD with nominally statistical significance, which allows hypothesizing about pathophysiological role of genes for the disease triggering and development.

HLA Class I Polymorphisms Influencing Both Peptide Binding and KIR Interactions Are Associated with Remission among Children with Atopic Dermatitis: A Longitudinal Study.

Atopic dermatitis (AD) is a disease of immune dysregulation and skin barrier dysfunction with a relapsing, remitting course and has been associated with several different genetic risk variants. HLA represent a highly variable set of genes that code for cell surface protein molecules involved in the Ag-specific immune response, including the regulation or functioning of T cells, NK cells, and APCs. The purpose of this study was to evaluate associations between HLA class I polymorphisms and the progression of AD over time. We evaluated the associations of AD symptoms and HLA class I polymorphisms based on high-resolution two-field typing in a longitudinal cohort of children with AD (up to 10 y of follow-up). Seven hundred and ninety-two children were evaluated every 6 mo, resulting in 12,752 AD evaluations. Using generalized estimating equations and corrected p values, B*44:02 was found to be associated with AD remission (1.83 [1.35, 2.47]; p = 0.0015). The HLA-B residues at position 116 (d-aspartate) and 80 (T-threonine) were associated with remission (1.42 [1.13, 1.76], p = 0.003; corrected p = 0.028) and (1.45 [1.17, 1.80], p = 0.0008; corrected p = 0.0024), respectively. B80T is a killer-cell Ig-like receptor (KIR) site. Our findings reveal that two axes of immune response (T cell and NK cell) may influence disease progression. Identifying binding pocket changes in addition to other factors (e.g., allergens) that increase the risk or severity of AD can improve our understanding of the immunologic mechanisms associated with AD and may lead to personalized therapies for improving patient care.

Immune Reconstitution Following TCRαß/CD19-Depleted Hematopoietic Cell Transplantation for Hematologic Malignancy in Pediatric Patients.

TCRαß/CD19-depleted HCT has been used with excellent outcomes in pediatric patients with hematologic malignancies, and several studies have demonstrated rapid immune reconstitution in the nonmalignant setting. However, immune recovery following TCRαß/CD19-depleted hematopoietic cell transplantation (HCT) for malignancy remains incompletely elucidated. Furthermore, the majority of studies to date have used haploidentical and matched unrelated donors. Here we report results of immune reconstitution following TCRαß/CD19-depleted HCT for hematologic malignancy in 51 pediatric patients with hematologic malignancies, the majority of whom received grafts from unrelated donors. Grafts were from matched unrelated (n = 20), mismatched unrelated (n = 20), and haploidentical (n = 11) donors. The median CD34+ cell dose was 10.2 × 106/kg (range, 4.54 to 20 × 106/kg), and the median TCRαß+ cell dose was 2.53 × 104/kg (range, 0 to 44.9 × 104/kg). Conditioning was myeloablative with either busulfan or total body irradiation, cyclophosphamide, and thiotepa. Thirty-three patients also received rabbit antithymocyte globulin. No prophylactic post-transplantation immune suppression was routinely given. Forty-three patients received rituximab on day +1 for recipient positive Epstein-Barr virus serology. Forty-nine patients (96%) engrafted with a median time to neutrophil recovery of 13 days (range, 8 to 30 days). Thirty-seven patients (73%) are alive at a median follow-up of 25 months (range, 6 to 50 months). Nine patients (18%) developed grade II-IV acute graft-versus-host disease (GVHD), and 5 patients (11%) developed extensive chronic GVHD. Twenty-six patients (51%) experienced viral reactivation. T cell reconstitution was rapid with significant numbers of CD3+, CD4+, and CD8+ T cells present on first assessment at 4 months post-HCT, and significant numbers of naïve CD4+ T cells were present by 8 months post-HCT. Chronic GVHD was associated with delayed T cell recovery; however, T cell reconstitution was not affected by underlying diagnosis, donor source, TCRαß+ T cell dose, conditioning regimen, or use of antithymocyte globulin. B cell recovery mirrored T cell recovery, and i.v. Ig was discontinued at a median of 8 months (range, 4 to 22 months) post-HCT in patients alive and relapse-free at last follow-up. Immune reconstitution is rapid following TCRαß/CD19-depleted HCT in pediatric patients with hematologic malignancies. Donor graft source, haploidentical or unrelated, did not affect immune reconstitution. Viral reactivation is common in the first 100 days post-HCT, indicating that improved T cell defense is needed in the early post-HCT period.

Next-generation sequencing technologies: An overview.

Since the days of Sanger sequencing, next-generation sequencing technologies have significantly evolved to provide increased data output, efficiencies, and applications. These next generations of technologies can be categorized based on read length. This review provides an overview of these technologies as two paradigms: short-read, or "second-generation," technologies, and long-read, or "third-generation," technologies. Herein, short-read sequencing approaches are represented by the most prevalent technologies, Illumina and Ion Torrent, and long-read sequencing approaches are represented by Pacific Biosciences and Oxford Nanopore technologies. All technologies are reviewed along with reported advantages and disadvantages. Until recently, short-read sequencing was thought to provide high accuracy limited by read-length, while long-read technologies afforded much longer read-lengths at the expense of accuracy. Emerging developments for third-generation technologies hold promise for the next wave of sequencing evolution, with the co-existence of longer read lengths and high accuracy.

Spatial constrains and information content of sub-genomic regions of the human genome.

Complexity metrics and machine learning (ML) models have been utilized to analyze the lengths of segmental genomic entities of DNA sequences (exonic, intronic, intergenic, repeat, unique) with the purpose to ask questions regarding the segmental organization of the human genome within the size distribution of these sequences. For this we developed an integrated methodology that is based upon the reconstructed phase space theorem, the non-extensive statistical theory of Tsallis, ML techniques, and a technical index, integrating the generated information, which we introduce and named complexity factor (COFA). Our analysis revealed that the size distribution of the genomic regions within chromosomes are not random but follow patterns with characteristic features that have been seen through its complexity character, and it is part of the dynamics of the whole genome. Finally, this picture of dynamics in DNA is recognized using ML tools for clustering, classification, and prediction with high accuracy.

Complex Linkage Disequilibrium Effects in HLA-DPB1 Expression and Molecular Mismatch Analyses of Transplantation Outcomes.

BACKGROUND: HLA molecular mismatch (MM) is a risk factor for de novo donor-specific antibody (dnDSA) development in solid organ transplantation. HLA expression differences have also been associated with adverse outcomes in hematopoietic cell transplantation. We sought to study both MM and expression in assessing dnDSA risk. METHODS: One hundred three HLA-DP-mismatched solid organ transplantation pairs were retrospectively analyzed. MM was computed using amino acids (aa), eplets, and, supplementarily, Grantham/Epstein scores. DPB1 alleles were classified as rs9277534-A (low-expression) or rs9277534-G (high-expression) linked. To determine the associations between risk factors and dnDSA, logistic regression, linkage disequilibrium (LD), and population-based analyses were performed. RESULTS: A high-risk AA:GX (recipient:donor) expression combination (X = A or G) demonstrated strong association with HLA-DP dnDSA (P = 0.001). MM was also associated with HLA-DP dnDSA when evaluated by itself (eplet P = 0.007, aa P = 0.003, Grantham P = 0.005, Epstein P = 0.004). When attempting to determine the relative individual effects of the risk factors in multivariable analysis, only AA:GX expression status retained a strong association (relative risk = 18.6, P = 0.007 with eplet; relative risk = 15.8, P = 0.02 with aa), while MM was no longer significant (eplet P = 0.56, aa P = 0.51). Importantly, these risk factors are correlated, due to LD between the expression-tagging single-nucleotide polymorphism and polymorphisms along HLA-DPB1. CONCLUSIONS: The MM and expression risk factors each appear to be strong predictors of HLA-DP dnDSA and to possess clinical utility; however, these two risk factors are closely correlated. These metrics may represent distinct ways of characterizing a common overlapping dnDSA risk profile, but they are not independent. Further, we demonstrate the importance and detailed implications of LD effects in dnDSA risk assessment and possibly transplantation overall.