An Expanded Role for HLA Genes: HLA-B Encodes a microRNA that Regulates IgA and Other Immune Response Transcripts

Phenome-wide association study on miRNA-related sequence variants: the UK Biobank

BACKGROUND

Genetic variants in the coding region could directly affect the structure and expression levels of genes and proteins. However, the importance of variants in the non-coding region, such as microRNAs (miRNAs), remain to be elucidated. Genetic variants in miRNA-related sequences could affect their biogenesis or functionality and ultimately affect disease risk. Yet, their implications and pleiotropic effects on many clinical conditions remain unknown.

METHODS

Here, we utilised genotyping and hospital records data in the UK Biobank (N = 423,419) to investigate associations between 346 genetic variants in miRNA-related sequences and a wide range of clinical diagnoses through phenome-wide association studies. Further, we tested whether changes in blood miRNA expression levels could affect disease risk through colocalisation and Mendelian randomisation analysis.

RESULTS

We identified 122 associations for six variants in the seed region of miRNAs, nine variants in the mature region of miRNAs, and 27 variants in the precursor miRNAs. These included associations with hypertension, dyslipidaemia, immune-related disorders, and others. Nineteen miRNAs were associated with multiple diagnoses, with six of them associated with multiple disease categories. The strongest association was reported between rs4285314 in the precursor of miR-3135b and celiac disease risk (odds ratio (OR) per effect allele increase = 0.37, P = 1.8 × 10-162). Colocalisation and Mendelian randomisation analysis highlighted potential causal role of miR-6891-3p in dyslipidaemia.

CONCLUSIONS

Our study demonstrates the pleiotropic effect of miRNAs and offers insights to their possible clinical importance.

Diagnostic utility of quantitative analysis of microRNA in bile samples obtained during endoscopic retrograde cholangiopancreatography for malignant biliary strictures

BACKGROUND

The sensitivity of bile cytology for malignant biliary strictures is not adequate. To overcome this limitation, we evaluated whether quantitative analysis of microRNAs (miRNAs) in bile can provide a precise diagnosis of malignant biliary strictures due to pancreatic cancer (PC) and biliary tract cancer (BTC).

METHODS

This was a retrospective evaluation of miRNA levels in stored bile samples of patients with PC, BTC or benign biliary stricture obtained during biliary drainage from April 2019 to December 2021 at our institution. A total of 113 patients (PC; n = 40, BTC; n = 38, control; n = 35) were enrolled. The miRNA candidates to be quantified were determined with microarray analysis from each 3 patients with PC, BTC and controls.

RESULTS

Using microarray analysis, we confirmed four significantly up-regulated miRNAs (miR-1275, miR-6891-5p, miR-7107-5p, miR-3197) in patients with PC and BTC compared to control patients. Quantitative PCR was then performed in 113 bile samples for these miRNAs. miR-1275 was significantly upregulated in PC (p = 0.003) and BTC (p = 0.049) compared to controls, miR-6891-5p was significantly upregulated in PC compared to controls (p = 0.025). In particular, a combination of bile cytology and miR-1275 in bile showed a sensitivity of 77.5% (95% CI, 70.7-77.5%), specificity of 100% (95% CI, 92.2-100%) and an area under the curve (AUC) of 0.93, and provided a significantly greater additional diagnostic effect than bile cytology alone (p = 0.014).

CONCLUSIONS

This study suggest that bile miRNAs could be potential biomarkers for pancreato-biliary diseases, particularly miR-1275 and miR-6891-5p may be helpful in the diagnosis of PC and BTC.

HLA Variation and SARS-CoV-2 Specific Antibody Response

Differences in SARS-CoV-2-specific immune responses have been observed between individuals following natural infection or vaccination. In addition to already known factors, such as age, sex, COVID-19 severity, comorbidity, vaccination status, hybrid immunity, and duration of infection, inter-individual variations in SARS-CoV-2 immune responses may, in part, be explained by structural differences brought about by genetic variation in the human leukocyte antigen (HLA) molecules responsible for the presentation of SARS-CoV-2 antigens to T effector cells. While dendritic cells present peptides with HLA class I molecules to CD8+ T cells to induce cytotoxic T lymphocyte responses (CTLs), they present peptides with HLA class II molecules to T follicular helper cells to induce B cell differentiation followed by memory B cell and plasma cell maturation. Plasma cells then produce SARS-CoV-2-specific antibodies. Here, we review published data linking HLA genetic variation or polymorphisms with differences in SARS-CoV-2-specific antibody responses. While there is evidence that heterogeneity in antibody response might be related to HLA variation, there are conflicting findings due in part to differences in study designs. We provide insight into why more research is needed in this area. Elucidating the genetic basis of variability in the SARS-CoV-2 immune response will help to optimize diagnostic tools and lead to the development of new vaccines and therapeutics against SARS-CoV-2 and other infectious diseases.

