Submitting Novel Full-Length HLA, MIC, and KIR Alleles with TypeLoader2.

The Immuno Polymorphism Database (IPD) plays a pivotal role for immunogenetics. Due to technical limitations, genotyping often focuses on specific key regions like the antigen recognition domain (ARD) for HLA genotyping, and the databases are populated accordingly. More recently, though, modern next generation sequencing (NGS) assays allow using larger gene segments or even complete genes for genotyping. It is therefore essential that the databases are updated with complete genetic reference sequences to fully serve current and future applications. However, the process of manually annotating and submitting full-length allele sequences to IPD is time-consuming and error-prone, which may discourage HLA-genotyping laboratories or researchers from submitting full-length sequences of novel alleles.Here, we detail the process of preparing and submitting novel HLA, MIC, and KIR alleles to ENA and IPD using TypeLoader2, a convenient software tool developed to streamline this process by automating the sequence annotation, the creation of all necessary files, as well as parts of the submission process itself. The software is freely available from GitHub ( https://github.com/DKMS-LSL/typeloader ).

Full-Length Characterization of Novel HLA-DRB1 Alleles for Reference Database Submission.

The prerequisite for successful HLA genotyping is the integrity of the large allele reference database IPD-IMGT/HLA. Consequently, it is in the laboratories' best interest that the data quality of submitted novel sequences is high. However, due to its long and variable length, the gene HLA-DRB1 presents the biggest challenge and as of today only 16% of the HLA-DRB1 alleles in the database are characterized in full length. To improve this situation, we developed a protocol for long-range PCR amplification of targeted HLA-DRB1 alleles. By subsequently combining both long-read and short-read sequencing technologies, our protocol ensures phased and error-corrected sequences of reference grade quality. This dual redundant reference sequencing (DR2S) approach is of particular importance for correctly resolving the challenging repeat regions of DRB1 intron 1. Until today, we used this protocol to characterize and submit 384 full-length HLA-DRB1 sequences to IPD-IMGT/HLA.

High population frequencies of MICA copy number variations originate from independent recombination events.

MICA is a stress-induced ligand of the NKG2D receptor that stimulates NK and T cell responses and was identified as a key determinant of anti-tumor immunity. The MICA gene is located inside the MHC complex and is in strong linkage disequilibrium with HLA-B. While an HLA-B*48-linked MICA deletion-haplotype was previously described in Asian populations, little is known about other MICA copy number variations. Here, we report the genotyping of more than two million individuals revealing high frequencies of MICA duplications (1%) and MICA deletions (0.4%). Their prevalence differs between ethnic groups and can rise to 2.8% (Croatia) and 9.2% (Mexico), respectively. Targeted sequencing of more than 70 samples indicates that these copy number variations originate from independent nonallelic homologous recombination events between segmental duplications upstream of MICA and MICB. Overall, our data warrant further investigation of disease associations and consideration of MICA copy number data in oncological study protocols.

DR2S: an integrated algorithm providing reference-grade haplotype sequences from heterozygous samples.

BACKGROUND: High resolution HLA genotyping of donors and recipients is a crucially important prerequisite for haematopoetic stem-cell transplantation and relies heavily on the quality and completeness of immunogenetic reference sequence databases of allelic variation. RESULTS: Here, we report on DR2S, an R package that leverages the strengths of two sequencing technologies-the accuracy of next-generation sequencing with the read length of third-generation sequencing technologies like PacBio's SMRT sequencing or ONT's nanopore sequencing-to reconstruct fully-phased high-quality full-length haplotype sequences. Although optimised for HLA and KIR genes, DR2S is applicable to all loci with known reference sequences provided that full-length sequencing data is available for analysis. In addition, DR2S integrates supporting tools for easy visualisation and quality control of the reconstructed haplotype to ensure suitability for submission to public allele databases. CONCLUSIONS: DR2S is a largely automated workflow designed to create high-quality fully-phased reference allele sequences for highly polymorphic gene regions such as HLA or KIR. It has been used by biologists to successfully characterise and submit more than 500 HLA alleles and more than 500 KIR alleles to the IPD-IMGT/HLA and IPD-KIR databases.

HLA-E diversity unfolded: Identification and characterization of 170 novel HLA-E alleles.

HLA-E is a member of the nonclassical HLA class Ib genes. Even though it is structurally highly similar to the classical HLA class Ia genes, it is less diverse and only 45 alleles and 12 proteins were known in December 2019 (IPD-IMGT/HLA, release 3.38.0). Since 2017, we have genotyped over 3 million voluntary stem cell donors for HLA-E by sequencing the most relevant allele-determining bases of exons 2 and 3. As expected, most donors harbor the two predominant alleles HLA-E*01:01 and/or HLA-E*01:03. However, in 1666 (0.05%) of our samples we detected 345 distinct novel HLA-E sequences. The most frequent one was identified in 162 samples and has by now been named HLA-E*01:114. To characterize these novel alleles in full-length, we used both short-read Illumina and long-read PacBio sequencing to obtain fully phased and highly accurate sequences. This resulted in 234 submissions to IPD-IMGT/HLA comprising 170 novel HLA-E alleles, which encode for 93 novel HLA-E proteins, as well as 64 confirmations or sequence extensions. Consequently, the number of HLA-E alleles in the database (release 3.42.0) has now increased to 256 HLA-E alleles and 110 HLA-E proteins.

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.

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.

Bioinformatics Strategies, Challenges, and Opportunities for Next Generation Sequencing-Based HLA Genotyping.

The advent of next generation sequencing (NGS) has altered the face of genotyping the human leukocyte antigen (HLA) system in clinical, stem cell donor registry, and research contexts. NGS has led to a dramatically increased sequencing throughput at high accuracy, while being more time and cost efficient than precursor technologies. This has led to a broader and deeper profiling of the key genes in the human immunogenetic make-up. The rapid evolution of sequencing technologies is evidenced by the development of varied short-read sequencing platforms with differing read lengths and sequencing capacities to long-read sequencing platforms capable of profiling full genes without fragmentation. Concomitantly, there has been development of a diverse set of computational analyses and software tools developed to deal with the various strengths and limitations of the sequencing data generated by the different sequencing platforms. This review surveys the different modalities involved in generating NGS HLA profiling sequence data. It systematically describes various computational approaches that have been developed to achieve HLA genotyping to different degrees of resolution. At each stage, this review enumerates the drawbacks and advantages of each of the platforms and analysis approaches, thus providing a comprehensive picture of the current state of HLA genotyping technologies.

Within-population variability in a moth sex pheromone blend, part 2: selection towards fixation.

To understand how variation in sexual communication systems evolves, the genetic architecture underlying sexual signals and responses needs to be identified. Especially in animals where mating signals are important for mate recognition, and signals and responses are governed by independently assorting genes, it is difficult to envision how signals and preferences can (co)evolve. Moths are a prime example of such animals. In the noctuid moth Heliothis virescens, we found within-population variation in the female pheromone. In previous selection experiments followed by quantitative trait locus (QTL) analysis and expression analysis of candidate desaturase genes, we developed a model involving a trans-acting repressor of the delta-11-desaturase. In our current study with new selection lines, we fixed the most extreme phenotype and found a single underlying mutation: a premature stop codon in the first coding exon of delta-11-desaturase, which we could trace back to its origin in the laboratory. Interestingly, we found no pleiotropic effects of this knock-out mutation on the male physiological or behavioural response, or on growth or fertility. This finding is in contrast to Drosophila melanogaster, where a single desaturase gene affects both female pheromone production and male behavioural response, but similar to other Lepidoptera where these traits are under independent genetic control. To our knowledge, this is the first time that a single point mutation has been identified that underlies the phenotypic variation in the pheromone signal of a moth.

TypeLoader2: Automated submission of novel HLA and killer-cell immunoglobulin-like receptor alleles in full length.

The Immuno Polymorphism Database (IPD) databases provide global, curated repositories for information regarding polymorphisms of genes of the immune system, thereby generating immense value for the research and clinical communities. The advent of high-throughput genotyping in immunogenetics has led to dramatically growing numbers of heretofore unknown HLA and lately also killer-cell immunoglobulin-like receptor (KIR) alleles, which are to be curated and deposited in the IPD-IMGT/HLA and IPD-KIR databases, respectively. It is highly desirable that these novel alleles are characterised and submitted in full length, and that known alleles are extended to cover the complete gene sequence. However, the manual annotation and submission of sequences to European Molecular Biology Laboratory's European Nucleotide Archive and the IPD-IMGT/HLA and IPD-KIR databases is time-consuming and error-prone. Here, we report the substantial extension of the HLA allele submission tool TypeLoader, which now also supports the annotation and submission of KIR alleles. To enable a more widespread use of this tool, we have made it available as a stand-alone application that can easily be installed on standard Windows or Linux computers. Furthermore, an internal SQLite database was added to store a wide range of metadata about each allele. This allows TypeLoader2 to be used as a lab's central information platform for the annotation, curation and submission of full-length HLA and KIR allele sequences. The software is freely available from GitHub (https://github.com/DKMS-LSL/typeloader). We hope that the increased convenience and scope of TypeLoader2 will foster the submission of more full-length sequences to the IPD-IMGT/HLA and IPD-KIR databases, ultimately promoting the use of full-length sequencing for genotyping both HLA and KIR.

Patterns of non-ARD variation in more than 300 full-length HLA-DPB1 alleles.

Our understanding of sequence variation in the HLA-DPB1 gene is largely restricted to the hypervariable antigen recognition domain (ARD) encoded by exon 2. Here, we employed a redundant sequencing strategy combining long-read and short-read data to accurately phase and characterise in full length the majority of common and well-documented (CWD) DPB1 alleles as well as alleles with an observed frequency of at least 0.0006% in our predominantly European sample set. We generated 664 DPB1 sequences, comprising 279 distinct allelic variants. This allows us to present the, to date, most comprehensive analysis of the nature and extent of DPB1 sequence variation. The full-length sequence analysis revealed the existence of two highly diverged allele clades. These clades correlate with the rs9277534 A → G variant, a known expression marker located in the 3'-UTR. The two clades are fully differentiated by 174 fixed polymorphisms throughout a 3.6 kb stretch at the 3'-end of DPB1. The region upstream of this differentiation zone is characterised by increasingly shared variation between the clades. The low-expression A clade comprises 59% of the distinct allelic sequences including the three by far most frequent DPB1 alleles, DPB1*04:01, DPB1*02:01 and DPB1*04:02. Alleles in the A clade show reduced nucleotide diversity with an excess of rare variants when compared to the high-expression G clade. This pattern is consistent with a scenario of recent proliferation of A-clade alleles. The full-length characterisation of all but the most rare DPB1 alleles will benefit the application of NGS for DPB1 genotyping and provides a helpful framework for a deeper understanding of high- and low-expression alleles and their implications in the context of unrelated haematopoietic stem-cell transplantation.

Allele-Level KIR Genotyping of More Than a Million Samples: Workflow, Algorithm, and Observations.

The killer-cell immunoglobulin-like receptor (KIR) genes regulate natural killer cell activity, influencing predisposition to immune mediated disease, and affecting hematopoietic stem cell transplantation (HSCT) outcome. Owing to the complexity of the KIR locus, with extensive gene copy number variation (CNV) and allelic diversity, high-resolution characterization of KIR has so far been applied only to relatively small cohorts. Here, we present a comprehensive high-throughput KIR genotyping approach based on next generation sequencing. Through PCR amplification of specific exons, our approach delivers both copy numbers of the individual genes and allelic information for every KIR gene. Ten-fold replicate analysis of a set of 190 samples revealed a precision of 99.9%. Genotyping of an independent set of 360 samples resulted in an accuracy of more than 99% taking into account consistent copy number prediction. We applied the workflow to genotype 1.8 million stem cell donor registry samples. We report on the observed KIR allele diversity and relative abundance of alleles based on a subset of more than 300,000 samples. Furthermore, we identified more than 2,000 previously unreported KIR variants repeatedly in independent samples, underscoring the large diversity of the KIR region that awaits discovery. This cost-efficient high-resolution KIR genotyping approach is now applied to samples of volunteers registering as potential donors for HSCT. This will facilitate the utilization of KIR as additional selection criterion to improve unrelated donor stem cell transplantation outcome. In addition, the approach may serve studies requiring high-resolution KIR genotyping, like population genetics and disease association studies.

Full-Length HLA Class I Genotyping with the MinION Nanopore Sequencer.

Nanopore sequencing, a paradigm change in sequencing technologies, offers a new cost-effective and scalable platform for HLA genotyping. Among the new generation of high-throughput sequencing technologies, the MinION nanopore sequencer is the first to offer a non-template-based direct DNA sensing sequencing technology. Oxford Nanopore Technologies (ONT) introduced the first version of the MinION in 2014; since then, the platform has gone through multiple iterations resulting in higher throughput and sequencing accuracy. The "what you put in is what you get" nature of the platform enables molecules to be sequenced without fragmentation. This results in ultra-long read lengths in the order of tens of kilobases enabling entire genes to be characterized with fully phased sequence information. With release R9.5, the MinION platform has reached a quality that enables HLA genotyping with minor shortcomings in long homopolymer regions. Within this chapter, we describe a protocol for sequencing and genotyping HLA Class I alleles using the MinION.

Genetic basis of allochronic differentiation in the fall armyworm.

BACKGROUND: Very little is known on how changes in circadian rhythms evolve. The noctuid moth Spodoptera frugiperda (Lepidoptera: Noctuidae) consists of two strains that exhibit allochronic differentiation in their mating time, which acts as a premating isolation barrier between the strains. We investigated the genetic basis of the strain-specific timing differences to identify the molecular mechanisms of differentiation in circadian rhythms. RESULTS: Through QTL analyses we identified one major Quantitative trait chromosome (QTC) underlying differentiation in circadian timing of mating activity. Using RADtags, we identified this QTC to be homologous to Bombyx mori C27, on which the clock gene vrille is located, which thus became the major candidate gene. In S. frugiperda, vrille showed strain-specific polymorphisms. Also, vrille expression differed significantly between the strains, with the rice-strain showing higher expression levels than the corn-strain. In addition, RT-qPCR experiments with the other main clock genes showed that pdp1, antagonist of vrille in the modulatory feedback loop of the circadian clock, showed higher expression levels in the rice-strain than in the corn-strain. CONCLUSIONS: Together, our results indicate that the allochronic differentiation in the two strains of S. frugiperda is associated with differential transcription of vrille or a cis-acting gene close to vrille, which contributes to the evolution of prezygotic isolation in S. frugiperda.

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.

Prediction of spurious HLA class II typing results using probabilistic classification.

While modern high-throughput sequence-based HLA genotyping methods generally provide highly accurate typing results, artefacts may nonetheless arise for numerous reasons, such as sample contamination, sequencing errors, read misalignments, or PCR amplification biases. To help detecting spurious typing results, we tested the performance of two probabilistic classifiers (binary logistic regression and random forest models) based on population-specific genotype frequencies. We trained the model using high-resolution typing results for HLA-DRB1, DQB1, and DPB1 from large samples of German, Polish and UK-based donors. The high predictive capacity of the best models replicated both in 10-fold cross-validation for each gene and in using independent evaluation data (AUC 0.820-0.893). While genotype frequencies alone provide enough predictive power to render the model generally useful for highlighting potentially spurious typing results, the inclusion of workflow-specific predictors substantially increases prediction specificity. Low initial DNA concentrations in combination with low-volume PCR reactions form a major source of stochastic error specific to the Fluidigm chip-based workflow at DKMS Life Science Lab. The addition of DNA concentrations as a predictor variable thus substantially increased AUC (0.947-0.959) over purely frequency-based models.

Population structure of the Indonesian giant tiger shrimp Penaeus monodon: a window into evolutionary similarities between paralogous mitochondrial DNA sequences and their genomes.

Here we used both microsatellites and mtCR (mitochondrial DNA control region) sequences as genetic markers to examine the genetic diversity and population structure of Penaeus monodon shrimp from six Indonesian regions. The microsatellite data showed that shrimp from the Indian and the Pacific Ocean were genetically distinct from each other. It has been reported previously that P. monodon mtCR sequences from the Indo-Pacific group into two major paralogous clades of unclear origin. Here we show that the population structure inferred from mtCR sequences matches the microsatellite-based population structure for one of these clades. This is consistent with the notion that this mtCR clade shares evolutionary history with nuclear DNA and may thus represent nuclear mitochondrial pseudogenes (Numts).

Expression of the toxin-antitoxin genes yefM(Lrh), yoeB(Lrh) in human Lactobacillus rhamnosus isolates.

Lactobacilli are important microorganisms in various activities, for example, diary products, meat ripening, bread and pickles, but, moreover, are associated directly with human skin and cavities (e.g., mouth, gut, or vagina). Some of them are used as probiotics. Therefore, the molecular biological investigation of these bacteria is important. Earlier we described several toxin antitoxin systems (type II) in lactobacilli. Here, we describe the structure and transcriptional regulation of genes, encoding TA system YefM-YoeB(Lrh) in three strains of Lactobacillus rhamnosus comparing stationary and exponential growth phases, the influence of stress factors and mRNA stability. The same TA system is responding to physiological and stress conditions differently in related strains. Using primer extension and RLM-RACE methods we determined three transcription start sites of RNAs in the operon. The promoter region of the operon is preceded by a conserved BOX element occurring at multiple positions in the genomes of L. rhamnosus strains. Downstream of and partially overlapping with the 3' end of the yoeB(Lrh) toxin gene, a divergently transcribed unexpected RNA was detected.

Chlamydia psittaci: new insights into genomic diversity, clinical pathology, host-pathogen interaction and anti-bacterial immunity.

The distinctive and unique features of the avian and mammalian zoonotic pathogen Chlamydia (C.) psittaci include the fulminant course of clinical disease, the remarkably wide host range and the high proportion of latent infections that are not leading to overt disease. Current knowledge on associated diseases is rather poor, even in comparison to other chlamydial agents. In the present paper, we explain and summarize the major findings of a national research network that focused on the elucidation of host-pathogen interactions in vitro and in animal models of C. psittaci infection, with the objective of improving our understanding of genomics, pathology, pathophysiology, molecular pathogenesis and immunology, and conceiving new approaches to therapy. We discuss new findings on comparative genome analysis, the complexity of pathophysiological interactions and systemic consequences, local immune response, the role of the complement system and antigen presentation pathways in the general context of state-of-the-art knowledge on chlamydial infections in humans and animals and single out relevant research topics to fill remaining knowledge gaps on this important yet somewhat neglected pathogen.