• European Marine Science
• Research dataclose
• Other research productsclose
• 2014-2023close
• European Commission close
• Wellcome Trust close
• EC|H2020 close
• EU close
• English close

To test associations between DNA methylation levels and gene expression levels in cis (cis eQTMs), we paired each Gene to CpGs closer than 500 kb from its TSS, either upstream or downstream. For each Gene, the TSS was defined based on HTA-2.0 annotation, using the start position for transcripts in the + strand, and the end position for transcripts in the - strand. CpGs position was obtained from Illumina 450K array annotation. Only CpGs in autosomal chromosomes (from chromosome 1 to 22) were tested. In the main analysis, we fitted for each CpG-Gene pair a linear regression model between gene expression and methylation levels adjusted for age, sex, cohort, and blood cell type composition. A second model was run without adjusting for blood cellular composition and it is only reported on the online web catalog, but not discussed in this manuscript. Although some of the unique associations of the unadjusted model might be real, others might be confounded by the large methylation and expression changes among blood cell types. To ensure that CpGs paired to a higher number of Genes do not have higher chances of being part of an eQTM, multiple-testing was controlled at the CpG level, following a procedure previously applied in the Genotype-Tissue Expression (GTEx) project (Gamazon et al., 2018). Briefly, our statistic used to test the hypothesis that a pair CpG-Gene is significantly associated is based on considering the lowest p-value observed for a given CpG and all its paired Gene (e.g., those in the 1 Mb window centered at the TSS). As we do not know the distribution of this statistic under the null, we used a permutation test. We generated 100 permuted gene expression datasets and ran our previous linear regression models obtaining 100 permuted p-values for each CpG-Gene pair. Then, for each CpG, we selected among all CpG-Gene pairs the minimum p-value in each permutation and fitted a beta distribution that is the distribution we obtain when dealing with extreme values (e.g. minimum) (Dudbridge and Gusnanto, 2008). Next, for each CpG, we took the minimum p-value observed in the real data and used the beta distribution to compute the probability of observing a lower p-value. We defined this probability as the empirical p-value of the CpG. Then, we considered as significant those CpGs with empirical p-values to be significant at 5% false discovery rate using Benjamini-Hochberg method. Finally, we applied a last step to identify all significant CpG-Gene pairs for all eCpGs. To do so, we defined a genome-wide empirical p-value threshold as the empirical p-value of the eCpG closest to the 5% false discovery rate threshold. We used this empirical p-value to calculate a nominal p-value threshold for each eCpG, based on the beta distribution obtained from the minimum permuted p-values. This nominal p-value threshold was defined as the value for which the inverse cumulative distribution of the beta distribution was equal to the empirical p-value. Then, for each eCpG, we considered as significant all eCpG-Gene variants with a p-value smaller than nominal p-value. For the meQTLs catalogue, we selected 9.9 M cis and trans meQTLs with a p-value <1e-7 in the ARIES dataset consisting of data from children of 7 years old (Gaunt et al., 2016). Then, we tested whether this subset of 9.9 M SNPs were also meQTLs in HELIX by running meQTL analyses using MatrixEQTL R package (Shabalin, 2012), adjusting for cohort, sex, age, blood cellular composition and the first 20 principal components (PCs) calculated from genome-wide genetic data of the GWAS variability. We confirmed 2.8 M meQTLs in HELIX (p-value <1e-7). Trans meQTLs represented <10% of the 2.8 M meQTLs. Enrichment of eCpGs for meQTLs was computed using a Chi-square test, using non eCpGs as background. Finally, we tested whether meQTLs were also eQTLs for the eGenes linked to the eCpGs. To this end, we run eQTL analyses (gene expression being the outcome and 2.8 M SNPs the predictors) with MatrixEQTL adjusting for cohort, sex, age, blood cellular composition and the first 20 GWAS PCs in HELIX. We considered as significant eQTLs the SNP-Gene pairs with p-value <1e-7 and with the direction of the effect consistent with the direction of the meQTL and the eQTM. Background: The identification of expression quantitative trait methylation (eQTMs), defined as associations between DNA methylation levels and gene expression, might help the biological interpretation of epigenome-wide association studies (EWAS). We aimed to identify autosomal cis eQTMs in children’s blood, using data from 832 children of the Human Early Life Exposome (HELIX) project. Methods: Blood DNA methylation and gene expression were measured with the Illumina 450K and the Affymetrix HTA v2 arrays, respectively. The relationship between methylation levels and expression of nearby genes (1 Mb window centered at the transcription start site, TSS) was assessed by fitting 13.6 M linear regressions adjusting for sex, age, cohort, and blood cell composition. Results: We identified 39,749 blood autosomal cis eQTMs, representing 21,966 unique CpGs (eCpGs, 5.7% of total CpGs) and 8,886 unique transcript clusters (eGenes, 15.3% of total transcript clusters, equivalent to genes). In 87.9% of these cis eQTMs, the eCpG was located at <250 kb from eGene’s TSS; and 58.8% of all eQTMs showed an inverse relationship between the methylation and expression levels. Only around half of the autosomal cis-eQTMs eGenes could be captured through annotation of the eCpG to the closest gene. eCpGs had less measurement error and were enriched for active blood regulatory regions and for CpGs reported to be associated with environmental exposures or phenotypic traits. 40.4% of eQTMs had at least one genetic variant associated with methylation and expression levels. The overlap of autosomal cis eQTMs in children’s blood with those described in adults was small (13.8%), and age-shared cis eQTMs tended to be proximal to the TSS and enriched for genetic variants. See HELIX_Blood_eQTM_READMEfile_20210205.xlsx.

