Haplotype-based analysis resolves missing heritability in oculocutaneous albinism type 1B.

Oculocutaneous albinism (OCA) is a rare disorder of pigment production. Affected individuals have variably decreased global pigmentation and visual-developmental changes that lead to low vision. OCA is notable for significant missing heritability, particularly among individuals with residual pigmentation. Tyrosinase (TYR) is the rate-limiting enzyme in melanin pigment biosynthesis and mutations that decrease enzyme function are one of the most common causes of OCA. We present the analysis of high-depth short-read TYR sequencing data for a cohort of 352 OCA probands, ∼50% of whom were previously sequenced without yielding a definitive diagnostic result. Our analysis identified 66 TYR single-nucleotide variants (SNVs) and small insertion/deletions (indels), 3 structural variants, and a rare haplotype comprised of two common frequency variants (p.Ser192Tyr and p.Arg402Gln) in cis-orientation, present in 149/352 OCA probands. We further describe a detailed analysis of the disease-causing haplotype, p.[Ser192Tyr; Arg402Gln] ("cis-YQ"). Haplotype analysis suggests that the cis-YQ allele arose by recombination and that multiple cis-YQ haplotypes are segregating in OCA-affected individuals and control populations. The cis-YQ allele is the most common disease-causing allele in our cohort, representing 19.1% (57/298) of TYR pathogenic alleles in individuals with type 1 (TYR-associated) OCA. Finally, among the 66 TYR variants, we found several additional alleles defined by a cis-oriented combination of minor, potentially hypomorph-producing alleles at common variant sites plus a second, rare pathogenic variant. Together, these results suggest that identification of phased variants for the full TYR locus are required for an exhaustive assessment for potentially disease-causing alleles.

Extreme phenotype sampling and next generation sequencing to identify genetic variants associated with tacrolimus in African American kidney transplant recipients.

African American (AA) kidney transplant recipients (KTRs) have poor outcomes, which may in-part be due to tacrolimus (TAC) sub-optimal immunosuppression. We previously determined the common genetic regulators of TAC pharmacokinetics in AAs which were CYP3A5 *3, *6, and *7. To identify low-frequency variants that impact TAC pharmacokinetics, we used extreme phenotype sampling and compared individuals with extreme high (n = 58) and low (n = 60) TAC troughs (N = 515 AA KTRs). Targeted next generation sequencing was conducted in these two groups. Median TAC troughs in the high group were 7.7 ng/ml compared with 6.3 ng/ml in the low group, despite lower daily doses of 5 versus 12 mg, respectively. Of 34,542 identified variants across 99 genes, 1406 variants were suggestively associated with TAC troughs in univariate models (p-value < 0.05), however none were significant after multiple testing correction. We suggest future studies investigate additional sources of TAC pharmacokinetic variability such as drug-drug-gene interactions and pharmacomicrobiome.

The Human Genome Organisation (HUGO) and a vision for Ecogenomics: the Ecological Genome Project.

BACKGROUND: The following outlines ethical reasons for widening the Human Genome Organisation's (HUGO) mandate to include ecological genomics. MAIN: The environment influences an organism's genome through ambient factors in the biosphere (e.g. climate and UV radiation), as well as the agents it comes into contact with, i.e. the epigenetic and mutagenic effects of inanimate chemicals and pollution, and pathogenic organisms. Emerging scientific consensus is that social determinants of health, environmental conditions and genetic factors work together to influence the risk of many complex illnesses. That paradigm can also explain the environmental and ecological determinants of health as factors that underlie the (un)healthy ecosystems on which communities rely. We suggest that The Ecological Genome Project is an aspirational opportunity to explore connections between the human genome and nature. We propose consolidating a view of Ecogenomics to provide a blueprint to respond to the environmental challenges that societies face. This can only be achieved by interdisciplinary engagement between genomics and the broad field of ecology and related practice of conservation. In this respect, the One Health approach is a model for environmental orientated work. The idea of Ecogenomics-a term that has been used to relate to a scientific field of ecological genomics-becomes the conceptual study of genomes within the social and natural environment. CONCLUSION: The HUGO Committee on Ethics, Law and Society (CELS) recommends that an interdisciplinary One Health approach should be adopted in genomic sciences to promote ethical environmentalism. This perspective has been reviewed and endorsed by the HUGO CELS and the HUGO Executive Board.

