Population-specific recombination sites within the human 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.

A cost-effective strategy for selection of third-party donors for a virus-specific T-cell bank for an Asian patient population

BACKGROUND AIMS

Third party virus-specific T cells (VST) has shown efficacy for opportunistic virus infection which do not have effective treatment or are drug-refractory. We describe our preparatory work in setting up a third-party VST bank for a multi-ethnic Asian population.

METHODS

Discarded white cells from regular blood bank plateletpheresis donors with known locally prevalent HLA antigens were cultured in small scale to generate VST against Adenovirus, BK virus, Cytomegalovirus, Epstein-Barr virus, and Human Herpes Virus 6. Multi-virus specific T cells (multi-VST) were also generated against all 5 viruses in single cultures. A strategy of allelic typing for donors with good and broad-spectrum cytotoxicity together with consideration on HLA restriction for the virus epitope was used to select combinations of VST lines for a hypothetical third party VST bank. The breadth of coverage based on these selection criteria was validated using our database of 100 post haematopoietic stem cell transplant patients.

RESULTS

We show that 50%, 42%, 56%, 56% and 42% of single VST cultures demonstrated specific cytotoxicity against AdV, BKV, CMV, EBV and HHV6 respectively. Twenty four of the 36 multi-VST lines showed activity against at least 2 of the 5 viruses studied. A carefully selected combination of just 6 VST lines can offer VST with at least 1 allelic match to 99% of potential recipients, while 92% can find 2 allelic matches and 79% can find 3 allelic matches.

CONCLUSIONS

This preparatory work confirms that a cost-effective strategy recruiting a small number of pre-characterized donors can generate VST lines with broad coverage for a multi-ethnic Asian patient population, thereby laying the foundation for setting up of a third party VST bank for Asian patients.

Recombination, selection and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.

A Comprehensive Review of HLA and Severe Cutaneous Adverse Drug Reactions: Implication for Clinical Pharmacogenomics and Precision Medicine

Human leukocyte antigen (HLA) encoded by the HLA gene is an important modulator for immune responses and drug hypersensitivity reactions as well. Genetic polymorphisms of HLA vary widely at population level and are responsible for developing severe cutaneous adverse drug reactions (SCARs) such as Stevens-Johnson syndrome (SJS), toxic epidermal necrolysis (TEN), drug reaction with eosinophilia and systemic symptoms (DRESS), maculopapular exanthema (MPE). The associations of different HLA alleles with the risk of drug induced SJS/TEN, DRESS and MPE are strongly supportive for clinical considerations. Prescribing guidelines generated by different national and international working groups for translation of HLA pharmacogenetics into clinical practice are underway and functional in many countries, including Thailand. Cutting edge genomic technologies may accelerate wider adoption of HLA screening in routine clinical settings. There are great opportunities and several challenges as well for effective implementation of HLA genotyping globally in routine clinical practice for the prevention of drug induced SCARs substantially, enforcing precision medicine initiatives.

TUBB8 Mutations Cause Female Infertility with Large Polar Body Oocyte and Fertilization Failure

Tubulin beta 8 class VIII (TUBB8) is a special β-tubulin isotype that mainly expressed in primate oocytes and early embryos and identified as the disease-causing gene of human oocyte maturation arrest. To identify the disease-causing genes in 2 patients with female infertility due to large polar body oocyte or fertilization failure, whole-exome sequencing was performed on the patients and available family members. We identified a novel heterozygous missense mutation c.817C>G (p.L273V) and a recently reported heterozygous missense mutation c.608A>G (p.D203G) in TUBB8 from two patients, respectively. We found oocyte with a large polar body in the patient who carried the p.D203G mutation in TUBB8. Bioinformatics analysis showed that these two mutations are harmful. The results of western blot and RT-PCR experiments showed that the D203G mutation caused a significant decrease in the expression of TUBB8, and immunostaining showed that the TUBB8 mutation caused abnormal microtubule morphology. These findings suggest that TUBB8 mutations resulted in oocyte with a large polar body and fertilization failure in patients.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

The changing landscape of HLA typing: Understanding how and when HLA typing data can be used with confidence from bench to bedside