Genome-wide association studies of COVID-19: Connecting the dots

Genome-wide association studies (GWASs) are a research approach used to identify genetic variants associated with common diseases, like COVID-19. The lead genetic variants (n = 41) reported by the eleven largest COVID-19 GWASs are mapped to 22 different chromosomal regions. The loci 3q21.31 (LZTFL1 and chemokine receptor genes) and 9q34.2 (ABO), associated with disease severity and susceptibility to infection, respectively, were the most replicated findings across studies. Genes involved with mucociliary clearance (CEP97, FOXP4), viral-entry (ACE2, SLC6A20) and mucosal immunity (MIR6891) are associated with the risk of SARS-CoV-2 infection while genes of antiviral immune response (IFNAR2, OAS1), leukocyte trafficking (CCR9, CXCR6) and lung injury (DPP9, NOTCH4) are associated with severe disease. The biological processes underlying the risk of infection occur prominently, but not exclusively, in the upper airways whereas the severe COVID-19-associated processes in alveolar-capillary interface. The COVID-19 GWASs has unraveled key genetic mechanisms of SARS-CoV-2 pathogenesis, although the genetic basis of other COVID-19 related phenotypes (long COVID and neurological impairment) remains to be elucidated.

X-linked genes exhibit miR6891-5p-regulated skewing in Sjögren's syndrome

Many autoimmune diseases exhibit a strikingly increased prevalence in females, with primary Sjögren's syndrome (pSS) being the most female-predominant example. However, the molecular basis underlying the female-bias in pSS remains elusive. To address this knowledge gap, we performed genome-wide, allele-specific profiling of minor salivary gland-derived mesenchymal stromal cells (MSCs) from pSS patients and control subjects, and detected major differences in the regulation of X-linked genes. In control female MSCs, X-linked genes were expressed from both paternal and maternal X chromosomes with a median paternal ratio of ~ 0.5. However, in pSS female MSCs, X-linked genes exhibited preferential expression from one of the two X chromosomes. Concomitantly, pSS MSCs showed decrease in XIST levels and reorganization of H3K27me3+ foci in the nucleus. Moreover, the HLA-locus-expressed miRNA miR6891-5p was decreased in pSS MSCs. miR6891-5p inhibition in control MSCs caused XIST dysregulation, ectopic silencing, and allelic skewing. Allelic skewing was accompanied by the mislocation of protein products encoded by the skewed genes, which was recapitulated by XIST and miR6891-5p disruption in control MSCs. Our data reveal X skewing as a molecular hallmark of pSS and highlight the importance of restoring X-chromosomal allelic balance for pSS treatment. KEY MESSAGES: X-linked genes exhibit skewing in primary Sjögren's syndrome (pSS). X skewing in pSS associates with alterations in H3K27me3 deposition. pSS MSCs show decreased levels of miR6891-5p, a HLA-expressed miRNA. miR6891-5p inhibition causes H3K27me3 dysregulation and allelic skewing.

Whole Exome Sequencing Reveals Genetic Variants in HLA Class II Genes Associated With Transplant-free Survival of Indeterminate Acute Liver Failure

INTRODUCTION

Indeterminate acute liver failure (IND-ALF) is a rare clinical syndrome with a high mortality rate. Lacking a known etiology makes rapid evaluation and treatment difficult, with liver transplantation often considered as the only therapeutic option. Our aim was to identify genetic variants from whole exome sequencing data that might be associated with IND-ALF clinical outcomes.

METHODS

Bioinformatics analysis was performed on whole exome sequencing data for 22 patients with IND-ALF. A 2-tier approach was used to identify significant single-nucleotide polymorphisms (SNPs) associated with IND-ALF clinical outcomes. Tier 1 identified the SNPs with a higher relative risk in the IND-ALF population compared with those identified in control populations. Tier 2 determined the SNPs connected to transplant-free survival and associated with model for end-stage liver disease serum sodium and Acute Liver Failure Study Group prognostic scores.

RESULTS

