Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

Genomic profiling of human vascular cells identifies TWIST1 as a causal gene for common vascular diseases.

Genome-wide association studies have identified multiple novel genomic loci associated with vascular diseases. Many of these loci are common non-coding variants that affect the expression of disease-relevant genes within coronary vascular cells. To identify such genes on a genome-wide level, we performed deep transcriptomic analysis of genotyped primary human coronary artery smooth muscle cells (HCASMCs) and coronary endothelial cells (HCAECs) from the same subjects, including splicing Quantitative Trait Loci (sQTL), allele-specific expression (ASE), and colocalization analyses. We identified sQTLs for TARS2, YAP1, CFDP1, and STAT6 in HCASMCs and HCAECs, and 233 ASE genes, a subset of which are also GTEx eGenes in arterial tissues. Colocalization of GWAS association signals for coronary artery disease (CAD), migraine, stroke and abdominal aortic aneurysm with GTEx eGenes in aorta, coronary artery and tibial artery discovered novel candidate risk genes for these diseases. At the CAD and stroke locus tagged by rs2107595 we demonstrate colocalization with expression of the proximal gene TWIST1. We show that disrupting the rs2107595 locus alters TWIST1 expression and that the risk allele has increased binding of the NOTCH signaling protein RBPJ. Finally, we provide data that TWIST1 expression influences vascular SMC phenotypes, including proliferation and calcification, as a potential mechanism supporting a role for TWIST1 in CAD.

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 elements 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 effects of genetic variation on 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 74 GWAS loci that have been associated with at least one complex phenotype. Our results reveal a 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.

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.

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.

SPAG7 deletion causes intrauterine growth restriction, resulting in adulthood obesity and metabolic dysfunction.

From a forward mutagenetic screen to discover mutations associated with obesity, we identified mutations in the Spag7 gene linked to metabolic dysfunction in mice. Here, we show that SPAG7 KO mice are born smaller and develop obesity and glucose intolerance in adulthood. This obesity does not stem from hyperphagia, but a decrease in energy expenditure. The KO animals also display reduced exercise tolerance and muscle function due to impaired mitochondrial function. Furthermore, SPAG7-deficiency in developing embryos leads to intrauterine growth restriction, brought on by placental insufficiency, likely due to abnormal development of the placental junctional zone. This insufficiency leads to loss of SPAG7-deficient fetuses in utero and reduced birth weights of those that survive. We hypothesize that a 'thrifty phenotype' is ingrained in SPAG7 KO animals during development that leads to adult obesity. Collectively, these results indicate that SPAG7 is essential for embryonic development and energy homeostasis later in life.

Obesity rates are climbing worldwide, leading to an increase in associated conditions such as type 2 diabetes. While new pharmaceutical approaches are available to help individuals manage their weight, many patients do not respond to them or experience prohibitive side effects. Identifying alternative treatments will likely require pinpointing the genes and molecular actors involved in the biological processes that control weight regulation. Previous research suggests that a protein known as SPAG7 could help shape how mice use and store the energy they extract from food. Flaherty et al. therefore set out to investigate the role this protein plays in the body. To do so, they created a line of mice born without SPAG7, which they monitored closely throughout life. These animals were underweight at birth and did not eat more than other mice, yet they were obese as adults. Their ability to exercise was reduced, their muscles were weaker and contained fibers with functional defects. The mice also exhibited biological changes associated with the onset of diabetes. Yet deleting SPAG7 during adulthood led to no such changes; these mice maintained normal muscle function and body weight. Closely examining how SPAG7-deficient mice developed in the womb revealed placental defects which likely caused these animals to receive fewer nutrients from their mother. Such early-life deprivation is known to be associated with the body shifting towards maximizing its use of resources and privileging fat storage, even into and throughout adulthood. By shedding light on the biological role of SPAG7, the work by Flaherty et al. helps to better understand how developmental events can increase the likelihood of obesity later in life. Further investigations are now needed to explore whether this knowledge could help design interventions relevant to human health.

Session Introduction: SALUD: Scalable Applications of cLinical risk Utility and preDiction.

This PSB 2023 session discusses challenges in clinical implication and application of risk prediction models, which includes but is not limited to: implementation of risk models, responsible use of polygenic risk scores (PGS), and other risk prediction strategies. We focus on the development and use of new, scalable methods for harmonizing and refining risk prediction models by incorporating genetic and non-genetic risk factors, applying new phenotyping strategies, and integrating clinical factors and biomarkers. Lastly, we will discuss innovation in expanding the utility of these prediction models to underrepresented populations. This session focuses on the overarching theme of enabling early diagnosis, and treatment and preventive measures related to complex diseases and comorbidities.

Evaluating the Contribution of Cell Type-Specific Alternative Splicing to Variation in Lipid Levels.

