Rapid Characterization of Complex Killer Cell Immunoglobulin-Like Receptor (KIR) Regions Using Cas9 Enrichment and Nanopore Sequencing

Complete Genomic Assembly of Mauritian Cynomolgus Macaque Killer Ig-like Receptor and Natural Killer Group 2 Haplotypes

Mauritian-origin cynomolgus macaques (MCMs) serve as a powerful nonhuman primate model in biomedical research due to their unique genetic homogeneity, which simplifies experimental designs. Despite their extensive use, a comprehensive understanding of crucial immune-regulating gene families, particularly killer Ig-like receptors (KIR) and NK group 2 (NKG2), has been hindered by the lack of detailed genomic reference assemblies. In this study, we employ advanced long-read sequencing techniques to completely assemble eight KIR and seven NKG2 genomic haplotypes, providing an extensive insight into the structural and allelic diversity of these immunoregulatory gene clusters. Leveraging these genomic resources, we prototype a strategy for genotyping KIR and NKG2 using short-read, whole-exome capture data, illustrating the potential for cost-effective multilocus genotyping at colony scale. These results mark a significant enhancement for biomedical research in MCMs and underscore the feasibility of broad-scale genetic investigations.

Advancements in HLA Typing Techniques and Their Impact on Transplantation Medicine

HLA typing serves as a standard practice in hematopoietic stem cell transplantation to ensure compatibility between donors and recipients, preventing the occurrence of allograft rejection and graft-versus-host disease. Conventional laboratory methods that have been widely employed in the past few years, including sequence-specific primer PCR and sequencing-based typing (SBT), currently face the risk of becoming obsolete. This risk stems not only from the extensive diversity within HLA genes but also from the rapid advancement of next-generation sequencing and third-generation sequencing technologies. Third-generation sequencing systems like single-molecule real-time (SMRT) sequencing and Oxford Nanopore (ONT) sequencing have the capability to analyze long-read sequences that span entire intronic-exonic regions of HLA genes, effectively addressing challenges related to HLA ambiguity and the phasing of multiple short-read fragments. The growing dominance of these advanced sequencers in HLA typing is expected to solidify further through ongoing refinements, cost reduction, and error rate minimization. This review focuses on hematopoietic stem cell transplantation (HSCT) and explores prospective advancements and application of HLA DNA typing techniques. It explores how the adoption of third-generation sequencing technologies can revolutionize the field by offering improved accuracy, reduced ambiguity, and enhanced assessment of compatibility in HSCT. Embracing these cutting-edge technologies is essential to advancing the success rates and outcomes of hematopoietic stem cell transplantation. This review underscores the importance of staying at the forefront of HLA typing techniques to ensure the best possible outcomes for patients undergoing HSCT.

The KIR repertoire of a West African chimpanzee population is characterized by limited gene, allele, and haplotype variation

Introduction

The killer cell immunoglobulin-like receptors (KIR) play a pivotal role in modulating the NK cell responses, for instance, through interaction with major histocompatibility complex (MHC) class I molecules. Both gene systems map to different chromosomes but co-evolved during evolution. The human KIR gene family is characterized by abundant allelic polymorphism and copy number variation. In contrast, our knowledge of the KIR repertoire in chimpanzees is limited to 39 reported alleles, with no available population data. Only three genomic KIR region configurations have been mapped, and seventeen additional ones were deduced by genotyping.

Methods

Previously, we documented that the chimpanzee MHC class I repertoire has been skewed due to an ancient selective sweep. To understand the depth of the sweep, we set out to determine the full-length KIR transcriptome - in our MHC characterized pedigreed West African chimpanzee cohort - using SMRT sequencing (PacBio). In addition, the genomic organization of 14 KIR haplotypes was characterized by applying a Cas9-mediated enrichment approach in concert with long-read sequencing by Oxford Nanopore Technologies.

Results

In the cohort, we discovered 35 undescribed and 15 already recorded Patr-KIR alleles, and a novel hybrid KIR gene. Some KIR transcripts are subject to evolutionary conserved alternative splicing events. A detailed insight on the KIR region dynamics (location and order of genes) was obtained, however, only five new KIR region configurations were detected. The population data allowed to investigate the distribution of the MHC-C1 and C2-epitope specificity of the inhibitory lineage III KIR repertoire, and appears to be skewed towards C2.

Discussion

Although the KIR region is known to evolve fast, as observed in other primate species, our overall conclusion is that the genomic architecture and repertoire in West African chimpanzees exhibit only limited to moderate levels of variation. Hence, the ancient selective sweep that affected the chimpanzee MHC class I region may also have impacted the KIR system.

