Next Article in Journal
The Effect of a Diet Enriched with Jerusalem artichoke, Inulin, and Fluoxetine on Cognitive Functions, Neurogenesis, and the Composition of the Intestinal Microbiota in Mice
Next Article in Special Issue
Potential Effects of High Temperature and Heat Wave on Nanorana pleskei Based on Transcriptomic Analysis
Previous Article in Journal
Significant Association of Estrogen Receptor-β Isoforms and Coactivators in Breast Cancer Subtypes
Previous Article in Special Issue
Pesticidal Toxicity of Phosphine and Its Interaction with Other Pest Control Treatments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

ParticleChromo3D+: A Web Server for ParticleChromo3D Algorithm for 3D Chromosome Structure Reconstruction

by
David Vadnais
and
Oluwatosin Oluwadare
*
Department of Computer Science, University of Colorado, Colorado Springs, CO 80918, USA
*
Author to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2023, 45(3), 2549-2560; https://doi.org/10.3390/cimb45030167
Submission received: 1 February 2023 / Revised: 9 March 2023 / Accepted: 11 March 2023 / Published: 17 March 2023
(This article belongs to the Collection Feature Papers in Current Issues in Molecular Biology)

Abstract

:
Understanding the three-dimensional (3D) structure of chromatin is invaluable for researching how it functions. One way to gather this information is the chromosome conformation capture (3C) technique and its follow-up technique Hi-C. Here, we present ParticleChromo3D+, a containerized web-based genome structure reconstruction server/tool that provides researchers with a portable and accurate tool for analyses. Additionally, ParticleChromo3D+ provides a more user-friendly way to access its capabilities via a graphical user interface (GUI). ParticleChromo3D+ can save time for researchers by increasing the accessibility of genome reconstruction, easing usage pain points, and offloading computational processing/installation time.

1. Introduction

DNA replication, gene regulation, and gene expression are just some of the areas that are affected by chromatin’s three-dimensional (3D) structure [1]. Traditionally, studying chromatin was done via microscopy tools, such as fluorescence in situ hybridization (FISH). FISH has been effective in showing that chromosomes organize in distant 3D territories [2]. Live-cell fluorescence microscopy has even been able to provide insight into the dynamic properties of living cell chromosomes [3]. Chromosome conformation capture (3C) was designed to allow for easier study of these 3D structures. 3C caused a paradigm shift in the way chromatins are studied because it often outperforms FISH [1,4,5]. 3C works through a biochemical procedure that binds and cuts areas of a genome and then measures the frequency at which two sections are bound together [4,6]. Each cut section is known as a chromatin bin, and each bin has a length that is specific to the number of base pairs in the bin [4,6]. Since the advent of 3C, many derivative technologies have been designed, such as Hi-C [7], chromosome conformation capture-on-chip (4C) [8], chromosome conformation capture carbon copy (5C) [9], tethered conformation capture (TCC) [10], and chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) [11]. In particular, Hi-C provides genome-wide chromatin information frequency (IF) data formatted as a contact matrix [7]. Hi-C uses the next-generation sequencing concepts of parallel sequencing and high-throughput sequencing to analyze read-pair interactions on an all-versus-all basis. This means that all of the genome’s read pairs can be sequenced at once [7]. In order to generate the IF data using Hi-C, the data must be pre-processed using data mapping and quality control [12,13]. This involves the process of mapping or aligning the read pairs to a reference genome using read-pair alignment algorithms, such as the Burrows–Wheeler aligner (BWA) [14] or the Bowtie 2 [15]. Following this process is a filtering of the alignment for quality control and, finally, a conversion of the read count data, based on the depth of the resolution, into contact matrix data or maps, which can be plain text, such as a two-dimensional (symmetric) or sparse (three-column) matrix, a binary .hic format [16], or the .cool/.mcool [17] file formats.
The 3D structure reconstruction field is currently focused on developing accurate and robust algorithms to predict the 3D chromosome and genome structure from the IF data produced from in situ Hi-C. In order for a researcher to interpret these Hi-C data, they often pre-process the data into the desired state of 3D positional data using one of a variety of tools [18,19,20,21,22,23,24,25,26,27,28,29,30]. Oluwadare O. et al. wrote a review paper about this in 2019, and, generally, 3D chromosome structure inference tools are based on either distance, contact, or probabilistic methodologies [31]. Distance methods work by first converting the IF data into relative distances and then using optimization techniques to optimize the structure into following the relative distances as closely as possible. Distance-based tools mostly differ from each other according to how the tools perform IF-to-distance conversion and which algorithm is used to optimize the result [31]. Contact methods use IF data directly to perform modeling [31,32]. Finally, probabilistic models measure the probabilistic measure that the chromatin bins will have contact [33]. These methods generally handle noise better than other methods, but they can be more time consuming [31].
In 2022, we published a paper in which we designed and tested a tool for chromosome reconstruction named ParticleChromo3D [18]. ParticleChromo3D is a distance-based tool that reconstructs 3D genome structures using particle swarm optimization (PSO) as its optimization algorithm. PSO is described in detail in our Section 2, but, in general, it is an optimization approach that uses many different particles to explore different configurations [34]. This exploration is completed by having the particles interact/iterate and move towards a local best and global best configuration simultaneously [34]. PSO can become stuck during local best optimizations and may not always find the global best solution [34]. During this, ParticleChromo3D was compared to ChromSDE [21], Chromosome3D [22], 3DMax [23], ShRec3D [24], LorDG [25], HSA [26], MOGEN [27], GEM [28], and PASTIS [29], and it performed well [18]. After creating ParticleChromo3D, we were left with some regrets in that it is not very portable, some knowledge of command line tools is needed to operate it, and it lacks usability and quality-of-life features. We created ParticleChromo3D+ as a containerized web-based implementation of ParticleChromo3D, designed to make it both easier to use externally and highly portable for offline usage.

