Next Article in Journal
Ex Vivo Human Histology Fractional Treatment with a New CO2 Scanner: A Potential Application on Deep Scarring
Previous Article in Journal
Angiogenesis and Hepatocellular Carcinoma: From Molecular Mechanisms to Systemic Therapies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detection of Nonsynonymous Single Variants in Human HLA-DRB1 Exon 2 Associated with Renal Transplant Rejection

1
Department of Hematology, Faculty of Medical Laboratory Sciences, National University, Khartoum 11111, Sudan
2
Department of Pharmaceutical Microbiology, Faculty of Pharmacy, International University of Africa, Khartoum 11111, Sudan
3
Department of Parasitology and Medical Entomology, Faculty of Medical Laboratory Sciences, National University, Khartoum11111, Sudan
4
Department of Chemical Pathology, Faculty of Medical Laboratory Sciences, National University, Khartoum 11111, Sudan
5
Molecular Biology Unit, Sirius Training and Research Centre, Khartoum 11111, Sudan
6
Department of Molecular Biology, National University Biomedical Research Institute, National University, Khartoum 11111, Sudan
7
Department of Pharmacognosy, Faculty of Pharmacy, Prince Sattam Bin Abdulaziz University, Al-kharj 11942, Saudi Arabia
8
Department of Pharmacognosy, Faculty of Pharmacy, Mansoura University, Mansoura 35516, Egypt
9
Department of Pharmacognosy, Faculty of Pharmacy, University of Khartoum, Al-Qasr Ave, Khartoum 11111, Sudan
10
Preparatory Year Program, Department of Chemistry, Batterjee Medical College, Jeddah 21442, Saudi Arabia
11
Department of Pharmacognosy, Faculty of Pharmacy, Assiut University, Assiut 71526, Egypt
12
Clinical Pharmacy Department, College of Pharmacy, Umm Al-Qura University, Makkah 24382, Saudi Arabia
13
Department of Pharmacology and Toxicology, College of Pharmacy, Taibah University, Al-Madinah Al-Munawwarah 30078, Saudi Arabia
14
Department of Natural Products and Alternative Medicine, Faculty of Pharmacy, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Medicina 2023, 59(6), 1116; https://doi.org/10.3390/medicina59061116
Submission received: 14 April 2023 / Revised: 1 June 2023 / Accepted: 6 June 2023 / Published: 9 June 2023
(This article belongs to the Section Urology & Nephrology)

Abstract

:
Background: HLA-DRB1 is the most polymorphic gene in the human leukocyte antigen (HLA) class II, and exon 2 is critical because it encodes antigen-binding sites. This study aimed to detect functional or marker genetic variants of HLA-DRB1 exon 2 in renal transplant recipients (acceptance and rejection) using Sanger sequencing. Methods: This hospital-based case-control study collected samples from two hospitals over seven months. The 60 participants were equally divided into three groups: rejection, acceptance, and control. The target regions were amplified and sequenced by PCR and Sanger sequencing. Several bioinformatics tools have been used to assess the impact of non-synonymous single-nucleotide variants (nsSNVs) on protein function and structure. The sequences data that support the findings of this study with accession numbers (OQ747803-OQ747862) are available in National Center for Biotechnology Information (GenBank database). Results: Seven SNVs were identified, two of which were novel (chr6(GRCh38.p12): 32584356C>A (K41N) and 32584113C>A (R122R)). Three of the seven SNVs were non-synonymous and found in the rejection group (chr6(GRCh38.p12): 32584356C>A (K41N), 32584304A>G (Y59H), and 32584152T>A (R109S)). The nsSNVs had varying effects on protein function, structure, and physicochemical parameters and could play a role in renal transplant rejection. The chr6(GRCh38.p12):32584152T>A variant showed the greatest impact. This is because of its conserved nature, main domain location, and pathogenic effects on protein structure, function, and stability. Finally, no significant markers were identified in the acceptance samples. Conclusion: Pathogenic variants can affect intramolecular/intermolecular interactions of amino acid residues, protein function/structure, and disease risk. HLA typing based on functional SNVs could be a comprehensive, accurate, and low-cost method for covering all HLA genes while shedding light on previously unknown causes in many graft rejection cases.

1. Introduction

