YoSon Park

Transposable elements are the primary source of novelty in primate gene regulation

Gene regulation shapes the evolution of phenotypic diversity. We investigated the evolution of liver promoters and enhancers in six primate species using ChIP-seq (H3K27ac and H3K4me1) to profile cis-regulatory elements (CREs) and using RNA-seq to characterize gene expression in the same individuals. To quantify regulatory divergence, we compared CRE activity across species by testing differential ChIP-seq read depths directly measured for orthologous sequences. We show that the primate regulatory landscape is largely conserved across the lineage, with 63% of the tested human liver CREs showing similar activity across species. Conserved CRE function is associated with sequence conservation, proximity to coding genes, cell-type specificity, and transcription factor binding. Newly evolved CREs are enriched in immune response and neurodevelopmental functions. We further demonstrate that conserved CREs bind master regulators, suggesting that while CREs contribute to species adaptation to the environment, core functions remain intact. Newly evolved CREs are enriched in young transposable elements (TEs), including Long-Terminal-Repeats (LTRs) and SINE-VNTR-Alus (SVAs), that significantly affect gene expression. Conversely, only 16% of conserved CREs overlap TEs. We tested the cis-regulatory activity of 69 TE subfamilies by luciferase reporter assays, spanning all major TE classes, and showed that 95.6% of tested TEs can function as either transcriptional activators or repressors. In conclusion, we demonstrated the critical role of TEs in primate gene regulation and illustrated potential mechanisms underlying evolutionary divergence among the primate species through the noncoding genome.

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.

Transposable element exaptation is the primary source of novelty in the primate gene regulatory landscape

Gene regulation plays a critical role in the evolution of phenotypic diversity. We investigated the evolution of liver promoters and enhancers in six primate species. We performed ChIP-seq for two histone modifications and RNA-seq to profile cis-regulatory element (CRE) activity and gene expression. The primate regulatory landscape is largely conserved across the lineage. Conserved CRE function is associated with sequence conservation, proximity to coding genes, cell type specificity of CRE function, and transcription factor binding. Newly evolved CREs are enriched in immune response and neurodevelopmental functions, while conserved CREs bind master regulators. Transposable elements (TEs) are the primary source of novelty in primate gene regulation. Newly evolved CREs are enriched in young TEs that affect gene expression. However, only 17% of conserved CREs overlap a TE, suggesting that target gene expression is under strong selection. Finally, we identified specific genomic features driving the functional recruitment of newly inserted TEs.

Genetic And Epigenetic Fine Mapping Of Complex Trait Associated Loci In The Human Liver

Deciphering the impact of genetic variation on gene regulation is fundamental to understanding common, complex human diseases. Although histone modifications are important markers of gene regulatory regions of the genome, any specific histone modification has not been assayed in more than a few individuals in the human liver. As a result, the impacts of genetic variation that direct histone modification states in the liver are poorly understood. Here, we generate the most comprehensive genome-wide dataset of two epigenetic marks, H3K4me3 and H3K27ac, and annotate thousands of putative regulatory elements in the human liver. We integrate these findings with genome-wide gene expression data collected from the same human liver tissues and high-resolution promoter-focused chromatin interaction maps collected from human liver-derived HepG2 cells. We demonstrate widespread functional consequences of natural genetic variation on putative regulatory element activity and gene expression levels. Leveraging these extensive datasets, we fine-map a total of 77 GWAS loci that have been associated with at least one complex phenotype. Our results contribute to the repertoire of genes and regulatory mechanisms governing complex disease development and further the basic understanding of genetic and epigenetic regulation of gene expression in the human liver tissue.

RNA-binding protein A1CF modulates plasma triglyceride levels through posttranscriptional regulation of stress-induced VLDL secretion

