Major Histocompatibility Complex Region Research Articles

Modeling of horizontal pleiotropy identifies possible causal gene expression in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a chronic autoimmune condition with complex causes involving genetic and environmental factors. While genome-wide association studies (GWASs) have identified genetic loci associated with SLE, the functional genomic elements responsible for disease development remain largely unknown. Mendelian Randomization (MR) is an instrumental variable approach to causal inference based on data from observational studies, where genetic variants are employed as instrumental variables (IVs). This study utilized a two-step strategy to identify causal genes for SLE. In the first step, the classical MR method was employed, assuming the absence of horizontal pleiotropy, to estimate the causal effect of gene expression on SLE. In the second step, advanced probabilistic MR methods (PMR-Egger, MRAID, and MR-MtRobin) were applied to the genes identified in the first step, considering horizontal pleiotropy, to filter out false positives. PMR-Egger and MRAID analyses utilized whole blood expression quantitative trait loci (eQTL) and SLE GWAS summary data, while MR-MtRobin analysis used an independent eQTL dataset from multiple immune cell types along with the same SLE GWAS data. The initial MR analysis identified 142 genes, including 43 outside of chromosome 6. Subsequently, applying the advanced MR methods reduced the number of genes with significant causal effects on SLE to 66. PMR-Egger, MRAID, and MR-MtRobin, respectively, identified 13, 7, and 16 non-chromosome 6 genes with significant causal effects. All methods identified expression of PHRF1 gene as causal for SLE. A comprehensive literature review was conducted to enhance understanding of the functional roles and mechanisms of the identified genes in SLE development. The findings from the three MR methods exhibited overlapping genes with causal effects on SLE, demonstrating consistent results. However, each method also uncovered unique genes due to different modelling assumptions and technical factors, highlighting the complementary nature of the approaches. Importantly, MRAID demonstrated a reduced percentage of causal genes from the Major Histocompatibility complex (MHC) region on chromosome 6, indicating its potential in minimizing false positive findings. This study contributes to unraveling the mechanisms underlying SLE by employing advanced probabilistic MR methods to identify causal genes, thereby enhancing our understanding of SLE pathogenesis.

A case of forensic genomics in Uganda reveals animal ownership and low exotic genetic introgression in indigenous cattle.

The cattle industry contributes to Uganda's agricultural output. It faces challenges that include theft and parentage ascertainment. These challenges can benefit from recent molecular genomics and bioinformatics technologies. We employed genomic analyses to establish potential ownership of a group of nine cattle that were being claimed by two farmers in Uganda. We investigated the genetic relationship of Ugandan cattle with regional indigenous breeds as well as exotic breeds that are currently present in Uganda. In addition, we investigated regions that are likely to be under selection in the Ugandan cattle. Hair samples were collected from seven and two animals from farmers A and B, respectively. They were genotyped for 53,218 Single Nucleotide Polymorphism markers. To establish genetic relationships between the sampled animals, we performed genomic analyses including, principal component analysis (PCA), hierarchical clustering analysis and identity by state/descent. We also performed admixture and runs of homozygosity analyses to assess the ancestry composition and identify regions potentially under selection in Ugandan cattle, respectively. The seven animals from Farmer A were genetically close to each other but showed minimal relationship with the disputed animals. The two animals from Farmer B were genetically distant from each other but showed greater similarity to four of the disputed animals. Four of the disputed animals showed great dissimilarity from the animals of both farmers. Comparison of these with the reference breeds revealed minimal European exotic genetic introgression into these animals, but rather high similarity to the Sheko. Results also revealed high homozygosity in the major histocompatibility complex regions. Our results demonstrate the use of currently available genomic tools to empirically establish the ownership of cattle; these could be scaled up as a resourceful and viable tool that could be employed to support conflict resolution where reliable livestock identification is unavailable.

Evolutionary trade-off between innate and acquired immune defences in birds

