Single-molecule studies of origin licensing reveal mechanisms ensuring bidirectional helicase loading.

Loading of the ring-shaped Mcm2-7 replicative helicase around DNA licenses eukaryotic origins of replication. During loading, Cdc6, Cdt1, and the origin-recognition complex (ORC) assemble two heterohexameric Mcm2-7 complexes into a head-to-head double hexamer that facilitates bidirectional replication initiation. Using multi-wavelength single-molecule fluorescence to monitor the events of helicase loading, we demonstrate that double-hexamer formation is the result of sequential loading of individual Mcm2-7 complexes. Loading of each Mcm2-7 molecule involves the ordered association and dissociation of distinct Cdc6 and Cdt1 proteins. In contrast, one ORC molecule directs loading of both helicases in each double hexamer. Based on single-molecule FRET, arrival of the second Mcm2-7 results in rapid double-hexamer formation that anticipates Cdc6 and Cdt1 release, suggesting that Mcm-Mcm interactions recruit the second helicase. Our findings reveal the complex protein dynamics that coordinate helicase loading and indicate that distinct mechanisms load the oppositely oriented helicases that are central to bidirectional replication initiation.

Eukaryotic origin-dependent DNA replication in vitro reveals sequential action of DDK and S-CDK kinases.

Proper eukaryotic DNA replication requires temporal separation of helicase loading from helicase activation and replisome assembly. Using an in vitro assay for eukaryotic origin-dependent replication initiation, we investigated the control of these events. After helicase loading, we found that the Dbf4-dependent Cdc7 kinase (DDK) but not S phase cyclin-dependent kinase (S-CDK) is required for the initial origin recruitment of Sld3 and the Cdc45 helicase-activating protein. Likewise, in vivo, DDK drives early-firing-origin recruitment of Cdc45 before activation of S-CDK. After S-CDK activation, a second helicase-activating protein (GINS) and the remainder of the replisome are recruited to the origin. Finally, recruitment of lagging but not leading strand DNA polymerases depends on Mcm10 and DNA unwinding. Our studies identify distinct roles for DDK and S-CDK during helicase activation and support a model in which the leading strand DNA polymerase is recruited prior to origin DNA unwinding and RNA primer synthesis.

Regulation of replication origin licensing by ORC phosphorylation reveals a two-step mechanism for Mcm2-7 ring closing.

Eukaryotic DNA replication must occur exactly once per cell cycle to maintain cell ploidy. This outcome is ensured by temporally separating replicative helicase loading (G1 phase) and activation (S phase). In budding yeast, helicase loading is prevented outside of G1 by cyclin-dependent kinase (CDK) phosphorylation of three helicase-loading proteins: Cdc6, the Mcm2-7 helicase, and the origin recognition complex (ORC). CDK inhibition of Cdc6 and Mcm2-7 is well understood. Here we use single-molecule assays for multiple events during origin licensing to determine how CDK phosphorylation of ORC suppresses helicase loading. We find that phosphorylated ORC recruits a first Mcm2-7 to origins but prevents second Mcm2-7 recruitment. The phosphorylation of the Orc6, but not of the Orc2 subunit, increases the fraction of first Mcm2-7 recruitment events that are unsuccessful due to the rapid and simultaneous release of the helicase and its associated Cdt1 helicase-loading protein. Real-time monitoring of first Mcm2-7 ring closing reveals that either Orc2 or Orc6 phosphorylation prevents Mcm2-7 from stably encircling origin DNA. Consequently, we assessed formation of the MO complex, an intermediate that requires the closed-ring form of Mcm2-7. We found that ORC phosphorylation fully inhibits MO complex formation and we provide evidence that this event is required for stable closing of the first Mcm2-7. Our studies show that multiple steps of helicase loading are impacted by ORC phosphorylation and reveal that closing of the first Mcm2-7 ring is a two-step process started by Cdt1 release and completed by MO complex formation.

Changing protein-DNA interactions promote ORC binding-site exchange during replication origin licensing.

