HLA Haplotypes In 250 Families: The Baylor Laboratory Results And A Perspective On A Core NGS Testing Model For The 17th International HLA And Immunogenetics Workshop.

Since their inception, the International HLA & Immunogenetics Workshops (IHIW) served as a collaborative platform for exchange of specimens, reference materials, experiences and best practices. In this report we present a subset of the results of human leukocyte antigen (HLA) haplotypes in families tested by next generation sequencing (NGS) under the 17th IHIW. We characterized 961 haplotypes in 921 subjects belonging to 250 families from 8 countries (Argentina, Austria, Egypt, Jamaica, Germany, Greece, Kuwait, and Switzerland). These samples were tested in a single core laboratory in a high throughput fashion using 6 different reagents/software platforms. Families tested included patients evaluated clinically as transplant recipients (kidney and hematopoietic cell transplant) and their respective family members. We identified 486 HLA alleles at the following loci HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQA1, -DQB1, -DPA1, -DPB1 (77, 115, 68, 69, 10, 6, 4, 44, 31, 20 and 42 alleles, respectively). We also identified nine novel alleles with polymorphisms in coding regions. This approach of testing samples from multiple laboratories across the world in different stages of technology implementation in a single core laboratory may be useful for future international workshops. Although data presented may not be reflective of allele and haplotype frequencies in the countries to which the families belong, they represent an extensive collection of 3rd and 4th field resolution level 11-locus haplotype associations of 486 alleles identified in families from 8 countries.

Assessment by Extended-Coverage Next-Generation Sequencing Typing of DPA1 and DPB1 Mismatches in Siblings Matching at HLA-A, -B, -C, -DRB1, and -DQ Loci.

Allogeneic hematopoietic stem cell transplant from an HLA matched sibling donor is usually the preferable choice. The use of next-generation sequencing (NGS) for HLA typing in clinical practice provides broader coverage and higher resolution of HLA genes. We evaluated the frequency of DPB1 crossing-over events among patients and potential related donors typed with NGS. From July 2016 to January 2018, 593 patients and 2385 siblings were typed. We evaluated sibling matching status in 546 patients, and 44.8% of these patients had siblings that matched at HLA-A, -B, -C, -DRB1, and -DQB1 loci. In 306 patient-HLA matched sibling pairs, we found 6 pairs (1.96%) with 1 DPB1 mismatch, and 5 of these pairs included an additional mismatch in DPA1. No additional mismatches were observed at the low expression loci. Using the T cell epitope algorithm, 4 of these DP mismatches were classified as permissive, 1 as nonpermissive in the host-versus-graft direction, and 1 as nonpermissive in the graft-versus-host direction. The frequency of DPB1 and DPA1 mismatches is low, and their impact in related donor transplants is not well established. Although DP typing in related transplants goes beyond guidelines, it is especially relevant for sensitized patients. NGS-based HLA typing provides full gene coverage, and its use in clinical practice can enable better donor selection.

HLA alleles and haplotypes observed in 263 US families.

The 17th International HLA and Immunogenetics Workshop (IHIW) conducted a project entitled "The Study of Haplotypes in Families by NGS HLA". We investigated the HLA haplotypes of 1017 subjects in 263 nuclear families sourced from five US clinical immunogenetics laboratories, primarily as part of the evaluation of related donor candidates for hematopoietic stem cell and solid organ transplantation. The parents in these families belonged to five broad groups - African (72 parents), Asian (115), European (210), Hispanic (118) and "Other" (11). High-resolution HLA genotypes were generated for each subject using next-generation sequencing (NGS) HLA typing systems. We identified the HLA haplotypes in each family using HaplObserve, software that builds haplotypes in families by reviewing HLA allele segregation from parents to children. We calculated haplotype frequencies within each broad group, by treating the parents in each family as unrelated individuals. We also calculated standard measures of global linkage disequilibrium (LD) and conditional asymmetric LD for each ethnic group, and used untruncated and two-field allele names to investigate LD patterns. Finally we demonstrated the utility of consensus DNA sequences in identifying novel variants, confirming them using HLA allele segregation at the DNA sequence level.

Next-generation HLA typing of 382 International Histocompatibility Working Group reference B-lymphoblastoid cell lines: Report from the 17th International HLA and Immunogenetics Workshop.

Extended molecular characterization of HLA genes in the IHWG reference B-lymphoblastoid cell lines (B-LCLs) was one of the major goals for the 17th International HLA and Immunogenetics Workshop (IHIW). Although reference B-LCLs have been examined extensively in previous workshops complete high-resolution typing was not completed for all the classical class I and class II HLA genes. To address this, we conducted a single-blind study where select panels of B-LCL genomic DNA samples were distributed to multiple laboratories for HLA genotyping by next-generation sequencing methods. Identical cell panels comprised of 24 and 346 samples were distributed and typed by at least four laboratories in order to derive accurate consensus HLA genotypes. Overall concordance rates calculated at both 2- and 4-field allele-level resolutions ranged from 90.4% to 100%. Concordance for the class I genes ranged from 91.7 to 100%, whereas concordance for class II genes was variable; the lowest observed at HLA-DRB3 (84.2%). At the maximum allele-resolution 78 B-LCLs were defined as homozygous for all 11 loci. We identified 11 novel exon polymorphisms in the entire cell panel. A comparison of the B-LCLs NGS HLA genotypes with the HLA genotypes catalogued in the IPD-IMGT/HLA Database Cell Repository, revealed an overall allele match at 68.4%. Typing discrepancies between the two datasets were mostly due to the lower-resolution historical typing methods resulting in incomplete HLA genotypes for some samples listed in the IPD-IMGT/HLA Database Cell Repository. Our approach of multiple-laboratory NGS HLA typing of the B-LCLs has provided accurate genotyping data. The data generated by the tremendous collaborative efforts of the 17th IHIW participants is useful for updating the current cell and sequence databases and will be a valuable resource for future studies.

