Search

Atherosclerotic disease is a major cause of acute cardiovascular events. A deeper understanding of its underlying mechanisms will allow advancing personalized and patient-centered healthcare. Transcriptomic research has proven to be a powerful tool for unravelling the complex molecular pathways that drive atherosclerosis. However, low reproducibility of research findings and lack of standardization of procedures pose significant challenges in this field. In this review, we discuss how transcriptomic research can help in understanding the different phenotypes of the atherosclerotic plaque that contribute to the development and progression of atherosclerosis. We highlight the methodological challenges that need to be addressed to improve research outputs, and emphasize the importance of research protocols harmonization. We also discuss recent advances in transcriptomic research, including bulk or single-cell sequencing, and their added value in plaque phenotyping. Finally, we explore how integrated multiomics data and machine learning improve understanding of atherosclerosis and provide directions for future research.

Background: Cladribine is a highly effective immunotherapy that is applied in two short-term courses over 2 years and reduces relapse rate and disease progression in patients with relapsing multiple sclerosis (MS). Despite the short treatment period, cladribine has a long-lasting effect on disease activity even after recovery of lymphocyte counts, suggesting a yet undefined long-term immune modulating effect. Objectives: Our aim was to provide a more profound understanding of the detailed effects of cladribine, also with regard to the patients’ therapy response. Design: We performed an open-labeled, explorative, prospective, single-arm study, in which we examined the detailed lymphocyte subset development of MS patients who received cladribine treatment over 2 years. Methods: We performed in-depth profiling of the effects of cladribine on peripheral blood lymphocytes by flow cytometry, bulk RNA sequencing of sorted CD4 + T cells, CD8 + T cells, and CD19 + B cells as well as single-cell RNA sequencing of peripheral blood mononuclear cells in a total of 23 MS patients before and at different time points up to 24 months after cladribine treatment. Data were correlated with clinical and cranial magnetic resonance imaging (MRI) disease activity. Results: Flow cytometry revealed a predominant and sustained reduction of memory B cells compared to other B cell subsets after cladribine treatment, whereas T cell subsets were slightly reduced in a more uniform pattern. The overall transcriptional profile of total blood B cells exhibited reduced expression of proinflammatory and T cell activating genes, while single-cell transcriptomics revealed that gene expression within each B cell cluster did not change over time. Stable patients displayed stronger reductions of selected memory B cell clusters as compared to patients with clinical or cerebral MRI disease activity. Conclusion: We describe a pronounced and sustained effect of cladribine on the memory B cell compartment, and the resulting change in B cell subset composition causes a significant alteration of B cell transcriptional profiles resulting in reduced proinflammatory and T cell activating capacities. The extent of reduction in selected memory B cell clusters by cladribine may predict treatment response.

ABSTRACTInfluenza A viruses (IAV) cause annual epidemics and occasional pandemics in humans. The most recent pandemic outbreak occurred in 2009 with H1N1pdm09. This virus, which most likely reassorted in swine before its transmission to humans, was reintroduced into the swine population and continues circulating ever since. In order to assess its potential to cause reassortants on a cellular level, human origin H1N1pdm09 and a recent Eurasian avian-like H1N1 swine IAV were (co-)passaged in the newly generated swine lung cell line C22. Co-infection with both viruses gave rise to numerous reassortants that additionally carry different mutations which can partially be found in nature as well. Reassortment most frequently affected the PB1, PA and NA segments with the swine IAV as recipient. These reassortants reached higher titers in swine lung cells and were able to replicate in genuine human lung tissue explants ex vivo, suggesting a possible zoonotic potential. Interestingly, reassortment and mutations in the viral ribonucleoprotein complex influence the viral polymerase activity in a cell type and species-specific manner. In summary, we demonstrate reassortment promiscuity of these viruses in a novel swine lung cell model and indicate a possible zoonotic potential of the reassortants.

The spatial and frontal QRS-T angles are electrocardiographic (ECG) predictors for arrhythmia. This work used genetic analyses to identify associated loci and pathways, and explore their relationships with other ECG traits and cardiovascular disease.

The QT interval is a heritable electrocardiographic measure associated with arrhythmia risk when prolonged. Here, the authors used a series of genetic analyses to identify genetic loci, pathways, therapeutic targets, and relationships with cardiovascular disease.

