Targeted and complete genomic sequencing of the Major Histocompatibility Complex in haplotypic form of individual heterozygous samples.

The human Major Histocompatibility Complex (MHC) is an approximately 4 Mb genomic segment on Chromosome 6 that plays a pivotal role in the immune response. Despite its importance in various traits and diseases, its complex nature makes it challenging to accurately characterize on a routine basis. We present a novel approach allowing targeted sequencing and de novo haplotypic assembly of the MHC region in heterozygous samples, using long-read sequencing technologies. Our approach is validated using two reference samples, two family trios, and an African-American sample. We achieved excellent coverage (96.6-99.9% with at least 30× depth) and high accuracy (99.89-99.99%) for the different haplotypes. This methodology offers a reliable and cost-effective method for sequencing and fully characterizing the MHC without the need for whole-genome sequencing, facilitating broader studies on this important genomic segment and having significant implications in immunology, genetics and medicine.

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.

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.

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.

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.

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.

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.

Partially CD3+-Depleted Unrelated and Haploidentical Donor Peripheral Stem Cell Transplantation Has Favorable Graft-versus-Host Disease and Survival Rates in Pediatric Hematologic Malignancy.

Most children who may benefit from stem cell transplantation lack a matched related donor. Alternative donor transplantations with an unrelated donor (URD) or a partially matched related donor (PMRD) carry an increased risk of graft-versus-host-disease (GVHD) and mortality compared with matched related donor transplantations. We hypothesized that a strategy of partial CD3+/CD19+ depletion for URD or PMRD peripheral stem cell transplantation (PSCT) would attenuate the risks of GVHD and mortality. We enrolled 84 pediatric patients with hematologic malignancies at the Children's Hospital of Philadelphia and the Children's Hospital of Wisconsin between April 2005 and February 2015. Two patients (2.4%) experienced primary graft failure. Relapse occurred in 23 patients (27.4%; cumulative incidence 26.3%), and 17 patients (20.2%) experienced nonrelapse mortality (NRM). Grade III-IV acute GVHD was observed in 18 patients (21.4%), and chronic GVHD was observed and graded as limited in 24 patients (35.3%) and extensive in 8 (11.7%). Three-year overall survival (OS) was 61.8% (95% confidence interval [CI], 50.2% to 71.4%) and event-free survival (EFS) was 52.0% (95% CI, 40.3% to 62.4%). Age ≥15 years was associated with decreased OS (P= .05) and EFS (P= .05). Relapse was more common in children in second complete remission (P = .03). Partially CD3+-depleted alternative donor PSCT NRM, OS, and EFS compare favorably with previously published studies of T cell-replete PSCT. Historically, T cell-replete PSCT has been associated with a higher incidence of extensive chronic GVHD compared with limited chronic GVHD, which may explain the comparatively low relapse and NRM rates in our study cohort despite similar overall rates of chronic GVHD. Partial T cell depletion may expand donor options for children with malignant transplantation indications lacking a matched related donor by mitigating, but not eliminating, chronic GVHD.

Assessing the utilization of high-resolution 2-field HLA typing in solid organ transplantation.

HLA typing in solid organ transplantation (SOT) is necessary for determining HLA-matching status between donor-recipient pairs and assessing patients' anti-HLA antibody profiles. Histocompatibility has traditionally been evaluated based on serologically defined HLA antigens. The evolution of HLA typing and antibody identification technologies, however, has revealed many limitations with using serologic equivalents for assessing compatibility in SOT. The significant improvements to HLA typing introduced by next-generation sequencing (NGS) require an assessment of the impact of this technology on SOT. We have assessed the role of high-resolution 2-field HLA typing (HR-2F) in SOT by retrospectively evaluating NGS-typed pre- and post-SOT cases. HR-2F typing was highly instructive or necessary in 41% (156/385) of the cases. Several pre- and posttransplant scenarios were identified as being better served by HR-2F typing. Five different categories are presented with specific case examples. The experience of another center (Temple University Hospital) is also included, whereby 21% of the cases required HR-2F typing by Sanger sequencing, as supported by other legacy methods, to properly address posttransplant anti-HLA antibody issues.

Novel findings with reassessment of exome data: implications for validation testing and interpretation of genomic data.

PurposeThe objective of this study was to assess the ability of our laboratory's exome-sequencing test to detect known and novel sequence variants and identify the critical factors influencing the interpretation of a clinical exome test.MethodsWe developed a two-tiered validation strategy: (i) a method-based approach that assessed the ability of our exome test to detect known variants using a reference HapMap sample, and (ii) an interpretation-based approach that assessed our relative ability to identify and interpret disease-causing variants, by analyzing and comparing the results of 19 randomly selected patients previously tested by external laboratories.ResultsWe demonstrate that this approach is reproducible with >99% analytical sensitivity and specificity for single-nucleotide variants and indels <10 bp. Our findings were concordant with the reference laboratories in 84% of cases. A new molecular diagnosis was applied to three cases, including discovery of two novel candidate genes.ConclusionWe provide an assessment of critical areas that influence interpretation of an exome test, including comprehensive phenotype capture, assessment of clinical overlap, availability of parental data, and the addressing of limitations in database updates. These results can be used to inform improvements in phenotype-driven interpretation of medical exomes in clinical and research settings.

