Immunogenomics of Killer Cell Immunoglobulin-Like Receptor (KIR) and HLA Class I: Coevolution and Consequences for Human Health.

Interactions of killer cell immunoglobin-like receptors (KIR) with human leukocyte antigens (HLA) class I regulate effector functions of key cytotoxic cells of innate and adaptive immunity. The extreme diversity of this interaction is genetically determined, having evolved in the ever-changing environment of pathogen exposure. Diversity of KIR and HLA genes is further facilitated by their independent segregation on separate chromosomes. That fetal implantation relies on many of the same types of immune cells as infection control places certain constraints on the evolution of KIR interactions with HLA. Consequently, specific inherited combinations of receptors and ligands may predispose to specific immune-mediated diseases, including autoimmunity. Combinatorial diversity of KIR and HLA class I can also differentiate success rates of immunotherapy directed to these diseases. Progress toward both etiopathology and predicting response to therapy is being achieved through detailed characterization of the extent and consequences of the combinatorial diversity of KIR and HLA. Achieving these goals is more tractable with the development of integrated analyses of molecular evolution, function, and pathology that will establish guidelines for understanding and managing risks. Here, we present what is known about the coevolution of KIR with HLA class I and the impact of their complexity on immune function and homeostasis.

Allele imputation for the killer cell immunoglobulin-like receptor KIR3DL1/S1.

Highly polymorphic interaction of KIR3DL1 and KIR3DS1 with HLA class I ligands modulates the effector functions of natural killer (NK) cells and some T cells. This genetically determined diversity affects severity of infections, immune-mediated diseases, and some cancers, and impacts the course of immunotherapies, including transplantation. KIR3DL1 is an inhibitory receptor, and KIR3DS1 is an activating receptor encoded by the KIR3DL1/S1 gene that has more than 200 diverse and divergent alleles. Determination of KIR3DL1/S1 genotypes for medical application is hampered by complex sequence and structural variation, requiring targeted approaches to generate and analyze high-resolution allele data. To overcome these obstacles, we developed and optimized a model for imputing KIR3DL1/S1 alleles at high-resolution from whole-genome SNP data. We designed the model to represent a substantial component of human genetic diversity. Our Global imputation model is effective at genotyping KIR3DL1/S1 alleles with an accuracy ranging from 88% in Africans to 97% in East Asians, with mean specificity of 99% and sensitivity of 95% for alleles >1% frequency. We used the established algorithm of the HIBAG program, in a modification named Pulling Out Natural killer cell Genomics (PONG). Because HIBAG was designed to impute HLA alleles also from whole-genome SNP data, PONG allows combinatorial diversity of KIR3DL1/S1 with HLA-A and -B to be analyzed using complementary techniques on a single data source. The use of PONG thus negates the need for targeted sequencing data in very large-scale association studies where such methods might not be tractable.

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.

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians.

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.

The combinatorial diversity of KIR and HLA class I allotypes in Peninsular Malaysia.

Killer cell immunoglobulin-like receptors (KIRs) interact with polymorphic human leucocyte antigen (HLA) class I molecules, modulating natural killer (NK) cell functions and affecting both the susceptibility and outcome of immune-mediated diseases. The KIR locus is highly diverse in gene content, copy number and allelic polymorphism within individuals and across geographical populations. To analyse currently under-represented Asian and Pacific populations, we investigated the combinatorial diversity of KIR and HLA class I in 92 unrelated Malay and 75 Malaysian Chinese individuals from the Malay Peninsula. We identified substantial allelic and structural diversity of the KIR locus in both populations and characterized novel variations at each analysis level. The Malay population is more diverse than Malay Chinese, likely representing a unique history including admixture with immigrating populations spanning several thousand years. Characterizing the Malay population are KIR haplotypes with large structural variants present in 10% individuals, and KIR and HLA alleles previously identified in Austronesian populations. Despite the differences in ancestries, the proportion of HLA allotypes that serve as KIR ligands is similar in each population. The exception is a significantly reduced frequency of interactions of KIR2DL1 with C2+ HLA-C in the Malaysian Chinese group, caused by the low frequency of C2+ HLA. One likely implication is a greater protection from preeclampsia, a pregnancy disorder associated with KIR2DL1, which shows higher incidence in the Malay than in the Malaysian Chinese. This first complete, high-resolution, characterization of combinatorial diversity of KIR and HLA in Malaysians will form a valuable reference for future clinical and population studies.