2. Materials and Methods

The web server allows easy access to our three-dimensional (3D) structure optimizer algorithm without having to install Python and the required script dependencies. The algorithm developed in our previous work utilizes a particle swarm optimization (PSO) algorithm to minimize the distance between chromatin bins with the goal of using the interaction frequency to create expected Euclidean distances and then incrementally changing the spatial positioning data to minimize the 3D structure’s expected distance from its current distance [18]. We built the server by running a Docker container that hosts our representational state transfer (REST)-based back-end API and front-end GUI.

2.1. Particle Swarm Optimization and ParticleChromo3D

PSO was designed by Kennedy J. and Eberhart R. based on the flight patterns of birds [34]. PSO works by creating a set of particles (in our case, a matrix of chromatin bin locations in [x, y, z]) and then having each particle adjust its position based on a velocity that is a combination of the particle’s history and the best swarm position [18,35]. In general, the position and velocity functions are shown below in Equations (1) and (2), respectively [36].
P n + 1 = P n + V n
V n + 1 = w × V n + c 1 × R 1 × ( Ψ n P n ) + c 2 × R 2 × ( G n P n )
where
  • P n and P n + 1 are the position at a time stamp n and the position at the next time stamp n + 1 , respectively.
  • V n and V n + 1 are the velocity at a time stamp n and the velocity at the next time stamp n + 1 , respectively.
  • w is used to create inertia for the velocity. This helps reduce or increase the contribution of the individual particle’s past velocity.
  • c 1 and c 2 are the local and global weights, respectively. These weights are used to tell each particle if it should prioritize its information or the swarm information.
  • R 1 and R 2 are randomized values that increase the breadth of the geometries that the particles explore.
  • Φ n contains the position of the chromatin bin for the optimum structure that the individual particle has found. Each particle has its own structure made up of each chromatin bin’s position.
  • G n contains the position of the chromatin bin for the optimum structure that the swarm has found.
Specifically, in ParticleChromo3D (Figure 1), each particle runs one of these update equations on every cell of a matrix at each iteration [18]. This matrix contains continuous numeric data with six features representing the 3D position and velocity of each chromatin bin instance [18]. Each particle’s structure is then compared (see Section 2.2 for details) at every time step to see if it is better than the current best-found structure at all time steps and all particles. If it is the best time structure, it will be saved as the new best structure, and the process will continue until a stop condition is reached.

2.2. Scoring Intermittent 3D Structures

In order to know which particle contains the global best structure, we need to be able to rank them. We perform this by comparing the Euclidean distances with a user-selected loss function. We provide the root mean squared error (RMSE) [31], the mean squared error (MSE) [37], the sum of squared errors (SSE) [38], and the Huber loss [39] as options from a drop-down menu. The four loss functions, the RMSE, MSE, SSE, and Huber Loss, are defined below in Equations (7), (4), (5), and (6), respectively.
R M S E = 1 n c o u n t × ( d i D i ) 2
M S E = 1 n c o u n t × ( d i D i ) 2
S S E = ( d i D i ) n c o u n t
H u b e r L o s s = 1 2 × ( d i D i ) 2 if ( d i D i ) α α × ( | d i D i | 1 2 × α ) if x = 0
where:
  • d i is the distance between two chromatin bins found by our particle’s structure.
  • D i is the expected distance between two chromatin bins based on the IF data.
  • n c o u n t is the total number of chromatin bins.
  • α is a positive real number for alternating between the top and bottom loss functions. We set α to 0.5 .