Many sensory systems use ribbon-type synapses to transmit their signals to downstream circuits. The properties of this synaptic transfer fundamentally dictate which aspects in the original stimulus will be accentuated or suppressed, thereby partially defining the detection limits of the circuit. Accordingly, sensory neurons have evolved a wide variety of ribbon geometries and vesicle pool properties to best support their diverse functional requirements. However, the need for diverse synaptic functions does not only arise across neuron types, but also within. Here we show that UV-cones, a single type of photoreceptor of the larval zebrafish eye, exhibit striking differences in their synaptic ultrastructure and consequent calcium to glutamate transfer function depending on their location in the eye. We arrive at this conclusion by combining serial section electron microscopy and simultaneous "dual-colour" 2-photon imaging of calcium and glutamate signals from the same synapse in vivo. We further use the functional dataset to fit a cascade-like model of the ribbon synapse with different vesicle pool sizes, transfer rates and other synaptic properties. Exploiting recent developments in simulation-based inference, we obtain full posterior estimates for the parameters and compare these across different retinal regions. The model enables us to extrapolate to new stimuli and to systematically investigate different response behaviours of various ribbon configurations. We also provide an interactive, easy-to-use version of this model as an online tool. Overall, we show that already on the synaptic level of single neuron types there exist highly specialized mechanisms which are advantageous for the encoding of different visual features.

