Skip to main content
Advertisement
  • Loading metrics

Agent-based modeling of nuclear chromosome ensembles identifies determinants of homolog pairing during meiosis

  • Ariana Chriss,

    Roles Conceptualization, Formal analysis, Investigation, Software, Writing – original draft, Writing – review & editing

    Current address: Department of Genetics and Genomic Sciences, Icahn School of Medicine, Mount Sinai, New York, New York, United States of America

    Affiliations Department of Mathematics and Statistics, Cleveland State University, Cleveland, Ohio, United States of America, Department of Biological, Geological, and Environmental Sciences, Cleveland State University, Cleveland, Ohio, United States of America

  • G. Valentin Börner ,

    Roles Conceptualization, Formal analysis, Funding acquisition, Project administration, Writing – original draft, Writing – review & editing

    g.boerner@csuohio.edu (GVB); s.d.ryan@csuohio.edu (SDR)

    Affiliations Department of Biological, Geological, and Environmental Sciences, Cleveland State University, Cleveland, Ohio, United States of America, Center for Gene Regulation in Health and Disease, Cleveland State University, Cleveland, Ohio, United States of America

  • Shawn D. Ryan

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing

    g.boerner@csuohio.edu (GVB); s.d.ryan@csuohio.edu (SDR)

    Affiliations Department of Mathematics and Statistics, Cleveland State University, Cleveland, Ohio, United States of America, Center for Applied Data Analysis and Modeling, Cleveland State University, Cleveland, Ohio, United States of America

Abstract

During meiosis, pairing of homologous chromosomes (homologs) ensures the formation of haploid gametes from diploid precursor cells, a prerequisite for sexual reproduction. Pairing during meiotic prophase I facilitates crossover recombination and homolog segregation during the ensuing reductional cell division. Mechanisms that ensure stable homolog alignment in the presence of an excess of non-homologous chromosomes have remained elusive, but rapid chromosome movements appear to play a role in the process. Apart from homolog attraction, provided by early intermediates of homologous recombination, dissociation of non-homologous associations also appears to contribute to homolog pairing, as suggested by the detection of stable non-homologous chromosome associations in pairing-defective mutants. Here, we have developed an agent-based model for homolog pairing derived from the dynamics of a naturally occurring chromosome ensemble. The model simulates unidirectional chromosome movements, as well as collision dynamics determined by attractive and repulsive forces arising from close-range physical interactions. Chromosome number and size as well as movement velocity and repulsive forces are identified as key factors in the kinetics and efficiency of homologous pairing in addition to homolog attraction. Dissociation of interactions between non-homologous chromosomes may contribute to pairing by crowding homologs into a limited nuclear area thus creating preconditions for close-range homolog attraction. Incorporating natural chromosome lengths, the model accurately recapitulates efficiency and kinetics of homolog pairing observed for wild-type and mutant meiosis in budding yeast, and can be adapted to nuclear dimensions and chromosome sets of other organisms.

Author summary

Pairing of homologous chromosomes is a key feature of multiple cellular processes including gene expression control, chromosome break repair, and chromosome segregation. Homolog pairing during meiosis is shared among all sexually reproducing eukaryotes. Mechanistic determinants of homology-specific chromosome alignment are presently unknown. We have developed an agent-based model where contributions of the entire chromosome set to the pairing process is taken into account, comprising both homologous and non-homologous chromosomal encounters. Predictions from the model are readily compared to experimental data from budding yeast, parameters can be adjusted to other cellular systems and predictions can be tested via experimental manipulation of the relevant chromosomal features.

Introduction

Double stranded DNA has an uncanny ability to find a homologous partner in DNA mixtures of staggering complexity. While homologs in somatic cells tend to occupy nuclear areas more distant than expected, somatic pairing nevertheless underlies important biological processes, including X chromosome inactivation and association of loci affected by genomic imprinting [1, 2]. During meiosis, homolog pairing is a key prerequisite for the separation of homologs to opposite spindle poles. When pairing is compromised, homologs fail to form crossovers, resulting in homolog nondisjunction and the formation of gametes with a surplus or deficit of one or several chromosomes. The resulting chromosomal imbalances are a leading cause for birth defects and still births [3].

In many organisms, meiotic pairing depends on recombination initiation via double strand breaks (DSBs), enzymatically induced by the Spo11 transesterase [4]. DSBs typically occur at a multitude of chromosomal sites, at different positions in different cells. DSB processing is closely associated with the homology search, a process whereby DNA breaks carrying 3’ single-stranded overhangs assess homology between nearby chromosomes at the DNA sequence level [5]. If matched, DSBs are processed via homologous recombination into crossovers as well as other recombination products [6, 7]. Crossovers involve the reciprocal exchange of chromosome arms between homologs at allelic positions. In addition to providing physical linkage between homologs and ensuring their attachment to opposite spindle poles, crossovers also increase genetic diversity [6, 7].

The timing and genetic requirements of homolog pairing have been extensively studied in several organisms [610]. In budding yeast, homologs are somatically paired in G1-arrested cells, unpair during premeiotic DNA replication and commence re-pairing as cells initiate homologous recombination [1013]. Around the time when homolog pairing is established, chromosomes also undergo rapid movements throughout prophase I, as a prerequisite for efficient pairing. Movements of chromosomes are mediated by motile cytoplasmic filaments, which drag chromosome ends (telomeres) through the semi-fluid nuclear envelope. Cytoplasmic motor proteins, actin in budding yeast and the dynein-microtubule complex in most other organisms, mediate nuclear chromosome movements due to the attachment of their ends via the conserved SUN-KASH protein complex where SUN proteins interact with chromosome ends and reach across the inner nuclear envelope whereas KASH proteins span the outer nuclear envelope providing a link between SUN proteins and the cytoplasmic filaments [1417]. Pairing is completed around the time when both DSB ends have undergone strand exchange, giving rise to double Holliday junctions, a critical precursor of crossovers [10, 13].

The synaptonemal complex (SC) is a proteinaceous structure that assembles between paired homologs and juxtaposes their axes closely at 100 nm, a distance conserved in most taxa [610]. Homologous chromosomes are considered as paired once they have associated at a distance less than or equal to 400 nm, corresponding to the distance between co-aligned homolog axes in absence of the SC [8].

Little is known about the molecular mechanism(s) of homolog pairing, in part due to experimental challenges. Fixated, surface-spread cells exhibit superior resolution, but they do not provide insights about chromosome dynamics during the pairing process and nuclear architecture may become distorted during sample preparation [10, 13]. Tracking of chromosome trajectories in live cells is hampered by limited resolution and potential effects of phototoxicity [18]. Importantly, such studies are limited to a small subset of chromosomes due to the necessity to fluorescently label individual chromosomes [10, 13, 19].

Several mathematical models have examined potential contributions of molecular processes to homolog pairing, including telomere attachment to the nuclear envelope, chromosome bending stiffness, and polymer chains exhibiting an excluded volume repulsive potential [18, 2023]. Moreover, a cellular automaton model was developed that examines random searching via chromosome shuffling [24]. Importantly, existing pairing models cannot be validated as they model interactions between individual chromosomes that are presently inaccessible to experimental analysis.

Here, we have developed an agent-based model (ABM), in which movements and interactions between individual chromosomes (i.e. agents) are simulated, thereby recapitulating the interaction dynamics of an entire nuclear chromosome ensemble. Rather than making assumptions about pairing dynamics along individual chromosomes, we use differential equations based on first principles that govern the movements of chromosomes as an outcome of interactions with all other chromosomes within the same nucleus. Our model allows for the analysis of individual trajectories for all chromosomes throughout meiotic prophase I, facilitating comparison with experimental data. Modifications of various chromosomal parameters, including chromosome number, size, and movement velocity as well as attractive and repulsive forces provide insights into the contributions of each of these factors to homolog pairing.

Materials and methods

Modeling approach

Our homologous pairing model considers three contributing processes, i.e., chromosome interactions (homologous and non-homologous collision dynamics), chromosome translation (directed movements), and random chromosome motion (thermal noise within the nucleus) (Fig 1). By deriving a model entirely from these first principles, we avoid the introduction of mathematical parameters not apparent from the underlying biological process. A key feature of our model is the inclusion of attractive and repulsive forces representing processes that stabilize homologous pairing or disrupt non-homologous interactions, respectively. Rather than focusing on a single homolog pair, our model captures trajectories of complete chromosome ensembles throughout prophase of meiosis I, facilitating comparison with experimental data and adjustments to species-specific features such as chromosome number and size as well as nuclear dimensions.

thumbnail
Fig 1. Determinants of chromosome dynamics during the meiotic homology search.

While searching for their homologous partner, chromosomes move about the cell nucleus and while performing a continuous random walk (dashed arrows) within the confines of the nuclear envelope (bold circles), with velocity and changes in direction determined by interactions with other chromosomes as well as thermal noise. (A) When homologs enter each other’s attractive radii (i.e., centers are 400 nm apart), they exert attractive forces (Eq (3); j hom) on each other which move them closer until they reach their respective exclusive radii (purple; here 50 nm) keeping them at a constant distance of 100 nm. They subsequently continue their effective “random walk” in a paired status moving as a single entity non-homologous with respect to all other chromosomes. (B) When a chromosome enters the repulsive radius of a non-homolog (red; here 400 nm), the repulsive force in (Eq (3); j nhom) diverts their movement at the angle of their encounter with new velocity proportional to the minimum interaction distance. For illustration, ring colors match the terms in the Morse potential governing the interactions between chromosomes in Eqs (1)–(3); also see Fig 2.

