Haplotype reconstruction for genetically complex regions with ambiguous genotype calls: Illustration by the KIR gene region.

Advances in DNA sequencing technologies have enabled genotyping of complex genetic regions exhibiting copy number variation and high allelic diversity, yet it is impossible to derive exact genotypes in all cases, often resulting in ambiguous genotype calls, that is, partially missing data. An example of such a gene region is the killer-cell immunoglobulin-like receptor (KIR) genes. These genes are of special interest in the context of allogeneic hematopoietic stem cell transplantation. For such complex gene regions, current haplotype reconstruction methods are not feasible as they cannot cope with the complexity of the data. We present an expectation-maximization (EM)-algorithm to estimate haplotype frequencies (HTFs) which deals with the missing data components, and takes into account linkage disequilibrium (LD) between genes. To cope with the exponential increase in the number of haplotypes as genes are added, we add three components to a standard EM-algorithm implementation. First, reconstruction is performed iteratively, adding one gene at a time. Second, after each step, haplotypes with frequencies below a threshold are collapsed in a rare haplotype group. Third, the HTF of the rare haplotype group is profiled in subsequent iterations to improve estimates. A simulation study evaluates the effect of combining information of multiple genes on the estimates of these frequencies. We show that estimated HTFs are approximately unbiased. Our simulation study shows that the EM-algorithm is able to combine information from multiple genes when LD is high, whereas increased ambiguity levels increase bias. Linear regression models based on this EM, show that a large number of haplotypes can be problematic for unbiased effect size estimation and that models need to be sparse. In a real data analysis of KIR genotypes, we compare HTFs to those obtained in an independent study. Our new EM-algorithm-based method is the first to account for the full genetic architecture of complex gene regions, such as the KIR gene region. This algorithm can handle the numerous observed ambiguities, and allows for the collapsing of haplotypes to perform implicit dimension reduction. Combining information from multiple genes improves haplotype reconstruction.

HLA-DP on Epithelial Cells Enables Tissue Damage by NKp44+ Natural Killer Cells in Ulcerative Colitis.

BACKGROUND & AIMS: Ulcerative colitis (UC) is characterized by severe inflammation and destruction of the intestinal epithelium, and is associated with specific risk single nucleotide polymorphisms in HLA class II. Given the recently discovered interactions between subsets of HLA-DP molecules and the activating natural killer (NK) cell receptor NKp44, genetic associations of UC and HLA-DP haplotypes and their functional implications were investigated. METHODS: HLA-DP haplotype and UC risk association analyses were performed (UC: n = 13,927; control: n = 26,764). Expression levels of HLA-DP on intestinal epithelial cells (IECs) in individuals with and without UC were quantified. Human intestinal 3-dimensional (3D) organoid cocultures with human NK cells were used to determine functional consequences of interactions between HLA-DP and NKp44. RESULTS: These studies identified HLA-DPA1∗01:03-DPB1∗04:01 (HLA-DP401) as a risk haplotype and HLA-DPA1∗01:03-DPB1∗03:01 (HLA-DP301) as a protective haplotype for UC in European populations. HLA-DP expression was significantly higher on IECs of individuals with UC compared with controls. IECs in human intestinal 3D organoids derived from HLA-DP401pos individuals showed significantly stronger binding of NKp44 compared with HLA-DP301pos IECs. HLA-DP401pos IECs in organoids triggered increased degranulation and tumor necrosis factor production by NKp44+ NK cells in cocultures, resulting in enhanced epithelial cell death compared with HLA-DP301pos organoids. Blocking of HLA-DP401-NKp44 interactions (anti-NKp44) abrogated NK cell activity in cocultures. CONCLUSIONS: We identified an UC risk HLA-DP haplotype that engages NKp44 and activates NKp44+ NK cells, mediating damage to intestinal epithelial cells in an HLA-DP haplotype-dependent manner. The molecular interaction between NKp44 and HLA-DP401 in UC can be targeted by therapeutic interventions to reduce NKp44+ NK cell-mediated destruction of the intestinal epithelium in UC.

Engagement of TRAIL triggers degranulation and IFNγ production in human natural killer cells.