Beyond assembly: the increasing flexibility of single-molecule sequencing technology

The maturation of high-throughput short-read sequencing technology over the past two decades has shaped the way genomes are studied. Recently, single-molecule, long-read sequencing has emerged as an essential tool in deciphering genome structure and function, including filling gaps in the human reference genome, measuring the epigenome and characterizing splicing variants in the transcriptome. With recent technological developments, these single-molecule technologies have moved beyond genome assembly and are being used in a variety of ways, including to selectively sequence specific loci with long reads, measure chromatin state and protein-DNA binding in order to investigate the dynamics of gene regulation, and rapidly determine copy number variation. These increasingly flexible uses of single-molecule technologies highlight a young and fast-moving part of the field that is leading to a more accessible era of nucleic acid sequencing.

HLA/MHC and KIR characterization in humans and non-human primates using Oxford Nanopore Technologies and Pacific Biosciences sequencing platforms

The gene products of the HLA/MHC and KIR multigene families are important modulators of the immune system and are associated with health and disease. Characterization of the genes encoding these receptors has been integrated into different biomedical applications, including transplantation and reproduction biology, immune therapies and in fundamental research into disease susceptibility or resistance. Conventional short-read sequencing strategies have shown their value in high throughput typing, but are insufficient to uncover the entire complexity of the highly polymorphic HLA/MHC and KIR gene systems. The implementation of single-molecule and real-time sequencing platforms, offered by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), revolutionized the fields of genomics and transcriptomics. Using fundamentally distinct principles, these platforms generate long-read data that can unwire the plasticity of the HLA/MHC and KIR genes, including high-resolution characterization of genes, alleles, phased haplotypes, transcription levels and epigenetics modification patterns. These insights might have profound clinical relevance, such as improved matching of donors and patients in clinical transplantation, but could also lift disease association studies to a higher level. Even more, a comprehensive characterization may refine animal models in preclinical studies. In this review, the different HLA/MHC and KIR characterization approaches using PacBio and ONT platforms are described and discussed.

Complete sequencing of a cynomolgus macaque major histocompatibility complex haplotype

Macaques provide the most widely used nonhuman primate models for studying the immunology and pathogenesis of human diseases. Although the macaque major histocompatibility complex (MHC) region shares most features with the human leukocyte antigen (HLA) region, macaques have an expanded repertoire of MHC class I genes. Although a chimera of two rhesus macaque MHC haplotypes was first published in 2004, the structural diversity of MHC genomic organization in macaques remains poorly understood owing to a lack of adequate genomic reference sequences. We used ultralong Oxford Nanopore and high-accuracy Pacific Biosciences (PacBio) HiFi sequences to fully assemble the ∼5.2-Mb M3 haplotype of an MHC-homozygous, Mauritian-origin cynomolgus macaque (Macaca fascicularis). The MHC homozygosity allowed us to assemble a single MHC haplotype unambiguously and avoid chimeric assemblies that hampered previous efforts to characterize this exceptionally complex genomic region in macaques. The high quality of this new assembly is exemplified by the identification of an extended cluster of six Mafa-AG genes that contains a recent duplication with a highly similar ∼48.5-kb block of sequence. The MHC class II region of this M3 haplotype is similar to the previously sequenced rhesus macaque haplotype and HLA class II haplotypes. The MHC class I region, in contrast, contains 13 MHC-B genes, four MHC-A genes, and three MHC-E genes (vs. 19 MHC-B, two MHC-A, and one MHC-E in the previously sequenced haplotype). These results provide an unambiguously assembled single contiguous cynomolgus macaque MHC haplotype with fully curated gene annotations that will inform infectious disease and transplantation research.

CRISPR-Cas9 Targeted Enrichment and Next-Generation Sequencing for Mutation Detection

Despite the rapid application of next-generation sequencing (NGS) technologies, target sequencing in regions of the genome is often required to diagnose many genetic diseases. Target enrichment can be an effective factor in reducing the cost of sequencing and the duration of sequencing. Recently, several clustered system regularly interspaced short palindromic repeats (CRISPR)-based methods (amplification-free sequencing) have been developed to target enrichment in combination with one of the NGS platforms. CRISPR-based target enrichment strategies act as an auxiliary tool to improve NGS analytical performance, thereby indirectly facilitating nucleic acid detection. The direct DNA cleavage approach by CRISPR-Cas at genome-specific sites enhances the possibility of separating native large fragments from disease-related genomic regions. The CRISPR-Cas can isolate the target region without any amplification; subsequently, long-read sequencing technologies were also implemented. These methods, as promising tools, have the ability to assess genetic and epigenetic composition for clinical application and treatment responses in cancer precision medicine. By modifying CRISPR-based enrichment protocols, it was possible to identify different types of mutations, including structural variants, short tandem repeats, fusion genes, and mobile elements. The Cas9 can specifically eliminate wild-type sequences, and it also enables the enrichment and detection of small amounts of tumor DNA fragments among the highly heterogeneous fragments of wild-type DNA.

