A genomic perspective on HLA evolution

HLA diversity and signatures of selection in the Maniq, a nomadic hunter-gatherer population in Southern Thailand

The Maniq, a small and isolated nomadic hunter-gatherer population from the rainforests of Southern Thailand, offer a unique context for investigating how demographic history, genetic drift, and pathogen-driven selection shape human leucocyte antigen (HLA) diversity. Using high-coverage whole-genome data from 21 individuals (12 unrelated), we genotyped HLA alleles with HLA-HD and T1K, identifying 32 alleles in classical and 14 in non-classical HLA genes. Although overall HLA diversity was comparatively low, a few alleles at each locus occurred at high frequency, mirroring patterns observed in other small, isolated populations. Principal-component analysis clustered the Maniq with other Austroasiatic-speaking Semang hunter-gatherers (Jehai, Kintaq) on the Malay Peninsula and, intriguingly, with the Austronesian-speaking Tao of Taiwan, indicating shared immunogenetic features across linguistic boundaries. Despite reduced diversity, multiple loci bore signatures of both long-term balancing and recent positive selection. Several SNPs under selection were in complete linkage disequilibrium with eQTLs known to influence responses to hepatitis B virus (HBV) and other pathogens, suggesting that pathogen-driven pressure-in particular HBV-may have contributed to recent HLA evolution in the Maniq. These findings provide critical insights into how demographic constraints and pathogen landscapes converge to shape HLA diversity and evolution. In light of increasing infectious disease burdens in indigenous communities, our results underscore the importance of studying small, isolated populations to better understand the adaptive significance of HLA genes.

Novel HLA-A, -B, -C and -DRB1 alleles identified in the Australian New South Wales tissue typing laboratory

HLA compatibility between patients and donors is a determining factor for the success of stem cell and solid organ transplantations. However, finding suitable donors remains challenging due to the highly polymorphic nature of HLA genes. Here we are describing 16 novel HLA-A, -B, -C and -DRB1 alleles identified over the period of 2 years. Fourteen differed from their closest reference allele sequence by single nucleotide substitutions detected in the coding regions and 13 contained polymorphism outside antigen recognition domains (exons 2 and 3 in class I and exon 2 in HLA class II). Two novel alleles, a non-sense substitution and a deletion of 2 bp resulted in two null alleles, while a substitution in a splice site resulted in an allele with questionable expression status. In HLA matching procedures, particularly in donor selection, it is important to determine alternatively expressed HLA alleles. Thirteen additional HLA-A, -B, -C and -DRB1 alleles detected as novel at the time of testing were already reported by other laboratories by the time of submission.

Genetic Diversity and Ethnic Tapestry of Kazakhstan as Inferred from HLA Polymorphism and Population Dynamics: A Comprehensive Review

Background: The human leukocyte antigen (HLA) system represents the most polymorphic segment within human DNA sequences and constitutes a core component of immune defense responses and in understanding population genetics. This research investigates the distribution of HLA class I and II polymorphisms across different ethnic groups in Kazakhstan, offering valuable insights into the genetic diversity and demographic evolution within this region. Methods: We performed an in-depth examination of HLA class I and II polymorphisms across diverse ethnic communities living in Kazakhstan, including Kazakhs, Russians, Uzbeks, Ukrainians, Germans, Tatars, and Koreans. Utilizing data from high-resolution HLA typing studies allowed us to assess allele frequencies alongside haplotype distributions while analyzing genetic interrelations between these populations. Additionally, we performed comparative assessments with global HLA databases to determine the genetic affiliations between these groups and their relationships with neighboring and more distant populations. Results: Our study revealed over 200 HLA alleles within the analyzed populations, and significant variations were observed in their allele and haplotype frequencies. Notably, the Kazakh group exhibited strong genetic ties to Asian and Siberian demographics; conversely, other ethnicities showed associations reflective of their historical roots. Notable alleles included HLA-A*02:01, B*07:02, C*07:02, DRB1*07:01, and DQB1*03:01, commonly observed across various groups. Linkage disequilibrium analysis revealed the presence of population-specific haplotypes, highlighting distinct genetic structures within these communities. Conclusions: The findings highlight the significant genetic diversity in Kazakhstan, influenced by its geographical location at the crossroads of Europe and Asia. These results are pertinent to immunogenetics, transplantation medicine, and personalized healthcare within Kazakhstan and adjacent regions. Future research should expand the sample size and explore disease associations to enhance our comprehension of HLA genetics across Central Asia.

Natural selection and adaptive traits in the Maniq, a nomadic hunter-gatherer society from Mainland Southeast Asia

Asia is home to diverse hunter-gatherer populations characterized by significant morphological, anthropological, cultural, and linguistic diversity. Despite their importance in understanding ancestral human subsistence, little is known about the essential genetic adaptations of these groups. This study investigates the evolutionary pressures shaping the genome of the Maniq population, a nomadic hunter-gatherer group inhabiting the rainforests of southern Thailand. Using genome-wide approaches, including iHS, xp-EHH, PBE, and beta statistics, we identified signatures of positive and balancing selection. Genes under positive selection were enriched in pathways related to immunity, metabolic regulation, structural adaptation, cardiovascular performance, and neuromodulatory traits. Several genes associated with the Southeast Asian 'negrito-like' phenotype were also under positive selection. Balancing selection was primarily detected in immune-related genes, particularly within the HLA region, underscoring the critical role of genetic diversity in surviving pathogen-rich environments. Additionally, balancing selection in olfactory receptor genes highlights their importance in environmental sensing and adaptation. These results reveal the intricate interplay of positive and balancing selection in shaping the genetic landscape of the Maniq population and highlight their adaptations to the ecological and lifestyle challenges of life in the rainforest. This study contributes to our understanding of human evolutionary processes in tropical environments and hunter-gatherer societies.

Old vs. New Local Ancestry Inference in HCHS/SOL: A Comparative Study

Hispanic/Latino populations are admixed, with genetic contributions from multiple ancestral populations. Studies of genetic association in these admixed populations often use methods such as admixture mapping, which relies on inferred counts of "local" ancestry, i.e., of the source ancestral population at a locus. Local ancestries are inferred using external reference panels that represent ancestral populations, making the choice of inference method and reference panel critical. This study used a dataset of Hispanic/Latino individuals from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) to evaluate the "old" local ancestry inference performed using the state-of-the-art inference method, RFMix, alongside "new" inferences performed using Fast Local Ancestry Estimation (FLARE), which also used an updated reference panel. We compared their performance in terms of global and local ancestry correlations, as well as admixture mapping-based associations. Overall, the old RFMix and new FLARE inferences were highly similar for both global and local ancestries, with FLARE-inferred datasets yielding admixture mapping results consistent with those computed from RFMix. However, in some genomic regions the old and new local ancestries have relatively lower correlations (Pearson R < 0.9). Most of these genomic regions (86.42%) were mapped to either ENCODE blacklist regions, or to gene clusters, compared to 7.67% of randomly-matched regions with high correlations (Pearson R > 0.97) between old and new local ancestries.