Thirty-one SNPs were found associated with a higher relative risk in the IND-ALF population compared with those in controls, of which 11 belong to the human leukocyte antigen (HLA) class II genes but none for the class I. Further analysis showed that 5 SNPs: rs796202376, rs139189937, and rs113473719 of HLA-DRB5; rs9272712 of HLA-DQA1; and rs747397929 of IDO1 were associated with a higher probability of IND-ALF transplant-free survival. Using 3 selected SNPs, a model for the polygenic risk score was developed to predict IND-ALF prognoses, which are comparable with those by model for end-stage liver disease serum sodium and Acute Liver Failure Study Group prognostic scores.

DISCUSSION

Certain gene variants in HLA-DRB5, HLA-DQA1, and IDO1 were found associated with IND-ALF transplant-free survival. Once validated, these identified SNPs may help elucidate the mechanism of IND-ALF and assist in its diagnosis and management.

Preliminary Investigation of the Biomarkers of Acute Renal Transplant Rejection Using Integrated Proteomics Studies, Gene Expression Omnibus Datasets, and RNA Sequencing

A kidney transplant is often the best treatment for end-stage renal disease. Although immunosuppressive therapy sharply reduces the occurrence of acute allograft rejection (AR), it remains the main cause of allograft dysfunction. We aimed to identify effective biomarkers for AR instead of invasive kidney transplant biopsy. We integrated the results of several proteomics studies related to AR and utilized public data sources. Gene ontology (GO) and pathway analyses were used to identify important biological processes and pathways. The performance of the identified proteins was validated using several public gene expression omnibus (GEO) datasets. Samples that performed well were selected for further validation through RNA sequencing of peripheral blood mononuclear cells of patients with AR (n = 16) and non-rejection (n = 19) from our medical center. A total of 25 differentially expressed proteins (DEPs) overlapped in proteomic studies of urine and blood samples. GO analysis showed that the DEPs were mainly involved in the immune system and blood coagulation. Pathway analysis showed that the complement and coagulation cascade pathways were well enriched. We found that immunoglobulin heavy constant alpha 1 (IGHA1) and immunoglobulin κ constant (IGKC) showed good performance in distinguishing AR from non-rejection groups validated with several GEO datasets. Through RNA sequencing, the combination of IGHA1, IGKC, glomerular filtration rate, and donor age showed good performance in the diagnosis of AR with ROC AUC 91.4% (95% CI: 82–100%). Our findings may contribute to the discovery of potential biomarkers for AR monitoring.

Serum Exosomes MicroRNAs Are Novel Non-Invasive Biomarkers of Intrahepatic Cholestasis of Pregnancy

Background

Intrahepatic cholestasis of pregnancy (ICP) is closely related to the occurrence of adverse outcomes. Currently, total bile acids (TBAs) are the only diagnostic index for ICP, and its sensitivity and specificity have certain limitations. In this study, we aimed to develop potential biomarkers for the diagnosis of ICP.

Methods

Sixty pregnant women diagnosed with ICP and 48 healthy pregnant controls were enrolled in this study. We used the Agilent microRNA (miRNA) array followed by quantitative reverse transcriptase polymerase chain reaction assays to identify and validate the serum exosome miRNA profiles in ICP and healthy pregnant controls. We employed bioinformatics to identify metabolic processes associated with differentially expressed serum exosome miRNAs.

Results

The expression levels of hsa-miR-4271, hsa-miR-1275, and hsa-miR-6891-5p in maternal serum exosomes were significantly lower in ICP patients compared to controls; the diagnostic accuracy of hsa-miR-4271, hsa-miR-1275, and hsa-miR-6891-5p was evaluated with the area under the receiver operating characteristic curve (AUC) values of 0.861, 0.886, and 0.838, respectively. Multiple logistic regression analysis showed that a combination of the levels of hsa-miR-4271and hsa-miR-1275 afforded a significantly higher AUC (0.982). The non-error rate of a combination of all three exosome miRNAs was the highest (95%), thus more reliable ICP diagnosis. The expression levels of all three exosome miRNAs were negatively associated with TBAs. Furthermore, according to bioinformatics analysis, the three exosome miRNAs were related to lipid metabolism, apoptosis, oxidative stress, and the Mitogen Activated Protein Kinase (MAPK) signaling pathway.

Conclusions

This study may identify the novel non-invasive biomarkers for ICP and provided new insights into the important role of the exosome miRNA regulation in ICP.

STXBP6 and B3GNT6 Genes are Associated With Selective IgA Deficiency