NK cells utilize a large array of receptors to screen their surroundings for aberrant or virus-infected cells. Given the vast diversity of receptors expressed on NK cells we seek to identify receptors involved in the recognition of HIV-1-infected cells. By combining an unbiased large-scale screening approach with a functional assay, we identify TRAIL to be associated with NK cell degranulation against HIV-1-infected target cells. Further investigating the underlying mechanisms, we demonstrate that TRAIL is able to elicit multiple effector functions in human NK cells independent of receptor-mediated induction of apoptosis. Direct engagement of TRAIL not only results in degranulation but also IFNγ production. Moreover, TRAIL-mediated NK cell activation is not limited to its cognate death receptors but also decoy receptor I, adding a new perspective to the perceived regulatory role of decoy receptors in TRAIL-mediated cytotoxicity. Based on these findings, we propose that TRAIL not only contributes to the anti-HIV-1 activity of NK cells but also possesses a multifunctional role beyond receptor-mediated induction of apoptosis, acting as a regulator for the induction of different effector functions.

Post-COVID-19 condition in the German working population: A cross-sectional study of 200,000 registered stem cell donors.

BACKGROUND: The SARS-CoV-2 pandemic has strained health systems worldwide, and infection numbers continue to rise. While previous data have already shown that many patients suffer from symptoms for months after an acute infection, data on risk factors and long-term outcomes are incomplete, particularly for the working population. OBJECTIVES: We aimed to provide information on the prevalence of post-COVID-19 conditions in a subset of the German working-age population (18-61 years old) and to analyze risk factors. METHODS: We conducted an online survey with a health questionnaire among registered potential stem cell donors with or without a self-reported history of polymerase chain reaction (PCR)-confirmed SARS-CoV-2 infection. Logistic regression models were used to examine the risks of severity of acute infection, sex, age, body mass index, diabetes mellitus, and arterial hypertension medication on post-COVID-19 symptoms. RESULTS: A total of 199,377 donors reported evaluable survey questionnaires-12,609 cases had a history of SARS-CoV-2 infection and 186,768 controls had none. Overall, cases reported physical, cognitive, and psychological complaints more frequently compared to controls. Increased rates of complaints persisted throughout 15 months postinfection, for example, 28.4%/19.3% of cases/controls reported fatigue (p <0.0001) and 9.5%/3.6% of cases/controls reported loss of concentration (p <0.0001). No significant differences were observed in the frequency of reported symptoms between 3 and 15 months postinfection. Multivariate analysis revealed a strong influence of the severity of the acute SARS-CoV-2 infection episode and age on the risk for post-COVID-19 conditions. CONCLUSION: We report the prevalence of post-COVID-19 conditions in mainly unvaccinated individuals with SARS-CoV-2 infections between February 2020 and August 2021. The severity of the acute course and age were major risk factors. Vaccinations may reduce the risk of post-COVID-19 conditions by reducing the risk of severe infections.

Analysis of Genomic DNA from Medieval Plague Victims Suggests Long-Term Effect of Yersinia pestis on Human Immunity Genes.

Pathogens and associated outbreaks of infectious disease exert selective pressure on human populations, and any changes in allele frequencies that result may be especially evident for genes involved in immunity. In this regard, the 1346-1353 Yersinia pestis-caused Black Death pandemic, with continued plague outbreaks spanning several hundred years, is one of the most devastating recorded in human history. To investigate the potential impact of Y. pestis on human immunity genes, we extracted DNA from 36 plague victims buried in a mass grave in Ellwangen, Germany in the 16th century. We targeted 488 immune-related genes, including HLA, using a novel in-solution hybridization capture approach. In comparison with 50 modern native inhabitants of Ellwangen, we find differences in allele frequencies for variants of the innate immunity proteins Ficolin-2 and NLRP14 at sites involved in determining specificity. We also observed that HLA-DRB1*13 is more than twice as frequent in the modern population, whereas HLA-B alleles encoding an isoleucine at position 80 (I-80+), HLA C*06:02 and HLA-DPB1 alleles encoding histidine at position 9 are half as frequent in the modern population. Simulations show that natural selection has likely driven these allele frequency changes. Thus, our data suggest that allele frequencies of HLA genes involved in innate and adaptive immunity responsible for extracellular and intracellular responses to pathogenic bacteria, such as Y. pestis, could have been affected by the historical epidemics that occurred in Europe.