Human leukocyte antigen (HLA) genes are extraordinary for their extreme diversity and widespread impact on human health and disease. More than 30,000 HLA alleles have been officially named and more alleles continue to be discovered at a rapid pace. HLA typing systems which have been developed to detect HLA diversity have advanced rapidly and are revolutionizing our understanding of HLA's clinical importance. However, continuous improvements in knowledge and technology have created challenges for clinicians and scientists. This review explains how differences in HLA typing systems can impact the HLA types that are assigned. The consequences of differences in laboratory testing methods and reference databases are described. The challenges of using HLA types that are not equivalent are illustrated. A fundamental understanding of the continual expansion of our understanding of HLA diversity and limitations in some of the typing data is essential for using typing data appropriately in clinical and research settings.

HLA-DQB1 6672G>C (rs113332494) is associated with clozapine-induced neutropenia and agranulocytosis in individuals of European ancestry

The atypical antipsychotic clozapine is the only effective medication for treatment-resistant schizophrenia. However, it can also induce serious adverse drug reactions, including agranulocytosis and neutropenia. The mechanism by which it does so is largely unknown, but there is evidence for contributing genetic factors. Several studies identified HLA-DQB1 variants and especially a polymorphism located in HLA-DQB1 (6672G>C, rs113332494) as associated with clozapine-induced agranulocytosis and neutropenia. We analysed the risk allele distribution of SNP rs113332494 in a sample of 1396 controls and 178 neutropenia cases of which 60 developed agranulocytosis. Absolute neutrophil counts of 500/mm³ and 1500/mm³ were used for defining agranulocytosis and neutropenia cases, respectively. We also performed association analyses and analysed local ancestry patterns in individuals of European ancestry, seeking replication and extension of earlier findings. HLA-DQB1 (6672G>C, rs113332494) was associated with neutropenia (OR = 6.20, P = 2.20E-06) and agranulocytosis (OR = 10.49, P = 1.83E-06) in individuals of European ancestry. The association signal strengthened after including local ancestry estimates (neutropenia: OR = 10.38, P = 6.05E-08; agranulocytosis: OR = 16.31, P = 1.39E-06), with effect sizes being considerably larger for agranulocytosis. Using local ancestry estimates for prediction, the sensitivity of rs113332494 increased from 11.28 to 55.64% for neutropenia and from 16.67 to 53.70% for agranulocytosis. Our study further strengthens the evidence implicating HLA-DQB1 in agranulocytosis and neutropenia, suggesting components of the immune system as contributing to this serious adverse drug reaction. Using local ancestry estimates might help in identifying risk variants and improve prediction of haematological adverse effects.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

Discovery of runs-of-homozygosity diplotype clusters and their associations with diseases in UK Biobank

Runs of homozygosity (ROH) segments, contiguous homozygous regions in a genome were traditionally linked to families and inbred populations. However, a growing literature suggests that ROHs are ubiquitous in outbred populations. Still, most existing genetic studies of ROH in populations are limited to aggregated ROH content across the genome, which does not offer the resolution for mapping causal loci. This limitation is mainly due to a lack of methods for efficient identification of shared ROH diplotypes. Here, we present a new method, ROH-DICE, to find large ROH diplotype clusters, sufficiently long ROHs shared by a sufficient number of individuals, in large cohorts. ROH-DICE identified over 1 million ROH diplotypes that span over 100 SNPs and shared by more than 100 UK Biobank participants. Moreover, we found significant associations of clustered ROH diplotypes across the genome with various self-reported diseases, with the strongest associations found between the extended HLA region and autoimmune disorders. We found an association between a diplotype covering the HFE gene and haemochromatosis, even though the well-known causal SNP was not directly genotyped nor imputed. Using genome-wide scan, we identified a putative association between carriers of an ROH diplotype in chromosome 4 and an increase of mortality among COVID-19 patients. In summary, our ROH-DICE method, by calling out large ROH diplotypes in a large outbred population, enables further population genetics into the demographic history of large populations. More importantly, our method enables a new genome-wide mapping approach for finding disease-causing loci with multi-marker recessive effects at population scale.

The Genetic Mechanisms Driving Diversification of the KIR Gene Cluster in Primates