Immunoglobulin A Deficiency (IgAD) is a polygenic primary immune deficiency, with a strong genetic association to the human leukocyte antigen (HLA) region. Previous genome-wide association studies (GWAS) have identified five non-HLA risk loci (IFIH1, PVT1, ATG13-AMBRA1, AHI1 and CLEC16A). In this study, we investigated the genetic interactions between different HLA susceptibility haplotypes and non-MHC genes in IgAD. To do this, we stratified IgAD subjects and healthy controls based on HLA haplotypes (N = 10,993), and then performed GWAS to identify novel genetic regions contributing to IgAD susceptibility. After replicating previously published HLA risk haplotypes, we compared individuals carrying at least one HLA risk allele (HLA-B*08:01-DRB1*03:01-DQB1*02:01 or HLA-DRB1*07:01-DQB1*02:02 or HLA-DRB1*01-DQB1*05:01) with individuals lacking an HLA risk allele. Subsequently, we stratified subjects based on the susceptibility alleles/haplotypes and performed gene-based association analysis using 572,856 SNPs and 24,125 genes. A significant genome-wide association in STXBP6 (rs4097492; p = 7.63 × 10^-9) was observed in the cohort carrying at least one MHC risk allele. We also identified a significant gene-based association for B3GNT6 (P _Gene = 2.1 × 10^-6) in patients not carrying known HLA susceptibility alleles. Our findings indicate that the etiology of IgAD differs depending on the genetic background of HLA susceptibility haplotypes.

The Role of HLA in the Association between IgA Deficiency and Celiac Disease

Selective IgA deficiency (SIgAD) is the most frequent primary immune defect. Since SIgAD is not characterized by relevant infectious issues in most cases, it is often diagnosed during the diagnostic work up of several and different autoimmune disorders, which are associated with this primary immune defect. The genetic background of SIgAD is complex and three HLA haplotypes resulted to be more frequently associated with it; in detail, two of them include HLA-DQB1^∗02 allelic variants, which are essential predisposing factors to develop Celiac Disease (CD). Here, we discuss the evidence regarding the role of HLA in the etiopathogenesis of SIgAD and its association with CD. Actually, the HLA region seems to play a modest role in the genetic predisposition to SIgAD and we may speculate that the association with the HLA-DQB1^∗02 alleles (or haplotypes including them) could derive from its link with CD. Indeed, SIgAD and some related immunological alterations are likely to predispose to several autoimmune diseases (with and despite different HLA backgrounds), including CD, which is relatively common and directly associated with the HLA-DQB1^∗02 allelic variants coding the DQ2 heterodimer. Further and specific studies are needed to make final conclusions in this regard.

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.

Viruses, cancer and non-self recognition

Virus-host interactions form an essential part of every aspect of life, and this review is aimed at looking at the balance between the host and persistent viruses with a focus on the immune system. The virus-host interaction is like a cat-and-mouse game and viruses have developed ingenious mechanisms to manipulate cellular pathways, most notably the major histocompatibility (MHC) class I pathway, to reside within infected cell while evading detection and destruction by the immune system. However, some of the signals sensing and responding to viral infection are derived from viruses and the fact that certain viruses can prevent the infection of others, highlights a more complex coexistence between the host and the viral microbiota. Viral immune evasion strategies also illustrate that processes whereby cells detect and present non-self genetic material to the immune system are interlinked with other cellular pathways. Immune evasion is a target also for cancer cells and a more detailed look at the interfaces between viral factors and components of the MHC class I peptide-loading complex indicates that these interfaces are also targets for cancer mutations. In terms of the immune checkpoint, however, viral and cancer strategies appear different.

Regulatory noncoding RNAs and the major histocompatibility complex

The Major Histocompatibility Complex (MHC) is a 4 Mbp genomic region located on the short arm of chromosome 6. The MHC region contains many key immune-related genes such as Human Leukocyte Antigens (HLAs). There has been a growing realization that, apart from MHC encoded proteins, RNAs derived from noncoding regions of the MHC-specifically microRNAs (miRNAs) and long noncoding RNAs (lncRNAs)-play a significant role in cellular regulation. Furthermore, regulatory noncoding RNAs (ncRNAs) derived from other parts of the genome fine-tune the expression of many immune-related MHC proteins. Although the field of ncRNAs of the MHC is a research area that is still in its infancy, ncRNA regulation of MHC genes has already been shown to be vital for immune function, healthy pregnancy and cellular homeostasis. Dysregulation of this intricate network of ncRNAs can lead to serious perturbations in homeostasis and subsequent disease.

