Lee, Vivian S.

Introduction: Genome-wide association studies (GWAS) have discovered >300 associations for CAD. Few loci are functionally characterized and may represent new mechanisms of disease. Hypothesis: Can an unbiased transcriptional screen for the >2000 candidate GWAS genes in endothelial cells (ECs) identify causal genes and disease-relevant pathways at CAD GWAS loci. Methods: We applied CRISPRi-Perturb-seq to knock down expression of genes within 500 Kb of coronary artery disease GWAS loci (2,300 genes in total) and measure their transcriptomic effects using single-cell RNA-seq. Validation of candidate causal genes was conducted in primary human ECs and zebrafish. Results: We identified 50 programs of co-expressed genes, which represent core cellular pathways (i.e. ribosome biogenesis) and 13 EC-specific pathways such as flow response and angiogenesis. The EC-specific programs show the greatest contribution to CAD heritability. 5 EC-specific programs have the greatest number of GWAS candidate genes. We find that 43 of 306 CAD GWAS signals converge onto these 5 gene programs, and all are regulated by genes in the cerebral cavernous malformations (CCM) pathway—which coordinates KLF2/4 transcriptional responses in endothelial cells, but has not been linked to CAD risk. The strongest regulator of these programs is TLNRD1 , which we show is a new CAD gene and novel regulator of the CCM pathway. TLNRD1 loss-of-function alters actin organization and barrier function in endothelial cells in vitro , and heart development in zebrafish in vivo . Conclusions: High-throughput functional analysis of 2,300 genes proximal to CAD GWAS loci prioritized pathways—such as angiogenesis and EC migration—that are regulated by multiple risk SNPs. Our study identifies new genes that likely influence risk for CAD, identifies convergence of CAD genes into certain pathways in endothelial cells, and prioritizes study of genes in the CCM-signaling pathway for contribution to CAD risk.

Peripheral artery disease (PAD) in the lower extremities often leads to intermittent claudication. In the present study, we proposed a low-dose DCE MRI protocol for quantifying calf muscle perfusion stimulated with plantar flexion and multiple new metrics for interpreting perfusion maps, including the ratio of gastrocnemius over soleus perfusion (G/S; for assessing the vascular redistribution between the two muscles) and muscle perfusion normalized by whole body perfusion (for quantifying the muscle’s active hyperemia). Twenty-eight human subjects participated in this Institutional Review Board-approved study, with 10 healthy subjects ( group A) for assessing interday reproducibility and 8 healthy subjects ( group B) for exploring the relationship between plantar-flexion load and induced muscle perfusion. In a pilot group of five elderly healthy subjects and five patients with PAD ( group C), we proposed a protocol that measured perfusion for a low-intensity exercise and for an exhaustion exercise in a single MRI session. In group A, perfusion estimates for calf muscles were highly reproducible, with correlation coefficients of 0.90–0.93. In group B, gastrocnemius perfusion increased linearly with the exercise workload ( P < 0.05). With the low-intensity exercise, patients with PAD in group C showed substantially lower gastrocnemius perfusion compared with elderly healthy subjects [43.4 (SD 23.5) vs. 106.7 (SD 73.2) ml·min−1·100 g−1]. With exhaustion exercise, G/S [1.0 (SD 0.4)] for patients with PAD was lower than both its low-intensity level [1.9 (SD 1.3)] and the level in elderly healthy subjects [2.7 (SD 2.1)]. In conclusion, the proposed MRI protocol and the new metrics are feasible for quantifying exercise-induced muscle hyperemia, a promising functional test of PAD.NEW & NOTEWORTHY To quantitatively map exercise-induced hyperemia in calf muscles, we proposed a high-resolution MRI method shown to be highly reproducible and sensitive to exercise load. With the use of low contrast, it is feasible to measure calf muscle hyperemia for both low-intensity and exhaustion exercises in a single MRI session. The newly proposed metrics for interpreting perfusion maps are promising for quantifying intermuscle vascular redistribution or a muscle’s active hyperemia.

Blood oxygen level-dependent (BOLD) MRI data of kidney, while indicative of tissue oxygenation level (Po2), is in fact influenced by multiple confounding factors, such as R2, perfusion, oxygen permeability, and hematocrit. We aim to explore the feasibility of extracting tissue Po2 from renal BOLD data. A method of two steps was proposed: first, a Monte Carlo simulation to estimate blood oxygen saturation (SHb) from BOLD signals, and second, an oxygen transit model to convert SHb to tissue Po2. The proposed method was calibrated and validated with 20 pigs (12 before and after furosemide injection) in which BOLD-derived tissue Po2 was compared with microprobe-measured values. The method was then applied to nine healthy human subjects (age: 25.7 ± 3.0 yr) in whom BOLD was performed before and after furosemide. For the 12 pigs before furosemide injection, the proposed model estimated renal tissue Po2 with errors of 2.3 ± 5.2 mmHg (5.8 ± 13.4%) in cortex and −0.1 ± 4.5 mmHg (1.7 ± 18.1%) in medulla, compared with microprobe measurements. After injection of furosemide, the estimation errors were 6.9 ± 3.9 mmHg (14.2 ± 8.4%) for cortex and 2.6 ± 4.0 mmHg (7.7 ± 11.5%) for medulla. In the human subjects, BOLD-derived medullary Po2 increased from 16.0 ± 4.9 mmHg (SHb: 31 ± 11%) at baseline to 26.2 ± 3.1 mmHg (SHb: 53 ± 6%) at 5 min after furosemide injection, while cortical Po2 did not change significantly at ∼58 mmHg (SHb: 92 ± 1%). Our proposed method, validated with a porcine model, appears promising for estimating tissue Po2 from renal BOLD MRI data in human subjects.