Moving Domain Computational Fluid Dynamics to Interface with an Embryonic Model of Cardiac Morphogenesis

Peristaltic contraction of the embryonic heart tube produces time- and spatial-varying wall shear stress (WSS) and pressure gradients (∇P) across the atrioventricular (AV) canal. Zebrafish (Danio rerio) are a genetically tractable system to investigate cardiac morphogenesis. The use of Tg(fli1a:EGFP)y1 transgenic embryos allowed for delineation and two-dimensional reconstruction of the endocardium. This time-varying wall motion was then prescribed in a two-dimensional moving domain computational fluid dynamics (CFD) model, providing new insights into spatial and temporal variations in WSS and ∇P during cardiac development. The CFD simulations were validated with particle image velocimetry (PIV) across the atrioventricular (AV) canal, revealing an increase in both velocities and heart rates, but a decrease in the duration of atrial systole from early to later stages. At 20-30 hours post fertilization (hpf), simulation results revealed bidirectional WSS across the AV canal in the heart tube in response to peristaltic motion of the wall. At 40-50 hpf, the tube structure undergoes cardiac looping, accompanied by a nearly 3-fold increase in WSS magnitude. At 110-120 hpf, distinct AV valve, atrium, ventricle, and bulbus arteriosus form, accompanied by incremental increases in both WSS magnitude and ∇P, but a decrease in bi-directional flow. Laminar flow develops across the AV canal at 20-30 hpf, and persists at 110-120 hpf. Reynolds numbers at the AV canal increase from 0.07±0.03 at 20-30 hpf to 0.23±0.07 at 110-120 hpf (p< 0.05, n=6), whereas Womersley numbers remain relatively unchanged from 0.11 to 0.13. Our moving domain simulations highlights hemodynamic changes in relation to cardiac morphogenesis; thereby, providing a 2-D quantitative approach to complement imaging analysis. © 2013 Lee et al

An atypical epigenetic mechanism affects uniparental expression of Pol IV-dependent siRNAs.

BACKGROUND: Small RNAs generated by RNA polymerase IV (Pol IV) are the most abundant class of small RNAs in flowering plants. In Arabidopsis thaliana Pol IV-dependent short interfering (p4-si)RNAs are imprinted and accumulate specifically from maternal chromosomes in the developing seeds. Imprinted expression of protein-coding genes is controlled by differential DNA or histone methylation placed in gametes. To identify epigenetic factors required for maternal-specific expression of p4-siRNAs we analyzed the effect of a series of candidate mutations, including those required for genomic imprinting of protein-coding genes, on uniparental expression of a representative p4-siRNA locus. RESULTS: Paternal alleles of imprinted genes are marked by DNA or histone methylation placed by DNA METHYLTRANSFERASE 1 or the Polycomb Repressive Complex 2. Here we demonstrate that repression of paternal p4-siRNA expression at locus 08002 is not controlled by either of these mechanisms. Similarly, loss of several chromatin modification enzymes, including a histone acetyltransferase, a histone methyltransferase, and two nucleosome remodeling proteins, does not affect maternal expression of locus 08002. Maternal alleles of imprinted genes are hypomethylated by DEMETER DNA glycosylase, yet expression of p4-siRNAs occurs irrespective of demethylation by DEMETER or related glycosylases. CONCLUSIONS: Differential DNA methylation and other chromatin modifications associated with epigenetic silencing are not required for maternal-specific expression of p4-siRNAs at locus 08002. These data indicate that there is an as yet unknown epigenetic mechanism causing maternal-specific p4-siRNA expression that is distinct from the well-characterized mechanisms associated with DNA methylation or the Polycomb Repressive Complex 2

Spatial and temporal variations in hemodynamic forces initiate cardiac trabeculation

Hemodynamic shear force has been implicated as modulating Notch signaling-mediated cardiac trabeculation. Whether the spatiotemporal variations in wall shear stress (WSS) coordinate the initiation of trabeculation to influence ventricular contractile function remains unknown. Using light-sheet fluorescent microscopy, we reconstructed the 4D moving domain and applied computational fluid dynamics to quantify 4D WSS along the trabecular ridges and in the groves. In WT zebrafish, pulsatile shear stress developed along the trabecular ridges, with prominent endocardial Notch activity at 3 days after fertilization (dpf), and oscillatory shear stress developed in the trabecular grooves, with epicardial Notch activity at 4 dpf. Genetic manipulations were performed to reduce hematopoiesis and inhibit atrial contraction to lower WSS in synchrony with attenuation of oscillatory shear index (OSI) during ventricular development. γ-Secretase inhibitor of Notch intracellular domain (NICD) abrogated endocardial and epicardial Notch activity. Rescue with NICD mRNA restored Notch activity sequentially from the endocardium to trabecular grooves, which was corroborated by observed Notch-mediated cardiomyocyte proliferations on WT zebrafish trabeculae. We also demonstrated in vitro that a high OSI value correlated with upregulated endothelial Notch-related mRNA expression. In silico computation of energy dissipation further supports the role of trabeculation to preserve ventricular structure and contractile function. Thus, spatiotemporal variations in WSS coordinate trabecular organization for ventricular contractile function

OPENMENDEL: A Cooperative Programming Project for Statistical Genetics