During origin licensing, the eukaryotic replicative helicase Mcm2-7 forms head-to-head double hexamers to prime origins for bidirectional replication. Recent single-molecule and structural studies revealed that one molecule of the helicase loader ORC (origin recognition complex) can sequentially load two Mcm2-7 hexamers to ensure proper head-to-head helicase alignment. To perform this task, ORC must release from its initial high-affinity DNA-binding site and "flip" to bind a weaker, inverted DNA site. However, the mechanism of this binding-site switch remains unclear. In this study, we used single-molecule Förster resonance energy transfer to study the changing interactions between DNA and ORC or Mcm2-7. We found that the loss of DNA bending that occurs during DNA deposition into the Mcm2-7 central channel increases the rate of ORC dissociation from DNA. Further studies revealed temporally controlled DNA sliding of helicase-loading intermediates and that the first sliding complex includes ORC, Mcm2-7, and Cdt1. We demonstrate that sequential events of DNA unbending, Cdc6 release, and sliding lead to a stepwise decrease in ORC stability on DNA, facilitating ORC dissociation from its strong binding site during site switching. In addition, the controlled sliding we observed provides insight into how ORC accesses secondary DNA-binding sites at different locations relative to the initial binding site. Our study highlights the importance of dynamic protein-DNA interactions in the loading of two oppositely oriented Mcm2-7 helicases to ensure bidirectional DNA replication.

Distinct RPA functions promote eukaryotic DNA replication initiation and elongation.

Replication protein A (RPA) binds single-stranded DNA (ssDNA) and serves critical functions in eukaryotic DNA replication, the DNA damage response, and DNA repair. During DNA replication, RPA is required for extended origin DNA unwinding and DNA synthesis. To determine the requirements for RPA during these processes, we tested ssDNA-binding proteins (SSBs) from different domains of life in reconstituted Saccharomyces cerevisiae origin unwinding and DNA replication reactions. Interestingly, Escherichia coli SSB, but not T4 bacteriophage Gp32, fully substitutes for RPA in promoting origin DNA unwinding. Using RPA mutants, we demonstrated that specific ssDNA-binding properties of RPA are required for origin unwinding but that its protein-interaction domains are dispensable. In contrast, we found that each of these auxiliary RPA domains have distinct functions at the eukaryotic replication fork. The Rfa1 OB-F domain negatively regulates lagging-strand synthesis, while the Rfa2 winged-helix domain stimulates nascent strand initiation. Together, our findings reveal a requirement for specific modes of ssDNA binding in the transition to extensive origin DNA unwinding and identify RPA domains that differentially impact replication fork function.

Mcm10 regulates DNA replication elongation by stimulating the CMG replicative helicase.

Activation of the Mcm2-7 replicative DNA helicase is the committed step in eukaryotic DNA replication initiation. Although Mcm2-7 activation requires binding of the helicase-activating proteins Cdc45 and GINS (forming the CMG complex), an additional protein, Mcm10, drives initial origin DNA unwinding by an unknown mechanism. We show that Mcm10 binds a conserved motif located between the oligonucleotide/oligosaccharide fold (OB-fold) and A subdomain of Mcm2. Although buried in the interface between these domains in Mcm2-7 structures, mutations predicted to separate the domains and expose this motif restore growth to conditional-lethal MCM10 mutant cells. We found that, in addition to stimulating initial DNA unwinding, Mcm10 stabilizes Cdc45 and GINS association with Mcm2-7 and stimulates replication elongation in vivo and in vitro. Furthermore, we identified a lethal allele of MCM10 that stimulates initial DNA unwinding but is defective in replication elongation and CMG binding. Our findings expand the roles of Mcm10 during DNA replication and suggest a new model for Mcm10 function as an activator of the CMG complex throughout DNA replication.

Putting two heads together to unwind DNA.

The loading of replicative helicases onto DNA is tightly regulated in all organisms, yet the molecular mechanisms for this event remain poorly defined. Remus et al. (2009) provide important insights into helicase loading in eukaryotes, showing that the Mcm2-7 replicative helicase encircles double-stranded DNA as head-to-head double hexamers.

The dynamics of eukaryotic replication initiation: origin specificity, licensing, and firing at the single-molecule level.