Over 800 million people worldwide (10% of the population) have chronic renal disease (CRD) [1]. CRD is more common in older people, women, and people with diabetes and high blood pressure [2]. Low- and middle-income countries face a significant burden of CRD [3,4]. Chronic renal disease is one of the leading causes of death worldwide [1]. During the period 1990–2017, CRD mortality increased by 41.5% globally [5]. Over one million patients were predicted to have end-stage renal disease (ESRD) globally two decades ago, with a 7% annual increase [6]. In Sudan, the prevalence of CRD ranges from 7.7 to 11%, with an estimated incidence of new cases of 70–140/million per year [7,8]. Further, 1000 new patients are diagnosed with ESRD each year, and the most common cause of renal failure in Sudanese people (53.6%) is unknown [9,10]. In general, the number of patients with ESRD (who require dialysis or renal transplant to survive) is increasing, and it is becoming a major public health concern worldwide [11].
In 2007, there were over 1.6 million dialysis patients and half a million renal transplant recipients worldwide [12]. Renal transplantation remains the most effective CRD/ESRD treatment, accounting for 28% of total renal therapy in Sudan [12,13]. Transplant rejection can be hyperacute (minutes to hours), acute (days to weeks), or chronic (months to years) [14]. The International Society of Nephrology analyzed data from 182 countries and reported a rejection rate of 59% [15]. Another study conducted in Iran identified the clinical causes of renal allograft nephrectomy, with chronic rejection (38%) being the most common cause [16,17]. Furthermore, a study performed between November 2011 and 2015 at Sharg El-Neel Hospital in Khartoum, Sudan, discovered that the rate of acute rejection was 10.4% [9]. Human leukocyte antigen (HLA) typing is an important step in transplantation, and a well-matched donor is critical for successful transplantation [17].
HLA genes are located on chromosome 6p (short arm) in the distal portion of the 21.3 band, one of the most polymorphic and gene-dense regions [18,19]. HLA complex genes and their protein products are divided into three classes based on their tissue distribution, structure, and function [20]. MHC class II antigens encoded by the HLA-DM, -DO, -DP, -DQ, and -DR loci, and their products are included in the immunoglobulin supergene family [21]. HLA-DR is a heterodimer comprising an alpha chain (DRA) and a beta chain (DRB) [19]. According to the IPD-IMGT/HLA database, HLA-DRB1 is the most polymorphic in class II of this system, with 3298 alleles in September 2022 (https://www.ebi.ac.uk/ipd/imgt/hla/about/statistics/) (accessed on 25 September 2022). The HLA-DRB1 gene is located in GRCh38.p12 (Genome Reference Consortium Human Build 38.p12) coordinates 32,578,775 to 32,589,848, has five introns, and is encoded by six exons [22]. Exon one encodes the leader peptide, exons two and three encode the two extracellular domains, exon four encodes the transmembrane domain, and exon five encodes the cytoplasmic tail (https://www.ncbi.nlm.nih.gov/gene/3123) (accessed on 25 September 2022 [23]. Many studies emphasized the importance of exon 2 as it encodes antigen-binding sites, contains the most pathological single nucleotide variants (SNVs), and is commonly included in high throughput HLA-typing commercial kits [24,25,26,27].
Several benefits are associated with HLA matching in organ transplants, such as kidneys, including improved graft function, reduced the incidence of acute or chronic rejection, extended graft survival, and the potential for reduced immunosuppression [28]. It was reported that patient-donor matching of HLA determinants lowers the risks of chronic and acute GVHD (graft-versus-host disease) [29]. Early studies indicated that HLA-DRB1 mismatch is a particular risk factor for rejection and is critical in the first six months after transplantation [24,27,30,31]. As a result of the realization that HLA plays a significant role in transplantation, the use of HLA typing in transplantation has seen numerous advancements [32,33]. Owing to this development, HLA typing has progressed from identifying HLA proteins to identifying HLA gene variations [31]. The 1000 Genomes Project provides an in-depth analysis of common genetic variations (single nucleotide polymorphisms (SNPs) and Insertions–deletions (indels)) in humans and their association with diseases [34]. HLA variants are strongly linked to various diseases and organ transplantation [22,31,35,36]. SNPs are single nucleotide variants (SNVs) in DNA sequences with a population allele frequency of 1% or higher [37,38]. SNVs can be found in both the coding and non-coding regions of the human genome [39]. Non-synonymous SNVs are a type of single variant that represent amino acid substitutions and protein variations [40]. Previous studies have indicated that nsSNPs account for approximately half of the mutations involved in various genetic diseases [41]. Indels are another type of significant genomic variant that are insertions or deletions of one or more DNA nucleotides [42]. The current study aimed to identify functional or marker genetic variants within HLA-DRB1 exon 2 in patients with renal transplant status (acceptance and rejection).

2. Materials and Methods

2.1. Study Design and Samples Information

This hospital-based case-control study was conducted at Ahmed Gasim and Ibn Sena Hospitals. Blood samples were collected from March to September 2021 using a convenience sampling method. Samples were collected from individuals of any age and sex who had received a renal transplant, regardless of whether they developed graft rejection (acute, hyperacute, or chronic) or acceptance within the first six months. Participants who had their renal transplants rejected because of medical errors or negligence were excluded. The total sample size was 60, divided equally into three groups. The first group included participants who had graft acceptance for more than the first six months. The second group of participants had graft rejection in the first six months, and the third group was the control group. The study was conducted in accordance with the guidelines of the Declaration of Helsinki and approved by the Institutional Ethics Committee of the National University-Sudan under approval No. NU170220219 (date of approval: 17 February 2021). Written informed consent was obtained from the patient(s) for their anonymized information to be published in this article.

2.2. DNA Isolation, Amplification, and Sequencing

Genomic DNA was extracted from the blood samples using QiaAmp blood extraction kits according to the manufacturer’s instructions (Qiagen, Hilden, Germany). The extracted DNA was tested for quality using a NanoDrop spectrophotometer (Implen, München, Germany) and stored at −20 °C until molecular analysis.
The following primers were used to amplify the human HLA-DRB1 gene (target region): Forward primer: 5′GTG CTC TCA GAA CTG CTT GC 3′, and reverse primer: 5′ CCT CAG GAA GAC GGA GGA TGA 3′. The PCR reaction mixture contained 1 µL of the extracted DNA added to 4 µL PCR master mix (Solis Biodyne, Tartu, Estonia) containing 1 U DNA polymerase, 12.5 mM MgCl2, and 4 mM dNTPs. The PCR thermal conditions were as follows: Initial denaturation at 94 °C for 5 min, followed by 35 cycles of denaturation at 94 °C for 30 s, annealing at 57 °C for 30 s, extension at 72 °C for 30 s, and a final extension step at 72 °C for 10 min. The thermal conditions were determined using a 2721 Thermocycler (Applied Biosystems, Thermo Fisher Scientific, Budapest, Hungary). Following PCR, the amplicons were visualized using 2% gel electrophoresis (Major Sciences, London, United Kingdom) by applying the PCR product to an electrical current adjusted to 100 V and 70 A for one h.
The amplified PCR amplicons were sequenced in duplicate based on both directions’ primers by the Sanger dideoxynucleotide chain-termination sequencing method using a 3730XL DNA analyzer (Applied Biosystems, Waltham, MA, USA) by Macrogen (Macrogen Inc., Amsterdam, The Netherlands).

2.3. Sequences and Variants Analysis Using Bioinformatics

DNA sequencing results of the 60 samples were obtained as AB1 files. Initially, the Chromatogram Explorer program (version 5.0.2.3) was used to assess the overall quality of the sequences (read Phred quality score), trim low-quality ends, and convert AB1 to FASTA formats [43]. Subsequently, the Basic Local Alignment Search Tool (BLAST) algorithm was used to check the specificity of these sequences by comparing them to the Homo sapiens genome (GRCh38.p12) using the Ensembl genome browser (https://www.ensembl.org/index.html) (accessed on 30 September 2022) [44]. Exon 2 regions were manually extracted from the high-quality sequences for further analysis. All sequences, including the reference (target region), were manually prepared and submitted for Multiple Sequence Alignment (MSA) using CLC genomics workbench program version 21.0.5. CLC is a Qiagen-bioinformatics commercial analysis and visualization product (https://www.qiagenbioinformatics.com/) (accessed on 5 October 2022). The gap cost parameters of the alignment algorithm were as follows: Gap open 10.0, gap extension 1.0, and end gap cost. To generate MSA, the CLC employs a progressive alignment algorithm [45]. Organizing sequence data in MSAs can reveal conserved and variable sites (variants or mutations) [46]. Variants were manually extracted from MSA and prepared for further analysis.
The chromosomal location of the detected variants was initially submitted to the Ensembl Variant Effect Predictor (VEP) [47]. VEP can annotate, analyze, and prioritize genomic variants in both coding and non-coding regions. VEP was used to determine the variants’ availability, frequency, and amino acid positions. Non-synonymous SNVs were then submitted sequentially to SIFT [48], PolyPhen-2 [49], PredictSNP [50], PANTHER [51], SNP&GO [52], SNAP2 [53], and PhD-SNP [54] tools to differentiate between functional (deleterious) and non-functional nsSNVs. The I-mutant server was used to determine whether nsSNVs affected protein stability [55]. The HOPE server was used to evaluate the effects of nsSNVs on protein structure [56]. The location of the domain and high evolutionary conservation were then determined using the InterPro and Consurf servers [57,58]. Moreover, the ProtParam server was used to assess the impact of nsSNVs on protein physicochemical parameters [59]. Finally, the STRING database (Version 11.5) was used to predict associations between HLA-DRB1 and most related proteins to construct a protein-protein network based on physical interactions and functional associations [60]. The study methodology is summarized in Figure 1.
The following are websites for the previous tools: VEP https://www.ensembl.org/info/docs/tools/vep/index.html, (accessed on 7 October 2022), SIFT https://sift.bii.a-star.edu.sg/, (accessed on 7 October 2022). PolyPhen-2 http://genetics.bwh.harvard.edu/pph2/, (accessed on 9 October 2022), PredictSNP tool https://loschmidt.chemi.muni.cz/predictsnp/, (accessed on 12 October 2022), PANTHER http://www.pantherdb.org/tools/csnpScoreForm.jsp, (accessed on 15 October 2022), SNP&GO https://snps.biofold.org/snps-and-go/snps-and-go.html, (accessed on 20 October 2022), SNAP2 https://www.rostlab.org/services/snap/, (accessed on 23 October 2022), PhD-SNP https://snps.biofold.org/phd-snp/phd-snp.html, (accessed on 27 October 2022), I-mutant v3.0 http://gpcr2.biocomp.unibo.it/cgi/predictors/I-Mutant3.0/I-Mutant3.0.cgi, (accessed on 1 November 2022), HOPE https://www3.cmbi.umcn.nl/hope/, (accessed on 5 November 2022), InterPro database https://www.ebi.ac.uk/interpro/, (accessed on 10 November 2022), Consurf https://consurf.tau.ac.il/, (accessed on 13 November 2022) ProtParam https://web.expasy.org/protparam/, (accessed on 17 November 2022), STRING database https://string-db.org/ (accessed on 20 November 2022).

3. Results

This study used DNA from 40 samples (acceptance and rejection) and 20 controls to target the HLA-DRB1 exon 2. BLAST revealed that all DNA quality-checked sequences showed high similarity (>99%) and specificity for the HLA-DRB1 target region. Multiple sequence alignments revealed seven SNVs in ten samples and controls. The Ensembl variant effect predictor was used to collect broad information on seven SNVs, two of which were novel (Table 1).
Three of the seven detected SNVs were non-synonymous and were only located in the rejected samples (R3, R9, and R16), whereas the remaining were synonymous (Table 1). To identify the deleterious effects of nsSNVs at the functional level, seven different tools (SIFT, PolyPhen, PredictSNP, PANTHER, SNP&GO, SNAP2, and PhD-SNP) with different prediction algorithms were used. Two nsSNVs (K41N and R109S) were predicted to be pathogenic by all seven tools, whereas the third (Y59H) was predicted by only four (Table 2). Furthermore, the I-mutant server predicted that all nsSNVs would affect the protein stability (Table 2).
At the structural level (using the HOPE server), the new residues differed in size, charge, and hydrophobicity. The new residues also influenced hydrogen bond formation, ionic interactions, multimer contacts or interactions, and the function of their region (Table 3 and Figure 2).
The Consurf server and InterPro database were used to predict the locations of variants in evolutionarily conserved and domain regions. The three nsSNVs K41N, Y59H, and R109S received scores of 5, 1, and 6, respectively, indicating that they were average, variable, and conserved, respectively (Figure 3).
Additionally, two of the three nsSNVs (Y59H and R109S) were discovered in the MHC II b N domain (MHC class II, beta chain, N-terminal) with the accession number IPR000353. All nsSNVs demonstrated changes in the overall protein physicochemical parameters. The properties altered by all three nsSNVs were the molecular weight, theoretical isoelectric point (pI), atomic composition, instability index, and GRAVY (Table 4).
HLA-DRB1 interacts with HLA-DRA, HLA-DMA, CD74, HLA-DMB, HLA-DPA1, CD4, BTNL2, and CD86, in that order, according to the protein-protein interaction network (Figure 4). Finally, no significant variants that had an impact on or worked as markers for transplant acceptance samples were found (Table 1).

4. Discussion

Chronic renal diseases in general and end-stage renal disease in particular are major health concerns in Sudan and around the world [1,6,7,8,9]. The growing number of ESRD patients places strain on both individuals (costs of dialysis or transplantation) and governments (increasing the financial burden of health care) [61]. Renal transplantation remains the most effective ESRD treatment, and HLA typing is the most important test in this process [13,17]. The HLA region is extremely diverse, and HLA-DRB1 is the most polymorphic in class II of this system [19,62]. HLA-DRB1 protein is significantly associated with graft survival, particularly in the first six months after transplantation [27,28]. HLA-DRB1 exon 2 is important because it encodes antigen-binding sites and contains the most pathological single nucleotide variants (SNVs) [24,25,26,27]. The present study aimed to identify functional or marker genetic variants in HLA-DRB1 exon 2 in renal transplant recipients (acceptance and rejection) using sequencing technology.
This study included 60 DNA samples (from various families) divided into three equal groups: Control, renal transplant rejection, and acceptance. The alignment algorithm identified seven SNVs at six locations, three of which were non-synonymous and could have functional consequences. Two of the three nsSNVs were found in the public domain archive of simple genetic polymorphisms (https://www.ncbi.nlm.nih.gov/snp/), (accessed on 20 December 2022), whereas the third was not found. Several algorithms and methods have been used to assess the impact of these three nsSNVs on protein function and stability indices. They discovered that two were high-impact, indicating deleterious effects by all tools, and the third had a lower impact. Previous (the preceding three decades) and recent studies have shown a close relationship between HLA mismatch and transplantation graft rejection [63,64]. Previous findings have indicated that R109S and K41N variants may have a greater impact on protein function than the Y59H variant. None of the three detected variants corresponded to an in silico study conducted in 2014 by Hassan et al., who identified the most deleterious nsSNVs in the human HLA-DRB1 gene [25]. The nsSNVs discovered in this study were submitted to dbSNP after Hassan et al.’s study, resulting in incompatibility. https://www.ncbi.nlm.nih.gov/snp/rs750986830#submissions, https://www.ncbi.nlm.nih.gov/snp/rs11554462#submissions (accessed on 5 April 2023).
At the structural level, all the variant amino acids showed differences in physical properties, bond formation, and a variety of interactions. The HOPE server predicted that the three nsSNVs are situated in a unique region called Beta-1. According to the Universal Protein Resource, the beta-1 region is a structural part of the peptide-binding cleft of HLA-DRB1, consisting of 94 amino acids. Additionally, beta-1 interacts with the T-cell receptors CDR2 and CDR3 (complementarity-determining regions 2 and 3) alpha domains through hydrogen bonds. https://www.uniprot.org/uniprotkb/P01911/entry (accessed on 15 December 2022). Differences in amino acid properties can disrupt this region (beta-1) and its function. According to evolutionary conservation, the K41N, Y59H, and R109S variants were average, variable, and conserved, respectively. Two of these (Y59H and R109S) were found in the MHC II b N domain. Compared to other variants, the R109S variant may have the greatest impact. This is because of its conserved nature, main domain location, and pathogenic impact on protein structure, function, and stability.
In terms of physicochemical properties, the HOPE tool, as previously mentioned, revealed differences in the residue levels, whereas ProtParam indicated that the variants caused changes in the entire protein. All the nsSNVs detected agreed to alter the molecular weight, theoretical pi (isoelectric point), atomic composition, and GRAVY of the protein. The extinction coefficients and total positive charge were altered in a few variants, but the aliphatic index and total negative charge remained the same. In general, all the nsSNVs had nearly equal effects on the overall physicochemical properties of the protein. Most proteins function consecutively with other proteins in living organisms, and protein-protein interaction studies provide crucial information for understanding the complicated biological processes that occur in live cells [65,66]. Thus, to gain a better understanding, a network of protein-protein interactions (PPIs) was constructed using the STRING database. Deleterious variants in the HLA-DRB1 protein could disrupt its interaction with confidence interaction proteins. Only one synonymous SNV (K41K) was found in the transplant acceptance samples in the current study. The K41K variant could not be identified as an acceptance marker because it appeared in both acceptance and rejection cases.

Author Contributions

Conceptualization, M.M.H.; Methodology, M.M.H., M.A.H., S.S.A., N.S.M., H.A., O.M. and M.A.M.; Formal analysis, M.M.H.; Investigation, M.M.H., N.S.M., H.A. and O.M.; Supervision, M.A.H., S.S.A., W.O. and M.A.M.; Validation, M.A.H., S.S.A., A.E.S., A.A. and M.A.M.; Visualization, M.M.H.; Writing—original draft, M.M.H., M.A.H. and N.S.M.; Writing—review & editing, S.R.M.I., K.F.G., S.F.M., G.A.M., A.A., A.E.S., W.O., S.S.A. and M.A.M.; Resources: S.R.M.I., K.F.G., S.F.M. and G.A.M.; Project administration, M.M.H.; Funding acquisition, W.O. All authors have read and agreed to the published version of the manuscript.

Funding

This study is supported via funding from Prince Sattam bin Abdulaziz University. Project number (PSAU/2023/R/1444).

Institutional Review Board Statement

The study was conducted in accordance with the guidelines of the Declaration of Helsinki and approved by the Institutional Ethics Committee of the National University-Sudan (NU170220219; date of approval: 17 February 2021).

Data Availability Statement

The sequences data that support the findings of this study are available in [National Center for Biotechnology Information-genebank database] at [https://www.ncbi.nlm.nih.gov/nucleotide/, accessed on 10 April 2023], accession numbers [OQ747803-OQ747862].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kovesdy, C.P. Epidemiology of chronic kidney disease: An update 2022. Kidney Int. Suppl. 2022, 12, 7–11. [Google Scholar] [CrossRef] [PubMed]
  2. Jager, K.J.; Kovesdy, C.; Langham, R.; Rosenberg, M.; Jha, V.; Zoccali, C. A single number for advocacy and communication-worldwide more than 850 million individuals have kidney diseases. Kidney Int. 2019, 96, 1048–1050. [Google Scholar] [CrossRef] [PubMed]
  3. Hill, N.R.; Fatoba, S.T.; Oke, J.L.; Hirst, J.A.; O’Callaghan, C.A.; Lasserson, D.S.; Hobbs, F.D. Global Prevalence of Chronic Kidney Disease—A Systematic Review and Meta-Analysis. PLoS ONE 2016, 11, e0158765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Xie, Y.; Bowe, B.; Mokdad, A.H.; Xian, H.; Yan, Y.; Li, T.; Maddukuri, G.; Tsai, C.Y.; Floyd, T.; Al-Aly, Z. Analysis of the Global Burden of Disease study highlights the global, regional, and national trends of chronic kidney disease epidemiology from 1990 to 2016. Kidney Int. 2018, 94, 567–581. [Google Scholar] [CrossRef] [Green Version]
  5. GBD Chronic Kidney Disease Collaboration. Global, regional, and national burden of chronic kidney disease, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet 2020, 395, 709–733. [Google Scholar] [CrossRef] [Green Version]
  6. Lysaght, M.J. Maintenance dialysis population dynamics: Current trends and long-term implications. J. Am. Soc. Nephrol. 2002, 13 (Suppl. S1), S37–S40. [Google Scholar] [CrossRef]
  7. Abu-Aishaa, H.; Elhassanb, E.A.M.; Khamisc, A.H.; Abu-Elmaali, A. Chronic Kidney Disease in Police Forces Households in Khartoum, Sudan: Pilot Report. Arab. J. Nephrol. Transpl. 2009, 2, 21–26. [Google Scholar] [CrossRef] [Green Version]
  8. Suliman, S.M.; Beliela, M.H.; Hamza, H. Dialysis and transplantation in Sudan. Saudi J. Kidney Dis. Transpl. 1995, 6, 312–314. [Google Scholar]
  9. Ahmed, I.A.M.; Musa, E.A.; Shaikh, Q.; Elsharif, M.E. Kidney Transplantation in Sudan. Transplant. 2018, 102, 1583–1585. [Google Scholar] [CrossRef]
  10. Elsharif, M.E.; Elsharif, E.G. Causes of end-stage renal disease in Sudan: A single-center experience. Saudi J. Kidney Dis. Transpl. 2011, 22, 373–376. [Google Scholar]
  11. Parmar, M.S. Chronic renal disease. BMJ 2002, 325, 85–90. [Google Scholar] [CrossRef] [PubMed]
  12. Elamin, S.; Obeid, W.; Abu-Aisha, H. Renal Replacement Therapy in Sudan, 2009. Arab. J. Nephrol. Transpl. 2010, 3, 31–36. [Google Scholar] [CrossRef] [Green Version]
  13. Banaga, A.; Mohammed, E.; Siddig, R.; Salama, D.; Elbashir, S.; Khojali, M.; Babiker, R.; Elmusharaf, K.; Homeida, M.M. Why Did Sudanese End Stage Renal Failure Patients Refuse Renal Transplantation? Open J. Nephrol. 2015, 5, 35–39. [Google Scholar] [CrossRef] [Green Version]
  14. Justiz Vaillant, A.A.; Mohseni, M. Chronic Transplantation Rejection. In StatPearls; StatPearls Publishing: Tampa, FL, USA, 2022. [Google Scholar]
  15. Mudiayi, D.; Shojai, S.; Okpechi, I.; Christie, E.A.; Wen, K.; Kamaleldin, M.; Elsadig, O.M.; Lunney, M.; Prasad, B.; Osman, M.A.; et al. Global Estimates of Capacity for Kidney Transplantation in World Countries and Regions. Transplantation 2022, 106, 1113–1122. [Google Scholar] [CrossRef] [PubMed]
  16. Panahi, A.; Bidaki, R.; Mirhosseini, S.M.; Mehraban, D. Renal allograft nephrectomy: Comparison between clinical and pathological diagnosis. Nephrourol. Mon. 2013, 5, 1001–1004. [Google Scholar] [CrossRef] [Green Version]
  17. Law, S.C.; Haigh, O.L.; Walpole, C.M.; Keane, C.; Miles, J.J.; Gandhi, M.K.; Radford, K.J.; Steptoe, R.J. Simple, rapid and inexpensive typing of common HLA class I alleles for immunological studies. J. Immunol. Methods 2019, 465, 72–76. [Google Scholar] [CrossRef] [Green Version]
  18. Choo, S.Y. The HLA system: Genetics, immunology, clinical testing, and clinical implications. Yonsei Med. J. 2007, 48, 11–23. [Google Scholar] [CrossRef] [Green Version]
  19. Mehra, N.K.; Kaur, G. Histocompatibility Antigen Complex of Man. In eLS; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2016; pp. 1–8. [Google Scholar] [CrossRef]
  20. Medhasi, S.; Chantratita, N. Human Leukocyte Antigen (HLA) System: Genetics and Association with Bacterial and Viral Infections. J. Immunol. Res. 2022, 2022, 9710376. [Google Scholar] [CrossRef]
  21. Shiina, T.; Hosomichi, K.; Inoko, H.; Kulski, J.K. The HLA genomic loci map: Expression, interaction, diversity and disease. J. Hum. Genet. 2009, 54, 15–39. [Google Scholar] [CrossRef] [Green Version]
  22. Hassan, M.; Mohamed, S.; Hussain, M.; Dowd, A. Deleterious Nonsynonymous SNP Found within HLA-DRB1 Gene Involved in Allograft Rejection in Sudanese Family: Using DNA Sequencing and Bioinformatics Methods. Open J. Immunol. 2015, 5, 222–232. [Google Scholar] [CrossRef] [Green Version]
  23. Sayers, E.W.; Bolton, E.E.; Brister, J.R.; Canese, K.; Chan, J.; Comeau, D.C.; Connor, R.; Funk, K.; Kelly, C.; Kim, S.; et al. Database resources of the national center for biotechnology information. Nucleic Acids. Res. 2022, 50, D20–D26. [Google Scholar] [CrossRef] [PubMed]
  24. Fürst, D.; Neuchel, C.; Tsamadou, C.; Schrezenmeier, H.; Mytilineos, J. HLA Matching in Unrelated Stem Cell Transplantation up to Date. Transfus. Med. Hemother. 2019, 46, 326–336. [Google Scholar] [CrossRef] [PubMed]
  25. Hassan, M.M.; Dowd, A.A.; Mohamed, A.H.; Mahalah, S.M.; Kaheel, H.H.; Mohamed, S.N.; Hassan, M.A. Computational analysis of deleterious nsSNPs within HLA-DRB1 and HLA-DQB1 genes responsible for Allograft rejection. Int. J. Comput. Bioinform. Silico Model. 2014, 3, 562–577. [Google Scholar]
  26. Baek, I.C.; Choi, E.J.; Shin, D.H.; Kim, H.J.; Choi, H.; Kim, T.G. Allele and haplotype frequencies of human leukocyte antigen-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 by next generation sequencing-based typing in Koreans in South Korea. PLoS ONE 2021, 16, e0253619. [Google Scholar] [CrossRef]
  27. Mahdi, B.M. A glow of HLA typing in organ transplantation. Clin. Transl. Med. 2013, 2, 2–6. [Google Scholar] [CrossRef] [Green Version]
  28. Zachary, A.A.; Leffell, M.S. HLA Mismatching Strategies for Solid Organ Transplantation—A Balancing Act. Front. Immunol. 2016, 7, 575. [Google Scholar] [CrossRef]
  29. Petersdorf, E.W. Which factors influence the development of GVHD in HLA-matched or mismatched transplants? Best Pract. Res. Clin. Haematol. 2017, 30, 333–335. [Google Scholar] [CrossRef]
  30. Vu, L.T.; Baxter-Lowe, L.A.; Garcia, J.; McEnhill, M.; Summers, P.; Hirose, R.; Lee, M.; Stock, P.G. HLA-DR matching in organ allocation: Balance between waiting time and rejection in pediatric kidney transplantation. Arch. Surg. 2011, 146, 824–829. [Google Scholar] [CrossRef] [Green Version]
  31. Edgerly, C.H.; Weimer, E.T. The Past, Present, and Future of HLA Typing in Transplantation. In Methods in Molecular Biology; Humana Press: New York, NY, USA, 2018; Volume 1802. [Google Scholar]
  32. Erlich, H. HLA DNA typing: Past, present, and future. Tissue Antigens 2012, 80, 1–11. [Google Scholar] [CrossRef]
  33. Cereb, N.; Kim, H.R.; Ryu, J.; Yang, S.Y. Advances in DNA sequencing technologies for high resolution HLA typing. Hum. Immunol. 2015, 76, 923–927. [Google Scholar] [CrossRef]
  34. 1000 Genomes Project Consortium. A global reference for human genetic variation. Nature 2015, 526, 68–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Dendrou, C.A.; Petersen, J.; Rossjohn, J.; Fugger, L. HLA variation and disease. Nat. Rev. Immunol. 2018, 18, 325–339. [Google Scholar] [CrossRef] [PubMed]
  36. Larjo, A.; Eveleigh, R.; Kilpeläinen, E.; Kwan, T.; Pastinen, T.; Koskela, S.; Partanen, J. Accuracy of Programs for the Determination of Human Leukocyte Antigen Alleles from Next-Generation Sequencing Data. Front. Immunol. 2017, 8, 1815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. He, Q.; He, Q.; Liu, X.; Wei, Y.; Shen, S.; Hu, X.; Li, Q.; Peng, X.; Wang, L.; Yu, L. Genome-wide prediction of cancer driver genes based on SNP and cancer SNV data. Am. J. Cancer Res. 2014, 4, 394–410. [Google Scholar] [PubMed]
  38. Zou, H.; Wu, L.X.; Tan, L.; Shang, F.F.; Zhou, H.H. Significance of Single-Nucleotide Variants in Long Intergenic Non-protein Coding RNAs. Front. Cell Dev. Biol. 2020, 8, 347. [Google Scholar] [CrossRef]
  39. Wu, J.; Li, Y.; Rendahl, A.; Bhargava, M. Novel Human FCGR1A Variants Affect CD64 Functions and Are Risk Factors for Sarcoidosis. Front. Immunol. 2022, 13, 841099. [Google Scholar] [CrossRef]
  40. Cline, M.S.; Karchin, R. Using bioinformatics to predict the functional impact of SNVs. Bioinformatics 2011, 27, 441–448. [Google Scholar] [CrossRef] [Green Version]
  41. Hossain, M.S.; Roy, A.S.; Islam, M.S. In silico analysis predicting effects of deleterious SNPs of human RASSF5 gene on its structure and functions. Sci. Rep. 2020, 10, 14542. [Google Scholar] [CrossRef]
  42. Lin, M.; Whitmire, S.; Chen, J.; Farrel, A.; Shi, X.; Guo, J.T. Effects of short indels on protein structure and function in human genomes. Sci. Rep. 2017, 7, 9313. [Google Scholar] [CrossRef] [Green Version]
  43. Alachiotis, N.; Vogiatzi, E.; Pavlidis, P.; Stamatakis, A. ChromatoGate: A Tool for Detecting Base Mis-Calls in Multiple Sequence Alignments by Semi-Automatic Chromatogram Inspection. Comput. Struct. Biotechnol. J. 2013, 6, e201303001. [Google Scholar] [CrossRef] [Green Version]
  44. Cunningham, F.; Allen, J.E.; Allen, J.; Alvarez-Jarreta, J.; Amode, M.R.; Armean, I.M.; Austine-Orimoloye, O.; Azov, A.G.; Barnes, I.; Bennett, R.; et al. Ensembl 2022. Nucleic. Acids Res. 2022, 50, D988–D995. [Google Scholar] [CrossRef] [PubMed]
  45. Feng, D.F.; Doolittle, R.F. Progressive sequence alignment as a prerequisite to correct phylogenetic trees. J. Mol. Evol. 1987, 25, 351–360. [Google Scholar] [CrossRef] [PubMed]
  46. Pirovano, W.; Heringa, J. Multiple sequence alignment. Methods Mol. Biol. 2008, 452, 143–161. [Google Scholar]
  47. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Sim, N.L.; Kumar, P.; Hu, J.; Henikoff, S.; Schneider, G.; Ng, P.C. SIFT web server: Predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012, 40, 452–457. [Google Scholar] [CrossRef] [PubMed]
  49. Adzhubei, I.A.; Schmidt, S.; Peshkin, L.; Ramensky, V.E.; Gerasimova, A.; Bork, P.; Kondrashov, A.S.; Sunyaev, S.R. A method and server for predicting damaging missense mutations. Nat. Methods 2010, 7, 248–249. [Google Scholar] [CrossRef] [Green Version]
  50. Bendl, J.; Stourac, J.; Salanda, O.; Pavelka, A.; Wieben, E.D.; Zendulka, J.; Brezovsky, J.; Damborsky, J. PredictSNP: Robust and accurate consensus classifier for prediction of disease-related mutations. PLoS Comput. Biol. 2014, 10, e1003440. [Google Scholar] [CrossRef]
  51. Tang, H.; Thomas, P.D. PANTHER-PSEP: Predicting disease-causing genetic variants using position-specific evolutionary preservation. Bioinformatics 2016, 32, 2230–2232. [Google Scholar] [CrossRef] [Green Version]
  52. Capriotti, E.; Calabrese, R.; Fariselli, P.; Martelli, P.L.; Altman, R.B.; Casadio, R. WS-SNPs&GO: A web server for predicting the deleterious effect of human protein variants using functional annotation. BMC Genom. 2013, 14 (Suppl. S3), S6. [Google Scholar]
  53. Hecht, M.; Bromberg, Y.; Rost, B. Better prediction of functional effects for sequence variants. BMC Genom. 2015, 16 (Suppl. S8), S1. [Google Scholar] [CrossRef] [Green Version]
  54. Capriotti, E.; Calabrese, R.; Casadio, R. Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information. Bioinformatics 2006, 22, 2729–2734. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Capriotti, E.; Fariselli, P.; Calabrese, R.; Casadio, R. Predicting protein stability changes from sequences using support vector machines. Bioinformatics 2005, 21 (Suppl. S2), ii54–ii58. [Google Scholar] [CrossRef] [PubMed]
  56. Venselaar, H.; Te Beek, T.A.; Kuipers, R.K.; Hekkelman, M.L.; Vriend, G. Protein structure analysis of mutations causing inheritable diseases. An e-Science approach with life scientist friendly interfaces. BMC Bioinform. 2010, 11, 548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Blum, M.; Chang, H.; Chuguransky, S.; Grego, T.; Kandasaamy, S.; Mitchell, A.; Nuka, G.; Paysan-Lafosse, T.; Qureshi, M.; Raj, S.; et al. The InterPro protein families and domains database: 20 years on. Nucleic Acids Res. 2021, 49, D344–D354. [Google Scholar] [CrossRef]
  58. Ashkenazy, H.; Abadi, S.; Martz, E.; Chay, O.; Mayrose, I.; Pupko, T.; Ben-Tal, N. ConSurf 2016: An improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res. 2016, 44, W344–W350. [Google Scholar] [CrossRef] [Green Version]
  59. Hassan, M.M.; Hussain, M.A.; Kambal, S.; Elshikh, A.A.; Gendeel, O.R.; Ahmed, S.A.; Altayeb, R.A.; Muhajir, A.M.; Mohamed, S.B. NeoCoV Is Closer to MERS-CoV than SARS-CoV. Infect. Dis. 2020, 13, 1178633720930711. [Google Scholar] [CrossRef]
  60. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic. Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef]
  61. Himmelfarb, J.; Vanholder, R.; Mehrotra, R.; Tonelli, M. The current and future landscape of dialysis. Nat. Rev. Nephrol. 2020, 16, 573–585. [Google Scholar] [CrossRef]
  62. Khan, T.; Rahman, M.; Ahmed, I.; Al Ali, F.; Jithesh, P.V.; Marr, N. Human leukocyte antigen class II gene diversity tunes antibody repertoires to common pathogens. Front. Immunol. 2022, 13, 856497. [Google Scholar] [CrossRef]
  63. Kobayashi, T.; Yokoyama, I.; Uchida, K.; Orihara, A.; Takagi, H. HLA-DRB1 matching as a recipient selection criterion in cadaveric renal transplantation. Transplantation 1993, 55, 1294–1297. [Google Scholar] [CrossRef]
  64. Kamoun, M.; McCullough, K.P.; Maiers, M.; Fernandez Vina, M.A.; Li, H.; Teal, V.; Leichtman, A.B.; Merion, R.M. HLA Amino Acid Polymorphisms and Kidney Allograft Survival. Transplantation 2017, 101, e170–e177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Speer, S.L.; Zheng, W.; Jiang, X.; Chu, I.T.; Guseman, A.J.; Liu, M.; Pielak, G.J.; Li, C. The intracellular environment affects protein-protein interactions. Proc. Natl. Acad. Sci. USA 2021, 118, e2019918118. [Google Scholar] [CrossRef] [PubMed]
  66. Glass, F.; Takenaka, M. The Yeast Three-Hybrid System for Protein Interactions. Methods Mol. Biol. 2018, 1794, 195–205. [Google Scholar] [PubMed]
Figure 1. Flowchart of the work methodology. The names that appear after the symbol “-” represent prediction tools.
Figure 1. Flowchart of the work methodology. The names that appear after the symbol “-” represent prediction tools.
Medicina 59 01116 g001
Figure 2. Structural alteration by HOPE server. The protein is shown in grey, the wild-type residue in green, and the variant residue in red. The structures from left to right represent R109S, Y59H, and K41N variants.
Figure 2. Structural alteration by HOPE server. The protein is shown in grey, the wild-type residue in green, and the variant residue in red. The structures from left to right represent R109S, Y59H, and K41N variants.
Medicina 59 01116 g002
Figure 3. Evolutionary conservancy of HLA-DRB1 produced by Consurf server.
Figure 3. Evolutionary conservancy of HLA-DRB1 produced by Consurf server.
Medicina 59 01116 g003
Figure 4. Protein-protein interaction network of HLA-DRB1 predicted by the STRING database. The current illustration depicts the ten proteins with the highest total probability score. The red, green, blue, purple, yellow, light blue, and black lines indicate the presence of fusion, neighborhood, co-occurrence, experimental, text mining, database, and co-expression evidence, respectively.
Figure 4. Protein-protein interaction network of HLA-DRB1 predicted by the STRING database. The current illustration depicts the ten proteins with the highest total probability score. The red, green, blue, purple, yellow, light blue, and black lines indicate the presence of fusion, neighborhood, co-occurrence, experimental, text mining, database, and co-expression evidence, respectively.
Medicina 59 01116 g004
Table 1. General information of detected SNVs.
Table 1. General information of detected SNVs.
Serial No.SMPs Code No.N.P Exon 2Chr. LocationVariantsVariants TypeSNV AVAIL.AA ChangeVariants Allele Frequencies ALL/African ClinVar
1000 GenomeGnomAD GenomesNCBI ALFA
1.R3/C1256:32584113C/ASNVNovelR122R----
2.R20146:32,584,122T/CSNVrs1,136,782T119T0.065/0.1290.001/0.0020.006/0.013-
3.R16446:32,584,152T/ASNVrs750,986,830R109S----
4.R91966:32,584,304A/GSNVrs11,554,462Y59H0.072/0.1670.169/0.3180.166/0.290-
5.R3/C102426:32,584,350C/TSNVrs17,885,011E43E0.063/0.0600.112/0.0690.096/0.071-
6.A7/R122486:32,584,356C/TSNVrs17,887,028K41K--0.000/0.000-
7.R32486:32,584,356C/ASNVNovelK41N----
SMPs code No.: The code number of samples containing genetic variant(s), which are as follows: Rejection (R), acceptance (A), and control (C). N.P Exon 2: Position of a nucleotide in Exon 2. Chr. Location: Chromosomes Location. SNV AVAIL.: SNV availability in the Single Nucleotide Polymorphism database (dbSNP). AA Change: Type and location of amino acid changed. ClinVar: ClinVar database clinical significance records (https://www.ncbi.nlm.nih.gov/clinvar/) (accessed on 1 December 2022). The symbol “-” refers to unavailable data.
Table 2. The functional effect and stability index of non-synonymous detected SNVs.
Table 2. The functional effect and stability index of non-synonymous detected SNVs.
Variant DescriptionSNV IDAA
Change
SIFTPolyPhen-2PredictSNPPANTHERSNP&GOSNAP2PhD-SNPI-Mutant
chr6(GRCh38.p12):32584152T>Ars750986830R109SDeleterious
(0.0)
Probably damaging
(0.999)
DeleteriousPossibly damagingDisease (0.908)Effect (47)DiseaseDecrease (−1.29)
chr6(GRCh38.p12):32584304A>Grs11554462Y59HTolerated
(0.3)
Possibly damaging
(0.833)
NeutralProbably damagingDisease (0.551)Effect (52)NeutralDecrease (−1.60)
chr6(GRCh38.p12):32584356C>ANovelK41NDeleterious
(0.01)
Probably damaging
(0.993)
DeleteriousPossibly damagingDisease (0.785)Effect (44)DiseaseDecrease (−0.34)
Variant Description: In variant call format (VCF). SNV ID: Accession number in dbSNP. AA Change: Type and location of amino acid changed. The value of the predicted score is represented by the numbers between the brackets.
Table 3. Structural effects of nsSNVs in a protein sequence using the HOPE server.
Table 3. Structural effects of nsSNVs in a protein sequence using the HOPE server.
Distinctions and Characteristicsrs750986830
R109S
rs11554462
Y59H
Novel
K41N
ArginineSerineTyrosineHistidineLysineAsparagine
Schematic structuresMedicina 59 01116 i001Medicina 59 01116 i002Medicina 59 01116 i003Medicina 59 01116 i004Medicina 59 01116 i005Medicina 59 01116 i006
SizelargeSmalllargeSmalllargeSmall
ChargePositiveNeutral--PositiveNeutral
Hydrophobicity-valueLess hydrophobicMore hydrophobicMore hydrophobicLess hydrophobic--
ContactsThe wild-type residue forms eight hydrogen bonds and one salt bridge with other residues.The mutant-type has an impact on the original’s hydrogen bond formation, binding site, and ionic interactions.The wild-type residue forms a hydrogen bond with eight residues.The mutant-type has an impact on the original’s hydrogen bond formation.The wild-type forms a hydrogen bond, a salt bridge with one residue, and is involved in multimer contacts.The mutant type affects hydrogen bond formation, ionic interaction, and the development of multimer interactions.
StructureThe mutation is located within a stretch of residues annotated in UniProt as a special region: Beta-1. The differences in amino acid properties can disturb this region and disturb its function.
UniProt: Universal database of protein (https://www.uniprot.org/) (accessed on 10 December 2022). The symbol “-” refers to unavailable data.
Table 4. The effect of nsSNVs on HLA-DRB1’ protein physicochemical parameters.
Table 4. The effect of nsSNVs on HLA-DRB1’ protein physicochemical parameters.
Reference & VariantsMolecular WeightTheoretical pIAtomic CompositionTotal −veTotal +veExtinction CoefficientsInstability IndexAliphatic IndexGRAVY
Reference29966.147.64C1342H2068N368O389S1225264128548.9277.93−0.207
R109S29897.037.00C1339H2061N365O390S1225254128548.2077.93−0.193
Y59H29940.107.66C1339H2066N370O388S1225263979549.5477.93−0.214
K41N29952.077.00C1340H2062N368O390S1225254128547.7077.93−0.205
The accession number for the reference sequence is P01911 (https://www.uniprot.org/) (accessed on 15 December 2022). Total –ve: Total negatively charged residues. Total +ve: Total positively charged residues. GRAVY: Grand average of hydropathicity index. The parameters that have been changed compared to the reference are highlighted in bold.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hassan, M.M.; Hussain, M.A.; Ali, S.S.; Mahdi, M.A.; Mohamed, N.S.; AbdElbagi, H.; Mohamed, O.; Sherif, A.E.; Osman, W.; Ibrahim, S.R.M.; et al. Detection of Nonsynonymous Single Variants in Human HLA-DRB1 Exon 2 Associated with Renal Transplant Rejection. Medicina 2023, 59, 1116. https://doi.org/10.3390/medicina59061116

AMA Style

Hassan MM, Hussain MA, Ali SS, Mahdi MA, Mohamed NS, AbdElbagi H, Mohamed O, Sherif AE, Osman W, Ibrahim SRM, et al. Detection of Nonsynonymous Single Variants in Human HLA-DRB1 Exon 2 Associated with Renal Transplant Rejection. Medicina. 2023; 59(6):1116. https://doi.org/10.3390/medicina59061116

Chicago/Turabian Style

Hassan, Mohamed M., Mohamed A. Hussain, Sababil S. Ali, Mohammed A. Mahdi, Nouh Saad Mohamed, Hanadi AbdElbagi, Osama Mohamed, Asmaa E. Sherif, Wadah Osman, Sabrin R. M. Ibrahim, and et al. 2023. "Detection of Nonsynonymous Single Variants in Human HLA-DRB1 Exon 2 Associated with Renal Transplant Rejection" Medicina 59, no. 6: 1116. https://doi.org/10.3390/medicina59061116

APA Style

Hassan, M. M., Hussain, M. A., Ali, S. S., Mahdi, M. A., Mohamed, N. S., AbdElbagi, H., Mohamed, O., Sherif, A. E., Osman, W., Ibrahim, S. R. M., Ghazawi, K. F., Miski, S. F., Mohamed, G. A., & Ashour, A. (2023). Detection of Nonsynonymous Single Variants in Human HLA-DRB1 Exon 2 Associated with Renal Transplant Rejection. Medicina, 59(6), 1116. https://doi.org/10.3390/medicina59061116

Article Metrics

Back to TopTop