Genomic signature of trait-­associated variants

Genome-­‐wide association studies have been used extensively to study hundreds of phenotypes and have determined thousands of associated SNPs whose underlying biology and causation is as yet largely unknown. Many previous studies attempted to clarify the causal biology by investigating overlaps of trait-­‐ associated variants with functional annotations, but lacked statistical rigor and examined incomplete subsets of available functional annotations. Additionally, it has been difficult to disentangle the relative contributions of different annotations that may show strong correlations with one another. In this thesis, we address these shortcomings and strengthen and extend the obtained results. Two methods, permutations and logistic regression, are applied in statistically rigorous analyses of genomic annotations and their observed enrichment or depletion of trait-­‐associated SNPs. The genomic annotations range from genic regions and regulatory features to measures of conservation and aspects of chromatin structure. Logistic regressions in a number of trait-­‐specific subsets identify genomic annotations influencing SNPs associated with both normal variation (e.g., eye or hair colour) and diseases, suggesting some generalities in the biological underpinnings of phenotypes. SNPs associated with phenotypes of the immune system are investigated and the results highlight the distinct aetiology for this subset. Despite the heterogeneity of the studied cancers, SNPs associated to different cancers are particularly enriched for conserved regions, unlike all other trait-­‐subsets. Nonetheless, chromatin states are, perhaps surprisingly, among the most influential genomic annotations in all trait-­‐ subsets. Evolutionary conserved regions are rarely within the top genomic annotations despite their widespread use in prioritisation methods for follow-­‐ up studies. We identify a common set of enriched or depleted genomic annotations that significantly influence all traits, but also highlight trait-­‐specific differences. These annotations may be used for the computational prioritisation of variants implicated in phenotypes of interest. The approaches developed for this thesis are further applied to studies of a specific human complex trait (height) and gene expression in atherosclerosis

The genomic signature of trait-associated variants

BACKGROUND: Genome-wide association studies have identified thousands of SNP variants associated with hundreds of phenotypes. For most associations the causal variants and the molecular mechanisms underlying pathogenesis remain unknown. Exploration of the underlying functional annotations of trait-associated loci has thrown some light on their potential roles in pathogenesis. However, there are some shortcomings of the methods used to date, which may undermine efforts to prioritize variants for further analyses. Here, we introduce and apply novel methods to rigorously identify annotation classes showing enrichment or depletion of trait-associated variants taking into account the underlying associations due to co-location of different functional annotations and linkage disequilibrium. RESULTS: We assessed enrichment and depletion of variants in publicly available annotation classes such as genic regions, regulatory features, measures of conservation, and patterns of histone modifications. We used logistic regression to build a multivariate model that identified the most influential functional annotations for trait-association status of genome-wide significant variants. SNPs associated with all of the enriched annotations were 8 times more likely to be trait-associated variants than SNPs annotated with none of them. Annotations associated with chromatin state together with prior knowledge of the existence of a local expression QTL (eQTL) were the most important factors in the final logistic regression model. Surprisingly, despite the widespread use of evolutionary conservation to prioritize variants for study we find only modest enrichment of trait-associated SNPs in conserved regions. CONCLUSION: We established odds ratios of functional annotations that are more likely to contain significantly trait-associated SNPs, for the purpose of prioritizing GWAS hits for further studies. Additionally, we estimated the relative and combined influence of the different genomic annotations, which may facilitate future prioritization methods by adding substantial information

Opening sequence:computational genomics in the era of high-throughput sequencing

A report on the 11th Cold Spring Harbor Laboratory/Wellcome Trust conference on Genome Informatics, Cold Spring Harbor Laboratories, New York, USA, November 2-5, 2011

Biomarkers for assessing pain and pain relief in the neonatal intensive care unit