Eukaryotic replication initiation is highly regulated and dynamic. It begins with the origin recognition complex (ORC) binding DNA sites called origins of replication. ORC, together with Cdc6 and Cdt1, mediate pre-replicative complex (pre-RC) assembly by loading a double hexamer of Mcm2-7: the core of the replicative helicase. Here, we use single-molecule imaging to directly visualize Saccharomyces cerevisiae pre-RC assembly and replisome firing in real time. We show that ORC can locate and stably bind origins within large tracts of non-origin DNA and that Cdc6 drives ordered pre-RC assembly. We further show that the dynamics of the ORC-Cdc6 interaction dictate Mcm2-7 loading specificity and that Mcm2-7 double hexamers form preferentially at a native origin sequence. Finally, we demonstrate that single Mcm2-7 hexamers propagate bidirectionally, monotonically, and processively as constituents of active replisomes.

Multiple functions for Mcm2-7 ATPase motifs during replication initiation.

The Mcm2-7 replicative helicase is central to all steps of eukaryotic DNA replication. The hexameric ring of Mcm subunits forms six essential ATPases whose contributions to replication initiation remain unclear. Mcm2-7 complexes containing ATPase-motif mutations showed Mcm2-7 ATP binding and hydrolysis are required for helicase loading. Loading-defective Mcm2-7 mutant complexes were defective in initial Mcm2-7 recruitment or Cdt1 release. Comparison with Cdc6 ATPase mutants showed that Cdc6 ATP hydrolysis is not required for helicase loading but instead drives removal of Mcm2-7 complexes that cannot complete loading. A subset of Mcm2-7 ATPase-site mutants completed helicase loading but could not initiate replication. Individual mutants were defective in distinct events during helicase activation, including maintenance of DNA association, recruitment of the GINS helicase activator, and DNA unwinding. Consistent with its heterohexameric structure, our findings show that the six Mcm2-7 ATPase active sites are specialized for different functions during helicase loading and activation.

Transcriptional repression of CDC6 and SLD2 during meiosis is associated with production of short heterogeneous RNA isoforms.

Execution of the meiotic and mitotic cell division programs requires distinct gene expression patterns. Unlike mitotic cells, meiotic cells reduce ploidy by following one round of DNA replication with two rounds of chromosome segregation (meiosis I and meiosis II). However, the mechanisms by which cells prevent DNA replication between meiosis I and meiosis II are not fully understood. Here, we show that transcriptional repression of two essential DNA replication genes, CDC6 and SLD2, is associated with production of shorter meiosis-specific RNAs containing the 3' end of both genes. Despite the short CDC6 RNA coding for a short protein (Cdc6short), this protein is not essential for meiosis and it does not have either a positive or negative impact on DNA replication. Production of CDC6short mRNA does not require the upstream CDC6 promoter (PCDC6) and is not a processed form of the full-length RNA. Instead, CDC6short depends on transcription initiation from within the ORF upon repression of PCDC6. Finally, using CDC6 genes from related yeast, we show that repression of full-length CDC6 mRNA is evolutionarily conserved and that this repression is consistently associated with production of unique short CDC6 RNAs. Together, these data demonstrate that meiotic cells transcriptionally repress full-length CDC6 and SLD2, and that inactivation of PCDC6 results in heterogeneous transcription initiation from within the CDC6 ORF.

CDK prevents Mcm2-7 helicase loading by inhibiting Cdt1 interaction with Orc6.

In Saccharomyces cerevisiae cells, B-type cyclin-dependent kinases (CDKs) target two origin recognition complex (ORC) subunits, Orc2 and Orc6, to inhibit helicase loading. We show that helicase loading by ORC is inhibited by two distinct CDK-dependent mechanisms. Independent of phosphorylation, binding of CDK to an "RXL" cyclin-binding motif in Orc6 sterically reduces the initial recruitment of the Cdt1/Mcm2-7 complex to ORC. CDK phosphorylation of Orc2 and Orc6 inhibits the same step in helicase loading. This phosphorylation of Orc6 is stimulated by the RXL motif and mediates the bulk of the phosphorylation-dependent inhibition of helicase loading. Direct binding experiments show that CDK phosphorylation specifically blocks one of the two Cdt1-binding sites on Orc6. Consistent with the inactivation of one Cdt1-binding site preventing helicase loading, CDK phosphorylation of ORC causes a twofold reduction of initial Cdt1/Mcm2-7 recruitment but results in nearly complete inhibition of Mcm2-7 loading. Intriguingly, in addition to being a target of both CDK inhibitory mechanisms, the Orc6 RXL/cyclin-binding motif plays a positive role in the initial recruitment of Cdt1/Mcm2-7 to the origin, suggesting that this motif is critical for the switch between active and inhibited ORC function at the G1-to-S-phase transition.