A specific amino acid motif of HLA-DRB1 mediates risk and interacts with smoking history in Parkinson's disease.

Parkinson's disease (PD) is a neurodegenerative disease in which genetic risk has been mapped to HLA, but precise allelic associations have been difficult to infer due to limitations in genotyping methodology. Mapping PD risk at highest possible resolution, we performed sequencing of 11 HLA genes in 1,597 PD cases and 1,606 controls. We found that susceptibility to PD can be explained by a specific combination of amino acids at positions 70-74 on the HLA-DRB1 molecule. Previously identified as the primary risk factor in rheumatoid arthritis and referred to as the "shared epitope" (SE), the residues Q/R-K/R-R-A-A at positions 70-74 in combination with valine at position 11 (11-V) is highly protective in PD, while risk is attributable to the identical epitope in the absence of 11-V. Notably, these effects are modified by history of cigarette smoking, with a strong protective effect mediated by a positive history of smoking in combination with the SE and 11-V (P = 10-4; odds ratio, 0.51; 95% confidence interval, 0.36-0.72) and risk attributable to never smoking in combination with the SE without 11-V (P = 0.01; odds ratio, 1.51; 95% confidence interval, 1.08-2.12). The association of specific combinations of amino acids that participate in critical peptide-binding pockets of the HLA class II molecule implicates antigen presentation in PD pathogenesis and provides further support for genetic control of neuroinflammation in disease. The interaction of HLA-DRB1 with smoking history in disease predisposition, along with predicted patterns of peptide binding to HLA, provide a molecular model that explains the unique epidemiology of smoking in PD.

High-resolution genomic profiles of breast cancer cell lines assessed by tiling BAC array comparative genomic hybridization.

A BAC-array platform for comparative genomic hybridization was constructed from a library of 32,433 clones providing complete genome coverage, and evaluated by screening for DNA copy number changes in 10 breast cancer cell lines (BT474, MCF7, HCC1937, SK-BR-3, L56Br-C1, ZR-75-1, JIMT1, MDA-MB-231, MDA-MB-361, and HCC2218) and one cell line derived from fibrocystic disease of the breast (MCF10A). These were also characterized by gene expression analysis and found to represent all five recently described breast cancer subtypes using the "intrinsic gene set" and centroid correlation. Three cell lines, HCC1937 and L56BrC1 derived from BRCA1 mutation carriers and MDA-MB-231, were of basal-like subtype and characterized by a high frequency of low-level gains and losses of typical pattern, including limited deletions on 5q. Four estrogen receptor positive cell lines were of luminal A subtype and characterized by a different pattern of aberrations and high-level amplifications, including ERBB2 and other 17q amplicons in BT474 and MDA-MB-361. SK-BR-3 cells, characterized by a complex genome including ERBB2 amplification, massive high-level amplifications on 8q and a homozygous deletion of CDH1 at 16q22, had an expression signature closest to luminal B subtype. The effects of gene amplifications were verified by gene expression analysis to distinguish targeted genes from silent amplicon passengers. JIMT1, derived from an ERBB2 amplified trastuzumab resistant tumor, was of the ERBB2 subtype. Homozygous deletions included other known targets such as PTEN (HCC1937) and CDKN2A (MDA-MB-231, MCF10A), but also new candidate suppressor genes such as FUSSEL18 (HCC1937) and WDR11 (L56Br-C1) as well as regions without known genes. The tiling BAC-arrays constitute a powerful tool for high-resolution genomic profiling suitable for cancer research and clinical diagnostics.

Genomic analysis of Drosophila melanogaster telomeres: full-length copies of HeT-A and TART elements at telomeres.

The repetitive nature of heterochromatin hampers its analysis in general genome-sequencing projects. Specific studies are needed to extend the sequence into telomeric and centromeric heterochromatin. Drosophila telomeres lack the telomerase-generated repeats that are characteristic of other eukaryotic chromosomes. Instead, they consist of tandem arrays of HeT-A and TART elements. Herein, we present the genomic organization of the telomeres in the isogenic strain (y; cn bw sp) that was used for the Drosophila melanogaster sequencing project. The data indicate that the canonical features of telomere organization are widely conserved in evolution. In addition, we have identified full-length elements, likely competent elements, for HeT-A and TART.

Genome-wide array CGH analysis of murine neuroblastoma reveals distinct genomic aberrations which parallel those in human tumors.

Neuroblastoma, the third most common tumor of childhood, is a complex disease in which few genetic mutations have been identified.Mice expressing a human MYCN oncogene driven by the rat tyrosine hydroxylase promoter (TH-MYCN) represent an animal model for this disorder. We performed microarray-based comparative genomic hybridization analysis on murine tumors, identifying gains on chromosomes 1, 3, 11, 14, 17, and 18 and losses on chromosomes 5, 9, and 16. Fluorescence in situ hybridization analysis confirmed an amplicon on chromosome 18 as the site of TH-MYCN transgene integration. Selected tumors with localized gains of chromosome 11 delineate a 15-Mb region orthologous to human chromosome 17q and help to narrow the minimal region gained in human tumors. We observed clustered loss of chromosomes 5, 9, and 16, orthologous to a similar pattern of combined loss of chromosomes 3p, 4p, and 11q in human tumors. These data demonstrate conservation of many genetic changes in murine and human neuroblastoma and suggest that further delineation of genetic abnormalities in murine tumors may identify genes important in human disease.