Capturing Differential Allele-Level Expression and Genotypes of All Classical HLA Loci and Haplotypes by a New Capture RNA-Seq Method

The highly polymorphic human major histocompatibility complex (MHC) also known as the human leukocyte antigen (HLA) encodes class I and II genes that are the cornerstone of the adaptive immune system. Their unique diversity (>25,000 alleles) might affect the outcome of any transplant, infection, and susceptibility to autoimmune diseases. The recent rapid development of new next-generation sequencing (NGS) methods provides the opportunity to study the influence/correlation of this high level of HLA diversity on allele expression levels in health and disease. Here, we describe the NGS capture RNA-Seq method that we developed for genotyping all 12 classical HLA loci (HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, HLA-DRA, HLA-DRB1, HLA-DRB3, HLA-DRB4, and HLA-DRB5) and assessing their allelic imbalance by quantifying their allele RNA levels. This is a target enrichment method where total RNA is converted to a sequencing-ready complementary DNA (cDNA) library and hybridized to a complex pool of RNA-specific HLA biotinylated oligonucleotide capture probes, prior to NGS. This method was applied to 161 peripheral blood mononuclear cells and 48 umbilical cord blood cells of healthy donors. The differential allelic expression of 10 HLA loci (except for HLA-DRA and HLA-DPA1) showed strong significant differences (P < 2.1 × 10-15). The results were corroborated by independent methods. This newly developed NGS method could be applied to a wide range of biological and medical questions including graft rejections and HLA-related diseases.

Role of microRNA 690 in Mediating Angiotensin II Effects on Inflammation and Endoplasmic Reticulum Stress

Overactivation of the renin–angiotensin system (RAS) during obesity disrupts adipocyte metabolic homeostasis and induces endoplasmic reticulum (ER) stress and inflammation; however, underlying mechanisms are not well known. We propose that overexpression of angiotensinogen (Agt), the precursor protein of RAS in adipose tissue or treatment of adipocytes with Angiotensin II (Ang II), RAS bioactive hormone, alters specific microRNAs (miRNA), that target ER stress and inflammation leading to adipocyte dysfunction. Epididymal white adipose tissue (WAT) from B6 wild type (Wt) and transgenic male mice overexpressing Agt (Agt-Tg) in adipose tissue and adipocytes treated with Ang II were used. Small RNA sequencing and microarray in WAT identified differentially expressed miRNAs and genes, out of which miR-690 and mitogen-activated protein kinase kinase 3 (MAP2K3) were validated as significantly up- and down-regulated, respectively, in Agt-Tg, and in Ang II-treated adipocytes compared to respective controls. Additionally, the direct regulatory role of miR-690 on MAP2K3 was confirmed using mimic, inhibitors and dual-luciferase reporter assay. Downstream protein targets of MAP2K3 which include p38, NF-κB, IL-6 and CHOP were all reduced. These results indicate a critical post-transcriptional role for miR-690 in inflammation and ER stress. In conclusion, miR-690 plays a protective function and could be a useful target to reduce obesity.

MHC Class I Regulation: The Origin Perspective

Viral-derived elements and non-coding RNAs that build up “junk DNA” allow for flexible and context-dependent gene expression. They are extremely dense in the MHC region, accounting for flexible expression of the MHC I, II, and III genes and adjusting the level of immune response to the environmental stimuli. This review brings forward the viral-mediated aspects of the origin and evolution of adaptive immunity and aims to link this perspective with the MHC class I regulation. The complex regulatory network behind MHC expression is largely controlled by virus-derived elements, both as binding sites for immune transcription factors and as sources of regulatory non-coding RNAs. These regulatory RNAs are imbalanced in cancer and associate with different tumor types, making them promising targets for diagnostic and therapeutic interventions.

Increased Ileal Immunoglobulin A Production and Immunoglobulin A-Coated Bacteria in Diarrhea-Predominant Irritable Bowel Syndrome

OBJECTIVES

Immune activation and intestinal microbial dysbiosis could induce diarrhea-predominant irritable bowel syndrome (IBS-D). We examined the roles of ileal immunoglobulin A (IgA) and IgA-coated bacteria in IBS-D pathogenesis.

METHODS

Peripheral blood, fecal samples, and ileal and cecal biopsies were collected from 32 healthy volunteers and 44 patients with IBS-D. Quantitative reverse transcriptase polymerase chain reaction was used to assess differential gene expression. IgA levels in the blood and fecal samples were quantified by an enzyme-linked immunosorbent assay. IgA cells were assessed by immunofluorescence imaging. Flow-cytometry-based IgA bacterial cell sorting and 16S rRNA gene sequencing (IgA-SEQ) was used to isolate and identify fecal IgA bacteria.

