Faculty Bibliography

In a previous study, heart xenografts from 10-gene-edited pigs transplanted into two human decedents did not show evidence of acute-onset cellular- or antibody-mediated rejection. Here, to better understand the detailed molecular landscape following xenotransplantation, we carried out bulk and single-cell transcriptomics, lipidomics, proteomics and metabolomics on blood samples obtained from the transplanted decedents every 6 h, as well as histological and transcriptomic tissue profiling. We observed substantial early immune responses in peripheral blood mononuclear cells and xenograft tissue obtained from decedent 1 (male), associated with downstream T cell and natural killer cell activity. Longitudinal analyses indicated the presence of ischemia reperfusion injury, exacerbated by inadequate immunosuppression of T cells, consistent with previous findings of perioperative cardiac xenograft dysfunction in pig-to-nonhuman primate studies. Moreover, at 42 h after transplantation, substantial alterations in cellular metabolism and liver-damage pathways occurred, correlating with profound organ-wide physiological dysfunction. By contrast, relatively minor changes in RNA, protein, lipid and metabolism profiles were observed in decedent 2 (female) as compared to decedent 1. Overall, these multi-omics analyses delineate distinct responses to cardiac xenotransplantation in the two human decedents and reveal new insights into early molecular and immune responses after xenotransplantation. These findings may aid in the development of targeted therapeutic approaches to limit ischemia reperfusion injury-related phenotypes and improve outcomes.

Optimizing immunologic compatibility in organ transplantation extends beyond the conventional approach of Human Leukocyte Antigen (HLA) antigen matching, which exhibits significant limitations. A broader comprehension of the roles of classical and non-classical HLA genes in transplantation is imperative for enhancing long-term graft survival. High-resolution molecular HLA genotyping, despite its inherent challenges, has emerged as the cornerstone for precise patient-donor compatibility assessment. Leveraging understanding of eplet biology and indirect immune activation, eplet mismatch calculators and the PIRCHE-II algorithm surpass traditional methods in predicting allograft rejection. Understanding minor histocompatibility antigens may also present an opportunity to personalize the compatibility process. While the application of molecular matching in deceased donor organ allocation presents multiple technical, logistical, and conceptual barriers, rendering it premature for mainstream use, several other areas of donor-recipient matching and post-transplant management are ready to incorporate molecular matching. Provision of molecular mismatch scores to physicians during potential organ offer evaluations could potentially amplify long-term outcomes. The implementation of molecular matching in living organ donation and kidney paired exchange programs is similarly viable. This article will explore the current understanding of immunologic matching in transplantation and the potential applications of epitope and non-epitope molecular biology and genetics in clinical transplantation.

Genome-wide association studies (GWAS) have yielded significant insights into the genetic architecture of myocardial infarction (MI), although studies in non-European populations are still lacking. Saudi Arabian cohorts offer an opportunity to discover novel genetic variants impacting disease risk due to a high rate of consanguinity. Genome-wide genotyping (GWG), imputation and GWAS followed by meta-analysis were performed based on two independent Saudi Arabian studies comprising 3950 MI patients and 2324 non-MI controls. Meta-analyses were then performed with these two Saudi MI studies and the CardioGRAMplusC4D and UK BioBank GWAS as controls. Meta-analyses of the two Saudi MI studies resulted in 17 SNPs with genome-wide significance. Meta-analyses of all 4 studies revealed 66 loci with genome-wide significance levels of p < 5 × 10^-8. All of these variants, except rs2764203, have previously been reported as MI-associated loci or to have high linkage disequilibrium with known loci. One SNP association in Shisa family member 5 (SHISA5) (rs11707229) was evident at a much higher frequency in the Saudi MI populations (> 12% MAF). In conclusion, our results replicated many MI associations, whereas in Saudi-only GWAS (meta-analyses), several new loci were implicated that require future validation and functional analyses.