Newborns admitted to the neonatal intensive care unit (NICU) regularly undergo painful procedures and may face various painful conditions such as postoperative pain. Optimal management of pain in these vulnerable preterm and term born neonates is crucial to ensure their comfort and prevent negative consequences of neonatal pain. This entails accurate and timely identification of pain, non-pharmacological pain treatment and if needed administration of analgesic therapy, evaluation of treatment effectiveness, and monitoring of adverse effects. Despite the widely recognized importance of pain management, pain assessment in neonates has thus far proven to be a challenge. As self-report, the gold standard for pain assessment, is not possible in neonates, other methods are needed. Several observational pain scales have been developed, but these often rely on snapshot and largely subjective observations and may fail to capture pain in certain conditions. Incorporation of biomarkers alongside observational pain scores holds promise in enhancing pain assessment and, by extension, optimizing pain treatment and neonatal outcomes. This review explores the possibilities of integrating biomarkers in pain assessment in the NICU.</p

Delineating morbidity patterns in preterm infants at near-term age using a data-driven approach

BACKGROUND: Long-term survival after premature birth is significantly determined by development of morbidities, primarily affecting the cardio-respiratory or central nervous system. Existing studies are limited to pairwise morbidity associations, thereby lacking a holistic understanding of morbidity co-occurrence and respective risk profiles. METHODS: Our study, for the first time, aimed at delineating and characterizing morbidity profiles at near-term age and investigated the most prevalent morbidities in preterm infants: bronchopulmonary dysplasia (BPD), pulmonary hypertension (PH), mild cardiac defects, perinatal brain pathology and retinopathy of prematurity (ROP). For analysis, we employed two independent, prospective cohorts, comprising a total of 530 very preterm infants: AIRR ("Attention to Infants at Respiratory Risks") and NEuroSIS ("Neonatal European Study of Inhaled Steroids"). Using a data-driven strategy, we successfully characterized morbidity profiles of preterm infants in a stepwise approach and (1) quantified pairwise morbidity correlations, (2) assessed the discriminatory power of BPD (complemented by imaging-based structural and functional lung phenotyping) in relation to these morbidities, (3) investigated collective co-occurrence patterns, and (4) identified infant subgroups who share similar morbidity profiles using machine learning techniques. RESULTS: First, we showed that, in line with pathophysiologic understanding, BPD and ROP have the highest pairwise correlation, followed by BPD and PH as well as BPD and mild cardiac defects. Second, we revealed that BPD exhibits only limited capacity in discriminating morbidity occurrence, despite its prevalence and clinical indication as a driver of comorbidities. Further, we demonstrated that structural and functional lung phenotyping did not exhibit higher association with morbidity severity than BPD. Lastly, we identified patient clusters that share similar morbidity patterns using machine learning in AIRR (n=6 clusters) and NEuroSIS (n=8 clusters). CONCLUSIONS: By capturing correlations as well as more complex morbidity relations, we provided a comprehensive characterization of morbidity profiles at discharge, linked to shared disease pathophysiology. Future studies could benefit from identifying risk profiles to thereby develop personalized monitoring strategies

Abundant local interactions in the 4p16.1 region suggest functional mechanisms underlying SLC2A9 associations with human serum uric acid

Human serum uric acid concentration (SUA) is a complex trait. A recent meta-analysis of multiple genome-wide association studies (GWAS) identified 28 loci associated with SUA jointly explaining only 7.7% of the SUA variance, with 3.4% explained by two major loci (SLC2A9 and ABCG2). Here we examined whether gene–gene interactions had any roles in regulating SUA using two large GWAS cohorts included in the meta-analysis [the Atherosclerosis Risk in Communities study cohort (ARIC) and the Framingham Heart Study cohort (FHS)]. We found abundant genome-wide significant local interactions in ARIC in the 4p16.1 region located mostly in an intergenic area near SLC2A9 that were not driven by linkage disequilibrium and were replicated in FHS. Taking the forward selection approach, we constructed a model of five SNPs with marginal effects and three epistatic SNP pairs in ARIC—three marginal SNPs were located within SLC2A9 and the remaining SNPs were all located in the nearby intergenic area. The full model explained 1.5% more SUA variance than that explained by the lead SNP alone, but only 0.3% was contributed by the marginal and epistatic effects of the SNPs in the intergenic area. Functional analysis revealed strong evidence that the epistatically interacting SNPs in the intergenic area were unusually enriched at enhancers active in ENCODE hepatic (HepG2, P = 4.7E−05) and precursor red blood (K562, P = 5.0E−06) cells, putatively regulating transcription of WDR1 and SLC2A9. These results suggest that exploring epistatic interactions is valuable in uncovering the complex functional mechanisms underlying the 4p16.1 region

