LoopSage: An energy-based Monte Carlo approach for the loop extrusion modeling of chromatin.

The connection between the patterns observed in 3C-type experiments and the modeling of polymers remains unresolved. This paper presents a simulation pipeline that generates thermodynamic ensembles of 3D structures for topologically associated domain (TAD) regions by loop extrusion model (LEM). The simulations consist of two main components: a stochastic simulation phase, employing a Monte Carlo approach to simulate the binding positions of cohesins, and a dynamical simulation phase, utilizing these cohesins' positions to create 3D structures. In this approach, the system's total energy is the combined result of the Monte Carlo energy and the molecular simulation energy, which are iteratively updated. The structural maintenance of chromosomes (SMC) protein complexes are represented as loop extruders, while the CCCTC-binding factor (CTCF) locations on DNA sequence are modeled as energy minima on the Monte Carlo energy landscape. Finally, the spatial distances between DNA segments from ChIA-PET experiments are compared with the computer simulations, and we observe significant Pearson correlations between predictions and the real data. LoopSage model offers a fresh perspective on chromatin loop dynamics, allowing us to observe phase transition between sparse and condensed states in chromatin.

Applicability of Scoring Calyculin A-Induced Premature Chromosome Condensation Objects for Dose Assessment Including for Radiotherapy Patients.

As an extension to a previous study, a linear calibration curve covering doses from 0 to 10 Gy was constructed and evaluated in the present study using calyculin A-induced premature chromosome condensation (PCC) by scoring excess PCC objects. The main aim of this study was to assess the applicability of this PCC assay for doses below 2 Gy that are critical for triage categorization. Two separate blind tests involving a total of 6 doses were carried out; 4 out of 6 dose estimates were within the 95% confidence limits (95% CL) with the other 2 just outside. In addition, blood samples from five cancer patients undergoing external beam radiotherapy (RT) were also analyzed, and the results showed whole-body dose estimates statistically comparable to the dicentric chromosome assay (DCA) results. This is the first time that calyculin A-induced PCC was used to analyze clinical samples by scoring excess objects. Although dose estimates for the pre-RT patient samples were found to be significantly higher than the mean value for the healthy donors and were also significantly higher than those obtained using DCA, all these pre-treatment patients fell into the same category as those who may have received a low dose (<1 Gy) and do not require immediate medical care during emergency triage. Additionally, for radiological accidents with unknown exposure scenario, PCC objects and rings can be scored in parallel for the assessment of both low- and high-dose exposures. In conclusion, scoring excess objects using calyculin A-induced PCC is confirmed to be another potential biodosimetry tool in radiological emergency particularly in mass casualty scenarios, even though the data need to be interpreted with caution when cancer patients are among the casualties.

Use of Atomic Force Microscopy in UVB-Induced Chromosome Damage Provides Important Bioinformation for Cell Damage Assessment.

The chromosomal structure derived from UVB-stimulated HaCaT cells was detected by atomic force microscopy (AFM) to evaluate the effect of UVB irradiation. The results showed that the higher the UVB irradiation dose, the more the cells that had chromosome aberration. At the same time, different representative types of chromosome structural aberrations were investigated. We also revealed damage to both DNA and cells under the corresponding irradiation doses. It was found that the degree of DNA damage was directly proportional to the irradiation dose. The mechanical properties of cells were also changed after UVB irradiation, suggesting that cells experienced a series of chain reactions from inside to outside after irradiation. The high-resolution imaging of chromosome structures by AFM after UVB irradiation enables us to relate the damage between chromosomes, DNA, and cells caused by UVB irradiation and provides specific information on genetic effects.

cudaMMC: GPU-enhanced multiscale Monte Carlo chromatin 3D modelling.

MOTIVATION: Investigating the 3D structure of chromatin provides new insights into transcriptional regulation. With the evolution of 3C next-generation sequencing methods like ChiA-PET and Hi-C, the surge in data volume has highlighted the need for more efficient chromatin spatial modelling algorithms. This study introduces the cudaMMC method, based on the Simulated Annealing Monte Carlo approach and enhanced by GPU-accelerated computing, to efficiently generate ensembles of chromatin 3D structures. RESULTS: The cudaMMC calculations demonstrate significantly faster performance with better stability compared to our previous method on the same workstation. cudaMMC also substantially reduces the computation time required for generating ensembles of large chromatin models, making it an invaluable tool for studying chromatin spatial conformation. AVAILABILITY AND IMPLEMENTATION: Open-source software and manual and sample data are freely available on https://github.com/SFGLab/cudaMMC.

