High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis.

Background: The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods. Methods: We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls. Results: We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk. Results: The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

The genotype list string code syntax for exchanging nomenclature-level genotyping results in clinical and research data management and analysis systems.

The nomenclatures used to describe HLA and killer-cell immunoglobulin-like receptor (KIR) alleles distinguish unique nucleotide and peptide sequences, and patterns of expression, but are insufficient for describing genotyping results, as description of ambiguities and relations across loci require terminology beyond allele names. The genotype list (GL) String grammar describes genotyping results for genetic systems with defined nomenclatures, like HLA and KIR, documenting what is known and unknown about a given genotyping result. However, the accuracy of a GL String is dependent on the reference database version under which it was generated. Here, we describe the GL string code (GLSC) system, which associates each GL String with meta-data describing the specific reference context in which the GL String was created, and in which it should be interpreted. GLSC is a defined syntax for exchanging GL Strings in the context of a specific gene-family namespace, allele-name code-system, and pertinent reference database version. GLSC allows HLA and KIR genotyping data to be transmitted, parsed and interpreted in the appropriate context, in an unambiguous manner, on modern data-systems, including Health Level 7 Fast Healthcare Interoperability Resource systems. Technical specification for GLSC can be found at https://glstring.org.

Deconvoluting complex correlates of COVID-19 severity with a multi-omic pandemic tracking strategy.

The SARS-CoV-2 pandemic has differentially impacted populations across race and ethnicity. A multi-omic approach represents a powerful tool to examine risk across multi-ancestry genomes. We leverage a pandemic tracking strategy in which we sequence viral and host genomes and transcriptomes from nasopharyngeal swabs of 1049 individuals (736 SARS-CoV-2 positive and 313 SARS-CoV-2 negative) and integrate them with digital phenotypes from electronic health records from a diverse catchment area in Northern California. Genome-wide association disaggregated by admixture mapping reveals novel COVID-19-severity-associated regions containing previously reported markers of neurologic, pulmonary and viral disease susceptibility. Phylodynamic tracking of consensus viral genomes reveals no association with disease severity or inferred ancestry. Summary data from multiomic investigation reveals metagenomic and HLA associations with severe COVID-19. The wealth of data available from residual nasopharyngeal swabs in combination with clinical data abstracted automatically at scale highlights a powerful strategy for pandemic tracking, and reveals distinct epidemiologic, genetic, and biological associations for those at the highest risk.

A new strategy for systematically classifying HLA alleles into serological specificities.

HLA serological specificities were defined by the reactivity of HLA molecules with sets of sera and monoclonal antibodies. Many recently identified alleles defined by molecular typing lack their serotype assignment. We surveyed the literature describing the correlation of the reactivity of serologic reagents with AA residues. 20 - 25 AA residues determining epitopes (DEP) that correlated with 82 WHO serologic specificities were identified for HLA class I loci. Thirteen DEP each located in the beta-1 domains that correlated with 24 WHO serologic specificities were identified for HLA-DRB1 and -DQB1 loci. The designation of possible HLA-DPB1, -DQA1, -DPA1, and additional serological specificities that result from epitopes defined by residues located at both -DQA1 and -DQB1 subunits were also examined. HATS software was developed for automated serotype assignments to HLA alleles in one of the three hierarchical matching criteria: (1) all DEP (FULL); (2) selected DEP specific to each serological specificity (SEROTYPE); (3) one AA mismatch with one or more SEROTYPES (INCOMPLETE). Results were validated by evaluating the alleles whose serotypes do not correspond to the first field of the allele name listed in the HLA dictionary. Additional 85 and 21 DEP patterns that do not correspond to any WHO serologic specificities for common HLA class I and DRB1 alleles were identified, respectively. A comprehensive antibody identification panel would allow for accurate unacceptable antigen listing and compatibility predictions in solid organ transplantation. We propose that antibody-screening panels should include all serologic specificities identified in this study.