The activity and function of natural killer (NK) cells are modulated through the interactions of multiple receptor families, of which some recognize MHC class I molecules. The high level of MHC class I polymorphism requires their ligands either to interact with conserved epitopes, as is utilized by the NKG2A receptor family, or to co-evolve with the MHC class I allelic variation, which task is taken up by the killer cell immunoglobulin-like receptor (KIR) family. Multiple molecular mechanisms are responsible for the diversification of the KIR gene system, and include abundant chromosomal recombination, high mutation rates, alternative splicing, and variegated expression. The combination of these genetic mechanisms generates a compound array of diversity as is reflected by the contraction and expansion of KIR haplotypes, frequent birth of fusion genes, allelic polymorphism, structurally distinct isoforms, and variegated expression, which is in contrast to the mainly allelic nature of MHC class I polymorphism in humans. A comparison of the thoroughly studied human and macaque KIR gene repertoires demonstrates a similar evolutionarily conserved toolbox, through which selective forces drove and maintained the diversified nature of the KIR gene cluster. This hypothesis is further supported by the comparative genetics of KIR haplotypes and genes in other primate species. The complex nature of the KIR gene system has an impact upon the education, activity, and function of NK cells in coherence with an individual’s MHC class I repertoire and pathogenic encounters. Although selection operates on an individual, the continuous diversification of the KIR gene system in primates might protect populations against evolving pathogens.

HLA Correlations with Clinical Phenotypes and Risk of Metabolic Comorbidities in Singapore Chinese Psoriasis Patients

INTRODUCTION

Psoriasis is a systemic, chronic inflammatory disease that not only afflicts the skin but is also associated with cardiovascular disease and metabolic syndrome. The strongest susceptibility loci for the disease is within the human leukocyte antigen (HLA) complex, though specific HLA allelic associations vary between populations.

OBJECTIVE

Our objective was to investigate HLA associations with clinical phenotypes of psoriasis and metabolic syndrome in Chinese psoriasis cases.

METHODS

We conducted an observational case-control study in Singapore with a cohort of psoriasis cases consecutively recruited from an outpatient specialist dermatological center (n = 120) compared with 130 healthy controls.

RESULTS

Significant HLA associations with psoriasis were observed with HLA-A*02:07, B*46:01, C*01:02, and C*06:02. The three-locus haplotype of A*02:07-C*01:02-B*46:01 was also significant (odds ratio [OR] 3.07; p = 9.47 × 10^-5). We also observed an association between nail psoriasis and HLA-A*02:07 carriers (OR 4.50; p = 0.002), whereas C*06:02 carriers were less prone to have nail involvement (OR 0.16; p = 0.004). HLA-A*02:07 was also identified as a possible risk allele for hypertension (OR 2.90; p < 0.05), and C*01:02 was a possible risk allele for dyslipidemia (OR 3.36; p < 0.05), both known to be common comorbidities in patients with psoriasis.

CONCLUSION

Our results demonstrate the growing importance of discerning population-specific clinical phenotypes and their association with certain HLA alleles in psoriasis.

Single haplotype admixture models using large scale HLA genotype frequencies to reproduce human admixture

The human leukocyte antigen (HLA) is the most polymorphic region in humans. Anthropologists use HLA to trace populations' migration and evolution. However, recent admixture between populations can mask the ancestral haplotype frequency distribution. We present a statistical method based on high-resolution HLA haplotype frequencies to resolve population admixture using a non-negative matrix factorization formalism and validated using haplotype frequencies from 56 world populations. The result is a minimal set of source components (SCs) decoding roughly 90% of the total variance in the studied admixtures. These SCs agree with the geographical distribution, phylogenies, and recent admixture events of the studied groups. With the growing population of multi-ethnic individuals, or individuals that do not report race/ethnic information, the HLA matching process for stem-cell and solid organ transplants is becoming more challenging. The presented algorithm provides a framework that facilitates the breakdown of highly admixed populations into SCs, which can be used to better match the rapidly growing population of multi-ethnic individuals worldwide.

MHC haplotype diversity in Icelandic horses determined by polymorphic microsatellites