External validation of models for KIR2DS1/KIR3DL1-informed selection of hematopoietic cell donors fails.

Several studies suggest that harnessing natural killer (NK) cell reactivity mediated through killer cell immunoglobulin-like receptors (KIRs) could reduce the risk of relapse after allogeneic hematopoietic cell transplantation. Based on one promising model, information on KIR2DS1 and KIR3DL1 and their cognate ligands can be used to classify donors as KIR-advantageous or KIR-disadvantageous. This study was aimed at externally validating this model in unrelated donor hematopoietic cell transplantation. The impact of the predictor on overall survival (OS) and relapse incidence was tested in a Cox regression model adjusted for patient age, a modified disease risk index, Karnofsky performance status, donor age, HLA match, sex match, cytomegalovirus match, conditioning intensity, type of T-cell depletion, and graft type. Data from 2222 patients with acute myeloid leukemia or myelodysplastic syndrome were analyzed. KIR genes were typed by using high-resolution amplicon-based next-generation sequencing. In univariable analyses and subgroup analyses, OS and the cumulative incidence of relapse of patients with a KIR-advantageous donor were comparable to patients with a KIR-disadvantageous donor. The adjusted hazard ratio from the multivariable Cox regression model was 0.99 (Wald test, P = .93) for OS and 1.04 (Wald test, P = .78) for relapse incidence. We also tested the impact of activating donor KIR2DS1 and inhibition by KIR3DL1 separately but found no significant impact on OS and the risk of relapse. Thus, our study shows that the proposed model does not universally predict NK-mediated disease control. Deeper knowledge of NK-mediated alloreactivity is necessary to predict its contribution to graft-versus-leukemia reactions and to eventually use KIR genotype information for donor selection.

High-throughput Interpretation of Killer-cell Immunoglobulin-like Receptor Short-read Sequencing Data with PING.

The killer-cell immunoglobulin-like receptor (KIR) complex on chromosome 19 encodes receptors that modulate the activity of natural killer cells, and variation in these genes has been linked to infectious and autoimmune disease, as well as having bearing on pregnancy and transplant outcomes. The medical relevance and high variability of KIR genes makes short-read sequencing an attractive technology for interrogating the region, providing a high-throughput, high-fidelity sequencing method that is cost-effective. However, because this gene complex is characterized by extensive nucleotide polymorphism, structural variation including gene fusions and deletions, and a high level of homology between genes, its interrogation at high resolution has been thwarted by bioinformatic challenges, with most studies limited to examining presence or absence of specific genes. Here, we present the PING (Pushing Immunogenetics to the Next Generation) pipeline, which incorporates empirical data, novel alignment strategies and a custom alignment processing workflow to enable high-throughput KIR sequence analysis from short-read data. PING provides KIR gene copy number classification functionality for all KIR genes through use of a comprehensive alignment reference. The gene copy number determined per individual enables an innovative genotype determination workflow using genotype-matched references. Together, these methods address the challenges imposed by the structural complexity and overall homology of the KIR complex. To determine copy number and genotype determination accuracy, we applied PING to European and African validation cohorts and a synthetic dataset. PING demonstrated exceptional copy number determination performance across all datasets and robust genotype determination performance. Finally, an investigation into discordant genotypes for the synthetic dataset provides insight into misaligned reads, advancing our understanding in interpretation of short-read sequencing data in complex genomic regions. PING promises to support a new era of studies of KIR polymorphism, delivering high-resolution KIR genotypes that are highly accurate, enabling high-quality, high-throughput KIR genotyping for disease and population studies.

Noninvasive Determination of CMV Serostatus From Dried Buccal Swab Samples: Assay Development, Validation, and Application to 1.2 Million Samples.