Characterizing substructure via mixture modeling in large-scale genetic summary statistics

Genetic summary data are broadly accessible and highly useful, including for risk prediction, causal inference, fine mapping, and incorporation of external controls. However, collapsing individual-level data into summary data, such as allele frequencies, masks intra- and inter-sample heterogeneity, leading to confounding, reduced power, and bias. Ultimately, unaccounted-for substructure limits summary data usability, especially for understudied or admixed populations. There is a need for methods to enable the harmonization of summary data where the underlying substructure is matched between datasets. Here, we present Summix2, a comprehensive set of methods and software based on a computationally efficient mixture model to enable the harmonization of genetic summary data by estimating and adjusting for substructure. In extensive simulations and application to public data, we show that Summix2 characterizes finer-scale population structure, identifies ascertainment bias, and scans for potential regions of selection due to local substructure deviation. Summix2 increases the robust use of diverse, publicly available summary data, resulting in improved and more equitable research.

Cell state-dependent allelic effects and contextual Mendelian randomization analysis for human brain phenotypes

Gene expression quantitative trait loci are widely used to infer relationships between genes and central nervous system (CNS) phenotypes; however, the effect of brain disease on these inferences is unclear. Using 2,348,438 single-nuclei profiles from 391 disease-case and control brains, we report 13,939 genes whose expression correlated with genetic variation, of which 16.7-40.8% (depending on cell type) showed disease-dependent allelic effects. Across 501 colocalizations for 30 CNS traits, 23.6% had a disease dependency, even after adjusting for disease status. To estimate the unconfounded effect of genes on outcomes, we repeated the analysis using nondiseased brains (n = 183) and reported an additional 91 colocalizations not present in the larger mixed disease and control dataset, demonstrating enhanced interpretation of disease-associated variants. Principled implementation of single-cell Mendelian randomization in control-only brains identified 140 putatively causal gene-trait associations, of which 11 were replicated in the UK Biobank, prioritizing candidate peripheral biomarkers predictive of CNS outcomes.

Human immune system: Exploring diversity across individuals and populations

The immune response is an intricate system that involves the complex connection of cellular and molecular components, each with distinct functional specialisations. It has a distinct capacity to adjust and mould the immune response in accordance with specific stimuli, influenced by both genetic and environmental factors. The presence of genetic diversity, particularly across different ethnic and racial groups, significantly contributes to the impact of incidence of diseases, disease susceptibility, autoimmune disorders, and cancer risks in specific regions and certain populations. Environmental factors, including geography and socioeconomic status, further modulate the variety of the immune system responses. These, in turn, affect the susceptibility to infectious diseases and development of autoimmune disorders. Despite the complexity of the relationship, there remains a gap in understanding the specificity of immune indices across races, immune reference ranges among populations, highlighting the need for deeper understanding of immune diversity for personalized approaches in diagnostics and therapeutics. This review systematically organizes these findings, with the goal of emphasizing the potential of targeted interventions to address health disparities and advance translational research, enabling a more comprehensive strategy. This approach promises significant advancements in identifying specific immunological conditions, focusing on personalized interventions, through both genetic and environmental factors.

Evaluating the Utilities of Foundation Models in Single-cell Data Analysis

Foundation Models (FMs) have made significant strides in both industrial and scientific domains. In this paper, we evaluate the performance of FMs for single-cell sequencing data analysis through comprehensive experiments across eight downstream tasks pertinent to single-cell data. Overall, the top FMs include scGPT, Geneformer, and CellPLM by considering model performances and user accessibility among ten single-cell FMs. However, by comparing these FMs with task-specific methods, we found that single-cell FMs may not consistently excel than task-specific methods in all tasks, which challenges the necessity of developing foundation models for single-cell analysis. In addition, we evaluated the effects of hyper-parameters, initial settings, and stability for training single-cell FMs based on a proposed scEval framework, and provide guidelines for pre-training and fine-tuning, to enhance the performances of single-cell FMs. Our work summarizes the current state of single-cell FMs, points to their constraints and avenues for future development, and offers a freely available evaluation pipeline to benchmark new models and improve method development.

Full-resolution HLA and KIR gene annotations for human genome assemblies

The human leukocyte antigen (HLA) genes and the killer cell immunoglobulin-like receptor (KIR) genes are critical to immune responses and are associated with many immune-related diseases. Located in highly polymorphic regions, it is difficult to study them with traditional short-read alignment-based methods. Although modern long-read assemblers can often assemble these genes, using existing tools to annotate HLA and KIR genes in these assemblies remains a nontrivial task. Here, we describe Immuannot, a new computation tool to annotate the gene structures of HLA and KIR genes and to type the allele of each gene. Applying Immuannot to 56 regional and 212 whole-genome assemblies from previous studies, we annotate 9931 HLA and KIR genes and found that almost half of these genes, 4068, have novel sequences compared with the current Immuno Polymorphism Database (IPD). These novel gene sequences are represented by 2664 distinct alleles, some of which contained nonsynonymous variations, resulting in 92 novel protein sequences. We demonstrate the complex haplotype structures at the two loci and report the linkage between HLA/KIR haplotypes and gene alleles. We anticipate that Immuannot will speed up the discovery of new HLA/KIR alleles and enable the association of HLA/KIR haplotype structures with clinical outcomes in the future.

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome: the biology of a neglected disease

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is a chronic, debilitating disease characterised by a wide range of symptoms that severely impact all aspects of life. Despite its significant prevalence, ME/CFS remains one of the most understudied and misunderstood conditions in modern medicine. ME/CFS lacks standardised diagnostic criteria owing to variations in both inclusion and exclusion criteria across different diagnostic guidelines, and furthermore, there are currently no effective treatments available. Moving beyond the traditional fragmented perspectives that have limited our understanding and management of the disease, our analysis of current information on ME/CFS represents a significant paradigm shift by synthesising the disease's multifactorial origins into a cohesive model. We discuss how ME/CFS emerges from an intricate web of genetic vulnerabilities and environmental triggers, notably viral infections, leading to a complex series of pathological responses including immune dysregulation, chronic inflammation, gut dysbiosis, and metabolic disturbances. This comprehensive model not only advances our understanding of ME/CFS's pathophysiology but also opens new avenues for research and potential therapeutic strategies. By integrating these disparate elements, our work emphasises the necessity of a holistic approach to diagnosing, researching, and treating ME/CFS, urging the scientific community to reconsider the disease's complexity and the multifaceted approach required for its study and management.