RESULTS

Fecal IgA, particularly IgA1, was upregulated in patients with IBS-D. IgA class switch and B cell-activating factor-receptor were increased in the terminal ileum of patients. The intestinal microbiota composition was altered in patients compared with that in controls. IgA-SEQ showed that the proportion of fecal IgA-coated bacteria was increased significantly in patients with IBS-D. IgA bacteria in patients with IBS-D showed higher abundances of Escherichia-Shigella, Granulicatella, and Haemophilus compared with healthy controls and IgA bacteria in patients with IBS-D. The Escherichia-Shigella IgA coating index was positively correlated with anxiety and depression. The Escherichia-Shigella relative abundance, luminal IgA activity, and some altered IgA-coated bacteria were positively associated with the clinical manifestations of IBS-D.

DISCUSSION

Microbial dysbiosis may promote the terminal ileal mucosa to produce higher levels of IgA, increasing the proportion of IgA-coated bacteria by activating IgA class switching, which might regulate local inflammation and clinical manifestations in IBS-D. IgA may mediate the effects of microbial dysbiosis on the pathogenesis of IBS-D.

Exploring the RNA Gap for Improving Diagnostic Yield in Primary Immunodeficiencies

Challenges in diagnosing primary immunodeficiency are numerous and diverse, with current whole-exome and whole-genome sequencing approaches only able to reach a molecular diagnosis in 25–60% of cases. We assess these problems and discuss how RNA-focused analysis has expanded and improved in recent years and may now be utilized to gain an unparalleled insight into cellular immunology. We review how investigation into RNA biology can give information regarding the differential expression, monoallelic expression, and alternative splicing—which have important roles in immune regulation and function. We show how this information can inform bioinformatic analysis pipelines and aid in the variant filtering process, expediting the identification of causal variants—especially those affecting splicing—and enhance overall diagnostic ability. We also demonstrate the challenges, which remain in the design of this type of investigation, regarding technological limitation and biological considerations and suggest potential directions for the clinical applications.

MHC class I polymorphic Alu insertion (POALIN) allele and haplotype frequencies in the Arabs of the United Arab Emirates and other world populations

Polymorphic Alu insertions (POALINs) are found throughout the human genome and have been used in various studies to infer geographic origin of human populations. The main aim of this study was to determine the allele and haplotype frequencies of five POALINs, AluHF, AluHG, AluHJ, AluTF and AluMICB, within the major histocompatibility complex (MHC) class I region of 95 UAE Arabs, and correlate their frequencies to those of the HLA-A, HLA-C and HLA-B class I allele lineages. Evolutionary relationships between the POALINs of the Arabs and those previously studied in populations of African, Asian and European descent were compared. At each of the five Alu loci (AluHF, AluHG, AluHJ, AluTF and AluMICB), Alu insertion was designated as Alu(locus)*02 and absence was Alu(locus)*01. The AluHG insertion (AluHG*02) had the highest frequency (0.332), followed by AluHF*02 (0.300), AluHJ*02 (0.263), AluMICB*02 (0.111) and AluTF*02 (0.058). Of the 270 Alu-HLA haplotypes pairs in the UAE Arabs, 110 had no Alu insertion, and 54 had an Alu insertion at >50% per haplotype. An Alu insertion >75% per haplotype was found between AluMICB*02 and HLA-B*14, HLA-B*22, HLA-B*44, HLA-B*55, HLA-B*57 and HLA-B*73, and with HLA-C*01 and HLA-C*18; AluHJ*02 with HLA-A*01, HLA-A*19, HLA-A*24 and HLA-A*32; AluHG*02 with HLA-A*02 and HLA-B*18; and AluHF*02 with HLA-A*10. The genotyped allele and haplotype frequencies of the MHC POALINs in UAE Arabs were compared with the results of 30 previously published Asian, European, American and African populations. Phylogenetic and multidimensional scaling (MDS) analysis of the relative MHC POALINs allele and haplotype frequencies revealed that the UAE Arabs have a similar lineage to Caucasians and the most distant genetic relationship to the Waorani native American population of Ecuador. The structure of both the phylogenetic tree and the MDS analysis supports the Out of Africa theory of human evolution. The nature of the clusters suggests the Arabian Middle East represents a crossroads from which human populations migrated towards Asia in the east and Europe to the north-west.

