Haplotype-Specific Expression Analysis of MHC Class II Genes in Healthy Individuals and Rheumatoid Arthritis Patients

The effect of type 1 diabetes protection and susceptibility associated HLA class II genotypes on DNA methylation in immune cells

The HLA region, especially HLA class I and II genes, which encode molecules for antigen presentation to T cells, plays a major role in the predisposition to autoimmune disorders. To clarify the mechanisms behind this association, we examined genome-wide DNA methylation by microarrays to cover over 850,000 CpG sites in the CD4+ T cells and CD19+ B cells of healthy subjects homozygous either for DRB1*15-DQA1*01-DQB1*06:02 (DR2-DQ6, n = 14), associated with a strongly decreased T1D risk, DRB1*03-DQA1*05-DQB1*02 (DR3-DQ2, n = 19), or DRB1*04:01-DQA1*03-DQB1*03:02 (DR4-DQ8, n = 17), associated with a moderately increased T1D risk. In total, we discovered 14 differentially methylated CpG probes, of which 10 were located in the HLA region and six in the HLA-DRB1 locus. The main differences were between the protective genotype DR2-DQ6 and the risk genotypes DR3-DQ2 and DR4-DQ8, where the DR2-DQ6 group was hypomethylated compared to the other groups in all but four of the differentially methylated probes. The differences between the risk genotypes DR3-DQ2 and DR4-DQ8 were small. Our results indicate that HLA variants have few systemic effects on methylation and that their effect on autoimmunity is conveyed directly by HLA molecules, possibly by differences in expression levels or function.

Impact of genome build on RNA-seq interpretation and diagnostics

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network and Genomics Research to Elucidate the Genetics of Rare Disease Consortium. Across six routinely collected biospecimens, 61% of quantified genes were not influenced by genome build. However, we identified 1,492 genes with build-dependent quantification, 3,377 genes with build-exclusive expression, and 9,077 genes with annotation-specific expression across six routinely collected biospecimens, including 566 clinically relevant and 512 known OMIM genes. Further, we demonstrate that between builds for a given gene, a larger difference in quantification is well correlated with a larger change in expression outlier calling. Combined, we provide a database of genes impacted by build choice and recommend that transcriptomics-guided analyses and diagnoses are cross referenced with these data for robustness.

A rigorous benchmarking of alignment-based HLA typing algorithms for RNA-seq data

Accurate identification of human leukocyte antigen (HLA) alleles is essential for various clinical and research applications, such as transplant matching and drug sensitivities. Recent advances in RNA-seq technology have made it possible to impute HLA types from sequencing data, spurring the development of a large number of computational HLA typing tools. However, the relative performance of these tools is unknown, limiting the ability for clinical and biomedical research to make informed choices regarding which tools to use. Here we report the study design of a comprehensive benchmarking of the performance of 12 HLA callers across 682 RNA-seq samples from 8 datasets with molecularly defined gold standard at 5 loci, HLA-A, -B, -C, -DRB1, and -DQB1. For each HLA typing tool, we will comprehensively assess their accuracy, compare default with optimized parameters, and examine for discrepancies in accuracy at the allele and loci levels. We will also evaluate the computational expense of each HLA caller measured in terms of CPU time and RAM. We also plan to evaluate the influence of read length over the HLA region on accuracy for each tool. Most notably, we will examine the performance of HLA callers across European and African groups, to determine discrepancies in accuracy associated with ancestry. We hypothesize that RNA-Seq HLA callers are capable of returning high-quality results, but the tools that offer a good balance between accuracy and computational expensiveness for all ancestry groups are yet to be developed. We believe that our study will provide clinicians and researchers with clear guidance to inform their selection of an appropriate HLA caller.

Design of a novel multi-epitopes vaccine against Escherichia fergusonii: a pan-proteome based in- silico approach

Escherichia fergusonii a gram-negative rod-shaped bacterium in the Enterobacteriaceae family, infect humans, causing serious illnesses such as urinary tract infection, cystitis, biliary tract infection, pneumonia, meningitis, hemolytic uremic syndrome, and death. Initially treatable with penicillin, antibiotic misuse led to evolving resistance, including resistance to colistin, a last-resort drug. With no licensed vaccine, the study aimed to design a multi-epitope vaccine against E. fergusonii. The study started with the retrieval of the complete proteome of all known strains and proceeded to filter the surface exposed virulent proteins. Seventeen virulent proteins (4 extracellular, 4 outer membranes, 9 periplasmic) with desirable physicochemical properties were identified from the complete proteome of known strains. Further, these proteins were processed for B-cell and T-cell epitope mapping. Obtained epitopes were evaluated for antigenicity, allergenicity, solubility, MHC-binding, and toxicity and the filtered epitopes were fused by specific linkers and an adjuvant into a vaccine construct. Structure of the vaccine candidate was predicted and refined resulting in 78.1% amino acids in allowed regions and VERIFY3D score of 81%. Vaccine construct was docked with TLR-4, MHC-I, and MHC-II, showing binding energies of -1040.8 kcal/mol, -871.4 kcal/mol, and -1154.6 kcal/mol and maximum interactions. Further, molecular dynamic simulation of the docked complexes was carried out resulting in a significant stable nature of the docked complexes (high B-factor and deformability values, lower Eigen and high variance values) in terms of intermolecular binding conformation and interactions. The vaccine was also reported to stimulate a variety of immunological pathways after administration. In short, the designed vaccine revealed promising predictions about its immune protective potential against E. fergusonii infections however experimental validation is needed to validate the results.

Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.

