Neural induction by the node and placode induction by head mesoderm share an initial state resembling neural plate border and ES cells.

Around the time of gastrulation in higher vertebrate embryos, inductive interactions direct cells to form central nervous system (neural plate) or sensory placodes. Grafts of different tissues into the periphery of a chicken embryo elicit different responses: Hensen's node induces a neural plate whereas the head mesoderm induces placodes. How different are these processes? Transcriptome analysis in time course reveals that both processes start by induction of a common set of genes, which later diverge. These genes are remarkably similar to those induced by an extraembryonic tissue, the hypoblast, and are normally expressed in the pregastrulation stage epiblast. Explants of this epiblast grown in the absence of further signals develop as neural plate border derivatives and eventually express lens markers. We designate this state as "preborder"; its transcriptome resembles embryonic stem cells. Finally, using sequential transplantation experiments, we show that the node, head mesoderm, and hypoblast are interchangeable to begin any of these inductions while the final outcome depends on the tissue emitting the later signals.

Cardiac circRNAs arise mainly from constitutive exons rather than alternatively spliced exons.

Circular RNAs (circRNAs) are a relatively new class of RNA molecules, and knowledge about their biogenesis and function is still in its infancy. It was recently shown that alternative splicing underlies the formation of circular RNAs (circRNA) arising from the Titin (TTN) gene. Since the main mechanism by which circRNAs are formed is still unclear, we hypothesized that alternative splicing, and in particular exon skipping, is a major driver of circRNA production. We performed RNA sequencing on human and mouse hearts, mapped alternative splicing events, and overlaid these with expressed circRNAs at exon-level resolution. In addition, we performed RNA sequencing on hearts of Rbm20 KO mice to address how important Rbm20-mediated alternative splicing is in the production of cardiac circRNAs. In human and mouse hearts, we show that cardiac circRNAs are mostly (∼90%) produced from constitutive exons and less (∼10%) from alternatively spliced exons. In Rbm20 KO hearts, we identified 38 differentially expressed circRNAs of which 12 were produced from the Ttn gene. Even though Ttn appeared the most prominent target of Rbm20 for circularization, we also detected Rbm20-dependent circRNAs arising from other genes including Fan1, Stk39, Xdh, Bcl2l13, and Sorbs1 Interestingly, only Ttn circRNAs seemed to arise from Rbm20-mediated skipped exons. In conclusion, cardiac circRNAs are mostly derived from constitutive exons, suggesting that these circRNAs are generated at the expense of their linear counterpart and that circRNA production impacts the accumulation of the linear mRNA.

RBM20 Mutations Induce an Arrhythmogenic Dilated Cardiomyopathy Related to Disturbed Calcium Handling.

BACKGROUND: Mutations in RBM20 (RNA-binding motif protein 20) cause a clinically aggressive form of dilated cardiomyopathy, with an increased risk of malignant ventricular arrhythmias. RBM20 is a splicing factor that targets multiple pivotal cardiac genes, such as Titin (TTN) and CAMK2D (calcium/calmodulin-dependent kinase II delta). Aberrant TTN splicing is thought to be the main determinant of RBM20-induced dilated cardiomyopathy, but is not likely to explain the increased risk of arrhythmias. Here, we investigated the extent to which RBM20 mutation carriers have an increased risk of arrhythmias and explore the underlying molecular mechanism. METHODS: We compared clinical characteristics of RBM20 and TTN mutation carriers and used our previously generated Rbm20 knockout (KO) mice to investigate downstream effects of Rbm20-dependent splicing. Cellular electrophysiology and Ca2+ measurements were performed on isolated cardiomyocytes from Rbm20 KO mice to determine the intracellular consequences of reduced Rbm20 levels. RESULTS: Sustained ventricular arrhythmias were more frequent in human RBM20 mutation carriers than in TTN mutation carriers (44% versus 5%, respectively, P=0.006). Splicing events that affected Ca2+- and ion-handling genes were enriched in Rbm20 KO mice, most notably in the genes CamkIIδ and RyR2. Aberrant splicing of CamkIIδ in Rbm20 KO mice resulted in a remarkable shift of CamkIIδ toward the δ-A isoform that is known to activate the L-type Ca2+ current ( ICa,L). In line with this, we found an increased ICa,L, intracellular Ca2+ overload and increased sarcoplasmic reticulum Ca2+ content in Rbm20 KO myocytes. In addition, not only complete loss of Rbm20, but also heterozygous loss of Rbm20 increased spontaneous sarcoplasmic reticulum Ca2+ releases, which could be attenuated by treatment with the ICa,L antagonist verapamil. CONCLUSIONS: We show that loss of Rbm20 disturbs Ca2+ handling and leads to more proarrhythmic Ca2+ releases from the sarcoplasmic reticulum. Patients that carry a pathogenic RBM20 mutation have more ventricular arrhythmias despite a similar left ventricular function, in comparison with patients with a TTN mutation. Our experimental data suggest that RBM20 mutation carriers may benefit from treatment with an ICa,L blocker to reduce their arrhythmia burden.