The Icelandic horse has been maintained as a closed population in its eponymous homeland for many generations, with no recorded introductions of new horses of any breed since the year 1000 CE. Here we determined the diversity of major histocompatibility complex (MHC) haplotypes in 156 Icelandic horses from two groups, based on a panel of 12 polymorphic intra-MHC microsatellites tested in families of various composition. We identified a total of 79 MHC haplotypes in these two groups, including one documented intra-MHC recombination event from a total of 147 observed meioses. None of these MHC haplotypes have been previously described in any other horse breed. Only one MHC homozygote was found in the entire population studied. These results indicate a very high level of MHC heterozygosity and haplotype diversity in the Icelandic horse. The environment in Iceland is remarkable for its lack of common agents of equine infectious disease, including equine herpesvirus type 1, influenza virus, and streptococcus equi. The driving forces for maintenance of MHC heterozygosity in Icelandic horses must thus be sought outside of these major horse pathogens. Based on our results, we propose that intra-MHC recombination may play a major role in the generation of novel haplotypes.

LDSplitDB: a database for studies of meiotic recombination hotspots in MHC using human genomic data

Background

Meiotic recombination happens during the process of meiosis when chromosomes inherited from two parents exchange genetic materials to generate chromosomes in the gamete cells. The recombination events tend to occur in narrow genomic regions called recombination hotspots. Its dysregulation could lead to serious human diseases such as birth defects. Although the regulatory mechanism of recombination events is still unclear, DNA sequence polymorphisms have been found to play crucial roles in the regulation of recombination hotspots.

Method

To facilitate the studies of the underlying mechanism, we developed a database named LDSplitDB which provides an integrative and interactive data mining and visualization platform for the genome-wide association studies of recombination hotspots. It contains the pre-computed association maps of the major histocompatibility complex (MHC) region in the 1000 Genomes Project and the HapMap Phase III datasets, and a genome-scale study of the European population from the HapMap Phase II dataset. Besides the recombination profiles, related data of genes, SNPs and different types of epigenetic modifications, which could be associated with meiotic recombination, are provided for comprehensive analysis. To meet the computational requirement of the rapidly increasing population genomics data, we prepared a lookup table of 400 haplotypes for recombination rate estimation using the well-known LDhat algorithm which includes all possible two-locus haplotype configurations.

Conclusion

To the best of our knowledge, LDSplitDB is the first large-scale database for the association analysis of human recombination hotspots with DNA sequence polymorphisms. It provides valuable resources for the discovery of the mechanism of meiotic recombination hotspots. The information about MHC in this database could help understand the roles of recombination in human immune system.

Database URL

http://histone.scse.ntu.edu.sg/LDSplitDB

Hidden genomic MHC disparity between HLA-matched sibling pairs in hematopoietic stem cell transplantation

Matching classical HLA alleles between donor and recipient is an important factor in avoiding adverse immunological effects in HSCT. Siblings with no differences in HLA alleles, either due to identical-by-state or identical-by-descent status, are considered to be optimal donors. We carried out a retrospective genomic sequence and SNP analysis of 336 fully HLA-A, -B, -DRB1 matched and 14 partially HLA-matched sibling HSCT pairs to determine the level of undetected mismatching within the MHC segment as well as to map their recombination sites. The genomic sequence of 34 genes locating in the MHC region revealed allelic mismatching at 1 to 8 additional genes in partially HLA-matched pairs. Also, fully matched pairs were found to have mismatching either at HLA-DPB1 or at non-HLA region within the MHC segment. Altogether, 3.9% of fully HLA-matched HSCT pairs had large genomic mismatching in the MHC segment. Recombination sites mapped to certain restricted locations. The number of mismatched nucleotides correlated with the risk of GvHD supporting the central role of full HLA matching in HSCT. High-density genome analysis revealed that fully HLA-matched siblings may not have identical MHC segments and even single allelic mismatching at any classical HLA gene often implies larger genomic differences along MHC.

Tailoring of recommendations to reduce serious cutaneous adverse drug reactions: a pharmacogenomics approach

The Health Sciences Authority launched a pharmacogenetics initiative in 2008 to facilitate evaluation of pharmacogenetics associations pertinent for Chinese, Malays and Indians in Singapore. The aim was to reduce the incidence and unpredictability of serious adverse drug reactions, with a focus on serious skin adverse drug reactions. This paper describes the gathering of evidence and weighing of factors that led to different genotyping recommendations for HLA-B*15:02 with carbamazepine and HLA-B*58:01 with allopurinol, despite both having strong genetic associations. Translation of pharmacogenomics at a national level requires careful deliberation of the prevalence of at-risk allele, strength of genetic associations, positive predictive value, cost-effectiveness and availability of alternative therapies. Our experience provides a perspective on translating genomic discoveries in advancing drug safety.

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.