The equation to find D i is [12,31]:
D i = 1 I F i , j β
where:
  • D i is the expected distance between two chromatin bins based on the IF data.
  • I F i , j is the information frequency between two chromatin bins i and j.
  • β is a conversion factor.

2.3. Scoring Final 3D Structures

When ParticleChromo3D+ is used, the results will be returned via email as a PDB file, the input parameters, the best optimization’s resulting Spearman correlation coefficient (SCC), and the best Pearson correlation coefficient (PCC). Both the SCC and PCC are defined below [31].
P C C = ( ( d i d ¯ ) × ( D i D ¯ ) ) ( d i d ¯ ) 2 × ( D i D ¯ ) 2
where
  • d i and D i are individual distances.
  • d ¯ and D ¯ are sample means of the distances.
S C C = ( x i x ¯ ) × ( y i y ¯ ) ( x i x ¯ ) 2 × ( y i y ¯ ) 2
where
  • x i and y i are the individual distances d i and D i converted into ranked variables.
  • x ¯ and y ¯ are sample means of the ranked distances.

2.4. Containerization

Containerization is a lightweight competitor to traditional virtual machines (VMs). Containers compete with VMs, reducing the load on the hypervisor [40]. Both containers and VMs, when implemented well, offer isolation, manageability, consolidation, and reliability [41]. Containers also offer better memory efficiency, short cloning times, and numerous continuous integration/continuous delivery (CI/CD) improvements over VMs when lightweight containers are used [42]. We chose to use containerization because we wanted our web server to be as lightweight and portable as possible. To try to gain the advantages of lightweight containers, we based our image on the Apache Tomcat official image sourced from Docker Hub. Our implementation ends up having an image size of roughly one gigabyte and was based on the containerization tool Docker.

2.5. Representational State Transfer

REST APIs are a mainstay in web-based development [43]. Alternatives to REST APIs include the simple object access protocol (SOAP), web service description language (WSDL), and gRPC (gRPC remote procedure calls) [43,44]. REST APIs assist in offloading storage and computation to cloud computers often through JavaScript object notation (JSON) or hypertext markup language (HTML) International Conference on Web Engineering. We chose REST for our API due to its widespread use, as shown in Figure 2. Figure 2 shows Google Trends’ search count from 13 February 2022 to 5 February 2023 and was sourced from Google Trends on 10 February 2023. An additional reason that we use REST over a non-HTML API is that our data volume does not currently seem to need an implementation such as gRPC, and these implementations would preclude the user from using non-HTML tools, such as Postman.

3. Results

ParticleChromo3D+ can be utilized in two ways. The first is to access our web server at http://particlechromo3d.online/ (accessed on 31 January 2023) using a web browser or through back-end services. For web browsers, we tested Firefox and Google Chrome. The second is to host a Docker server and run an instance of the ParticleChromo3D+ image as a container or pod.

3.1. Usage

Whether our server or a local container was used, the user can now either upload an IF matrix or use one of the example files therein. This can be done through the front-end GUI or a back-end REST interface via the command line. These implementations default to ports 8080 and 5001, respectively. The GUI is written using HTML and CSS. The REST interface is written using the Python Flask web framework, and it wraps a Python-based script generated for previous research [18].

3.1.1. Front-End Access

Upon submission of a job to the web server, it will run the particle swarm optimization-based algorithm developed in earlier research through the following process (Figure 3): First, the user can tweak the run parameters. We provide access to the parameters of swarm size, maximum iterations, change threshold, initialization value range, loss function, output file name, and email recipient upon job completion. We provided the default parameters of 15, 30,000, 0.000001, 1.0, and the root mean squared error based on past research. The user must provide a valid email address to begin processing data. Second, the user can choose the IF data they want to process. Third, the data are then sent to the back end and evaluated upon the user pressing the submit button. Finally, the results are emailed to the user and maintained online at a provided download URL.
The GUI is broken into three parts (Figure 4). First, the optimization parameters are available to be changed. This is labeled as section 1 in Figure 4. Here the user can change or use the provided default values for the Swarm Size, iteration count, error threshold, the random range for the creation of the initial or stating x, y, z coordinate of the 3D chromosome structure, the loss function, or the output file name. The only required field from the user is the provision of their email address to receive the output results. Once the desired optimization parameters have been selected, the user can either upload an IF matrix to process (Figure 4 section 2a) or use the provided 1 mb Hi-C contact map file for the GM12878 cell Hi-C dataset from Rao et al. [45] (Figure 4 section 2b).

