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.

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.

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.

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.

Tobacco transcription repressors NtJAZ: Potential involvement in abiotic stress response and glandular trichome induction.

Members of the Jasmonate ZIM domain (JAZ) proteins act as transcriptional repressors in the jasmonate (JA) hormonal response. To characterize the potential roles of JAZ gene family in plant development and abiotic stress response, fifteen JAZs were identified based on the genome of Nicotiana tabacum. Structural analysis confirmed the presence of single Jas and TIFY motif. Tissue expression pattern analysis indicated that NtJAZ-2, -3, -5, and -10 were highly expressed in roots and NtJAZ-11 was expressed only in the cotyledons. The transcript level of NtJAZ-3, -5, -9, and -10 in the stem epidermis was higher than that in the stem without epidermis. Dynamic expression of NtJAZs exposed to abiotic stress and phytohormone indicated that the expression of most NtJAZs was activated by salicylic acid, methyl jasmonate, gibberellic acid, cold, salt, and heat stresses. With abscisic acid treatment, NtJAZ-1, -2, and -3 were not activated; NtJAZ-4, -5, and -6 were up-regulated; and the remaining NtJAZ genes were inhibited. With drought stress, the expression of NtJAZ-1, -2, -3, -4, -5, -6, -7, and -8 was up-regulated, whereas the transcript of the remaining genes was inhibited. Moreover, high concentration MeJA (more than 1 mM MeJA) had an effect on secreting trichome induction, but inhabited the plant growth. Nine NtJAZs may play important role in secreting trichome induction. These results indicate that the JAZ proteins are convergence points for various phytohormone signal networks, which are involved in abiotic stress responses.

A Polymer Model for the Quantitative Reconstruction of Chromosome Architecture from HiC and GAM Data.

It is widely believed that the folding of the chromosome in the nucleus has a major effect on genetic expression. For example, coregulated genes in several species have been shown to colocalize in space despite being far away on the DNA sequence. In this manuscript, we present a new, to our knowledge, method to model the three-dimensional structure of the chromosome in live cells based on DNA-DNA interactions measured in high-throughput chromosome conformation capture experiments and genome architecture mapping. Our approach incorporates a polymer model and directly uses the contact probabilities measured in high-throughput chromosome conformation capture experiments and genome architecture mapping experiments rather than estimates of average distances between genomic loci. Specifically, we model the chromosome as a Gaussian polymer with harmonic interactions and extract the coupling coefficients best reproducing the experimental contact probabilities. In contrast to existing methods, we give an exact expression of the contact probabilities at thermodynamic equilibrium. The Gaussian effective model reconstructed with our method reproduces experimental contacts with high accuracy. We also show how Brownian dynamics simulations of our reconstructed Gaussian effective model can be used to study chromatin organization and possibly give some clue about its dynamics.

Digestion-ligation-only Hi-C is an efficient and cost-effective method for chromosome conformation capture.

Chromosome conformation capture (3C) technologies can be used to investigate 3D genomic structures. However, high background noise, high costs, and a lack of straightforward noise evaluation in current methods impede the advancement of 3D genomic research. Here we developed a simple digestion-ligation-only Hi-C (DLO Hi-C) technology to explore the 3D landscape of the genome. This method requires only two rounds of digestion and ligation, without the need for biotin labeling and pulldown. Non-ligated DNA was efficiently removed in a cost-effective step by purifying specific linker-ligated DNA fragments. Notably, random ligation could be quickly evaluated in an early quality-control step before sequencing. Moreover, an in situ version of DLO Hi-C using a four-cutter restriction enzyme has been developed. We applied DLO Hi-C to delineate the genomic architecture of THP-1 and K562 cells and uncovered chromosomal translocations. This technology may facilitate investigation of genomic organization, gene regulation, and (meta)genome assembly.

Flexible and dynamic nucleosome fiber in living mammalian cells.