Steroid-tacrolimus drug-drug interaction and the effect of CYP3A genotypes.

AIMS: Tacrolimus, metabolized by CYP3A4 and CYP3A5 enzymes, is susceptible to drug-drug interactions (DDI). Steroids induce CYP3A genes to increase tacrolimus clearance, but the effect is variable. We hypothesized that the extent of the steroid-tacrolimus DDI differs by CYP3A4/5 genotypes. METHODS: Kidney transplant recipients (n = 2462) were classified by the number of loss of function alleles (LOF) (CYP3A5*3, *6 and *7 and CYP3A4*22) and steroid use at each tacrolimus trough in the first 6 months post-transplant. A population pharmacokinetic analysis was performed by nonlinear mixed-effect modelling (NONMEM) and stepwise covariate modelling to define significant covariates affecting tacrolimus clearance. A stochastic simulation was performed and translated into a Shiny application with the mrgsolve and Shiny packages in R. RESULTS: Steroids were associated with modestly higher (3%-11.8%) tacrolimus clearance. Patients with 0-LOF alleles receiving steroids showed the greatest increase (11.8%) in clearance compared to no steroids, whereas those with 2-LOFs had a negligible increase (2.6%) in the presence of steroids. Steroid use increased tacrolimus clearance by 5% and 10.3% in patients with 1-LOF and 3/4-LOFs, respectively. CONCLUSIONS: Steroids increase the clearance of tacrolimus but vary slightly by CYP3A genotype. This is important in individuals of African ancestry who are more likely to carry no LOF alleles, may more commonly receive steroid treatment, and will need higher tacrolimus doses.

Limited Sampling Strategies Fail to Accurately Predict Mycophenolic Acid Area Under the Curve in Kidney Transplant Recipients and the Impact of Enterohepatic Recirculation.

BACKGROUND: Therapeutic drug monitoring for mycophenolic acid (MPA) is challenging due to difficulties in measuring the area under the curve (AUC). Limited sampling strategies (LSSs) have been developed for MPA therapeutic drug monitoring but come with risk of unacceptable performance. The authors hypothesized that the poor predictive performance of LSSs were due to the variability in MPA enterohepatic recirculation (EHR). This study is the first to evaluate LSSs models performance in the context of EHR. METHODS: Adult kidney transplant recipients (n = 84) receiving oral mycophenolate mofetil underwent intensive MPA pharmacokinetic sampling. MPA AUC0-12hr and EHR were determined. Published MPA LSSs in kidney transplant recipients receiving tacrolimus were evaluated for their predictive performance in estimating AUC0-12hr in our full cohort and separately in individuals with high and low EHR. RESULTS: None of the evaluated LSS models (n = 12) showed good precision or accuracy in predicting MPA AUC0-12hr in the full cohort. In the high EHR group, models with late timepoints had better accuracy but low precision, except for 1 model with late timepoints at 6 and 10 hours postdose, which had marginally acceptable precision. For all models, the good guess of predicted AUC0-12hr (±15% of observed AUC0-12hr) was highly variable (range, full cohort = 19%-61.9%; high EHR = 4.5%-65.9%; low EHR = 27.5%-62.5%). CONCLUSIONS: The predictive performance of the LSS models varied according to EHR status. Timepoints ≥5 hours postdose in LSS models are essential to capture EHR. Models and strategies that incorporate EHR during development are required to accurately ascertain MPA exposure.

Higher number of tacrolimus dose adjustments in kidney transplant recipients who are extensive and intermediate CYP3A5 metabolizers.