Statistical methods for genomewide association studies (GWAS) continue to improve. However, the increasing volume and variety of genetic and genomic data make computational speed and ease of data manipulation mandatory in future software. In our view, a collaborative effort of statistical geneticists is required to develop open source software targeted to genetic epidemiology. Our attempt to meet this need is called the OPENMENDELproject (https://openmendel.github.io). It aims to (1) enable interactive and reproducible analyses with informative intermediate results, (2) scale to big data analytics, (3) embrace parallel and distributed computing, (4) adapt to rapid hardware evolution, (5) allow cloud computing, (6) allow integration of varied genetic data types, and (7) foster easy communication between clinicians, geneticists, statisticians, and computer scientists. This article reviews and makes recommendations to the genetic epidemiology community in the context of the OPENMENDEL project.Comment: 16 pages, 2 figures, 2 table

Structural basis for potency differences between GDF8 and GDF11.

BACKGROUND: Growth/differentiation factor 8 (GDF8) and GDF11 are two highly similar members of the transforming growth factor β (TGFβ) family. While GDF8 has been recognized as a negative regulator of muscle growth and differentiation, there are conflicting studies on the function of GDF11 and whether GDF11 has beneficial effects on age-related dysfunction. To address whether GDF8 and GDF11 are functionally identical, we compared their signaling and structural properties. RESULTS: Here we show that, despite their high similarity, GDF11 is a more potent activator of SMAD2/3 and signals more effectively through the type I activin-like receptor kinase receptors ALK4/5/7 than GDF8. Resolution of the GDF11:FS288 complex, apo-GDF8, and apo-GDF11 crystal structures reveals unique properties of both ligands, specifically in the type I receptor binding site. Lastly, substitution of GDF11 residues into GDF8 confers enhanced activity to GDF8. CONCLUSIONS: These studies identify distinctive structural features of GDF11 that enhance its potency, relative to GDF8; however, the biological consequences of these differences remain to be determined

Spatial and temporal variations in hemodynamic forces initiate cardiac trabeculation

Hemodynamic shear force has been implicated as modulating Notch signaling-mediated cardiac trabeculation. Whether the spatiotemporal variations in wall shear stress (WSS) coordinate the initiation of trabeculation to influence ventricular contractile function remains unknown. Using light-sheet fluorescent microscopy, we reconstructed the 4D moving domain and applied computational fluid dynamics to quantify 4D WSS along the trabecular ridges and in the groves. In WT zebrafish, pulsatile shear stress developed along the trabecular ridges, with prominent endocardial Notch activity at 3 days after fertilization (dpf), and oscillatory shear stress developed in the trabecular grooves, with epicardial Notch activity at 4 dpf. Genetic manipulations were performed to reduce hematopoiesis and inhibit atrial contraction to lower WSS in synchrony with attenuation of oscillatory shear index (OSI) during ventricular development. γ-Secretase inhibitor of Notch intracellular domain (NICD) abrogated endocardial and epicardial Notch activity. Rescue with NICD mRNA restored Notch activity sequentially from the endocardium to trabecular grooves, which was corroborated by observed Notch-mediated cardiomyocyte proliferations on WT zebrafish trabeculae. We also demonstrated in vitro that a high OSI value correlated with upregulated endothelial Notch-related mRNA expression. In silico computation of energy dissipation further supports the role of trabeculation to preserve ventricular structure and contractile function. Thus, spatiotemporal variations in WSS coordinate trabecular organization for ventricular contractile function

The United States COVID-19 Forecast Hub dataset

Academic researchers, government agencies, industry groups, and individuals have produced forecasts at an unprecedented scale during the COVID-19 pandemic. To leverage these forecasts, the United States Centers for Disease Control and Prevention (CDC) partnered with an academic research lab at the University of Massachusetts Amherst to create the US COVID-19 Forecast Hub. Launched in April 2020, the Forecast Hub is a dataset with point and probabilistic forecasts of incident cases, incident hospitalizations, incident deaths, and cumulative deaths due to COVID-19 at county, state, and national, levels in the United States. Included forecasts represent a variety of modeling approaches, data sources, and assumptions regarding the spread of COVID-19. The goal of this dataset is to establish a standardized and comparable set of short-term forecasts from modeling teams. These data can be used to develop ensemble models, communicate forecasts to the public, create visualizations, compare models, and inform policies regarding COVID-19 mitigation. These open-source data are available via download from GitHub, through an online API, and through R packages

Variance Component Selection With Applications to Microbiome Taxonomic Data.

High-throughput sequencing technology has enabled population-based studies of the role of the human microbiome in disease etiology and exposure response. Microbiome data are summarized as counts or composition of the bacterial taxa at different taxonomic levels. An important problem is to identify the bacterial taxa that are associated with a response. One method is to test the association of specific taxon with phenotypes in a linear mixed effect model, which incorporates phylogenetic information among bacterial communities. Another type of approaches consider all taxa in a joint model and achieves selection via penalization method, which ignores phylogenetic information. In this paper, we consider regression analysis by treating bacterial taxa at different level as multiple random effects. For each taxon, a kernel matrix is calculated based on distance measures in the phylogenetic tree and acts as one variance component in the joint model. Then taxonomic selection is achieved by the lasso (least absolute shrinkage and selection operator) penalty on variance components. Our method integrates biological information into the variable selection problem and greatly improves selection accuracies. Simulation studies demonstrate the superiority of our methods versus existing methods, for example, group-lasso. Finally, we apply our method to a longitudinal microbiome study of Human Immunodeficiency Virus (HIV) infected patients. We implement our method using the high performance computing language Julia. Software and detailed documentation are freely available at https://github.com/JingZhai63/VCselection