CRISPR/Cas9-Mediated Enrichment Coupled to Nanopore Sequencing Provides a Valuable Tool for the Precise Reconstruction of Large Genomic Target Regions

Complete and accurate identification of genetic variants associated with specific phenotypes can be challenging when there is a high level of genomic divergence between individuals in a study and the corresponding reference genome. We have applied the Cas9-mediated enrichment coupled to nanopore sequencing to perform a targeted de novo assembly and accurately reconstruct a genomic region of interest. This approach was used to reconstruct a 250-kbp target region on chromosome 5 of the common bean genome (Phaseolus vulgaris) associated with the shattering phenotype. Comparing a non-shattering cultivar (Midas) with the reference genome revealed many single-nucleotide variants and structural variants in this region. We cut five 50-kbp tiled sub-regions of Midas genomic DNA using Cas9, followed by sequencing on a MinION device and de novo assembly, generating a single contig spanning the whole 250-kbp region. This assembly increased the number of Illumina reads mapping to genes in the region, improving their genotypability for downstream analysis. The Cas9 tiling approach for target enrichment and sequencing is a valuable alternative to whole-genome sequencing for the assembly of ultra-long regions of interest, improving the accuracy of downstream genotype-phenotype association analysis.

Long-read assembly of major histocompatibility complex and killer cell immunoglobulin-like receptor genome regions in cynomolgus macaque

BACKGROUND

The major histocompatibility complex (MHC) and the killer cell immunoglobulin-like receptors (KIR) are key regulators of immune responses. The cynomolgus macaque, an Old World monkey species, can be applied as an important preclinical model for studying human diseases, including coronavirus disease 2019 (COVID-19). Several MHC-KIR combinations have been associated with either a poor or good prognosis. Therefore, macaques with a well-characterized immunogenetic profile may improve drug evaluation and speed up vaccine development. At present, a complete overview of the MHC and KIR haplotype organizations in cynomolgus macaques is lacking, and characterization by conventional techniques is hampered by the extensive expansion of the macaque MHC-B region that complicates the discrimination between genes and alleles.

METHODS

We assembled complete MHC and KIR genomic regions of cynomolgus macaque using third-generation long-read sequencing approach. We identified functional Mafa-B loci at the transcriptome level using locus-specific amplification in a cohort of 33 Vietnamese cynomolgus macaques.

RESULTS

This is the first physical mapping of complete MHC and KIR gene regions in a Vietnamese cynomolgus macaque. Furthermore, we identified four functional Mafa-B loci (B2, B3, B5, and B6) and showed that alleles of the Mafa-I*01, -B*056, -B*034, and -B*001 functional lineages, respectively, are highly frequent in the Vietnamese cynomolgus macaque population.

CONCLUSION

The insights into the MHC and KIR haplotype organizations and the level of diversity may refine the selection of animals with specific genetic markers for future medical research.

An efficient CRISPR-Cas9 enrichment sequencing strategy for characterizing complex and highly duplicated genomic regions. A case study in the Prunus salicina LG3-MYB10 genes cluster

BACKGROUND

Genome complexity is largely linked to diversification and crop innovation. Examples of regions with duplicated genes with relevant roles in agricultural traits are found in many crops. In both duplicated and non-duplicated genes, much of the variability in agronomic traits is caused by large as well as small and middle scale structural variants (SVs), which highlights the relevance of the identification and characterization of complex variability between genomes for plant breeding.

RESULTS

Here we improve and demonstrate the use of CRISPR-Cas9 enrichment combined with long-read sequencing technology to resolve the MYB10 region in the linkage group 3 (LG3) of Japanese plum (Prunus salicina). This region, which has a length from 90 to 271 kb according to the P. salicina genomes available, is associated with fruit color variability in Prunus species. We demonstrate the high complexity of this region, with homology levels between Japanese plum varieties comparable to those between Prunus species. We cleaved MYB10 genes in five plum varieties using the Cas9 enzyme guided by a pool of crRNAs. The barcoded fragments were then pooled and sequenced in a single MinION Oxford Nanopore Technologies (ONT) run, yielding 194 Mb of sequence. The enrichment was confirmed by aligning the long reads to the plum reference genomes, with a mean read on-target value of 4.5% and a depth per sample of 11.9x. From the alignment, 3261 SNPs and 287 SVs were called and phased. A de novo assembly was constructed for each variety, which also allowed detection, at the haplotype level, of the variability in this region.

