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.

BACKGROUND:Xenotransplantation of genetically engineered porcine organs has the potential to address the challenge of organ donor shortage. Two cases of porcine-to-human kidney xenotransplantation were performed, yet the physiological effects on the xenografts and the recipients' immune responses remain largely uncharacterized. METHODS:We performed single-cell RNA sequencing (scRNA-seq) and longitudinal RNA-seq analyses of the porcine kidneys to dissect xenotransplantation-associated cellular dynamics and xenograft-recipient interactions. We additionally performed longitudinal scRNA-seq of the peripheral blood mononuclear cells (PBMCs) to detect recipient immune responses across time. FINDINGS/RESULTS:Although no hyperacute rejection signals were detected, scRNA-seq analyses of the xenografts found evidence of endothelial cell and immune response activation, indicating early signs of antibody-mediated rejection. Tracing the cells' species origin, we found human immune cell infiltration in both xenografts. Human transcripts in the longitudinal bulk RNA-seq revealed that human immune cell infiltration and the activation of interferon-gamma-induced chemokine expression occurred by 12 and 48 h post-xenotransplantation, respectively. Concordantly, longitudinal scRNA-seq of PBMCs also revealed two phases of the recipients' immune responses at 12 and 48-53 h. Lastly, we observed global expression signatures of xenotransplantation-associated kidney tissue damage in the xenografts. Surprisingly, we detected a rapid increase of proliferative cells in both xenografts, indicating the activation of the porcine tissue repair program. CONCLUSIONS:Longitudinal and single-cell transcriptomic analyses of porcine kidneys and the recipient's PBMCs revealed time-resolved cellular dynamics of xenograft-recipient interactions during xenotransplantation. These cues can be leveraged for designing gene edits and immunosuppression regimens to optimize xenotransplantation outcomes. FUNDING/BACKGROUND:This work was supported by NIH RM1HG009491 and DP5OD033430.

BACKGROUND:Hydrocephalus is a highly heterogeneous multifactorial disease that arises from genetic and environmental factors. Familial genetic studies of hydrocephalus have elucidated four robustly associated hydrocephalus associated loci. This study aims to identify potential genetic causation in cases of hydrocephalus, with or without spina bifida and Dandy Walker Syndrome (DWS), using family-based rare variant association analysis of whole exome sequencing. METHODS:We performed whole exome sequencing in 143 individuals across 48 families where at least one offspring was affected with hydrocephalus (N.=27), with hydrocephalus with spina bifida (N.=21) and with DWS (N.=3), using Illumina HiSeq 2500 instrument. RESULTS:No pathogenic or putative pathogenic single-nucleotide variants were evident in the four known hydrocephalus loci in our subjects. However, after examining 73 known hydrocephalus genes previously identified from literature, we identified three potentially impactful variants from the cohort. Using a gene panel comprising variants in known neural tube defects loci, we identified a total of 1024 potentially deleterious variants, of which 797 were missense variants and 191 were frameshift variants, 36 were stop gain/loss variants. A small portion of our family pedigree analyses yielded putative genetic signals which may be responsible for hydrocephaly elated phenotypes, however the low diagnostic yield may be due to lack of capture of genetic variants in the exonic regions i.e. structural variants may only be evident from whole genome sequencing. CONCLUSIONS:We identified three potentially impactful variants from our cohort in 73 known hydrocephalus genes previously identified in literature.

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.

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.

INTRODUCTION/UNASSIGNED:Genome-wide association studies have discovered common polymorphisms in regions associated with schizophrenia. No genome-wide analyses have been performed in Saudi schizophrenia subjects. METHODS/UNASSIGNED:Genome-wide genotyping data from 136 Saudi schizophrenia cases and 97 Saudi controls in addition to 4,625 American were examined for copy number variants (CNVs). A hidden Markov model approach was used to call CNVs. RESULTS/UNASSIGNED: = 0.04). The analyses focused on extremely large >250 kilobases CNVs or homozygous deletions of any size. One extremely large deletion was noted in a single case (16.5 megabases on chromosome 10). Two cases had an 814 kb duplication of chromosome 7 spanning a cluster of genes, including circadian-related loci, and two other cases had 277 kb deletions of chromosome 9 encompassing an olfactory receptors gene family. CNVs were also seen in loci previously associated with schizophrenia, namely a 16p11 proximal duplication and two 22q11.2 deletions. DISCUSSION/UNASSIGNED:Runs of homozygosity (ROHs) were analyzed across the genome to investigate correlation with schizophrenia risk. While rates and sizes of these ROHs were similar in cases and controls, we identified 10 regions where multiple cases had ROHs and controls did not.

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.

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.