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.

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.

Wet and dry biological data are potentially complementary. By visually integrating the initiation and developmental processes of organisms, we might reveal new causalities in biological data. Here we present an integrated visualization system for a causality network constructed from phenotypic developmental characters and their related scientific literature. To obtain the phenotypic characters, we applied bio-imaging informatics techniques to the data of wet experiments. The phenotypic character network was visually rendered in the CausalNet system, which provides both explanatory and verification visualization functions. Statistical analysis and scientific literature mining proved useful for determining the mechanisms underlying the phenotypic trait network. The validity of the system was confirmed in an application example and expert feedback on the developmental process of the nematode Caenorhabditis elegans. The discussed methodology is applicable to other multicellular organisms.

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.

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: [email protected].

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/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

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.

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.