CRISPR-based targeted haplotype-resolved assembly of a megabase region

Constructing high-quality haplotype-resolved genome assemblies has substantially improved the ability to detect and characterize genetic variants. A targeted approach providing readily access to the rich information from haplotype-resolved genome assemblies will be appealing to groups of basic researchers and medical scientists focused on specific genomic regions. Here, using the 4.5 megabase, notoriously difficult-to-assemble major histocompatibility complex (MHC) region as an example, we demonstrated an approach to construct haplotype-resolved assembly of the targeted genomic region with the CRISPR-based enrichment. Compared to the results from haplotype-resolved genome assembly, our targeted approach achieved comparable completeness and accuracy with reduced computing complexity, sequencing cost, as well as the amount of starting materials. Moreover, using the targeted assembled personal MHC haplotypes as the reference both improves the quantification accuracy for sequencing data and enables allele-specific functional genomics analyses of the MHC region. Given its highly efficient use of resources, our approach can greatly facilitate population genetic studies of targeted regions, and may pave a new way to elucidate the molecular mechanisms in disease etiology.

HLA allele-specific expression: Methods, disease associations, and relevance in hematopoietic stem cell transplantation

Varying HLA allele-specific expression levels are associated with human diseases, such as graft versus host disease (GvHD) in hematopoietic stem cell transplantation (HSCT), cytotoxic T cell response and viral load in HIV infection, and the risk of Crohn’s disease. Only recently, RNA-based next generation sequencing (NGS) methodologies with accompanying bioinformatics tools have emerged to quantify HLA allele-specific expression replacing the quantitative PCR (qPCR) -based methods. These novel NGS approaches enable the systematic analysis of the HLA allele-specific expression changes between individuals and between normal and disease phenotypes. Additionally, analyzing HLA allele-specific expression and allele-specific expression loss provide important information for predicting efficacies of novel immune cell therapies. Here, we review available RNA sequencing-based approaches and computational tools for NGS to quantify HLA allele-specific expression. Moreover, we explore recent studies reporting disease associations with differential HLA expression. Finally, we discuss the role of allele-specific expression in HSCT and how considering the expression quantification in recipient-donor matching could improve the outcome of HSCT.

A systematic review and meta-analysis of HLA class II associations in patients with IgG4 autoimmunity

Autoimmune diseases caused by pathogenic IgG4 subclass autoantibodies (IgG4-AID) include diseases like MuSK myasthenia gravis, pemphigus vulgaris or thrombotic thrombocytopenic purpura. Their etiology is still unknown. Polymorphisms in the human leukocyte antigen (HLA) gene locus, particularly in HLA-DRB1, are known genetic susceptibility factors for autoimmune diseases. We hypothesized a similar role for HLA polymorphisms in IgG4-AID and conducted a systematic review and meta-analysis with case-control studies on IgG4-AID based on MOOSE/ HuGENet guidelines. Genotype (G) and allele (A) frequencies of HLA-DQB1*05 (G: OR 3.8; 95% CI 2.44-5.9; p < 0.00001; A: OR 2.54; 95% CI 1.82-3.55; p < 0.00001) and HLA-DRB1*14 (G: OR 4.31; 95% CI 2.82-6.59; p < 0.00001; A: OR 4.78; 95% CI 3.52-6.49; p < 0.00001) and the HLA-DRB1*14-DQB1*05 haplotype (OR 6.3; 95% CI 3.28-12.09; p < 0.00001/OR 4.98; 95% CI 3.8-6.53; p < 0.00001) were increased while HLA-DRB1*13 (G: OR 0.48; 95% CI 0.34-0.68; p < 0.0001; A: OR 0.46; 95% CI 0.34-0.62; p < 0.00001) was decreased in IgG4-AID patients. In conclusion, the HLA-DQB1*05, HLA-DRB1*14 alleles and the HLA-DQB1*05-DRB1*14 haplotype could be genetic risk factors that predispose for the production of pathogenic IgG4 autoantibodies and the HLA-DRB1*13 allele may protect from IgG4 autoimmunity.

Bioinformatics and System Biology Approach to Identify the Influences of COVID-19 on Rheumatoid Arthritis

Background

Severe coronavirus disease 2019 (COVID -19) has led to a rapid increase in mortality worldwide. Rheumatoid arthritis (RA) was a high-risk factor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, whereas the molecular mechanisms underlying RA and CVOID-19 are not well understood. The objectives of this study were to analyze potential molecular mechanisms and identify potential drugs for the treatment of COVID-19 and RA using bioinformatics and a systems biology approach.

Methods

Two Differentially expressed genes (DEGs) sets extracted from GSE171110 and GSE1775544 datasets were intersected to generate common DEGs, which were used for functional enrichment, pathway analysis, and candidate drugs analysis.

Results

A total of 103 common DEGs were identified in the two datasets between RA and COVID-19. A protein-protein interaction (PPI) was constructed using various combinatorial statistical methods and bioinformatics tools. Subsequently, hub genes and essential modules were identified from the PPI network. In addition, we performed functional analysis and pathway analysis under ontological conditions and found that there was common association between RA and progression of COVID-19 infection. Finally, transcription factor-gene interactions, protein-drug interactions, and DEGs-miRNAs coregulatory networks with common DEGs were also identified in the datasets.

Conclusion

We successfully identified the top 10 hub genes that could serve as novel targeted therapy for COVID-19 and screened out some potential drugs useful for COVID-19 patients with RA.