BackgroundThe development, maintenance, and use of immune defences are costly. Therefore, animals face trade-offs in terms of resource allocation within their immune system and between their immune system and other physiological processes. To maximize fitness, evolution may favour investment in one immunological defence or subsystem over another in a way that matches a species broader life history strategy. Here, we used phylogenetically-informed comparative analyses to test for relationships between two immunological components. Natural antibodies and complement were used as proxies for the innate branch; structural complexity of the major histocompatibility complex (MHC) region was used for the acquired branch.ResultsWe found a negative association between the levels of natural antibodies (i.e., haemagglutination titre) and the total MHC gene copy number across the avian phylogeny, both at the species and family level. The family-level analysis indicated that this association was apparent for both MHC-I and MHC-II, when copy numbers within these two MHC regions were analysed separately. The association remained significant after controlling for basic life history components and for ecological traits commonly linked to pathogen exposure.ConclusionOur results provide the first phylogenetically robust evidence for an evolutionary trade-off within the avian immune system, with a more developed acquired immune system (i.e., more complex MHC architecture) in more derived bird lineages (e.g., passerines) being accompanied by an apparent downregulation of the innate immune system.

Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons.

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

Genome-wide association study identifies risk loci within the major histocompatibility complex region for Hunner-type interstitial cystitis

Hunner-type interstitial cystitis (HIC) is a rare, chronic inflammatory disease of the urinary bladder with unknown etiology and genetic background. Here, we conduct a genome-wide association study of 144 patients with HIC and 41,516 controls of Japanese ancestry. The genetic variant, rs1794275, in the major histocompatibility complex (MHC) region (chromosome 6p21.3) is associated with HIC risk (odds ratio [OR]= 2.32; p= 3.4×10-9). The association is confirmed in a replication set of 26 cases and 1,026 controls (p=0.014). Fine mapping demonstrates the contribution to the disease risk of a completely linked haplotype of three human leukocyte antigen HLA-DQβ1 amino acid positions, 71, 74, and 75 (OR= 1.94; p= 5×10-8) and of HLA-DPβ1 amino acid position 178, which tags HLA-DPB1∗04:02 (OR= 2.35; p= 7.5×10-8). The three HLA-DQβ1 amino acid positions are located together at the peptide binding groove, suggesting their functional importance in antigen presentation. Our study reveals genetic contributions to HIC risk that may be associated with class II MHC molecule antigen presentation.

Race/ethnicity-stratified fine-mapping of the MHC locus reveals genetic variants associated with late-onset asthma.

Introduction: Asthma is a chronic disease of the airways that impairs normal breathing. The etiology of asthma is complex and involves multiple factors, including the environment and genetics, especially the distinct genetic architecture associated with ancestry. Compared to early-onset asthma, little is known about genetic predisposition to late-onset asthma. We investigated the race/ethnicity-specific relationship among genetic variants within the major histocompatibility complex (MHC) region and late-onset asthma in a North Carolina-based multiracial cohort of adults. Methods: We stratified all analyses by self-reported race (i.e., White and Black) and adjusted all regression models for age, sex, and ancestry. We conducted association tests within the MHC region and performed fine-mapping analyses conditioned on the race/ethnicity-specific lead variant using whole-genome sequencing (WGS) data. We applied computational methods to infer human leukocyte antigen (HLA) alleles and residues at amino acid positions. We replicated findings in the UK Biobank. Results: The lead signals, rs9265901 on the 5' end of HLA-B, rs55888430 on HLA-DOB, and rs117953947 on HCG17, were significantly associated with late-onset asthma in all, White, and Black participants, respectively (OR = 1.73, 95%CI: 1.31 to 2.14, p = 3.62 × 10-5; OR = 3.05, 95%CI: 1.86 to 4.98, p = 8.85 × 10-6; OR = 19.5, 95%CI: 4.37 to 87.2, p = 9.97 × 10-5, respectively). For the HLA analysis, HLA-B*40:02 and HLA-DRB1*04:05, HLA-B*40:02, HLA-C*04:01, and HLA-DRB1*04:05, and HLA-DRB1*03:01 and HLA-DQB1 were significantly associated with late-onset asthma in all, White, and Black participants. Conclusion: Multiple genetic variants within the MHC region were significantly associated with late-onset asthma, and the associations were significantly different by race/ethnicity group.

Identifying potential risk genes and pathways for neuropsychiatric and substance use disorders using intermediate molecular mediator information.