Facile and phase-defined determination of HLA alleles with morpholino-functionalized nanoparticle probes

A number of human leukocyte antigen (HLA) gene alleles have been found to be genetic risk markers for immunologically mediated drug hypersensitivity. Clinical adoption of HLA pharmacogenomics requires facile and accurate allele screening assays. As HLA genes are highly polymorphic, currently available methods are usually labor-intensive and liable to generate false positives. Herein we report a general strategy for screening HLA alleles with nanoparticle probes. Specific HLA alleles can be identified by gauging three to five sequence variants. Single-polymerase chain reaction (PCR) and dual-PCR methods have been proposed to achieve phase-defined determination of the sequence variants. Morpholino-functionalized gold nanoparticle probes allow for colorimetric and highly specific detection. Assays for HLA-B*58:01 and HLA-B*15:02 have been developed and validated with 49 selected human genomic DNA samples. The facile nanoparticle probe-based assays can be implemented easily in molecular diagnostic laboratories for accurate and cost-effective screening of HLA alleles.

Conserved 33-kb haplotype in the MHC class III region regulates chronic arthritis

Genome-wide association studies have revealed many genetic loci associated with complex autoimmune diseases. In rheumatoid arthritis (RA), the MHC gene HLA-DRB1 is the strongest candidate predicting disease development. It has been suggested that other immune-regulating genes in the MHC contribute to the disease risk, but this contribution has been difficult to show because of the strong linkage disequilibrium within the MHC. We isolated genomic regions in the form of congenic fragments in rats to test whether there are additional susceptibility loci in the MHC. By both congenic mapping in inbred strains and SNP typing in wild rats, we identified a conserved, 33-kb large haplotype Ltab-Ncr3 in the MHC-III region, which regulates the onset, severity, and chronicity of arthritis. The Ltab-Ncr3 haplotype consists of five polymorphic immunoregulatory genes: Lta (lymphotoxin-α), Tnf, Ltb (lymphotoxin-β), Lst1 (leukocyte-specific transcript 1), and Ncr3 (natural cytotoxicity-triggering receptor 3). Significant correlation in the expression of the Ltab-Ncr3 genes suggests that interaction of these genes may be important in keeping these genes clustered together as a conserved haplotype. We studied the arthritis association and the spliceo-transcriptome of four different Ltab-Ncr3 haplotypes and showed that higher Ltb and Ncr3 expression, lower Lst1 expression, and the expression of a shorter splice variant of Lst1 correlate with reduced arthritis severity in rats. Interestingly, patients with mild RA also showed higher NCR3 expression and lower LST1 expression than patients with severe RA. These data demonstrate the importance of a conserved haplotype in the regulation of complex diseases such as arthritis.

Actinic Prurigo in Singaporean Chinese: A Positive Association with HLA-DRB1*03:01

Studies have reported the association of human leukocyte antigen (HLA) genes with susceptibility to develop actinic prurigo (AP) in Caucasians, but there were no studies in Asian populations, including the Chinese. Our study was performed to determine if AP is associated with susceptibility or protective HLA alleles or haplotypes in Singaporean Chinese. All Chinese patients diagnosed with AP at National Skin Center, Singapore, from January 2002 to April 2015 were invited to participate in the study. Clinical data and phototesting results were collated, and HLA typing was performed. Among 14 patients included, 11 were male and the mean age was 49.6 (37.9-61.3) years. All patients did not have a family history of AP and none had mucosal involvement, as such these clinical features differed from Caucasian AP patients. The frequency of DRB1*03:01 in AP patients was significantly higher compared to healthy controls (43% vs 16%, P = 0.022, odds ratio (OR) 3.89). Concurrently, the frequency of HLA-B*58:01-DRB1*03:01 haplotype was also significantly increased (25% vs 7%, P = 0.004, OR 4.23). In conclusion, HLA-DRB1*03:01 was associated with AP in Singaporean Chinese patients. This novel allelic association may possibly be utilized as a biological marker to aid in the diagnosis of AP in Chinese patients.