Challenges for the standardized reporting of NGS HLA genotyping: Surveying gaps between clinical and research laboratories.

Next generation sequencing (NGS) is being applied for HLA typing in research and clinical settings. NGS HLA typing has made it feasible to sequence exons, introns and untranslated regions simultaneously, with significantly reduced labor and reagent cost per sample, rapid turnaround time, and improved HLA genotype accuracy. NGS technologies bring challenges for cost-effective computation, data processing and exchange of NGS-based HLA data. To address these challenges, guidelines and specifications such as Genotype List (GL) String, Minimum Information for Reporting Immunogenomic NGS Genotyping (MIRING), and Histoimmunogenetics Markup Language (HML) were proposed to streamline and standardize reporting of HLA genotypes. As part of the 17th International HLA and Immunogenetics Workshop (IHIW), we implemented standards and systems for HLA genotype reporting that included GL String, MIRING and HML, and found that misunderstanding or misinterpretations of these standards led to inconsistencies in the reporting of NGS HLA genotyping results. This may be due in part to a historical lack of centralized data reporting standards in the histocompatibility and immunogenetics community. We have worked with software and database developers, clinicians and scientists to address these issues in a collaborative fashion as part of the Data Standard Hackathons (DaSH) for NGS. Here we report several categories of challenges to the consistent exchange of NGS HLA genotyping data we have observed. We hope to address these challenges in future DaSH for NGS efforts.

High Resolution Haplotype Analyses of Classical HLA Genes in Families With Multiple Sclerosis Highlights the Role of HLA-DP Alleles in Disease Susceptibility.

Multiple sclerosis (MS) susceptibility shows strong genetic associations with HLA alleles and haplotypes. We genotyped 11 HLA genes in 477 non-Hispanic European MS patients and their 954 unaffected parents using a validated next-generation sequencing (NGS) methodology. HLA haplotypes were assigned unequivocally by tracing HLA allele transmissions. We explored HLA haplotype/allele associations with MS using the genotypic transmission disequilibrium test (gTDT) and multiallelic TDT (mTDT). We also conducted a case-control (CC) study with all patients and 2029 healthy unrelated ethnically matched controls. We performed separate analyses of 54 extended multi-case families by reviewing transmission of haplotype blocks. The haplotype fragment including DRB5*01:01:01~DRB1*15:01:01:01 was significantly associated with predisposition (gTDT: p < 2.20e-16; mTDT: p =1.61e-07; CC: p < 2.22e-16) as reported previously. A second risk allele, DPB1*104:01 (gTDT: p = 3.69e-03; mTDT: p = 2.99e-03; CC: p = 1.00e-02), independent from the haplotype bearing DRB1*15:01 was newly identified. The allele DRB1*01:01:01 showed significant protection (gTDT: p = 8.68e-06; mTDT: p = 4.50e-03; CC: p = 1.96e-06). Two DQB1 alleles, DQB1*03:01 (gTDT: p = 2.86e-03; mTDT: p = 5.56e-02; CC: p = 4.08e-05) and DQB1*03:03 (gTDT: p = 1.17e-02; mTDT: p = 1.16e-02; CC: p = 1.21e-02), defined at two-field level also showed protective effects. The HLA class I block, A*02:01:01:01~C*03:04:01:01~B*40:01:02 (gTDT: p = 5.86e-03; mTDT: p = 3.65e-02; CC: p = 9.69e-03) and the alleles B*27:05 (gTDT: p = 6.28e-04; mTDT: p = 2.15e-03; CC: p = 1.47e-02) and B*38:01 (gTDT: p = 3.20e-03; mTDT: p = 6.14e-03; CC: p = 1.70e-02) showed moderately protective effects independently from each other and from the class II associated factors. By comparing statistical significance of 11 HLA loci and 19 haplotype segments with both untruncated and two-field allele names, we precisely mapped MS candidate alleles/haplotypes while eliminating false signals resulting from 'hitchhiking' alleles. We assessed genetic burden for the HLA allele/haplotype identified in this study. This family-based study including the highest-resolution of HLA alleles proved to be powerful and efficient for precise identification of HLA genotypes associated with both, susceptibility and protection to development of MS.