Burden of disease of spinal muscular atrophy linked to chromosome 5q (5q-SMA) in Colombia.

OBJECTIVE: This article estimates the disease burden of 5q-SMA in Colombia by using the disability-adjusted life years (DALYs) metric. METHODS: Epidemiological data were obtained from local databases and medical literature and were adjusted in the DisMod II tool. DALYs were obtained by adding years of life lost due to premature death (YLL) and years lived with disability (YLD). RESULTS: The modeled prevalence of 5q-SMA in Colombia was 0.74 per 100,000 population. The fatality rate for all types was 14.1%. The disease burden of 5q-SMA was estimated at 4,421 DALYs (8.6 DALYs/100,000), corresponding to 4,214 (95.3%) YLLs and 207 (4.7%) YLDs. Most of the DALYs were accounted in the 2-17 age group. Of the total burden, 78% correspond to SMA type 1, 18% to type 2, and 4% to type 3. CONCLUSIONS: Although 5q-SMA is a rare disease, it is linked to a significant disease burden due to premature mortality and severe sequelae. The estimates shown in this article are important inputs to inform public policy decisions on how to ensure adequate health service provision for patients with 5q-SMA.

Myoid gonadal stromal tumours are characterised by recurrent chromosome-level copy number gains: molecular assessment of a multi-institutional series.

Myoid gonadal stromal tumours (MGST) represent a rare type of testicular sex cord-stromal tumour that has recently been recognised as a distinct entity by the World Health Organization (WHO) classification of genitourinary tumours. MGSTs affect adult men and have been reported to behave in an indolent fashion. Histologically, MGSTs are pure spindle cell neoplasms that coexpress SMA and S100 protein. Given that the molecular features of these neoplasms remain largely undescribed, we evaluated a multi-institutional series of MGSTs using DNA and RNA sequencing. This study included 12 tumours from 12 patients aged 28 to 57 years. Tumour sizes ranged from 0.6 to 4.3 cm. Aggressive histologic features, such as vascular invasion, necrosis, invasive growth, and atypical mitoses were invariably absent. Mitotic activity was low, with a median of less than 1 mitosis per 10 high power fields (HPF; maximum: 3 mitoses per 10 HPF). Molecular analyses did not identify recurrent mutations or gene fusions. All cases with interpretable copy number variant data (9/10 cases sequenced successfully) demonstrated a consistent pattern of chromosome arm-level and whole-chromosome-level copy number gains indicative of ploidy shifts, with recurrent gains involving chromosomes 3, 6, 7, 8, 9, 11, 12, 14q, 15q, 17, 18q, 20, and 21q. Similar findings have also been recognised in pure spindle cell and spindle-cell predominant sex cord-stromal tumours without S100 protein expression. MGSTs are characterised by ploidy shifts and may be part of a larger spectrum of spindle cell-predominant sex cord-stromal tumours, including cases without S100 protein expression.

Application of chromosome microarray analysis and karyotyping in diagnostic assessment of abnormal Down syndrome screening results.

BACKGROUND: Down syndrome (DS) is the most common congenital cause of intellectual disability and also leads to numerous metabolic and structural problems. This study aims to explore the application value of chromosomal microarray analysis (CMA) and karyotyping in prenatal diagnosis for pregnant women with abnormal DS screening results. METHODS: The study recruited 1452 pregnant women with abnormal DS screening results including 493 with an enlarged nuchal translucency thickness (NT ≥ 2.5 mm) and 959 with an abnormal second-trimester maternal serum biomarker screening results. They underwent amniocentesis to obtain amniotic fluid for CMA and karyotyping. RESULTS: CMA identified 74/1452 abnormal results, which was more efficient than karyotyping (51/1452, P < 0.05.) CMA is equivalent to traditional karyotyping for identifying aneuploidies. Compared to karyotyping CMA identified 1.90% more copy number variants (CNVs) ranging from 159Kb to 6496Kb. However, 34.4% of them were recurrent pathogenic CNVs associated with risk of neurodevelopmental disorders. CMA identified 13 variants of uncertain significance (VUS) results and 1 maternal uniparental disomy (UPD) of chromosome 7. Karyotyping identified 3 mosaic sex chromosome aneuploidy and 4 balanced translocation which could not be identified by CMA. In enlarged NT group, karyotyping identified 80.9% abnormal results while in serum screening group karyotyping identified 35.7%. However, the incidence of pathogenic/likely pathogenic (P/LP) CNVs was nearly the same in both groups. That was because aneuploidies and gross duplication/deletion were previously screened out by NT scan. CONCLUSIONS: CMA and karyotyping have both advantages and disadvantages in prenatal diagnosis of pregnant women with abnormal DS screening results. However, there was not enough evidence to support routine CMA in pregnant women with abnormal DS screening results.