DNA methylation is a reversible epigenetic mechanism of transcriptional regulation. Changes in DNA methylation patterns are influenced by internal and external factors, and have been associated with various human diseases. We hypothesized that specific methylation changes occur in depressed individuals, and that these changes are associated with disease progression and responses to pharmacological treatments. To test our hypothesis, we analyzed DNA methylation-wide changes in the BiDirect Study's depression and control cohorts. Over the course of two consecutive follow-up visits in 3-year intervals from baseline (i.e. 3 and 6 years after initial assessment), we identified a set of changes in methylation sites within gene enhancer and promoter regions that consistently changed only in depression patients or only in control individuals, respectively. Longitudinal changes in these methylation sites correlated with changes in circulating markers of inflammation and were enriched in biological processes related to innate and adaptive immunity. Currently, we are starting to explore the translational value of these findings by investigating their relationships with depressive symptoms and immunity. Moreover, by applying different machine learning algorithms, we test the utility of a prioritized subset of the identified methylation sites to predict disease status, disease progression, chronic inflammation status and antidepressant response. Overall, we expect to leverage longitudinal DNA methylation information to gain insights into the role of immunity in depression and generate computational tools with potential applications in the clinic.

IntroductionThe measurement of neurofilament light chain (NfL) in blood is a promising biomarker of neurological injury and disease. We investigated the genetic factors that underlie serum NfL levels (sNfL) of individuals without neurological conditions.MethodsWe performed a discovery genome-wide association study (GWAS) of sNfL in participants of the German BiDirect Study (N = 1,899). A secondary GWAS for meta-analysis was performed in a small Austrian cohort (N = 287). Results from the meta-analysis were investigated in relation with several clinical variables in BiDirect.ResultsOur discovery GWAS identified 12 genomic loci at the suggestive threshold ((p < 1 × 10−5). After meta-analysis, 7 loci were suggestive of an association with sNfL. Genotype-specific differences in sNfL were observed for the lead variants of meta-analysis loci (rs34523114, rs114956339, rs529938, rs73198093, rs34372929, rs10982883, and rs1842909) in BiDirect participants. We identified potential associations in meta-analysis loci with markers of inflammation and renal function. At least 6 protein-coding genes (ACTG2, TPRKB, DMXL1, COL23A1, NAT1, and RIMS2) were suggested as genetic factors contributing to baseline sNfL levels.DiscussionOur findings suggest that polygenic regulation of neuronal processes, inflammation, metabolism and clearance modulate the variability of NfL in the circulation. These could aid in the interpretation of sNfL measurements in a personalized manner.

Abstract Coronary artery disease (CAD) is a long-lasting inflammatory disease characterized by monocyte migration into the vessel wall leading to clinical events like myocardial infarction (MI). However, the role of monocyte subsets, especially their miRNA-driven differentiation in this scenario is still in its infancy. Here, we characterized monocyte subsets in controls and disease phenotypes of CAD and MI patients using flow cytometry and miRNA and mRNA expression profiling using RNA sequencing. We observed major differences in the miRNA profiles between the classical (CD14++CD16−) and nonclassical (CD14+CD16++) monocyte subsets irrespective of the disease phenotype suggesting the Cyclin-dependent Kinase 6 (CDK6) to be an important player in monocyte maturation. Between control and MI patients, we found a set of miRNAs to be differentially expressed in the nonclassical monocytes and targeting CCND2 (Cyclin D2) that is able to enhance myocardial repair. Interestingly, miRNAs as miR-125b playing a role in vascular calcification were differentially expressed in the classical subset in patients suffering from CAD and not MI in comparison to control samples. In conclusion, our study describes specific peculiarities of monocyte subset miRNA expression in control and diseased samples and provides basis to further functional analysis and to identify new cardiovascular disease treatment targets.

Abstract Background The mTOR-PI3K-Akt pathway influences cell metabolism and (malignant) cell growth. We generated sex-specific polygenic risk scores capturing natural variation in 7 out of 10 top-ranked genes in this pathway. We studied the scores directly and in interaction with energy balance-related factors (body mass index (BMI), trouser/skirt size, height, physical activity, and early life energy restriction) in relation to colorectal cancer (CRC) risk in the Netherlands Cohort Study (NLCS) (n=120,852). The NLCS has a case-cohort design and 20.3 years of follow-up. Participants completed a baseline questionnaire on diet and cancer in 1986 when 55–69 years old. ~75% of the cohort returned toenail clippings used for DNA isolation and genotyping (n subcohort=3,793, n cases=3,464). To generate the scores, the dataset was split in two and risk alleles were defined and weighted based on sex-specific associations with CRC risk in the other dataset half, because there were no SNPs in the top-ranked genes associated with CRC risk in previous genome-wide association studies at a significance level p

Background: The amygdala is crucial for emotional cognitive processing. Affective or emotional states can bias cognitive processes, including attention, memory, and decision-making. This can result in optimistic or pessimistic behaviors that are partially driven by the activation of the amygdala. The resulting emotional cognitive bias is a common feature of anxiety and mood disorders, both of which are interactively influenced by genetic and environmental factors. It is also known that emotional cognitive biases can be influenced by environmental factors. However, little is known about the effects of genetics and/or gene-environment interactions on emotional cognitive biases. We investigated the effects of the genetic background and environmental enrichment on the transcriptional profiles of the mouse amygdala following a well-established cognitive bias test.Methods: Twenty-four female C57BL/6J and B6D2F1N mice were housed either in standard (control) conditions or in an enriched environment. After appropriate training, the cognitive bias test was performed on 19 mice that satisfactorily completed the training scheme to assess their responses to ambiguous cues. This allowed us to calculate an “optimism score” for each mouse. Subsequently, we dissected the anterior and posterior portions of the amygdala to perform RNA-sequencing for differential expression and other statistical analyses.Results: In general, we found only minor changes in the amygdala’s transcriptome associated with the levels of optimism in our mice. In contrast, we observed wide molecular effects of the genetic background in both housing environments. The C57BL/6J animals showed more transcriptional changes in response to enriched environments than the B6D2F1N mice. We also generally found more dysregulated genes in the posterior than in the anterior portion of the amygdala. Gene set overrepresentation analyses consistently implicated cellular metabolic responses and immune processes in the differences observed between mouse strains, while processes favoring neurogenesis and neurotransmission were implicated in the responses to environmental enrichment. In a correlation analysis, lipid metabolism in the anterior amygdala was suggested to influence the levels of optimism.Conclusions: Our observations underscore the importance of selecting appropriate animal models when performing molecular studies of affective conditions or emotional states, and suggest an important role of immune and stress responses in the genetic component of emotion regulation.

BackgroundDuring early phases of life, such as prenatal or early postnatal development and adolescence, an organism's phenotype can be shaped by the environmental conditions it experiences. According to the Match-Mismatch hypothesis (MMH), changes to this environment during later life stages can result in a mismatch between the individual's adaptations and the prevailing environmental conditions. Thus, negative consequences in welfare and health can occur. We aimed to test the MMH in the context of food availability, assuming adolescence as a sensitive period of adaptation.MethodsWe have previously reported a study of the physiological and behavioral effects of match and mismatch conditions of high (ad libitum) and low (90% of ad libitum intake) food availability from adolescence to early adulthood in female C57BL/6J mice (n = 62). Here, we performed RNA-sequencing of the livers of a subset of these animals (n = 16) to test the effects of match and mismatch feeding conditions on the liver transcriptome.ResultsIn general, we found no effect of the match-mismatch situations. Contrarily, the amount of food available during early adulthood (low vs. high) drove the differences we observed in final body weight and gene expression in the liver, regardless of the amount of food available to the animals during adolescence. Many of the differentially expressed genes and the corresponding biological processes found to be overrepresented overlapped, implicating common changes in various domains. These included metabolism, homeostasis, cellular responses to diverse stimuli, transport of bile acids and other molecules, cell differentiation, major urinary proteins, and immunity and inflammation.ConclusionsOur previous and present observations found no support for the MMH in the context of low vs high food availability from adolescence to early adulthood in female C57BL/6J mice. However, even small differences of approximately 10% in food availability during early adulthood resulted in physiological and molecular changes with potential beneficial implications for metabolic diseases.

The proinflammatory alarmins S100A8 and S100A9 are among the most abundant proteins in neutrophils and monocytes but are completely silenced after differentiation to macrophages. The molecular mechanisms of the extraordinarily dynamic transcriptional regulation of S100a8 and S100a9 genes, however, are only barely understood. Using an unbiased genome-wide CRISPR/Cas9 knockout (KO)-based screening approach in immortalized murine monocytes, we identified the transcription factor C/EBPδ as a central regulator of S100a8 and S100a9 expression. We showed that S100A8/A9 expression and thereby neutrophil recruitment and cytokine release were decreased in C/EBPδ KO mice in a mouse model of acute lung inflammation. S100a8 and S100a9 expression was further controlled by the C/EBPδ antagonists ATF3 and FBXW7. We confirmed the clinical relevance of this regulatory network in subpopulations of human monocytes in a clinical cohort of cardiovascular patients. Moreover, we identified specific C/EBPδ-binding sites within S100a8 and S100a9 promoter regions, and demonstrated that C/EBPδ-dependent JMJD3-mediated demethylation of H3K27me3 is indispensable for their expression. Overall, our work uncovered C/EBPδ as a novel regulator of S100a8 and S100a9 expression. Therefore, C/EBPδ represents a promising target for modulation of inflammatory conditions that are characterized by S100a8 and S100a9 overexpression.

Background: Pathogens or trauma-derived danger signals induced maturation and activation of plasmacytoid dendritic cells (pDCs) is a pivotal step in pDC-dependent host defense. Exposure of pDC to cardiometabolic disease-associated lipids and proteins may well influence critical signaling pathways, thereby compromising immune responses against endogenous, bacterial and viral pathogens. In this study, we have addressed if hyperlipidemia impacts human pDC activation, cytokine response and capacity to prime CD4+ T cells. METHODS AND RESULTS: We show that exposure to pro-atherogenic oxidized low-density lipoproteins (oxLDL) led to pDC lipid accumulation, which in turn ablated a Toll-like receptor (TLR) 7 and 9 dependent up-regulation of pDC maturation markers CD40, CD83, CD86 and HLA-DR. Moreover, oxLDL dampened TLR9 activation induced the production of pro-inflammatory cytokines in a NUR77/IRF7 dependent manner and impaired the capacity of pDCs to prime and polarize CD4+ T helper (Th) cells. CONCLUSION: Our findings reveal profound effects of dyslipidemia on pDC responses to pathogen-derived signals.

BACKGROUND Genomic and experimental studies suggest a role for PITX2 in atrial fibrillation (AF). To assess if this association is relevant for recurrent AF in patients, we tested whether left atrial PITX2 affects recurrent AF after AF ablation.METHODS mRNA concentrations of PITX2 and its cardiac isoform, PITX2c, were quantified in left atrial appendages (LAAs) from patients undergoing thoracoscopic AF ablation, either in whole LAA tissue (n = 83) or in LAA cardiomyocytes (n = 52), and combined with clinical parameters to predict AF recurrence. Literature suggests that BMP10 is a PITX2-repressed, atrial-specific, secreted protein. BMP10 plasma concentrations were combined with 11 cardiovascular biomarkers and clinical parameters to predict recurrent AF after catheter ablation in 359 patients.RESULTS Reduced concentrations of cardiomyocyte PITX2, but not whole LAA tissue PITX2, were associated with AF recurrence after thoracoscopic AF ablation (16% decreased recurrence per 2–(ΔΔCt) increase in PITX2). RNA sequencing, quantitative PCR, and Western blotting confirmed that BMP10 is one of the most PITX2-repressed atrial genes. Left atrial size (HR per mm increase [95% CI], 1.055 [1.028, 1.082]); nonparoxysmal AF (HR 1.672 [1.206, 2.318]), and elevated BMP10 (HR 1.339 [CI 1.159, 1.546] per quartile increase) were predictive of recurrent AF. BMP10 outperformed 11 other cardiovascular biomarkers in predicting recurrent AF.CONCLUSIONS Reduced left atrial cardiomyocyte PITX2 and elevated plasma concentrations of the PITX2-repressed, secreted atrial protein BMP10 identify patients at risk of recurrent AF after ablation.TRIAL REGISTRATION ClinicalTrials.gov NCT01091389, NL50069.018.14, Dutch National Registry of Clinical Research Projects EK494-16.FUNDING British Heart Foundation, European Union (H2020), Leducq Foundation.

We recently reported a family-based genome wide association study (GWAS) for pediatric stroke pointing our attention to two significantly associated genes of the ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs) gene family ADAMTS2 (rs469568, p = 8x10-6) and ADAMTS12 (rs1364044, p = 2.9x10-6). To further investigate these candidate genes, we applied a targeted resequencing approach on 48 discordant sib-pairs for pediatric stroke followed by genotyping of the detected non-synonymous variants in the full cohort of 270 offspring trios and subsequent fine mapping analysis. We identified eight non-synonymous SNPs in ADAMTS2 and six in ADAMTS12 potentially influencing the respective protein function. These variants were genotyped within a cohort of 270 affected offspring trios, association analysis revealed the ADAMTS12 variant rs77581578 to be significantly under-transmitted (p = 6.26x10-3) to pediatric stroke patients. The finding was validated in a pediatric venous thromboembolism (VTE) cohort of 189 affected trios. Subsequent haplotype analysis of ADAMTS12 detected a significantly associated haplotype comprising the originally identified GWAS variant. Several ADAMTS genes such as ADAMTS13 are involved in thromboembolic disease process. Here, we provide further evidence for ADAMTS12 to likely play a role in pediatric stroke. Further functional studies are warranted to assess the functional role of ADAMTS12 in the pathogenesis of stroke.