BACKGROUND: Buccal swab sampling constitutes an attractive noninvasive alternative to blood drawings for antibody serostatus assays. Here we describe a method to determine the cytomegalovirus immunoglobulin G (CMV IgG) serostatus from dried buccal swab samples. METHODS: Upon solubilization, CMV IgG is determined by an ELISA assay specifically adapted to cope with low IgG concentrations. The derived CMV titer is normalized against the total protein concentration to adjust for incorrectly or less efficiently sampled buccal swabs. Assay parameters were optimized on a set of 713 samples. RESULTS: Validation with 1784 samples revealed distinct results for > 80% of samples with 98.6% specificity and 99.1% sensitivity. Based on the analysis of 1.2 million samples we derived age- and sex-stratified CMV prevalence statistics for Germany, Poland, United Kingdom, and Chile. To confirm accuracy of the assay in routine operation, the CMV status of 6518 donors was reassessed by independent laboratories based on conventional blood samples revealing 96.9% specificity and 97.4% sensitivity. CONCLUSIONS: The assay accurately delivers the CMV IgG serostatus from dried buccal swab samples for > 80% of the participants. Thereby it provides a noninvasive alternative to plasma-based CMV monitoring for nondiagnostic purposes such as hematopoietic stem cell transplantation donor screening or population studies.

Estimating HLA haplotype frequencies from homozygous individuals - A Technical Report.

We estimated HLA haplotype frequencies based on individuals homozygous for 4, 5 or 6 loci. Validation of our approach using a sample of over 3.4 million German individuals was successful. Compared to an expectation-maximization algorithm, the errors were larger. However, our approach allows the unequivocal detection of rare haplotypes.

Marriage does not relate to major histocompatibility complex: a genetic analysis based on 3691 couples.

Optimization of chances for healthy offspring is thought to be one of the factors driving mate choice and compatibility of the major histocompatibility complex (MHC) is assumed to determine the offspring's fitness. While humans have been claimed to be able to perceive information of MHC compatibility via the olfactory channel, it remains unknown whether humans use such information for mate choice. By investigation of 3691 married couples, we observed that the high polymorphism of MHC leads to a low chance for homozygous offspring. MHC similarity between couples did not differ from chance, we hence observed no MHC effect in married couples. Hormonal contraception at the time of relationship initiation had no significant effect towards enhanced similarity. A low variety of alleles within a postcode area led to a higher likelihood of homozygous offspring. Based on this data, we conclude that there is no pattern of MHC dis-assortative mating in a genetically diverse Western society. We discuss the question of olfactory mate preference, in-group mating bias and the high polymorphism as potential explanations.

Immunogenetics in stem cell donor registry work: The DKMS example (Part 1).

Currently, stem cell donor registries include more than 35 million potential donors worldwide to provide HLA-matched stem cell products for patients in need of an unrelated donor transplant. DKMS is a leading stem cell donor registry with more than 9 million donors from Germany, Poland, the United States, the United Kingdom, India and Chile. DKMS donors have donated hematopoietic stem cells more than 80,000 times. Many aspects of donor registry work are closely related to topics from immunogenetics or population genetics. In this two-part review article, we describe, analyse and discuss these areas of donor registry work by using the example of DKMS. Part 1 of the review gives a general overview on DKMS and includes typical donor registry activities with special focus on the HLA system: high-throughput HLA typing of potential stem cell donors, HLA haplotype frequencies and resulting matching probabilities, and donor file optimization with regard to HLA diversity.

Immunogenetics in stem cell donor registry work: The DKMS example (Part 2).

DKMS is a leading stem cell donor registry with more than 9 million donors. Donor registry activities share many touch points with topics from immunogenetics or population genetics. In this two-part review article, we deal with these aspects of donor registry work by using the example of DKMS. In the second part of the review, we focus on donor typing of non-HLA genes, the impact of donor age, gender and CMV serostatus on donation probabilities, the identification of novel HLA, KIR and MIC alleles by high-throughput donor typing, the activities of the Collaborative Biobank and pharmacogenetics in the donor registry context.

Correction: HLA-C*07 in axial spondyloarthritis: data from the German Spondyloarthritis Inception Cohort and the Spondyloarthritis Caught Early cohort.

The original version of this Article contained an error in the spelling of the author Denis Poddubnyy, which was incorrectly given as Denis Podubbnyy. This has now been corrected in both the PDF and HTML versions of the Article.

HLA-C*07 in axial spondyloarthritis: data from the German Spondyloarthritis Inception Cohort and the Spondyloarthritis Caught Early cohort.