Neuropsychiatric and substance use disorders (NPSUDs) have a complex etiology that includes environmental and polygenic risk factors with significant cross-trait genetic correlations. Genome-wide association studies (GWAS) of NPSUDs yield numerous association signals. However, for most of these regions, we do not yet have a firm understanding of either the specific risk variants or the effects of these variants. Post-GWAS methods allow researchers to use GWAS summary statistics and molecular mediators (transcript, protein, and methylation abundances) infer the effect of these mediators on risk for disorders. One group of post-GWAS approaches is commonly referred to as transcriptome/proteome/methylome-wide association studies, which are abbreviated as T/P/MWAS (or collectively as XWAS). Since these approaches use biological mediators, the multiple testing burden is reduced to the number of genes (∼20,000) instead of millions of GWAS SNPs, which leads to increased signal detection. In this work, our aim is to uncover likely risk genes for NPSUDs by performing XWAS analyses in two tissues-blood and brain. First, to identify putative causal risk genes, we performed an XWAS using the Summary-data-based Mendelian randomization, which uses GWAS summary statistics, reference xQTL data, and a reference LD panel. Second, given the large comorbidities among NPSUDs and the shared cis-xQTLs between blood and the brain, we improved XWAS signal detection for underpowered analyses by performing joint concordance analyses between XWAS results i) across the two tissues and ii) across NPSUDs. All XWAS signals i) were adjusted for heterogeneity in dependent instruments (HEIDI) (non-causality) p-values and ii) used to test for pathway enrichment. The results suggest that there were widely shared gene/protein signals within the major histocompatibility complex region on chromosome 6 (BTN3A2 and C4A) and elsewhere in the genome (FURIN, NEK4, RERE, and ZDHHC5). The identification of putative molecular genes and pathways underlying risk may offer new targets for therapeutic development. Our study revealed an enrichment of XWAS signals in vitamin D and omega-3 gene sets. So, including vitamin D and omega-3 in treatment plans may have a modest but beneficial effect on patients with bipolar disorder.

Benchmark of tools for in silico prediction of MHC class I and class II genotypes from NGS data

BackgroundThe Human Leukocyte Antigen (HLA) genes are a group of highly polymorphic genes that are located in the Major Histocompatibility Complex (MHC) region on chromosome 6. The HLA genotype affects the presentability of tumour antigens to the immune system. While knowledge of these genotypes is of utmost importance to study differences in immune responses between cancer patients, gold standard, PCR-derived genotypes are rarely available in large Next Generation Sequencing (NGS) datasets. Therefore, a variety of methods for in silico NGS-based HLA genotyping have been developed, bypassing the need to determine these genotypes with separate experiments. However, there is currently no consensus on the best performing tool.ResultsWe evaluated 13 MHC class I and/or class II HLA callers that are currently available for free academic use and run on either Whole Exome Sequencing (WES) or RNA sequencing data. Computational resource requirements were highly variable between these tools. Three orthogonal approaches were used to evaluate the accuracy on several large publicly available datasets: a direct benchmark using PCR-derived gold standard HLA calls, a correlation analysis with population-based allele frequencies and an analysis of the concordance between the different tools. The highest MHC-I calling accuracies were found for Optitype (98.0%) and arcasHLA (99.4%) on WES and RNA sequencing data respectively, while for MHC-II HLA-HD was the most accurate tool for both data types (96.2% and 99.4% on WES and RNA data respectively).ConclusionThe optimal strategy for HLA genotyping from NGS data depends on the availability of either WES or RNA data, the size of the dataset and the available computational resources. If sufficient resources are available, we recommend Optitype and HLA-HD for MHC-I and MHC-II genotype calling respectively.

MHC architecture in amphibians - ancestral reconstruction, gene rearrangements and duplication patterns.

The hypervariable major histocompatibility complex (MHC) is a crucial component of vertebrate adaptive immunity, but large-scale studies on MHC macroevolution in non-model vertebrates have long been constrained by methodological limitations. Here, we used rapidly accumulating genomic data to reconstruct macroevolution of the MHC region in amphibians. We retrieved contigs containing the MHC region from genome assemblies of 32 amphibian species and examined major structural rearrangements, duplication patterns and gene structure across the amphibian phylogeny. Based on the few available caecilian and urodele genomes we showed that the structure of ancestral MHC region in amphibians was probably relatively simple and compact, with a close physical linkage between MHC-I and MHC-II regions. This ancestral MHC architecture was generally conserved in anurans, although the evolution of class I subregion proceeded towards more extensive duplication and rapid expansion of gene copy number, providing evidence for dynamic evolutionary trajectories. Although in anurans we recorded tandems of duplicated MHC-I genes outside the core subregion, our phylogenetic analyses of MHC-I sequences provided little support for an expansion of nonclassical MHC-Ib genes across amphibian families. Finally, we found that intronic regions of amphibian classical MHC genes were much longer when compared to other tetrapod lineages (birds and mammals), which could partly be driven by the expansion of genome size. Our study reveals novel evolutionary patterns of the MHC region in amphibians and provides a comprehensive framework for further studies on the MHC macroevolution across vertebrates.

