4,777 research outputs found
3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries
Recent advances in electron microscopy have enabled the imaging of single
cells in 3D at nanometer length scale resolutions. An uncharted frontier for in
silico biology is the ability to simulate cellular processes using these
observed geometries. Enabling such simulations requires watertight meshing of
electron micrograph images into 3D volume meshes, which can then form the basis
of computer simulations of such processes using numerical techniques such as
the Finite Element Method. In this paper, we describe the use of our recently
rewritten mesh processing software, GAMer 2, to bridge the gap between poorly
conditioned meshes generated from segmented micrographs and boundary marked
tetrahedral meshes which are compatible with simulation. We demonstrate the
application of a workflow using GAMer 2 to a series of electron micrographs of
neuronal dendrite morphology explored at three different length scales and show
that the resulting meshes are suitable for finite element simulations. This
work is an important step towards making physical simulations of biological
processes in realistic geometries routine. Innovations in algorithms to
reconstruct and simulate cellular length scale phenomena based on emerging
structural data will enable realistic physical models and advance discovery at
the interface of geometry and cellular processes. We posit that a new frontier
at the intersection of computational technologies and single cell biology is
now open.Comment: 39 pages, 14 figures. High resolution figures and supplemental movies
available upon reques
Generative models of the human connectome
The human connectome represents a network map of the brain's wiring diagram
and the pattern into which its connections are organized is thought to play an
important role in cognitive function. The generative rules that shape the
topology of the human connectome remain incompletely understood. Earlier work
in model organisms has suggested that wiring rules based on geometric
relationships (distance) can account for many but likely not all topological
features. Here we systematically explore a family of generative models of the
human connectome that yield synthetic networks designed according to different
wiring rules combining geometric and a broad range of topological factors. We
find that a combination of geometric constraints with a homophilic attachment
mechanism can create synthetic networks that closely match many topological
characteristics of individual human connectomes, including features that were
not included in the optimization of the generative model itself. We use these
models to investigate a lifespan dataset and show that, with age, the model
parameters undergo progressive changes, suggesting a rebalancing of the
generative factors underlying the connectome across the lifespan.Comment: 38 pages, 5 figures + 19 supplemental figures, 1 tabl
Thermalization, Error-Correction, and Memory Lifetime for Ising Anyon Systems
We consider two-dimensional lattice models that support Ising anyonic
excitations and are coupled to a thermal bath. We propose a phenomenological
model for the resulting short-time dynamics that includes pair-creation,
hopping, braiding, and fusion of anyons. By explicitly constructing topological
quantum error-correcting codes for this class of system, we use our
thermalization model to estimate the lifetime of the quantum information stored
in the encoded spaces. To decode and correct errors in these codes, we adapt
several existing topological decoders to the non-Abelian setting. We perform
large-scale numerical simulations of these two-dimensional Ising anyon systems
and find that the thresholds of these models range between 13% to 25%. To our
knowledge, these are the first numerical threshold estimates for quantum codes
without explicit additive structure.Comment: 34 pages, 9 figures; v2 matches the journal version and corrects a
misstatement about the detailed balance condition of our Metropolis
simulations. All conclusions from v1 are unaffected by this correctio
Inferring Geodesic Cerebrovascular Graphs: Image Processing, Topological Alignment and Biomarkers Extraction
A vectorial representation of the vascular network that embodies quantitative features - location, direction, scale, and bifurcations - has many potential neuro-vascular applications. Patient-specific models support computer-assisted surgical procedures in neurovascular interventions, while analyses on multiple subjects are essential for group-level studies on which clinical prediction and therapeutic inference ultimately depend. This first motivated the development of a variety of methods to segment the cerebrovascular system. Nonetheless, a number of limitations, ranging from data-driven inhomogeneities, the anatomical intra- and inter-subject variability, the lack of exhaustive ground-truth, the need for operator-dependent processing pipelines, and the highly non-linear vascular domain, still make the automatic inference of the cerebrovascular topology an open problem. In this thesis, brain vessels’ topology is inferred by focusing on their connectedness. With a novel framework, the brain vasculature is recovered from 3D angiographies by solving a connectivity-optimised anisotropic level-set over a voxel-wise tensor field representing the orientation of the underlying vasculature. Assuming vessels joining by minimal paths, a connectivity paradigm is formulated to automatically determine the vascular topology as an over-connected geodesic graph. Ultimately, deep-brain vascular structures are extracted with geodesic minimum spanning trees. The inferred topologies are then aligned with similar ones for labelling and propagating information over a non-linear vectorial domain, where the branching pattern of a set of vessels transcends a subject-specific quantized grid. Using a multi-source embedding of a vascular graph, the pairwise registration of topologies is performed with the state-of-the-art graph matching techniques employed in computer vision. Functional biomarkers are determined over the neurovascular graphs with two complementary approaches. Efficient approximations of blood flow and pressure drop account for autoregulation and compensation mechanisms in the whole network in presence of perturbations, using lumped-parameters analog-equivalents from clinical angiographies. Also, a localised NURBS-based parametrisation of bifurcations is introduced to model fluid-solid interactions by means of hemodynamic simulations using an isogeometric analysis framework, where both geometry and solution profile at the interface share the same homogeneous domain. Experimental results on synthetic and clinical angiographies validated the proposed formulations. Perspectives and future works are discussed for the group-wise alignment of cerebrovascular topologies over a population, towards defining cerebrovascular atlases, and for further topological optimisation strategies and risk prediction models for therapeutic inference. Most of the algorithms presented in this work are available as part of the open-source package VTrails
- …