Host methylation predicts SARS-CoV-2 infection and clinical outcome.

BACKGROUND: Since the onset of the SARS-CoV-2 pandemic, most clinical testing has focused on RT-PCR1. Host epigenome manipulation post coronavirus infection2-4 suggests that DNA methylation signatures may differentiate patients with SARS-CoV-2 infection from uninfected individuals, and help predict COVID-19 disease severity, even at initial presentation. METHODS: We customized Illumina's Infinium MethylationEPIC array to enhance immune response detection and profiled peripheral blood samples from 164 COVID-19 patients with longitudinal measurements of disease severity and 296 patient controls. RESULTS: Epigenome-wide association analysis revealed 13,033 genome-wide significant methylation sites for case-vs-control status. Genes and pathways involved in interferon signaling and viral response were significantly enriched among differentially methylated sites. We observe highly significant associations at genes previously reported in genetic association studies (e.g. IRF7, OAS1). Using machine learning techniques, models built using sparse regression yielded highly predictive findings: cross-validated best fit AUC was 93.6% for case-vs-control status, and 79.1%, 80.8%, and 84.4% for hospitalization, ICU admission, and progression to death, respectively. CONCLUSIONS: In summary, the strong COVID-19-specific epigenetic signature in peripheral blood driven by key immune-related pathways related to infection status, disease severity, and clinical deterioration provides insights useful for diagnosis and prognosis of patients with viral infections.

Host methylation predicts SARS-CoV-2 infection and clinical outcome.

BACKGROUND: Since the onset of the SARS-CoV-2 pandemic, most clinical testing has focused on RT-PCR1. Host epigenome manipulation post coronavirus infection2-4 suggests that DNA methylation signatures may differentiate patients with SARS-CoV-2 infection from uninfected individuals, and help predict COVID-19 disease severity, even at initial presentation. METHODS: We customized Illumina's Infinium MethylationEPIC array to enhance immune response detection and profiled peripheral blood samples from 164 COVID-19 patients with longitudinal measurements of disease severity and 296 patient controls. RESULTS: Epigenome-wide association analysis revealed 13,033 genome-wide significant methylation sites for case-vs-control status. Genes and pathways involved in interferon signaling and viral response were significantly enriched among differentially methylated sites. We observe highly significant associations at genes previously reported in genetic association studies (e.g. IRF7, OAS1). Using machine learning techniques, models built using sparse regression yielded highly predictive findings: cross-validated best fit AUC was 93.6% for case-vs-control status, and 79.1%, 80.8%, and 84.4% for hospitalization, ICU admission, and progression to death, respectively. CONCLUSIONS: In summary, the strong COVID-19-specific epigenetic signature in peripheral blood driven by key immune-related pathways related to infection status, disease severity, and clinical deterioration provides insights useful for diagnosis and prognosis of patients with viral infections.

Viral infections affect the body in many ways, including via changes to the epigenome, the sum of chemical modifications to an individual's collection of genes that affect gene activity. Here, we analyzed the epigenome in blood samples from people with and without COVID-19 to determine whether we could find changes consistent with SARS-CoV-2 infection. Using a combination of statistical and machine learning techniques, we identify markers of SARS-CoV-2 infection as well as of severity and progression of COVID-19 disease. These signals of disease progression were present from the initial blood draw when first walking into the hospital. Together, these approaches demonstrate the potential of measuring the epigenome for monitoring SARS-CoV-2 status and severity.

Natural selection contributed to immunological differences between hunter-gatherers and agriculturalists.

The shift from a hunter-gatherer to an agricultural mode of subsistence is believed to have been associated with profound changes in the burden and diversity of pathogens across human populations. Yet, the extent to which the advent of agriculture affected the evolution of the human immune system remains unknown. Here we present a comparative study of variation in the transcriptional responses of peripheral blood mononuclear cells to bacterial and viral stimuli between Batwa rainforest hunter-gatherers and Bakiga agriculturalists from Uganda. We observed increased divergence between hunter-gatherers and agriculturalists in the early transcriptional response to viruses compared with that for bacterial stimuli. We demonstrate that a significant fraction of these transcriptional differences are under genetic control and we show that positive natural selection has helped to shape population differences in immune regulation. Across the set of genetic variants underlying inter-population immune-response differences, however, the signatures of positive selection were disproportionately observed in the rainforest hunter-gatherers. This result is counter to expectations on the basis of the popularized notion that shifts in pathogen exposure due to the advent of agriculture imposed radically heightened selective pressures in agriculturalist populations.