Objective: We employed Mendelian Randomization to explore whether the effects of blood pressure (BP) and BP lowering through different antihypertensive drug classes on stroke risk vary by stroke etiology. Methods: We selected genetic variants associated with systolic and diastolic BP and BP-lowering variants in genes encoding antihypertensive drug targets from a GWAS on 757,601 individuals. Applying two-sample Mendelian randomization, we examined associations with any stroke (67,162 cases; 454,450 controls), ischemic stroke and its subtypes (large artery, cardioembolic, small vessel stroke), intracerebral hemorrhage (ICH, deep and lobar), and the related small vessel disease phenotype of WMH. Results: Genetic predisposition to higher systolic and diastolic BP was associated with higher risk of any stroke, ischemic stroke, and ICH. We found associations between genetically determined BP and all ischemic stroke subtypes with a higher risk of large artery and small vessel stroke compared to cardioembolic stroke, as well as associations with deep, but not lobar ICH. Genetic proxies for calcium channel blockers, but not beta blockers, were associated with lower risk of any stroke and ischemic stroke. Proxies for CCBs showed particularly strong associations with small vessel stroke and the related radiological phenotype of WMH. Conclusions: This study supports a causal role of hypertension in all major stroke subtypes except lobar ICH. We find differences in the effects of BP and BP lowering through antihypertensive drug classes between stroke subtypes and identify calcium channel blockade as a promising strategy for preventing manifestations of cerebral small vessel disease. Objective: We employed Mendelian Randomization to explore whether the effects of blood pressure (BP) and BP lowering through different antihypertensive drug classes on stroke risk vary by stroke etiology. Methods: We selected genetic variants associated with systolic and diastolic BP and BP-lowering variants in genes encoding antihypertensive drug targets from a GWAS on 757,601 individuals. Applying two-sample Mendelian randomization, we examined associations with any stroke (67,162 cases; 454,450 controls), ischemic stroke and its subtypes (large artery, cardioembolic, small vessel stroke), intracerebral hemorrhage (ICH, deep and lobar), and the related small vessel disease phenotype of WMH. Results: Genetic predisposition to higher systolic and diastolic BP was associated with higher risk of any stroke, ischemic stroke, and ICH. We found associations between genetically determined BP and all ischemic stroke subtypes with a higher risk of large artery and small vessel stroke compared to cardioembolic stroke, as well as associations with deep, but not lobar ICH. Genetic proxies for calcium channel blockers, but not beta blockers, were associated with lower risk of any stroke and ischemic stroke. Proxies for CCBs showed particularly strong associations with small vessel stroke and the related radiological phenotype of WMH. Conclusions: This study supports a causal role of hypertension in all major stroke subtypes except lobar ICH. We find differences in the effects of BP and BP lowering through antihypertensive drug classes between stroke subtypes and identify calcium channel blockade as a promising strategy for preventing manifestations of cerebral small vessel disease. 1

3D-printing networks of droplets connected by interface bilayers is a powerful platform to build synthetic tissues, in which functionality relies on precisely ordered structures. However, the structural precision and consistency in assembling these structures is currently limited, which restricts intricate designs and the complexity of functions performed by synthetic tissues. Here, we report that the equilibrium contact angle (θDIB) between a pair of droplets is a key parameter that dictates the tessellation and precise positioning of hundreds of picolitre droplets within 3D-printed, multi-layer networks. When θDIB approximates the geometrically-derived critical angle (θc) of 35.3º, the resulting networks of droplets arrange in regular hexagonally close-packed (hcp) lattices with the least fraction of defects. With this improved control over droplet packing, we can 3D-print functional synthetic tissues with single-droplet-wide conductive pathways. Our new insights into 3D droplet packing permit the fabrication of complex synthetic tissues, where precisely positioned compartments perform coordinated tasks. The files are labelled in the format: SILXX_POPCYY_NetNumber_Process.tif Where: "SILXX_POPCYY" indicates the volume fraction of silicone oil (φSIL) and molar fraction of POPC (xPOPC) at which the network was printed (e.g. "SIL55_POPC13" corresponds to φSIL = 0.55 and xPOPC = 0.13) "NetNumber" indicates the repeat number of the specific printed network (e.g. "net3" indicates the 3rd network printed at a specific condition) "Process" indicates the type of processing visualised in the image. This can be: " " : Raw confocal image "Network": Segmented image showing automatic identification of the lipid bilayers in droplet networks "Network_Corrected": Manually corrected image after automatic segmentation "LinkClasses": Segmented image showing classification of lipid bilayers and lipid monolayers "Overlay": Overlay of the confocal image and corresponding packing classification based on Delaunay triangulation.

Axon size plays a crucial role in determining conductance velocity and, consequently, in the the timing and synchronization of neural activation. Noninvasive measurement of axon radii could have significant impact on the understanding of healthy and diseased neural processes. However, until now, accurate axon radius mapping has eluded in vivo neuroimaging, mainly due to a lack of sensitivity of the MRI signal to micron-sized axons. Here, we show how -- when confounding factors such as extra-axonal water and axonal orientation dispersion are eliminated -- heavily diffusion-weighted MRI signals becomes sensitive to axon radii. However, diffusion MRI is only capable of estimating a single metric representing the entire axon radius distribution within a voxel that emphasizes the largest axons. Our findings, both in rodents and humans, enable noninvasive mapping of critical information on axon radii, as well as resolve the long-standing debate on whether axon radii can be quantified. The dataset is comprised of three main items: 1. Human diffusion MRI 2. Rodent diffusion MRI 3. Rodent confocal microscopy Please consider the information and technical details provided in the file "Documentation.pdf" when using the data. The data acquisition is decribed in the manuscript. Please contact the authors for additional information.