Next-Generation Sequencing Based HLA Typing: Deciphering Immunogenetic Aspects of Sarcoidosis

Unraveling of the HLA-related immunogenetic basis of several immune disorders is complex due to the extensive HLA polymorphism and strong linkage-disequilibrium between HLA loci. A lack of in phase sequence information, a relative deficiency of high resolution genotyping including non-coding regions and ambiguous haplotype assignment make it difficult to compare findings across association studies and to attribute a causal role to specific HLA alleles/haplotypes in disease susceptibility and modification of disease phenotypes. Earlier, historical antibody and DNA-based methods of HLA typing, primarily of low resolution at antigen/alellic group levels, yielded "indicative" findings which were partially improved by high-resolution DNA-based typing. Only recently, next-generation sequencing (NGS) approaches based on deep-sequencing of the complete HLA genes combined with bioinformatics tools began to provide the access to complete information at an allelic level. Analyzing HLA with NGS approaches, therefore, promises to provide further insight in the etiopathogenesis of several immune disorders in which HLA associations have been implicated. These range from coeliac disease and rheumatological conditions to even more complex disorders, such as type-1 diabetes, systemic lupus erythematosus and sarcoidosis. A systemic disease of unknown etiology, sarcoidosis has previously been associated with numerous HLA variants and also other gene polymorphisms, often in linkage with the HLA region. To date, the biological significance of these associations has only partially been defined. Therefore, more precise assignments of HLA alleles/haplotypes using NGS approaches could help to elucidate the exact role of HLA variation in the multifaceted etiopathogenesis of sarcoidosis, including epigenetic mechanisms. NGS-based HLA analyses may be also relevant for defining variable clinical phenotypes and for predicting the disease course or the response to current/plausible novel therapies.

Reference Grade Characterization of Polymorphisms in Full-Length HLA Class I and II Genes With Short-Read Sequencing on the ION PGM System and Long-Reads Generated by Single Molecule, Real-Time Sequencing on the PacBio Platform

Although NGS technologies fuel advances in high-throughput HLA genotyping methods for identification and classification of HLA genes to assist with precision medicine efforts in disease and transplantation, the efficiency of these methods are impeded by the absence of adequately-characterized high-frequency HLA allele reference sequence databases for the highly polymorphic HLA gene system. Here, we report on producing a comprehensive collection of full-length HLA allele sequences for eight classical HLA loci found in the Japanese population. We augmented the second-generation short read data generated by the Ion Torrent technology with long amplicon spanning consensus reads delivered by the third-generation SMRT sequencing method to create reference grade high-quality sequences of HLA class I and II gene alleles resolved at the genomic coding and non-coding level. Forty-six DNAs were obtained from a reference set used previously to establish the HLA allele frequency data in Japanese subjects. The samples included alleles with a collective allele frequency in the Japanese population of more than 99.2%. The HLA loci were independently amplified by long-range PCR using previously designed HLA-locus specific primers and subsequently sequenced using SMRT and Ion PGM sequencers. The mapped long and short-reads were used to produce a reference library of consensus HLA allelic sequences with the help of the reference-aware software tool LAA for SMRT Sequencing. A total of 253 distinct alleles were determined for 46 healthy subjects. Of them, 137 were novel alleles: 101 SNVs and/or indels and 36 extended alleles at a partial or full-length level. Comparing the HLA sequences from the perspective of nucleotide diversity revealed that HLA-DRB1 was the most divergent among the eight HLA genes, and that the HLA-DPB1 gene sequences diverged into two distinct groups, DP2 and DP5, with evidence of independent polymorphisms generated in exon 2. We also identified two specific intronic variations in HLA-DRB1 that might be involved in rheumatoid arthritis. In conclusion, full-length HLA allele sequencing by third-generation and second-generation technologies has provided polymorphic gene reference sequences at a genomic allelic resolution including allelic variations assigned up to the field-4 level for a stronger foundation in precision medicine and HLA-related disease and transplantation studies.

Epigenetic Regulation of Autophagy: A Path to the Control of Autoimmunity

Autoimmune diseases are a significant cause of debilitation and mortality globally and are in need of cost-effective therapeutics. Autophagy is a cellular pathway that facilitates immune modulation involved in both pathogen control and autoimmunity. Regulation is multifactorial and includes a number of epigenetic pathways which can involve modification of DNA-binding histones to induce autophagy-related mRNA synthesis or microRNA and decapping-associated mRNA degradation which results in autophagy suppression. Appreciation of epigenetic-based pathways involved in autophagy and autoimmunity may facilitate application of a burgeoning group of epigenetic pharmaceuticals to these important diseases.

