Computational assessment of miRNA binding to low and high expression HLA-DPB1 allelic sequences

HLA-DPB1 genotype variants predict DP molecule cell surface expression and DP donor specific antibody binding capacity

The contribution of alloresponses to mismatched HLA-DP in solid organ transplantation and hematopoietic stem cell transplantation (HCT) has been well documented. Exploring the regulatory mechanisms of DPB1 alleles has become an important question to be answered. In this study, our initial investigation focused on examining the correlation between the rs9277534G/A SNP and DPB1 mRNA expression. The result showed that there was a significant increase in DPB1 mRNA expression in B lymphoblastoid cell lines (BLCLs) with the rs9277534GG genotype compared to rs9277534AA genotype. In addition, B cells with the rs9277534GG exhibited significantly higher DP protein expression than those carrying the rs9277534AA genotype in primary B cells. Furthermore, we observed a significant upregulation of DP expression in B cells following treatment with Interleukin 13 (IL-13) compared to untreated B cells carrying rs9277534GG-linked DPB1 alleles. Fluorescence in situ hybridization (FISH) analysis of DPB1 in BLCL demonstrated significant differences in both the cytoplasmic (p=0.0003) and nuclear (p=0.0001) localization of DP mRNA expression comparing DPB1*04:01 (rs9277534AA) and DPB1*05:01 (rs9277534GG) homozygous cells. The study of the correlation between differential DPB1 expression and long non-coding RNAs (lncRNAs) showed that lnc-HLA-DPB1-13:1 is strongly associated with DP expression (r=0.85), suggesting the potential involvement of lncRNA in regulating DP expression. The correlation of DP donor specific antibody (DSA) with B cell flow crossmatch (B-FCXM) results showed a better linear correlation of DP DSA against GG and AG donor cells (R2 = 0.4243, p=0.0025 and R2 = 0.6172, p=0.0003, respectively), compared to DSA against AA donor cells (R2 = 0.0649, p=0.4244). This explained why strong DP DSA with a low expression DP leads to negative B-FCXM. In conclusion, this study provides evidence supporting the involvement of lncRNA in modulating HLA-DP expression, shedding lights on the intricate regulatory mechanisms of DP, particularly under inflammatory conditions in transplantation.

In silico prioritisation of microRNA-associated common variants in multiple sclerosis

BACKGROUND

Genome-wide association studies (GWAS) have highlighted over 200 autosomal variants associated with multiple sclerosis (MS). However, variants in non-coding regions such as those encoding microRNAs have not been explored thoroughly, despite strong evidence of microRNA dysregulation in MS patients and model organisms. This study explores the effect of microRNA-associated variants in MS, through the largest publicly available GWAS, which involved 47,429 MS cases and 68,374 controls.

METHODS

We identified SNPs within the coordinates of microRNAs, ± 5-kb microRNA flanking regions and predicted 3'UTR target-binding sites using miRBase v22, TargetScan 7.0 RNA22 v2.0 and dbSNP v151. We established the subset of microRNA-associated SNPs which were tested in the summary statistics of the largest MS GWAS by intersecting these datasets. Next, we prioritised those microRNA-associated SNPs which are among known MS susceptibility SNPs, are in strong linkage disequilibrium with the former or meet a microRNA-specific Bonferroni-corrected threshold. Finally, we predicted the effects of those prioritised SNPs on their microRNAs and 3'UTR target-binding sites using TargetScan v7.0, miRVaS and ADmiRE.

RESULTS

We have identified 30 candidate microRNA-associated variants which meet at least one of our prioritisation criteria. Among these, we highlighted one microRNA variant rs1414273 (MIR548AC) and four 3'UTR microRNA-binding site variants within SLC2A4RG (rs6742), CD27 (rs1059501), MMEL1 (rs881640) and BCL2L13 (rs2587100). We determined changes to the predicted microRNA stability and binding site recognition of these microRNA and target sites.

CONCLUSIONS

We have systematically examined the functional, structural and regulatory effects of candidate MS variants among microRNAs and 3'UTR targets. This analysis allowed us to identify candidate microRNA-associated MS SNPs and highlights the value of prioritising non-coding RNA variation in GWAS. These candidate SNPs could influence microRNA regulation in MS patients. Our study is the first thorough investigation of both microRNA and 3'UTR target-binding site variation in multiple sclerosis using GWAS summary statistics.

Human Leukocyte Antigen Association Study Reveals DRB1*04:02 Effects Additional to DRB1*07:01 in Anti-LGI1 Encephalitis

BACKGROUND AND OBJECTIVES

To study human leukocyte antigen (HLA) allele associations in anti-leucine-rich glioma-inactivated 1 (LGI1) encephalitis.

METHODS

A multiethnic cohort of 269 patients with anti-LGI1 encephalitis and 1,359 controls was included. Four-digit HLA sequencing and genome wide association single-nucleotide polymorphism typing imputation (0.99 concordance) were used for HLA typing. Significance of primary and secondary associations was tested using χ2, Fisher exact tests, or logistic regression with the control of population stratification covariates when applicable.

