Prototyping an Ontological Framework for Cellular Senescence Mechanisms: A Homeostasis Imbalance Perspective

Although cellular senescence is a key factor in organismal aging, with both positive and negative effects on individuals, its mechanisms remain largely unknown. Thus, integrating knowledge is essential to explain how cellular senescence manifests in tissue damage and age-related diseases. Here, we propose an ontological model that organizes knowledge of cellular senescence in a computer-readable form. We manually annotated and defined cellular senescence processes, molecules, anatomical structures, phenotypes, and other entities based on the Homeostasis Imbalance Process ontology (HOIP). We described the mechanisms as causal relationships of processes and modelled a homeostatic imbalance between stress and stress response in cellular senescence for a unified framework. HOIP was assessed formally, and the relationships between cellular senescence and diseases were inferred for higher-order knowledge processing. We visualized cellular senescence processes to support knowledge utilization. Our study provides a knowledge base to help elucidate mechanisms linking cellular and organismal aging.

Unifying community-wide whole-brain imaging datasets enables robust automated neuron identification and reveals determinants of neuron positioning in C. elegans

We develop a data harmonization approach for C. elegans volumetric microscopy data, still or video, consisting of a standardized format, data pre-processing techniques, and a set of human-in-the-loop machine learning based analysis software tools. We unify a diverse collection of 118 whole-brain neural activity imaging datasets from 5 labs, storing these and accompanying tools in an online repository called WormID ( wormid.org ). We use this repository to generate a statistical atlas that, for the first time, enables accurate automated cellular identification that generalizes across labs, approaching human performance in some cases. We mine this repository to identify factors that influence the developmental positioning of neurons. To facilitate communal use of this repository, we created open-source software, code, web-based tools, and tutorials to explore and curate datasets for contribution to the scientific community. This repository provides a growing resource for experimentalists, theorists, and toolmakers to investigate neuroanatomical organization and neural activity across diverse experimental paradigms, develop and benchmark algorithms for automated neuron detection, segmentation, cell identification, tracking, and activity extraction, and inform models of neurobiological development and function.

Harmonizing the Generation and Pre-publication Stewardship of FAIR Image data

