Boxiang Liu

Local genetic effects on gene expression across 44 human tissues

Expression quantitative trait locus (eQTL) mapping provides a powerful means to identify functional variants influencing gene expression and disease pathogenesis. We report the identification of cis-eQTLs from 7,051 post-mortem samples representing 44 tissues and 449 individuals as part of the Genotype-Tissue Expression (GTEx) project. We find a cis-eQTL for 88% of all annotated protein-coding genes, with one-third having multiple independent effects. We identify numerous tissue-specific cis-eQTLs, highlighting the unique functional impact of regulatory variation in diverse tissues. By integrating large-scale functional genomics data and state-of-the-art fine-mapping algorithms, we identify multiple features predictive of tissue-specific and shared regulatory effects. We improve estimates of cis-eQTL sharing and effect sizes using allele specific expression across tissues. Finally, we demonstrate the utility of this large compendium of cis-eQTLs for understanding the tissue-specific etiology of complex traits, including coronary artery disease. The GTEx project provides an exceptional resource that has improved our understanding of gene regulation across tissues and the role of regulatory variation in human genetic diseases.

TCF21 and the environmental sensor aryl-hydrocarbon receptor cooperate to activate a pro-inflammatory gene expression program in coronary artery smooth muscle cells

Both environmental factors and genetic loci have been associated with coronary artery disease (CAD), however gene-gene and gene-environment interactions that might identify molecular mechanisms of risk are not easily studied by human genetic approaches. We have previously identified the transcription factor TCF21 as the causal CAD gene at 6q23.2 and characterized its downstream transcriptional network that is enriched for CAD GWAS genes. Here we investigate the hypothesis that TCF21 interacts with a downstream target gene, the aryl hydrocarbon receptor (AHR), a ligand-activated transcription factor that mediates the cellular response to environmental contaminants, including dioxin and polycyclic aromatic hydrocarbons (e.g., tobacco smoke). Perturbation of TCF21 expression in human coronary artery smooth muscle cells (HCASMC) revealed that TCF21 promotes expression of AHR, its heterodimerization partner ARNT, and cooperates with these factors to upregulate a number of inflammatory downstream disease related genes including IL1A, MMP1, and CYP1A1. TCF21 was shown to bind in AHR, ARNT and downstream target gene loci, and co-localization was noted for AHR-ARNT and TCF21 binding sites genome-wide in regions of HCASMC open chromatin. These regions of co-localization were found to be enriched for GWAS signals associated with cardio-metabolic as well as chronic inflammatory disease phenotypes. Finally, we show that similar to TCF21, AHR gene expression is increased in atherosclerotic lesions in mice in vivo using laser capture microdissection, and AHR protein is localized in human carotid atherosclerotic lesions where it is associated with protein kinases with a critical role in innate immune response. These data suggest that TCF21 can cooperate with AHR to activate an inflammatory gene expression program that is exacerbated by environmental stimuli, and may contribute to the overall risk for CAD.

Genetic Regulatory Mechanisms of Smooth Muscle Cells Map to Coronary Artery Disease Risk Loci

Coronary artery disease (CAD) is the leading cause of death globally. Genome-wide association studies (GWASs) have identified more than 95 independent loci that influence CAD risk, most of which reside in non-coding regions of the genome. To interpret these loci, we generated transcriptome and whole-genome datasets using human coronary artery smooth muscle cells (HCASMCs) from 52 unrelated donors, as well as epigenomic datasets using ATAC-seq on a subset of 8 donors. Through systematic comparison with publicly available datasets from GTEx and ENCODE projects, we identified transcriptomic, epigenetic, and genetic regulatory mechanisms specific to HCASMCs. We assessed the relevance of HCASMCs to CAD risk using transcriptomic and epigenomic level analyses. By jointly modeling eQTL and GWAS datasets, we identified five genes (SIPA1, TCF21, SMAD3, FES, and PDGFRA) that may modulate CAD risk through HCASMCs, all of which have relevant functional roles in vascular remodeling. Comparison with GTEx data suggests that SIPA1 and PDGFRA influence CAD risk predominantly through HCASMCs, while other annotated genes may have multiple cell and tissue targets. Together, these results provide tissue-specific and mechanistic insights into the regulation of a critical vascular cell type associated with CAD in human populations.

Ocular disease mechanisms elucidated by genetics of human fetal retinal pigment epithelium gene expression

The eye is an intricate organ with limited representation in large-scale functional genomics datasets. The retinal pigment epithelium (RPE) serves vital roles in ocular development and retinal homeostasis. We interrogated the genetics of gene expression of cultured human fetal RPE (fRPE) cells under two metabolic conditions. Genes with disproportionately high fRPE expression are enriched for genes related to inherited ocular diseases. Variants near these fRPE-selective genes explain a larger fraction of risk for both age-related macular degeneration (AMD) and myopia than variants near genes enriched in 53 other human tissues. Increased mitochondrial oxidation of glutamine by fRPE promoted expression of lipid synthesis genes implicated in AMD. Expression and splice quantitative trait loci (e/sQTL) analysis revealed shared and metabolic condition-specific loci of each type and several eQTL not previously described in any tissue. Fine mapping of fRPE e/sQTL across AMD and myopia genome-wide association data suggests new candidate genes, and mechanisms by which the same common variant of RDH5 contributes to both increased AMD risk and decreased myopia risk. Our study highlights the unique transcriptomic characteristics of fRPE and provides a resource to connect e/sQTL in a critical ocular cell type to monogenic and complex eye disorders.