Full-length extension of HLA allele sequences by HLA allele-specific hemizygous Sanger sequencing (SSBT)

The gold standard for typing at the allele level of the highly polymorphic Human Leucocyte Antigen (HLA) gene system is sequence based typing. Since sequencing strategies have mainly focused on identification of the peptide binding groove, full-length sequence information is lacking for >90% of the HLA alleles. One of the goals of the 17th IHIWS workshop is to establish full-length sequences for as many HLA alleles as possible. In our component "Extension of HLA sequences by full-length HLA allele-specific hemizygous Sanger sequencing" we have used full-length hemizygous Sanger Sequence Based Typing to achieve this goal. We selected samples of which full length sequences were not available in the IPD-IMGT/HLA database. In total we have generated the full-length sequences of 48 HLA-A, 45 -B and 31 -C alleles. For HLA-A extended alleles, 39/48 showed no intron differences compared to the first allele of the corresponding allele group, for HLA-B this was 26/45 and for HLA-C 20/31. Comparing the intron sequences to other alleles of the same allele group revealed that in 5/48 HLA-A, 16/45 HLA-B and 8/31 HLA-C alleles the intron sequence was identical to another allele of the same allele group. In the remaining 10 cases, the sequence either showed polymorphism at a conserved nucleotide or was the result of a gene conversion event. Elucidation of the full-length sequence gives insight in the polymorphic content of the alleles and facilitates the identification of its evolutionary origin.

Novel and Haplotype Specific MicroRNAs Encoded by the Major Histocompatibility Complex

The MHC is recognized for its importance in human health and disease. However, many disease-associated variants throughout the region remain of unknown significance, residing predominantly within non-coding regions of the MHC. The characterization of non-coding RNA transcripts throughout the MHC is thus central to understanding the genetic contribution of these variants. Therefore, we characterize novel miRNA transcripts throughout the MHC by performing deep RNA sequencing of two B lymphoblastoid cell lines with completely characterized MHC haplotypes. Our analysis identifies 89 novel miRNA transcripts, 48 of which undergo Dicer-dependent biogenesis and are loaded onto the Argonaute silencing complex. Several of the identified mature miRNA and pre-miRNA transcripts are unique to specific MHC haplotypes and overlap common SNPs. Furthermore, 43 of the 89 identified novel miRNA transcripts lie within linkage disequilibrium blocks that contain a disease-associated SNP. These disease associated SNPs are associated with 65 unique disease phenotypes, suggesting that these transcripts may play a role in the etiology of numerous diseases associated with the MHC. Additional in silico analysis reveals the potential for thousands of putative pre-miRNA encoding loci within the MHC that may be expressed by different cell types and at different developmental stages.

Peroxiredoxin-2 and zinc-alpha-2-glycoprotein as potentially combined novel salivary biomarkers for early detection of oral squamous cell carcinoma using proteomic approaches

No effective screening method is available for oral squamous cell carcinoma (OSCC) that is recognized to influence by environmental factors as well as human papillomavirus (HPV) and Epstein-Barr virus (EBV). Therefore, we sought to identify salivary biomarkers for screening of OSCC with or without HPV and/or EBV infection. Saliva, lesion and oral exfoliated cells were collected from OSCC patients and cancer-free controls (CFCs) and grouped depending on their HPV- and EBV-infection status. Salivary protein was precipitated and subjected to 2-dimensional gel electrophoresis. Differential expression of proteins was identified by mass spectrometry and validated by Western blotting. Distinctive expression patterns of salivary proteins were detected in OSCC as compared with CFCs. Levels of peroxiredoxin-2 (PRDX-2) and zinc-alpha-2-glycoprotein (ZAG) were significantly up-regulated in OSCC cases (p<0.001) relative to CFCs. Similarly, these proteins were also up-regulated in lesion cells compared with oral exfoliated cells (p<0.001). However, the expression patterns of these proteins were not significantly influenced by patient histories (risk factors). In combination, these proteins yielded the highest discriminatory power (AUC=0.999), sensitivity (100%), and specificity (98.77%) in distinguishing the early stages of OSCC. The detection of PRDX-2 combining with ZAG protein could potentially be used as salivary biomarkers for early screening of OSCC.

SIGNIFICANCE

Our findings demonstrate a useful of combined detection of PRDX-2 and ZAG as a salivary biomarker for the early detection of OSCC.