We aimed to assess the mRNA expression of MHC class 1-related molecules in ankylosing spondylitis (AS) patients vs healthy controls (HCs) and, subsequently, if the absence of HLA-C*07 is associated with genetic susceptibility to axial spondyloarthritis (axSpA). HLA-C*07 was assessed in (a) an exploratory cohort of 24 AS patients vs 40 HCs, (b) a confirmatory cohort of 113 AS patients and 83 non-radiographic axSpA patients from the GErman SPondyloarthritis Inception Cohort (GESPIC) vs 134,528 German potential stem cell donors, and (c) an early back pain cohort with 94 early axSpA patients vs 216 chronic back pain (CBP) patients from the SPondyloArthritis Caught Early (SPACE) cohort. In the exploratory cohort, 79% of the AS patients were HLA-C*07 negative compared to 35% of the HCs (p < 0.001). This difference was confirmed in GESPIC with 73% of AS patients being HLA-C*07 negative compared to 50% of the controls (p < 0.0001); 59% of the nr-axSpA patients were HLA-C*07 negative. In the SPACE cohort, 70% of the axSpA patients were HLA-C*07 negative compared to 44% of CBP patients (p < 0.0001); the association between HLA-C*07 negativity and a diagnosis of axSpA was independent from HLA-B*27. In conclusion, the absence of HLA-C*07 is associated with genetic susceptibility to axSpA.

Human Leukocyte Antigen similarity decreases partners' and strangers' body odor attractiveness for women not using hormonal contraception.

The Human Leukocyte Antigen (HLA) is a gene complex that encodes important elements of the human immune system. HLA profile is communicated via olfaction and interindividual diversity is assumed to be advantageous for mate choice. Additionally, HLA diversity appears to enhance satisfaction and sexual attraction in existing romantic partnerships. However, whether this effect is transmitted via body odors and whether it results in an attraction towards HLA-dissimilar individuals and/or an avoidance of HLA-similar ones remains unclear. In the present study, we genotyped couples and asked each participant to rate a body odor sample from their partner and from three strangers of the opposite sex who expressed a similar or dissimilar HLA-B and HLA-C genotype. We found no statistically significant preference for HLA similarity or dissimilarity in men. Among women, the observed effects differed depending on hormonal contraception status. Like men, women on hormonal contraception did not exhibit significant HLA-related preferences. However, for women not using hormonal contraceptives, odors from HLA-B and HLA-C similar donors were significantly less attractive than those from HLA-dissimilar donors, regardless of whether the samples were from a partner or a stranger. Our findings support the hypothesis that HLA similarity is perceived via body odors and that such similarity affects human attraction. This mechanism may serve an evolutionarily adaptive function in preventing prospective offspring from having decreased immunocompetence, or decreasing the chance of kin mating.

Hapl-o-Mat: open-source software for HLA haplotype frequency estimation from ambiguous and heterogeneous data.

BACKGROUND: Knowledge of HLA haplotypes is helpful in many settings as disease association studies, population genetics, or hematopoietic stem cell transplantation. Regarding the recruitment of unrelated hematopoietic stem cell donors, HLA haplotype frequencies of specific populations are used to optimize both donor searches for individual patients and strategic donor registry planning. However, the estimation of haplotype frequencies from HLA genotyping data is challenged by the large amount of genotype data, the complex HLA nomenclature, and the heterogeneous and ambiguous nature of typing records. RESULTS: To meet these challenges, we have developed the open-source software Hapl-o-Mat. It estimates haplotype frequencies from population data including an arbitrary number of loci using an expectation-maximization algorithm. Its key features are the processing of different HLA typing resolutions within a given population sample and the handling of ambiguities recorded via multiple allele codes or genotype list strings. Implemented in C++, Hapl-o-Mat facilitates efficient haplotype frequency estimation from large amounts of genotype data. We demonstrate its accuracy and performance on the basis of artificial and real genotype data. CONCLUSIONS: Hapl-o-Mat is a versatile and efficient software for HLA haplotype frequency estimation. Its capability of processing various forms of HLA genotype data allows for a straightforward haplotype frequency estimation from typing records usually found in stem cell donor registries.

2.7 million samples genotyped for HLA by next generation sequencing: lessons learned.