https://doi.org/10.1371/journal.pcbi.1011416.g001

To compute the net motion of meiotic chromosome sets, we consider pairwise interactions between chromosomes at short distances as provided by a semi-dilute solution [25, 26]. A two-dimensional model is developed to match the available two-dimensional experimental measurements from fixated, surface-spread yeast nuclei used to calibrate the model [13]. For the dynamics of the center of mass, , a coupled system of N equations can be proposed as (1) where N gives the number of chromosomes and Fj is the force between two chromosomes, which is obtained as the negative gradient of the potential energy U (2) The potential energy is defined using the Morse potential [28] (3) Three components in Eq (1) are noteworthy. First, there is a translational velocity term (v0di) defining a chromosome’s present straight-line motion as one component of its total velocity. We refer to this as the chromosome’s translational velocity because it is a straight-line motion in the direction the chromosome is already moving. The velocity (v0) is multiplied by a chromosome’s dimensionless orientation di = vi/‖vi‖, a unit vector in the direction of the chromosome’s velocity. di can be changed by collisions with other chromosomes or with the nuclear envelope as well as by thermal noise. In living cells, directed chromosome movement is typically generated by cytoplasmic motile filaments which are connected to chromosome ends (telomeres) via a protein complex that traverses the semi-fluid nuclear envelope (see Introduction) [17, 27]. In our model, chromosome end attachment to the nuclear envelope is reflected by the two-dimensional features of the nuclear volume simulating movement along an invariable z-surface.

Second, interaction forces between chromosomes are captured by Eqs (2) and (3). The force F is derived from the potential energy defined by the Morse potential U and is singularly determined by the distance between homologous or non-homologous pairs of chromosomes ‖xixj‖ [28]. A Morse potential is defined as a difference of Yukawa potential energies, assigning to homologous chromosome pairs a weak attractive (i.e. associative) force strength Ca (blue ring in Fig 1A; here a = 400 nm) and a weak repulsive (i.e. excluded volume) force strength Cr1 acting at very short distances and ensuring that paired homologs are kept at a fixed center-to-center distance rather than overlapping (purple ring in Fig 1; here r1 = 50 nm; see also [23]). Accordingly, interactions between homologs are determined by their distance: namely, attractive with increasing strength as distance decreases from 400 to 50 nm, and repulsive at or below 50 nm. For pairs of non-homologous chromosomes, the potential energy is provided by a purely repulsive force strength Cr that may model dissociation of strand exchange (red ring, in Fig 1B; here r = 400 nm), and also the excluded volume constraint Cr1 acting at a distance of r1 = 50 nm (see Table 1). Following a non-homologous encounter, chromosomes move away from each other while maintaining their orientation with the new velocity proportional to the interaction distance. Mathematically, the Yukawa potential has a simple formula for its derivative defined in Eq (2) which gives the corresponding interaction force and allows for straightforward numerical computation.

thumbnail
Table 1. Table of parameters used in the simulations.

Movement speed for spo11 hypomorph inferred from [15].

https://doi.org/10.1371/journal.pcbi.1011416.t001

The form of the interactions as a Yukawa decaying exponential ensures that forces between chromosomes become rapidly negligible beyond the effective interaction radius a or r respectively, see Eq (3). The three forces are at their respective maxima when chromosomal centers of mass overlap, and decrease exponentially from there, falling e.g., to ∼ 1/3 of their maxima when their distance has reached the effective radius. Beyond this effective radius, forces are essentially negligible compared to chromosome motion because these contributions of velocity are significantly smaller than the translation velocity, v0, or the random motion term. For example, when chromosome centers of mass are separated by three effective length radii, forces decrease by an order of magnitude from their maxima. Thus, all computed forces represent short-range interactions that do not affect the motion of either homologous or non-homologous chromosomes beyond one effective radius scale, as further evident from the Supplemental Movies [45] of representative simulation runs (see Fig 2; below).

thumbnail
Fig 2. Strength of exponentially decaying forces between non-homologous (short-range repulsive) and homologous (short-range attractive) chromosomes.

Vertical black lines represent the centers of mass, xi, and the colored regions indicate each chromosome’s force modeled as a decaying exponential. The first column shows the strength of the forces. The second column shows a cartoon illustration of the corresponding chromosome dynamics when they meet, and the third column contains representative frames from Supplemental Movie 1 [45]. The solid circles on the righthand side indicate the effective interaction radius (see Table 1) beyond which the force is negligible. (A-C) Illustration of exponentially decaying forces when two non-homologous chromosomes meet. The repulsive force felt by a non-homologous chromosome corresponds to where its center of mass crosses into the other chromosome’s repulsive region (denoted by * in A). (D-F) Similar illustrations but for an attractive homologous interaction.

https://doi.org/10.1371/journal.pcbi.1011416.g002

Within living cells, attractive and repulsive forces may result from a combination of diverse interactions. These interactions primarily entail homology search and strand exchange between a resected double-strand break (DSB) and intact template DNA. They may also encompass interactions between intact double-stranded DNA molecules. Attractive forces manifest during the formation and elongation of a heteroduplex between a DSB and a homologous template. This process is initiated by microhomologies, typically consisting of 8 base pairs [29]. Repulsive forces involve the dissociation of heteroduplexes containing internal or flanking mismatches, a phenomenon facilitated by ATP hydrolysis [30]. Factors influencing the success of pairing and strand exchange further include the degree of coiling in both the template and invading DNA [31] and, at closer proximity, electrostatic interactions [32]. Proteins, such as histones [33], RecA orthologs [30], mediator proteins [34], diverse ATPases [35], and the mismatch repair machinery [36, 37], can all modulate both attraction and dissociation processes.

Third, there is a term in Eq (1) accounting for thermal noise within the nucleus. This is a small white-noise process where can be modeled as a normally distributed random variable with mean zero and variance proportional to the time step. While the chromosomes are tethered to the nuclear envelope, their motion is not linear and thermal noise plays an active role at this small scale. This is best modeled as a random walk. Although the net motion is not that of a random walk, this component contributes biologically-relevant overall motion. This formulation has been successfully used to capture thermal noise in other models describing microscale biological processes [3844].

Our approach replaces microscopic details of the physical chromosome shape with a representation of the dynamics of its center of mass, which resembles the green fluorescent protein (GFP) dot utilized in the experimental setting. This allows for chromosome movement to be the direct readout of interaction frequencies between homologous and non-homologous chromosomes throughout the nucleus. Our model is simplifying the nucleus from three to two dimensions, providing a direct correspondence to the two-dimensional nature of the experimental data set used to calibrate the model where biological sample preparation involved fixation and flattening of the nucleus [13]. Furthermore, our modeling framework treats homolog pairing as an endpoint and does not consider the ensuing homolog segregation during meiosis I, again facilitating comparison to experimental data where cells were arrested at the pairing stage due to absence of the meiotic progression factor Ndt80 [13, 19].

Parameter choice

Table 1 shows the parameter values used in the initial simulations. While our model is adaptable to a wide range of parameters, it is important to identify realistic values for comparison with specific experimental observations. Meiosis in the budding yeast S. cerevisiae (the organism which provided the experimental data set used for model calibration) involves 16 homolog pairs or 32 individual chromosomes, corresponding to 32 coupled equations, each simulating the dynamics of the center of mass for an individual chromosome. Chromosomes may come in close contact with the nuclear boundary without ever crossing it due to an appropriate value set for the boundary repulsion strength and characteristic length scale using a purely repulsive Morse potential (similar to the pure repulsion case in Eq 3). While the intact yeast nucleus is roughly 2 μm in diameter, it is set here to 6.5 μm to account for its flattening and spreading during experimental data collection [5, 13]. Flat nuclei in the simulation result in chromosome motion exclusively within the two-dimensional cross-section shown in Fig 1. Experimentally determined distances between GFP-tagged chromosomes range between 200 nm (the resolution limit of light microscopy) and 6.0 μm [13].

While attractive and repulsive interactions between homologous and non-homologous chromosomes, respectively, can presently not be measured directly, their length scales can be estimated from features that have been determined in several biological systems (above). Many of these interactions involve evolutionarily conserved components of the homologous recombination machinery for which measurements are available [6, 7]. We use a “symmetric” model where attractive and repulsive forces act at the same length. The attractive length radius (a) is set to 400 nm corresponding to approximately 800 nucleotides of fully extended, single stranded DNA (assuming 0.5 nm per nucleotide; [46, 47]). During the homology search, such 5’ resected DSB tentacles may reach out from the broken chromosome to other chromosomes to probe homology [9]. The exponential decay of both attractive and repulsive forces further reflects single stranded DNA at meiotic DSBs that may range from 500–1500 nucleotides [47, 48] Thus, two chromosomes with centers of mass within 400 nm can initiate the homology search. The minimum homology length requirement for strand capture entails 8 nucleotides of homology; such initial contacts then are extended in nucleotide triplet steps [29]. Thus, a relatively short DSB tip region appears to be sufficient to make the first contact followed by homolog movement toward one another by heteroduplex extension. Once homologs are paired, they remain so henceforth and assume a single equation of motion ensuring they do not break apart, with the very short distance repulsive force keeping them at a constant distance of 50 nm which approximately corresponds to the conserved 100 nm width of the synaptonemal complex [6, 7].