Genome-wide analyses reveal novel opioid use disorder loci and genetic overlap with schizophrenia, bipolar disorder, and major depression.

Opioid use disorder (OUD) and mental disorders are often comorbid, with increased morbidity and mortality. The causes underlying this relationship are poorly understood. Although these conditions are highly heritable, their shared genetic vulnerabilities remain unaccounted for. We applied the conditional/conjunctional false discovery rate (cond/conjFDR) approach to analyse summary statistics from independent genome wide association studies of OUD, schizophrenia (SCZ), bipolar disorder (BD) and major depression (MD) of European ancestry. Next, we characterized the identified shared loci using biological annotation resources. OUD data were obtained from the Million Veteran Program, Yale-Penn and Study of Addiction: Genetics and Environment (SAGE) (15 756 cases, 99 039 controls). SCZ (53 386 cases, 77 258 controls), BD (41 917 cases, 371 549 controls) and MD (170 756 cases, 329 443 controls) data were provided by the Psychiatric Genomics Consortium. We discovered genetic enrichment for OUD conditional on associations with SCZ, BD, MD and vice versa, indicating polygenic overlap with identification of 14 novel OUD loci at condFDR < 0.05 and 7 unique loci shared between OUD and SCZ (n = 2), BD (n = 2) and MD (n = 7) at conjFDR < 0.05 with concordant effect directions, in line with estimated positive genetic correlations. Two loci were novel for OUD, one for BD and one for MD. Three OUD risk loci were shared with more than one psychiatric disorder, at DRD2 on chromosome 11 (BD and MD), at FURIN on chromosome 15 (SCZ, BD and MD) and at the major histocompatibility complex region (SCZ and MD). Our findings provide new insights into the shared genetic architecture between OUD and SCZ, BD and MD, indicating a complex genetic relationship, suggesting overlapping neurobiological pathways.

Newly identified variability of the antigen binding site coding sequences of the equine major histocompatibility complex class I and class II genes.

The major histocompatibility complex (MHC) with its class I and II genes plays a crucial role in the immune response to pathogens by presenting oligopeptide antigens to various immune response effector cells. In order to counteract the vast variability of infectious agents, MHC class I and II genes usually retain high levels of SNPs mainly concentrated in the exons encoding the antigen binding sites. The aim of the study was to reveal new variability of selected MHC genes with a special focus on MHC class I physical haplotypes. Long-range NGS to was used to identify exon 2-exon 3 alleles in three genetically distinct horse breeds. A total of 116 allelic variants were found in the MHC class I genes Eqca-1, Eqca-2, Eqca-7 and Eqca-Ψ, 112 of which were novel. The MHC class II DRA locus was confirmed to comprise five exon 2 alleles, and no new sequences were observed. Additional variability in terms of 15 novel exon 2 alleles was identified in the DQA1 locus. Extensive overall variability across the entire MHC region was confirmed by an analysis of MHC-linked microsatellite loci. Both diversifying and purifying selection were detected within the MHC class I and II loci analyzed.

Exploring the causal role of the immune response to varicella-zoster virus on multiple traits: a phenome-wide Mendelian randomization study