Replication origin-flanking roadblocks reveal origin-licensing dynamics and altered sequence dependence.

In eukaryotes, DNA replication initiates from multiple origins of replication for timely genome duplication. These sites are selected by origin licensing, during which the core enzyme of the eukaryotic DNA replicative helicase, the Mcm2-7 (minichromosome maintenance) complex, is loaded at each origin. This origin licensing requires loading two Mcm2-7 helicases around origin DNA in a head-to-head orientation. Current models suggest that the origin-recognition complex (ORC) and cell-division cycle 6 (Cdc6) proteins recognize and encircle origin DNA and assemble an Mcm2-7 double-hexamer around adjacent double-stranded DNA. To test this model and assess the location of Mcm2-7 initial loading, we placed DNA-protein roadblocks at defined positions adjacent to the essential ORC-binding site within Saccharomyces cerevisiae origin DNA. Roadblocks were made either by covalent cross-linking of the HpaII methyltransferase to DNA or through binding of a transcription activator-like effector (TALE) protein. Contrary to the sites of Mcm2-7 recruitment being precisely defined, only single roadblocks that inhibited ORC-DNA binding showed helicase loading defects. We observed inhibition of helicase loading without inhibition of ORC-DNA binding only when roadblocks were placed on both sides of the origin to restrict sliding of a helicase-loading intermediate. Consistent with a sliding helicase-loading intermediate, when either one of the flanking roadblocks was eliminated, the remaining roadblock had no effect on helicase loading. Interestingly, either origin-flanking nucleosomes or roadblocks resulted in helicase loading being dependent on an additional origin sequence known to be a weaker ORC-DNA-binding site. Together, our findings support a model in which sliding helicase-loading intermediates increase the flexibility of the DNA sequence requirements for origin licensing.

Mec1 is one of multiple kinases that prime the Mcm2-7 helicase for phosphorylation by Cdc7.

Activation of the eukaryotic replicative DNA helicase, the Mcm2-7 complex, requires phosphorylation by Cdc7/Dbf4 (Dbf4-dependent kinase or DDK), which, in turn, depends on prior phosphorylation of Mcm2-7 by an unknown kinase (or kinases). We identified DDK phosphorylation sites on Mcm4 and Mcm6 and found that phosphorylation of either subunit suffices for cell proliferation. Importantly, prior phosphorylation of either S/T-P or S/T-Q motifs on these subunits is required for DDK phosphorylation of Mcm2-7 and for normal S phase passage. Phosphomimetic mutations of DDK target sites bypass both DDK function and mutation of the priming phosphorylation sites. Mrc1 facilitates Mec1 phosphorylation of the S/T-Q motifs of chromatin-bound Mcm2-7 during S phase to activate replication. Genetic interactions between priming site mutations and MRC1 or TOF1 deletion support a role for these modifications in replication fork stability. These findings identify regulatory mechanisms that modulate origin firing and replication fork assembly during cell cycle progression.

Conserved nucleosome positioning defines replication origins.

The origin recognition complex (ORC) specifies replication origin location. The Saccharomyces cerevisiae ORC recognizes the ARS (autonomously replicating sequence) consensus sequence (ACS), but only a subset of potential genomic sites are bound, suggesting other chromosomal features influence ORC binding. Using high-throughput sequencing to map ORC binding and nucleosome positioning, we show that yeast origins are characterized by an asymmetric pattern of positioned nucleosomes flanking the ACS. The origin sequences are sufficient to maintain a nucleosome-free origin; however, ORC is required for the precise positioning of nucleosomes flanking the origin. These findings identify local nucleosomes as an important determinant for origin selection and function.

Incorporation into the prereplicative complex activates the Mcm2-7 helicase for Cdc7-Dbf4 phosphorylation.