Background A recent human exome-chip study on plasma lipids identified a missense mutation in the A1CF (APOBEC1 complementation factor) gene that is associated with elevated triglyceride (TG) levels, but how A1CF, an RNA binding protein, influences plasma TG is unknown. Methods We generated A1cf knockout (A1cf −/−) mice and knock-in mice homozygous for the TG-associated Gly398Ser mutation (A1cfGS/GS), determined lipid phenotypes, and assessed TG physiology through measurements of clearance and secretion. We further identified A1CF’s RNA binding targets using enhanced cross-linking and immunoprecipitation sequencing of cultured HepG2 cells and investigated pathways enriched for these targets. Transcriptomic effects of A1CF deficiency were evaluated through RNA sequencing and analyses for differential expression, alternative splicing, and RNA editing. Results Both A1cf −/−and A1cfGS/GS mice exhibited increased fasting plasma TG, establishing that the TG phenotype is due to A1CF loss of function. In vivo TG secretion and clearance studies revealed increased TG secretion without changes in clearance in A1cf −/−mice. Increased VLDL-apoB secretion was also seen in A1cf −/−rat hepatoma cells, but no increase in apoB synthesis was observed. This phenotype was seen without significant shifts in apoB-100/apoB-48 in A1CF deficiency. To discover novel pathways for A1CF’s role in TG metabolism, we identified A1CF’s RNA binding targets, which were enriched for pathways related to proteasomal catabolism and endoplasmic reticulum (ER) stress. Indeed, proteasomal inhibition led to increased cellular stress in A1cf −/−cells, and higher expression of ER-stress protein GRP78 was observed in resting A1cf −/−cells. RNA-seq of whole livers from wild-type and A1cf −/−mice revealed that pro-inflammatory, not lipogenesis, genes were upregulated as a secondary effect of A1CF deficiency. Differential alternative splicing (AS) analysis and RNA editing analysis revealed that genes involved in cellular stress and metabolism underwent differential changes in A1CF deficiency, and top A1CF binding target proteins with relevance to intracellular stress were differentially expressed on the protein but not mRNA level, implicating multiple mechanisms by which A1CF influences TG secretion. Conclusions These data suggest an important role for A1CF in mediating VLDL-TG secretion through regulating intracellular stress.

Renal compartment–specific genetic variation analyses identify new pathways in chronic kidney disease

Chronic kidney disease (CKD), a condition in which the kidneys are unable to clear waste products, affects 700 million people globally. Genome-wide association studies (GWASs) have identified sequence variants for CKD; however, the biological basis of these GWAS results remains poorly understood. To address this issue, we created an expression quantitative trait loci (eQTL) atlas for the glomerular and tubular compartments of the human kidney. Through integrating the CKD GWAS with eQTL, single-cell RNA sequencing and regulatory region maps, we identified novel genes for CKD. Putative causal genes were enriched for proximal tubule expression and endolysosomal function, where DAB2, an adaptor protein in the TGF-β pathway, formed a central node. Functional experiments confirmed that reducing Dab2 expression in renal tubules protected mice from CKD. In conclusion, compartment-specific eQTL analysis is an important avenue for the identification of novel genes and cellular pathways involved in CKD development and thus potential new opportunities for its treatment.Kidney compartment–specific eQTL analysis goes beyond GWAS to reveal causal genes and pathways involved in renal disease development.

The YscE/YscG chaperone and YscF N-terminal sequences target YscF to the Yersinia pestis type III secretion apparatus

The needle structures of type III secretion (T3S) systems are formed by the secretion and polymerization of a needle subunit protein, YscF in Yersinia pestis. A subset of T3S systems employ unique heterodimeric chaperones, YscE and YscG in Y. pestis, to prevent the polymerization of needle subunits within the bacterial cell. We demonstrate that the YscE/YscG chaperone is also required for stable YscF expression and for secretion of YscF. Overexpression of a functional maltose-binding protein (MBP)-YscG hybrid protein stabilized cytoplasmic YscF but YscF was not secreted in the absence of YscE. Furthermore, a YscE mutant protein was identified that functioned with YscG to stabilize cytosolic YscF; however, YscF was not secreted. These findings confirm a role for the YscE/YscG chaperone in YscF secretion and suggest that YscE may have a specific role in this process. Recent studies have shown that YscF deleted of its N-terminal 15 residues is still secreted and functional, suggesting that YscF may not require an N-terminal secretion signal. However, we demonstrate that YscF contains an N-terminal secretion signal and that a functional N-terminal signal is required for YscF secretion.