Standard reference sequences for submission of HLA genotyping for the 18th International HLA and Immunogenetics Workshop.

The International human leukocyte antigen (HLA) and Immunogenetics Workshops (IHIWs) have fostered international collaborations of researchers and experts in the fields of HLA, histocompatibility and immunology. These IHIW collaborations have comprised many projects focused on achieving a variety of specific goals. The international and collaborative nature of these projects necessitates the collection and analysis of complex data generated in multiple laboratories, often using multiple methods of acquisition. Collection and storage of these data in a consistent way adds value to IHIW projects, which can be extended to future work. DNA-based genotyping data, especially HLA genotyping data, can be transmitted in the form of a Histoimmunogenetics Markup Language (HML) document. HML facilitates clear communication of a genotype and supporting metadata, such as, sequencing platform, laboratory assays, consensus sequence, and interpretation. Sequence information can be reported relative to known reference sequences, which add meaning and context to genotypes. Selecting the correct reference sequence for a given allele sequence is nuanced, and guidelines have emerged through collaborative community efforts such as Data Standards Hackathons. Here, we describe the guidelines established for the selection of reference sequences to be used in transmission of HLA (and MICA/MICB) genotyping data for the 18th IHIW.

High-throughput SARS-CoV-2 and host genome sequencing from single nasopharyngeal swabs.

During COVID19 and other viral pandemics, rapid generation of host and pathogen genomic data is critical to tracking infection and informing therapies. There is an urgent need for efficient approaches to this data generation at scale. We have developed a scalable, high throughput approach to generate high fidelity low pass whole genome and HLA sequencing, viral genomes, and representation of human transcriptome from single nasopharyngeal swabs of COVID19 patients.

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.

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.

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.

Quality control project of NGS HLA genotyping for the 17th International HLA and Immunogenetics Workshop.

The 17th International HLA and Immunogenetics Workshop (IHIW) organizers conducted a Pilot Study (PS) in which 13 laboratories (15 groups) participated to assess the performance of the various sequencing library preparation protocols, NGS platforms and software in use prior to the workshop. The organizers sent 50 cell lines to each of the 15 groups, scored the 15 independently generated sets of NGS HLA genotyping data, and generated "consensus" HLA genotypes for each of the 50 cell lines. Proficiency Testing (PT) was subsequently organized using four sets of 24 cell lines, selected from 48 of 50 PS cell lines, to validate the quality of NGS HLA typing data from the 34 participating IHIW laboratories. Completion of the PT program with a minimum score of 95% concordance at the HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1 loci satisfied the requirements to submit NGS HLA typing data for the 17th IHIW projects. Together, these PS and PT efforts constituted the 17th IHIW Quality Control project. Overall PT concordance rates for HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, HLA-DRB1, HLA-DRB3, HLA-DRB4 and HLA-DRB5 were 98.1%, 97.0% and 98.1%, 99.0%, 98.6%, 98.8%, 97.6%, 96.0%, 99.1%, 90.0% and 91.7%, respectively. Across all loci, the majority of the discordance was due to allele dropout. The high cost of NGS HLA genotyping per experiment likely prevented the retyping of initially failed HLA loci. Despite the high HLA genotype concordance rates of the software, there remains room for improvement in the assembly of more accurate consensus DNA sequences by NGS HLA genotyping software.

High-resolution characterization of allelic and haplotypic HLA frequency distribution in a Spanish population using high-throughput next-generation sequencing.