Correction: Novel findings with reassessment of exome data: implications for validation testing and interpretation of genomic data.

In the published version of this article, the name of the 18th author was misspelled as Minjie Lou. The correct name is Minjie Luo. The authors regret the error.

Correction: Novel findings with reassessment of exome data: implications for validation testing and interpretation of genomic data.

In the published version of this article, the degree of author Bo Zhang was incorrectly listed as PhD. The correct degree is BS.

The impact of next-generation sequencing in immunogenetics: current status and future directions.

PURPOSE OF REVIEW: Next-generation sequencing (NGS) has now been established, and widely recognized, to be the preferred choice for human leukocyte antigen (HLA) typing. This transformation is based upon the many scientific, operational and economic benefits this technology affords. In this report, we review the major advantages, existing limitations and significant promise derived from adopting this technology in immunogenetics. RECENT FINDINGS: Significant benefits have emerged from the usage of NGS in a relatively short period, whereby we realize that this technology not only helps addressing the technical and operational problems we have had with the legacy methods for HLA typing, but equally important, it also allows for creative applications in stem cell and organ transplantation, new ways to investigate associations of the major histocompatibility complex (MHC) with many diseases and enhance our understanding regarding the MHC and non-MHC genomic interactions. The emerging picture is one of significant benefits in the diagnostic sphere of immunogenetics and transplantation and one of interconnectivity, integrating the many biological pathways controlled and affected by this unique genomic region. SUMMARY: NGS has revolutionized the science and practice of immunogenetics. In this article, we identify the still unresolved issues, the current benefits to transplantation and the potential for dissecting the complexity of the MHC, one of the most fascinating regions of the human genome. Using current trends, an attempt is made to predict future directions and outcomes.

The next generation of target capture technologies - large DNA fragment enrichment and sequencing determines regional genomic variation of high complexity.

BACKGROUND: The ability to capture and sequence large contiguous DNA fragments represents a significant advancement towards the comprehensive characterization of complex genomic regions. While emerging sequencing platforms are capable of producing several kilobases-long reads, the fragment sizes generated by current DNA target enrichment technologies remain a limiting factor, producing DNA fragments generally shorter than 1 kbp. The DNA enrichment methodology described herein, Region-Specific Extraction (RSE), produces DNA segments in excess of 20 kbp in length. Coupling this enrichment method to appropriate sequencing platforms will significantly enhance the ability to generate complete and accurate sequence characterization of any genomic region without the need for reference-based assembly. RESULTS: RSE is a long-range DNA target capture methodology that relies on the specific hybridization of short (20-25 base) oligonucleotide primers to selected sequence motifs within the DNA target region. These capture primers are then enzymatically extended on the 3'-end, incorporating biotinylated nucleotides into the DNA. Streptavidin-coated beads are subsequently used to pull-down the original, long DNA template molecules via the newly synthesized, biotinylated DNA that is bound to them. We demonstrate the accuracy, simplicity and utility of the RSE method by capturing and sequencing a 4 Mbp stretch of the major histocompatibility complex (MHC). Our results show an average depth of coverage of 164X for the entire MHC. This depth of coverage contributes significantly to a 99.94 % total coverage of the targeted region and to an accuracy that is over 99.99 %. CONCLUSIONS: RSE represents a cost-effective target enrichment method capable of producing sequencing templates in excess of 20 kbp in length. The utility of our method has been proven to generate superior coverage across the MHC as compared to other commercially available methodologies, with the added advantage of producing longer sequencing templates amenable to DNA sequencing on recently developed platforms. Although our demonstration of the method does not utilize these DNA sequencing platforms directly, our results indicate that the capture of long DNA fragments produce superior coverage of the targeted region.

Generation of Full-Length Class I Human Leukocyte Antigen Gene Consensus Sequences for Novel Allele Characterization.