Genomic DNA is organized three dimensionally within cells as chromatin and is searched and read by various proteins by an unknown mechanism; this mediates diverse cell functions. Recently, several pieces of evidence, including our cryomicroscopy and synchrotron X-ray scattering analyses, have demonstrated that chromatin consists of irregularly folded nucleosome fibers without a 30-nm chromatin fiber (i.e., a polymer melt-like structure). This melt-like structure implies a less physically constrained and locally more dynamic state, which may be crucial for protein factors to scan genomic DNA. Using a combined approach of fluorescence correlation spectroscopy, Monte Carlo computer simulations, and single nucleosome imaging, we demonstrated the flexible and dynamic nature of the nucleosome fiber in living mammalian cells. We observed local nucleosome fluctuation (~50 nm movement/30 ms) caused by Brownian motion. Our in vivo/in silico results suggest that local nucleosome dynamics facilitate chromatin accessibility and play a critical role in the scanning of genome information.

Estimating divergence time and ancestral effective population size of Bornean and Sumatran orangutan subspecies using a coalescent hidden Markov model.

Due to genetic variation in the ancestor of two populations or two species, the divergence time for DNA sequences from two populations is variable along the genome. Within genomic segments all bases will share the same divergence-because they share a most recent common ancestor-when no recombination event has occurred to split them apart. The size of these segments of constant divergence depends on the recombination rate, but also on the speciation time, the effective population size of the ancestral population, as well as demographic effects and selection. Thus, inference of these parameters may be possible if we can decode the divergence times along a genomic alignment. Here, we present a new hidden Markov model that infers the changing divergence (coalescence) times along the genome alignment using a coalescent framework, in order to estimate the speciation time, the recombination rate, and the ancestral effective population size. The model is efficient enough to allow inference on whole-genome data sets. We first investigate the power and consistency of the model with coalescent simulations and then apply it to the whole-genome sequences of the two orangutan sub-species, Bornean (P. p. pygmaeus) and Sumatran (P. p. abelii) orangutans from the Orangutan Genome Project. We estimate the speciation time between the two sub-species to be thousand years ago and the effective population size of the ancestral orangutan species to be , consistent with recent results based on smaller data sets. We also report a negative correlation between chromosome size and ancestral effective population size, which we interpret as a signature of recombination increasing the efficacy of selection.

Spatial and topological organization of DNA chains induced by gene co-localization.

Transcriptional activity has been shown to relate to the organization of chromosomes in the eukaryotic nucleus and in the bacterial nucleoid. In particular, highly transcribed genes, RNA polymerases and transcription factors gather into discrete spatial foci called transcription factories. However, the mechanisms underlying the formation of these foci and the resulting topological order of the chromosome remain to be elucidated. Here we consider a thermodynamic framework based on a worm-like chain model of chromosomes where sparse designated sites along the DNA are able to interact whenever they are spatially close by. This is motivated by recurrent evidence that there exist physical interactions between genes that operate together. Three important results come out of this simple framework. First, the resulting formation of transcription foci can be viewed as a micro-phase separation of the interacting sites from the rest of the DNA. In this respect, a thermodynamic analysis suggests transcription factors to be appropriate candidates for mediating the physical interactions between genes. Next, numerical simulations of the polymer reveal a rich variety of phases that are associated with different topological orderings, each providing a way to increase the local concentrations of the interacting sites. Finally, the numerical results show that both one-dimensional clustering and periodic location of the binding sites along the DNA, which have been observed in several organisms, make the spatial co-localization of multiple families of genes particularly efficient.

Double duty for nuclear proteins--the price of more open forms of mitosis.

During cell division, eukaryotic cells pass on their genetic material to the next generation by undergoing mitosis, which segregates their chromosomes. During mitosis, the nuclear envelope, nuclear pore complexes and nucleolus must also be segregated. Cells achieve this in a range of different forms of mitosis, from closed, in which these nuclear structures remain intact, to open, in which these nuclear structures are disassembled. In between lies a smorgasbord of intermediate forms of mitosis, displaying varying degrees of nuclear disassembly. Gathering evidence is revealing links between the extent of nuclear disassembly and the evolution of new roles for nuclear proteins during mitosis. We propose that proteins with such double duties help coordinate reassembly of the nucleus with chromosomal segregation.