RBM20 Regulates Circular RNA Production From the Titin Gene.

RATIONALE: RNA-binding motif protein 20 (RBM20) is essential for normal splicing of many cardiac genes, and loss of RBM20 causes dilated cardiomyopathy. Given its role in splicing, we hypothesized an important role for RBM20 in forming circular RNAs (circRNAs), a novel class of noncoding RNA molecules. OBJECTIVE: To establish the role of RBM20 in the formation of circRNAs in the heart. METHODS AND RESULTS: Here, we performed circRNA profiling on ribosomal depleted RNA from human hearts and identified the expression of thousands of circRNAs, with some of them regulated in disease. Interestingly, we identified 80 circRNAs to be expressed from the titin gene, a gene that is known to undergo highly complex alternative splicing. We show that some of these circRNAs are dynamically regulated in dilated cardiomyopathy but not in hypertrophic cardiomyopathy. We generated RBM20-null mice and show that they completely lack these titin circRNAs. In addition, in a cardiac sample from an RBM20 mutation carrier, titin circRNA production was severely altered. Interestingly, the loss of RBM20 caused only a specific subset of titin circRNAs to be lost. These circRNAs originated from the RBM20-regulated I-band region of the titin transcript. CONCLUSIONS: We show that RBM20 is crucial for the formation of a subset of circRNAs that originate from the I-band of the titin gene. We propose that RBM20, by excluding specific exons from the pre-mRNA, provides the substrate to form this class of RBM20-dependent circRNAs.

LXRα improves myocardial glucose tolerance and reduces cardiac hypertrophy in a mouse model of obesity-induced type 2 diabetes.

AIMS/HYPOTHESIS: Diabetic cardiomyopathy is a myocardial disease triggered by impaired insulin signalling, increased fatty acid uptake and diminished glucose utilisation. Liver X receptors (LXRs) are key transcriptional regulators of metabolic homeostasis. However, their effect in the diabetic heart is largely unknown. METHODS: We cloned murine Lxrα (also known as Nr1h3) behind the α-myosin heavy chain (αMhc; also known as Myh6) promoter to create transgenic (Lxrα-Tg) mice and transgene-negative littermates (wild-type [WT]). A mouse model of type 2 diabetes was induced by a high-fat diet (HFD, 60% energy from fat) over 16 weeks and compared with a low-fat diet (10% energy from fat). A mouse model of type 1 diabetes was induced via streptozotocin injection over 12 weeks. RESULTS: HFD manifested comparable increases in body weight, plasma triacylglycerol and insulin resistance per OGTT in Lxrα-Tg and WT mice. HFD significantly increased left ventricular weight by 21% in WT hearts, but only by 5% in Lxrα-Tg. To elucidate metabolic effects in the heart, microPET (positron emission tomography) imaging revealed that cardiac glucose uptake was increased by 1.4-fold in WT mice on an HFD, but further augmented by 1.7-fold in Lxrα-Tg hearts, in part through 5' adenosine monophosphate-activated protein kinase (AMPK) phosphorylation and restoration of glucose transporter 4 (GLUT4). By contrast, streptozotocin-induced ablation of insulin signalling diminished cardiac glucose uptake levels and caused cardiac dysfunction, indicating that insulin may be important in LXRα-mediated glucose uptake. Chromatin immunoprecipitation assays identified natriuretic peptides, atrial natriuretic peptide (ANP) and B-type natriuretic peptide (BNP), as potential direct targets of cardiac LXRα overexpression. CONCLUSIONS/INTERPRETATION: Cardiac-specific LXRα overexpression ameliorates the progression of HFD-induced left ventricular hypertrophy in association with increased glucose reliance and natriuretic peptide signalling during the early phase of diabetic cardiomyopathy. These findings implicate a potential protective role for LXR in targeting metabolic disturbances underlying diabetes.