Next-generation sequencing (NGS) at the HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1 and -DRB3/4/5 loci was performed on 282 healthy unrelated individuals from different major regions of Spain. High-resolution HLA genotypes defined by full sequencing of class I loci and extended coverage of class II loci were obtained to determine allele frequencies and also to estimate extended haplotype frequencies. HLA alleles were typed at the highest resolution level (4-field level, 4FL); with exception of a minor deviation in HLA-DPA1, no statistically significant deviations from expected Hardy Weinberg Equilibrium (HWE) proportions were observed for all other HLA loci. This study provides new 4FL-allele and -haplotype frequencies estimated for the first time in the Spanish population. Furthermore, our results describe extended haplotypes (including the less frequently typed HLA-DPA1 and HLA-DQA1 loci) and show distinctive haplotype associations found at 4FL-allele definition in this Spanish population study. The distinctive allelic and haplotypic diversity found at the 4FL reveals the high level of heterozygosity and specific haplotypic associations displayed that were not apparent at 2-field level (2FL). Overall, these results may contribute as a useful reference source for future population studies, for HLA-disease association studies as a healthy control group dataset and for improving donor recruitment strategies of bone marrow registries. HLA genotyping data of this Spanish population cohort was also included in the 17th International Histocompatibility and Immunogenetics Workshop (IHIW) as part of the study of HLA diversity in unrelated worldwide populations using NGS.

Tools for building, analyzing and evaluating HLA haplotypes from families.

The highly polymorphic classical human leukocyte antigen (HLA) genes display strong linkage disequilibrium (LD) that results in conserved multi-locus haplotypes. For unrelated individuals in defined populations, HLA haplotype frequencies can be estimated using the expectation-maximization (EM) method. Haplotypes can also be constructed using HLA allele segregation from nuclear families. It is straightforward to identify many HLA genotyping inconsistencies by visually reviewing HLA allele segregation in family members. It is also possible to identify potential crossover events when two or more children are available in a nuclear family. This process of visual inspection can be unwieldy, and we developed the "HaplObserve" program to standardize the process and automatically build haplotypes using family-based HLA allele segregation. HaplObserve facilitates systematically building haplotypes, and reporting potential crossover events. HLA Haplotype Validator (HLAHapV) is a program originally developed to impute chromosomal phase from genotype data using reference haplotype data. We updated and adapted HLAHapV to systematically compare observed and estimated haplotypes. We also used HLAHapV to identify haplotypes when uninformative HLA genotypes are present in families. Finally, we developed "pould", an R package that calculates haplotype frequencies, and estimates standard measures of global (locus-level) LD from both observed and estimated haplotypes.

High resolution HLA analysis reveals independent class I haplotypes and amino-acid motifs protective for multiple sclerosis.

We investigated association between HLA class I and class II alleles and haplotypes, and KIR loci and their HLA class I ligands, with multiple sclerosis (MS) in 412 European American MS patients and 419 ethnically matched controls, using next-generation sequencing. The DRB1*15:01~DQB1*06:02 haplotype was highly predisposing (odds ratio (OR) = 3.98; 95% confidence interval (CI) = 3-5.31; p-value (p) = 2.22E-16), as was DRB1*03:01~DQB1*02:01 (OR = 1.63; CI = 1.19-2.24; p = 1.41E-03). Hardy-Weinberg (HW) analysis in MS patients revealed a significant DRB1*03:01~DQB1*02:01 homozyote excess (15 observed; 8.6 expected; p = 0.016). The OR for this genotype (5.27; CI = 1.47-28.52; p = 0.0036) suggests a recessive MS risk model. Controls displayed no HW deviations. The C*03:04~B*40:01 haplotype (OR = 0.27; CI = 0.14-0.51; p = 6.76E-06) was highly protective for MS, especially in haplotypes with A*02:01 (OR = 0.15; CI = 0.04-0.45; p = 6.51E-05). By itself, A*02:01 is moderately protective, (OR = 0.69; CI = 0.54-0.87; p = 1.46E-03), and haplotypes of A*02:01 with the HLA-B Thr80 Bw4 variant (Bw4T) more so (OR = 0.53; CI = 0.35-0.78; p = 7.55E-04). Protective associations with the Bw4 KIR ligand resulted from linkage disequilibrium (LD) with DRB1*15:01, but the Bw4T variant was protective (OR = 0.64; CI = 0.49-0.82; p = 3.37-04) independent of LD with DRB1*15:01. The Bw4I variant was not associated with MS. Overall, we find specific class I HLA polymorphisms to be protective for MS, independent of the strong predisposition conferred by DRB1*15:01.