Polygenic adaptation and convergent evolution on growth and cardiac genetic pathways in African and Asian rainforest hunter-gatherers.

Different human populations facing similar environmental challenges have sometimes evolved convergent biological adaptations, for example, hypoxia resistance at high altitudes and depigmented skin in northern latitudes on separate continents. The "pygmy" phenotype (small adult body size), characteristic of hunter-gatherer populations inhabiting both African and Asian tropical rainforests, is often highlighted as another case of convergent adaptation in humans. However, the degree to which phenotypic convergence in this polygenic trait is due to convergent versus population-specific genetic changes is unknown. To address this question, we analyzed high-coverage sequence data from the protein-coding portion of the genomes of two pairs of populations: Batwa rainforest hunter-gatherers and neighboring Bakiga agriculturalists from Uganda and Andamanese rainforest hunter-gatherers and Brahmin agriculturalists from India. We observed signatures of convergent positive selection between the rainforest hunter-gatherers across the set of genes with "growth factor binding" functions ([Formula: see text]). Unexpectedly, for the rainforest groups, we also observed convergent and population-specific signatures of positive selection in pathways related to cardiac development (e.g., "cardiac muscle tissue development"; [Formula: see text]). We hypothesize that the growth hormone subresponsiveness likely underlying the adult small body-size phenotype may have led to compensatory changes in cardiac pathways, in which this hormone also plays an essential role. Importantly, in the agriculturalist populations, we did not observe similar patterns of positive selection on sets of genes associated with growth or cardiac development, indicating our results most likely reflect a history of convergent adaptation to the similar ecology of rainforests rather than a more general evolutionary pattern.

Development and validation of an arthropod maceration protocol for zoonotic pathogen detection in mosquitoes and fleas.

Arthropod-borne diseases remain a pressing international public health concern. While progress has been made in the rapid detection of arthropod-borne pathogens via quantitative real-time (qPCR), or even hand-held detection devices, a simple and robust maceration and nucleic acid extraction method is necessary to implement biosurveillance capabilities. In this study, a comparison of maceration techniques using five types of beads followed by nucleic acid extraction and detection were tested using two morphologically disparate arthropods, the Aedes aegypti mosquito and Xenopsylla spp. flea, to detect the zoonotic diseases dengue virus serotype-1 and Yersinia pestis. Post-maceration nucleic acid extraction was carried out using the 1-2-3 Platinum-Path-Sample-Purification (PPSP) kit followed by qPCR detection using the Joint Biological Agent Identification and Diagnostic System (JBAIDS). We found that the 5mm stainless steel beads added to the beads provided in the PPSP kit were successful in macerating the exoskeleton for both Ae. aegypti and Xenopsylla spp. Replicates in the maceration/extraction/detection protocol were increased in a stepwise fashion until a final 128 replicates were obtained. For dengue virus detection there was a 99% positivity rate and for Y. pestis detection there was a 95% positive detection rate. In the examination of both pathogens, there were no significant differences between qPCR instruments, days ran, time of day ran, or operators.

Plasmodium-specific molecular assays produce uninterpretable results and non-Plasmodium spp. sequences in field-collected Anopheles vectors.

The Malaria Research and Reference Reagent Resource-recommended PLF/UNR/VIR polymerase chain reaction (PCR) was used to detect Plasmodium vivax in Anopheles spp. mosquitoes collected in South Korea. Samples that were amplified were sequenced and compared with known Plasmodium spp. by using the PlasmoDB.org Basic Local Alignment Search Tool/n and the National Center for Biotechnology Information Basic Local Alignment Search Tool/n tools. Results show that the primers PLF/UNR/VIR used in this PCR can produce uninterpretable results and non-specific sequences in field-collected mosquitoes. Three additional PCRs (PLU/VIV, specific for 18S small subunit ribosomal DNA; Pvr47, specific for a nuclear repeat; and GDCW/PLAS, specific for the mitochondrial marker, cytB) were then used to find a more accurate and interpretable assay. Samples that were amplified were again sequenced. The PLU/VIV and Pvr47 assays showed cross-reactivity with non-Plasmodium spp. and an arthropod fungus (Zoophthora lanceolata). The GDCW/PLAS assay amplified only Plasmodium spp. but also amplified the non-human specific parasite P. berghei from an Anopheles belenrae mosquito. Detection of P. berghei in South Korea is a new finding.