Long non-coding RNAs (lncRNAs) have been increasingly studied during the past decade. This led to an immense number of annotated transcripts, out of which many were linked to a diverse range of biological mechanisms and diseases. Due to the variety of their regulatory potential, they are seen as an important link in understanding complex epigenetic mechanisms. Prominent examples of lncRNAs in the cardiovascular system are ANRIL, Braveheart, MALAT1 and HOTAIR which have been excessively studied. But despite the impressive number of described transcripts, only a few examples are characterized functionally. One way to do this is to identify accessible structural domains in the RNA secondary structure which have the ability to bind to DNA, RNA or proteins. Through recent improvements in computational as well as experimental methods, this exploration of secondary structure became not only more efficient than traditional methods like crystallization, but also feasible to investigate whole genome RNA structures. The purpose of this review is to highlight the recent advances in secondary structure probing methods and how these can be applied in order to investigate the functional roles of lncRNAs in the cardiovascular system.

Abstract Background Genetic variation is an important determinant of RNA transcription and splicing, which in turn contributes to variation in human traits, including cardiovascular diseases. Results Here we report the first in-depth survey of heart transcriptome variation using RNA-sequencing in 97 patients with dilated cardiomyopathy and 108 non-diseased controls. We reveal extensive differences of gene expression and splicing between dilated cardiomyopathy patients and controls, affecting known as well as novel dilated cardiomyopathy genes. Moreover, we show a widespread effect of genetic variation on the regulation of transcription, isoform usage, and allele-specific expression. Systematic annotation of genome-wide association SNPs identifies 60 functional candidate genes for heart phenotypes, representing 20% of all published heart genome-wide association loci. Focusing on the dilated cardiomyopathy phenotype we found that eQTL variants are also enriched for dilated cardiomyopathy genome-wide association signals in two independent cohorts. Conclusions RNA transcription, splicing, and allele-specific expression are each important determinants of the dilated cardiomyopathy phenotype and are controlled by genetic factors. Our results represent a powerful resource for the field of cardiovascular genetics.