Likewise, the repulsive length (r) is set to a 400 nm center of mass radius under consideration that invading resected DSB tentacles also generate a dissociative (repulsive) force. Dissociation of strand exchange due to insufficient homology likely involves the same forces as attraction with net movement in the opposite direction. Again both attractive as well as repulsive forces are significant only over short distances, i.e., when chromosome centers of mass are less than 400 nm apart, corresponding to less than 1/15 of the maximum distance provided by the 6.5 μm diameter of the nuclear area (Fig 3).

thumbnail
Fig 3. Effects of chromosome number on pairing efficiencies and kinetics.

Pairing frequencies for the indicated number of homolog pairs in 200 realizations are combined to compute the percent of homologs paired at the indicated time points. Average pairing levels of the indicated number of equally sized homolog pairs. Error bars are omitted for clarity.

https://doi.org/10.1371/journal.pcbi.1011416.g003

Our choice for the strengths of homologous and non-homologous interaction forces is motivated by kinetic considerations. Non-homologous interaction strength is directly related to the time spent exploring that interaction by a given chromosome. Thus, doubling the repulsive strength cuts the time spent in close proximity in half for non-homologous chromosomes (see Fig 2). By contrast, interactions between homologous chromosomes requires a longer time period for the assessment of homology. Hence, parameters need to reflect features of the probing process whether or not two chromosomes are homologous, but also not result in non-homologous? chromosomes jumping apart. It is further noteworthy that a given chromosome has an attractive interaction with one homologous chromosome, but a repulsive interaction with up to 30 non-homologous chromosomes, even though in the semi-dilute scenario, interactions typically occur with only one chromosome at a time.

Apart from collisions, the homology search involves the following processes. Directed chromosome movements increase the probability of chromosomal encounters allowing them to probe homologous and non-homologous DNA partners. The velocity of chromosome movements at 300 nm/s is taken from live-cell imaging of meiotic chromosome movements which range between 200 and 500 nm/s [1517]. Importantly, prophase I chromosome movements enhance chromosomal encounters. Without movements, the homology search for a given chromosome with a 400 nm attractive radius would be limited to approximately 5% of the nuclear area [49]. The parameter for chromosome motion due to thermal noise is provided from the diffusion constant of interphase chromatin which ranges between 102 and 103 nm2/sec [50].

At the start of each simulation, chromosomes are assigned unit velocities and random initial placements throughout the nucleus using a basic exclusion algorithm where a chromosome is placed and then its center of mass distance to all other chromosomes previously placed is computed [20]. If any initial distance is below 400 nm, another random location is chosen for the respective chromosome to avoid overlap in the initial chromosome placement, repeating this process until all chromosomes have been placed. Notably, the model is easily adaptable to alternative initial placement exclusion distances.

Results

Initial model and comparison with experimental data

In the reference experiment, homologous and non-homologous chromosome interactions, were monitored by retrieving cell aliquots from semi-synchronous meiotic cultures which either carry two GFP-tagged copies of chromosome III or single copies of GFP-tagged chromosomes III and II, respectively [13]. Pairing was inferred from the frequency of cells with GFP dots separated by less than 400 nm in the strain carrying GFP-tags in homologs adjusted for fortuitous co-localization as derived from measurements in cultures carrying GFP-tags on non-homologous chromosomes pairs. This analysis suggested that homolog pairing is at a minimum around t = 3h after transfer to meiosis medium, corresponding to the time of pre-meiotic S-phase, and reaches maximum levels at t = 7.5h, when most cells have entered the pachytene stage [13]. To model the complete homology search process throughout prophase I, we simulated chromosomal dynamics using Eqs (1)–(3). The numerical approach was designed in MATLAB [51]. A base code for the wild-type dynamics is available at [52] and corresponding Supplemental Movies are available at [45]. Importantly, the simulation allows concurrent tracking of all 16 homolog pairs and comparing their pairing dynamics to appropriately matched non-homologous partners within the same nucleus. Running Monte Carlo simulations with up to N = 32 chromosomes, typically 200 realizations were carried out and averaged.

We began assessing the validity of our model under the simplification that all chromosomes are of uniform length corresponding to the measured length of a mid-sized yeast chromosome at the pachytene stage [16] and a uniform reach corresponding to 400 nm. Parameters in Table 1 are chosen for the best correspondence with experimental data following exploration of the parameter space. We initially modeled pairing with a minimum set of two pairs of homologous chromosomes and subsequently increased the number of chromosome pairs to 16. Homologs were considered as paired when they became stably juxtaposed at 400 nm (see Modeling Approach), exploring pairing dynamics in dilute (less then 4 pairs) to semi-dilute (4 or more pairs) conditions. This study revealed unexpected effects of chromosome numbers on the efficiency and kinetics of pairing.

The model was started at the t = 3h time point, corresponding to minimum pairing levels in the experimental data set (above). For the scenario with 2 homolog pairs only 75% of homolog pairs completed pairing by t = 7h indicating that all homologs are paired in as few as 50% of cells, and increases at later times are negligible. An increase of homolog pairs to 4 increased the efficiency of pairing somewhat and this was increased further in the presence of 8 or 16 homolog pairs per nucleus. This reveals the qualitative difference between the dilute (2 pairs) and semi-dilute regime (4 or more pairs) in chromosome dynamics. Accordingly, with at least 8 homolog pairs, essentially all homologs have completed pairing by t = 9h. Thus, in the dilute scenario, isolated movements without chromosome interactions are predominant, whereas in the semi-dilute scenario, chromosome interactions dominate motion and such interactions are a prerequisite for successful homolog pairing.

Having established that chromosome number affects pairing efficiency, we next explored additional features of the scenario with 16 homolog pairs corresponding to the chromosome set in diploid budding yeast. In this analysis, homolog pairs were arranged based on their initial distances in ascending order. Homolog distances were then recorded over time and plotted at four time points corresponding to those examined experimentally for chromosome III (Fig 4A). Each data point represents average distances from 200 realizations. Initial distances between homologs in the simulation range from 0.7 μm to ∼ 6 μm (Fig 4A). Modeling suggests that chromosomes initially separated by 0.7 μm to 4 μm (chromosome indices 1 to 12) have completed pairing within 2 hours, whereas chromosomes separated by a larger distance require up to 4 hours to complete pairing (Fig 4A). These results suggest a nonlinear relationship between initial homolog distance and pairing kinetics.

thumbnail
Fig 4. The pairing model captures experimentally determined association kinetics between homologous and non-homologous chromosomes during wild-type meiosis.

(A) Modeling of pairing between 16 uniformly-sized partner chromosomes. For each modeling run, chromosome distances at t = 3h are used to uniquely index homologous chromosome pairs, which are arranged along the x-axis based on ascending initial distances. The y-axis indicates average distances from 200 realizations of the agent-based model Eqs (1)–(3) at the time points indicated by the color code. Error bars indicate standard deviations. The inset shows experimentally determined distances between a single pair of GFP-labeled yeast chromosome III measured in a synchronized meiotic culture in fixated nuclei (n > 100) (data from [13]). Nuclei are arranged based on distances between homologous GFP signals at a given time point. Note that x-axes are different in the inset due to the fact that in the experiment cells were observed at a given time point and then discarded, whereas in the simulation the same cell was tracked over time. (B) Results from modeling of distances between non-homologous chromosome pairs in the same nuclei analyzed in (A). For non-homologous pairs, each homolog partner is matched with a non-homologous chromosome that at t = 3h exhibits a distance optimally matched to that with its cognate homolog partner. Inset (B) shows experimentally determined distances between non-homologous GFP-tagged budding yeast chromosomes II and III (error bars indicate SD). For details on experimental conditions see [13]. The pairing distance is highlighted in gray at 400 nm. *Note the experimental data do not include information for t = 9h.

https://doi.org/10.1371/journal.pcbi.1011416.g004

Distances between non-homologous chromosome partners were monitored in the same simulation. Data points in this case represent any non-homologous chromosome initially placed at the same distance as the homologous partner chromosome in the same nucleus (Fig 4B). Accordingly, initial distances between non-homologous and homologous pairs are similar, yet non-homologous distances remain unchanged over time whereas homologous distances become progressively smaller. This confirms that the simulation is indeed specific for homologous pairing.

Comparison between simulation and experiment indicates good qualitative and quantitative convergence for homologous and non-homologous association kinetics (see insets Fig 4). Importantly, non-homologous GFP signals remain at similar distances throughout meiosis in both experiment and simulation, with little changes in inter-chromosomal distances. We conclude that the modeling framework efficiently recapitulates the chromosome dynamics during homologous pairing throughout prophase I. Notably, however, modeled pairing kinetics are derived from the entire set of 16 homolog pairs within the same nucleus whereas experimental data represent distances between a single GFP-tagged chromosome with its homologous partner or a single representative non-homologous partner measured in different nuclei and in different cultures.