3.1.2. Back-End API

The back end can be reached at port 5001 by default. The minimal expected use case is to hit the upload and process endpoints. An sample bash script with examples of many helpful curl commands is provided on GitHub. First, a POST request must be made to the upload endpoint with the desired IF matrix (Table 1). Then, a GET request must be made to the process endpoint with all of its parameters defined (Table 2).
Additionally, we provide access to three more endpoints named uploaded, download, and convert. The uploaded endpoint allows the user to retrieve an HTML formatted list of all of the available IF files for processing. An expected use case for this endpoint is verifying that the desired IF file has been successfully uploaded (accessed on 31 January 2023).
http://biomlearn.uccs.edu:5001/uploaded
The download endpoint allows the user to download the contents of a known name PDB file.
biomlearn.uccs.edu:5001/download?ofname=${filename}
Lastly, the convert endpoint allows the user to convert 3xN matrices to square matrices so long as they are uploaded in a tab-separated values format. This endpoint can be used to format the IF data for use with ParticleChromo3D+ or any other square-matrix-based solution.
biomlearn.uccs.edu:5001/convert?filename=${filename}
Examples of accessing all endpoints via curl can be found at in our GitHub repository under help/exampleCurlScript.bash.

3.2. Consistency

We tested our server on the GM12878 cell Hi-C dataset, GEO Accession number GSE63525 [45]. The normalized contact matrix was downloaded from the GSDB database with GSDB ID: OO7429SF [46]. We timed our web server’s runtime on Chromosomes 1, 10, and 20 ten times each. While doing this, we recorded the runtime, the best Pearson correlation coefficient, and the best Spearman correlation coefficient from processing start time to the time of email send off completion. The reason that we did not use the time until email received is because the results ended up varying too greatly due to external factors. As expected, the runtime did increase with chromosome size (Figure 5). The maximum runtime was for Chromosome 1 at 97.729 s. The SCC values had ranges of [0.0099, 0.0026, 0.0140] for Chromosomes 1, 10, and 20, respectively. In the same order, the ranges of the PCC values were [0.0060, 0.0068, 0.0188]. Both the SCC and PCC values remained consistent run after run (Figure 6 and Figure 7).

3.3. Starting/Extending a Local ParticleChromo3D+ Server

A user may desire to run a local instance of ParticleChromo3D+. Some reasons a user may desire a local instance are for data confidentiality (our research team will never see locally stored data) or offline development environments. If the user needs to run locally, we provide a parameterized way for them to run an instance of ParticleChromo3D+ without worrying about dependencies/builds through GitHub workflows and containerization. Additionally, the user may want to extend our container. A user may want to extend this in order to modify our algorithm or modify our process. By using a standard and open container base, we make it easier to extend our build directly in the Dockerfile or through the FROM notation. We developed ParticleChromo3D+ on a Windows 10 operating system and were able to seamlessly deploy the image on our Linux servers. This supports our goal of portability for the container.

3.3.1. Installation

Our Docker images can be found at https://github.com/OluwadareLab/ParticleChromo3D_Plus/actions (accessed on 31 January 2023) by selecting the desired build and then clicking on the file named particlechromo3d_image.tar.gz under the artifacts section. This image is built by GitHub workflows automatically. Now, the user must load the image with:
 docker image load -i particlechromo3d_image.tar.gz
 docker tag \${IMAGE_ID} particlechromo3d:latest
Once the image has been loaded, it can be run ephemerally or persistently with the instructions on our GitHub homepage. The basic command to run our image is:
docker run -d \
  -p 5001:5001 -p 8080:8080 \
  -e SERVICE_EMAIL=${YOUR_SVC_EMAIL} \
  -e HOSTNAME_BE=${YOUR_URL} \
  -e SERVICE_EMAIL_KEY=${KEY} \
  particlechromo3d:latest
where:
  • -d runs the container detached.
  • -p [external:internal] sets the ports.
  • -e SERVICE_EMAIL is the email address that will be used by the server.
  • -e HOSTNAME_BE is the DNS name of the server.
  • -e SERVICE_EMAIL_KEY is the password to the service email.

3.3.2. Extending the Image

If the user wants to extend the image, they can create a new Docker file and use FROM and our image name/tag, or they can build the image and extend our services. To build the image, download or clone our Git repository and then path a Docker CLI into the top-level directory where the Dockerfile is located. In this directory, the user will build the image by running:
 docker build -t particlechromo3D:latest.