Unraveling the genetic susceptibility of complex diseases such as chronic kidney disease remains challenging. Here, we used inbred rat models of kidney damage associated with elevated blood pressure for the comprehensive analysis of a major albuminuria susceptibility locus detected in these models. We characterized its genomic architecture by congenic substitution mapping, targeted next-generation sequencing, and compartment-specific RNA sequencing analysis in isolated glomeruli. This led to prioritization of transmembrane protein Tmem63c as a novel potential target. Tmem63c is differentially expressed in glomeruli of allele-specific rat models during onset of albuminuria. Patients with focal segmental glomerulosclerosis exhibited specific TMEM63C loss in podocytes. Functional analysis in zebrafish revealed a role for tmem63c in mediating the glomerular filtration barrier function. Our data demonstrate that integrative analysis of the genomic architecture of a complex trait locus is a powerful tool for identification of new targets such as Tmem63c for further translational investigation.

With the rising interest in the regulatory functions of long non-coding RNAs (lncRNAs) in complex human diseases such as cardiovascular diseases, there is an increasing need in public databases offering comprehensive and integrative data for all aspects of these versatile molecules. Recently, a variety of public data repositories that specialized in lncRNAs have been developed, which make use of huge high-throughput data particularly from next-generation sequencing (NGS) approaches. Here, we provide an overview of current lncRNA databases covering basic and functional annotation, lncRNA expression and regulation, interactions with other biomolecules, and genomic variants influencing the structure and function of lncRNAs. The prominent lncRNA antisense noncoding RNA in the INK4 locus (ANRIL), which has been unequivocally associated with coronary artery disease through genome-wide association studies (GWAS), serves as an example to demonstrate the features of each individual database.