Kidney transplant recipients carrying the CYP3A5*1 allele have lower tacrolimus troughs, and higher dose requirements compared to those with the CYP3A5*3/*3 genotype. However, data on the effect of CYP3A5 alleles on post-transplant tacrolimus management are lacking. The effect of CYP3A5 metabolism phenotypes on the number of tacrolimus dose adjustments and troughs in the first 6 months post-transplant was evaluated in 78 recipients (64% Caucasians). Time to first therapeutic concentration, percentage of time in therapeutic range (TTR), and estimated glomerular filtration rate (eGFR) were also evaluated. Fifty-five kidney transplant recipients were CYP3A5 poor metabolizers (PM), 17 were intermediate metabolizers (IM), and 6 were extensive metabolizers (EM). Compared to PMs, EMs/IMs had significantly more dose adjustments (6.1 vs. 8.1, p = .015). Overall, 33.82% of trough measurements resulted in a dose change. There was no difference in the number of tacrolimus trough measurements between PMs and EM/IMs. The total daily tacrolimus dose requirements were higher in EMs and IMs compared to PMs (<.001). TTR was ∼50% in the PMs and EMs/IMs groups. CYP3A5 EM/IM metabolizers have more tacrolimus dose changes and higher dose requirements which increases clinical management complexity. Larger studies are needed to assess the cost and benefits of including genotyping data to improve clinical management.

A custom capture sequence approach for oculocutaneous albinism identifies structural variant alleles at the OCA2 locus.

Oculocutaneous albinism (OCA) is a heritable disorder of pigment production that manifests as hypopigmentation and altered eye development. Exon sequencing of known OCA genes is unsuccessful in producing a complete molecular diagnosis for a significant number of affected individuals. We sequenced the DNA of individuals with OCA using short-read custom capture sequencing that targeted coding, intronic, and noncoding regulatory regions of known OCA genes, and genome-wide association study-associated pigmentation loci. We identified an OCA2 complex structural variant (CxSV), defined by a 143 kb inverted segment reintroduced in intron 1, upstream of the native location. The corresponding CxSV junctions were observed in 11/390 probands screened. The 143 kb CxSV presents in one family as a copy number variant duplication for the 143 kb region. In the remaining 10/11 families, the 143 kb CxSV acquired an additional 184 kb deletion across the same region, restoring exons 3-19 of OCA2 to a copy-number neutral state. Allele-associated haplotype analysis found rare SNVs rs374519281 and rs139696407 are linked with the 143 kb CxSV in both OCA2 alleles. For individuals in which customary molecular evaluation does not reveal a biallelic OCA diagnosis, we recommend preliminary screening for these haplotype-associated rare variants, followed by junction-specific validation for the OCA2 143 kb CxSV.

Pharmacogenomics in kidney transplant recipients and potential for integration into practice.

WHAT IS KNOWN AND OBJECTIVE: Pharmacogenomic biomarkers are now used in many clinical care settings and represent one of the successes of precision medicine. Genetic variants are associated with pharmacokinetic and pharmacodynamic changes leading to medication adverse effects and changes in clinical response. Actionable pharmacogenomic variants are common in transplant recipients and have implications for medications used in transplant, but yet are not broadly incorporated into practice. METHODS: From the Clinical Pharmacogenetics Implementation Consortium and Dutch Pharmacogenetics Working Group guidelines, and PharmGKB databases, 12 pharmacogenomic genes with 30 variants were selected and used to create diplotypes and actionable pharmacogenomic phenotypes. A total of 853 kidney allograft recipients who had genomic information available from a genome-wide association study were included. RESULTS: Each recipient had at least one actionable pharmacogenomic diplotype/phenotype, whereas the majority (58%) had three or four actionable diplotypes/phenotypes and 17.4% had five or more among the 12 genes. The participants carried actionable diplotypes/phenotypes for multiple medications, including tacrolimus, azathioprine, clopidogrel, warfarin, simvastatin, voriconazole, antidepressants and proton-pump inhibitors. WHAT IS NEW AND CONCLUSION: Pharmacogenomic variants are common in transplant recipients, and transplant recipients receive medications that have actionable variants. CLINICAL TRIAL: Genomics of Transplantation, clinicaltrials.gov (NCT01714440).

Tacrolimus troughs and genetic determinants of metabolism in kidney transplant recipients: A comparison of four ancestry groups.

Tacrolimus trough and dose requirements vary dramatically between individuals of European and African American ancestry. These differences are less well described in other populations. We conducted an observational, prospective, multicenter study from which 2595 kidney transplant recipients of European, African, Native American, and Asian ancestry were studied for tacrolimus trough, doses, and genetic determinants of metabolism. We studied the well-known variants and conducted a CYP3A4/5 gene-wide analysis to identify new variants. Daily doses, and dose-normalized troughs were significantly different between the four groups (P < .001). CYP3A5*3 (rs776746) was associated with higher dose-normalized tacrolimus troughs in all groups but occurred at different allele frequencies and had differing effect sizes. The CYP3A5*6 (rs10264272) and *7 (rs413003343) variants were only present in African Americans. CYP3A4*22 (rs35599367) was not found in any of the Asian ancestry samples. We identified seven suggestive variants in the CYP3A4/5 genes associated with dose-normalized troughs in Native Americans (P = 1.1 × 10-5 -8.8 × 10-6 ) and one suggestive variant in Asian Americans (P = 5.6 × 10-6 ). Tacrolimus daily doses and dose-normalized troughs vary significantly among different ancestry groups. We identified potential new variants important in Asians and Native Americans. Studies with larger populations should be conducted to assess the importance of the identified suggestive variants.

The impact of donor and recipient common clinical and genetic variation on estimated glomerular filtration rate in a European renal transplant population.

Genetic variation across the human leukocyte antigen loci is known to influence renal-transplant outcome. However, the impact of genetic variation beyond the human leukocyte antigen loci is less clear. We tested the association of common genetic variation and clinical characteristics, from both the donor and recipient, with posttransplant eGFR at different time-points, out to 5 years posttransplantation. We conducted GWAS meta-analyses across 10 844 donors and recipients from five European ancestry cohorts. We also analyzed the impact of polygenic risk scores (PRS), calculated using genetic variants associated with nontransplant eGFR, on posttransplant eGFR. PRS calculated using the recipient genotype alone, as well as combined donor and recipient genotypes were significantly associated with eGFR at 1-year posttransplant. Thirty-two percent of the variability in eGFR at 1-year posttransplant was explained by our model containing clinical covariates (including weights for death/graft-failure), principal components and combined donor-recipient PRS, with 0.3% contributed by the PRS. No individual genetic variant was significantly associated with eGFR posttransplant in the GWAS. This is the first study to examine PRS, composed of variants that impact kidney function in the general population, in a posttransplant context. Despite PRS being a significant predictor of eGFR posttransplant, the effect size of common genetic factors is limited compared to clinical variables.

Identification of genetic variants associated with tacrolimus metabolism in kidney transplant recipients by extreme phenotype sampling and next generation sequencing.

An extreme phenotype sampling (EPS) model with targeted next-generation sequencing (NGS) identified genetic variants associated with tacrolimus (Tac) metabolism in subjects from the Deterioration of Kidney Allograft Function (DeKAF) Genomics cohort which included 1,442 European Americans (EA) and 345 African Americans (AA). This study included 48 subjects separated into 4 groups of 12 (AA high, AA low, EA high, EA low). Groups were selected by the extreme phenotype of dose-normalized Tac trough concentrations after adjusting for common genetic variants and clinical factors. NGS spanned > 3 Mb of 28 genes and identified 18,661 genetic variants (3961 previously unknown). A group of 125 deleterious variants, by SIFT analysis, were associated with Tac troughs in EAs (burden test, p = 0.008), CYB5R2 was associated with Tac troughs in AAs (SKAT, p = 0.00079). In CYB5R2, rs61733057 (increased allele frequency in AAs) was predicted to disrupt protein function by SIFT and PolyPhen2 analysis. The variants merit further validation.

NPHP1 (Nephrocystin-1) Gene Deletions Cause Adult-Onset ESRD.

Background Nephronophthisis (NPH) is the most prevalent genetic cause for ESRD in children. However, little is known about the prevalence of NPH in adult-onset ESRD. Homozygous full gene deletions of the NPHP1 gene encoding nephrocystin-1 are a prominent cause of NPH. We determined the prevalence of NPH in adults by assessing homozygous NPHP1 full gene deletions in adult-onset ESRD.Methods Adult renal transplant recipients from five cohorts of the International Genetics and Translational Research in Transplantation Network (iGeneTRAiN) underwent single-nucleotide polymorphism genotyping. After quality control, we determined autosomal copy number variants (such as deletions) on the basis of median log2 ratios and B-allele frequency patterns. The findings were independently validated in one cohort. Patients were included in the analysis if they had adult-onset ESRD, defined as start of RRT at ≥18 years old.Results We included 5606 patients with adult-onset ESRD; 26 (0.5%) showed homozygous NPHP1 deletions. No donor controls showed homozygosity for this deletion. Median age at ESRD onset was 30 (range, 18-61) years old for patients with NPH, with 54% of patients age ≥30 years old. Notably, only three (12%) patients were phenotypically classified as having NPH, whereas most patients were defined as having CKD with unknown etiology (n=11; 42%).Conclusions Considering that other mutation types in NPHP1 or mutations in other NPH-causing genes were not analyzed, NPH is a relatively frequent monogenic cause of adult-onset ESRD. Because 88% of patients had not been clinically diagnosed with NPH, wider application of genetic testing in adult-onset ESRD may be warranted.

Urinary microbiome associated with chronic allograft dysfunction in kidney transplant recipients.

BACKGROUND: We performed a study to identify differences in the urinary microbiome associated with chronic allograft dysfunction (CAD) and compared the urinary microbiome of male and female transplant recipients with CAD. METHODS: This case-control study enrolled 67 patients within the Deterioration of Kidney Allograft Function (DeKAF) Genomics cohort at two transplant centers. CAD was defined as a greater than 25% rise in serum creatinine relative to a 3 month post-transplant baseline. Urine samples from patients with and without CAD were analyzed using 16S V4 bacterial ribosomal DNA sequences. RESULTS: Corynebacterium was more prevalent in female and male patients with CAD compared to non-CAD female patients (P = 0.0005). A total 21 distinct Operational Taxonomic Unit (OTUs) were identified as significantly different when comparing CAD and non-CAD patients using Kruskal-Wallis (P < 0.01). A subset analysis of female patients with CAD compared to non-CAD females identified similar differentially abundant OTUs, including the genera Corynebacterium and Staphylococcus (Kruskal-Wallis; P = 0.01; P = 0.004, respectively). Male CAD vs female CAD analysis showed greater abundance of phylum Proteobacteria in males. CONCLUSION: There were differences in the urinary microbiome when comparing female and male CAD patients with their female non-CAD counterparts and these differences persisted in the subset analysis limited to female patients only.

Tacrolimus trough and dose intra-patient variability and CYP3A5 genotype: Effects on acute rejection and graft failure in European American and African American kidney transplant recipients.

BACKGROUND: Suboptimal immunosuppression after kidney transplantation contributes to toxicity and loss of efficacy. Little is known regarding the impact of intra-patient variability of tacrolimus (TAC) doses and troughs in the early post-transplant period or the influence of genetic variants on variability. METHODS: Coefficients of variation (CV) of TAC troughs and doses of 1226 European American (EA) and 246 African American (AA) adult recipients enrolled in DeKAF Genomics were compared for association with acute rejection and graft failure. Additionally, the influence of recipients' number of CYP3A5 loss-of-function alleles was assessed. RESULTS: Acute rejection was associated with greater CV of dose in AA (P < 0.001) and EA recipients (P = 0.012). Graft failure was associated with a greater CV of dose (P = 0.022) and trough (P < 0.001) in AA, and higher CV of trough (P = 0.024) in EA recipients. In EA, CYP3A5 loss-of-function alleles were associated with decreased CV of trough (P = 0.0042) and increased CV of dose (P < 0.0001). CONCLUSION: CYP3A5 loss-of-function alleles influence intra-patient TAC trough and dose variability. High variability of TAC dose increases risk of acute rejection. High variability of TAC trough increases risk of graft failure. Early clinical recognition of TAC dose and trough variability may improve patient management and outcomes.

Genetics of acute rejection after kidney transplantation.

Treatment of acute rejection (AR) following kidney transplantation has improved in recent years, but there are still limitations to successful outcomes. This review article covers literature in regard to recipient and donor genetics of AR kidney and secondarily of liver allografts. Many candidate gene and some genome-wide association studies (GWASs) have been conducted for AR in kidney transplantation. Genetic associations with AR in kidney and liver are mostly weak, and in most cases, the associations have not been reproducible. A limitation in the study of AR is the lack of sufficiently large populations that account for population stratification to study the AR phenotype which in this era occurs in <10% of transplants. Furthermore, the AR phenotype has been difficult to define and the definitions of classifications have evolved over time. Literature related to the pharmacogenomics of tacrolimus is robust and has been validated in many studies. Associations between gene expression and AR are emerging as markers of outcomes and AR classification. In the future, combinations of pretransplant genotype for AR risk prediction, genotype-based immune suppressant dosing, and pharmacogenomic markers to select AR maintenance or treatment and expression markers from biopsies may provide valuable clinical tools for guiding treatment.

CRISPR/Cas9 Genetic Modification of CYP3A5 *3 in HuH-7 Human Hepatocyte Cell Line Leads to Cell Lines with Increased Midazolam and Tacrolimus Metabolism.

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 engineering of the CYP3A5 *3 locus (rs776746) in human liver cell line HuH-7 (CYP3A5 *3/*3) has led to three CYP3A5 *1 cell lines by deletion of the exon 3B splice junction or point mutation. Cell lines CYP3A5 *1/*3 sd (single deletion), CYP3A5 *1/*1 dd (double deletion), or CYP3A5 *1/*3 pm (point mutation) expressed the CYP3A5 *1 mRNA and had elevated CYP3A5 mRNA (P < 0.0005 for all engineered cell lines) and protein expression compared with HuH-7. In metabolism assays, HuH-7 had less tacrolimus (all P < 0.05) or midazolam (MDZ) (all P < 0.005) disappearance than all engineered cell lines. HuH-7 had less 1-OH MDZ (all P < 0.0005) or 4-OH (all P < 0.005) production in metabolism assays than all bioengineered cell lines. We confirmed CYP3A5 metabolic activity with the CYP3A4 selective inhibitor CYP3CIDE. This is the first report of genomic CYP3A5 bioengineering in human cell lines with drug metabolism analysis.

Clinical Interpretation of Variants from Next-Generation Sequencing: The 2016 Scientific Meeting of the Human Genome Variation Society.

The 2016 scientific meeting of the Human Genome Variation Society (HGVS; http://www.hgvs.org) was held on the 20th of May in Barcelona, Spain, with the theme of "Clinical Interpretation of Variants from Next-Generation Sequencing."

Clinically relevant variants - identifying, collecting, interpreting, and disseminating: the 2013 annual scientific meeting of the Human Genome Variation Society.

The dramatic advances in genetic sequencing technologies used in research laboratories are now entering the clinic, and applications of whole-genome and whole-exome sequencing to disease diagnosis, predisposition, and treatment will soon be commonplace. However, the standards and methods for identifying clinically relevant variants are currently being debated and defined. Multiple agencies worldwide have recognized that we have reached an exciting and critical transition point into the clinic, and many important issues are being discussed that impact how genetic variation data in the clinic will be interpreted and used. The 2013 annual scientific meeting of the Human Genome Variation Society (HGVS) had as its main theme the discovery, interpretation, and dissemination of clinically relevant DNA variants. The meeting featured the continuously developing technology of databasing genetic variation and computational tools for allelic variant discovery. Attention was given to curating and integrating these data with clinical findings, including approaches to distinguish between functional alleles underlying clinical phenotypes and benign sequence variants and making data sources interoperable and functional for clinical diagnostic utility, citing examples in specific diseases.

Gamma-aminobutyric acid system genes--no evidence for a role in alcohol use and abuse in a community-based sample.

BACKGROUND: While twin and adoption studies point to substantial genetic influence upon alcohol use, dependence, and other alcohol-related phenotypes, few of the genes underlying variation in these phenotypes have been identified. Markers in genes related to GABAergic activity-a system integral to many of alcohol's biological effects-have been implicated in alcohol use and alcohol-related psychopathology in linkage and association studies. METHODS: Using multiple methods, we conducted a comprehensive examination of the effects of markers in γ-aminobutyric acid (GABA) system genes in a community-based sample of 7,224 individuals assessed in early and middle adulthood. In addition to testing the effect of individual single nucleotide polymorphism (SNP) markers on alcohol-related phenotypes, we computed a polygenic score reflecting the aggregated effects of multiple GABA system SNPs. We also estimated the variance in alcohol-related phenotypes attributable to all GABA system markers considered simultaneously and conducted gene-based association tests. RESULTS: No method produced results indicative of an effect of GABA system variants on measures of alcohol use or misuse. CONCLUSIONS: These results reflect alcohol-related behaviors in a population-representative sample, many of whom are still in adolescence, and in which the incidence of heavy drinking and alcohol-related symptomatology are relatively low. Contrasted with existing studies of the association between alcohol use and GABA system genes, our results suggest that the relationship may be limited to particular contexts, such as when accompanied by polysubstance abuse or a familial history of alcoholism.

Polygenic risk score for acute rejection based on donor-recipient non-HLA genotype mismatch.

BACKGROUND: Acute rejection (AR) after kidney transplantation is an important allograft complication. To reduce the risk of post-transplant AR, determination of kidney transplant donor-recipient mismatching focuses on blood type and human leukocyte antigens (HLA), while it remains unclear whether non-HLA genetic mismatching is related to post-transplant complications. METHODS: We carried out a genome-wide scan (HLA and non-HLA regions) on AR with a large kidney transplant cohort of 784 living donor-recipient pairs of European ancestry. An AR polygenic risk score (PRS) was constructed with the non-HLA single nucleotide polymorphisms (SNPs) filtered by independence (r2 < 0.2) and P-value (< 1×10-3) criteria. The PRS was validated in an independent cohort of 352 living donor-recipient pairs. RESULTS: By the genome-wide scan, we identified one significant SNP rs6749137 with HR = 2.49 and P-value = 2.15×10-8. 1,307 non-HLA PRS SNPs passed the clumping plus thresholding and the PRS exhibited significant association with the AR in the validation cohort (HR = 1.54, 95% CI = (1.07, 2.22), p = 0.019). Further pathway analysis attributed the PRS genes into 13 categories, and the over-representation test identified 42 significant biological processes, the most significant of which is the cell morphogenesis (GO:0000902), with 4.08 fold of the percentage from homo species reference and FDR-adjusted P-value = 8.6×10-4. CONCLUSIONS: Our results show the importance of donor-recipient mismatching in non-HLA regions. Additional work will be needed to understand the role of SNPs included in the PRS and to further improve donor-recipient genetic matching algorithms. Trial registry: Deterioration of Kidney Allograft Function Genomics (NCT00270712) and Genomics of Kidney Transplantation (NCT01714440) are registered on ClinicalTrials.gov.