RESULTS

DRB1*07:01 and DQA1*02:01, 2 alleles in strong linkage disequilibrium, were associated with the disease (90% vs 24%, OR = 27.8, p < 10e-50) across ethnicity independent of variation at DRB3 and DQB1, 2 flanking HLA loci. DRB1*07:01 homozygosity was associated with a doubling of risk (OR = 2.1, p = 0.010), suggesting causality. DRB1*07:01 negative subjects were younger (p = 0.003) and more frequently female (p = 0.015). Three patients with malignant thymomas did not carry DRB1*07:01, whereas patients with other tumors had high DRB1*07:01 frequency, suggesting that the presence of tumors other than thymomas may be coincidental and not causal. In both DRB1*07:01 heterozygous individuals and DRB1*07:01 negative subjects, DRB1*04:02 was associated with anti-LGI1 encephalitis, indicating an independent effect of this allele (OR = 6.85, p = 4.57 × 10-6 and OR = 8.93, p = 2.50 × 10-3, respectively). DRB1*04:02 was also independently associated with younger age at onset (β = -6.68, p = 9.78 × 10-3). Major histocompatibility complex peptide-binding predictions using LGI1-derived peptides revealed divergent binding propensities for DRB1*04:02 and DRB1*07:01 alleles, suggesting independent pathogenic mechanisms.

DISCUSSION

In addition to the established primary DRB1*07:01 association in anti-LGI1 encephalitis, we observe a secondary effect of DRB1*04:02 with lower age at onset. Our study provides evidence for secondary effects within HLA locus that correlate with clinical phenotypes in anti-LGI1 encephalitis.

Next-generation sequencing and clinical histocompatibility testing

Histocompatibility testing is essential for donor identification and risk assessment in solid organ and hematopoietic stem cell transplant. Additionally, it is useful for identifying donor specific alleles for monitoring donor specific antibodies in post-transplant patients. Next-generation sequence (NGS) based human leukocyte antigen (HLA) typing has improved many aspects of histocompatibility testing in hematopoietic stem cell and solid organ transplant. HLA disease association testing and research has also benefited from the advent of NGS technologies. In this review we discuss the current impact and future applications of NGS typing on clinical histocompatibility testing for transplant and non-transplant purposes.

Complex Linkage Disequilibrium Effects in HLA-DPB1 Expression and Molecular Mismatch Analyses of Transplantation Outcomes

BACKGROUND

HLA molecular mismatch (MM) is a risk factor for de novo donor-specific antibody (dnDSA) development in solid organ transplantation. HLA expression differences have also been associated with adverse outcomes in hematopoietic cell transplantation. We sought to study both MM and expression in assessing dnDSA risk.

METHODS

One hundred three HLA-DP-mismatched solid organ transplantation pairs were retrospectively analyzed. MM was computed using amino acids (aa), eplets, and, supplementarily, Grantham/Epstein scores. DPB1 alleles were classified as rs9277534-A (low-expression) or rs9277534-G (high-expression) linked. To determine the associations between risk factors and dnDSA, logistic regression, linkage disequilibrium (LD), and population-based analyses were performed.

RESULTS

A high-risk AA:GX (recipient:donor) expression combination (X = A or G) demonstrated strong association with HLA-DP dnDSA (P = 0.001). MM was also associated with HLA-DP dnDSA when evaluated by itself (eplet P = 0.007, aa P = 0.003, Grantham P = 0.005, Epstein P = 0.004). When attempting to determine the relative individual effects of the risk factors in multivariable analysis, only AA:GX expression status retained a strong association (relative risk = 18.6, P = 0.007 with eplet; relative risk = 15.8, P = 0.02 with aa), while MM was no longer significant (eplet P = 0.56, aa P = 0.51). Importantly, these risk factors are correlated, due to LD between the expression-tagging single-nucleotide polymorphism and polymorphisms along HLA-DPB1.

CONCLUSIONS

The MM and expression risk factors each appear to be strong predictors of HLA-DP dnDSA and to possess clinical utility; however, these two risk factors are closely correlated. These metrics may represent distinct ways of characterizing a common overlapping dnDSA risk profile, but they are not independent. Further, we demonstrate the importance and detailed implications of LD effects in dnDSA risk assessment and possibly transplantation overall.

APOC3 rs2070667 Associates with Serum Triglyceride Profile and Hepatic Inflammation in Nonalcoholic Fatty Liver Disease