The most frequent HLA alleles around the world: A fundamental synopsis

A comprehensive knowledge of human leukocyte antigen (HLA) molecular variation worldwide is essential in human population genetics research and disease association studies and is also indispensable for clinical applications such as allogeneic hematopoietic cell transplantation, where ensuring HLA compatibility between donors and recipients is paramount. Enormous progress has been made in this field thanks to several decades of HLA population studies allowing the development of helpful databases and bioinformatics tools. However, it is still difficult to appraise the global HLA population diversity in a synthetic way. We thus introduce here a novel approach, based on approximately 2000 data sets, to assess this complexity by providing a fundamental synopsis of the most frequent HLA alleles observed in different regions of the world. This new knowledge will be useful not only as a fundamental reference for basic research, but also as an efficient guide for clinicians working in the field of transplantation.

Characterizing substructure via mixture modeling in large-scale genetic summary statistics

Genetic summary data are broadly accessible and highly useful including for risk prediction, causal inference, fine mapping, and incorporation of external controls. However, collapsing individual-level data into groups masks intra- and inter-sample heterogeneity, leading to confounding, reduced power, and bias. Ultimately, unaccounted substructure limits summary data usability, especially for understudied or admixed populations. Here, we present Summix2, a comprehensive set of methods and software based on a computationally efficient mixture model to estimate and adjust for substructure in genetic summary data. In extensive simulations and application to public data, Summix2 characterizes finer-scale population structure, identifies ascertainment bias, and identifies potential regions of selection due to local substructure deviation. Summix2 increases the robust use of diverse publicly available summary data resulting in improved and more equitable research.

Two new HLA-DQB1*03:02:01 variants with substitutions in intron 2: HLA-DQB1*03:02:01:19 and -DQB1*03:02:01:21

Two different single nucleotide substitutions in intron 2 give rise to novel HLA-DQB1*03:02:01 alleles.

The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

Immune-driven clonal cell selection at the intersection among cancer, infections, autoimmunity and senescence

Immune surveillance mechanisms play a crucial role in maintaining lifelong immune homeostasis in response to pathologic stimuli and aberrant cell states. However, their persistence, especially in the context of chronic antigenic exposure, can create a fertile ground for immune evasion. These escaping cell phenotypes, harboring a variety of genomic and transcriptomic aberrances, chiefly in human leukocyte antigen (HLA) and antigen presentation machinery genes, may survive and proliferate, featuring a scenario of clonal cell expansion with immune failure characteristics. While well characterized in solid and, to some extent, hematological malignancies, little is known about their occurrence and significance in other disease contexts. Historical literature highlights the role for escaping HLA-mediated recognition as a strategy adopted by virus to evade from the immune system, hinting at the potential for immune aberrant cell expansion in the context of chronic infections. Additionally, unmasked in idiopathic aplastic anemia as a mechanism able to rescue failing hematopoiesis, HLA clonal escape may operate in autoimmune disorders, particularly in tissues targeted by aberrant immune responses. Furthermore, senescent cell status emerging as immunogenic phenotypes stimulating T cell responses, may act as a bottleneck for the selection of such immune escaping clones, blurring the boundaries between neoplastic transformation, aging and inflammation. Here we provide a fresh overview and perspective on this immune-driven clonal cell expansion, linking pathophysiological features of neoplastic, autoimmune, infectious and senescence processes exposed to immune surveillance.

Full resolution HLA and KIR genes annotation for human genome assemblies

The HLA (Human Leukocyte Antigen) genes and the KIR (Killer cell Immunoglobulin-like Receptor) genes are critical to immune responses and are associated with many immune-related diseases. Located in highly polymorphic regions, they are hard to be studied with traditional short-read alignment-based methods. Although modern long-read assemblers can often assemble these genes, using existing tools to annotate HLA and KIR genes in these assemblies remains a non-trivial task. Here, we describe Immuannot, a new computation tool to annotate the gene structures of HLA and KIR genes and to type the allele of each gene. Applying Immuannot to 56 regional and 212 whole-genome assemblies from previous studies, we annotated 9,931 HLA and KIR genes and found that almost half of these genes, 4,068, had novel sequences compared to the current Immuno Polymorphism Database (IPD). These novel gene sequences were represented by 2,664 distinct alleles, some of which contained non-synonymous variations resulting in 92 novel protein sequences. We demonstrated the complex haplotype structures at the two loci and reported the linkage between HLA/KIR haplotypes and gene alleles. We anticipate that Immuannot will speed up the discovery of new HLA/KIR alleles and enable the association of HLA/KIR haplotype structures with clinical outcomes in the future.

HLA Genes: A Hallmark of Functional Genetic Variation and Complex Evolution

The major histocompatibility complex (MHC) with its highly polymorphic HLA genes represents one of the most intensely studied genomic regions in the genome. MHC proteins play a key role in antigen-specific immunity and are associated with a wide range of complex diseases. Despite decades of research and many advances in the field, the characterization and interpretation of its genetic and genomic variability remain challenging. Here an overview is provided of the MHC, the nature of its exceptional variability, and the complex evolutionary processes assumed to drive this variability. Highlighted are also recent advances in the field that promise to improve our understanding of the variability in the MHC and in antigen-specific immunity more generally.

Association of Dental Caries in Primary Teeth With HLA Class II Gene Among Syrian Children

Background and objectives Numerous studies have shown that there is evidence that genetic factors contribute in different ways to the occurrence of caries lesions, and the multiple differences in the human leucocytes antigen (HLA) gene patterns play an important role in the body's immune response. Therefore, this study aimed to evaluate the association between some HLA class II alleles (HLA-DR4, HLA-DQ2, HLA-DQ4) and early childhood caries (ECC) occurrence in Syrian children. Materials and methods The study included 80 children aged three to six, divided into two groups: Group 1 (n = 40) being severe early childhood caries (S-ECC) children (deft ≥ 10) and Group 2 (n = 40) being free caries children (deft = 0). The genomic DNAs were extracted and collected by taking a buccal swab using a sterile metal strip and were amplified by polymerase chain reaction-single specific primer (PCR-SSP) and then HLA-typing was performed for all alleles. Results There were no statistically significant differences in the frequency of occurrence of S-ECC, HLA-DR4, and HLA-DQ2 alleles (p = 0.626, 0.256 respectively), while the incidence of S-ECC was associated with the DQ4 allele (p = 0.012). Conclusion HLA class II molecules may play a crucial role in predisposing to ECC, since positive HLA-DQ4 may enhance the chance of developing ECC. However, no association were found between (HLA-DR4 and HLA-DQ2) and ECC.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Design and deep learning of synthetic B-cell-specific promoters