BackgroundThe immune response to infections could be largely driven by the individual’s genes, especially in the major histocompatibility complex (MHC) region. Varicella-zoster virus (VZV) is a highly communicable pathogen. In addition to infection, the reactivations of VZV can be a potential causal factor for multiple traits. Identification of VZV immune response-related health conditions can therefore help elucidate the aetiology of certain diseases.MethodsA phenome-wide Mendelian randomization (MR) study of anti-VZV immunoglobulin G (IgG) levels with 1370 traits was conducted to explore the potential causal role of VZV-specific immune response on multiple traits using the UK Biobank cohort. For the robustness of the results, we performed MR analyses using five different methods. To investigate the impact of the MHC region on MR results, the analyses were conducted using instrumental variables (IVs) inside (IVmhc) and outside (IVno.mhc) the MHC region or all together (IVfull).ResultsForty-nine single nucleotide polymorphisms (IVfull) were associated with anti-VZV IgG levels, of which five (IVmhc) were located in the MHC region and 44 (IVno.mhc) were not. Statistical evidence (false discovery rate < 0.05 in at least three of the five MR methods) for a causal effect of anti-VZV IgG levels was found on 22 traits using IVmhc, while no evidence was found when using IVno.mhc or IVfull. The reactivations of VZV increased the risk of Dupuytren disease, mononeuropathies of the upper limb, sarcoidosis, coeliac disease, teeth problems and earlier onset of allergic rhinitis, which evidence was concordant with the literature. Suggestive causal evidence (P < 0.05 in at least three of five MR methods) using IVfull, IVmhc and IVno.mhc was detected in 92, 194 and 56 traits, respectively. MR results from IVfull correlated with those from IVmhc or IVno.mhc. However, the results between IVmhc and IVno.mhc were noticeably different, as evidenced by causal associations in opposite directions between anti-VZV IgG and ten traits.ConclusionsIn this exploratory study, anti-VZV IgG was causally associated with multiple traits. IVs in the MHC region might have a substantial impact on MR, and therefore, could be potentially considered in future MR studies.

Abstract 766: A genome-wide association study identifies the association of common variants in the major histocompatibility complex with cervical cancer in black South African women