Cardiovascular diseases (CVDs) are of multifactorial origin and can be attributed to several genetic and environmental components. CVDs are the leading cause of mortality worldwide and they primarily damage the heart and the vascular system. Non-coding RNA (ncRNA) refers to functional RNA molecules, which have been transcribed into DNA but do not further get translated into proteins. Recent transcriptomic studies have identified the presence of thousands of ncRNA molecules across species. In humans, less than 2% of the total genome represents the protein-coding genes. While the role of many ncRNAs is yet to be ascertained, some long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) have been associated with disease progression, serving as useful diagnostic and prognostic biomarkers. A plethora of data repositories specialized in ncRNAs have been developed over the years using publicly available high-throughput data from next-generation sequencing and other approaches, that cover various facets of ncRNA research like basic and functional annotation, expressional profile, structural and molecular changes, and interaction with other biomolecules. Here, we provide a compendium of the current ncRNA databases relevant to cardiovascular research.

The EU-CardioRNA Cooperation in Science and Technology (COST) Action is a European-wide consortium established in 2018 with 31 European country members and four associate member countries to build bridges between translational researchers from academia and industry who conduct research on non-coding RNAs, cardiovascular diseases and similar research areas. EU-CardioRNA comprises four core working groups (WG1–4). In the first year since its launch, EU-CardioRNA met biannually to exchange and discuss recent findings in related fields of scientific research, with scientific sessions broadly divided up according to WG. These meetings are also an opportunity to establish interdisciplinary discussion groups, brainstorm ideas and make plans to apply for joint research grants and conduct other scientific activities, including knowledge transfer. Following its launch in Brussels in 2018, three WG meetings have taken place. The first of these in Lisbon, Portugal, the second in Istanbul, Turkey, and the most recent in Maastricht, The Netherlands. Each meeting includes a scientific session from each WG. This meeting report briefly describes the highlights and key take-home messages from each WG session in this first successful year of the EU-CardioRNA COST Action.