Association of Altered Plasma Lipidome with Disease Severity in COVID-19 Patients

The severity of COVID-19 is linked to an imbalanced immune response. The dysregulated metabolism of small molecules and bioactive lipids has also been associated with disease severity. To promote understanding of the disease biochemistry and provide targets for intervention, we applied a range of LC-MS platforms to analyze over 100 plasma samples from patients with varying COVID-19 severity and with detailed clinical information on inflammatory responses (&gt;30 immune markers). This is the third publication in a series, and it reports the results of comprehensive lipidome profiling using targeted LC-MS/MS. We identified 1076 lipid features across 25 subclasses, including glycerophospholipids, sterols, glycerolipids, and sphingolipids, among which 531 lipid features were dramatically changed in the plasma of intensive care unit (ICU) patients compared to patients in the ward. Patients in the ICU showed 1.3–57-fold increases in ceramides, (lyso-)glycerophospholipids, diglycerides, triglycerides, and plasmagen phosphoethanolamines, and 1.3–2-fold lower levels of a cyclic lysophosphatidic acid, sphingosine-1-phosphates, sphingomyelins, arachidonic acid-containing phospholipids, lactosylceramide, and cholesterol esters compared to patients in the ward. Specifically, phosphatidylinositols (PIs) showed strong fatty acid saturation-dependent behavior, with saturated fatty acid (SFA)- and monosaturated fatty acid (MUFA)-derived PI decreasing and polystaturated (PUFA)-derived PI increasing. We also found ~4000 significant Spearman correlations between lipids and multiple clinical markers of immune response with |R| ≥ 0.35 and FDR corrected Q &lt; 0.05. Except for lysophosphatidic acid, lysophospholipids were positively associated with the CD4 fraction of T cells, and the cytokines IL-8 and IL-18. In contrast, sphingosine-1-phosphates were negatively correlated with innate immune markers such as CRP and IL-6. Further indications of metabolic changes in moderate COVID-19 disease were demonstrated in recovering ward patients compared to those at the start of hospitalization, where 99 lipid species were altered (6 increased by 30–62%; 93 decreased by 1.3–2.8-fold). Overall, these findings support and expand on early reports that dysregulated lipid metabolism is involved in COVID-19.</p

Loss of adipose triglyceride lipase is associated with human cancer and induces mouse pulmonary neoplasia

Metabolic reprogramming is a hallmark of cancer. Understanding cancer metabolism is instrumental to devise innovative therapeutic approaches. Anabolic metabolism, including the induction of lipogenic enzymes, is a key feature of proliferating cells. Here, we report a novel tumor suppressive function for adipose triglyceride lipase (ATGL), the rate limiting enzyme in the triglyceride hydrolysis cascade. In immunohistochemical analysis, non-small cell lung cancers, pancreatic adenocarcinoma as well as leiomyosarcoma showed significantly reduced levels of ATGL protein compared to corresponding normal tissues. The ATGL gene was frequently deleted in various forms of cancers. Low levels of ATGL mRNA correlated with significantly reduced survival in patients with ovarian, breast, gastric and non-small cell lung cancers. Remarkably, pulmonary neoplasia including invasive adenocarcinoma developed spontaneously in mice lacking ATGL pointing to an important role for this lipase in controlling tumor development. Loss of ATGL, as detected in several forms of human cancer, induces spontaneous development of pulmonary neoplasia in a mouse model. Our results, therefore, suggest a novel tumor suppressor function for ATGL and contribute to the understanding of cancer metabolism. We propose to evaluate loss of ATGL protein expression for the diagnosis of malignant tumors. Finally, modulation of the lipolytic pathway may represent a novel therapeutic approach in the treatment of human cancer