Synthetic biology and deep learning synergistically revolutionize our ability for decoding and recoding DNA regulatory grammar. The B-cell-specific transcriptional regulation is intricate, and unlock the potential of B-cell-specific promoters as synthetic elements is important for B-cell engineering. Here, we designed and pooled synthesized 23 640 B-cell-specific promoters that exhibit larger sequence space, B-cell-specific expression, and enable diverse transcriptional patterns in B-cells. By MPRA (Massively parallel reporter assays), we deciphered the sequence features that regulate promoter transcriptional, including motifs and motif syntax (their combination and distance). Finally, we built and trained a deep learning model capable of predicting the transcriptional strength of the immunoglobulin V gene promoter directly from sequence. Prediction of thousands of promoter variants identified in the global human population shows that polymorphisms in promoters influence the transcription of immunoglobulin V genes, which may contribute to individual differences in adaptive humoral immune responses. Our work helps to decipher the transcription mechanism in immunoglobulin genes and offers thousands of non-similar promoters for B-cell engineering.

Human adipose tissue as a major reservoir of cytomegalovirus-reactive T cells

Introduction

Cytomegalovirus (CMV) is a common herpesvirus with a high prevalence worldwide. After the acute infection phase, CMV can remain latent in several tissues. CD8 T cells in the lungs and salivary glands mainly control its reactivation control. White adipose tissue (WAT) contains a significant population of memory T cells reactive to viral antigens, but CMV specificity has mainly been studied in mouse WAT. Therefore, we obtained blood, omental WAT (oWAT), subcutaneous WAT (sWAT), and liver samples from 11 obese donors to characterize the human WAT adaptive immune landscape from a phenotypic and immune receptor specificity perspective.

Methods

We performed high-throughput sequencing of the T cell receptor (TCR) locus to analyze tissue and blood TCR repertoires of the 11 donors. The presence of TCRs specific to CMV epitopes was tested through ELISpot assays. Moreover, phenotypic characterization of T cells was carried out through flow cytometry.

Results

High-throughput sequencing analyses revealed that tissue TCR repertoires in oWAT, sWAT, and liver samples were less diverse and dominated by hyperexpanded clones when compared to blood samples. Additionally, we predicted the presence of TCRs specific to viral epitopes, particularly from CMV, which was confirmed by ELISpot assays. Remarkably, we found that oWAT has a higher proportion of CMV-reactive T cells than blood or sWAT. Finally, flow cytometry analyses indicated that most WAT-infiltrated lymphocytes were tissue-resident effector memory CD8 T cells.

Discussion

Overall, these findings postulate human oWAT as a major reservoir of CMV-specific T cells, presumably for latent viral reactivation control. This study enhances our understanding of the adaptive immune response in human WAT and highlights its potential role in antiviral defense.

Joint multi-ancestry and admixed GWAS reveals the complex genetics behind human cranial vault shape

The cranial vault in humans is highly variable, clinically relevant, and heritable, yet its genetic architecture remains poorly understood. Here, we conduct a joint multi-ancestry and admixed multivariate genome-wide association study on 3D cranial vault shape extracted from magnetic resonance images of 6772 children from the ABCD study cohort yielding 30 genome-wide significant loci. Follow-up analyses indicate that these loci overlap with genomic risk loci for sagittal craniosynostosis, show elevated activity cranial neural crest cells, are enriched for processes related to skeletal development, and are shared with the face and brain. We present supporting evidence of regional localization for several of the identified genes based on expression patterns in the cranial vault bones of E15.5 mice. Overall, our study provides a comprehensive overview of the genetics underlying normal-range cranial vault shape and its relevance for understanding modern human craniofacial diversity and the etiology of congenital malformations.

HLA-DQ and alcohol in the pathogenesis of irritable bowel syndrome in college students: a case-control study

Some researchers have shown that genetics contribute to the incidence of IBS. However, no research has focused on the interaction between HLA-DQ and living habits in the pathogenesis of IBS. The present study explored the risk factors for IBS in college students of Guangxi Han nationality and explored the interaction between HLA-DQ and living habits on the pathogenesis of IBS. Univariate and multivariate analyses were used to determine the risk factors for IBS. Logistic interaction analysis and the Excel table made by Andersson were used to explore the interaction between genes and living habits in the context of IBS. We found that low expression of HLA-DQ2 and DQ8 were associated with the pathogenesis of IBS, while mild to moderate alcohol consumption was associated with the occurrence of IBS symptoms. Only the HLA-DQ8 gene and alcohol consumption had additive interactions in the context of the occurrence of IBS. In other words, for college students of Guangxi Han nationality, HLA-DQ2 and HLA-DQ8 might be protective against IBS, while alcohol consumption might be an independent risk factor. There was an additive interaction between HLA-DQ8 and alcohol consumption in the occurrence of IBS.

Comparative physiological anthropogeny: exploring molecular underpinnings of distinctly human phenotypes

Anthropogeny is a classic term encompassing transdisciplinary investigations of the origins of the human species. Comparative anthropogeny is a systematic comparison of humans and other living nonhuman hominids (so-called "great apes"), aiming to identify distinctly human features in health and disease, with the overall goal of explaining human origins. We begin with a historical perspective, briefly describing how the field progressed from the earliest evolutionary insights to the current emphasis on in-depth molecular and genomic investigations of "human-specific" biology and an increased appreciation for cultural impacts on human biology. While many such genetic differences between humans and other hominids have been revealed over the last two decades, this information remains insufficient to explain the most distinctive phenotypic traits distinguishing humans from other living hominids. Here we undertake a complementary approach of "comparative physiological anthropogeny," along the lines of the preclinical medical curriculum, i.e., beginning with anatomy and considering each physiological system and in each case considering genetic and molecular components that are relevant. What is ultimately needed is a systematic comparative approach at all levels from molecular to physiological to sociocultural, building networks of related information, drawing inferences, and generating testable hypotheses. The concluding section will touch on distinctive considerations in the study of human evolution, including the importance of gene-culture interactions.

The landscape of somatic mutations in lymphoblastoid cell lines

Somatic mutations have important biological ramifications while exerting substantial rate, type, and genomic location heterogeneity. Yet, their sporadic occurrence makes them difficult to study at scale and across individuals. Lymphoblastoid cell lines (LCLs), a model system for human population and functional genomics, harbor large numbers of somatic mutations and have been extensively genotyped. By comparing 1,662 LCLs, we report that the mutational landscape of the genome varies across individuals in terms of the number of mutations, their genomic locations, and their spectra; this variation may itself be modulated by somatic trans-acting mutations. Mutations attributed to the translesion DNA polymerase η follow two different modes of formation, with one mode accounting for the hypermutability of the inactive X chromosome. Nonetheless, the distribution of mutations along the inactive X chromosome appears to follow an epigenetic memory of the active form.