Assessment of mitochondrial DNA viability ratio in day-4 biopsied embryos as an add-in to select euploid embryos for single embryo transfer.

The aim of this study was to assess mitochondrial DNA analysis as a predictor of the pregnancy potential of biopsied preimplantation embryos. The study included 78 blastomeres biopsied from day 4 cleavage stage euploid embryos. The embryo karyotype was confirmed by 24-chromosome preimplantation genetic testing for aneuploidies using the Illumina Next-Generation Sequencing (NGS) system. Mitochondria viability ratios (mtV) were determined from BAM files subjected to the web-based genome-analysis tool Galaxy. From this cohort of patients, 30.4% of patients (n = 34) failed to establish pregnancy. The mean mtV ratio [mean = 1.51 ± 1.25-1.77 (95% CI)] for this group was significantly (P < 0.01) lower compared with the embryo population that resulted in established pregnancies [mean = 2.5 ± 1.82-2.68 (95% CI)]. mtV multiple of mean (MoM) values were similarly significantly (P < 0.01) lower in blastocysts failing to establish pregnancy. At a 0.5 MoM cut-off, the sensitivity of mtV quantitation was 35.3% and specificity was 78.2%. The positive predictive value for an mtV value > 0.5 MoM was 41.4%. This study demonstrates the clinical utility of preimplantation quantification of viable mitochondrial DNA in biopsied blastomeres as a prognosticator of pregnancy potential.

Chromosomal-level genome of velvet bean (Mucuna pruriens) provides resources for L-DOPA synthetic research and development.

Mucuna pruriens, commonly called velvet bean, is the main natural source of levodopa (L-DOPA), which has been marketed as a psychoactive drug for the clinical management of Parkinson's disease and dopamine-responsive dystonia. Although velvet bean is a very important plant species for food and pharmaceutical manufacturing, the lack of genetic and genomic information about this species severely hinders further molecular research thereon and biotechnological development. Here, we reported the first velvet bean genome, with a size of 500.49 Mb and 11 chromosomes encoding 28,010 proteins. Genomic comparison among legume species indicated that velvet bean speciated ∼29 Ma from soybean clade, without specific genome duplication. Importantly, we identified 21 polyphenol oxidase coding genes that catalyse l-tyrosine to L-DOPA in velvet bean, and two subfamilies showing tandem expansion on Chr3 and Chr7 after speciation. Interestingly, disease-resistant and anti-pathogen gene families were found contracted in velvet bean, which might be related to the expansion of polyphenol oxidase. Our study generated a high-quality genomic reference for velvet bean, an economically important agricultural and medicinal plant, and the newly reported L-DOPA biosynthetic genes could provide indispensable information for the biotechnological and sustainable development of an environment-friendly L-DOPA biosynthesis processing method.

Geometrical Properties of the Nucleus and Chromosome Intermingling Are Possible Major Parameters of Chromosome Aberration Formation.

Ionizing radiation causes chromosome aberrations, which are possible biomarkers to assess space radiation cancer risks. Using the Monte Carlo codes Relativistic Ion Tracks (RITRACKS) and Radiation-Induced Tracks, Chromosome Aberrations, Repair and Damage (RITCARD), we investigated how geometrical properties of the cell nucleus, irradiated with ion beams of linear energy transfer (LET) ranging from 0.22 keV/µm to 195 keV/µm, influence the yield of simple and complex exchanges. We focused on the effect of (1) nuclear volume by considering spherical nuclei of varying radii; (2) nuclear shape by considering ellipsoidal nuclei of varying thicknesses; (3) beam orientation; and (4) chromosome intermingling by constraining or not constraining chromosomes in non-overlapping domains. In general, small nuclear volumes yield a higher number of complex exchanges, as compared to larger nuclear volumes, and a higher number of simple exchanges for LET < 40 keV/µm. Nuclear flattening reduces complex exchanges for high-LET beams when irradiated along the flattened axis. The beam orientation also affects yields for ellipsoidal nuclei. Reducing chromosome intermingling decreases both simple and complex exchanges. Our results suggest that the beam orientation, the geometry of the cell nucleus, and the organization of the chromosomes within are important parameters for the formation of aberrations that must be considered to model and translate in vitro results to in vivo risks.