BACKGROUND: Genome-wide association studies have identified hundreds of loci associated with lipid levels. However, the genetic mechanisms underlying most of these loci are not well-understood. Recent work indicates that changes in the abundance of alternatively spliced transcripts contribute to complex trait variation. Consequently, identifying genetic loci that associate with alternative splicing in disease-relevant cell types and determining the degree to which these loci are informative for lipid biology is of broad interest. METHODS: We analyze gene splicing in 83 sample-matched induced pluripotent stem cell (iPSC) and hepatocyte-like cell lines (n=166), as well as in an independent collection of primary liver tissues (n=96) to perform discovery of splicing quantitative trait loci (sQTLs). RESULTS: We observe that transcript splicing is highly cell type specific, and the genes that are differentially spliced between iPSCs and hepatocyte-like cells are enriched for metabolism pathway annotations. We identify 1384 hepatocyte-like cell sQTLs and 1455 iPSC sQTLs at a false discovery rate of <5% and find that sQTLs are often shared across cell types. To evaluate the contribution of sQTLs to variation in lipid levels, we conduct colocalization analysis using lipid genome-wide association data. We identify 19 lipid-associated loci that colocalize either with an hepatocyte-like cell expression quantitative trait locus or sQTL. Only 2 loci colocalize with both a sQTL and expression quantitative trait locus, indicating that sQTLs contribute information about genome-wide association studies loci that cannot be obtained by analysis of steady-state gene expression alone. CONCLUSIONS: These results provide an important foundation for future efforts that use iPSC and iPSC-derived cells to evaluate genetic mechanisms influencing both cardiovascular disease risk and complex traits in general.

Identification and Development of Therapeutics for COVID-19.

After emerging in China in late 2019, the novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spread worldwide, and as of mid-2021, it remains a significant threat globally. Only a few coronaviruses are known to infect humans, and only two cause infections similar in severity to SARS-CoV-2: Severe acute respiratory syndrome-related coronavirus, a species closely related to SARS-CoV-2 that emerged in 2002, and Middle East respiratory syndrome-related coronavirus, which emerged in 2012. Unlike the current pandemic, previous epidemics were controlled rapidly through public health measures, but the body of research investigating severe acute respiratory syndrome and Middle East respiratory syndrome has proven valuable for identifying approaches to treating and preventing novel coronavirus disease 2019 (COVID-19). Building on this research, the medical and scientific communities have responded rapidly to the COVID-19 crisis and identified many candidate therapeutics. The approaches used to identify candidates fall into four main categories: adaptation of clinical approaches to diseases with related pathologies, adaptation based on virological properties, adaptation based on host response, and data-driven identification (ID) of candidates based on physical properties or on pharmacological compendia. To date, a small number of therapeutics have already been authorized by regulatory agencies such as the Food and Drug Administration (FDA), while most remain under investigation. The scale of the COVID-19 crisis offers a rare opportunity to collect data on the effects of candidate therapeutics. This information provides insight not only into the management of coronavirus diseases but also into the relative success of different approaches to identifying candidate therapeutics against an emerging disease. IMPORTANCE The COVID-19 pandemic is a rapidly evolving crisis. With the worldwide scientific community shifting focus onto the SARS-CoV-2 virus and COVID-19, a large number of possible pharmaceutical approaches for treatment and prevention have been proposed. What was known about each of these potential interventions evolved rapidly throughout 2020 and 2021. This fast-paced area of research provides important insight into how the ongoing pandemic can be managed and also demonstrates the power of interdisciplinary collaboration to rapidly understand a virus and match its characteristics with existing or novel pharmaceuticals. As illustrated by the continued threat of viral epidemics during the current millennium, a rapid and strategic response to emerging viral threats can save lives. In this review, we explore how different modes of identifying candidate therapeutics have borne out during COVID-19.

Identification and Development of Therapeutics for COVID-19.

After emerging in China in late 2019, the novel coronavirus SARS-CoV-2 spread worldwide and as of mid-2021 remains a significant threat globally. Only a few coronaviruses are known to infect humans, and only two cause infections similar in severity to SARS-CoV-2: Severe acute respiratory syndrome-related coronavirus, a closely related species of SARS-CoV-2 that emerged in 2002, and Middle East respiratory syndrome-related coronavirus, which emerged in 2012. Unlike the current pandemic, previous epidemics were controlled rapidly through public health measures, but the body of research investigating severe acute respiratory syndrome and Middle East respiratory syndrome has proven valuable for identifying approaches to treating and preventing novel coronavirus disease 2019 (COVID-19). Building on this research, the medical and scientific communities have responded rapidly to the COVID-19 crisis to identify many candidate therapeutics. The approaches used to identify candidates fall into four main categories: adaptation of clinical approaches to diseases with related pathologies, adaptation based on virological properties, adaptation based on host response, and data-driven identification of candidates based on physical properties or on pharmacological compendia. To date, a small number of therapeutics have already been authorized by regulatory agencies such as the Food and Drug Administration (FDA), while most remain under investigation. The scale of the COVID-19 crisis offers a rare opportunity to collect data on the effects of candidate therapeutics. This information provides insight not only into the management of coronavirus diseases, but also into the relative success of different approaches to identifying candidate therapeutics against an emerging disease.