The impact of HLA polymorphism on herpesvirus infection and disease

Human Leukocyte Antigens (HLA) are cell surface molecules, central in coordinating innate and adaptive immune responses, that are targets of strong diversifying natural selection by pathogens. Of these pathogens, human herpesviruses have a uniquely ancient relationship with our species, where coevolution likely has reciprocating impact on HLA and viral genomic diversity. Consistent with this notion, genetic variation at multiple HLA loci is strongly associated with modulating immunity to herpesvirus infection. Here, we synthesize published genetic associations of HLA with herpesvirus infection and disease, both from case/control and genome-wide association studies. We analyze genetic associations across the eight human herpesviruses and identify HLA alleles that are associated with diverse herpesvirus-related phenotypes. We find that whereas most HLA genetic associations are virus- or disease-specific, HLA-A*01 and HLA-A*02 allotypes may be more generally associated with immune susceptibility and control, respectively, across multiple herpesviruses. Connecting genetic association data with functional corroboration, we discuss mechanisms by which diverse HLA and cognate receptor allotypes direct variable immune responses during herpesvirus infection and pathogenesis. Together, this review examines the complexity of HLA-herpesvirus interactions driven by differential T cell and Natural Killer cell immune responses.

HLA allele and haplotype frequencies of registered stem cell donors in Chile

Patients in need of hematopoietic stem cell transplantation often rely on unrelated stem cell donors matched in certain human leukocyte antigen (HLA) genes. Donor search is complicated by the extensive allelic variability of the HLA system. Therefore, large registries of potential donors are maintained in many countries worldwide. Population-specific HLA characteristics determine the registry benefits for patients and also the need for further regional donor recruitment. In this work, we analyzed HLA allele and haplotype frequencies of donors of DKMS Chile, the first Chilean donor registry, with self-assessed "non-Indigenous" (n=92,788) and "Mapuche" (n=1,993) ancestry. We identified HLA alleles that were distinctly more abundant in the Chilean subpopulations than in worldwide reference populations, four of them particularly characteristic for the Mapuche subpopulation, namely B*39:09g, B*35:09, DRB1*04:07g, and DRB1*16:02g. Both population subsamples carried haplotypes of both Native American and European origin at high frequencies, reflecting Chile's complex history of admixture and immigration. Matching probability analysis revealed limited benefits for Chilean patients (both non-Indigenous and Mapuche) from donor registries of non-Chilean donors, thus indicating a need for ongoing significant donor recruitment efforts in Chile.

A rarefaction approach for measuring population differences in rare and common variation

In studying allele-frequency variation across populations, it is often convenient to classify an allelic type as "rare," with nonzero frequency less than or equal to a specified threshold, "common," with a frequency above the threshold, or entirely unobserved in a population. When sample sizes differ across populations, however, especially if the threshold separating "rare" and "common" corresponds to a small number of observed copies of an allelic type, discreteness effects can lead a sample from one population to possess substantially more rare allelic types than a sample from another population, even if the two populations have extremely similar underlying allele-frequency distributions across loci. We introduce a rarefaction-based sample-size correction for use in comparing rare and common variation across multiple populations whose sample sizes potentially differ. We use our approach to examine rare and common variation in worldwide human populations, finding that the sample-size correction introduces subtle differences relative to analyses that use the full available sample sizes. We introduce several ways in which the rarefaction approach can be applied: we explore the dependence of allele classifications on subsample sizes, we permit more than two classes of allelic types of nonzero frequency, and we analyze rare and common variation in sliding windows along the genome. The results can assist in clarifying similarities and differences in allele-frequency patterns across populations.

Identifying signatures of positive selection in human populations from North Africa

Because of its location, North Africa (NA) has witnessed continuous demographic movements with an impact on the genomes of present-day human populations. Genomic data describe a complex scenario with varying proportions of at least four main ancestry components: Maghrebi, Middle Eastern-, European-, and West-and-East-African-like. However, the footprint of positive selection in NA has not been studied. Here, we compile genome-wide genotyping data from 190 North Africans and individuals from surrounding populations, investigate for signatures of positive selection using allele frequencies and linkage disequilibrium-based methods and infer ancestry proportions to discern adaptive admixture from post-admixture selection events. Our results show private candidate genes for selection in NA involved in insulin processing (KIF5A), immune function (KIF5A, IL1RN, TLR3), and haemoglobin phenotypes (BCL11A). We also detect signatures of positive selection related to skin pigmentation (SLC24A5, KITLG), and immunity function (IL1R1, CD44, JAK1) shared with European populations and candidate genes associated with haemoglobin phenotypes (HPSE2, HBE1, HBG2), other immune-related (DOCK2) traits, and insulin processing (GLIS3) traits shared with West and East African populations. Finally, the SLC8A1 gene, which codifies for a sodium-calcium exchanger, was the only candidate identified under post-admixture selection in Western NA.

Is HLA type a possible cancer risk modifier in Lynch syndrome?

Lynch syndrome (LS) is the most common inherited cancer syndrome. It is inherited via a monoallelic germline variant in one of the DNA mismatch repair (MMR) genes. LS carriers have a broad 30% to 80% risk of developing various malignancies, and more precise, individual risk estimations would be of high clinical value, allowing tailored cancer prevention and surveillance. Due to MMR deficiency, LS cancers are characterized by the accumulation of frameshift mutations leading to highly immunogenic frameshift peptides (FSPs). Thus, immune surveillance is proposed to inhibit the outgrowth of MMR-deficient cell clones. Recent studies have shown that immunoediting during the evolution of MMR-deficient cancers leads to a counter-selection of highly immunogenic antigens. The immunogenicity of FSPs is dependent on the antigen presentation. One crucial factor determining antigen presentation is the HLA genotype. Hence, a LS carrier's HLA genotype plays an important role in the presentation of FSP antigens to the immune system, and may influence the likelihood of progression from precancerous lesions to cancer. To address the challenge of clarifying this possibility including diverse populations with different HLA types, we have established the INDICATE initiative (Individual cancer risk by HLA type, http://indicate-lynch.org/), an international network aiming at a systematic evaluation of the HLA genotype as a possible cancer risk modifier in LS. Here we summarize the current knowledge on the role of HLA type in cancer risk and outline future research directions to delineate possible association in the scenario of LS with genetically defined risk population and highly immunogenic tumors.