Together with the molecular knowledge of genes and proteins, biological images promise to significantly enhance the scientific understanding of complex cellular systems and to advance predictive and personalized therapeutic products for human health. For this potential to be realized, quality-assured image data must be shared among labs at a global scale to be compared, pooled, and reanalyzed, thus unleashing untold potential beyond the original purpose for which the data was generated. There are two broad sets of requirements to enable image data sharing in the life sciences. One set of requirements is articulated in the companion White Paper entitled "Enabling Global Image Data Sharing in the Life Sciences," which is published in parallel and addresses the need to build the cyberinfrastructure for sharing the digital array data (arXiv:2401.13023 [q-bio.OT], https://doi.org/10.48550/arXiv.2401.13023). In this White Paper, we detail a broad set of requirements, which involves collecting, managing, presenting, and propagating contextual information essential to assess the quality, understand the content, interpret the scientific implications, and reuse image data in the context of the experimental details. We start by providing an overview of the main lessons learned to date through international community activities, which have recently made considerable progress toward generating community standard practices for imaging Quality Control (QC) and metadata. We then provide a clear set of recommendations for amplifying this work. The driving goal is to address remaining challenges, and democratize access to common practices and tools for a spectrum of biomedical researchers, regardless of their expertise, access to resources, and geographical location.

Single-molecule tracking of Nanog and Oct4 in living mouse embryonic stem cells uncovers a feedback mechanism of pluripotency maintenance

Nanog and Oct4 are core transcription factors that form part of a gene regulatory network to regulate hundreds of target genes for pluripotency maintenance in mouse embryonic stem cells (ESCs). To understand their function in the pluripotency maintenance, we visualised and quantified the dynamics of single molecules of Nanog and Oct4 in a mouse ESCs during pluripotency loss. Interestingly, Nanog interacted longer with its target loci upon reduced expression or at the onset of differentiation, suggesting a feedback mechanism to maintain the pluripotent state. The expression level and interaction time of Nanog and Oct4 correlate with their fluctuation and interaction frequency, respectively, which in turn depend on the ESC differentiation status. The DNA viscoelasticity near the Oct4 target locus remained flexible during differentiation, supporting its role either in chromatin opening or a preferred binding to uncondensed chromatin regions. Based on these results, we propose a new negative feedback mechanism for pluripotency maintenance via the DNA condensation state-dependent interplay of Nanog and Oct4.

OME-Zarr: a cloud-optimized bioimaging file format with international community support

A growing community is constructing a next-generation file format (NGFF) for bioimaging to overcome problems of scalability and heterogeneity. Organized by the Open Microscopy Environment (OME), individuals and institutes across diverse modalities facing these problems have designed a format specification process (OME-NGFF) to address these needs. This paper brings together a wide range of those community members to describe the cloud-optimized format itself-OME-Zarr-along with tools and data resources available today to increase FAIR access and remove barriers in the scientific process. The current momentum offers an opportunity to unify a key component of the bioimaging domain-the file format that underlies so many personal, institutional, and global data management and analysis tasks.

Centromere/kinetochore is assembled through CENP-C oligomerization

Kinetochore is an essential protein complex required for accurate chromosome segregation. The constitutive centromere-associated network (CCAN), a subcomplex of the kinetochore, associates with centromeric chromatin and provides a platform for the kinetochore assembly. The CCAN protein CENP-C is thought to be a central hub for the centromere/kinetochore organization. However, the role of CENP-C in CCAN assembly needs to be elucidated. Here, we demonstrate that both the CCAN-binding domain and the C-terminal region that includes the Cupin domain of CENP-C are necessary and sufficient for chicken CENP-C function. Structural and biochemical analyses reveal self-oligomerization of the Cupin domains of chicken and human CENP-C. We find that the CENP-C Cupin domain oligomerization is vital for CENP-C function, centromeric localization of CCAN, and centromeric chromatin organization. These results suggest that CENP-C facilitates the centromere/kinetochore assembly through its oligomerization.

OME-Zarr: a cloud-optimized bioimaging file format with international community support

A growing community is constructing a next-generation file format (NGFF) for bioimaging to overcome problems of scalability and heterogeneity. Organized by the Open Microscopy Environment (OME), individuals and institutes across diverse modalities facing these problems have designed a format specification process (OME-NGFF) to address these needs. This paper brings together a wide range of those community members to describe the cloud-optimized format itself -- OME-Zarr -- along with tools and data resources available today to increase FAIR access and remove barriers in the scientific process. The current momentum offers an opportunity to unify a key component of the bioimaging domain -- the file format that underlies so many personal, institutional, and global data management and analysis tasks.

Condensed but liquid-like domain organization of active chromatin regions in living human cells

In eukaryotes, higher-order chromatin organization is spatiotemporally regulated as domains, for various cellular functions. However, their physical nature in living cells remains unclear (e.g., condensed domains or extended fiber loops; liquid-like or solid-like). Using novel approaches combining genomics, single-nucleosome imaging, and computational modeling, we investigated the physical organization and behavior of early DNA replicated regions in human cells, which correspond to Hi-C contact domains with active chromatin marks. Motion correlation analysis of two neighbor nucleosomes shows that nucleosomes form physically condensed domains with ~150-nm diameters, even in active chromatin regions. The mean-square displacement analysis between two neighbor nucleosomes demonstrates that nucleosomes behave like a liquid in the condensed domain on the ~150 nm/~0.5 s spatiotemporal scale, which facilitates chromatin accessibility. Beyond the micrometers/minutes scale, chromatin seems solid-like, which may contribute to maintaining genome integrity. Our study reveals the viscoelastic principle of the chromatin polymer; chromatin is locally dynamic and reactive but globally stable.

Systematic analysis of cell morphodynamics in C. elegans early embryogenesis

The invariant cell lineage of Caenorhabditis elegans allows unambiguous assignment of the identity for each cell, which offers a unique opportunity to study developmental dynamics such as the timing of cell division, dynamics of gene expression, and cell fate decisions at single-cell resolution. However, little is known about cell morphodynamics, including the extent to which they are variable between individuals, mainly due to the lack of sufficient amount and quality of quantified data. In this study, we systematically quantified the cell morphodynamics in 52 C. elegans embryos from the two-cell stage to mid-gastrulation at the high spatiotemporal resolution, 0.5 μm thickness of optical sections, and 30-second intervals of recordings. Our data allowed systematic analyses of the morphological features. We analyzed sphericity dynamics and found a significant increase at the end of metaphase in every cell, indicating the universality of the mitotic cell rounding. Concomitant with the rounding, the volume also increased in most but not all cells, suggesting less universality of the mitotic swelling. Combining all features showed that cell morphodynamics was unique for each cell type. The cells before the onset of gastrulation could be distinguished from all the other cell types. Quantification of reproducibility in cell-cell contact revealed that variability in division timings and cell arrangements produced variability in contacts between the embryos. However, the area of such contacts occupied less than 5% of the total area, suggesting the high reproducibility of spatial occupancies and adjacency relationships of the cells. By comparing the morphodynamics of identical cells between the embryos, we observed diversity in the variability between cells and found it was determined by multiple factors, including cell lineage, cell generation, and cell-cell contact. We compared the variabilities of cell morphodynamics and cell-cell contacts with those in ascidian Phallusia mammillata embryos. The variabilities were larger in C. elegans, despite smaller differences in embryo size and number of cells at each developmental stage.

PHi-C2: interpreting Hi-C data as the dynamic 3D genome state

Summary

High-throughput chromosome conformation capture (Hi-C) is a widely used assay for studying the three-dimensional (3D) genome organization across the whole genome. Here, we present PHi-C2, a Python package supported by mathematical and biophysical polymer modeling that converts input Hi-C matrix data into the polymer model’s dynamics, structural conformations and rheological features. The updated optimization algorithm for regenerating a highly similar Hi-C matrix provides a fast and accurate optimal solution compared to the previous version by eliminating the factors underlying the inefficiency of the optimization algorithm in the iterative optimization process. In addition, we have enabled a Google Colab workflow to run the algorithm, wherein users can easily change the parameters and check the results in the notebook. Overall, PHi-C2 represents a valuable tool for mining the dynamic 3D genome state embedded in Hi-C data.

Availability and implementation

PHi-C2 as the phic Python package is freely available under the GPL license and can be installed from the Python package index. The source code is available from GitHub at https://github.com/soyashinkai/PHi-C2. Moreover, users do not have to prepare a Python environment because PHi-C2 can run on Google Colab (https://bit.ly/3rlptGI).

Supplementary information

Supplementary data are available at Bioinformatics online.

Single-nucleosome imaging reveals steady-state motion of interphase chromatin in living human cells

Dynamic chromatin behavior plays a critical role in various genome functions. However, it remains unclear how chromatin behavior changes during interphase, where the nucleus enlarges and genomic DNA doubles. While the previously reported chromatin movements varied during interphase when measured using a minute or longer time scale, we unveil that local chromatin motion captured by single-nucleosome imaging/tracking on a second time scale remained steady throughout G₁, S, and G₂ phases in live human cells. This motion mode appeared to change beyond this time scale. A defined genomic region also behaved similarly. Combined with Brownian dynamics modeling, our results suggest that this steady-state chromatin motion was mainly driven by thermal fluctuations. Steady-state motion temporarily increased following a DNA damage response. Our findings support the viscoelastic properties of chromatin. We propose that the observed steady-state chromatin motion allows cells to conduct housekeeping functions, such as transcription and DNA replication, under similar environments during interphase.

The PAF1 complex cell‐autonomously promotes oogenesis in Caenorhabditis elegans

The RNA polymerase II‐associated factor 1 complex (PAF1C) is a protein complex that consists of LEO1, RTF1, PAF1, CDC73, and CTR9, and has been shown to be involved in RNA polymerase II‐mediated transcriptional and chromatin regulation. Although it has been shown to regulate a variety of biological processes, the precise role of the PAF1C during germ line development has not been clarified. In this study, we found that reduction in the function of the PAF1C components, LEO‐1, RTFO‐1, PAFO‐1, CDC‐73, and CTR‐9, in Caenorhabditis elegans affects oogenesis. Defects in oogenesis were also confirmed using an oocyte maturation marker, OMA‐1::GFP. While four to five OMA‐1::GFP‐positive oocytes were observed in wild‐type animals, their numbers were significantly decreased in pafo‐1 mutant and leo‐1(RNAi), pafo‐1(RNAi), and cdc‐73(RNAi) animals. Expression of a functional PAFO‐1::mCherry transgene in the germline significantly rescued the oogenesis‐defective phenotype of the pafo‐1 mutants, suggesting that expression of the PAF1C in germ cells is required for oogenesis. Notably, overexpression of OMA‐1::GFP partially rescued the oogenesis defect in the pafo‐1 mutants. Based on our findings, we propose that the PAF1C promotes oogenesis in a cell‐autonomous manner by positively regulating the expression of genes involved in oocyte maturation.

NeuroGT: A brain atlas of neurogenic tagging CreER drivers for birthdate-based classification and manipulation of mouse neurons

Neuronal birthdate is one of the major determinants of neuronal phenotypes. However, most birthdating methods are retrospective in nature, allowing very little experimental access to the classified neuronal subsets. Here, we introduce four neurogenic tagging mouse lines, which can assign CreER-loxP recombination to neuron subsets that share the same differentiation timing in living animals and enable various experimental manipulations of the classified subsets. We constructed a brain atlas of the neurogenic tagging mouse lines (NeuroGT), which includes holistic image data of the loxP-recombined neurons and their processes across the entire brain that were tagged on each single day during the neurodevelopmental period. This image database, which is open to the public, offers investigators the opportunity to find specific neurogenic tagging driver lines and the stages of tagging appropriate for their own research purposes.

Improvement in image quality via the pseudo confocal effect in multidirectional digital scanned laser light-sheet microscopy

Multidirectional digital scanned laser light-sheet microscopy (mDSLM) cannot be used with the current pseudo confocal system to reduce blurring and background signals. The multiline scanning for light-sheet illumination and the simple image construction proposed in this study are alternative to the pseudo confocal system. We investigate the effectiveness of our pseudo confocal method combined with mDSLM on artificial phantoms and biological samples. The results indicate that image quality from mDSLM can be improved by the confocal effect; their combination is effective and can be applied to biological investigations.

QUAREP-LiMi: A community-driven initiative to establish guidelines for quality assessment and reproducibility for instruments and images in light microscopy

A modern day light microscope has evolved from a tool devoted to making primarily empirical observations to what is now a sophisticated , quantitative device that is an integral part of both physical and life science research. Nowadays, microscopes are found in nearly every experimental laboratory. However, despite their prevalent use in capturing and quantifying scientific phenomena, neither a thorough understanding of the principles underlying quantitative imaging techniques nor appropriate knowledge of how to calibrate, operate and maintain microscopes can be taken for granted. This is clearly demonstrated by the well-documented and widespread difficulties that are routinely encountered in evaluating acquired data and reproducing scientific experiments. Indeed, studies have shown that more than 70% of researchers have tried and failed to repeat another scientist's experiments, while more than half have even failed to reproduce their own experiments. One factor behind the reproducibility crisis of experiments published in scientific journals is the frequent underreporting of imaging methods caused by a lack of awareness and/or a lack of knowledge of the applied technique. Whereas quality control procedures for some methods used in biomedical research, such as genomics (e.g. DNA sequencing, RNA-seq) or cytometry, have been introduced (e.g. ENCODE), this issue has not been tackled for optical microscopy instrumentation and images. Although many calibration standards and protocols have been published, there is a lack of awareness and agreement on common standards and guidelines for quality assessment and reproducibility. In April 2020, the QUality Assessment and REProducibility for instruments and images in Light Microscopy (QUAREP-LiMi) initiative was formed. This initiative comprises imaging scientists from academia and industry who share a common interest in achieving a better understanding of the performance and limitations of microscopes and improved quality control (QC) in light microscopy. The ultimate goal of the QUAREP-LiMi initiative is to establish a set of common QC standards, guidelines, metadata models and tools, including detailed protocols, with the ultimate aim of improving reproducible advances in scientific research. This White Paper (1) summarizes the major obstacles identified in the field that motivated the launch of the QUAREP-LiMi initiative; (2) identifies the urgent need to address these obstacles in a grassroots manner, through a community of stakeholders including, researchers, imaging scientists, bioimage analysts, bioimage informatics developers, corporate partners, funding agencies, standards organizations, scientific publishers and observers of such; (3) outlines the current actions of the QUAREP-LiMi initiative and (4) proposes future steps that can be taken to improve the dissemination and acceptance of the proposed guidelines to manage QC. To summarize, the principal goal of the QUAREP-LiMi initiative is to improve the overall quality and reproducibility of light microscope image data by introducing broadly accepted standard practices and accurately captured image data metrics.

C. elegans spermatozoa lacking spe-45 are incapable of fusing with the oocyte plasma membrane

C. elegans spe-9 class genes encode sperm proteins with indispensable roles during fertilization. We have previously reported that spe-45 belongs to the spe-9 class, based on the finding that self-sperm of spe-45(tm3715) hermaphrodites were not consumed by fertilization. In this study, we directly observed live fertilization in the spermatheca of fem-1(hc17) females after mating with spe-45(tm3715) males. As expected, it was clearly shown that spe-45 mutant spermatozoa failed to fuse with the oocyte plasma membrane. Thus, our live imaging system for C. elegans fertilization seems to be useful for evaluation of the functions of male and female gametes.

Statistical image processing quantifies the changes in cytoplasmic texture associated with aging in Caenorhabditis elegans oocytes

Background

Oocyte quality decreases with aging, thereby increasing errors in fertilization, chromosome segregation, and embryonic cleavage. Oocyte appearance also changes with aging, suggesting a functional relationship between oocyte quality and appearance. However, no methods are available to objectively quantify age-associated changes in oocyte appearance.

Results

We show that statistical image processing of Nomarski differential interference contrast microscopy images can be used to quantify age-associated changes in oocyte appearance in the nematode Caenorhabditis elegans. Max–min value (mean difference between the maximum and minimum intensities within each moving window) quantitatively characterized the difference in oocyte cytoplasmic texture between 1- and 3-day-old adults (Day 1 and Day 3 oocytes, respectively). With an appropriate parameter set, the gray level co-occurrence matrix (GLCM)-based texture feature Correlation (COR) more sensitively characterized this difference than the Max–min Value. Manipulating the smoothness of and/or adding irregular structures to the cytoplasmic texture of Day 1 oocyte images reproduced the difference in Max–min Value but not in COR between Day 1 and Day 3 oocytes. Increasing the size of granules in synthetic images recapitulated the age-associated changes in COR. Manual measurements validated that the cytoplasmic granules in oocytes become larger with aging.

Conclusions

The Max–min value and COR objectively quantify age-related changes in C. elegans oocyte in Nomarski DIC microscopy images. Our methods provide new opportunities for understanding the mechanism underlying oocyte aging.

Toward understanding the dynamic state of 3D genome

The three-dimensional (3D) genome organization and its role in biological activities have been investigated for over a decade in the field of cell biology. Recent studies using live-imaging and polymer simulation have suggested that the higher-order chromatin structures are dynamic; the stochastic fluctuations of nucleosomes and genomic loci cannot be captured by bulk-based chromosome conformation capture techniques (Hi-C). In this review, we focus on the physical nature of the 3D genome architecture. We first describe how to decode bulk Hi-C data with polymer modeling. We then introduce our recently developed PHi-C method, a computational tool for modeling the fluctuations of the 3D genome organization in the presence of stochastic thermal noise. We also present another new method that analyzes the dynamic rheology property (represented as microrheology spectra) as a measure of the flexibility and rigidity of genomic regions over time. By applying these methods to real Hi-C data, we highlighted a temporal hierarchy embedded in the 3D genome organization; chromatin interaction boundaries are more rigid than the boundary interior, while functional domains emerge as dynamic fluctuations within a particular time interval. Our methods may bridge the gap between live-cell imaging and Hi-C data and elucidate the nature of the dynamic 3D genome organization.

BD5: An open HDF5-based data format to represent quantitative biological dynamics data

BD5 is a new binary data format based on HDF5 (hierarchical data format version 5). It can be used for representing quantitative biological dynamics data obtained from bioimage informatics techniques and mechanobiological simulations. Biological Dynamics Markup Language (BDML) is an XML (Extensible Markup Language)-based open format that is also used to represent such data; however, it becomes difficult to access quantitative data in BDML files when the file size is large because parsing XML-based files requires large computational resources to first read the whole file sequentially into computer memory. BD5 enables fast random (i.e., direct) access to quantitative data on disk without parsing the entire file. Therefore, it allows practical reuse of data for understanding biological mechanisms underlying the dynamics.

Microrheology for Hi-C Data Reveals the Spectrum of the Dynamic 3D Genome Organization

The one-dimensional information of genomic DNA is hierarchically packed inside the eukaryotic cell nucleus and organized in a three-dimensional (3D) space. Genome-wide chromosome conformation capture (Hi-C) methods have uncovered the 3D genome organization and revealed multiscale chromatin domains of compartments and topologically associating domains (TADs). Moreover, single-nucleosome live-cell imaging experiments have revealed the dynamic organization of chromatin domains caused by stochastic thermal fluctuations. However, the mechanism underlying the dynamic regulation of such hierarchical and structural chromatin units within the microscale thermal medium remains unclear. Microrheology is a way to measure dynamic viscoelastic properties coupling between thermal microenvironment and mechanical response. Here, we propose a new, to our knowledge, microrheology for Hi-C data to analyze the dynamic compliance property as a measure of rigidness and flexibility of genomic regions along with the time evolution. Our method allows the conversion of an Hi-C matrix into the spectrum of the dynamic rheological property along the genomic coordinate of a single chromosome. To demonstrate the power of the technique, we analyzed Hi-C data during the neural differentiation of mouse embryonic stem cells. We found that TAD boundaries behave as more rigid nodes than the intra-TAD regions. The spectrum clearly shows the dynamic viscoelasticity of chromatin domain formation at different timescales. Furthermore, we characterized the appearance of synchronous and liquid-like intercompartment interactions in differentiated cells. Together, our microrheology data derived from Hi-C data provide physical insights into the dynamics of the 3D genome organization.

PHi-C: deciphering Hi-C data into polymer dynamics

Genomes are spatiotemporally organized within the cell nucleus. Genome-wide chromosome conformation capture (Hi-C) technologies have uncovered the 3D genome organization. Furthermore, live-cell imaging experiments have revealed that genomes are functional in 4D. Although computational modeling methods can convert 2D Hi-C data into population-averaged static 3D genome models, exploring 4D genome nature based on 2D Hi-C data remains lacking. Here, we describe a 4D simulation method, PHi-C (polymer dynamics deciphered from Hi-C data), that depicts 4D genome features from 2D Hi-C data by polymer modeling. PHi-C allows users to interpret 2D Hi-C data as physical interaction parameters within single chromosomes. The physical interaction parameters can then be used in the simulations and analyses to demonstrate dynamic characteristics of genomic loci and chromosomes as observed in live-cell imaging experiments. PHi-C is available at https://github.com/soyashinkai/PHi-C.

Biologically constrained optimization based cell membrane segmentation in C. elegans embryos

BACKGROUND

Recent advances in bioimaging and automated analysis methods have enabled the large-scale systematic analysis of cellular dynamics during the embryonic development of Caenorhabditis elegans. Most of these analyses have focused on cell lineage tracing rather than cell shape dynamics. Cell shape analysis requires cell membrane segmentation, which is challenging because of insufficient resolution and image quality. This problem is currently solved by complicated segmentation methods requiring laborious and time consuming parameter adjustments.

RESULTS

Our new framework BCOMS (Biologically Constrained Optimization based cell Membrane Segmentation) automates the extraction of the cell shape of C. elegans embryos. Both the segmentation and evaluation processes are automated. To automate the evaluation, we solve an optimization problem under biological constraints. The performance of BCOMS was validated against a manually created ground truth of the 24-cell stage embryo. The average deviation of 25 cell shape features was 5.6%. The deviation was mainly caused by membranes parallel to the focal planes, which either contact the surfaces of adjacent cells or make no contact with other cells. Because segmentation of these membranes was difficult even by manual inspection, the automated segmentation was sufficiently accurate for cell shape analysis. As the number of manually created ground truths is necessarily limited, we compared the segmentation results between two adjacent time points. Across all cells and all cell cycles, the average deviation of the 25 cell shape features was 4.3%, smaller than that between the automated segmentation result and ground truth.

CONCLUSIONS

BCOMS automated the accurate extraction of cell shapes in developing C. elegans embryos. By replacing image processing parameters with easily adjustable biological constraints, BCOMS provides a user-friendly framework. The framework is also applicable to other model organisms. Creating the biological constraints is a critical step requiring collaboration between an experimentalist and a software developer.

In vivo Live Imaging of Calcium Waves and Other Cellular Processes during Fertilization in Caenorhabditis elegans

Fertilization calcium waves are a conserved trigger for animal development; however, genetic analysis of these waves has been limited due to the difficulty of imaging in vivo fertilization. Here we describe a protocol to image calcium dynamics during in vivo fertilization in the genetic animal model Caenorhabditis elegans. This protocol consists of germline microinjection of a chemical calcium indicator, worm immobilization, live imaging, and image processing that quantifies the calcium fluorescence in the oocyte region moving in the field-of-view during ovulation. This imaging protocol can also be used to image other cellular processes during in vivo fertilization in C. elegans, such as membrane fusion and cytoskeletal dynamics.

Phenotype Analysis Method for Identification of Gene Functions Involved in Asymmetric Division of Caenorhabditis elegans

In gene function analysis, it is arduous to identify gene function individually, and the way to screen out all involved genes according to a particular phenotype or disease usually shows us little information for a specific problem. We present a data-driven analysis system based on wild type (WT) embryos to study the concrete function of each gene associated with certain category of abnormal phenotypes. It can be applied to genes with very few RNAi embryos. Instead of presupposing the particular function of a gene, its function is confirmed by the statistical testing of built models. The scheme includes the following five: first, verify the to be detected genes and determine related recognized features according to the given category; second, compute the value of each feature based on WT embryos and merge them by principal component analysis (PCA); third, for each of the selected components of PCA, build a normal distribution and verify its normality; fourth, project the RNAi embryos to each component and probe them; and finally, analyze the more detailed functions of each gene based on the physical or biological meaning of each component. Choosing the first-round asymmetric division process of Caenorhabditis elegans as the phenotype, experimental results show that on the different aspects of the asymmetric division process, par-2, par-3, and let-754 are related to scalar differences; dcn-1 and mcm-5 are associated with the divergences of scalar variation, which may reflect the disaccord in development; and dcn-1, par-2, and par-3 are involved with morphological discrepancies.

SSBD: a database of quantitative data of spatiotemporal dynamics of biological phenomena

Motivation: Rapid advances in live-cell imaging analysis and mathematical modeling have produced a large amount of quantitative data on spatiotemporal dynamics of biological objects ranging from molecules to organisms. There is now a crucial need to bring these large amounts of quantitative biological dynamics data together centrally in a coherent and systematic manner. This will facilitate the reuse of this data for further analysis.

Results: We have developed the Systems Science of Biological Dynamics database (SSBD) to store and share quantitative biological dynamics data. SSBD currently provides 311 sets of quantitative data for single molecules, nuclei and whole organisms in a wide variety of model organisms from Escherichia coli to Mus musculus. The data are provided in Biological Dynamics Markup Language format and also through a REST API. In addition, SSBD provides 188 sets of time-lapse microscopy images from which the quantitative data were obtained and software tools for data visualization and analysis.

Availability and Implementation: SSBD is accessible at http://ssbd.qbic.riken.jp.

Contact:sonami@riken.jp

Biological Dynamics Markup Language (BDML): an open format for representing quantitative biological dynamics data

Motivation: Recent progress in live-cell imaging and modeling techniques has resulted in generation of a large amount of quantitative data (from experimental measurements and computer simulations) on spatiotemporal dynamics of biological objects such as molecules, cells and organisms. Although many research groups have independently dedicated their efforts to developing software tools for visualizing and analyzing these data, these tools are often not compatible with each other because of different data formats.

Results: We developed an open unified format, Biological Dynamics Markup Language (BDML; current version: 0.2), which provides a basic framework for representing quantitative biological dynamics data for objects ranging from molecules to cells to organisms. BDML is based on Extensible Markup Language (XML). Its advantages are machine and human readability and extensibility. BDML will improve the efficiency of development and evaluation of software tools for data visualization and analysis.

Availability and implementation: A specification and a schema file for BDML are freely available online at http://ssbd.qbic.riken.jp/bdml/.

Contact:sonami@riken.jp

Supplementary Information:Supplementary data are available at Bioinformatics online.

Evaluation of the effectiveness of simple nuclei-segmentation methods on Caenorhabditis elegans embryogenesis images

BACKGROUND

For the analysis of spatio-temporal dynamics, various automated processing methods have been developed for nuclei segmentation. These methods tend to be complex for segmentation of images with crowded nuclei, preventing the simple reapplication of the methods to other problems. Thus, it is useful to evaluate the ability of simple methods to segment images with various degrees of crowded nuclei.

RESULTS

Here, we selected six simple methods from various watershed based and local maxima detection based methods that are frequently used for nuclei segmentation, and evaluated their segmentation accuracy for each developmental stage of the Caenorhabditis elegans. We included a 4D noise filter, in addition to 2D and 3D noise filters, as a pre-processing step to evaluate the potential of simple methods as widely as possible. By applying the methods to image data between the 50- to 500-cell developmental stages at 50-cell intervals, the error rate for nuclei detection could be reduced to ≤ 2.1% at every stage until the 350-cell stage. The fractions of total errors throughout the stages could be reduced to ≤ 2.4%. The error rates improved at most of the stages and the total errors improved when a 4D noise filter was used. The methods with the least errors were two watershed-based methods with 4D noise filters. For all the other methods, the error rate and the fraction of errors could be reduced to ≤ 4.2% and ≤ 4.1%, respectively. The minimum error rate for each stage between the 400- to 500-cell stages ranged from 6.0% to 8.4%. However, similarities between the computational and manual segmentations measured by volume overlap and Hausdorff distance were not good. The methods were also applied to Drosophila and zebrafish embryos and found to be effective.

CONCLUSIONS

The simple segmentation methods were found to be useful for detecting nuclei until the 350-cell stage, but not very useful after the 400-cell stage. The incorporation of a 4D noise filter to the simple methods could improve their performances. Error types and the temporal biases of errors were dependent on the methods used. Combining multiple simple methods could also give good segmentations.

WDDD: Worm Developmental Dynamics Database

During animal development, cells undergo dynamic changes in position and gene expression. A collection of quantitative information about morphological dynamics under a wide variety of gene perturbations would provide a rich resource for understanding the molecular mechanisms of development. Here, we created a database, the Worm Developmental Dynamics Database (http://so.qbic.riken.jp/wddd/), which stores a collection of quantitative information about cell division dynamics in early Caenorhabditis elegans embryos with single genes silenced by RNA-mediated interference. The information contains the three-dimensional coordinate values of the outlines of nuclear regions and the dynamics of the outlines over time. The database provides free access to 50 sets of quantitative data for wild-type embryos and 136 sets of quantitative data for RNA-mediated interference embryos corresponding to 72 of the 97 essential embryonic genes on chromosome III. The database also provides sets of four-dimensional differential interference contrast microscopy images on which the quantitative data were based. The database will provide a novel opportunity for the development of computational methods to obtain fresh insights into the mechanisms of development. The quantitative information and microscopy images can be synchronously viewed through a web browser, which is designed for easy access by experimental biologists.

Cell-to-cell heterogeneity in cortical tension specifies curvature of contact surfaces in Caenorhabditis elegans embryos

In the two-cell stage embryos of Caenorhabditis elegans, the contact surface of the two blastomeres forms a curve that bulges from the AB blastomere to the P₁ blastomere. This curve is a consequence of the high intracellular hydrostatic pressure of AB compared with that of P₁. However, the higher pressure in AB is intriguing because AB has a larger volume than P₁. In soap bubbles, which are a widely used model of cell shape, a larger bubble has lower pressure than a smaller bubble. Here, we reveal that the higher pressure in AB is mediated by its higher cortical tension. The cell fusion experiments confirmed that the curvature of the contact surface is related to the pressure difference between the cells. Chemical and genetic interferences showed that the pressure difference is mediated by actomyosin. Fluorescence imaging indicated that non-muscle myosin is enriched in the AB cortex. The cell killing experiments provided evidence that AB but not P₁ is responsible for the pressure difference. Computer simulation clarified that the cell-to-cell heterogeneity of cortical tensions is indispensable for explaining the pressure difference. This study demonstrates that heterogeneity in surface tension results in significant deviations of cell behavior compared to simple soap bubble models, and thus must be taken into consideration in understanding cell shape within embryos.

Estimating intercellular surface tension by laser-induced cell fusion

Intercellular surface tension is a key variable in understanding cellular mechanics. However, conventional methods are not well suited for measuring the absolute magnitude of intercellular surface tension because these methods require determination of the effective viscosity of the whole cell, a quantity that is difficult to measure. In this study, we present a novel method for estimating the intercellular surface tension at single-cell resolution. This method exploits the cytoplasmic flow that accompanies laser-induced cell fusion when the pressure difference between cells is large. Because the cytoplasmic viscosity can be measured using well-established technology, this method can be used to estimate the absolute magnitudes of tension. We applied this method to two-cell-stage embryos of the nematode Caenorhabditis elegans and estimated the intercellular surface tension to be in the 30-90 µN m(-1) range. Our estimate was in close agreement with cell-medium surface tensions measured at single-cell resolution.

A proof of the DBRF-MEGN method, an algorithm for deducing minimum equivalent gene networks

Background

We previously developed the DBRF-MEGN (difference-based regulation finding-minimum equivalent gene network) method, which deduces the most parsimonious signed directed graphs (SDGs) consistent with expression profiles of single-gene deletion mutants. However, until the present study, we have not presented the details of the method's algorithm or a proof of the algorithm.

Results

We describe in detail the algorithm of the DBRF-MEGN method and prove that the algorithm deduces all of the exact solutions of the most parsimonious SDGs consistent with expression profiles of gene deletion mutants.

Conclusions

The DBRF-MEGN method provides all of the exact solutions of the most parsimonious SDGs consistent with expression profiles of gene deletion mutants.

Local cortical pulling-force repression switches centrosomal centration and posterior displacement in C. elegans

Centrosome positioning is actively regulated by forces acting on microtubules radiating from the centrosomes. Two mechanisms, center-directed and polarized cortical pulling, are major contributors to the successive centering and posteriorly displacing migrations of the centrosomes in single-cell-stage Caenorhabditis elegans. In this study, we analyze the spatial distribution of the forces acting on the centrosomes to examine the mechanism that switches centrosomal migration from centering to displacing. We clarify the spatial distribution of the forces using image processing to measure the micrometer-scale movements of the centrosomes. The changes in distribution show that polarized cortical pulling functions during centering migration. The polarized cortical pulling force directed posteriorly is repressed predominantly in the lateral regions during centering migration and is derepressed during posteriorly displacing migration. Computer simulations show that this local repression of cortical pulling force is sufficient for switching between centering and displacing migration. Local regulation of cortical pulling might be a mechanism conserved for the precise temporal regulation of centrosomal dynamic positioning.

Genome-wide prediction of genetic interactions in a metazoan

Genetic interactions provide information about genes and processes with overlapping functions in biological systems. For Saccharomyces cerevisiae, computational integration of multiple types of functional genomic data is used to generate genome-wide predictions of genetic interactions. However, this methodology cannot be applied to the vastly more complex genome of metazoans, and only recently has the first metazoan genome-wide prediction of genetic interactions been reported. The prediction for Caenorhabditis elegans was generated by computationally integrating functional genomic data from S. cerevisiae, C. elegans and Drosophila melanogaster. This achievement is an important step toward system-level understanding of biological systems and human diseases.

Computer simulations and image processing reveal length-dependent pulling force as the primary mechanism for C. elegans male pronuclear migration

A male pronucleus migrates toward the center of an egg to reach the female pronucleus for zygote formation. This migration depends on microtubules growing from two centrosomes associated with the male pronucleus. Two mechanisms were previously proposed for this migration: a "pushing mechanism," which uses the pushing force resulting from microtubule polymerization, and a "pulling mechanism," which uses the length-dependent pulling force generated by minus-end-directed motors anchored throughout the cytoplasm. We combined two computer-assisted analyses to examine the relative contribution of these mechanisms to male pronuclear migration. Computer simulation revealed an intrinsic difference in migration behavior of the male pronucleus between the pushing and pulling mechanisms. In vivo measurements using image processing showed that the actual migration behavior in Caenorhabditis elegans confirms the pulling mechanism. A male pronucleus having a single centrosome migrated toward the single aster. We propose that the pulling mechanism is the primary mechanism for male pronuclear migration.

Detection of nuclei in 4D Nomarski DIC microscope images of early Caenorhabditis elegans embryos using local image entropy and object tracking

Background

The ability to detect nuclei in embryos is essential for studying the development of multicellular organisms. A system of automated nuclear detection has already been tested on a set of four-dimensional (4D) Nomarski differential interference contrast (DIC) microscope images of Caenorhabditis elegans embryos. However, the system needed laborious hand-tuning of its parameters every time a new image set was used. It could not detect nuclei in the process of cell division, and could detect nuclei only from the two- to eight-cell stages.

Results

We developed a system that automates the detection of nuclei in a set of 4D DIC microscope images of C. elegans embryos. Local image entropy is used to produce regions of the images that have the image texture of the nucleus. From these regions, those that actually detect nuclei are manually selected at the first and last time points of the image set, and an object-tracking algorithm then selects regions that detect nuclei in between the first and last time points. The use of local image entropy makes the system applicable to multiple image sets without the need to change its parameter values. The use of an object-tracking algorithm enables the system to detect nuclei in the process of cell division. The system detected nuclei with high sensitivity and specificity from the one- to 24-cell stages.

Conclusion

A combination of local image entropy and an object-tracking algorithm enabled highly objective and productive detection of nuclei in a set of 4D DIC microscope images of C. elegans embryos. The system will facilitate genomic and computational analyses of C. elegans embryos.

Objective measurement of spindle orientation in early Caenorhabditis elegans embryo

The spindle orientation is a crucial piece of information to understand the development of embryo. The spindle forms during cell division and the cell divides along the spindle axis. Spindle orientation was measured in many different mutant embryos of Caenrohabditis elegans. However, the objectivity and the productivity of these measurements were low because these measurements were made manually. Here we present a system that automates the measurement of the spindle orientation in C. elegans embryo. Automation increases the objectivity and productivity of the measurement. We confirmed the applicability of the system by applying it to spindles during the second divisions in wild-type and mutant C. elegans embryos.

DBRF-MEGN method: an algorithm for deducing minimum equivalent gene networks from large-scale gene expression profiles of gene deletion mutants

MOTIVATION

Large-scale gene expression profiles measured in gene deletion mutants are invaluable sources for identifying gene regulatory networks. Signed directed graph (SDG) is the most common representation of gene networks in genetics and cell biology. However, no practical procedure that deduces SDGs consistent with such profiles has been developed.

RESULTS

We developed the DBRF-MEGN (difference-based regulation finding-minimum equivalent gene network) method in which an algorithm deduces the most parsimonious SDGs consistent with expression profiles of gene deletion mutants. Positive (or negative) directed edges representing positive (or negative) gene regulations are deduced by comparing the gene expression level between the wild-type and mutant. The most parsimonious SDGs are deduced using graph theoretical procedures. Compensation for excess removal of edges by restoring a minimum number of edges makes the method applicable to cyclic gene networks. Use of independent groups of edges greatly reduces the computational cost, thus making the method applicable to large-scale expression profiles. We confirmed the applicability of our method by applying it to the gene expression profiles of 265 Saccharomyces cerevisiae deletion mutants, and we confirmed our method's validity by comparing the pheromone response pathway, general amino acid control system, and copper and iron homeostasis system deduced by our method with those reported in the literature. Interpretation of the gene network deduced from the S. cerevisiae expression profiles by using our method led to the prediction of 132 transcriptional targets and modulators of transcriptional activity of 18 transcriptional regulators.

AVAILABILITY

The software is available on request.

A gene network inference method from continuous-value gene expression data of wild-type and mutants

In this paper we introduce a new inference method of a gene regulatory network from steady-state gene expression data. Our method determines a regulatory structure consistent with an observed set of steady-state expression profiles, each generated from wild-type and single deletion mutant of the target network. Our method derives the regulatory relationships in the network using a graph theoretic approach. The advantage of our method is to be able to deal with continuous values of steady-state data, while most of the methods proposed in past use a Boolean network model with binary data. Performance of our method is evaluated on simulated networks with varying the size of networks, indegree of each gene, and the data characteristics (continuous-value/binary), and is compared with that of predictor method proposed by Ideker et al. As a result, we show the superiority of using continuous values to binary values, and the performance of our method is much better than that of the predictor method.

Adrenocorticotropin-independent macronodular adrenocortical hyperplasia associated with multiple colon adenomas/carcinomas which showed a point mutation in the APC gene

We report a male Japanese with corticotropin (ACTH)-independent macronodular adrenocortical hyperplasia (AIMAH) associated with multiple colon adenomas/carcinomas. The plasma cortisol level was elevated with no diurnal rhythm and was not suppressed with dexamethasone. Basal plasma ACTH was unmeasurable but subnormally increased after administration of metyrapone or corticotropin releasing hormone. Both adrenals were resected and weighed 90g; the histopathologic findings were similar to those of AIMAH as previously reported. At least 21 colon lesions which were adenomas or carcinomas, were resected endoscopically or surgically. This is the second reported case of the association of AIMAH with multiple colon polyps. An APC gene point mutation was detected in the colon cancer tissue by polymerase chain reaction (PCR)-single strand conformation polymorphism (SSCP)/direct sequencing analysis at the putative splice acceptor site consensus sequence. However, no abnormality of APC gene was detected in the adrenocortical hyperplastic tissue. The possible etiological coexistence of these two diseases is discussed.

Splenic macrophages can modify steroidogenesis of Leydig cells

A large amount of LH/hCG treatment given to male rats is known to suppress the enzyme activity of cytochrome P450c17 in Leydig cells for 48 h. A high dose LH/hCG injection is also known to allow immunocytes, such as macrophages, to migrate into the testicular interstitial compartment. It has not been known, however, whether these cells play a role in that suppression. In this study, we examined if splenic macrophages have any effects on testosterone secretion from Leydig cells by culturing rat testicular interstitial cells (TIC). Splenic macrophages co-cultured with TIC significantly suppressed testosterone secretion. Macrophages co-cultured reduced both progesterone to testosterone conversion and the amount of cytochrome P450c17 mRNA. The conditioned medium (SMCM), prepared by culturing macrophages for 12 h, significantly reduced either testosterone secretion from TIC or progesterone to testosterone conversion by TIC. These results indicate that splenic macrophages suppress testosterone secretion from Leydig cells by suppressing the cytochrome P450c17 enzyme in vitro, and that this effect is mediated at least in part by some soluble factors secreted from macrophages. Splenic macrophages migrating into the testis after LH/hCG stimulation could play a role in suppressing cytochrome P450c17 in Leydig cells.

[Experimental studies on endocrine therapy for rat hepatocellular carcinoma]

Androgen receptor (AnR) and estrogen receptor (ER) are known to exist in human hepatocellular carcinoma (HCC), but the low binding capacity casts doubts on the efficacy of endocrinotherapy. However, we focussed on the favorable dissociation constant and displacement of AnR and ER. Efficacy of endocrinotherapy for HCC was investigated using rat HCC cell line (AH66F) resembling the properties of human HCC and sex hormone receptors. In rat HCC AH66F, we confirmed that the AnR and ER were both positive and were mentioned binding capacity, dissociation constant and displacement resembled those of human HCC. In rat HCC AH66F, administrations of Tamoxifen converted AnR and ER responses to negative. Rat HCC AH66F was transplanted intraperitoneally to Donryu rats, various endocrinotherapies administered and the number of survival days compared with a control group. In male and female rats, the number of survival days was both in the orchidectomied (p < 0.01) and the Tamoxifen treated (p < 0.001) group significantly prolonged. However, in the group treated with medroxyprogesterone acetate no significant differences were observed. Also, in experiments with AnR (-) and ER (-) rat HCC cell line AH60C all endocrinotherapies were ineffective. Above results confirmed the efficacy of endocrinotherapy for rat HCC with positive sex hormone receptors.