If the user wants to make code changes, they should complete this before building the image.

4. Discussion

The ParticleChromo3D+ web server provides a GUI and API that give it additional use cases beyond traditional CLI tools. The GUI/API provide ease-of-use for users, and the containerization adds a portable way to convert IF data into 3D chromosome structures. Additionally, the API allows users to procedurally script calls to our server to generate the 3D structures. Through both the API and the front-end web server, researchers can offload the computational load from their servers and place it on ours. Researchers can instantaneously begin processing IF files without having to spend time setting up the algorithm environment. This is achieved while still being comparable to other algorithms in terms of the 3D chromosome structure’s accuracy. ParticleChromo3D+ also has the advantages of offloading work to our servers and removing the need to install required languages and libraries.

5. Conclusions

In this work, we present ParticleChromo3D+, a user-friendly web server at http://particlechromo3d.online/ (accessed on 31 January 2023) for predicting 3D chromosome structures from user-uploaded IF data using a particle swarm optimization algorithm [18]. ParticleChromo3D+ reduces the setup time, increases accessibility/usability, and offloads the computational load from research environments. We believe ParticleChromo3D+ is a valuable tool for accelerating genetic research through increasing the accessibility to 3D structure prediction.

Author Contributions

O.O. conceived the project. D.V. designed the algorithm. D.V. implemented the algorithm. D.V. performed the statistical and simulation analyses. O.O. and D.V. evaluated the results and wrote the manuscript. All authors reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by start-up funding from the University of Colorado, Colorado Springs (to O.O.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All real Hi-C data files are available from the GSDB database (accession number(s) OO7429SF). The generated models, all datasets used for all performed analyses, and the source code for ParticleChromo3D are available at https://github.com/OluwadareLab/ParticleChromo3D_Plus, (accessed on 31 January 2023).

Conflicts of Interest

The authors declare that they have no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CLICommand Line Interface
etc.et cetera
FISHfluorescence in situ hybridization
gRPCgRPC Remote Procedure Calls
HTMLHypertext Markup Language
i.e.id est
IFInteraction Frequency
JSONJavaScript Object Notation
MSEMean Squared Error
PCCPearson Correlation Coefficient
RESTREpresentational State Transfer
RMSERoot Mean Squared Error
SCCSpearman Correlation Coefficient
SOAPSimple Object Access Protocol
SSESum of Squared Errors
VMsvirtual machines
WSDLWeb Service Description Language
3CChromosome Conformation Capture

References

  1. Sati, S.; Cavalli, G. Chromosome conformation capture technologies and their impact in understanding genome function. Chromosoma 2017, 126, 33–44. [Google Scholar] [CrossRef] [Green Version]
  2. Bolzer, A.; Kreth, G.; Solovei, I.; Koehler, D.; Saracoglu, K.; Fauth, C.; Müller, S.; Eils, R.; Cremer, C.; Speicher, M.R.; et al. Three-Dimensional Maps of All Chromosomes in Human Male Fibroblast Nuclei and Prometaphase Rosettes. PLoS Biol. 2005, 3, E157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Masui, O.; Bonnet, I.; Le Baccon, P.; Brito, I.; Pollex, T.; Murphy, N.; Hupé, P.; Barillot, E.; Belmont, A.S.; Heard, E. Live-cell chromosome dynamics and outcome of X chromosome pairing events during ES cell differentiation. Cell 2011, 145, 447–458. [Google Scholar] [CrossRef] [Green Version]
  4. Dekker, J.; Rippe, K.; Dekker, M.; Kleckner, N. Capturing chromosome conformation. Science 2002, 295, 1306–1311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. De Wit, E.; De Laat, W. A decade of 3C technologies: Insights into nuclear organization. Genes Dev. 2012, 26, 11–24. [Google Scholar] [CrossRef] [Green Version]
  6. Denker, A.; de Laat, W. The second decade of 3C technologies: Detailed insights into nuclear organization. Genes Dev. 2016, 30, 1357–1382. [Google Scholar] [CrossRef] [Green Version]
  7. Lieberman-Aiden, E.; Van Berkum, N.L.; Williams, L.; Imakaev, M.; Ragoczy, T.; Telling, A.; Amit, I.; Lajoie, B.R.; Sabo, P.J.; Dorschner, M.O.; et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science 2009, 326, 289–293. [Google Scholar] [CrossRef] [Green Version]
  8. Simonis, M.; Klous, P.; Splinter, E.; Moshkin, Y.; Willemsen, R.; De Wit, E.; Van Steensel, B.; De Laat, W. Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture–on-chip (4C). Nat. Genet. 2006, 38, 1348–1354. [Google Scholar] [CrossRef]
  9. Dostie, J.; Richmond, T.A.; Arnaout, R.A.; Selzer, R.R.; Lee, W.L.; Honan, T.A.; Rubio, E.D.; Krumm, A.; Lamb, J.; Nusbaum, C.; et al. Chromosome Conformation Capture Carbon Copy (5C): A massively parallel solution for mapping interactions between genomic elements. Genome Res. 2006, 16, 1299–1309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Kalhor, R.; Tjong, H.; Jayathilaka, N.; Alber, F.; Chen, L. Genome architectures revealed by tethered chromosome confor-mation capture and population-based modeling. Nat. Biotechnol. 2012, 30, 90–98. [Google Scholar] [CrossRef] [Green Version]
  11. Li, G.; Fullwood, M.J.; Xu, H.; Mulawadi, F.H.; Velkov, S.; Vega, V.; Ariyaratne, P.N.; Mohamed, Y.B.; Ooi, H.S.; Tenna-koon, C.; et al. ChIA-PET tool for comprehensive chromatin interaction analysis with paired-end tag sequencing. Genome Biol. 2010, 11, R22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Pal, K.; Forcato, M.; Ferrari, F. Hi-C analysis: From data generation to integration. Biophys. Rev. 2019, 11, 67–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. MacKay, K.; Kusalik, A. Computational methods for predicting 3D genomic organization from high-resolution chromo-some conformation capture data. Briefings Funct. Genom. 2020, 19, 292–308. [Google Scholar] [CrossRef] [PubMed]
  14. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  15. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome 2009, 10, R25. [Google Scholar] [CrossRef] [Green Version]
  16. Durand, N.C.; Shamim, M.S.; Machol, I.; Rao, S.S.; Huntley, M.H.; Lander, E.S.; Aiden, E.L. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 2016, 3, 95–98. [Google Scholar] [CrossRef] [Green Version]
  17. Abdennur, N.; Mirny, L.A. Cooler: Scalable storage for Hi-C data and other genomically labeled arrays. Bioinformatics 2020, 36, 311–316. [Google Scholar] [CrossRef]
  18. Vadnais, D.; Middleton, M.; Oliwadare, O. ParticleChromo3D: A Particle Swarm Optimization algorithm for chromo-some 3D structure prediction from Hi-C data. BioData Mining. 2022, 15, 19. [Google Scholar] [CrossRef] [PubMed]
  19. Peng, C.; Fu, L.Y.; Dong, P.F.; Deng, Z.L.; Li, J.X.; Wang, X.T.; Zhang, H.Y. The sequencing bias relaxed characteristics of Hi-C derived data and implications for chromatin 3D modeling. Nucleic Acids Res. 2013, 41, e183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Collins, B.; Oluwadare, O.; Brown, P. ChromeBat: A Bio-Inspired Approach to 3D Genome Reconstruction. Genes 2021, 12, 1757. [Google Scholar] [CrossRef]
  21. Zhang, Z.; Li, G.; Toh, K.C.; Sung, W.K. Inference of spatial organizations of chromosomes using semi-definite embed-ding approach and Hi-C data. In Proceedings of the Annual International Conference on Research in Computational Molecular Biology, Beijing, China, 7–10 April 2013; Springer: Berlin/Heidelberg, Germany, 2013; pp. 317–332. [Google Scholar]
  22. Adhikari, B.; Trieu, T.; Cheng, J. Chromosome3D: Reconstructing three-dimensional chromosomal structures from Hi-C interaction frequency data using distance geometry simulated annealing. BMC Genom. 2016, 17, 886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Oluwadare, O.; Zhang, Y.; Cheng, J. A maximum likelihood algorithm for reconstructing 3D structures of human chro-mosomes from chromosomal contact data. BMC Genom. 2018, 19, 161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lesne, A.; Riposo, J.; Roger, P.; Cournac, A.; Mozziconacci, J. 3D genome reconstruction from chromosomal contacts. Nat. Methods 2014, 11, 1141–1143. [Google Scholar] [CrossRef]
  25. Trieu, T.; Cheng, J. 3D genome structure modeling by Lorentzian objective function. Nucleic Acids Res. 2017, 45, 1049–1058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zou, C.; Zhang, Y.; Ouyang, Z. HSA: Integrating multi-track Hi-C data for genome-scale reconstruction of 3D chromatin structure. Genome Biol. 2016, 17, 40. [Google Scholar] [CrossRef] [Green Version]
  27. Trieu, T.; Cheng, J. MOGEN: A tool for reconstructing 3D models of genomes from chromosomal conformation capturing data. Bioinformatics 2016, 32, 1286–1292. [Google Scholar] [CrossRef] [Green Version]
  28. Zhu, G.; Deng, W.; Hu, H.; Ma, R.; Zhang, S.; Yang, J.; Peng, J.; Kaplan, T.; Zeng, J. Reconstructing spatial organizations of chromosomes through manifold learning. Nucleic Acids Res. 2018, 46, e50. [Google Scholar] [CrossRef] [Green Version]
  29. Varoquaux, N.; Ay, F.; Noble, W.S.; Vert, J.P. A statistical approach for inferring the 3D structure of the genome. Bioinformatics 2014, 30, i26–i33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Di Pierro, M.; Zhang, B.; Aiden, E.L.; Wolynes, P.G.; Onuchic, J.N. Transferable model for chromosome architecture. Proc. Natl. Acad. Sci. USA 2016, 113, 12168–12173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Oluwadare, O.; Highsmith, M.; Cheng, J. An Overview of Methods for Reconstructing 3-D Chromosome and Genome Structures from Hi-C Data. Biol. Proced. Online 2019, 21, 7. [Google Scholar] [CrossRef]
  32. Paulsen, J.; Sekelja, M.; Oldenburg, A.R.; Barateau, A.; Briand, N.; Delbarre, E.; Shah, A.; Sørensen, A.L.; Vigouroux, C.; Buendia, B.; et al. Chrom3D: Three-dimensional genome modeling from Hi-C and nuclear lamin-genome contacts. Genome Biol. 2017, 18, 21. [Google Scholar] [CrossRef]
  33. Yaffe, E.; Tanay, A. Probabilistic modeling of Hi-C contact maps eliminates systematic biases to characterize global chromosomal architecture. Nat. Genet. 2011, 43, 1059–1065. [Google Scholar] [CrossRef] [PubMed]
  34. Kennedy, J.; Eberhart, R.C. Particle swarm optimization. In Proceedings of the 1995 IEEE International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  35. Mohapatra, R.; Saha, S.; Dhavala, S.S. Adaswarm: A novel pso optimization method for the mathematical equivalence of error gradients. arXiv 2020, arXiv:2006.09875. [Google Scholar]
  36. Bonyadi, M.R.; Michalewicz, Z. Particle swarm optimization for single objective continuous space problems: A review. Evol. Comput. 2017, 25, 1–54. [Google Scholar] [CrossRef]
  37. Gupta, V.; Kling, H.; Yilmaz, K.; Martinez, G. Decomposition of the mean squared error and NSE performance criteria: Implications for improving hydrological modelling. J. Hydrol. 2009, 344, 80–91. [Google Scholar] [CrossRef] [Green Version]
  38. Nainggolan, R.; Perangin-angin, R.; Simarmata, E.; Tarigan, F. Improved the Performance of the K-Means Cluster Using the Sum of Squared Error (SSE) optimized by using the Elbow Method. J. Phys. Conf. Ser. 2018, 1361, 012015. [Google Scholar] [CrossRef]
  39. Huber, P.J. A Robust Version of the Probability Ratio Test. Ann. Math. Stat. 1965, 36, 1753–1758. [Google Scholar] [CrossRef]
  40. Li, Z.; Kihl, M.; Lu, Q.; Andersson, J.A. Performance Overhead Comparison between Hypervisor and Container Based Virtualization. In Proceedings of the IEEE 31st International Conference on Advanced Information Networking and Applications (AINA), Taipei, Taiwan, 27–29 March 2017; pp. 955–962. [Google Scholar]
  41. Parashar, M.; AbdelBaky, M.; Rodero, I.; Devarakonda, A. Cloud Paradigms and Practices for Computational and Da-ta-Enabled Science and Engineering. Comput. Sci. Eng. 2013, 15, 10–18. [Google Scholar] [CrossRef]
  42. Nogueira, A.F.; Ribeiro, C.; Rela, M.; Craske, A. Improving La Redoute’s CI/CD Pipeline and DevOps Processes by Applying Machine Learning Techniques. In Proceedings of the 11th International Conference on the Quality of Information and Communications Technology (QUATIC), Coimbra, Portugal, 4–7 September 2018; pp. 282–286. [Google Scholar]
  43. Rodríguez, C.; Báez, M.; Daniel, F.; Casati, F.; Trabucco, J.C.; Canali, L.; Percannella, G. REST APIs: A Large-Scale Analysis of Compliance with Principles and Best Practices. In Proceedings of the International Conference on Web Engineering, Lugano, Switzerland, 6–9 June 2016. [Google Scholar]
  44. Giretti, A. Understanding the gRPC Specification. In Beginning gRPC with ASP.NET Core 6; Apress: Berkeley, CA, USA, 2022. [Google Scholar]
  45. Rao, S.S.; Huntley, M.H.; Durand, N.C.; Stamenova, E.K.; Bochkov, I.D.; Robinson, J.T.; Sanborn, A.L.; Machol, I.; Omer, A.D.; Lander, E.S.; et al. A 3D Map of the Human Genome at Kilobase Resolution Reveals Principles of Chromatin Looping. Cell 2014, 159, 1665–1680. [Google Scholar] [CrossRef] [Green Version]
  46. Oluwadare, O.; Highsmith, M.; Turner, D.; Lieberman-Aiden, E.; Cheng, J. GSDB: A database of 3D chromosome and genome structures reconstructed from Hi-C data. BMC Mol. Cell Biol. 2020, 21, 30. [Google Scholar] [CrossRef] [PubMed]
Figure 1. ParticleChromo3D [18] Process Diagram.
Figure 1. ParticleChromo3D [18] Process Diagram.
Cimb 45 00167 g001
Figure 2. Comparison of REST, SOAP, gRPC, and WSDL searches worldwide.
Figure 2. Comparison of REST, SOAP, gRPC, and WSDL searches worldwide.
Cimb 45 00167 g002
Figure 3. Standard expected use case.
Figure 3. Standard expected use case.
Cimb 45 00167 g003
Figure 4. Graphical user interface screenshot for the ParticleChromo3D+ webserver.
Figure 4. Graphical user interface screenshot for the ParticleChromo3D+ webserver.
Cimb 45 00167 g004
Figure 5. Runtime box plots by chromosome.
Figure 5. Runtime box plots by chromosome.
Cimb 45 00167 g005
Figure 6. SCC box plots by chromosome.
Figure 6. SCC box plots by chromosome.
Cimb 45 00167 g006
Figure 7. PCC box plots by chromosome.
Figure 7. PCC box plots by chromosome.
Cimb 45 00167 g007
Table 1. POST http://biomlearn.uccs.edu:5001/upload, (accessed on 31 January 2023).
Table 1. POST http://biomlearn.uccs.edu:5001/upload, (accessed on 31 January 2023).
ParameterTypeDescription
fileText FileThis is your IF matrix
Table 2. GET http://biomlearn.uccs.edu:5001/process, (accessed on 31 January 2023).
Table 2. GET http://biomlearn.uccs.edu:5001/process, (accessed on 31 January 2023).
ParameterTypeDescription
ifnameStringName of the target IF matrix text file.
ssIntegerNumber of swarms desired.
ittIntegerIteration Count. This is an exit condition for the optimization script that defines the maximum iterations before stopping.
thresholdDoubleThis is an exit condition for the optimization script that defines the threshold of the error minimum before stopping.
randRangeDoubleRange of values between which the (x, y, z) initial position can be randomly assigned for each chromatin bin.
lfInteger [0, 3]Bit mask for choosing the desired loss function.
outfileStringAllows the user to tag the out put PDB file’s filename.
emailStringDefines the email address to be messaged upon job completion.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vadnais, D.; Oluwadare, O. ParticleChromo3D+: A Web Server for ParticleChromo3D Algorithm for 3D Chromosome Structure Reconstruction. Curr. Issues Mol. Biol. 2023, 45, 2549-2560. https://doi.org/10.3390/cimb45030167

AMA Style

Vadnais D, Oluwadare O. ParticleChromo3D+: A Web Server for ParticleChromo3D Algorithm for 3D Chromosome Structure Reconstruction. Current Issues in Molecular Biology. 2023; 45(3):2549-2560. https://doi.org/10.3390/cimb45030167

Chicago/Turabian Style

Vadnais, David, and Oluwatosin Oluwadare. 2023. "ParticleChromo3D+: A Web Server for ParticleChromo3D Algorithm for 3D Chromosome Structure Reconstruction" Current Issues in Molecular Biology 45, no. 3: 2549-2560. https://doi.org/10.3390/cimb45030167

APA Style

Vadnais, D., & Oluwadare, O. (2023). ParticleChromo3D+: A Web Server for ParticleChromo3D Algorithm for 3D Chromosome Structure Reconstruction. Current Issues in Molecular Biology, 45(3), 2549-2560. https://doi.org/10.3390/cimb45030167

Article Metrics

Back to TopTop