Impact of Functional Polymorphisms on Drug Survival of Biological Therapies in Patients with Moderate-to-Severe Psoriasis

Biological therapies (BTs) indicated for psoriasis are highly effective; however, not all patients obtain good results, and loss of effectiveness is the main reason for switching. Genetic factors may be involved. The objective of this study was to evaluate the influence of single-nucleotide polymorphisms (SNPs) on the drug survival of tumor necrosis factor inhibitors (anti-TNF) medications and ustekinumab (UTK) in patients diagnosed with moderate-to-severe psoriasis. We conducted an ambispective observational cohort study that included 379 lines of treatment with anti-TNF (n = 247) and UTK (132) in 206 white patients from southern Spain and Italy. The genotyping of the 29 functional SNPs was carried out using real-time polymerase chain reaction (PCR) with TaqMan probes. Drug survival was evaluated with Cox regression and Kaplan-Meier curves. The multivariate analysis showed that the HLA-C rs12191877-T (hazard ratio [HR] = 0.560; 95% CI = 0.40-0.78; p = 0.0006) and TNF-1031 (rs1799964-C) (HR = 0.707; 95% CI = 0.50-0.99; p = 0.048) polymorphisms are associated with anti-TNF drug survival, while TLR5 rs5744174-G (HR = 0.589; 95% CI = 0.37-0.92; p = 0.02), CD84 rs6427528-GG (HR = 0.557; 95% CI = 0.35-0.88; p = 0.013) and PDE3A rs11045392-T together with SLCO1C1 rs3794271-T (HR = 0.508; 95% CI = 0.32-0.79; p = 0.002) are related to UTK survival. The limitations are the sample size and the clustering of anti-TNF drugs; we used a homogeneous cohort of patients from 2 hospitals only. In conclusion, SNPs in the HLA-C, TNF, TLR5, CD84, PDE3A, and SLCO1C1 genes may be useful as biomarkers of drug survival of BTs indicated for psoriasis, making it possible to implement personalized medicine that will reduce financial healthcare costs, facilitate medical decision-making and improve patient quality of life. However, further pharmacogenetic studies need to be conducted to confirm these associations.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation, and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here, we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues, using personalized reference genomes to mitigate technical confounding. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B, and T cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

New insights from GWAS on BMI-related growth traits in a longitudinal cohort of admixed children with Native American and European ancestry

Body-mass index (BMI) is a hallmark of adiposity. In contrast with adulthood, the genetic architecture of BMI during childhood is poorly understood. The few genome-wide association studies (GWAS) on children have been performed almost exclusively in Europeans and at single ages. We performed cross-sectional and longitudinal GWAS for BMI-related traits on 904 admixed children with mostly Mapuche Native American and European ancestries. We found regulatory variants of the immune gene HLA-DQB3 strongly associated with BMI at 1.5 - 2.5 years old. A variant in the sex-determining gene DMRT1 was associated with the age at adiposity rebound (Age-AR) in girls (P = 9.8 × 10 - 9 ). BMI was significantly higher in Mapuche than in Europeans between 5.5 and 16.5 years old. Finally, Age-AR was significantly lower (P = 0.004 ) by 1.94 years and BMI at AR was significantly higher (P = 0.04 ) by 1.2 kg/m 2 , in Mapuche children compared with Europeans.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

Trans-ethnic Polygenic Risk Scores for Body Mass Index: An International Hundred K+ Cohorts Consortium Study

Background

While polygenic risk scores hold significant promise in estimating an individual's risk of developing a complex trait such as obesity, their application in the clinic has, to date, been limited by a lack of data from non-European populations. As a collaboration model of the International Hundred K+ Cohorts Consortium (IHCC), we endeavored to develop a globally applicable trans-ethnic PRS for body mass index (BMI) through this relatively new international effort.

Methods

The PRS model was developed trained and tested at the Center for Applied Genomics (CAG) of The Children's Hospital of Philadelphia (CHOP) based on a BMI meta-analysis from the GIANT consortium. The validated PRS models were subsequently disseminated to the participating sites. Scores were generated by each site locally on their cohorts and summary statistics returned to CAG for final analysis.

Results

We show that in the absence of a well powered trans-ethnic GWAS from which to derive SNPs and effect estimates, trans-ethnic scores can be generated from European ancestry GWAS using Bayesian approaches such as LDpred to adjust the summary statistics using trans-ethnic linkage disequilibrium reference panels. The ported trans-ethnic scores outperform population specific-PRS across all non-European ancestry populations investigated including East Asians and three-way admixed Brazilian cohort.

Conclusions

Widespread use of PRS in the clinic is hampered by a lack of genotyping data in individuals of non-European ancestry for the vast majority of traits. Here we show that for a truly polygenic trait such as BMI adjusting the summary statistics of a well powered European ancestry study using trans-ethnic LD reference results in a score that is predictive across a range of ancestries including East Asians and three-way admixed Brazilians.

Evaluation of Everolimus Activity against Mycobacterium tuberculosis Using In Vitro Models of Infection

Even though Everolimus has been investigated in a phase II randomized trial as a host-directed therapy (HDT) to treat tuberculosis (TB), an oncological patient treated with Everolimus for a neuroendocrine pancreatic neoplasia developed active TB twice and a non-tuberculous mycobacterial (NTM) infection in a year and a half time span. To investigate this interesting case, we isolated and genotypically characterized the Mycobacterium tuberculosis (Mtb) clinical strain from the patient and tested the effect of Everolimus on its viability in an axenic culture and in a peripheral blood mononuclear cell (PBMCs) infection model. To exclude strain-specific resistance, we tested the activity of Everolimus against Mtb strains of ancient and modern lineages. Furthermore, we investigated the Everolimus effect on ROS production and autophagy modulation during Mtb infection. Everolimus did not have a direct effect on mycobacteria viability and a negligible effect during Mtb infection in host cells, although it stimulated autophagy and ROS production. Despite being a biologically plausible HDT against TB, Everolimus does not exert a direct or indirect activity on Mtb. This case underlines the need for a careful approach to drug repurposing and implementation and the importance of pre-clinical experimental studies.

The evolutionary and functional significance of germline immunoglobulin gene variation

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

Immunogenetic Predisposition to SARS-CoV-2 Infection