Targeted proteomics and metabolomics for biomarker discovery in abdominal aortic aneurysm and post-EVAR sac volume

BACKGROUND AND AIMS: Abdominal aortic aneurysm (AAA) patients undergo uniform surveillance programs both leading up to, and following surgery. Circulating biomarkers could play a pivotal role in individualizing surveillance. We applied a multi-omics approach to identify relevant biomarkers and gain pathophysiological insights. MATERIALS AND METHODS: In this cross-sectional study, 108 AAA patients and 200 post-endovascular aneurysm repair (post-EVAR) patients were separately investigated. We performed partial least squares regression and ingenuity pathway analysis on circulating concentrations of 96 proteins (92 Olink Cardiovascular-III panel, 4 ELISA-assays) and 199 metabolites (measured by LC-TQMS), and their associations with CT-based AAA/sac volume. RESULTS: The median (25th-75th percentile) maximal diameter was 50.0 mm (46.0, 53.0) in the AAA group, and 55.4 mm (45.0, 64.2) in the post-EVAR group. Correcting for clinical characteristics in AAA patients, the aneurysm volume Z-score differed 0.068 (95 %CI: (0.042, 0.093)), 0.066 (0.047, 0.085) and -0.051 (-0.064, -0.038) per Z-score valine, leucine and uPA, respectively. After correcting for clinical characteristics and orthogonalization in the post-EVAR group, the sac volume Z-score differed 0.049 (0.034, 0.063) per Z-score TIMP-4, -0.050 (-0.064, -0.037) per Z-score LDL-receptor, -0.051 (-0.062, -0.040) per Z-score 1-OG/2-OG and -0.056 (-0.066, -0.045) per Z-score 1-LG/2-LG. CONCLUSIONS: The branched-chain amino acids and uPA were related to AAA volume. For post-EVAR patients, LDL-receptor, monoacylglycerols and TIMP-4 are potential biomarkers for sac volume. Additionally, distinct markers for sac change were identified.</p

Targeted proteomics and metabolomics for biomarker discovery in abdominal aortic aneurysm and post-EVAR sac volume

BACKGROUND AND AIMS: Abdominal aortic aneurysm (AAA) patients undergo uniform surveillance programs both leading up to, and following surgery. Circulating biomarkers could play a pivotal role in individualizing surveillance. We applied a multi-omics approach to identify relevant biomarkers and gain pathophysiological insights. MATERIALS AND METHODS: In this cross-sectional study, 108 AAA patients and 200 post-endovascular aneurysm repair (post-EVAR) patients were separately investigated. We performed partial least squares regression and ingenuity pathway analysis on circulating concentrations of 96 proteins (92 Olink Cardiovascular-III panel, 4 ELISA-assays) and 199 metabolites (measured by LC-TQMS), and their associations with CT-based AAA/sac volume. RESULTS: The median (25th-75th percentile) maximal diameter was 50.0 mm (46.0, 53.0) in the AAA group, and 55.4 mm (45.0, 64.2) in the post-EVAR group. Correcting for clinical characteristics in AAA patients, the aneurysm volume Z-score differed 0.068 (95 %CI: (0.042, 0.093)), 0.066 (0.047, 0.085) and -0.051 (-0.064, -0.038) per Z-score valine, leucine and uPA, respectively. After correcting for clinical characteristics and orthogonalization in the post-EVAR group, the sac volume Z-score differed 0.049 (0.034, 0.063) per Z-score TIMP-4, -0.050 (-0.064, -0.037) per Z-score LDL-receptor, -0.051 (-0.062, -0.040) per Z-score 1-OG/2-OG and -0.056 (-0.066, -0.045) per Z-score 1-LG/2-LG. CONCLUSIONS: The branched-chain amino acids and uPA were related to AAA volume. For post-EVAR patients, LDL-receptor, monoacylglycerols and TIMP-4 are potential biomarkers for sac volume. Additionally, distinct markers for sac change were identified.</p