Several additional differences between experimental data and modeling are noteworthy. In the simulation, all chromosome pairs are placed to an unpaired starting position and released into the moving chromosome ensemble. The experiment, by nature represents a more complex situation: (i) Due to the incomplete synchrony within a meiotic culture, a cell sample taken at a given time point comprises a mixture of cells that have progressed to various degrees in the pairing process. (ii) Moreover, prior to entry into meiosis, homologs are paired somatically in G1-arrested cells, potentially as a special case of yeast meiosis [10, 11, 13, 19], meaning that both laggards that have not yet completed unpairing and those undergoing pairing ahead of the population average will artificially increase the pairing frequency. Accordingly, a substantial subset of cells exhibits homologous pairing even at the time of minimum pairing levels (see insets of Fig 4A), either because they are still somatically paired or because they have already progressed to post-replicative pairing. (iii) While the model places unpaired chromosomes randomly on an idealized two-dimensional nuclear area, the experimental “area” is created through spreading and/or squashing of a three-dimensional cell in an imperfectly controlled manner potentially compressing homologs that are separated along the z-axes, a complication experimentally addressed via comparison with cells harboring GFP signals on non-homologous chromosomes. Thus, unpaired homolog pairs may end up in close vicinity, even though this close vicinity may not have existed when the nucleus was in its original 3D state. This experimental artifact is further evident from substantial frequencies of pairing between non-homologous chromosomes in the experiment (see inset Fig 4B).

Modeling predicts delayed pairing for the three smallest yeast chromosomes

The 16 budding yeast chromosomes range in length from 250 kilobasepairs (kbp) to 1,500 kbp. Chromosome size affects several meiotic processes, including the timing and/or density of initiating DSBs and crossovers [53, 54]. Chromosome length may further affect the cumulative range of attractive and repulsive forces as well as the available space for other chromosomes to move in [55]. We therefore took the length of yeast chromosomes into account by scaling all attractive and repulsive radii in Table 1 by each chromosome’s length in relation to the average chromosome length. The larger effective radius takes into account the higher number of DSBs engaged in the homology search along a larger chromosome, even though the reach of an individual DSB would likely remain unchanged. The inner repulsive radius enforcing binding distance between homologs was kept at 50 nm, independent of chromosome length, consistent with the uniform width of the synaptonemal complex.