Single-nucleotide polymorphisms (SNPs) of apolipoprotein C3 (APOC3) play important role in lipid metabolism, and dyslipidemia underlies nonalcoholic fatty liver disease (NAFLD). But the correlation of serum lipidomics, APOC3 SNPs, and NAFLD remains limited understood. Enrolling thirty-four biopsy-proven NAFLD patients from Tianjin, Shanghai, Fujian, we investigated their APOC3 genotype and serum lipid profile by DNA sequencing and ultraperformance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS), respectively. Scoring of hepatocyte steatosis, ballooning, lobular inflammation, and liver fibrosis was then performed to reveal the role of lipidomics-affecting APOC3 SNPs in NAFLD-specific pathological alterations. Here, we reported that APOC3 SNPs (rs4225, rs4520, rs5128, rs2070666, and rs2070667) intimately correlated to serum lipidomics in NAFLD patients. A allele instead of G allele at rs2070667, which dominated the SNPs underlying lipidomic alteration, exhibited downregulatory effect on triacylglycerols (TGs: TG 54 : 7, TG 54 : 8, and TG 56 : 9) containing polyunsaturated fatty acid (PUFA). Moreover, subjects with low-level PUFA-containing TGs were predisposed to high-grade lobular inflammation (TG 54 : 7, rho = -0.454 and P = 0.007; TG 54 : 8, rho = -0.411 and P =0.016; TG 56 : 9, rho = -0.481 and P = 0.004). The significant correlation of APOC3 rs2070667 and inflammation grading [G/G vs. G/A+A/A: 0.00 (0.00 and 1.00) vs. 1.50 (0.75 and 2.00), P = 0.022] further confirmed its pathological action on the basis of lipidomics-impacting activity. These findings suggest an inhibitory effect of A allele at APOC3 rs2070667 on serum levels of PUFA-containing TGs, which are associated with high-grade lobular inflammation in NAFLD patients.

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.

Factors affecting HLA expression: A review

The detection and semiquantitative measurement of circulating human leucocyte antigen (HLA)-specific antibodies is essential for the management of patients before and after transplantation. In addition, the pretransplant cross-match to assess the reactivity of recipient HLA antibody against donor lymphocytes has long been the gold standard to prevent hyperacute rejection. Whilst both of these tests assume that recipient HLA-specific antibody is the only variable in the assessment of transplant risk, this is not the case. Transplant immunologists recognize that some HLA antigens are expressed at levels a magnitude lower than others (e.g., HLA-C, HLA-DQ), but within loci, and between different cell types there are many factors that influence HLA expression in both resting and activated cells. HLA is not usually expressed without the specific promoter proteins NLRC5, for HLA class I, and CIITA, for class II. The quantity of HLA protein production is then affected by factors including promoter region polymorphisms, alternative exon splice sites, methylation and microRNA-directed degradation. Different loci are influenced by multiple combinations of these control mechanisms making prediction of HLA regulation difficult, but an ability to measure the cellular expression of each HLA antigen, in conjunction with knowledge of circulating HLA-specific antibody, would lead to a more informed algorithm to assess transplant risk.

Using Whole Genome Sequencing in an African Subphenotype of Myasthenia Gravis to Generate a Pathogenetic Hypothesis

Myasthenia gravis (MG) is a rare, treatable antibody-mediated disease which is characterized by muscle weakness. The pathogenic antibodies are most frequently directed at the acetylcholine receptors (AChRs) at the skeletal muscle endplate. An ophthalmoplegic subphenotype of MG (OP-MG), which is characterized by treatment resistant weakness of the extraocular muscles (EOMs), occurs in a proportion of myasthenics with juvenile symptom onset and African genetic ancestry. Since the pathogenetic mechanism(s) underlying OP-MG is unknown, the aim of this study was to use a hypothesis-generating genome-wide analysis to identify candidate OP-MG susceptibility genes and pathways. Whole genome sequencing (WGS) was performed on 25 AChR-antibody positive myasthenic individuals of African genetic ancestry sampled from the phenotypic extremes: 15 with OP-MG and 10 individuals with control MG (EOM treatment-responsive). Variants were called according to the Genome Analysis Toolkit (GATK) best practice guidelines using the hg38 reference genome. In addition to single variant association analysis, variants were mapped to genes (±200 kb) using VEGAS2 to calculate gene-based test statistics and HLA allele group assignment was inferred through "best-match" alignment of reads against the IMGT/HLA database. While there were no single variant associations that reached genome-wide significance in this exploratory sample, several genes with significant gene-based test statistics and known to be expressed in skeletal muscle had biological functions which converge on muscle atrophy signaling and myosin II function. The closely linked HLA-DPA1 and HLA-DPB1 genes were associated with OP-MG subjects (gene-based p < 0.05) and the frequency of a functional A > G SNP (rs9277534) in the HLA-DPB1 3'UTR, which increases HLA-DPB1 expression, differed between the two groups (G-allele 0.30 in OP-MG vs. 0.60 in control MG; p = 0.04). Furthermore, we show that rs9277534 is an HLA-DBP1 expression quantitative trait locus in patient-derived myocytes (p < 1 × 10^-3). The application of a SNP to gene to pathway approach to this exploratory WGS dataset of African myasthenic individuals, and comparing dichotomous subphenotypes, resulted in the identification of candidate genes and pathways that may contribute to OP-MG susceptibility. Overall, the hypotheses generated by this work remain to be verified by interrogating candidate gene and pathway expression in patient-derived extraocular muscle.