Herein, we included 527 individuals from two Hospitals, Chemnitz and University-Hospital Leipzig. In total, 199 were negative for PCR and 328 were positive upon first admission. We used next generation sequencing for HLA-A, B, C, DRB1, DRB345, DQA1, DQB1, DPA1, and DPB1, and in some cases, HLA-E, F, G, and H. Furthermore, we molecularly defined 22 blood group systems comprising 26 genes and 5 platelet antigen genes. We observed a significant enrichment of homozygosity for DQA/DQB in the positive group. Within the negative subjects, HLA-B*57:01, HLA-B*55:01, DRB1*13:01, and DRB1*01:01 were enriched, and in the positive group, homozygosity for DQA/DQB, DRB1*09:01, and DRB1*15:01 was observed. DQA1*01:01, DQA1*02:01, and DQA1*01:03 were enriched in the negative group. HLA-DQB1*06:02 was enriched in the positive group, and HLA-DQB1*05:01 and HLA-DQB1*06:03 were enriched in the negative group. For the blood group systems MNS, RH, LE, FY, JK, YT, DO, and KN, enrichment was seen in both groups, depending on the antigen under observation. Homozygosity for D-positive RHD alleles, as well as the phenotypes M-N+ of the MNS blood group system and Yk(a-) of the KN system, were enriched in the positive group. All of these significances disappeared upon correction. Subjects who carried homozygous HPA-1a were more frequent in the negative group, contrasting with the finding that HPA-1ab was enriched in the positive group.

Human leukocyte antigen super-locus: nexus of genomic supergenes, SNPs, indels, transcripts, and haplotypes

The human Major Histocompatibility Complex (MHC) or Human Leukocyte Antigen (HLA) super-locus is a highly polymorphic genomic region that encodes more than 140 coding genes including the transplantation and immune regulatory molecules. It receives special attention for genetic investigation because of its important role in the regulation of innate and adaptive immune responses and its strong association with numerous infectious and/or autoimmune diseases. In recent years, MHC genotyping and haplotyping using Sanger sequencing and next-generation sequencing (NGS) methods have produced many hundreds of genomic sequences of the HLA super-locus for comparative studies of the genetic architecture and diversity between the same and different haplotypes. In this special issue on 'The Current Landscape of HLA Genomics and Genetics', we provide a short review of some of the recent analytical developments used to investigate the SNP polymorphisms, structural variants (indels), transcription and haplotypes of the HLA super-locus. This review highlights the importance of using reference cell-lines, population studies, and NGS methods to improve and update our understanding of the mechanisms, architectural structures and combinations of human MHC genomic alleles (SNPs and indels) that better define and characterise haplotypes and their association with various phenotypes and diseases.

The immunogenetics of COVID-19

The worldwide coronavirus disease 2019 pandemic was sparked by the severe acute respiratory syndrome caused by coronavirus 2 (SARS-CoV-2) that first surfaced in December 2019 (COVID-19). The effects of COVID-19 differ substantially not just between patients individually but also between populations with different ancestries. In humans, the human leukocyte antigen (HLA) system coordinates immune regulation. Since HLA molecules are a major component of antigen-presenting pathway, they play an important role in determining susceptibility to infectious disease. It is likely that differential susceptibility to SARS-CoV-2 infection and/or disease course in COVID-19 in different individuals could be influenced by the variations in the HLA genes which are associated with various immune responses to SARS-CoV-2. A growing number of studies have identified a connection between HLA variation and diverse COVID-19 outcomes. Here, we review research investigating the impact of HLA on individual responses to SARS-CoV-2 infection and/or progression, also discussing the significance of MHC-related immunological patterns and its use in vaccine design.

Heterogeneity in disease resistance and the impact of antibiotics in the US

We hypothesize that the impact of antibiotics is moderated by a population's inherent (genetic) resistance to infectious disease. Using the introduction of sulfa drugs in 1937, we show that US states that are more genetically susceptible to infectious disease saw larger declines in their bacterial mortality rates following the introduction of sulfa drugs in 1937. This suggests area-level genetic endowments of disease resistance and the discovery of medical technologies have acted as substitutes in determining levels of health across the US. We also document immediate effects of sulfa drug exposure to the age of the workforce and cumulative effects on educational attainment for cohorts exposed to sulfa drugs in early life.

Individual HLA heterogeneity and its implications for cellular immune evasion in cancer and beyond

Structural and functional variability of human leukocyte antigen (HLA) is the foundation for competent adaptive immune responses against pathogen and tumor antigens as it assures the breadth of the presented immune-peptidome, theoretically sustaining an efficient and diverse T cell response. This variability is presumably the result of the continuous selection by pathogens, which over the course of evolution shaped the adaptive immune system favoring the assortment of a hyper-polymorphic HLA system able to elaborate efficient immune responses. Any genetic alteration affecting this diversity may lead to pathological processes, perturbing antigen presentation capabilities, T-cell reactivity and, to some extent, natural killer cell functionality. A highly variable germline HLA genotype can convey immunogenetic protection against infections, be associated with tumor surveillance or influence response to anti-neoplastic treatments. In contrast, somatic aberrations of HLA loci, rearranging the original germline configuration, theoretically decreasing its variability, can facilitate mechanisms of immune escape that promote tumor growth and immune resistance. The purpose of the present review is to provide a unified and up-to-date overview of the pathophysiological consequences related to the perturbations of the genomic heterogeneity of HLA complexes and their impact on human diseases, with a special focus on cancer.

Identifying signatures of natural selection in Indian populations

In this study, we present the results of a genome-wide scan for signatures of positive selection using data from four tribal groups (Kokana, Warli, Bhil, and Pawara) and two caste groups (Deshastha Brahmin and Kunbi Maratha) from West of the Maharashtra State In India, as well as two samples of South Asian ancestry from the 1KG project (Gujarati Indian from Houston, Texas and Indian Telugu from UK). We used an outlier approach based on different statistics, including PBS, xpEHH, iHS, CLR, Tajima's D, as well as two recently developed methods: Graph-aware Retrieval of Selective Sweeps (GRoSS) and Ascertained Sequentially Markovian Coalescent (ASMC). In order to minimize the risk of false positives, we selected regions that are outliers in all the samples included in the study using more than one method. We identified putative selection signals in 107 regions encompassing 434 genes. Many of the regions overlap with only one gene. The signals observed using microarray-based data are very consistent with our analyses using high-coverage sequencing data, as well as those identified with a novel coalescence-based method (ASMC). Importantly, at least 24 of these genomic regions have been identified in previous selection scans in South Asian populations or in other population groups. Our study highlights genomic regions that may have played a role in the adaptation of anatomically modern humans to novel environmental conditions after the out of Africa migration.

Type 1 diabetes in diverse ancestries and the use of genetic risk scores