BACKGROUND:Cross-species immunological incompatibilities have hampered pig-to-human xenotransplantation, but porcine genome engineering recently enabled the first successful experiments. However, little is known about the immune response after the transplantation of pig kidneys to human recipients. We aimed to precisely characterise the early immune responses to the xenotransplantation using a multimodal deep phenotyping approach. METHODS:We did a complete phenotyping of two pig kidney xenografts transplanted to decedent humans. We used a multimodal strategy combining morphological evaluation, immunophenotyping (IgM, IgG, C4d, CD68, CD15, NKp46, CD3, CD20, and von Willebrand factor), gene expression profiling, and whole-transcriptome digital spatial profiling and cell deconvolution. Xenografts before implantation, wild-type pig kidney autografts, as well as wild-type, non-transplanted pig kidneys with and without ischaemia-reperfusion were used as controls. FINDINGS:cells. Both xenografts showed increased expression of genes biologically related to a humoral response, including monocyte and macrophage activation, natural killer cell burden, endothelial activation, complement activation, and T-cell development. Whole-transcriptome digital spatial profiling showed that antibody-mediated injury was mainly located in the glomeruli of the xenografts, with significant enrichment of transcripts associated with monocytes, macrophages, neutrophils, and natural killer cells. This phenotype was not observed in control pig kidney autografts or in ischaemia-reperfusion models. INTERPRETATION:Despite favourable short-term outcomes and absence of hyperacute injuries, our findings suggest that antibody-mediated rejection in pig-to-human kidney xenografts might be occurring. Our results suggest specific therapeutic targets towards the humoral arm of rejection to improve xenotransplantation results. FUNDING:OrganX and MSD Avenir.

Genetically modified xenografts are one of the most promising solutions to the discrepancy between the numbers of available human organs for transplantation and potential recipients. To date, a porcine heart has been implanted into only one human recipient. Here, using 10-gene-edited pigs, we transplanted porcine hearts into two brain-dead human recipients and monitored xenograft function, hemodynamics and systemic responses over the course of 66 hours. Although both xenografts demonstrated excellent cardiac function immediately after transplantation and continued to function for the duration of the study, cardiac function declined postoperatively in one case, attributed to a size mismatch between the donor pig and the recipient. For both hearts, we confirmed transgene expression and found no evidence of cellular or antibody-mediated rejection, as assessed using histology, flow cytometry and a cytotoxic crossmatch assay. Moreover, we found no evidence of zoonotic transmission from the donor pigs to the human recipients. While substantial additional work will be needed to advance this technology to human trials, these results indicate that pig-to-human heart xenotransplantation can be performed successfully without hyperacute rejection or zoonosis.

BACKGROUND:Crohn's diseases and ulcerative colitis, both of which are chronic immune-mediated disorders of the gastrointestinal tract are major contributors to the overarching Inflammatory bowel diseases. It has become increasingly evident that the pathological processes of IBDs results from interactions between genetic and environmental factors, which can skew immune responses against normal intestinal flora. METHODS:The aim of this study is to assess and analyze the taxa diversity and relative abundances in CD and UC in the Saudi population. We utilized a sequencing strategy that targets all variable regions in the 16 S rRNA gene using the Swift Amplicon 16 S rRNA Panel on Illumina NovaSeq 6000. RESULTS:The composition of stool 16 S rRNA was analyzed from 219 patients with inflammatory bowel disease and from 124 healthy controls. We quantified the abundance of microbial communities to examine any significant differences between subpopulations of samples. At the genus level, two genera in particular, Veillonella and Lachnoclostridium showed significant association with CD versus controls. There were significant differences between subjects with CD versus UC, with the top differential genera spanning Akkermansia, Harryflintia, Maegamonas and Phascolarctobacterium. Furthermore, statistically significant taxa diversity in microbiome composition was observed within the UC and CD groups. CONCLUSIONS:In conclusion we have shown that there are significant differences in gut microbiota between UC, CD and controls in a Saudi Arabian inflammatory bowel disease cohort. This reinforces the need for further studies in large populations that are ethnically and geographically diverse. In addition, our results show the potential to develop classifiers that may have add additional richness of context to clinical diagnosis of UC and CD with larger inflammatory bowel disease cohorts.

BACKGROUND:Endomyocardial biopsy (EMB) is currently considered the gold standard for diagnosing cardiac allograft rejection. However, significant limitations related to histological interpretation variability are well-recognized. We sought to develop a methodology to evaluate EMB solely based on gene expression, without relying on histology interpretation. METHODS:Sixty-four EMBs were obtained from 47 post-heart transplant recipients, who were evaluated for allograft rejection. EMBs were subjected to mRNA sequencing, in which an unsupervised classification algorithm was used to identify the molecular signatures that best classified the EMBs. Cytokine and natriuretic peptide peripheral blood profiling was also performed. Subsequently, we performed gene network analysis to identify the gene modules and gene ontology to understand their biological relevance. We correlated our findings with the unsupervised and histological classifications. RESULTS:Our algorithm classifies EMBs into three categories based solely on clusters of gene expression: unsupervised classes 1, 2, and 3. Unsupervised and histological classifications were closely related, with stronger gene module-phenotype correlations for the unsupervised classes. Gene ontology enrichment analysis revealed processes impacting on the regulation of cardiac and mitochondrial function, immune response, and tissue injury response. Significant levels of cytokines and natriuretic peptides were detected following the unsupervised classification. CONCLUSION/CONCLUSIONS:We have developed an unsupervised algorithm that classifies EMBs into three distinct categories, without relying on histology interpretation. These categories were highly correlated with mitochondrial, immune, and tissue injury response. Significant cytokine and natriuretic peptide levels were detected within the unsupervised classification. If further validated, the unsupervised classification could offer a more objective EMB evaluation.