Molecular characterization of classical and nonclassical MHC class I genes from the golden pheasant (Chrysolophus pictus)

Classical major histocompatibility complex (MHC) class I allelic polymorphism is essential for competent antigen presentation. To improve the genotyping efforts in the golden pheasant, it is necessary to differentiate more accurately between classical and nonclassical class I molecules. In our study, all MHC class I genes were isolated from one golden pheasant based on two overlapping PCR amplifications. In total, six full-length class I nucleotide sequences (A-F) were identified, and four were novel. Two (A and C) belonged to the IA1 gene, two (B and D) were alleles derived from the IA2 gene through transgene amplification, and two (E and F) comprised a third novel locus, IA3 that was excluded from the core region of the golden pheasant MHC-B. IA1 and IA2 exhibited the broad expression profiles characteristic of classical loci, while IA3 showed no expression in multiple tissues and was therefore defined as a nonclassical gene. Phylogenetic analysis indicated that the three IA genes in the golden pheasant share a much closer evolutionary relationship than the corresponding sequences in other galliform species. This observation was consistent with high sequence similarity among them, which likely arises from the homogenizing effect of recombination. Our careful distinction between the classical and nonclassical MHC class I genes in the golden pheasant lays the foundation for developing locus-specific genotyping and establishing a good molecular marker system of classical MHC I loci.

Intrahaplotypic Variants Differentiate Complex Linkage Disequilibrium within Human MHC Haplotypes

Distinct regions of long-range genetic fixation in the human MHC region, known as conserved extended haplotypes (CEHs), possess unique genomic characteristics and are strongly associated with numerous diseases. While CEHs appear to be homogeneous by SNP analysis, the nature of fine variations within their genomic structure is unknown. Using multiple, MHC-homozygous cell lines, we demonstrate extensive sequence conservation in two common Asian MHC haplotypes: A33-B58-DR3 and A2-B46-DR9. However, characterization of phase-resolved MHC haplotypes revealed unique intra-CEH patterns of variation and uncovered 127 single nucleotide variants (SNVs) which are missing from public databases. We further show that the strong linkage disequilibrium structure within the human MHC that typically confounds precise identification of genetic features can be resolved using intra-CEH variants, as evidenced by rs3129063 and rs448489, which affect expression of ZFP57, a gene important in methylation and epigenetic regulation. This study demonstrates an improved strategy that can be used towards genetic dissection of diseases.

Tumor Necrosis Factor (TNF) -308G>A, Nitric Oxide Synthase 3 (NOS3) +894G>T Polymorphisms and Migraine Risk: A Meta-Analysis

Background and Objective

Conflicting data have been reported on the association between tumor necrosis factor (TNF) –308G>A and nitric oxide synthase 3 (NOS3) +894G>T polymorphisms and migraine. We performed a meta-analysis of case-control studies to evaluate whether the TNF –308G>A and NOS3 +894G>T polymorphisms confer genetic susceptibility to migraine.

Method

We performed an updated meta-analysis for TNF –308G>A and a meta-analysis for NOS3 +894G>T based on studies published up to July 2014. We calculated study specific odds ratios (OR) and 95% confidence intervals (95% CI) assuming allele contrast, dominant model, recessive model, and co-dominant model as pooled effect estimates.

Results

Eleven studies in 6682 migraineurs and 22591 controls for TNF –308G>A and six studies in 1055 migraineurs and 877 controls for NOS3 +894G>T were included in the analysis. Neither indicated overall associations between gene polymorphisms and migraine risk. Subgroup analyses suggested that the “A” allele of the TNF –308G>A variant increases the risk of migraine among non-Caucasians (dominant model: pooled OR = 1.82; 95% CI 1.15 – 2.87). The risk of migraine with aura (MA) was increased among both Caucasians and non-Caucasians. Subgroup analyses suggested that the “T” allele of the NOS3 +894G>T variant increases the risk of migraine among non-Caucasians (co-dominant model: pooled OR = 2.10; 95% CI 1.14 – 3.88).

Conclusions