Abstract Cervical cancer is the second most common female gynaecological malignancy in sub-Saharan Africa (SSA). Previous genome-wide association studies (GWASs) identified genetic variants associated with cervical cancer in women of European and Asian ancestry, particularly within the major histocompatibility complex (MHC). However, owing to the complexity of linkage disequilibrium and varied environmental exposures in different ancestral populations, regional GWASs are needed to investigate the genetic contribution to the pathogenesis of cervical cancer in continental African populations. A total of 1,637 cervical cancer patients and 761 controls in women of African ancestry were identified from the Johannesburg Cancer Study and the AWI-Gen study, respectively. DNA samples were genotyped on the 2.3 million H3Africa GWAS SNP array and genotyped SNPs imputed to a total of 13,937,308 SNPs. The association between genetic variants and cervical cancer was tested using linear mixed models, adjusting for five principal components and cryptic relatedness. We identified one locus within the MHC region associated with cervical cancer at genome-wide significance, led by rs3129965 (odds ratio [OR] = 0.89, 95% confidence interval [CI] = 0.87-0.92, P = 2.292 x 10-11) near HLA-DRA. Four other loci were strongly associated with cervical cancer (P&amp;lt;5x10-7), including rs575105445 (OR=0.77, 95%CI=0.70-0.85, P = 2.188 x 10-07) near HLA-C,rs184053148 (OR =0.79, 95% CI = 0.72-0.86, P = 4.19 x 10-07) in TRIM31,rs12496921 (OR=0.88, 95%CI = 0.84-0.92, P = 1.32 x 10-07) in ROPN1B and rs114031308 (OR= 1.27, 95% CI =1.16-1.39, P =3.06 x 10-07) near TTC26. The genes in these loci are implicated in numerous pathways including antiviral immunity, DNA repair, the regulation of p53, apoptosis, cell growth and differentiation. These discoveries demonstrate the important contribution of common genetic variants to cervical cancer susceptibility and provide insights into its pathogenesis in women of African ancestry. Citation Format: Abram Bunya Kamiza, Wenlong Carl Chen, Jean-Tristan Brandenberg, Michele Ramsay, Christopher G. Mathew. A genome-wide association study identifies the association of common variants in the major histocompatibility complex with cervical cancer in black South African women [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 766.

Complete sequences of six major histocompatibility complex haplotypes, including all the major MHC class II structures.

Accurate and comprehensive immunogenetic reference panels are key to the successful implementation of population-scale immunogenomics. The 5Mbp Major Histocompatibility Complex (MHC) is the most polymorphic region of the human genome and associated with multiple immune-mediated diseases, transplant matching and therapy responses. Analysis of MHC genetic variation is severely complicated by complex patterns of sequence variation, linkage disequilibrium and a lack of fully resolved MHC reference haplotypes, increasing the risk of spurious findings on analyzing this medically important region. Integrating Illumina, ultra-long Nanopore, and PacBio HiFi sequencing as well as bespoke bioinformatics, we completed five of the alternative MHC reference haplotypes of the current (GRCh38/hg38) build of the human reference genome and added one other. The six assembled MHC haplotypes encompass the DR1 and DR4 haplotype structures in addition to the previously completed DR2 and DR3, as well as six distinct classes of the structurally variable C4 region. Analysis of the assembled haplotypes showed that MHC class II sequence structures, including repeat element positions, are generally conserved within the DR haplotype supergroups, and that sequence diversity peaks in three regions around HLA-A, HLA-B+C, and the HLA class II genes. Demonstrating the potential for improved short-read analysis, the number of proper read pairs recruited to the MHC was found to be increased by 0.06%-0.49% in a 1000 Genomes Project read remapping experiment with seven diverse samples. Furthermore, the assembled haplotypes can serve as references for the community and provide the basis of a structurally accurate genotyping graph of the complete MHC region.

MICA and MICB allele assortment in Finland.

Genetic variation in the MICA and MICB genes located within the major histocompatibility complex region has been reported to be associated with transplantation outcome and susceptibility to autoimmune diseases and infections. Only limited data of polymorphism in these genes in different populations are available. We here report allelic variation at 2-field resolution and the haplotypes of the MICA and MICB genes in Finland (n= 1032 individuals), a north European population with historical bottleneck and founder effects. Altogether 24 MICA and 18 MICB alleles were found, forming 70 estimated MICA-MICB haplotypes. As compared to other populations frequency differences were found, for example, MICA*010:01 was found to be at an allele frequency of 0.133 in Finland which is higher than in other European populations (0.021-0.077), but close to Asian populations (0.151-0.220). Three novel alleles with amino acid change are described. The results demonstrate a relatively high level of polymorphism and population differences in MICA and MICB allele distribution.

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.

Polygenic risk score prediction of multiple sclerosis in individuals of South Asian ancestry.

Polygenic risk scores aggregate an individual's burden of risk alleles to estimate the overall genetic risk for a specific trait or disease. Polygenic risk scores derived from genome-wide association studies of European populations perform poorly for other ancestral groups. Given the potential for future clinical utility, underperformance of polygenic risk scores in South Asian populations has the potential to reinforce health inequalities. To determine whether European-derived polygenic risk scores underperform at multiple sclerosis prediction in a South Asian-ancestry population compared with a European-ancestry cohort, we used data from two longitudinal genetic cohort studies: Genes & Health (2015-present), a study of ∼50 000 British-Bangladeshi and British-Pakistani individuals, and UK Biobank (2006-present), which is comprised of ∼500 000 predominantly White British individuals. We compared individuals with and without multiple sclerosis in both studies (Genes & Health: N Cases = 42, N Control = 40 490; UK Biobank: N Cases = 2091, N Control = 374 866). Polygenic risk scores were calculated using clumping and thresholding with risk allele effect sizes obtained from the largest multiple sclerosis genome-wide association study to date. Scores were calculated with and without the major histocompatibility complex region, the most influential locus in determining multiple sclerosis risk. Polygenic risk score prediction was evaluated using Nagelkerke's pseudo-R 2 metric adjusted for case ascertainment, age, sex and the first four genetic principal components. We found that, as expected, European-derived polygenic risk scores perform poorly in the Genes & Health cohort, explaining 1.1% (including the major histocompatibility complex) and 1.5% (excluding the major histocompatibility complex) of disease risk. In contrast, multiple sclerosis polygenic risk scores explained 4.8% (including the major histocompatibility complex) and 2.8% (excluding the major histocompatibility complex) of disease risk in European-ancestry UK Biobank participants. These findings suggest that polygenic risk score prediction of multiple sclerosis based on European genome-wide association study results is less accurate in a South Asian population. Genetic studies of ancestrally diverse populations are required to ensure that polygenic risk scores can be useful across ancestries.

Identification of the BTN3A3 Gene as a Molecule Implicated in Generalized Pustular Psoriasis in a Chinese Population

The discovery of pathogenic variants provided biological insight into the role of host genetic factors in generalized pustular psoriasis (GPP). However, not all those affected by GPP carry variants in the reported genes. To comprehensively explore the molecular pathogenesis of GPP, whole-exome sequencing was performed, and two loci were identified with exome-wide significance through single variant association analysis: rs148755083 in the IL36RN gene (Pcombined= 1.19× 10-18, OR= 8.26) and HLA-C∗06:02 within the major histocompatibility complex region (Pcombined= 8.38× 10-12, OR= 2.98). Gene burden testing revealed that BTN3A3 correlated with GPP (Pcombined= 1.14× 10-10, OR= 5.59). Subtype analysis showed that IL36RN and BTN3A3 were both significantly associated with GPP alone and GPP with psoriasis vulgaris, whereas a correlation with HLA-C∗06:02 was only observed in GPP with psoriasis vulgaris. Functional analysis revealed that BTN3A3 regulated cell proliferation and inflammatory balance in GPP. In particular, loss of function of BTN3A3 activated NF-κB and promoted the production of inflammatory cytokines by inhibiting IL-36Ra expression to disturb the IL-1/IL-36 inflammatory axis and enhance the TNF-α-mediated pathway. Our findings identify BTN3A3 as, to our knowledge, a previously unreported pathogenic determinant, expanding our understanding of the genetic basis of GPP.

Exploring the overlap between alopecia areata and major depressive disorder: Epidemiological and genetic perspectives.

Research suggests that Alopecia areata (AA) and Major Depressive Disorder (MDD) show substantial comorbidity. To date, no study has investigated the hypothesis that this is attributable to shared genetic aetiology. To investigate AA-MDD comorbidity on the epidemiological and molecular genetic levels. First, epidemiological analyses were performed using data from a cohort of adult German health insurance beneficiaries (n=1.855million) to determine the population-based prevalence of AA-MDD comorbidity. Second, analyses were performed to determine the prevalence of MDD in a clinical AA case-control sample with data on psychiatric phenotypes, stratifying for demographic factors to identify possible contributing factors to AA-MDD comorbidity. Third, the genetic overlap between AA and MDD was investigated using a polygenic risk score (PRS) approach and linkage disequilibrium score (LDSC) regression. For PRS, summary statistics from a large MDD GWAS meta-analysis (PGC-MD2) were used as the training sample, while a Central European AA cohort, including the above-mentioned AA patients, and an independent replication US-AA cohort were used as target samples. LDSC was performed using summary statistics of PGC-MD2 and the largest AA meta-analysis to date. High levels of AA-MDD comorbidity were reported in the population-based (MDD in 24% of AA patients), and clinical samples (MDD in 44% of AA patients). MDD-PRS explained a modest proportion of variance in AA case-control status (R2 =1%). This signal was limited to the major histocompatibility complex (MHC) region on chromosome 6. LDSC regression (excluding MHC) revealed no significant genetic correlation between AA and MDD. As in previous research, AA patients showed an increased prevalence of MDD. The present analyses suggest that genetic overlap may be confined to the MHC region, which is implicated in immune function. More detailed investigation is required to refine understanding of how the MHC is involved in the development of AA and MDD comorbidity.

LncRNA HCG11 suppresses cell proliferation in hormone receptor-positive breast cancer via SRSF1/β-catenin.

Hormone receptor positive (HR-positive) breast cancer (BC) is the most common subtype of breast cancer. Despite adjuvant endocrine therapy and chemotherapy-based treatment, the therapeutic response is often not satisfactory in HR-positive BC patients. Therefore, elucidating the mechanisms that regulate the progression of HR-positive BC is urgently required to identify new therapeutic targets. Previously, HLA Complex Group 11 (HCG11), located on the major histocompatibility complex (MHC) region, was found to be abnormally expressed in a variety of tumor cells. However, the role of HCG11 in HR-positive BC cells has not been explored to date. In the current study, we found that HCG11 is downregulated in HR-positive BC tissues and cell lines. Both in vitro and in vivo, HCG11 acts as a tumor suppressor in HR-positive BC cells. Furthermore, the mechanistic details unraveled that HCG11 recruits Serine/arginine-rich splicing factor 1 (SRSF1) to target β-catenin mRNA for promoting the translation of β-catenin. Our study emphasizes the potential of HCG11 as a novel intervention target for HR-positive BC treatment.