As before, distances between homologous chromosomes were plotted as a function of time, but results were sorted by increasing chromosome lengths rather than by initial chromosome distance (Fig 5A). For one of the realizations of the entire homology search of size-adjusted homolog pairs, see [45] (or https://zenodo.org/records/10246589). Again, initial chromosome distances are on average 3μm for all chromosomes, independent of length, corresponding to half of the diameter of the cell nucleus. As time progresses, chromosomes longer than 400 kbp have completed homologous pairing by t = 7h in essentially all nuclei, whereas the three shortest homolog pairs ranging between 250 and 320 kbp remain unpaired, even at t = 9h, as suggested by their final distances of > 400 nm (Fig 5). Thus, the three shortest chromosomes take substantially longer to pair than the rest of chromosomes, revealing a nonlinear effect of chromosome length on pairing kinetics. Importantly, these distinct kinetics are distinct for small homolog pairs, yet absent for the same chromosomes and their equidistant non-homologous chromosomes in the same nuclei (Fig 5B). Hence, our model predicts that smaller chromosomes exhibit unique pairing requirements. For a direct comparison of Chromosome III dynamics between simulation and experiment see S6 Fig.

thumbnail
Fig 5. Effects of chromosome size on pairing kinetics and efficiency.

(A) Average distances over time between homolog pairs where repulsive and attractive radii are adjusted proportionally to the sizes of actual yeast chromosomes. Note that the three smallest homolog pairs are markedly slower in achieving pairing than all other chromosomes. (B) Non-homologous chromosomes equidistant to each of the two homologs were identified at t = 3h and their distances were monitored throughout the simulation in the same set of model nuclei monitored in (a). (200 realizations, error bars indicate SD). The gray rectangle highlights pairing distances at or below 400 nm.

https://doi.org/10.1371/journal.pcbi.1011416.g005

Several additional features of the model are evident from inspection of the Supplemental Movie Files [45]: Both attractive and repulsive forces affect the path of other chromosomes only at short distances, as the path of a given chromosome is only altered when their effective radii overlap. Moreover, even when homologous chromosomes approach each other, they may remain closely aligned for extended periods, but fail to complete pairing at that time and separate again, e.g., due to repulsive interaction with a non-homologous chromosome. Such dynamics may explain the mixed association previously observed during experimental live-cell imaging studies where GFP-tagged homologs approached and separated without completing pairing [18].

Effect of chromosome movement velocity on pairing

We next used our model to examine the role of chromosome translation velocity, v0, on pairing efficiencies and kinetics. Even though faster chromosome movements would be expected to uniformly accelerate homolog pairing, these simulations identified a velocity threshold for pairing, which particularly impacts mid-sized and large chromosomes. In the simulations described previously, the velocity parameter in Eq (1) was set to v0 = 300 nm/s (see Table 1). Reducing the chromosome translation velocity by 50% (to 150 nm/s) essentially eliminates pairing when average pairing frequencies for size-adjusted chromosomes were monitored over time (red line; Fig 6A). Increasing chromosome velocity in 30 nm/s increments to 210 nm/s results only in minor increases in pairing efficiencies, with on average only 30% of homolog pairs achieving pairing by t = 9h (Fig 6A, blue). Surprisingly, the four smallest chromosomes (sized below 500 kbp) disproportionately contribute to pairing at velocities at or below 210 nm/s, whereas pairing is essentially abrogated for mid-sized and larger chromosomes (Fig 6B). When chromosome translation velocity is further increased by 30 nm/s to 240 nm/s, this results in a dramatic improvement of pairing efficiency for all chromosomes, with now ∼ 50% of homolog pairing completed by t = 9h, regardless of chromosome size (orange; Fig 6A and 6B). Pairing efficiency and kinetics can be further improved for all chromosomes by an increase of velocity to 270 nm/s (light blue), whereas neither a further 30 nm/s incremental increase nor a doubling in chromosome velocity to 600 nm/s has a substantial effect (Fig 6A and 6B; black and purple).

thumbnail
Fig 6. Modeling identifies a critical threshold of chromosome movement velocity for efficient homolog pairing.

(A) Dots indicate the average pairing levels over time of all 32 true-sized chromosomes. Black indicates the velocity of chromosome movements in the wild-type model in Figs 4 and 5 (300 nm/s). Chromosomes fail to pair at velocities around 150 nm/s and only increases above 240 nm/s have a substantial impact. For a more detailed analysis of pairing kinetics of (a) see S1 Fig. (B) Effects of movement velocity for actual chromosome sizes. The graph shows pairing levels of chromosomes of increasing sizes at t = 9h, indicating that pairing efficiencies of larger chromosomes are more dramatically affected by changes in velocity than those of smaller chromosomes.

https://doi.org/10.1371/journal.pcbi.1011416.g006

Thus, instead of a linear relationship between chromosome velocity and pairing, our model predicts a threshold effect where velocities at or above 240 nm/s dramatically improve pairing, a threshold that particularly impacts pairing of mid-sized and larger chromosomes. Note that at higher velocities, chromosomes move further into the repulsive radii of their non-homologous partners and are therefore repelled faster, resulting in an accelerated homology search. Together, these findings suggest that at velocities below 240 nm/s, contributions from random diffusion are more pronounced, interfering with directional chromosome movements. As random motion becomes dominant over directed movement, chromosomes explore less nuclear area, making homologous encounters less likely. Moreover, at lower translational velocities, repulsive interactions between non-homologous chromosomes are less frequent, diminishing effects of excluded regions defined by the presence of non-homologous chromosomes.

Contributions of attractive and repulsive force strengths to homolog pairing

Next, we examined relative contributions of attractive and repulsive forces to pairing dynamics by appropriate parameter changes in Eqs. (1) to (3). The base model (Fig 7A), comprising size adjusted chromosomes with equal strengths and reach of attractive as well as repulsive forces and movements at 300 nm/s, was adjusted by setting either the attractive force strength Ca = 0 (Fig 7B) or by inversely setting the repulsive force strength Cr = 0 (Fig 7C). These simulations indicate that neither attractive nor repulsive forces alone achieve pairing, at least on the time scale examined here, yet with some important differences. With repulsion only, the repulsive forces progressively drive homologs together via volume exclusion, with slower pairing kinetics affecting only smaller chromosomes (Fig 7B). In contrast, with attractive forces only, homologous chromosomes are drawn together only when they enter the local proximity of each other, and pairing proceeds exceedingly slowly, but with little or no effect of chromosome size (Fig 7C).

thumbnail
Fig 7. Contributions of attractive and repulsive forces on pairing efficiencies and kinetics.

(A) Pairing wild-type model using true chromosome lengths and a standard translational movement velocity of 300 nm/s as primarily studied herein. In (B) the attractive strength Ca = 0 in the WT model is used to highlight the effect of repulsive interactions alone. In (C) the reverse is true, Cr = 0 in the WT model is used to highlight the effect of attraction alone. (200 realizations, error bars indicate SD). The pairing distance is highlighted by a gray rectangle.

https://doi.org/10.1371/journal.pcbi.1011416.g007

We also examined whether a dominant attractive force could result in homolog pairing by increasing the attractive force Ca by an order of magnitude compared to the base model with an equivalent reduction of the repulsive force Cr. Indeed, this parameter adjustment results in almost instantaneous pairing of all but the smallest chromosomes (see S5 Fig). We note, however, that this scenario is unrealistic, as it is equivalent to a 3-fold increase in the effective radius, resulting in an effective radius of roughly 1,200 nm for average-sized chromosomes, further corresponding to one-third of the nuclear area and to an average resection length of 2.4kb of fully extended single stranded DNA (S5 Fig). While symmetric attractive and repulsive forces more faithfully capture the dynamics observed in the reference experiment (Fig 5), this permutation of the model demonstrates how systems with faster or slower pairing can be captured by adjusting the magnitude of both contributing interaction forces.

Contributions of chromosome elasticity and orientation to homolog pairing

To better account for the rod-like shape of condensed prophase I chromosomes, we next extended the basic modeling framework to account for the presumed polymer properties of chromosomes. Following an approach previously developed for linear active polymers [56], we modeled each chromosome as an active dumbbell where two beads, as defined by the base model, are connected by a bond force Fb = ∇Ub. Here the bond potential and resulting force are expressed as (4) where λ is a Lagrange multiplier, kB is the Boltzmann constant, T is the temperature, γT is a friction coefficient of translational motion and r2 = (r2r1)2 the bond vector whose average length enforces the chromosome size, 〈r2〉 = 2. Following this approach we choose kbT/γT = D for the diffusion coefficient. This replaces the radial model above with an elastic dumbbell formed by two connected beads each with its own attractive and repulsive forces, but now explicitly adding the additional features of elasticity and orientation of the chromosome (see Fig 8A). In the equations of motion the bond force term is added as an additional contribution to the motion of each chromosome, but only affects the relative position of each bead in the active dumbbell. In terms of attractive and repulsive interactions, each of the beads interacts with the two beads on each of the other chromosomes, more faithfully representing the elongated physical structure of chromosomes.

thumbnail
Fig 8. Polymer chain model for chromosomes incorporates flexibility and orientation.

(A) Polymer chain model as an active dumbbell where each end is represented by the single bead model (1) with an additional term to ensure they stay together (4). (B) Representative still from a simulation of the active dumbbell movie, see Supplemental Movie 3 [45].

https://doi.org/10.1371/journal.pcbi.1011416.g008

Simulations were carried out and described with a still of the nucleus and all 32 active dumbbell chromosomes (see Fig 9B and Supplemental Movie 2 at [45]). Even though pairing is completed somewhat faster in the model with dumbbells compared to those with circular chromosomes with a single center of mass, both models provide similar results. Thus, the computational complexity added by the active dumbbell approach does not significantly alter the results, at least with the parameters chosen here. However, the active dumbbell version of the simulation may prove useful in situations, e.g., when interactions with unusually large chromosomes are under investigation.

thumbnail
Fig 9. Polymer chain model as an active dumbbell.

(A) Average distances over time between homolog pairs where repulsive and attractive radii are adjusted proportionally to the sizes of actual yeast chromosomes. Note that the three smallest homolog pairs are markedly slower in achieving pairing than all other chromosomes. (B) Non-homologous chromosomes equidistant to each of the two homologs were identified at t = 3h and their distances were monitored throughout the simulation in the same set of model nuclei monitored in (A). (200 realizations, error bars indicate SD). The gray rectangle highlights pairing distances at or below 400 nm.

https://doi.org/10.1371/journal.pcbi.1011416.g009

Adaptation of the pairing model to a mutant with reduced DSBs

Building on the expanded model that takes chromosome sizes and optimized velocities into account, we next investigated the case of a meiotic mutant for which experimental data are available [13]. Absence of meiotic DSBs, e.g., in a spo11 null mutant, essentially abrogates homolog pairing, consistent with a central role of recombination intermediates in establishing and/or stabilizing homolog pairing [10, 19, 57]. A decrease of initiating DSBs to around 30% of wild type in a hypomorphic spo11 mutant (spo11-HA/spo11-HA-Y135F; hereafter spo11-HA/yf) results in delayed, though largely efficient pairing [13]. Model parameters were adjusted to accommodate the fact that reduced DSB abundance would likely reduce both the cumulative attractive and repulsive forces exerted by these intermediates. Furthermore, chromosome translation velocity was reduced consistent with experimental data that indicate a reduction of the average chromosome velocity by 20% in a spo11 null mutant (to 110 nm/s from 140 nm/s in wild type; see Fig 4A in [15]). Accordingly, for the spo11 hypomorph, in one simulation we reduced the chromosome translation velocity from 300 nm/s to 230 nm/s.

Actual chromosome sizes were used for this simulation, yet results were plotted based on initial chromosome distances to facilitate comparison with the experimental data set in both wild-type SPO11 (Fig 10A) and spo11 hypomorph (Fig 10B and 10C). We considered two parameter sets to model the spo11 hypomorph. In Scenario I, the movement speed was reduced to 77% of the wild type, and attractive strength was reduced to one-third to represent the mutant strain’s decrease in DSBs available for the homology search (Fig 10B). In Scenario II, we considered exclusively the DSB reduction in spo11-HA/yf by decreasing to one-third both the attractive and repulsive force strengths, but kept the chromosome movement speed at wild-type levels. For a specific realization of the homology search in mutant spo11, see Supplemental Movie 3 [45] (or https://zenodo.org/records/10246589). Reducing both movement speed and attractive forces delays pairing of widely separated homolog pairs indefinitely (see Fig 10C), whereas all chromosome distances are affected more similarly when both attractive and repulsive contributions are reduced (Fig 10B). Importantly, these distinct responses demonstrate the versatility of our model to predict potential contributions of various parameters to pairing dynamics.

thumbnail
Fig 10. Modeling of pairing kinetics at reduced DSB abundance.

(a) Wild-type model with true-sized chromosomes. Results differ from Fig 4A) due to use of true sized rather than uniformly sized chromosomes. The inset shows experimental data for yeast chromosome III in hypomorphic spo11-HA/yf. Note the x-axis in the inset shows different nuclei derived from aliquots at the indicated time points, whereas the same cells were tracked through time during simulations. (b) Scenario I for hypomorphic spo11 with 3-fold reductions of the wild-type levels of both attractive and repulsive forces. For non-homologous chromosome distances see S3 Fig. (c) Scenario II for spo11 hypomorphic mutant with reduced movement speed and attractive force strength, with experimental observations for spo11-HA/yf shown in the inset. Translational movement velocity is reduced to 77% of wild type levels (230 nm/s), consistent with slower chromosome movements observed in spo11 [15], and attractive force is reduced 3-fold representing decreased attractive forces exerted by fewer DSBs (see Table 1). The pairing distance is highlighted in gray at 400 nm. (200 realizations, error bars indicate SD).

https://doi.org/10.1371/journal.pcbi.1011416.g010

Discussion

How homologs pair during meiosis in the presence of an excess of non-homologous chromosomes is presently unknown. To explore contributions to the pairing process of the entire nuclear chromosome ensemble, we have developed an agent-based mathematical model derived from first principles that takes into account both attractive forces between homologs and dissociative forces between non-homologous chromosomes (see Methods section). In many organisms, homologs enter meiosis separated by distances that far exceed the reach of resected DSBs that could assess homology. This necessitates a process that ensures initial homolog co-localization potentially provided by chromosome movements together with non-homologous repulsive forces. Results from our simulations suggest that repulsive forces together with chromosome movements are indeed a key determinant of bringing homologs into close vicinity. Both attractive and repulsive forces exert their effects over the same short distances that are within the range of the single stranded region of resected DSBs. Repulsive interactions between non-homologous chromosomes create excluded regions within the nucleus driving homologs into close vicinity, thereby facilitating close-range attractive pairing interactions. Our simulations further demonstrate that repulsive forces are most effective when chromosome numbers rise above a certain threshold, likely by reducing the available area per chromosome. Attractive forces come into play once distances between homologs are sufficiently small.

A repulsive force that shortens the time spent in non-homologous interactions is a key feature of our pairing model. It represents molecular processes that dissociate non-homologous chromosomes from each other. While pronounced contributions of a dissociative/repulsive force to homolog pairing may appear counter-intuitive, several mutant phenotypes are consistent with the existence of molecular processes that contribute to the dissociation of non-homologous interactions, which in the model are captured as a repulsive force.

First, when the heterodimeric Hop2/Mnd1 protein complex is defective, chromosomes undergo stable non-homologous synapsis in yeast, mammals and plants [34, 5861]. The Hop2/Mnd1 complex also mediates homologous strand exchange [62], yet non-homologous chromosome associations are unlikely to be an indirect consequence of strand exchange as non-homologous synapsis is not a shared feature of mutants defective for strand exchange [34, 58]. Importantly, involvement of a single protein complex in homologous strand exchange and dissociation of non-homologous chromosome interactions is consistent with both forces acting upon the same molecular intermediate, most likely displacement loops between DNA segments with extensive or limited sequence similarity, respectively. Second, mutation of the Ph1 locus in allopolyploid wheat results in erroneous stabilization of interactions between non-homologous chromosomes, suggesting that Ph1 normally mediates dissociation of such interactions [63]. Third, the mismatch repair machinery disrupts recombination between DNA segments with limited sequence similarity via ejection of the invading strand [30, 34, 36, 64]. In our model, this type of heteroduplex rejection is simulated by a repulsive force that minimizes the association time of non-homologous chromosome partners.

Repulsive interactions may also be the underlying cause for intermittent separation of paired homologs observed in live-cell imaging studies [18]. Such separations have been interpreted as spatially restricted sub-diffusion involving fully paired homologs [18]. Movie animations of our model indicate that similar disruptions could arise due to collisions with non-homologous chromosomes during incipient pairing interactions of homolog pairs (e.g., see [45]). Importantly, results from our model suggest that during this exploratory phase, homolog pairs are susceptible to becoming dislodged due to collisions with non-homologous chromosomes.

The agent-based modeling framework developed here is new to chromosome dynamics but has been useful in modeling other multicomponent systems that involve attractive and repulsive forces, ranging from molecular to macroscopic components [20, 21, 26, 6575]. Earlier modeling approaches have focused on interactions between individual homologs which likely play an important role in completion of the pairing process, but they have not considered the effects of interactions between non-homologous chromosomes [18, 2022].

Our model has identified unexpected threshold effects for several parameters where minor changes result in major nonlinear effects. Accordingly, for the current settings, homolog pairing levels and/or kinetics are disproportionately increased when chromosomes are initially separated by less than 4 μm (Fig 4), when chromosome size exceeds 400 kbp (Fig 5), and when chromosome velocity is increased from 210 nm/s to 240 nm/s (Fig 6). Such discontinuities are likely related to a critical threshold of non-homologous chromosome encounters that needs to be crossed for homologous chromosomes to become confined to the same nuclear areas, thereby facilitating homolog encounters. Accordingly, the frequency of non-homologous interactions is reduced and crowding of homologous chromosomes into the same nuclear area occurs at low frequencies when chromosomes are present in lower numbers, exhibit smaller sizes or fail to achieve unidirectional movement due to disturbance by Brownian motion.

Rapid chromosome movements during prophase of meiosis I have been observed in species from yeast to mouse [1517, 76]. Intriguingly, in different taxa, chromosome movements are mediated by different cytoskeleton components potentially resulting in a wide range of movement speeds [17]. Our model indicates that chromosome movements must occur above a certain velocity threshold, likely determined by the number of chromosomes and the dimensions of the nucleus, providing a potential reason why chromosomes move at distinct speeds in different organisms. Moreover, our model predicts that the three smallest yeast chromosomes would be slower in completing homolog pairing (Fig 5). Notably, the same chromosomes exhibit increased DSB and crossover frequencies [53, 54], features that may specifically compensate for size-related disadvantages in pairing [81].

A key feature of the current agent-based model is that the entire ensemble of participating entities is included in the simulations, capturing both typical behavior and deviations thereof. In contrast, experimental analysis of homolog pairing is limited by the availability of distinct tags for individual chromosomes. Accordingly, observations from a small number of homolog pairs have been extrapolated to the entire chromosome complement. Moreover, in cases where the pairing status is monitored in fixated cells, the progression of pairing must be inferred from different cells retrieved from the culture at different times. In contrast, our model captures pairing efficiencies and kinetics of all chromosomes in the same cell over time. This has already allowed us to predict different sensitivities to translation velocity thresholds of small, medium-sized, and very large chromosomes that would have eluded a population-based approach (see Fig 6B).

One of the key features of our modeling framework is that it is easy to explore various aspects of the parameter space. The model is readily scaled to different biological settings and may provide predictions for a multitude of cellular scenarios. For example, meiosis in the Indian muntjac involves only three very large homolog pairs [77], whereas in some insects and plants between 200 [78] and 600 homolog pairs need to complete pairing [79]. Such complexities are inaccessible to current experimentation but become analyzable by the current model. In future work, our model could be used to explore, e.g., how pairing dynamics are affected by different nucleus sizes and shapes, or by a nucleus represented by a 3-dimensional volume rather than a 2 dimensional surface area. As an example, we have already started exploring the effect of nucleus size on pairing dynamics, indicating that a smaller nuclear area accelerates pairing likely because each chromosome interacts with several non-homologous neighbors at the same time. Conversely, pairing becomes inefficient when the nuclear area is increased above a certain size (e.g., S4 Fig). Other extensions of our model might include chromosome-size dependent velocities, and temporary changes in effective nuclear volume, as in the case of directed chromosome movements during the horsetail stage in S. pombe where the entire chromosome complement temporarily becomes confined to small regions within the nuclear volume [82].

In summary, the modeling approach developed here suggests that homolog pairing is achieved by two mechanistically distinct, yet temporally coinciding processes: Homologs become confined to a nuclear area due to dissociation of interactions with the entire non-homologous chromosome set achieved via a repulsive force. Confinement to smaller areas enables homologs to assess similarities modeled as an attractive force. Importantly, both types of interactions involve close-range physical DNA interactions. Our model makes specific predictions about contributions of chromosome dimensions and movement velocity in combination with chromosome numbers that may further be affected by nuclear dimension. Chromosome numbers, nuclear dimensions, and movement speeds vary widely among different organisms and may affect pairing requirements. Some of these parameters are accessible to experimental manipulation [80], rendering predictions from our model testable in appropriate experiments.

Supporting information

S1 Appendix. Collection of additional supporting information.

This Appendix contains additional details and biological discussion regarding the supplemental figures to this manuscript.

https://doi.org/10.1371/journal.pcbi.1011416.s001

(PDF)

S1 Fig. Effects of chromosome velocity on pairing kinetics.

(A) Reproduction of Fig 6a, each dot indicates the average pairing levels of all 32 chromosomes of true sizes at the indicated time points. Black (300 nm/s) indicates pairing levels at the velocity of chromosome movements in the model in Figs 3 and 5 Chromosomes fail to pair at velocities around 150 nm/s, likely due to the effect of thermal noise. Increases in 30 nm/s increments increases pairing efficiencies at t = 9h ∼ 3-fold, with more modest gains above 240 nm/s where essentially all 16 homologs pair efficiently. (B) Results from S1A Fig were normalized by maximum pairing levels. With increased movement velocities, 50% pairing levels are achieved at progressively earlier time points.

https://doi.org/10.1371/journal.pcbi.1011416.s002

(EPS)

S2 Fig. Effects of movement velocity on homologous pairing kinetics.

All results for 16 homolog pairs using true chromosome lengths. Chromosomes are arranged according to their initial distance to facilitate comparison with the experimental data set. Chromosome movement velocity is (A) 180 nm/s, (B) 210 nm/s, (C) 240 nm/s, and (D) 270 nm/s. Note the sharp transition between 210 nm/s and 240 nm/s, as also shown in Fig 6A. The pairing distance is highlighted in gray at 400 nm. (200 realizations, error bars indicate SD).

https://doi.org/10.1371/journal.pcbi.1011416.s003

(EPS)

S3 Fig. Non-homologous chromosome distances for wild type and spo11 hypomorph, with experimental observations in the insets.

(A) Reproduction of Fig 4B for non-homologous distances in the wild-type simulations. (B) Non-homologous pairs in spo11 mutants where one of the two homolog partners is matched with a non-homologous chromosome exhibiting an optimally matched initial distance with its cognate homolog partner at t = 3h. Inset in (B) shows experimentally determined distances between non-homologous GFP-tagged chromosomes II and III. For details on the experimental conditions see [13]. The pairing distance is highlighted in gray at 400 nm. *Note the experimental data do not include information for t = 9h. (200 realizations, error bars indicate SD).

https://doi.org/10.1371/journal.pcbi.1011416.s004

(EPS)

S4 Fig. Effects of nucleus size on pairing kinetics and efficiency.

Average distances between homolog pairs for actual sizes of yeast chromosomes as a function of nucleus size. Observe that pairing occurs much less frequently as the nucleus radius increases consistent with the roles of confinement and repulsive non-homologous interactions in pairing kinetics. In a high-density environment of a small nucleus, the repulsive interactions have an even larger effect because each chromosome is interacting with many non-homologous neighbors in close proximity. This drives the pairing of homologs by quickly filling the nucleus space with excluded regions due to repulsion. Results over 200 realizations. Error bars indicate SD.

https://doi.org/10.1371/journal.pcbi.1011416.s005

(EPS)

S5 Fig. Contributions of attractive and repulsive forces on pairing efficiencies and kinetics.

(A) Pairing wild-type model using true chromosome lengths and a standard translational movement velocity of 300 nm/s as primarily studied herein. In (B) the repulsive strength is increased by an order of magnitude Cr = 0.05 and the attractive strength is decreased by an order of magnitude Ca = 0.0005. In (C) the reverse is true Cr = 0.0005 and Ca = 0.05 (200 realizations, error bars indicate SD). The pairing distance is highlighted by a gray rectangle.

https://doi.org/10.1371/journal.pcbi.1011416.s006

(EPS)

S6 Fig. Dynamics of chromosome III.

(A) Chromosome III in each of 200 realizations are sorted by their initial distance at t = 3h to their homologous chromosome partner for comparison with experiment. (B) The average distance between Chromosome III and Chromosome II in each of the 200 realizations sorted by their current distance at four time points. (C) Average homologous (Chromosome III with itself) and non-homologous (Chromosome III with Chromosome II) distances as a function of time.

https://doi.org/10.1371/journal.pcbi.1011416.s007

(EPS)

Acknowledgments

The authors thank Sebastian Sensale Rodriguez for useful discussions.

References

  1. 1. Xu N, Tsai CL, Lee JT. Transient homologous chromosome pairing marks the onset of X inactivation. Science. 2006;311(5764):1149–1152. pmid:16424298
  2. 2. LaSalle JM, Lalande M. Homologous association of oppositely imprinted chromosomal domains. Science. 1996;272(5262):725–728. pmid:8614834
  3. 3. Gruhn JR, Hoffmann ER. Errors of the Egg: The Establishment and Progression of Human Aneuploidy Research in the Maternal Germline. Annual Review of Genetics. 2022;56:369–390. pmid:36055648
  4. 4. Lam I, Keeney S. Mechanism and regulation of meiotic recombination initiation. Cold Spring Harbor Perspectives In Biology. 2015;7(1):a016634.
  5. 5. Barzel A, Kupiec M. Finding a match: how do homologous sequences get together for recombination? Nature Reviews Genetics. 2008;9(1):27–37. pmid:18040271
  6. 6. Zickler D, Kleckner N. Recombination, pairing, and synapsis of homologs during meiosis. Cold Spring Harbor Perspectives In Biology. 2015;7(6):a016626. pmid:25986558
  7. 7. Börner GV, Hochwagen A, MacQueen AJ. Meiosis in budding yeast. Genetics. 2023;225(2):iyad125. pmid:37616582
  8. 8. Sym M, Engebrecht J, Roeder GS. ZIP1 is a synaptonemal complex protein required for meiotic chromosome synapsis. Cell. 1993;72(3):365–378. pmid:7916652
  9. 9. Storlazzi A, Gargano S, Ruprich-Robert G, Falque M, David M, Kleckner N, et al. Recombination proteins mediate meiotic spatial chromosome organization and pairing. Cell. 2010;141(1):94–106. pmid:20371348
  10. 10. Weiner BM, Kleckner N. Chromosome pairing via multiple interstitial interactions before and during meiosis in yeast. Cell. 1994;77(7):977–991. pmid:8020104
  11. 11. Burgess SM, Kleckner N, Weiner BM. Somatic pairing of homologs in budding yeast: existence and modulation. Genes & Development. 1999;13(12):1627–1641. pmid:10385630
  12. 12. Peoples-Holst TL, Burgess SM. Multiple branches of the meiotic recombination pathway contribute independently to homolog pairing and stable juxtaposition during meiosis in budding yeast. Genes & Development. 2005;19(7):863–874. pmid:15805472
  13. 13. Sandhu R, Neria FM, Neria JM, Chen X, Hollingsworth NM, Börner GV. DNA helicase Mph1FANCM ensures meiotic recombination between parental chromosomes by dissociating precocious displacement loops. Developmental Cell. 2020;53(4):458–472. pmid:32386601
  14. 14. Brown MS, Zanders S, Alani E. Sustained and rapid chromosome movements are critical for chromosome pairing and meiotic progression in budding yeast. Genetics. 2011;188(1):21–32.
  15. 15. Conrad MN, Lee CY, Chao G, Shinohara M, Kosaka H, Shinohara A, et al. Rapid telomere movement in meiotic prophase is promoted by NDJ1, MPS3, and CSM4 and is modulated by recombination. Cell. 2008;133(7):1175–1187. pmid:18585352
  16. 16. Koszul R, Kim K, Prentiss M, Kleckner N, Kameoka S. Meiotic chromosomes move by linkage to dynamic actin cables with transduction of force through the nuclear envelope. Cell. 2008;133(7):1188–1201. pmid:18585353
  17. 17. Alleva B, Smolikove S. Moving and stopping: Regulation of chromosome movement to promote meiotic chromosome pairing and synapsis. Nucleus. 2017;8(6):613–624. pmid:28892406
  18. 18. Newman TA, Beltran B, McGehee JM, Elnatan D, Cahoon CK, Paddy MR, et al. Diffusion and distal linkages govern interchromosomal dynamics during meiotic prophase. Proceedings of the National Academy of Sciences. 2022;119(12):e2115883119. pmid:35302885
  19. 19. Brar GA, Hochwagen A, Ee LsS, Amon A. The multiple roles of cohesin in meiotic chromosome morphogenesis and pairing. Molecular Biology of the Cell. 2009;20(3):1030–1047. pmid:19073884
  20. 20. Marshall WF, Fung JC. Modeling meiotic chromosome pairing: nuclear envelope attachment, telomere-led active random motion, and anomalous diffusion. Physical Biology. 2016;13(2):026003. pmid:27046097
  21. 21. Marshall WF, Fung JC. Modeling meiotic chromosome pairing: a tug of war between telomere forces and a pairing-based Brownian ratchet leads to increased pairing fidelity. Physical Biology. 2019;16(4):046005. pmid:30943453
  22. 22. Penfold CA, Brown PE, Lawrence ND, Goldman AS. Modeling meiotic chromosomes indicates a size dependent contribution of telomere clustering and chromosome rigidity to homologue juxtaposition. PLoS Computational Biology. 2012;8(5):e1002496. pmid:22570605
  23. 23. Navarro EJ, Marshall WF, Fung JC. Modeling cell biological features of meiotic chromosome pairing to study interlock resolution. PLoS Computational Biology. 2022;18(6):e1010252. pmid:35696428
  24. 24. Dorninger D, Karigl G, Loidl J. Simulation of chromosomal homology searching in meiotic pairing. Journal of Theoretical Biology. 1995;176(2):247–260. pmid:7475113
  25. 25. Ryan SD, Berlyand L, Haines BM, Karpeev D. A kinetic model for semidilute bacterial suspensions. Multiscale Modeling & Simulation. 2013;11(4):1176–1196.
  26. 26. Ryan SD. Role of hydrodynamic interactions in chemotaxis of bacterial populations. Physical Biology. 2019;17(1):016003. pmid:31726435
  27. 27. Spichal M, Fabre E. The emerging role of the cytoskeleton in chromosome dynamics. Frontiers in Genetics. 2017;8:60. pmid:28580009
  28. 28. D’Orsogna MR, Chuang YL, Bertozzi AL, Chayes LS. Self-propelled particles with soft-core interactions: patterns, stability, and collapse. Physical Review Letters. 2006;96(10):104302. pmid:16605738
  29. 29. Qi Z, Redding S, Lee JY, Gibb B, Kwon YH, Niu H, Gaines WA, Sung P, Greene EC. DNA sequence alignment by microhomology sampling during homologous recombination. Cell. 2015; 160(5):856–869 pmid:25684365
  30. 30. Danilowicz C, Hermans L, Coljee V, Prévost C, Prentiss M. ATP hydrolysis provides functions that promote rejection of pairings between different copies of long repeated sequences. Nucleic Acids Research. 2017; 45(14):8448–8462 pmid:28854739
  31. 31. Forget AL, Kowalczykowski SC. Single-molecule imaging of DNA pairing by RecA reveals a three-dimensional homology search. Nature. 2012; 482(7385): 423–427 pmid:22318518
  32. 32. Kornyshev AA, Leikin S. Sequence recognition in the pairing of DNA duplexes. Physical Review Letters. 2001; 86(16): 3666 pmid:11328049
  33. 33. Gebala M, Johnson SL, Narlikar GJ, Herschlag D. Ion counting demonstrates a high electrostatic field generated by the nucleosome eLife. 2019; 8: e44993 pmid:31184587
  34. 34. Tsubouchi H, Roeder GS. The importance of genetic recombination for fidelity of chromosome pairing in meiosis. Developmental Cell. 2003;5(6):915–925. pmid:14667413
  35. 35. Tavares EM, Wright WD, Heyer W-D, Le Cam E, Dupaigne P. In vitro role of Rad54 in Rad51-ssDNA filament-dependent homology search and synaptic complexes formation Nature Communications. 2019; 10(1): 4058 pmid:31492866
  36. 36. Anand R, Beach A, Li K, Haber J. Rad51-mediated double-strand break repair and mismatch correction of divergent substrates Nature. 2017; 544(7650): 377–380 pmid:28405019
  37. 37. Hunter N, Chambers SR, Louis EJ, Borts RH. The mismatch repair system contributes to meiotic sterility in an interspecific yeast hybrid. The EMBO Journal. 1996;15(7):1726–1733. pmid:8612597
  38. 38. Ariel G, Sidortsov M, Ryan SD, Heidenreich S, Bär M, Be’Er A. Collective dynamics of two-dimensional swimming bacteria: Experiments and models. Physical Review E. 2018;98(3):032415.
  39. 39. Wang F, Yuan J, Berg HC. Switching dynamics of the bacterial flagellar motor near zero load. Proceedings of the National Academy of Sciences. 2014;111(44):15752–15755. pmid:25331864
  40. 40. Brosseau Q, Usabiaga FB, Lushi E, Wu Y, Ristroph L, Zhang J, et al. Relating rheotaxis and hydrodynamic actuation using asymmetric gold-platinum phoretic rods. Physical Review Letters. 2019;123(17):178004. pmid:31702241
  41. 41. Drescher K, Dunkel J, Cisneros LH, Ganguly S, Goldstein RE. Fluid dynamics and noise in bacterial cell–cell and cell–surface scattering. Proceedings of the National Academy of Sciences. 2011;108(27):10940–10945. pmid:21690349
  42. 42. Issa MW, Baumgartner NR, Kalil MA, Ryan SD, Wirth CL. Charged nanoparticles quench the propulsion of active janus colloids. ACS Omega. 2019;4(8):13034–13041. pmid:31460430
  43. 43. Chepizhko O, Saintillan D, Peruani F. Revisiting the emergence of order in active matter. Soft Matter. 2021;17(11):3113–3120. pmid:33599237
  44. 44. Lawley SD, Miles CE. Diffusive search for diffusing targets with fluctuating diffusivity and gating. Journal of Nonlinear Science. 2019;29(6):2955–2985.
  45. 45. Ryan SD. Zenodo Repository: Supplemental Movie Files for “Agent-Based Modeling of a Nuclear Chromosome Ensemble Identifies Determinants of Homolog Pairing During Meiosis” by Chriss et al.; 2024. Available from: https://zenodo.org/records/10246589.
  46. 46. Smith SB, Cui Y, Bustamante C. Overstretching B-DNA: the elastic response of individual double-stranded and single-stranded DNA molecules. Science. 1996;271(5250):795–799. pmid:8628994
  47. 47. Mimitou EP, Yamada S, Keeney S. A global view of meiotic double-strand break end resection. Science. 2017;355(6320):40–45. pmid:28059759
  48. 48. Zakharyevich K, Ma Y, Tang S, Hwang PY-H, Boiteux S, Hunter N. Temporally and biochemically distinct activities of Exo1 during meiosis: Double-Strand Break resection and resolution of double Holliday junctions Molecular Cell. 2010; 40(6): 1001–1015 pmid:21172664
  49. 49. Lee CY, Conrad MN, Dresser ME. Meiotic chromosome pairing is promoted by telomere-led chromosome movements independent of bouquet formation. PLoS Genetics. 2012;8(5):e1002730. pmid:22654677
  50. 50. Marshall W, Straight A, Marko J, Swedlow J, Dernburg A, Belmont A, et al. Interphase chromosomes undergo constrained diffusional motion in living cells. Current Biology. 1997;7(12):930–939. pmid:9382846
  51. 51. MATLAB. version 9.11.0.1809720 (R2021b). Natick, Massachusetts: The MathWorks Inc.; 2021.
  52. 52. Ryan SD. GitHub Repository: Chromosome Dynamics in Prophase I; 2024. Available from: https://github.com/sdryan/ChromosomeDynamicsProphase1.
  53. 53. Murakami H, Lam I, Huang PC, Song J, van Overbeek M, Keeney S. Multilayered mechanisms ensure that short chromosomes recombine in meiosis. Nature. 2020;582(7810):124–128. pmid:32494071
  54. 54. Kaback DB, Guacci V, Barber D, Mahon JW. Chromosome size-dependent control of meiotic recombination. Science. 1992;256(5054):228–232. pmid:1566070
  55. 55. Therizols P, Duong T, Dujon B, Zimmer C, Fabre E. Chromosome arm length and nuclear constraints determine the dynamic relationship of yeast subtelomeres. Proceedings of the National Academy of Sciences. 2010;107(5):2025–2030. pmid:20080699
  56. 56. Winkler RG. Dynamics of flexible active Brownian dumbbells in the absence and the presence of shear flow. Soft Matter. 2016;12(16):3737–3749 pmid:26980630
  57. 57. Neale MJ, Ramachandran M, Trelles-Sticken E, Scherthan H, Goldman AS. Wild-type levels of Spo11-induced DSBs are required for normal single-strand resection during meiosis. Molecular Cell. 2002;9(4):835–846. pmid:11983174
  58. 58. Leu JY, Chua PR, Roeder GS. The meiosis-specific Hop2 protein of S. cerevisiae ensures synapsis between homologous chromosomes. Cell. 1998;94(3):375–386. pmid:9708739
  59. 59. Zierhut C, Berlinger M, Rupp C, Shinohara A, Klein F. Mnd1 is required for meiotic interhomolog repair. Current Biology. 2004;14(9):752–762. pmid:15120066
  60. 60. Petukhova GV, Romanienko PJ, Camerini-Otero RD. The Hop2 protein has a direct role in promoting interhomolog interactions during mouse meiosis. Developmental Cell. 2003;5(6):927–936. pmid:14667414
  61. 61. Kerzendorfer C, Vignard J, Pedrosa-Harand A, Siwiec T, Akimcheva S, Jolivet S, et al. The Arabidopsis thaliana MND1 homologue plays a key role in meiotic homologous pairing, synapsis and recombination. Journal of Cell Science. 2006;119(12):2486–2496. pmid:16763194
  62. 62. Petukhova GV, Pezza RJ, Vanevski F, Ploquin M, Masson JY, Camerini-Otero RD. The Hop2 and Mnd1 proteins act in concert with Rad51 and Dmc1 in meiotic recombination. Nature Structural & Molecular Biology. 2005;12(5):449–453. pmid:15834424
  63. 63. Griffiths S, Sharp R, Foote TN, Bertin I, Wanous M, Reader S, et al. Molecular characterization of Ph1 as a major chromosome pairing locus in polyploid wheat. Nature. 2006;439(7077):749–752. pmid:16467840
  64. 64. Bozdag GO, Ono J, Denton JA, Karakoc E, Hunter N, Leu JY, et al. Breaking a species barrier by enabling hybrid recombination. Current Biology. 2021;31(4):R180–R181. pmid:33621502
  65. 65. Wong H, Arbona JM, Zimmer C. How to build a yeast nucleus. Nucleus. 2013;4(5):361–366. pmid:23974728
  66. 66. Marko JF, Siggia ED. Polymer models of meiotic and mitotic chromosomes. Molecular Biology of the Cell. 1997;8(11):2217–2231. pmid:9362064
  67. 67. Cheng TM, Heeger S, Chaleil RA, Matthews N, Stewart A, Wright J, et al. A simple biophysical model emulates budding yeast chromosome condensation. eLife. 2015;4:e05565. pmid:25922992
  68. 68. Ariel G, Ayali A. Locust collective motion and its modeling. PLoS Computational Biology. 2015;11(12):e1004522. pmid:26656851
  69. 69. Boissard E, Degond P, Motsch S. Trail formation based on directed pheromone deposition. Journal of Mathematical Biology. 2013;66:1267–1301. pmid:22526837
  70. 70. Ryan SD. A model for collective dynamics in ant raids. Journal of Mathematical Biology. 2016;72:1579–1606. pmid:26304617
  71. 71. Baumgartner NR II, Ryan SD. Interaction of red crabs with yellow crazy ants during migration on Christmas Island. Mathematical Biosciences. 2020;330:108486.
  72. 72. Ryan SD, McCarthy Z, Potomkin M. Motor protein transport along inhomogeneous microtubules. Bulletin of Mathematical Biology. 2021;83:1–29. pmid:33415532
  73. 73. Ryan SD, Haines BM, Berlyand L, Ziebert F, Aranson IS. Viscosity of bacterial suspensions: Hydrodynamic interactions and self-induced noise. Physical Review E. 2011;83(5):050904. pmid:21728480
  74. 74. Ryan SD, Ariel G, Be’er A. Anomalous fluctuations in the orientation and velocity of swarming bacteria. Biophysical Journal. 2016;111(1):247–255. pmid:27410751
  75. 75. Peled S, Ryan SD, Heidenreich S, Bär M, Ariel G, Be’Er A. Heterogeneous bacterial swarms with mixed lengths. Physical Review E. 2021;103(3):032413. pmid:33862716
  76. 76. Scherthan H, Weich S, Schwegler H, Heyting C, Härle M, Cremer T. Centromere and telomere movements during early meiotic prophase of mouse and man are associated with the onset of chromosome pairing. The Journal of Cell Biology. 1996;134(5):1109–1125. pmid:8794855
  77. 77. Wurster DH, Benirschke K. Indian muntjac, Muntiacus muntjak: a deer with a low diploid chromosome number. Science, 168 (1970), pp. 1364–1366 pmid:5444269
  78. 78. Lukhtanov VA. The blue butterfly Polyommatus (Plebicula) atlanticus (Lepidoptera, Lycaenidae) holds the record of the highest number of chromosomes in the non-polyploid eukaryotic organisms. Comparative Cytogenetics. 2015;9(4):683. pmid:26753083
  79. 79. Ninan CA. Studies on the cytology and phylogeny of the Pteridophytes. VI. Observations on the Ophioglossaceae. Cytologia. 1958;23:291–316.
  80. 80. Richardson SM, Mitchell LA, Stracquadanio G, Yang K, Dymond JS, DiCarlo JE, et al. Design of a synthetic yeast genome. Science. 2017;355(6329):1040–1044. pmid:28280199
  81. 81. Rosin L. F., Gil J. Jr, Drinnenberg I. A., & Lei E. P. (2021). Oligopaint DNA FISH reveals telomere-based meiotic pairing dynamics in the silkworm, Bombyx mori. PLoS genetics. 2021, 17(7): e1009700.
  82. 82. Chikashige Y., Ding D. Q., Funabiki H., Haraguchi T., Mashiko S., Yanagida M., & Hiraoka Y. Telomere-led premeiotic chromosome movement in fission yeast. Science. 1994, 264(5156): 270–273.