Our findings appear to support the hypothesis that the TNF –308G>A polymorphism may act as a genetic susceptibility factor for migraine among non-Caucasians and that the NOS3 +894G>T polymorphism may modulate the risk of migraine among non-Caucasians.

Human leukocyte antigen typing using a knowledge base coupled with a high-throughput oligonucleotide probe array analysis

Human leukocyte antigens (HLA) are important biomarkers because multiple diseases, drug toxicity, and vaccine responses reveal strong HLA associations. Current clinical HLA typing is an elimination process requiring serial testing. We present an alternative in situ synthesized DNA-based microarray method that contains hundreds of thousands of probes representing a complete overlapping set covering 1,610 clinically relevant HLA class I alleles accompanied by computational tools for assigning HLA type to 4-digit resolution. Our proof-of-concept experiment included 21 blood samples, 18 cell lines, and multiple controls. The method is accurate, robust, and amenable to automation. Typing errors were restricted to homozygous samples or those with very closely related alleles from the same locus, but readily resolved by targeted DNA sequencing validation of flagged samples. High-throughput HLA typing technologies that are effective, yet inexpensive, can be used to analyze the world's populations, benefiting both global public health and personalized health care.

The immunodominant influenza A virus M158-66 cytotoxic T lymphocyte epitope exhibits degenerate class I major histocompatibility complex restriction in humans

Cytotoxic T lymphocytes recognizing conserved peptide epitopes are crucial for protection against influenza A virus (IAV) infection. The CD8 T cell response against the M158-66 (GILGFVFTL) matrix protein epitope is immunodominant when restricted by HLA-A*02, a major histocompatibility complex (MHC) molecule expressed by approximately half of the human population. Here we report that the GILGFVFTL peptide is restricted by multiple HLA-C*08 alleles as well. We observed that M158-66 was able to elicit cytotoxic T lymphocyte (CTL) responses in both HLA-A*02- and HLA-C*08-positive individuals and that GILGFVFTL-specific CTLs in individuals expressing both restriction elements were distinct and not cross-reactive. The crystal structure of GILGFVFTL-HLA-C*08:01 was solved at 1.84 Å, and comparison with the known GILGFVFTL-HLA-A*02:01 structure revealed that the antigen bound both complexes in near-identical conformations, accommodated by binding pockets shaped from shared as well as unique residues. This discovery of degenerate peptide presentation by both HLA-A and HLA-C allelic variants eliciting unique CTL responses to IAV infection contributes fundamental knowledge with important implications for vaccine development strategies.

IMPORTANCE

The presentation of influenza A virus peptides to elicit immunity is thought to be narrowly restricted, with a single peptide presented by a specific HLA molecule. In this study, we show that the same influenza A virus peptide can be more broadly presented by both HLA-A and HLA-C molecules. This discovery may help to explain the differences in immunity to influenza A virus between individuals and populations and may also aid in the design of vaccines.

Differential recombination dynamics within the MHC of macaque species

A panel of 15 carefully selected microsatellites (short tandem repeats, STRs) has allowed us to study segregation and haplotype stability in various macaque species. The STRs span the major histocompatibility complex (MHC) region and map in more detail from the centromeric part of the Mhc-A to the DR region. Two large panels of Indian rhesus and Indonesian/Indochinese cynomolgus macaques have been subjected to pedigree analysis, allowing the definition of 161 and 36 different haplotypes and the physical mapping of 10 and 5 recombination sites, respectively. Although most recombination sites within the studied section of the Indian rhesus monkey MHC are situated between the Mhc-A and Mhc-B regions, the resulting recombination rate for this genomic segment is low and similar to that in humans. In contrast, in Indonesian/Indochinese macaques, two recombination sites, which appear to be absent in rhesus macaques, map between the class III and II regions. As a result, the mean recombination frequency of the core MHC, Mhc-A to class II, is higher in Indonesian/Indochinese cynomolgus than in Indian rhesus macaques, but as such is comparable to that in humans. The present communication demonstrates that the dynamics of recombination 'hot/cold spots' in the MHC, as well as their frequencies, may differ substantially between highly related macaque species.

Intertypic modular exchanges of genomic segments by homologous recombination at universally conserved segments in human adenovirus species D