BACKGROUND: Routine, high-resolution human leukocyte antigen (HLA) genotyping by next generation sequencing within clinical immunogenetics laboratories can now provide the full-length gene sequence characterization of fully phased HLA alleles. This powerful technique provides insights into HLA variation beyond the traditionally characterized antigen recognition domain, providing sequence annotation across the entire gene including untranslated and intronic regions and may be used to characterize novel alleles from massively parallel sequencing runs. METHODS: We evaluated the utility of the Omixon Holotype HLA assay to generate credible, fully phased full-length gene consensus sequences for 50 individuals at major histocompatibility complex, class I, A (HLA-A), HLA-B, and HLA-C loci (300 genotyped alleles in total) to identify and characterize novel class I HLA alleles using our downstream analytical pipeline. RESULTS: Our analysis revealed that 7.7% (23/300) of genotyped class I HLA alleles contain novel polymorphisms. Interestingly, all of the novel alleles identified by our analysis were found to harbor sequence variations within intronic regions of the respective locus. In total our analysis identified 17 unique novel class I HLA alleles from 23 of the 300 genotyped alleles and generated full-length gene sequence annotations for 9 previously incompletely annotated HLA class I allele sequences derived from 14 of the 300 genotyped alleles. CONCLUSIONS: The demonstrated utility of the Omixon Holotype HLA assay in combination with our downstream analytical framework to generate fully phased, full-length gene consensus sequences for the identification and characterization of novel HLA alleles, facilitates the study of HLA polymorphism beyond the antigen recognition domain in human health and disease.

Aneurysmal lesions of patients with abdominal aortic aneurysm contain clonally expanded T cells.

Abdominal aortic aneurysm (AAA) is a common disease with often life-threatening consequences. This vascular disorder is responsible for 1-2% of all deaths in men aged 65 years or older. Autoimmunity may be responsible for the pathogenesis of AAA. Although it is well documented that infiltrating T cells are essentially always present in AAA lesions, little is known about their role in the initiation and/or progression of the disease. To determine whether T cells infiltrating AAA lesions contain clonally expanded populations of T cells, we amplified ß-chain TCR transcripts by the nonpalindromic adaptor-PCR/Vß-specific PCR and/or Vß-specific PCR, followed by cloning and sequencing. We report in this article that aortic abdominal aneurysmal lesions from 8 of 10 patients with AAA contained oligoclonal populations of T cells. Multiple identical copies of ß-chain TCR transcripts were identified in these patients. These clonal expansions are statistically significant. These results demonstrate that αß TCR(+) T lymphocytes infiltrating aneurysmal lesions of patients with AAA have undergone proliferation and clonal expansion in vivo at the site of the aneurysmal lesion, in response to unidentified self- or nonself Ags. This evidence supports the hypothesis that AAA is a specific Ag-driven T cell disease.

A genome-wide meta-analysis of six type 1 diabetes cohorts identifies multiple associated loci.

Diabetes impacts approximately 200 million people worldwide, of whom approximately 10% are affected by type 1 diabetes (T1D). The application of genome-wide association studies (GWAS) has robustly revealed dozens of genetic contributors to the pathogenesis of T1D, with the most recent meta-analysis identifying in excess of 40 loci. To identify additional genetic loci for T1D susceptibility, we examined associations in the largest meta-analysis to date between the disease and ∼2.54 million SNPs in a combined cohort of 9,934 cases and 16,956 controls. Targeted follow-up of 53 SNPs in 1,120 affected trios uncovered three new loci associated with T1D that reached genome-wide significance. The most significantly associated SNP (rs539514, P = 5.66×10⁻¹¹) resides in an intronic region of the LMO7 (LIM domain only 7) gene on 13q22. The second most significantly associated SNP (rs478222, P = 3.50×10⁻9 resides in an intronic region of the EFR3B (protein EFR3 homolog B) gene on 2p23; however, the region of linkage disequilibrium is approximately 800 kb and harbors additional multiple genes, including NCOA1, C2orf79, CENPO, ADCY3, DNAJC27, POMC, and DNMT3A. The third most significantly associated SNP (rs924043, P = 8.06×10⁻9 lies in an intergenic region on 6q27, where the region of association is approximately 900 kb and harbors multiple genes including WDR27, C6orf120, PHF10, TCTE3, C6orf208, LOC154449, DLL1, FAM120B, PSMB1, TBP, and PCD2. These latest associated regions add to the growing repertoire of gene networks predisposing to T1D.

Coupled nanopores for single-molecule detection.

Rapid sensing of molecules is increasingly important in many studies and applications, such as DNA sequencing and protein identification. Here, beyond atomically thin 2D nanopores, we conceptualize, simulate and experimentally demonstrate coupled, guiding and reusable bilayer nanopore platforms, enabling advanced ultrafast detection of unmodified molecules. The bottom layer can collimate and decelerate the molecule before it enters the sensing zone, and the top 2D pore (~2 nm) enables position sensing. We varied the number of pores in the bottom layer from one to nine while fixing one 2D pore in the top layer. When the number of pores in the bottom layer is reduced to one, sensing is performed by both layers, and distinct T- and W-shaped translocation signals indicate the precise position of molecules and are sensitive to fragment lengths. This is uniquely enabled by microsecond resolution capabilities and precision nanofabrication. Coupled nanopores represent configurable multifunctional systems with inter- and intralayer structures for improved electromechanical control and prolonged dwell times in a 2D sensing zone.

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.