The marine silicon cycle is intrinsically linked with carbon cycling in the oceans via biological production of silica by a wide range of organisms. The stable silicon isotopic composition (denoted by δ30Si) of siliceous microfossils extracted from sediment cores can be used as an archive of past oceanic silicon cycling. However, the silicon isotopic composition of biogenic silica has only been measured in diatoms, sponges and radiolarians, and isotopic fractionation relative to seawater is entirely unknown for many other silicifiers. Furthermore, the biochemical pathways and mechanisms that determine isotopic fractionation during biosilicification remain poorly understood. Here, we present the first measurements of the silicon isotopic fractionation during biosilicification by loricate choanoflagellates, a group of protists closely related to animals. We cultured two species of choanoflagellates, Diaphanoeca grandis and Stephanoeca diplocostata, which showed consistently greater isotopic fractionation (approximately −5 ‰ to −7 ‰) than cultured diatoms (−0.5 ‰ to −2.1 ‰). Instead, choanoflagellate silicon isotopic fractionation appears to be more similar to sponges grown under similar dissolved silica concentrations. Our results highlight that there is a taxonomic component to silicon isotope fractionation during biosilicification, possibly via a shared or related biochemical transport pathway. These findings have implications for the use of biogenic silica δ30Si produced by different silicifiers as proxies for past oceanic change.

Perception adapts to mismatching multisensory information, both when different cues appear simultaneously and when they appear sequentially. While both multisensory integration and adaptive trial-by-trial recalibration are central for behavior, it remains unknown whether they are mechanistically linked and arise from a common neural substrate. To relate the neural underpinnings of sensory integration and recalibration, we measured whole-brain magnetoencephalography while human participants performed an audio-visual ventriloquist task. Using single-trial multivariate analysis, we localized the perceptually-relevant encoding of multisensory information within and between trials. While we found neural signatures of multisensory integration within temporal and parietal regions, only medial superior parietal activity encoded past and current sensory information and mediated the perceptual recalibration within and between trials. These results highlight a common neural substrate of sensory integration and perceptual recalibration, and reveal a role of medial parietal regions in linking present and previous multisensory evidence to guide adaptive behavior. PARK_KAYSER_RecalMEG_2019Please see README file.

Dogs were present in the Americas prior to the arrival of European colonists, but the origin and fate of these pre-contact dogs are largely unknown. We sequenced 71 mitochondrial and seven nuclear genomes from ancient North American and Siberian dogs spanning ~9,000 years. Our analysis indicates that American dogs were not domesticated from North American wolves. Instead, American dogs form a monophyletic lineage that likely originated in Siberia and dispersed into the Americas alongside people. After the arrival of Europeans, native American dogs almost completely disappeared, leaving a minimal genetic legacy in modern dog populations. Remarkably, the closest detectable extant lineage to pre-contact American dogs is the canine transmissible venereal tumor, a contagious cancer clone derived from an individual dog that lived up to 8,000 years ago. Mitochondrial DNA FASTA fileFASTA file containing 1166 dog mtDNA genomes used in this studyfull_mtDNA_alignment.fastaNEXUS treeMaximum likelihood tree (RAxML) of 1166 dogs mtDNA genomes used in this studyfull_mtDNA_alignment.treExcel sheetPublication source of the 1166 mtDNA genomes used in this studyfull_mtDNA_alignment.xlsxPlink (bed) fileContains genotype for dogs 54 dogsfull_data.bedPlink file (bim)Contains genotype for 54 dogsfull_data.bimPlink file (fam)Contains genotype for 54 dogsfull_data.famNJ tree in Figure 2bNJ tree in Figure 2b (see Table S2 for more info)Figure_b.treNexus fileNexus file used for producing Figure S12 (MKV model in MrBayes)Binary_char_MKV.nexNEXUS treeBayesian tree in Figure S12 (see Table S2 for more info)Figure_S12.tre