Correction: High resolution HLA analysis reveals independent class I haplotypes and amino-acid motifs protective for multiple sclerosis.

Since the publication of this article, the authors have found that the numbers of patients and controls were reversed. This study included 412 MS patients and 419 controls. This correction applies to the Abstract, the final paragraph of the Introduction, and the first paragraph of the Materials and Methods. This was entirely a reporting error and does not impact the Results or Conclusions.

Collection and storage of HLA NGS genotyping data for the 17th International HLA and Immunogenetics Workshop.

For over 50 years, the International HLA and Immunogenetics Workshops (IHIW) have advanced the fields of histocompatibility and immunogenetics (H&I) via community sharing of technology, experience and reagents, and the establishment of ongoing collaborative projects. Held in the fall of 2017, the 17th IHIW focused on the application of next generation sequencing (NGS) technologies for clinical and research goals in the H&I fields. NGS technologies have the potential to allow dramatic insights and advances in these fields, but the scope and sheer quantity of data associated with NGS raise challenges for their analysis, collection, exchange and storage. The 17th IHIW adopted a centralized approach to these issues, and we developed the tools, services and systems to create an effective system for capturing and managing these NGS data. We worked with NGS platform and software developers to define a set of distinct but equivalent NGS typing reports that record NGS data in a uniform fashion. The 17th IHIW database applied our standards, tools and services to collect, validate and store those structured, multi-platform data in an automated fashion. We have created community resources to enable exploration of the vast store of curated sequence and allele-name data in the IPD-IMGT/HLA Database, with the goal of creating a long-term community resource that integrates these curated data with new NGS sequence and polymorphism data, for advanced analyses and applications.

De Novo and Rare Variants at Multiple Loci Support the Oligogenic Origins of Atrioventricular Septal Heart Defects.

Congenital heart disease (CHD) has a complex genetic etiology, and recent studies suggest that high penetrance de novo mutations may account for only a small fraction of disease. In a multi-institutional cohort surveyed by exome sequencing, combining analysis of 987 individuals (discovery cohort of 59 affected trios and 59 control trios, and a replication cohort of 100 affected singletons and 533 unaffected singletons) we observe variation at novel and known loci related to a specific cardiac malformation the atrioventricular septal defect (AVSD). In a primary analysis, by combining developmental coexpression networks with inheritance modeling, we identify a de novo mutation in the DNA binding domain of NR1D2 (p.R175W). We show that p.R175W changes the transcriptional activity of Nr1d2 using an in vitro transactivation model in HUVEC cells. Finally, we demonstrate previously unrecognized cardiovascular malformations in the Nr1d2tm1-Dgen knockout mouse. In secondary analyses we map genetic variation to protein-interaction networks suggesting a role for two collagen genes in AVSD, which we corroborate by burden testing in a second replication cohort of 100 AVSDs and 533 controls (p = 8.37e-08). Finally, we apply a rare-disease inheritance model to identify variation in genes previously associated with CHD (ZFPM2, NSD1, NOTCH1, VCAN, and MYH6), cardiac malformations in mouse models (ADAM17, CHRD, IFT140, PTPRJ, RYR1 and ATE1), and hypomorphic alleles of genes causing syndromic CHD (EHMT1, SRCAP, BBS2, NOTCH2, and KMT2D) in 14 of 59 trios, greatly exceeding variation in control trios without CHD (p = 9.60e-06). In total, 32% of trios carried at least one putatively disease-associated variant across 19 loci,suggesting that inherited and de novo variation across a heterogeneous group of loci may contribute to disease risk.