The essential S-phase kinase Cdc7-Dbf4 acts at eukaryotic origins of replication to trigger a cascade of protein associations that activate the Mcm2-7 replicative helicase. Also known as Dbf4-dependent kinase (DDK), this kinase preferentially targets chromatin-associated Mcm2-7 complexes that are assembled on the DNA during prereplicative complex (pre-RC) formation. Here we address the mechanisms that control the specificity of DDK action. We show that incorporation of Mcm2-7 into the pre-RC increased the level and changes the specificity of DDK phosphorylation of this complex. In the context of the pre-RC, DDK preferentially targets a conformationally distinct subpopulation of Mcm2-7 complexes that is tightly linked to the origin DNA. This targeting requires DDK to tightly associate with Mcm2-7 complexes in a Dbf4-dependent manner. Importantly, we find that DDK association with and phosphorylation of origin-linked Mcm2-7 complexes require prior phosphorylation of the pre-RC. Our findings provide insights into the mechanisms that ensure that DDK action is spatially and temporally restricted to the origin-bound Mcm2-7 complexes that will drive replication fork movement during S phase and suggest new mechanisms to regulate origin activity.

Multiple Cdt1 molecules act at each origin to load replication-competent Mcm2-7 helicases.

Eukaryotic origins of replication are selected by loading a head-to-head double hexamer of the Mcm2-7 replicative helicase around origin DNA. Cdt1 plays an essential but transient role during this event; however, its mechanism of action is unknown. Through analysis of Cdt1 mutations, we demonstrate that Cdt1 performs multiple functions during helicase loading. The C-terminus of Cdt1 binds Mcm2-7, and this interaction is required for efficient origin recruitment of both proteins. We show that origin recognition complex (ORC) and Cdc6 recruit multiple Cdt1 molecules to the origin during helicase loading, and disruption of this multi-Cdt1 intermediate prevents helicase loading. Although dispensable for loading Mcm2-7 double hexamers that are topologically linked to DNA, the essential N-terminal domain of Cdt1 is required to load Mcm2-7 complexes that are competent for association with the Cdc45 and GINS helicase-activating proteins and replication initiation. Our data support a model in which origin-bound ORC and Cdc6 recruit two Cdt1 molecules to initiate double-hexamer formation prior to helicase loading and demonstrate that Cdt1 influences the replication competence of loaded Mcm2-7 helicases.

Separation of DNA replication from the assembly of break-competent meiotic chromosomes.

The meiotic cell division reduces the chromosome number from diploid to haploid to form gametes for sexual reproduction. Although much progress has been made in understanding meiotic recombination and the two meiotic divisions, the processes leading up to recombination, including the prolonged pre-meiotic S phase (meiS) and the assembly of meiotic chromosome axes, remain poorly defined. We have used genome-wide approaches in Saccharomyces cerevisiae to measure the kinetics of pre-meiotic DNA replication and to investigate the interdependencies between replication and axis formation. We found that replication initiation was delayed for a large number of origins in meiS compared to mitosis and that meiotic cells were far more sensitive to replication inhibition, most likely due to the starvation conditions required for meiotic induction. Moreover, replication initiation was delayed even in the absence of chromosome axes, indicating replication timing is independent of the process of axis assembly. Finally, we found that cells were able to install axis components and initiate recombination on unreplicated DNA. Thus, although pre-meiotic DNA replication and meiotic chromosome axis formation occur concurrently, they are not strictly coupled. The functional separation of these processes reveals a modular method of building meiotic chromosomes and predicts that any crosstalk between these modules must occur through superimposed regulatory mechanisms.

A noncanonical GTPase signaling mechanism controls exit from mitosis in budding yeast.

In the budding yeast Saccharomyces cerevisiae , exit from mitosis is coupled to spindle position to ensure successful genome partitioning between mother and daughter cell. This coupling occurs through a GTPase signaling cascade known as the mitotic exit network (MEN). The MEN senses spindle position via a Ras-like GTPase Tem1 which localizes to the spindle pole bodies (SPBs, yeast equivalent of centrosomes) during anaphase and signals to its effector protein kinase Cdc15. How Tem1 couples the status of spindle position to MEN activation is not fully understood. Here, we show that Cdc15 has a relatively weak preference for Tem1 GTP and Tem1's nucleotide state does not change upon MEN activation. Instead, we find that Tem1's nucleotide cycle establishes a localization-based concentration difference in the cell where only Tem1 GTP is recruited to the SPB, and spindle position regulates the MEN by controlling Tem1 localization. SPB localization of Tem1 primarily functions to promote Tem1-Cdc15 interaction for MEN activation by increasing the effective concentration of Tem1. Consistent with this model, we demonstrate that artificially tethering Tem1 to the SPB or concentrating Tem1 in the cytoplasm with genetically encoded multimeric nanoparticles could bypass the requirement of Tem1 GTP and correct spindle position for MEN activation. This localization/concentration-based GTPase signaling mechanism for Tem1 differs from the canonical Ras-like GTPase signaling paradigm and is likely relevant to other localization-based signaling scenarios.