Over 75 genetic loci within and outside of the HLA region influence type 1 diabetes risk. Genetic risk scores (GRS), which facilitate the integration of complex genetic information, have been developed in type 1 diabetes and incorporated into models and algorithms for classification, prognosis, and prediction of disease and response to preventive and therapeutic interventions. However, the development and validation of GRS across different ancestries is still emerging, as is knowledge on type 1 diabetes genetics in populations of diverse genetic ancestries. In this Review, we provide a summary of the current evidence on the evolutionary genetic variation in type 1 diabetes and the racial and ethnic differences in type 1 diabetes epidemiology, clinical characteristics, and preclinical course. We also discuss the influence of genetics on type 1 diabetes with differences across ancestries and the development and validation of GRS in various populations.

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.

Admixture Mapping of Alzheimer's disease in Caribbean Hispanics identifies a new locus on 22q13.1

Late-onset Alzheimer's disease (LOAD) is significantly more frequent in Hispanics than in non-Hispanic Whites. Ancestry may explain these differences across ethnic groups. To this end, we studied a large cohort of Caribbean Hispanics (CH, N = 8813) and tested the association between Local Ancestry (LA) and LOAD ("admixture mapping") to identify LOAD-associated ancestral blocks, separately for ancestral components (European [EUR], African [AFR], Native American[NA]) and jointly (AFR + NA). Ancestral blocks significant after permutation were fine-mapped employing multi-ethnic whole-exome sequencing (WES) to identify rare variants associated with LOAD (SKAT-O) and replicated in the UK Biobank WES dataset. Candidate genes were validated studying (A) protein expression in human LOAD and control brains; (B) two animal AD models, Drosophila and Zebrafish. In the joint AFR + NA model, we identified four significant ancestral blocks located on chromosomes 1 (p value = 8.94E-05), 6 (p value = 8.63E-05), 21 (p value = 4.64E-05) and 22 (p value = 1.77E-05). Fine-mapping prioritized the GCAT gene on chromosome 22 (SKAT-O p value = 3.45E-05) and replicated in the UK Biobank (SKAT-O p value = 0.05). In LOAD brains, a decrease of 28% in GCAT protein expression was observed (p value = 0.038), and GCAT knockdown in Amyloid-β42 Drosophila exacerbated rough eye phenotype (68% increase, p value = 4.84E-09). In zebrafish, gcat expression increased after acute amyloidosis (34%, p value = 0.0049), and decreased upon anti-inflammatory Interleukin-4 (39%, p value = 2.3E-05). Admixture mapping uncovered genomic regions harboring new LOAD-associated loci that might explain the observed different frequency of LOAD across ethnic groups. Our results suggest that the inflammation-related activity of GCAT is a response to amyloid toxicity, and reduced GCAT expression exacerbates AD pathology.

Human Leukocyte Antigen (HLA) System: Genetics and Association with Bacterial and Viral Infections

The human leukocyte antigen (HLA) system is one of the most crucial host factors influencing disease progression in bacterial and viral infections. This review provides the basic concepts of the structure and function of HLA molecules in humans. Here, we highlight the main findings on the associations between HLA class I and class II alleles and susceptibility to important infectious diseases such as tuberculosis, leprosy, melioidosis, Staphylococcus aureus infection, human immunodeficiency virus infection, coronavirus disease 2019, hepatitis B, and hepatitis C in populations worldwide. Finally, we discuss challenges in HLA typing to predict disease outcomes in clinical implementation. Evaluation of the impact of HLA variants on the outcome of bacterial and viral infections would improve the understanding of pathogenesis and identify those at risk from infectious diseases in distinct populations and may improve the individual treatment.

Characterizing the diversity of MHC conserved extended haplotypes using families from the United Arab Emirates

Aside from its anthropological relevance, the characterization of the allele frequencies of genes in the human Major Histocompatibility Complex (MHC) and the combination of these alleles that make up MHC conserved extended haplotypes (CEHs) is necessary for histocompatibility matching in transplantation as well as mapping disease association loci. The structure and content of the MHC region in Middle Eastern populations remain poorly characterized, posing challenges when establishing disease association studies in ethnic groups that inhabit the region and reducing the capacity to translate genetic research into clinical practice. This study was conceived to address a gap of knowledge, aiming to characterize CEHs in the United Arab Emirates (UAE) population through segregation analysis of high-resolution, pedigree-phased, MHC haplotypes derived from 41 families. Twenty per cent (20.5%) of the total haplotype pool derived from this study cohort were identified as putative CEHs in the UAE population. These consisted of CEHs that have been previously detected in other ethnic groups, including the South Asian CEH 8.2 [HLA- C*07:02-B*08:01-DRB1*03:01-DQA1*05:01-DQB1*02:01 (H.F. 0.094)] and the common East Asian CEH 58.1 [HLA- C*03:02-B*58:01-DRB1*03:01- DQA1*05:01-DQB1*02:01 (H.F. 0.024)]. Additionally, three novel CEHs were identified in the current cohort, including HLA- C*15:02-B*40:06-DRB1*16:02-DQB1*05:02 (H.F. 0.035), HLA- C*16:02-B*51:01-DRB1*16:01-DQA1*01:02-DQB1*05:02 (H.F. 0.029), and HLA- C*03:02-B*58:01-DRB1*16:01-DQA1*01:02-DQB1*05:02 (H.F. 0.024). Overall, the results indicate a substantial gene flow with neighbouring ethnic groups in the contemporary UAE population including South Asian, East Asian, African, and European populations. Importantly, alleles and haplotypes that have been previously associated with autoimmune diseases (e.g., Type 1 Diabetes) were also present. In this regard, this study emphasizes that an appreciation for ethnic differences can provide insights into subpopulation-specific disease-related polymorphisms, which has remained a difficult endeavour.

How HLA diversity is apportioned: influence of selection and relevance to transplantation

In his 1972 paper ‘The apportionment of human diversity’, Lewontin showed that, when averaged over loci, genetic diversity is predominantly attributable to differences among individuals within populations. However, selection can alter the apportionment of diversity of specific genes or genomic regions. We examine genetic diversity at the human leucocyte antigen (HLA) loci, located within the major histocompatibility complex (MHC) region. HLA genes code for proteins that are critical to adaptive immunity and are well-documented targets of balancing selection. The single-nucleotide polymorphisms (SNPs) within HLA genes show strong signatures of balancing selection on large timescales and are broadly shared among populations, displaying low F_ST values. However, when we analyse haplotypes defined by these SNPs (which define ‘HLA alleles’), we find marked differences in frequencies between geographic regions. These differences are not reflected in the F_ST values because of the extreme polymorphism at HLA loci, illustrating challenges in interpreting F_ST. Differences in the frequency of HLA alleles among geographic regions are relevant to bone-marrow transplantation, which requires genetic identity at HLA loci between patient and donor. We discuss the case of Brazil's bone marrow registry, where a deficit of enrolled volunteers with African ancestry reduces the chance of finding donors for individuals with an MHC region of African ancestry.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.