Identification of rare-disease genes in diverse undiagnosed cases using whole blood transcriptome sequencing and large control cohorts

RNA sequencing (RNA-seq) is a complementary approach for Mendelian disease diagnosis for patients in whom exome-sequencing is not informative. For both rare neuromuscular and mitochondrial disorders, its application has improved diagnostic rates. However, the generalizability of this approach to diverse Mendelian diseases has yet to be evaluated. We sequenced whole blood RNA from 56 cases with undiagnosed rare diseases spanning 11 diverse disease categories to evaluate the general application of RNA-seq to Mendelian disease diagnosis. We developed a robust approach to compare rare disease cases to existing large sets of RNA-seq controls (N=1,594 external and N=31 family-based controls) and demonstrated the substantial impacts of gene and variant filtering strategies on disease gene identification when combined with RNA-seq. Across our cohort, we observed that RNA-seq yields a 8.5% diagnostic rate. These diagnoses included diseases where blood would not intuitively reflect evidence of disease. We identified RARS2 as an under-expression outlier containing compound heterozygous pathogenic variants for an individual exhibiting profound global developmental delay, seizures, microcephaly, hypotonia, and progressive scoliosis. We also identified a new splicing junction in KCTD7 for an individual with global developmental delay, loss of milestones, tremors and seizures. Our study provides a broad evaluation of blood RNA-seq for the diagnosis of rare disease.

Coronary artery disease genes SMAD3 and TCF21 promote opposing interactive genetic programs that regulate smooth muscle cell differentiation and disease risk

Although numerous genetic loci have been associated with coronary artery disease (CAD) with genome wide association studies, efforts are needed to identify the causal genes in these loci and link them into fundamental signaling pathways. Recent studies have investigated the disease mechanism of CAD associated gene SMAD3, a central transcription factor (TF) in the TGFβ pathway, investigating its role in smooth muscle biology. In vitro studies in human coronary artery smooth muscle cells (HCASMC) revealed that SMAD3 modulates cellular phenotype, promoting expression of differentiation marker genes while inhibiting proliferation. RNA sequencing and chromatin immunoprecipitation sequencing studies in HCASMC identified downstream genes that reside in pathways which mediate vascular development and atherosclerosis processes in this cell type. HCASMC phenotype, and gene expression patterns promoted by SMAD3 were noted to have opposing direction of effect compared to another CAD associated TF, TCF21. At sites of SMAD3 and TCF21 colocalization on DNA, SMAD3 binding was inversely correlated with TCF21 binding, due in part to TCF21 locally blocking chromatin accessibility at the SMAD3 binding site. Further, TCF21 was able to directly inhibit SMAD3 activation of gene expression in transfection reporter gene studies. In contrast to TCF21 which is protective toward CAD, SMAD3 expression in HCASMC was shown to be directly correlated with disease risk. We propose that the pro-differentiation action of SMAD3 inhibits dedifferentiation that is required for HCASMC to expand and stabilize disease plaque as they respond to vascular stresses, counteracting the protective dedifferentiating activity of TCF21 and promoting disease risk.

Recurrently Mutated Genes Differ between Leptomeningeal and Solid Lung Cancer Brain Metastases

&NA; When compared with solid brain metastases from NSCLC, leptomeningeal disease (LMD) has unique growth patterns and is rapidly fatal. Patients with LMD do not undergo surgical resection, limiting the tissue available for scientific research. In this study we performed whole exome sequencing on eight samples of LMD to identify somatic mutations and compared the results with those for 26 solid brain metastases. We found that taste 2 receptor member 31 gene (TAS2R31) and phosphodiesterase 4D interacting protein gene (PDE4DIP) were recurrently mutated among LMD samples, suggesting involvement in LMD progression. Together with a retrospective review of the charts of an additional 44 patients with NSCLC LMD, we discovered a surprisingly low number of KRAS mutations (n = 4 [7.7%]) but a high number of EGFR mutations (n = 33 [63.5%]). The median interval for development of LMD from NSCLC was shorter in patients with mutant EGFR (16.3 months) than in patients with wild‐type EGFR (23.9 months) (p = 0.017). Targeted analysis of recurrent mutations thus presents a useful complement to the existing diagnostic tool kit, and correlations of EGFR in LMD and KRAS in solid metastases suggest that molecular distinctions or systemic treatment pressure underpin the differences in growth patterns within the brain.

Ubiquitination of ABCE1 by NOT4 in Response to Mitochondrial Damage Links Co-translational Quality Control to PINK1-Directed Mitophagy

Translation of mRNAs is tightly regulated and constantly surveyed for errors. Aberrant translation can trigger co-translational protein and RNA quality control processes, impairments of which cause neurodegeneration by still poorly understood mechanism(s). Here we show that quality control of translation of mitochondrial outer membrane (MOM)-localized mRNA intersects with the turnover of damaged mitochondria, both orchestrated by the mitochondrial kinase PINK1. Mitochondrial damage causes stalled translation of complex-I 30 kDa subunit (C-I30) mRNA on MOM, triggering the recruitment of co-translational quality control factors Pelo, ABCE1, and NOT4 to the ribosome/mRNA-ribonucleoprotein complex. Damage-induced ubiquitination of ABCE1 by NOT4 generates poly-ubiquitin signals that attract autophagy receptors to MOM to initiate mitophagy. In the Drosophila PINK1 model, these factors act synergistically to restore mitophagy and neuromuscular tissue integrity. Thus ribosome-associated co-translational quality control generates an early signal to trigger mitophagy. Our results have broad therapeutic implications for the understanding and treatment of neurodegenerative diseases.