Computational tools and resources for prediction and analysis of gene regulatory regions in the chick genome.

The discovery of cis-regulatory elements is a challenging problem in bioinformatics, owing to distal locations and context-specific roles of these elements in controlling gene regulation. Here we review the current bioinformatics methodologies and resources available for systematic discovery of cis-acting regulatory elements and conserved transcription factor binding sites in the chick genome. In addition, we propose and make available, a novel workflow using computational tools that integrate CTCF analysis to predict putative insulator elements, enhancer prediction, and TFBS analysis. To demonstrate the usefulness of this computational workflow, we then use it to analyze the locus of the gene Sox2 whose developmental expression is known to be controlled by a complex array of cis-acting regulatory elements. The workflow accurately predicts most of the experimentally verified elements along with some that have not yet been discovered. A web version of the CTCF tool, together with instructions for using the workflow can be accessed from http://toolshed.g2.bx.psu.edu/view/mkhan1980/ctcf_analysis. For local installation of the tool, relevant Perl scripts and instructions are provided in the directory named "code" in the supplementary materials.

A gene regulatory network for neural induction.

During early vertebrate development, signals from a special region of the embryo, the organizer, can redirect the fate of non-neural ectoderm cells to form a complete, patterned nervous system. This is called neural induction and has generally been imagined as a single signalling event, causing a switch of fate. Here, we undertake a comprehensive analysis, in very fine time course, of the events following exposure of competent ectoderm of the chick to the organizer (the tip of the primitive streak, Hensen's node). Using transcriptomics and epigenomics we generate a gene regulatory network comprising 175 transcriptional regulators and 5614 predicted interactions between them, with fine temporal dynamics from initial exposure to the signals to expression of mature neural plate markers. Using in situ hybridization, single-cell RNA-sequencing, and reporter assays, we show that the gene regulatory hierarchy of responses to a grafted organizer closely resembles the events of normal neural plate development. The study is accompanied by an extensive resource, including information about conservation of the predicted enhancers in other vertebrates.

Blood Pressure-Associated Genetic Variants in the Natriuretic Peptide Receptor 1 Gene Modulate Guanylate Cyclase Activity.

BACKGROUND: Human genetic variation in the NPR1 (natriuretic peptide receptor 1 gene, encoding NPR-A, atrial natriuretic peptide receptor 1) was recently shown to affect blood pressure (BP). NPR-A catalyzes the intracellular conversion of guanosine triphosphate to cGMP (cyclic 3',5'-guanosine monophosphate) on binding of ANP, BNP (atrial or brain natriuretic peptide). Increased levels of cGMP decrease BP by inducing natriuresis, diuresis, and vasodilation. METHODS: We performed a meta-analysis of low-frequency and rare NPR1 variants for BP association in up to 491 584 unrelated individuals. To examine whether the identified BP-associated variants affect NPR-A function, the cGMP response to ANP and BNP was measured in cells expressing wild-type NPR1 and cells expressing the NPR1 variants. RESULTS: In this study, we identified BP associations of 3 amino acid altering variants of NPR1. The minor alleles of rs35479618 (p.E967K, gnomAD non-Finnish European allele frequency 0.017) and rs116245325 (p.L1034F, allele frequency 0.0007) were associated with higher BP (P=4.0×10-25 and P=9.9×10-8, respectively), while the minor allele of rs61757359 (p.G541S, allele frequency 0.003) was associated with lower BP (P=1.8×10-9). Cells transiently expressing 967K or 1034F NPR-A displayed decreased cGMP production in response to ANP and BNP (all P<10-6), while cells expressing 541S NPR-A produced more cGMP compared with cells expressing wild-type NPR-A (P≤4.13×10-5 for ANP and P≤4.24×10-3 for BNP). CONCLUSIONS: In summary, the loss or gain of guanylate cyclase activity for these NPR1 allelic variants could explain the higher or lower BP observed for carriers in large population-based studies.

A novel murine model for arrhythmogenic cardiomyopathy points to a pathogenic role of Wnt signalling and miRNA dysregulation.

AIMS: Arrhythmogenic cardiomyopathy (AC) is one of the most common inherited cardiomyopathies, characterized by progressive fibro-fatty replacement in the myocardium. Clinically, AC manifests itself with ventricular arrhythmias, syncope, and sudden death and shows wide inter- and intra-familial variability. Among the causative genes identified so far, those encoding for the desmosomal proteins plakophilin-2 (PKP2), desmoplakin (DSP), and desmoglein-2 (DSG2) are the most commonly mutated. So far, little is known about the molecular mechanism(s) behind such a varied spectrum of phenotypes, although it has been shown that the causative mutations not only lead to structural abnormalities but also affect the miRNA profiling of cardiac tissue. Here, we aimed at studying the pathogenic effects of a nonsense mutation of the desmoglein-2 gene, both at the structural level and in terms of miRNA expression pattern. METHODS AND RESULTS: We generated transgenic mice with cardiomyocyte-specific overexpression of a FLAG-tagged human desmoglein-2 harbouring the Q558* nonsense mutation found in an AC patient. The hearts of these mice showed signs of fibrosis, decrease in desmosomal size and number, and reduction of the Wnt/ß-catenin signalling. Genome-wide RNA-Seq performed in Tg-hQ hearts and non-transgenic hearts revealed that 24 miRNAs were dysregulated in transgenic animals. Further bioinformatic analyses for selected miRNAs suggested that miR-217-5p, miR-499-5p, and miR-708-5p might be involved in the pathogenesis of the disease. CONCLUSION: Down-regulation of the canonical Wnt/ß-catenin signalling might be considered a common key event in the AC pathogenesis. We identified the miRNA signature in AC hearts, with miR-708-5p and miR-217-5p being the most up-regulated and miR-499-5p the most down-regulated miRNAs. All of them were predicted to be involved in the regulation of the Wnt/ß-catenin pathway and might reveal the potential pathophysiology mechanisms of AC, as well as be useful as therapeutic targets for the disease.

Serially measured circulating microRNAs and adverse clinical outcomes in patients with acute heart failure.

AIMS: Previous studies have identified candidate circulating microRNAs (circmiRs) as biomarkers for heart failure (HF) using relatively insensitive arrays, validated in small cohorts. The present study used RNA sequencing to identify novel candidate circmiRs and compared these with previously identified circmiRs in a large, prospective cohort of patients with acute HF (AHF). METHODS AND RESULTS: RNA sequencing of plasma from instrumented pigs was used to identify circmiRs produced by myocardium. Production of known myomiRs and microRNA (miR)-1306-5p was identified. The prognostic values of this and 11 other circmiRs were tested in a prospective cohort of 496 AHF patients, from whom blood samples were collected at up to seven time-points during the study's 1-year follow-up. The primary endpoint was the composite of all-cause mortality and HF rehospitalization. In the prospective AHF cohort, 188 patients reached the primary endpoint, and higher values of repeatedly measured miR-1306-5p were positively associated with risk for reaching the primary endpoint at the same time-point [hazard ratio (HR) 4.69, 95% confidence interval (CI) 2.18-10.06], independent of clinical characteristics and NT-proBNP. Baseline miR-1306-5p did not improve model discrimination/reclassification significantly compared with NT-proBNP. For miR-320a, miR-378a-3p, miR-423-5p and miR-1254, associations with the primary endpoint were present after adjustment for age and sex (HR 1.38, 95% CI 1.12-1.70; HR 1.35, 95% CI 1.04-1.74; HR 1.45, 95% CI 1.10-1.92; HR 1.22, 95% CI 1.00-1.50, respectively). Rates of detection of myomiRs miR-208a-3p and miR-499a-5p were very low. CONCLUSIONS: Repeatedly measured miR-1306-5p was positively associated with adverse clinical outcome in AHF, even after multivariable adjustment including NT-proBNP. However, baseline miR-1306-5p did not add significant discriminatory value to NT-proBNP. Low-abundance, heart-enriched myomiRs are often undetectable, which mandates the development of more sensitive assays.

Biomarker Profiles in Heart Failure Patients With Preserved and Reduced Ejection Fraction.

BACKGROUND: Biomarkers may help us to unravel differences in the underlying pathophysiology between heart failure (HF) patients with a reduced ejection fraction (HFrEF) and a preserved ejection fraction (HFpEF). Therefore, we compared biomarker profiles to characterize pathophysiological differences between patients with HFrEF and HFpEF. METHODS AND RESULTS: We retrospectively analyzed 33 biomarkers from different pathophysiological domains (inflammation, oxidative stress, remodeling, cardiac stretch, angiogenesis, arteriosclerosis, and renal function) in 460 HF patients (21% HFpEF, left ventricular ejection fraction ≥45%) measured at discharge after hospitalization for acute HF. The association between these markers and the occurrence of all-cause mortality and/or HF-related rehospitalizations at 18 months was compared between patients with HFrEF and HFpEF. Patients were 70.6±11.4 years old and 37.4% were female. Patients with HFpEF were older, more often female, and had a higher systolic blood pressure. Levels of high-sensitive C-reactive protein were significantly higher in HFpEF, while levels of pro-atrial-type natriuretic peptide and N-terminal pro-brain natriuretic peptide were higher in HFrEF. Linear regression followed by network analyses revealed prominent inflammation and angiogenesis-associated interactions in HFpEF and mainly cardiac stretch-associated interactions in HFrEF. The angiogenesis-specific marker, neuropilin and the remodeling-specific marker, osteopontin were predictive for all-cause mortality and/or HF-related rehospitalizations at 18 months in HFpEF, but not in HFrEF (P for interaction <0.05). CONCLUSIONS: In HFpEF, inflammation and angiogenesis-mediated interactions are predominantly observed, while stretch-mediated interactions are found in HFrEF. The remodeling marker osteopontin and the angiogenesis marker neuropilin predicted outcome in HFpEF, but not in HFrEF.

Biomarker Profiles of Acute Heart Failure Patients With a Mid-Range Ejection Fraction.

OBJECTIVES: In this study, the authors used biomarker profiles to characterize differences between patients with acute heart failure with a midrange ejection fraction (HFmrEF) and compare them with patients with a reduced (heart failure with a reduced ejection fraction [HFrEF]) and preserved (heart failure with a preserved ejection fraction [HFpEF]) ejection fraction. BACKGROUND: Limited data are available on biomarker profiles in acute HFmrEF. METHODS: A panel of 37 biomarkers from different pathophysiological domains (e.g., myocardial stretch, inflammation, angiogenesis, oxidative stress, hematopoiesis) were measured at admission and after 24 h in 843 acute heart failure patients from the PROTECT trial. HFpEF was defined as left ventricular ejection fraction (LVEF) of ≥50% (n = 108), HFrEF as LVEF of <40% (n = 607), and HFmrEF as LVEF of 40% to 49% (n = 128). RESULTS: Hemoglobin and brain natriuretic peptide levels (300 pg/ml [HFpEF]; 397 pg/ml [HFmrEF]; 521 pg/ml [HFrEF]; ptrend <0.001) showed an upward trend with decreasing LVEF. Network analysis showed that in HFrEF interactions between biomarkers were mostly related to cardiac stretch, whereas in HFpEF, biomarker interactions were mostly related to inflammation. In HFmrEF, biomarker interactions were both related to inflammation and cardiac stretch. In HFpEF and HFmrEF (but not in HFrEF), remodeling markers at admission and changes in levels of inflammatory markers across the first 24 h were predictive for all-cause mortality and rehospitalization at 60 days (pinteraction <0.05). CONCLUSIONS: Biomarker profiles in patients with acute HFrEF were mainly related to cardiac stretch and in HFpEF related to inflammation. Patients with HFmrEF showed an intermediate biomarker profile with biomarker interactions between both cardiac stretch and inflammation markers. (PROTECT-1: A Study of the Selective A1 Adenosine Receptor Antagonist KW-3902 for Patients Hospitalized With Acute HF and Volume Overload to Assess Treatment Effect on Congestion and Renal Function; NCT00328692).

A network analysis to compare biomarker profiles in patients with and without diabetes mellitus in acute heart failure.

AIMS: It is unclear whether distinct pathophysiological processes are present among patients with acute heart failure (AHF), with and without diabetes. Network analysis of biomarkers may identify correlative associations that reflect different pathophysiological pathways. METHODS AND RESULTS: We analysed a panel of 48 circulating biomarkers measured within 24 h of admission for AHF in a subset of patients enrolled in the PROTECT trial. In patients with and without diabetes, we performed a network analysis to identify correlations between measured biomarkers. Compared with patients without diabetes (n = 1111), those with diabetes (n = 922) had a higher prevalence of ischaemic heart disease and traditional coronary risk factors. After multivariable adjustment, patients with and without diabetes had significantly different levels of biomarkers across a spectrum of pathophysiological domains, including inflammation (TNFR-1a, periostin), cardiomyocyte stretch (BNP), angiogenesis (VEGFR, angiogenin), and renal function (NGAL, KIM-1) (adjusted P-value <0.05). Among patients with diabetes, network analysis revealed that periostin strongly clustered with C-reactive protein and interleukin-6. Furthermore, renal markers (creatinine and NGAL) closely associated with potassium and glucose. These findings were not seen among patients without diabetes. CONCLUSION: Patients with AHF and diabetes, compared with those without diabetes, have distinct biomarker profiles. Network analysis suggests that cardiac remodelling, inflammation, and fibrosis are closely associated with each other in patients with diabetes. Furthermore, potassium levels may be sensitive to changes in renal function as reflected by the strong renal-potassium-glucose correlation. These findings were not seen among patients without diabetes and may suggest distinct pathophysiological processes among AHF patients with diabetes.

A strategy to discover new organizers identifies a putative heart organizer.

Organizers are regions of the embryo that can both induce new fates and impart pattern on other regions. So far, surprisingly few organizers have been discovered, considering the number of patterned tissue types generated during development. This may be because their discovery has relied on transplantation and ablation experiments. Here we describe a new approach, using chick embryos, to discover organizers based on a common gene expression signature, and use it to uncover the anterior intestinal portal (AIP) endoderm as a putative heart organizer. We show that the AIP can induce cardiac identity from non-cardiac mesoderm and that it can pattern this by specifying ventricular and suppressing atrial regional identity. We also uncover some of the signals responsible. The method holds promise as a tool to discover other novel organizers acting during development.

Use of biomarkers to establish potential role and function of circulating microRNAs in acute heart failure.

BACKGROUND: Circulating microRNAs (miRNAs) emerge as potential heart failure biomarkers. We aimed to identify associations between acute heart failure (AHF)-specific circulating miRNAs and well-known heart failure biomarkers. METHODS: Associations between 16 biomarkers predictive for 180day mortality and the levels of 12 AHF-specific miRNAs were determined in 100 hospitalized AHF patients, at baseline and 48hours. Patients were divided in 4 pre-defined groups, based on clinical parameters during hospitalization. Correlation analyses between miRNAs and biomarkers were performed and complemented by miRNA target prediction and pathway analysis. RESULTS: No significant correlations were found at hospital admission. However, after 48hours, 7 miRNAs were significantly negatively correlated to biomarkers indicative for a worse clinical outcome in the patient group with the most unfavorable in-hospital course (n=21); miR-16-5p was correlated to C-reactive protein (R=-0.66, p-value=0.0027), miR-106a-5p to creatinine (R=-0.68, p-value=0.002), miR-223-3p to growth differentiation factor 15 (R=-0.69, p-value=0.0015), miR-652-3p to soluble ST-2 (R=-0.77, p-value<0.001), miR-199a-3p to procalcitonin (R=-0.72, p-value<0.001) and galectin-3 (R=-0.73, p-value<0.001) and miR-18a-5p to procalcitonin (R=-0.68, p-value=0.002). MiRNA target prediction and pathway analysis identified several pathways related to cardiac diseases, which could be linked to some of the miRNA-biomarker correlations. CONCLUSIONS: The majority of correlations between circulating AHF-specific miRNAs were related to biomarkers predictive for a worse clinical outcome in a subgroup of worsening heart failure patients at 48hours of hospitalization. The selective findings suggest a time-dependent effect of circulating miRNAs and highlight the susceptibility to individual patient characteristics influencing potential relations between miRNAs and biomarkers.

Low levels of vitamin D are associated with multimorbidity: results from the LifeLines Cohort Study.

BACKGROUND: The prevalence of multimorbidity (≥ 1 disease within an individual) is rapidly increasing. So far, studies on the relationship between vitamin D and morbidity are mainly focusing on effects on single disease domains only, while vitamin D biology is associated with several diseases throughout the human body. METHODS: We studied 8,726 participants from the LifeLines Cohort Study (a cross-sectional, population-based cohort study) and used the self-developed composite morbidity score to study the association between vitamin D levels and multimorbidity. RESULTS: Study participants (mean age 45 ± 13 years, 73% females) had a mean plasma vitamin D level of 59 ± 22 nmol/L. In participants aged between 50 and 60 years, 58% had ≥ 2 affected disease domains, while morbidity score increased with age (70-80 years: 82% morbidity score > 1; > 80 years: 89% morbidity score > 1). Each incremental reduction by 1 standard deviation (SD) of vitamin D level was associated with an 8% higher morbidity score (full model OR 0.92, 95% CI 0.88-0.97, P = 0.001). Participants with vitamin D levels < 25 nmol/L were at highest risk for increasing morbidity prevalence (versus > 80 nmol/L, OR 1.34, 95% CI 1.07-1.67, P = 0.01). CONCLUSIONS: Low levels of vitamin D are associated with higher prevalence of multimorbidity, especially in participants with vitamin D levels < 25 nmol/L. Collectively, our results favor a general, rather than an organ-specific, approach when assessing the impact of vitamin D deficiency.

The transcription factor Pitx2 positions the embryonic axis and regulates twinning.

Embryonic polarity of invertebrates, amphibians and fish is specified largely by maternal determinants, which fixes cell fates early in development. In contrast, amniote embryos remain plastic and can form multiple individuals until gastrulation. How is their polarity determined? In the chick embryo, the earliest known factor is cVg1 (homologous to mammalian growth differentiation factor 1, GDF1), a transforming growth factor beta (TGFß) signal expressed posteriorly before gastrulation. A molecular screen to find upstream regulators of cVg1 in normal embryos and in embryos manipulated to form twins now uncovers the transcription factor Pitx2 as a candidate. We show that Pitx2 is essential for axis formation, and that it acts as a direct regulator of cVg1 expression by binding to enhancers within neighbouring genes. Pitx2, Vg1/GDF1 and Nodal are also key actors in left-right asymmetry, suggesting that the same ancient polarity determination mechanism has been co-opted to different functions during evolution.