Assessment of Multiway Interactions with Tri-C.

Tri-C is a chromosome conformation capture (3C) approach that can efficiently identify multiway chromatin interactions with viewpoints of interest. As opposed to pair-wise interactions identified in methods such as Hi-C, 4C, and Capture-C, the detection of multiway interactions allows researchers to investigate how multiple cis-regulatory elements interact together in higher-order structures in single nuclei and address questions regarding structural cooperation between these elements. Here, we describe the procedure for designing and performing a Tri-C experiment.

A chromosome-level genome of Brachymystax tsinlingensis provides resources and insights into salmonids evolution.

Brachymystax tsinlingensis Li, 1966 is an endangered freshwater fish with economic, ecological, and scientific values. Study of the genome of B. tsinlingensis might be particularly insightful given that this is the only Brachymystax species with genome. We present a high-quality chromosome-level genome assembly and protein-coding gene annotation for B. tsinlingensis with Illumina short reads, Nanopore long reads, Hi-C sequencing reads, and RNA-seq reads from 5 tissues/organs. The final chromosome-level genome size is 2,031,709,341 bp with 40 chromosomes. We found that the salmonids have a unique GC content and codon usage, have a slower evolutionary rate, and possess specific positively selected genes. We also confirmed the salmonids have undergone a whole-genome duplication event and a burst of transposon-mediated repeat expansion, and lost HoxAbß Hox cluster, highly expressed genes in muscle may partially explain the migratory habits of B. tsinlingensis. The high-quality B. tsinlingensis assembled genome could provide a valuable reference for the study of other salmonids as well as aid the conservation of this endangered species.

TADeus2: a web server facilitating the clinical diagnosis by pathogenicity assessment of structural variations disarranging 3D chromatin structure.

In recent years great progress has been made in identification of structural variants (SV) in the human genome. However, the interpretation of SVs, especially located in non-coding DNA, remains challenging. One of the reasons stems in the lack of tools exclusively designed for clinical SVs evaluation acknowledging the 3D chromatin architecture. Therefore, we present TADeus2 a web server dedicated for a quick investigation of chromatin conformation changes, providing a visual framework for the interpretation of SVs affecting topologically associating domains (TADs). This tool provides a convenient visual inspection of SVs, both in a continuous genome view as well as from a rearrangement's breakpoint perspective. Additionally, TADeus2 allows the user to assess the influence of analyzed SVs within flaking coding/non-coding regions based on the Hi-C matrix. Importantly, the SVs pathogenicity is quantified and ranked using TADA, ClassifyCNV tools and sampling-based P-value. TADeus2 is publicly available at https://tadeus2.mimuw.edu.pl.

Polymer Folding Simulations from Hi-C Data.

In the absence of a clear molecular understanding of the mechanism that stabilizes specific contacts in interphasic chromatin, we resort to the principle of maximum entropy to build a polymeric model based on the Hi-C data of the specific system one wants to study. The interactions are set by an iterative Monte Carlo algorithm to reproduce the average contacts summarized by the Hi-C map. The study of the ensemble of conformations generated by the algorithm can report a much richer set of information than the experimental map alone, including colocalization of multiple sites, fluctuations of the contacts, and kinetical properties.

Clinical utility and cost-effectiveness analysis of chromosome testing concomitant with chromosomal microarray of patients with constitutional disorders in a U.S. academic medical center.