BACKGROUND: At the DKMS Life Science Lab, Next Generation Sequencing (NGS) has been used for ultra-high-volume high-resolution genotyping of HLA loci for the last three and a half years. Here, we report on our experiences in genotyping the HLA, CCR5, ABO, RHD and KIR genes using a direct amplicon sequencing approach on Illumina MiSeq and HiSeq 2500 instruments. RESULTS: Between January 2013 and June 2016, 2,714,110 samples largely from German, Polish and UK-based potential stem cell donors have been processed. 98.9% of all alleles for the targeted HLA loci (HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1) were typed at high resolution or better. Initially a simple three-step workflow based on nanofluidic chips in conjunction with 4-primer amplicon tagging was used. Over time, we found that this setup results in PCR artefacts such as primer dimers and PCR-mediated recombination, which may necessitate repeat typing. Split workflows for low- and high-DNA-concentration samples helped alleviate these problems and reduced average per-locus repeat rates from 3.1 to 1.3%. Further optimisations of the workflow included the use of phosphorothioate oligos to reduce primer degradation and primer dimer formation, and employing statistical models to predict read yield from initial template DNA concentration to avoid intermediate quantification of PCR products. Finally, despite the populations typed at DKMS Life Science Lab being relatively homogenous genetically, an analysis of 1.4 million donors processed between January 2015 and May 2016 led to the discovery of 1,919 distinct novel HLA alleles. CONCLUSIONS: Amplicon-based NGS HLA genotyping workflows have become the workhorse in high-volume tissue typing of registry donors. The optimisation of workflow practices over multiple years has led to insights and solutions that improve the efficiency and robustness of short amplicon based genotyping workflows.

ABO allele-level frequency estimation based on population-scale genotyping by next generation sequencing.

BACKGROUND: The characterization of the ABO blood group status is vital for blood transfusion and solid organ transplantation. Several methods for the molecular characterization of the ABO gene, which encodes the alleles that give rise to the different ABO blood groups, have been described. However, the application of those methods has so far been restricted to selected samples and not been applied to population-scale analysis. RESULTS: We describe a cost-effective method for high-throughput genotyping of the ABO system by next generation sequencing. Sample specific barcodes and sequencing adaptors are introduced during PCR, rendering the products suitable for direct sequencing on Illumina MiSeq or HiSeq instruments. Complete sequence coverage of exons 6 and 7 enables molecular discrimination of the ABO subgroups and many alleles. The workflow was applied to ABO genotype more than a million samples. We report the allele group frequencies calculated on a subset of more than 110,000 sampled individuals of German origin. Further we discuss the potential of the workflow for high resolution genotyping taking the observed allele group frequencies into account. Finally, sequence analysis revealed 287 distinct so far not described alleles of which the most abundant one was identified in 174 samples. CONCLUSIONS: The described workflow delivers high resolution ABO genotyping at low cost enabling population-scale molecular ABO characterization.

Cost-efficient high-throughput HLA typing by MiSeq amplicon sequencing.

BACKGROUND: A close match of the HLA alleles between donor and recipient is an important prerequisite for successful unrelated hematopoietic stem cell transplantation. To increase the chances of finding an unrelated donor, registries recruit many hundred thousands of volunteers each year. Many registries with limited resources have had to find a trade-off between cost and resolution and extent of typing for newly recruited donors in the past. Therefore, we have taken advantage of recent improvements in NGS to develop a workflow for low-cost, high-resolution HLA typing. RESULTS: We have established a straightforward three-step workflow for high-throughput HLA typing: Exons 2 and 3 of HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1 are amplified by PCR on Fluidigm Access Array microfluidic chips. Illumina sequencing adapters and sample specific tags are directly incorporated during PCR. Upon pooling and cleanup, 384 samples are sequenced in a single Illumina MiSeq run. We developed "neXtype" for streamlined data analysis and HLA allele assignment. The workflow was validated with 1140 samples typed at 6 loci. All neXtype results were concordant with the Sanger sequences, demonstrating error-free typing of more than 6000 HLA loci. Current capacity in routine operation is 12,000 samples per week. CONCLUSIONS: The workflow presented proved to be a cost-efficient alternative to Sanger sequencing for high-throughput HLA typing. Despite the focus on cost efficiency, resolution exceeds the current standards of Sanger typing for donor registration.