The transcription factor ISL1 is thought to be key for conveying the multipotent and proliferative properties of cardiac precursor cells. Here, we investigate its function upon cardiac induction of human embryonic stem cells. We find that ISL1 does not stabilize the transient cardiac precursor cell state but rather serves to accelerate cardiomyocyte differentiation. Conversely, ISL1 depletion delays cardiac differentiation and respecifies nascent cardiomyocytes from a ventricular to an atrial identity. Mechanistic analyses integrate this unrecognized anti-atrial function of ISL1 with known and newly identified atrial inducers. In this revised view, ISL1 is antagonized by retinoic acid signaling via a novel player, MEIS2. Conversely, ISL1 competes with the retinoic acid pathway for prospective cardiomyocyte fate, which converges on the atrial specifier NR2F1. This study reveals a core regulatory network putatively controlling human heart chamber formation and also bears implications for the subtype-specific production of human cardiomyocytes with enhanced functional properties.

Cardiovascular disease (CVD) remains the leading cause of death worldwide and, despite continuous advances, better diagnostic and prognostic tools, as well as therapy, are needed. The human transcriptome, which is the set of all RNA produced in a cell, is much more complex than previously thought and the lack of dialogue between researchers and industrials and consensus on guidelines to generate data make it harder to compare and reproduce results. This European Cooperation in Science and Technology (COST) Action aims to accelerate the understanding of transcriptomics in CVD and further the translation of experimental data into usable applications to improve personalized medicine in this field by creating an interdisciplinary network. It aims to provide opportunities for collaboration between stakeholders from complementary backgrounds, allowing the functions of different RNAs and their interactions to be more rapidly deciphered in the cardiovascular context for translation into the clinic, thus fostering personalized medicine and meeting a current public health challenge. Thus, this Action will advance studies on cardiovascular transcriptomics, generate innovative projects, and consolidate the leadership of European research groups in the field. COST (European Cooperation in Science and Technology) is a funding organization for research and innovation networks (www.cost.eu).

Cardiovascular diseases (CVDs) affect the heart and the vascular system with a high prevalence and place a huge burden on society as well as the healthcare system. These complex diseases are often the result of multiple genetic and environmental risk factors and pose a great challenge to understanding their etiology and consequences. With the advent of next generation sequencing, many non-coding RNA transcripts, especially long non-coding RNAs (lncRNAs), have been linked to the pathogenesis of CVD. Despite increasing evidence, the proper functional characterization of most of these molecules is still lacking. The exploration of conservation of sequences across related species has been used to functionally annotate protein coding genes. In contrast, the rapid evolutionary turnover and weak sequence conservation of lncRNAs make it difficult to characterize functional homologs for these sequences. Recent studies have tried to explore other dimensions of interspecies conservation to elucidate the functional role of these novel transcripts. In this review, we summarize various methodologies adopted to explore the evolutionary conservation of cardiovascular non-coding RNAs at sequence, secondary structure, syntenic, and expression level.

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25-55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38-75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.

