Prediction of causal genes at GWAS loci with pleiotropic gene regulatory effects using sets of correlated instrumental variables.

Multivariate Mendelian randomization (MVMR) is a statistical technique that uses sets of genetic instruments to estimate the direct causal effects of multiple exposures on an outcome of interest. At genomic loci with pleiotropic gene regulatory effects, that is, loci where the same genetic variants are associated to multiple nearby genes, MVMR can potentially be used to predict candidate causal genes. However, consensus in the field dictates that the genetic instruments in MVMR must be independent (not in linkage disequilibrium), which is usually not possible when considering a group of candidate genes from the same locus. Here we used causal inference theory to show that MVMR with correlated instruments satisfies the instrumental set condition. This is a classical result by Brito and Pearl (2002) for structural equation models that guarantees the identifiability of individual causal effects in situations where multiple exposures collectively, but not individually, separate a set of instrumental variables from an outcome variable. Extensive simulations confirmed the validity and usefulness of these theoretical results. Importantly, the causal effect estimates remained unbiased and their variance small even when instruments are highly correlated, while bias introduced by horizontal pleiotropy or LD matrix sampling error was comparable to standard MR. We applied MVMR with correlated instrumental variable sets at genome-wide significant loci for coronary artery disease (CAD) risk using expression Quantitative Trait Loci (eQTL) data from seven vascular and metabolic tissues in the STARNET study. Our method predicts causal genes at twelve loci, each associated with multiple colocated genes in multiple tissues. We confirm causal roles for PHACTR1 and ADAMTS7 in arterial tissues, among others. However, the extensive degree of regulatory pleiotropy across tissues and the limited number of causal variants in each locus still require that MVMR is run on a tissue-by-tissue basis, and testing all gene-tissue pairs with cis-eQTL associations at a given locus in a single model to predict causal gene-tissue combinations remains infeasible. Our results show that within tissues, MVMR with dependent, as opposed to independent, sets of instrumental variables significantly expands the scope for predicting causal genes in disease risk loci with pleiotropic regulatory effects. However, considering risk loci with regulatory pleiotropy that also spans across tissues remains an unsolved problem.

In Vitro Development of Human iPSC-Derived Functional Neuronal Networks on Laser-Fabricated 3D Scaffolds.

Neural progenitor cells generated from human induced pluripotent stem cells (hiPSCs) are the forefront of ″brain-on-chip″ investigations. Viable and functional hiPSC-derived neuronal networks are shaping powerful in vitro models for evaluating the normal and abnormal formation of cortical circuits, understanding the underlying disease mechanisms, and investigating the response to drugs. They therefore represent a desirable instrument for both the scientific community and the pharmacological industry. However, culture conditions required for the full functional maturation of individual neurons and networks are still unidentified. It has been recognized that three-dimensional (3D) culture conditions can better emulate in vivo neuronal tissue development compared to 2D cultures and thus provide a more desirable in vitro approach. In this paper, we present the design and implementation of a 3D scaffold platform that supports and promotes intricate neuronal network development. 3D scaffolds were produced through direct laser writing by two-photon polymerization (2PP), a high-resolution 3D laser microstructuring technology, using the biocompatible and nondegradable photoreactive resin Dental LT Clear (DClear). Neurons developed and interconnected on a 3D environment shaped by vertically stacked scaffold layers. The developed networks could support different cell types. Starting at the day 50 of 3D culture, neuronal progenitor cells could develop into cortical projection neurons (CNPs) of all six layers, different types of inhibitory neurons, and glia. Additionally and in contrast to 2D conditions, 3D scaffolds supported the long-term culturing of neuronal networks over the course of 120 days. Network health and functionality were probed through calcium imaging, which revealed a strong spontaneous neuronal activity that combined individual and collective events. Taken together, our results highlight advanced microstructured 3D scaffolds as a reliable platform for the 3D in vitro modeling of neuronal functions.

Comparison between instrumental variable and mediation-based methods for reconstructing causal gene networks in yeast.