Chromosomal microarray (CMA) is now widely used as first-tier testing for the detection of copy number variants (CNVs) and absence of heterozygosity (AOH) in patients with multiple congenital anomalies (MCA), autism spectrum disorder (ASD), developmental delay (DD), and/or intellectual disability (ID). Chromosome analysis is commonly used to complement CMA in the detection of balanced genomic aberrations. However, the cost-effectiveness and the impact on clinical management of chromosome analysis concomitant with CMA were not well studied, and there is no consensus on how to best utilize these two tests. To assess the clinical utility and cost-effectiveness of chromosome analysis concomitant with CMA in patients with MCA, ASD, DD, and/or ID, we retrospectively analyzed 3,360 postnatal cases for which CMA and concomitant chromosome analysis were performed in the Colorado Genetic Laboratory (CGL) at the University Of Colorado School Of Medicine. Chromosome analysis alone yielded a genetic diagnosis in two patients (0.06%) and contributed additional information to CMA results in 199 (5.92%) cases. The impact of abnormal chromosome results on patient management was primarily related to counseling for reproductive and recurrence risks assessment (101 cases, 3.01%) while a few (5 cases, 0.15%) led to changes in laboratory testing and specialist referral (25 cases, 0.74%). The incremental cost-effectiveness ratio (ICER) of combined testing demonstrated the cost of each informative chromosome finding was significantly higher for patients with clinically insignificant (CI) CMA findings versus clinically significant (CS) CMA results. Our results suggest that a stepwise approach with CMA testing with reflex to chromosome analysis on cases with CS CMA findings is a more cost-effective testing algorithm for patients with MCA, ASD, and/or DD/ID.

The economy of chromosomal distances in bacterial gene regulation.

In the transcriptional regulatory network (TRN) of a bacterium, the nodes are genes and a directed edge represents the action of a transcription factor (TF), encoded by the source gene, on the target gene. It is a condensed representation of a large number of biological observations and facts. Nonrandom features of the network are structural evidence of requirements for a reliable systemic function. For the bacterium Escherichia coli we here investigate the (Euclidean) distances covered by the edges in the TRN when its nodes are embedded in the real space of the circular chromosome. Our work is motivated by 'wiring economy' research in Computational Neuroscience and starts from two contradictory hypotheses: (1) TFs are predominantly employed for long-distance regulation, while local regulation is exerted by chromosomal structure, locally coordinated by the action of structural proteins. Hence long distances should often occur. (2) A large distance between the regulator gene and its target requires a higher expression level of the regulator gene due to longer reaching times and ensuing increased degradation (proteolysis) of the TF and hence will be evolutionarily reduced. Our analysis supports the latter hypothesis.

Quantitative assessment reveals the dominance of duplicated sequences in germline-derived extrachromosomal circular DNA.

Extrachromosomal circular DNA (eccDNA) originates from linear chromosomal DNA in various human tissues under physiological and disease conditions. The genomic origins of eccDNA have largely been investigated using in vitro-amplified DNA. However, in vitro amplification obscures quantitative information by skewing the total population stoichiometry. In addition, the analyses have focused on eccDNA stemming from single-copy genomic regions, leaving eccDNA from multicopy regions unexamined. To address these issues, we isolated eccDNA without in vitro amplification (naïve small circular DNA, nscDNA) and assessed the populations quantitatively by integrated genomic, molecular, and cytogenetic approaches. nscDNA of up to tens of kilobases were successfully enriched by our approach and were predominantly derived from multicopy genomic regions including segmental duplications (SDs). SDs, which account for 5% of the human genome and are hotspots for copy number variations, were significantly overrepresented in sperm nscDNA, with three times more sequencing reads derived from SDs than from the entire single-copy regions. SDs were also overrepresented in mouse sperm nscDNA, which we estimated to comprise 0.2% of nuclear DNA. Considering that eccDNA can be integrated into chromosomes, germline-derived nscDNA may be a mediator of genome diversity.

Efficiency and equity in origin licensing to ensure complete DNA replication.

The cell division cycle must be strictly regulated during both development and adult maintenance, and efficient and well-controlled DNA replication is a key event in the cell cycle. DNA replication origins are prepared in G1 phase of the cell cycle in a process known as origin licensing which is essential for DNA replication initiation in the subsequent S phase. Appropriate origin licensing includes: (1) Licensing enough origins at adequate origin licensing speed to complete licensing before G1 phase ends; (2) Licensing origins such that they are well-distributed on all chromosomes. Both aspects of licensing are critical for replication efficiency and accuracy. In this minireview, we will discuss recent advances in defining how origin licensing speed and distribution are critical to ensure DNA replication completion and genome stability.

Assessment of 3D Interactions Between Promoters and Distal Regulatory Elements with Promoter Capture Hi-C (PCHi-C).

Chromosome conformation capture and its variants interrogate population-average chromatin structure at a higher resolution and throughput than microscopic methods. Capture Hi-C is a variant tailored for the simultaneous assessment of all interactions with thousands of specific bait sequences, so is particularly suited to genome-wide studies of promoter interactions with distal regulatory elements, such as enhancers. We present the principles and methods for Promoter Capture Hi-C (PCHi-C), from experimental design to data analysis.