Reactivation of fetal gene expression patterns has been implicated in common cardiac diseases in adult life including left ventricular (LV) hypertrophy (LVH) in arterial hypertension. Thus, increased wall stress and neurohumoral activation are discussed to induce the return to expression of fetal genes after birth in LVH. We therefore aimed to identify novel potential candidates for LVH by analyzing fetal-adult cardiac gene expression in a genetic rat model of hypertension, i.e. the stroke-prone spontaneously hypertensive rat (SHRSP). To this end we performed genome-wide transcriptome analysis in SHRSP to identify differences in expression patterns between day 20 of fetal development (E20) and adult animals in week 14 in comparison to a normotensive rat strain with contrasting low LV mass, i.e. Fischer (F344). 15232 probes were detected as expressed in LV tissue obtained from rats at E20 and week 14 (p < 0.05) and subsequently screened for differential expression. We identified 24 genes with SHRSP specific up-regulation and 21 genes with down-regulation as compared to F344. Further bioinformatic analysis presented Efcab6 as a new candidate for LVH that showed only in the hypertensive SHRSP rat differential expression during development (logFC = 2.41, p < 0.001) and was significantly higher expressed in adult SHRSP rats compared with adult F344 (+ 76%) and adult normotensive Wistar-Kyoto rats (+ 82%). Thus, it represents an interesting new target for further functional analyses and the elucidation of mechanisms leading to LVH. Here we report a new approach to identify candidate genes for cardiac hypertrophy by combining the analysis of gene expression differences between strains with a contrasting cardiac phenotype with a comparison of fetal-adult cardiac expression patterns.

Recent studies reported altered DNA methylation in failing human hearts. This may suggest a role for de novo DNA methylation in the development of heart failure. Here, we tested whether cardiomyocyte-specific loss of de novo DNA methyltransferases Dnmt3a and Dnmt3b altered cardiac function and remodeling after chronic left ventricular pressure overload.Mice with specific ablation of Dnmt3a and Dnmt3b expression in cardiomyocytes were generated by crossing floxed Dnmt3afl and Dnmt3bfl alleles with mice expressing Cre recombinase under control of the atrial myosin light chain gene promoter. The efficacy of combined Dnmt3a/3b ablation (DKO) was characterized on cardiomyocyte-specific genomic DNA and mRNA levels. Cardiac phenotyping was carried out without (sham) or with left ventricular pressure overload induced by transverse aortic constriction (TAC). Under similar conditions, cardiac genome-wide transcriptional profiling was performed and DNA methylation levels of promoters of differentially regulated genes were assessed by pyrosequencing.DKO cardiomyocytes showed virtual absence of targeted Dnmt3a and Dnmt3b mRNA transcripts. Cardiac phenotyping revealed no significant differences between DKO and control mice under sham and TAC conditions. Transcriptome analyses identified upregulation of 44 and downregulation of 9 genes in DKO as compared with control sham mice. TAC mice showed similar changes with substantial overlap of regulated genes compared to sham. Promoters of upregulated genes were largely unmethylated in DKO compared to control mice.The absence of cardiac pathology in the presence of the predicted molecular phenotype suggests that de novo DNA methylation in cardiomyocytes is dispensable for adaptive mechanisms after chronic cardiac pressure overload.

Duplicate genes emerge as copy-number variations (CNVs) at the population level, and remain copy-number polymorphic until they are fixed or lost. The successful establishment of such structural polymorphisms in the genome plays an important role in evolution by promoting genetic diversity, complexity and innovation. To characterize the early evolutionary stages of duplicate genes and their potential adaptive benefits, we combine comparative genomics with population genomics analyses to evaluate the distribution and impact of CNVs across natural populations of an eco-genomic model, the three-spined stickleback. With whole genome sequences of 66 individuals from populations inhabiting three distinct habitats, we find that CNVs generally occur at low frequencies and are often only found in one of the 11 populations surveyed. A subset of CNVs, however, displays copy-number differentiation between populations, showing elevated within-population frequencies consistent with local adaptation. By comparing teleost genomes to identify lineage-specific genes and duplications in sticklebacks, we highlight rampant gene content differences among individuals in which over 30% of young duplicate genes are CNVs. These CNV genes are evolving rapidly at the molecular level and are enriched with functional categories associated with environmental interactions, depicting the dynamic early copy-number polymorphic stage of genes during population differentiation.

Vertebrate innate immunity is the first line of defense against an invading pathogen and has long been assumed to be largely unspecific with respect to parasite/pathogen species. However, recent phenotypic evidence suggests that immunogenetic variation, i.e. allelic variability in genes associated with the immune system, results in host-parasite genotype-by-genotype interactions and thus specific innate immune responses. Immunogenetic variation is common in all vertebrate taxa and this reflects an effective immunological function in complex environments. However, the underlying variability in host gene expression patterns as response of innate immunity to within-species genetic diversity of macroparasites in vertebrates is unknown. We hypothesized that intra-specific variation among parasite genotypes must be reflected in host gene expression patterns. Here we used high-throughput RNA-sequencing to examine the effect of parasite genotypes on gene expression patterns of a vertebrate host, the three-spined stickleback (Gasterosteus aculeatus). By infecting naïve fish with distinct trematode genotypes of the species Diplostomum pseudospathaceum we show that gene activity of innate immunity in three-spined sticklebacks depended on the identity of an infecting macroparasite genotype. In addition to a suite of genes indicative for a general response against the trematode we also find parasite-strain specific gene expression, in particular in the complement system genes, despite similar infection rates of single clone treatments. The observed discrepancy between infection rates and gene expression indicates the presence of alternative pathways which execute similar functions. This suggests that the innate immune system can induce redundant responses specific to parasite genotypes.