Ischemic stroke represents a significant societal burden across the globe. Rare high penetrant monogenic variants and less pathogenic common single nucleotide polymorphisms (SNPs) have been described as being associated with risk of diseases. Genetic studies in Saudi Arabian patients offer a greater opportunity to detect rare high penetrant mutations enriched in these consanguineous populations. We performed whole exome sequencing on 387 ischemic stroke subjects from Saudi Arabian hospital networks with up to 20,230 controls from the Saudi Human Genome Project and performed gene burden analyses of variants in 177 a priori loci derived from knowledge-driven curation of monogenic and genome-wide association studies of stroke. Using gene-burden analyses, we observed significant associations in numerous loci under autosomal dominant and/or recessive modelling. Stroke subjects with modified Rankin Scale (mRSs) above 3 were found to carry greater cumulative polygenic risk score (PRS) from rare variants in stroke genes (standardized PRS mean > 0) compared to the population average (standardized PRS mean = 0). However, patients with mRS of 3 or lower had lower cumulative genetic risk from rare variants in stroke genes (OR (95%CI) = 1.79 (1.29-2.49), p = 0.0005), with the means of standardized PRS at or lower than 0. In conclusion, gene burden testing in Saudi stroke populations reveals a number of statistically significant signals under different disease inheritance models. However, interestingly, stroke subjects with mRS of 3 or lower had lower cumulative genetic risk from rare variants in stroke genes and therefore, determining the potential mRS cutoffs to use for clinical significance may allow risk stratification of this population.

Background Cardiometabolic diseases are highly comorbid, but their relationship with female-specific or overwhelmingly female-predominant health conditions (breast cancer, endometriosis, pregnancy complications) is understudied. This study aimed to estimate the cross-trait genetic overlap and influence of genetic burden of cardiometabolic traits on health conditions unique to women. Methods and Results Using electronic health record data from 71 008 ancestrally diverse women, we examined relationships between 23 obstetrical/gynecological conditions and 4 cardiometabolic phenotypes (body mass index, coronary artery disease, type 2 diabetes, and hypertension) by performing 4 analyses: (1) cross-trait genetic correlation analyses to compare genetic architecture, (2) polygenic risk score-based association tests to characterize shared genetic effects on disease risk, (3) Mendelian randomization for significant associations to assess cross-trait causal relationships, and (4) chronology analyses to visualize the timeline of events unique to groups of women with high and low genetic burden for cardiometabolic traits and highlight the disease prevalence in risk groups by age. We observed 27 significant associations between cardiometabolic polygenic scores and obstetrical/gynecological conditions (body mass index and endometrial cancer, body mass index and polycystic ovarian syndrome, type 2 diabetes and gestational diabetes, type 2 diabetes and polycystic ovarian syndrome). Mendelian randomization analysis provided additional evidence of independent causal effects. We also identified an inverse association between coronary artery disease and breast cancer. High cardiometabolic polygenic scores were associated with early development of polycystic ovarian syndrome and gestational hypertension. Conclusions We conclude that polygenic susceptibility to cardiometabolic traits is associated with elevated risk of certain female-specific health conditions.

BACKGROUND:Loss-of-Function (LoF) variants in human genes are important due to their impact on clinical phenotypes and frequent occurrence in the genomes of healthy individuals. The association of LoF variants with complex diseases and traits may lead to the discovery and validation of novel therapeutic targets. Current approaches predict high-confidence LoF variants without identifying the specific genes or the number of copies they affect. Moreover, there is a lack of methods for detecting knockout genes caused by compound heterozygous (CH) LoF variants. RESULTS:We have developed the Loss-of-Function ToolKit (LoFTK), which allows efficient and automated prediction of LoF variants from genotyped, imputed and sequenced genomes. LoFTK enables the identification of genes that are inactive in one or two copies and provides summary statistics for downstream analyses. LoFTK can identify CH LoF variants, which result in LoF genes with two copies lost. Using data from parents and offspring we show that 96% of CH LoF genes predicted by LoFTK in the offspring have the respective alleles donated by each parent. CONCLUSIONS:LoFTK is a command-line based tool that provides a reliable computational workflow for predicting LoF variants from genotyped and sequenced genomes, identifying genes that are inactive in 1 or 2 copies. LoFTK is an open software and is freely available to non-commercial users at https://github.com/CirculatoryHealth/LoFTK .