Elevated rates of force development and MgATP binding in F764L and S532P myosin mutations causing dilated cardiomyopathy.

Dilated cardiomyopathy (DCM) is a disease characterized by dilation of the ventricular chambers and reduced contractile function. We examined the contractile performance of chemically-skinned ventricular strips from two heterozygous murine models of DCM-causing missense mutations of myosin, F764L/+ and S532P/+, in an α-myosin heavy chain (MyHC) background. In Ca(2+)-activated skinned myocardial strips, the maximum developed tension in F764L/+ was only ~50% that of litter-mate controls (+/+). The F764L/+ also exhibited significantly reduced rigor stiffness, loaded shortening velocity and power output. Corresponding indices for S532P/+ strips were not different from controls. Manipulation of MgATP concentration in conjunction with measures of viscoelasticity, which provides estimates of myosin detachment rate 2πc, allowed us to probe the molecular basis of changes in crossbridge kinetics that occur with the myosin mutations. By examining the response of detachment rate to varying MgATP we found the rate of MgADP release was unaffected by the myosin mutations. However, MgATP binding rate was higher in the DCM groups compared to controls (422±109mM(-1)·s(-1) in F764L/+, 483±74mM(-1)·s(-1) in S532P/+ and 303±18mM(-1)·s(-1) in +/+). In addition, the rate constant of force development, 2πb, was significantly higher in DCM groups compared to controls (at 5mM MgATP: 36.9±4.9s(-1) in F764L/+, 32.9±4.5s(-1) in S532P/+ and 18.2±1.7s(-1) in +/+). These results suggest that elevated rates of force development and MgATP binding are features of cardiac myofilament function that underlie the development of DCM.

Identifying cellular mechanisms of zinc-induced relaxation in isolated cardiomyocytes.

We tested several molecular and cellular mechanisms of cardiomyocyte contraction-relaxation function that could account for the reduced systolic and enhanced diastolic function observed with exposure to extracellular Zn(2+). Contraction-relaxation function was monitored in isolated rat and mouse cardiomyocytes maintained at 37°C, stimulated at 2 or 6 Hz, and exposed to 32 µM Zn(2+) or vehicle. Intracellular Zn(2+) detected using FluoZin-3 rose to a concentration of ∼13 nM in 3-5 min. Peak sarcomere shortening was significantly reduced and diastolic sarcomere length was elongated after Zn(2+) exposure. Peak intracellular Ca(2+) detected by Fura-2FF was reduced after Zn(2+) exposure. However, the rate of cytosolic Ca(2+) decline reflecting sarcoplasmic reticulum (SR) Ca(2+)-ATPase (SERCA2a) activity and the rate of Na(+)/Ca(2+) exchanger activity evaluated by rapid Na(+)-induced Ca(2+) efflux were unchanged by Zn(2+) exposure. SR Ca(2+) load evaluated by rapid caffeine exposure was reduced by ∼50%, and L-type calcium channel inward current measured by whole cell patch clamp was reduced by ∼70% in cardiomyocytes exposed to Zn(2+). Furthermore, ryanodine receptor (RyR) S2808 and phospholamban (PLB) S16/T17 were markedly dephosphorylated after perfusing hearts with 50 µM Zn(2+). Maximum tension development and thin-filament Ca(2+) sensitivity in chemically skinned cardiac muscle strips were not affected by Zn(2+) exposure. These findings suggest that Zn(2+) suppresses cardiomyocyte systolic function and enhances relaxation function by lowering systolic and diastolic intracellular Ca(2+) concentrations due to a combination of competitive inhibition of Ca(2+) influx through the L-type calcium channel, reduction of SR Ca(2+) load resulting from phospholamban dephosphorylation, and lowered SR Ca(2+) leak via RyR dephosphorylation. The use of the low-Ca(2+)-affinity Fura-2FF likely prevented the detection of changes in diastolic Ca(2+) and SERCA2a function. Other strategies to detect diastolic Ca(2+) in the presence of Zn(2+) are essential for future work.