Causal gene networks model the flow of information within a cell. Reconstructing causal networks from omics data is challenging because correlation does not imply causation. When genomics and transcriptomics data from a segregating population are combined, genomic variants can be used to orient the direction of causality between gene expression traits. Instrumental variable methods use a local expression quantitative trait locus (eQTL) as a randomized instrument for a gene's expression level, and assign target genes based on distal eQTL associations. Mediation-based methods additionally require that distal eQTL associations are mediated by the source gene. A detailed comparison between these methods has not yet been conducted, due to the lack of a standardized implementation of different methods, the limited sample size of most multi-omics datasets, and the absence of ground-truth networks for most organisms. Here we used Findr, a software package providing uniform implementations of instrumental variable, mediation, and coexpression-based methods, a recent dataset of 1012 segregants from a cross between two budding yeast strains, and the Yeastract database of known transcriptional interactions to compare causal gene network inference methods. We found that causal inference methods result in a significant overlap with the ground-truth, whereas coexpression did not perform better than random. A subsampling analysis revealed that the performance of mediation saturates at large sample sizes, due to a loss of sensitivity when residual correlations become significant. Instrumental variable methods on the other hand contain false positive predictions, due to genomic linkage between eQTL instruments. Instrumental variable and mediation-based methods also have complementary roles for identifying causal genes underlying transcriptional hotspots. Instrumental variable methods correctly predicted STB5 targets for a hotspot centred on the transcription factor STB5, whereas mediation failed due to Stb5p auto-regulating its own expression. Mediation suggests a new candidate gene, DNM1, for a hotspot on Chr XII, whereas instrumental variable methods could not distinguish between multiple genes located within the hotspot. In conclusion, causal inference from genomics and transcriptomics data is a powerful approach for reconstructing causal gene networks, which could be further improved by the development of methods to control for residual correlations in mediation analyses, and for genomic linkage and pleiotropic effects from transcriptional hotspots in instrumental variable analyses.

Impact of Physical Obstacles on the Structural and Effective Connectivity of in silico Neuronal Circuits.

Scaffolds and patterned substrates are among the most successful strategies to dictate the connectivity between neurons in culture. Here, we used numerical simulations to investigate the capacity of physical obstacles placed on a flat substrate to shape structural connectivity, and in turn collective dynamics and effective connectivity, in biologically-realistic neuronal networks. We considered µ-sized obstacles placed in mm-sized networks. Three main obstacle shapes were explored, namely crosses, circles and triangles of isosceles profile. They occupied either a small area fraction of the substrate or populated it entirely in a periodic manner. From the point of view of structure, all obstacles promoted short length-scale connections, shifted the in- and out-degree distributions toward lower values, and increased the modularity of the networks. The capacity of obstacles to shape distinct structural traits depended on their density and the ratio between axonal length and substrate diameter. For high densities, different features were triggered depending on obstacle shape, with crosses trapping axons in their vicinity and triangles funneling axons along the reverse direction of their tip. From the point of view of dynamics, obstacles reduced the capacity of networks to spontaneously activate, with triangles in turn strongly dictating the direction of activity propagation. Effective connectivity networks, inferred using transfer entropy, exhibited distinct modular traits, indicating that the presence of obstacles facilitated the formation of local effective microcircuits. Our study illustrates the potential of physical constraints to shape structural blueprints and remodel collective activity, and may guide investigations aimed at mimicking organizational traits of biological neuronal circuits.

Effect of salt on the H-bond symmetrization in ice.

The richness of the phase diagram of water reduces drastically at very high pressures where only two molecular phases, proton-disordered ice VII and proton-ordered ice VIII, are known. Both phases transform to the centered hydrogen bond atomic phase ice X above about 60 GPa, i.e., at pressures experienced in the interior of large ice bodies in the universe, such as Saturn and Neptune, where nonmolecular ice is thought to be the most abundant phase of water. In this work, we investigate, by Raman spectroscopy up to megabar pressures and ab initio simulations, how the transformation of ice VII in ice X is affected by the presence of salt inclusions in the ice lattice. Considerable amounts of salt can be included in ice VII structure under pressure via rock-ice interaction at depth and processes occurring during planetary accretion. Our study reveals that the presence of salt hinders proton order and hydrogen bond symmetrization, and pushes ice VII to ice X transformation to higher and higher pressures as the concentration of salt is increased.