CONCLUSIONS

CRISPR-Cas9 enrichment is a versatile and powerful tool for long-read targeted sequencing even on highly duplicated and/or polymorphic genomic regions, being especially useful when a reference genome is not available. Potential uses of this methodology as well as its limitations are further discussed.

Dynamic evolution of Mhc haplotypes in cynomolgus macaques of different geographic origins

The major histocompatibility complex (MHC) plays a key role in immune defense, and the Mhc genes of cynomolgus macaque display a high degree of polymorphism. Based on their geographic distribution, different populations of cynomolgus macaques are recognized. Here we present the characterization of the Mhc class I and II repertoire of a large pedigreed group of cynomolgus macaques originating from the mainland north of the isthmus of Kra (N = 42). Segregation analyses resulted in the definition of 81 unreported Mafa-A/B/DRB/DQ/DP haplotypes, which include 32 previously unknown DRB regions. In addition, we report 13 newly defined Mafa-A/B/DRB/DQ/DP haplotypes in a group of cynomolgus macaques originating from the mainland south of the isthmus of Kra/Maritime Southeast Asia (N = 16). A relatively high level of sharing of Mafa-A (51%) and Mafa-B (40%) lineage groups is observed between the populations native to the north and the south of isthmus of Kra. At the allelic level, however, the Mafa-A/B haplotypes seem to be characteristic of a population. An overall comparison of all currently known data revealed that each geographic population has its own specific combinations of Mhc class I and II haplotypes. This illustrates the dynamic evolution of the cynomolgus macaque Mhc region, which was most likely generated by recombination and maintained by selection due to the differential pathogenic pressures encountered in different geographic areas.

Comparative genetics of KIR haplotype diversity in humans and rhesus macaques: the balancing act

The role of natural killer (NK) cells is tightly modulated by interactions of killer cell immunoglobulin-like receptors (KIR) with their ligands of the MHC class I family. Several characteristics of the KIR gene products are conserved in primate evolution, like the receptor structures and the variegated expression pattern. At the genomic level, however, the clusters encoding the KIR family display species-specific diversity, reflected by differential gene expansions and haplotype architecture. The human KIR cluster is extensively studied in large cohorts from various populations, which revealed two KIR haplotype groups, A and B, that represent more inhibitory and more activating functional profiles, respectively. So far, genomic KIR analyses in large outbred populations of non-human primate species are lacking. In this study, we roughly quadrupled the number of rhesus macaques studied for their KIR transcriptome (n = 298). Using segregation analysis, we defined 112 unique KIR region configurations, half of which display a more inhibitory profile, whereas the other half has a more activating potential. The frequencies and functional potential of these profiles might mirror the human KIR haplotype groups. However, whereas the human group A and B KIR haplotypes are confined to largely fixed organizations, the haplotypes in macaques feature highly variable gene content. Moreover, KIR homozygosity was hardly encountered in this panel of macaques. This study exhibits highly diverse haplotype architectures in humans and macaques, which nevertheless might have an equivalent effect on the modulation of NK cell activity.

High-resolution human KIR genotyping

Killer immunoglobulin-like receptors (KIR) regulate the function of natural killer cells through interactions with various ligands on the surface of cells, thereby determining whether natural killer (NK) cells are to be activated or inhibited from killing the cell being interrogated. The genes encoding these proteins display extensive variation through variable gene content, copy number and allele polymorphism. The combination of KIR genes and their ligands is implicated in various clinical settings including haematopoietic stem cell and solid organ transplant and infectious disease progression. The determination of KIR genes has been used as a factor in the selection of optimal stem cell donors with haplotype variations in recipient and donor giving differential clinical outcomes. Methods to determine KIR genes have primarily involved ascertaining the presence or absence of genes in an individual. With the more recent introduction of massively parallel clonal next-generation sequencing and single molecule very long read length third-generation sequencing, high-resolution determination of KIR alleles has become feasible. Determining the extent and functional impact of allele variation has the potential to lead to further optimisation of clinical outcomes as well as a deeper understanding of the functional properties of the receptors and their interactions with ligands. This review summarizes recently published high-resolution KIR genotyping methods and considers the various advantages and disadvantages of the approaches taken. In addition the application of allele level genotyping in the setting of transplantation and infectious disease control is discussed.