Human adenovirus species D (HAdV-D), which is composed of clinically and epidemiologically important pathogens worldwide, contains more taxonomic "types" than any other species of the genus Mastadenovirus, although the mechanisms accounting for the high level of diversity remain to be disclosed. Recent studies of known and new types of HAdV-D have indicated that intertypic recombination between distant types contributes to the increasing diversity of the species. However, such findings raise the question as to how homologous recombination events occur between diversified types since homologous recombination is suppressed as nucleotide sequences diverge. In order to address this question, we investigated the distribution of the recombination boundaries in comparison with the landscape of intergenomic sequence conservation assessed according to the synonymous substitution rate (dS). The results revealed that specific genomic segments are conserved between even the most distantly related genomes; we call these segments "universally conserved segments" (UCSs). These findings suggest that UCSs facilitate homologous recombination, resulting in intergenomic segmental exchanges of UCS-flanking genomic regions as recombination modules. With the aid of such a mechanism, the haploid genomes of HAdV-Ds may have been reshuffled, resulting in chimeric genomes out of diversified repertoires in the HAdV-D population analogous to the MHC region reshuffled via crossing over in vertebrates. In addition, some HAdVs with chimeric genomes may have had the opportunity to avoid host immune responses thereby causing epidemics.

Building pharmacogenetics into a pharmacovigilance program in Singapore: using serious skin rash as a pilot study

To study the possible genetic associations with adverse drug reactions (ADR), the Singapore Health Sciences Authority (HSA) has piloted a program to collect DNA and phenotype data of ADR cases as part of its pharmacovigilance program. Between 2009 and 2012, HSA screened 158 cases of Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN). To assess the association between HLA-B*1502 and carbamazepine (CBZ)-induced SJS/TEN, 13 cases and 26 drug-tolerant controls were analyzed. All 13 CBZ-SJS/TEN cases and 3/26 controls were HLA-B*1502 positive (odds ratio 181, 95% confidence interval: 8.7-3785, P=6.9 × 10(-8)). Discussions of the finding with the Ministry of Health and an expert panel led to the decision to make HLA-B*1502 testing the standard of care prior to first use of CBZ in Asians and to subsidize the genotyping test at public hospitals. This program illustrates the role of a regulatory authority in advancing the use of pharmacogenetics for drug safety.

Diversity of extended HLA-DRB1 haplotypes in the Finnish population

The Major Histocompatibility Complex (MHC, 6p21) codes for traditional HLA and other host response related genes. The polymorphic HLA-DRB1 gene in MHC Class II has been associated with several complex diseases. In this study we focus on MHC haplotype structures in the Finnish population. We explore the variability of extended HLA-DRB1 haplotypes in relation to the other traditional HLA genes and a selected group of MHC class III genes. A total of 150 healthy Finnish individuals were included in the study. Subjects were genotyped for HLA alleles (HLA-A, -B, -DRB1, -DQB1, and -DPB1). The polymorphism of TNF, LTA, C4, BTNL2 and HLA-DRA genes was studied with 74 SNPs (single nucleotide polymorphism). The C4A and C4B gene copy numbers and a 2-bp silencing insertion at exon 29 in C4A gene were analysed with quantitative genomic realtime-PCR. The allele frequencies for each locus were calculated and haplotypes were constructed using both the traditional HLA alleles and SNP blocks. The most frequent Finnish A∼B∼DR -haplotype, uncommon in elsewhere in Europe, was A*03∼B*35∼DRB1*01∶01. The second most common haplotype was a common European ancestral haplotype AH 8.1 (A*01∼B*08∼DRB1*03∶01). Extended haplotypes containing HLA-B, TNF block, C4 and HLA-DPB1 strongly increased the number of HLA-DRB1 haplotypes showing variability in the extended HLA-DRB1 haplotype structures. On the contrary, BTNL2 block and HLA-DQB1 were more conserved showing linkage with the HLA-DRB1 alleles. We show that the use of HLA-DRB1 haplotypes rather than single HLA-DRB1 alleles is advantageous when studying the polymorphisms and LD patters of the MHC region. For disease association studies the HLA-DRB1 haplotypes with various MHC markers allows us to cluster haplotypes with functionally important gene variants such as C4 deficiency and cytokines TNF and LTA, and provides hypotheses for further assessment. Our study corroborates the importance of studying population-specific MHC haplotypes.