We present a comprehensive toolkit for post-processing, visualization and advanced analysis of GWAS results. In the spirit of comparable tools for gene-expression analysis, we attempt to unify and simplify several procedures that are essential for the interpretation of GWAS results. This includes the generation of advanced Manhattan and regional association plots including rare variant display as well as novel interaction network analysis tools for the investigation of systems-biology aspects. Our package supports virtually all model organisms and represents the first cohesive implementation of such tools for the popular language R. Previous software of that range is dispersed over a wide range of platforms and mostly not adaptable for custom work pipelines. We demonstrate the utility of this package by providing an example workflow on a publicly available dataset.

BACKGROUND: The distribution of human disease-associated mutations is not random across the human genome. Despite the fact that natural selection continually removes disease-associated mutations, an enrichment of these variants can be observed in regions of low recombination. There are a number of mechanisms by which such a clustering could occur, including genetic perturbations or demographic effects within different populations. Recent genome-wide association studies (GWAS) suggest that single nucleotide polymorphisms (SNPs) associated with complex disease traits are not randomly distributed throughout the genome, but tend to cluster in regions of low recombination. PRINCIPAL FINDINGS: Here we investigated whether deleterious mutations have accumulated in regions of low recombination due to the impact of recent positive selection and genetic hitchhiking. Using publicly available data on common complex diseases and population demography, we observed an enrichment of hitchhiked disease associations in conserved gene clusters subject to selection pressure. Evolutionary analysis revealed that these conserved gene clusters arose by multiple concerted rearrangements events across the vertebrate lineage. We observed distinct clustering of disease-associated SNPs in evolutionary rearranged regions of low recombination and high gene density, which harbor genes involved in immunity, that is, the interleukin cluster on 5q31 or RhoA on 3p21. CONCLUSIONS: Our results suggest that multiple lineage specific rearrangements led to a physical clustering of functionally related and linked genes exhibiting an enrichment of susceptibility loci for complex traits. This implies that besides recent evolutionary adaptations other evolutionary dynamics have played a role in the formation of linked gene clusters associated with complex disease traits.

Elevated levels of acute phase proteins (APP) are often found in patients with cardiovascular diseases. In a previous study, we demonstrated the importance of the IL-6-gp130 axis -as a key regulator of inflammatory acute phase signaling in hepatocytes-for the development of atherosclerosis.Gp130-dependent gene expression was analyzed in a previously established hepatocyte-specific gp130 knockout mouse model. We performed whole transcriptome analysis in isolated hepatocytes to measure tissue specific responses after proinflammatory stimulus with IL-6 across different time points. Our analyses revealed an unexpected small gene cluster that requires IL-6 stimulus for early activation. Several of the genes in this cluster are involved in different cell defense mechanisms. Thus, stressors that trigger both general stress and inflammatory responses lead to activation of a stereotypic innate cellular defense response. Furthermore, we identified a potential biomarker Lipocalin (LCN) 2 for the gp130 dependent early inflammatory response.Our findings suggest a complex network of tightly linked genes involved in the early activation of different parts of the innate immune response including acute phase proteins, complement and coagulation cascade.

Crohn disease (CD), a sub-entity of inflammatory bowel disease (IBD), is a complex polygenic disorder. Although recent studies have successfully identified CD-associated genetic variants, these susceptibility loci explain only a fraction of the heritability of the disease. Here, we report on a multi-stage genome-wide scan of 393 German CD cases and 399 controls. Among the 116,161 single-nucleotide polymorphisms tested, an association with the known CD susceptibility gene NOD2, the 5q31 haplotype, and the recently reported CD locus at 5p13.1 was confirmed. In addition, SNP rs1793004 in the gene encoding nel-like 1 precursor (NELL1, chromosome 11p15.1) showed a consistent disease-association in independent German population- and family-based samples (942 cases, 1082 controls, 375 trios). Subsequent fine mapping and replication in an independent sample of 454 French/Canadian CD trios supported the authenticity of the NELL1 association. Further confirmation in a large German ulcerative colitis (UC) sample indicated that NELL1 is a ubiquitous IBD susceptibility locus (combined p