Data-Based Optimal Switching and Control With Admissibility Guaranteed Q -Learning.

This article addresses the data-based optimal switching and control codesign for discrete-time nonlinear switched systems via a two-stage approximate dynamic programming (ADP) algorithm. Through offline policy improvement and policy evaluation, the proposed algorithm iteratively determines the optimal hybrid control policy using system input/output data. Moreover, a strict proof of the convergence is given for the two-stage ADP algorithm. Admissibility, an essential property of the hybrid control policy must be ensured for practical application. To this end, the properties of the hybrid control policies are analyzed and an admissibility criterion is obtained. To realize the proposed Q -learning algorithm, an actor-critic neural network (NN) structure that employs multiple NNs to approximate the Q -functions and control policies for different subsystems is adopted. By applying the proposed admissibility criterion, the obtained hybrid control policy is guaranteed to be admissible. Finally, two numerical simulations verify the effectiveness of the proposed algorithm.

Novel Evaluation Method for Flushing Efficiency Based on the Principle of Wall Shear Rate Equality under High Temperature and High Pressure.

An excellent flushing efficiency evaluation device is indispensable for studying the cementing flushing fluid. However, the existing flushing efficiency evaluation device cannot simulate the formation and flushing process of the downhole mud cake. Therefore, this paper proposes a novel flushing efficiency evaluation device that can simulate the formation and flushing of downhole mud cakes. The rotational speed of the device during flushing and the rotor's diameter is deduced based on the principle of equal wall shear rate. The evaluation device and method can be used to quantitatively evaluate the flushing efficiency of the flushing liquid on the mud cake under high temperatures and high pressures. This device analyzed the effects of the annular gap, temperature, construction displacement, and flushing fluids on flushing efficiency. The results show that the smaller the annular gap, the higher the temperature, the larger the displacement, the higher the scouring efficiency, and the higher the shear bond strength. Fiber flushing has the dual function of mechanical and chemical flushing, so its flushing efficiency is 14.75% higher than that of heavy flushing fluid at 10 min. The surfactant in the emulsified rinse leads to a sudden increase in rinse efficiency in the middle and late stages. The reduced flushing efficiency of the flushing fluid is due to the reduced ζ potential.

Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.).

Powdery mildew (PM), caused by the fungus Oidium lini in flax, can cause defoliation and reduce seed yield and quality. To date, one major dominant gene (Pm1) and three quantitative trait loci (QTL) on chromosomes 1, 7 and 9 have been reported for PM resistance. To fully dissect the genetic architecture of PM resistance and identify QTL, a diverse flax core collection of 372 accessions augmented with an additional 75 breeding lines were sequenced, and PM resistance was evaluated in the field for eight years (2010-2017) in Morden, Manitoba, Canada. Genome-wide association studies (GWAS) were performed using two single-locus and seven multi-locus statistical models with 247,160 single nucleotide polymorphisms (SNPs) and the phenotypes of the 447 individuals for each year separately as well as the means over years. A total of 349 quantitative trait nucleotides (QTNs) were identified, of which 44 large-effect QTNs (R2 = 10-30%) were highly stable over years. The total number of favourable alleles per accession was significantly correlated with PM resistance (r = 0.74), and genomic selection (GS) models using all identified QTNs generated significantly higher predictive ability (r = 0.93) than those constructed using the 247,160 genome-wide random SNP (r = 0.69), validating the overall reliability of the QTNs and showing the additivity of PM resistance in flax. The QTNs were clustered on the distal ends of all 15 chromosomes, especially on chromosome 5 (0.4-5.6 Mb and 9.4-16.9 Mb) and 13 (4.7-5.2 Mb). To identify candidate genes, a dataset of 3230 SNPs located in resistance gene analogues (RGAs) was used as input for GWAS, from which an additional 39 RGA-specific QTNs were identified. Overall, 269 QTN loci harboured 445 RGAs within the 200 Kb regions spanning the QTNs, including 45 QTNs located within the RGAs. These RGAs supported by significant QTN/SNP allele effects were mostly nucleotide binding site and leucine-rich repeat receptors (NLRs) belonging to either coiled-coil (CC) NLR (CNL) or toll interleukin-1 (TIR) NLR (TNL), receptor-like kinase (RLK), receptor-like protein kinase (RLP), transmembrane-coiled-coil (TM-CC), WRKY, and mildew locus O (MLO) genes. These results constitute an important genomic tool for resistance breeding and gene cloning for PM in flax.

Aegilops tauschii genome assembly Aet v5.0 features greater sequence contiguity and improved annotation.

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and an important genetic resource. The reference-quality genome sequence Aet v4.0 for Ae. tauschii acc. AL8/78 was therefore an important milestone for wheat biology and breeding. Further advances in sequencing acc. AL8/78 and release of the Aet v5.0 sequence assembly are reported here. Two new optical maps were constructed and used in the revision of pseudomolecules. Gaps were closed with Pacific Biosciences long-read contigs, decreasing the gap number by 38,899. Transposable elements and protein-coding genes were reannotated. The number of annotated high-confidence genes was reduced from 39,635 in Aet v4.0 to 32,885 in Aet v5.0. A total of 2245 biologically important genes, including those affecting plant phenology, grain quality, and tolerance of abiotic stresses in wheat, was manually annotated and disease-resistance genes were annotated by a dedicated pipeline. Disease-resistance genes encoding nucleotide-binding site domains, receptor-like protein kinases, and receptor-like proteins were preferentially located in distal chromosome regions, whereas those encoding transmembrane coiled-coil proteins were dispersed more evenly along the chromosomes. Discovery, annotation, and expression analyses of microRNA (miRNA) precursors, mature miRNAs, and phasiRNAs are reported, including miRNA target genes. Other small RNAs, such as hc-siRNAs and tRFs, were characterized. These advances enhance the utility of the Ae. tauschii genome sequence for wheat genetics, biotechnology, and breeding.

Developing Chloroplast Genomic Resources from 25 Avena Species for the Characterization of Oat Wild Relative Germplasm.

Chloroplast (cp) genomics will play an important role in the characterization of crop wild relative germplasm conserved in worldwide gene banks, thanks to the advances in genome sequencing. We applied a multiplexed shotgun sequencing procedure to sequence the cp genomes of 25 Avena species with variable ploidy levels. Bioinformatics analysis of the acquired sequences generated 25 de novo genome assemblies ranging from 135,557 to 136,006 bp. The gene annotations revealed 130 genes and their duplications, along with four to six pseudogenes, for each genome. Little differences in genome structure and gene arrangement were observed across the 25 species. Polymorphism analyses identified 1313 polymorphic sites and revealed an average of 277 microsatellites per genome. Greater nucleotide diversity was observed in the short single-copy region. Genome-wide scanning of selection signals suggested that six cp genes were under positive selection on some amino acids. These research outputs allow for a better understanding of oat cp genomes and evolution, and they form an essential set of cp genomic resources for the studies of oat evolutionary biology and for oat wild relative germplasm characterization.

Fine-Tuning of MiR528 Accumulation Modulates Flowering Time in Rice.

In plants, microRNA (miRNA) functions in the post-transcriptional repression of target mRNAs have been well explored. However, the mechanisms regulating the accumulation of miRNAs remain poorly understood. Here, we report that distinct mechanisms regulate accumulation of a monocot-specific miRNA, rice (Oryza sativa) miR528. At the transcriptional level, miR528 accumulated to higher levels in older plants than in young seedlings and exhibited aging-modulated gradual accumulation and diurnal rhythms in leaves; at the post-transcriptional level, aging also modulated miR528 levels by enhancing pri-miR528 alternative splicing. We found that miR528 promotes rice flowering under long-day conditions by targeting RED AND FAR-RED INSENSITIVE 2 (OsRFI2). Moreover, natural variations in the MIR528 promoter region caused differences in miR528 expression among rice varieties, which are correlated with their different binding affinities with the transcription factor OsSPL9 that activates the expression of miR528. Taken together, our findings reveal rice plants have evolved sophisticated modes fine-tuning miR528 levels and provide insight into the mechanisms that regulate MIRNA expression in plants.

Evaluation of Genomic Prediction for Pasmo Resistance in Flax.

Pasmo (Septoria linicola) is a fungal disease causing major losses in seed yield and quality and stem fibre quality in flax. Pasmo resistance (PR) is quantitative and has low heritability. To improve PR breeding efficiency, the accuracy of genomic prediction (GP) was evaluated using a diverse worldwide core collection of 370 accessions. Four marker sets, including three defined by 500, 134 and 67 previously identified quantitative trait loci (QTL) and one of 52,347 PR-correlated genome-wide single nucleotide polymorphisms, were used to build ridge regression best linear unbiased prediction (RR-BLUP) models using pasmo severity (PS) data collected from field experiments performed during five consecutive years. With five-fold random cross-validation, GP accuracy as high as 0.92 was obtained from the models using the 500 QTL when the average PS was used as the training dataset. GP accuracy increased with training population size, reaching values >0.9 with training population size greater than 185. Linear regression of the observed PS with the number of positive-effect QTL in accessions provided an alternative GP approach with an accuracy of 0.86. The results demonstrate the GP models based on marker information from all identified QTL and the 5-year PS average is highly effective for PR prediction.

Genome-Wide Association Study and Selection Signatures Detect Genomic Regions Associated with Seed Yield and Oil Quality in Flax.

A genome-wide association study (GWAS) was performed on a set of 260 lines which belong to three different bi-parental flax mapping populations. These lines were sequenced to an averaged genome coverage of 19× using the Illumina Hi-Seq platform. Phenotypic data for 11 seed yield and oil quality traits were collected in eight year/location environments. A total of 17,288 single nucleotide polymorphisms were identified, which explained more than 80% of the phenotypic variation for days to maturity (DTM), iodine value (IOD), palmitic (PAL), stearic, linoleic (LIO) and linolenic (LIN) acid contents. Twenty-three unique genomic regions associated with 33 quantitative trait loci (QTL) for the studied traits were detected, thereby validating four genomic regions previously identified. The 33 QTL explained 48⁻73% of the phenotypic variation for oil content, IOD, PAL, LIO and LIN but only 8⁻14% for plant height, DTM and seed yield. A genome-wide selective sweep scan for selection signatures detected 114 genomic regions that accounted for 7.82% of the flax pseudomolecule and overlapped with the 11 GWAS-detected genomic regions associated with 18 QTL for 11 traits. The results demonstrate the utility of GWAS combined with selection signatures for dissection of the genetic structure of traits and for pinpointing genomic regions for breeding improvement.

Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax.

Genomes of varying sizes have been sequenced with next-generation sequencing platforms. However, most reference sequences include draft unordered scaffolds containing chimeras caused by mis-scaffolding. A BioNano genome (BNG) optical map was constructed to improve the previously sequenced flax genome (Linum usitatissimum L., 2n = 30, about 373 Mb), which consisted of 3852 scaffolds larger than 1 kb and totalling 300.6 Mb. The high-resolution BNG map of cv. CDC Bethune totalled 317 Mb and consisted of 251 BNG contigs with an N50 of 2.15 Mb. A total of 622 scaffolds (286.6 Mb, 94.9%) aligned to 211 BNG contigs (298.6 Mb, 94.2%). Of those, 99 scaffolds, diagnosed to contain assembly errors, were refined into 225 new scaffolds. Using the newly refined scaffold sequences and the validated bacterial artificial chromosome-based physical map of CDC Bethune, the 211 BNG contigs were scaffolded into 94 super-BNG contigs (N50 of 6.64 Mb) that were further assigned to the 15 flax chromosomes using the genetic map. The pseudomolecules total about 316 Mb, with individual chromosomes of 15.6 to 29.4 Mb, and cover 97% of the annotated genes. Evidence from the chromosome-scale pseudomolecules suggests that flax has undergone palaeopolyploidization and mesopolyploidization events, followed by rearrangements and deletions or fusion of chromosome arms from an ancient progenitor with a haploid chromosome number of eight.

Genome-Wide Association Studies for Pasmo Resistance in Flax (Linum usitatissimum L.).

Pasmo is one of the most widespread diseases threatening flax production. To identify genetic regions associated with pasmo resistance (PR), a genome-wide association study was performed on 370 accessions from the flax core collection. Evaluation of pasmo severity was performed in the field from 2012 to 2016 in Morden, MB, Canada. Genotyping-by-sequencing has identified 258,873 single nucleotide polymorphisms (SNPs) distributed on all 15 flax chromosomes. Marker-trait associations were identified using ten different statistical models. A total of 692 unique quantitative trait nucleotides (QTNs) associated with 500 putative quantitative trait loci (QTL) were detected from six phenotypic PR datasets (five individual years and average across years). Different QTNs were identified with various statistical models and from individual PR datasets, indicative of the complementation between analytical methods and/or genotype × environment interactions of the QTL effects. The single-locus models tended to identify large-effect QTNs while the multi-loci models were able to detect QTNs with smaller effects. Among the putative QTL, 67 had large effects (3-23%), were stable across all datasets and explained 32-64% of the total variation for PR in the various datasets. Forty-five of these QTL spanned 85 resistance gene analogs including a large toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (TNL) type gene cluster on chromosome 8. The number of QTL with positive-effect or favorite alleles (NPQTL) in accessions was significantly correlated with PR (R 2 = 0.55), suggesting that these QTL effects are mainly additive. NPQTL was also significantly associated with morphotype (R 2 = 0.52) and major QTL with positive effect alleles were present in the fiber type accessions. The 67 large effect QTL are suited for marker-assisted selection and the 500 QTL for effective genomic prediction in PR molecular breeding.

Genotyping-by-sequencing data of 272 crested wheatgrass (Agropyron cristatum) genotypes.

Crested wheatgrass [Agropyron cristatum L. (Gaertn.)] is an important cool-season forage grass widely used for early spring grazing. However, the genomic resources for this non-model plant are still lacking. Our goal was to generate the first set of next generation sequencing data using the genotyping-by-sequencing technique. A total of 272 crested wheatgrass plants representing seven breeding lines, five cultivars and five geographically diverse accessions were sequenced with an Illumina MiSeq instrument. These sequence datasets were processed using different bioinformatics tools to generate contigs for diploid and tetraploid plants and SNPs for diploid plants. Together, these genomic resources form a fundamental basis for genomic studies of crested wheatgrass and other wheatgrass species. The raw reads were deposited into Sequence Read Archive (SRA) database under NCBI accession SRP115373 (https://www.ncbi.nlm.nih.gov/sra?term=SRP115373) and the supplementary datasets are accessible in Figshare (10.6084/m9.figshare.5345092).

Genome sequence of the progenitor of the wheat D genome Aegilops tauschii.

Aegilops tauschii is the diploid progenitor of the D genome of hexaploid wheat (Triticum aestivum, genomes AABBDD) and an important genetic resource for wheat. The large size and highly repetitive nature of the Ae. tauschii genome has until now precluded the development of a reference-quality genome sequence. Here we use an array of advanced technologies, including ordered-clone genome sequencing, whole-genome shotgun sequencing, and BioNano optical genome mapping, to generate a reference-quality genome sequence for Ae. tauschii ssp. strangulata accession AL8/78, which is closely related to the wheat D genome. We show that compared to other sequenced plant genomes, including a much larger conifer genome, the Ae. tauschii genome contains unprecedented amounts of very similar repeated sequences. Our genome comparisons reveal that the Ae. tauschii genome has a greater number of dispersed duplicated genes than other sequenced genomes and its chromosomes have been structurally evolving an order of magnitude faster than those of other grass genomes. The decay of colinearity with other grass genomes correlates with recombination rates along chromosomes. We propose that the vast amounts of very similar repeated sequences cause frequent errors in recombination and lead to gene duplications and structural chromosome changes that drive fast genome evolution.

ROS accumulation and antiviral defence control by microRNA528 in rice.

MicroRNAs (miRNAs) are key regulators of plant-pathogen interactions. Modulating miRNA function has emerged as a new strategy to produce virus resistance traits1-5. However, the miRNAs involved in antiviral defence and the underlying mechanisms remain largely elusive. We previously demonstrated that sequestration by Argonaute (AGO) proteins plays an important role in regulating miRNA function in antiviral defence pathways6. Here we reveal that cleavage-defective AGO18 complexes sequester microRNA528 (miR528) upon viral infection. We show that miR528 negatively regulates viral resistance in rice by cleaving L-ascorbate oxidase (AO) messenger RNA, thereby reducing AO-mediated accumulation of reactive oxygen species. Upon viral infection, miR528 becomes preferentially associated with AGO18, leading to elevated AO activity, higher basal reactive oxygen species accumulation and enhanced antiviral defence. Our findings reveal a mechanism in which antiviral defence is boosted through suppression of an miRNA that negatively regulates viral resistance. This mechanism could be manipulated to engineer virus-resistant crop plants.

Pseudogenes and Their Genome-Wide Prediction in Plants.

Pseudogenes are paralogs generated from ancestral functional genes (parents) during genome evolution, which contain critical defects in their sequences, such as lacking a promoter, having a premature stop codon or frameshift mutations. Generally, pseudogenes are functionless, but recent evidence demonstrates that some of them have potential roles in regulation. The majority of pseudogenes are generated from functional progenitor genes either by gene duplication (duplicated pseudogenes) or retro-transposition (processed pseudogenes). Pseudogenes are primarily identified by comparison to their parent genes. Bioinformatics tools for pseudogene prediction have been developed, among which PseudoPipe, PSF and Shiu's pipeline are publicly available. We compared these three tools using the well-annotated Arabidopsis thaliana genome and its known 924 pseudogenes as a test data set. PseudoPipe and Shiu's pipeline identified ~80% of A. thaliana pseudogenes, of which 94% were shared, while PSF failed to generate adequate results. A need for improvement of the bioinformatics tools for pseudogene prediction accuracy in plant genomes was thus identified, with the ultimate goal of improving the quality of genome annotation in plants.

RGAugury: a pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants.

BACKGROUND: Resistance gene analogs (RGAs), such as NBS-encoding proteins, receptor-like protein kinases (RLKs) and receptor-like proteins (RLPs), are potential R-genes that contain specific conserved domains and motifs. Thus, RGAs can be predicted based on their conserved structural features using bioinformatics tools. Computer programs have been developed for the identification of individual domains and motifs from the protein sequences of RGAs but none offer a systematic assessment of the different types of RGAs. A user-friendly and efficient pipeline is needed for large-scale genome-wide RGA predictions of the growing number of sequenced plant genomes. RESULTS: An integrative pipeline, named RGAugury, was developed to automate RGA prediction. The pipeline first identifies RGA-related protein domains and motifs, namely nucleotide binding site (NB-ARC), leucine rich repeat (LRR), transmembrane (TM), serine/threonine and tyrosine kinase (STTK), lysin motif (LysM), coiled-coil (CC) and Toll/Interleukin-1 receptor (TIR). RGA candidates are identified and classified into four major families based on the presence of combinations of these RGA domains and motifs: NBS-encoding, TM-CC, and membrane associated RLP and RLK. All time-consuming analyses of the pipeline are paralleled to improve performance. The pipeline was evaluated using the well-annotated Arabidopsis genome. A total of 98.5, 85.2, and 100 % of the reported NBS-encoding genes, membrane associated RLPs and RLKs were validated, respectively. The pipeline was also successfully applied to predict RGAs for 50 sequenced plant genomes. A user-friendly web interface was implemented to ease command line operations, facilitate visualization and simplify result management for multiple datasets. CONCLUSIONS: RGAugury is an efficiently integrative bioinformatics tool for large scale genome-wide identification of RGAs. It is freely available at Bitbucket: https://bitbucket.org/yaanlpc/rgaugury .

Synteny analysis in Rosids with a walnut physical map reveals slow genome evolution in long-lived woody perennials.

BACKGROUND: Mutations often accompany DNA replication. Since there may be fewer cell cycles per year in the germlines of long-lived than short-lived angiosperms, the genomes of long-lived angiosperms may be diverging more slowly than those of short-lived angiosperms. Here we test this hypothesis. RESULTS: We first constructed a genetic map for walnut, a woody perennial. All linkage groups were short, and recombination rates were greatly reduced in the centromeric regions. We then used the genetic map to construct a walnut bacterial artificial chromosome (BAC) clone-based physical map, which contained 15,203 exonic BAC-end sequences, and quantified with it synteny between the walnut genome and genomes of three long-lived woody perennials, Vitis vinifera, Populus trichocarpa, and Malus domestica, and three short-lived herbs, Cucumis sativus, Medicago truncatula, and Fragaria vesca. Each measure of synteny we used showed that the genomes of woody perennials were less diverged from the walnut genome than those of herbs. We also estimated the nucleotide substitution rate at silent codon positions in the walnut lineage. It was one-fifth and one-sixth of published nucleotide substitution rates in the Medicago and Arabidopsis lineages, respectively. We uncovered a whole-genome duplication in the walnut lineage, dated it to the neighborhood of the Cretaceous-Tertiary boundary, and allocated the 16 walnut chromosomes into eight homoeologous pairs. We pointed out that during polyploidy-dysploidy cycles, the dominant tendency is to reduce the chromosome number. CONCLUSION: Slow rates of nucleotide substitution are accompanied by slow rates of synteny erosion during genome divergence in woody perennials.

Disease Resistance Gene Analogs (RGAs) in Plants.

Plants have developed effective mechanisms to recognize and respond to infections caused by pathogens. Plant resistance gene analogs (RGAs), as resistance (R) gene candidates, have conserved domains and motifs that play specific roles in pathogens' resistance. Well-known RGAs are nucleotide binding site leucine rich repeats, receptor like kinases, and receptor like proteins. Others include pentatricopeptide repeats and apoplastic peroxidases. RGAs can be detected using bioinformatics tools based on their conserved structural features. Thousands of RGAs have been identified from sequenced plant genomes. High-density genome-wide RGA genetic maps are useful for designing diagnostic markers and identifying quantitative trait loci (QTL) or markers associated with plant disease resistance. This review focuses on recent advances in structures and mechanisms of RGAs, and their identification from sequenced genomes using bioinformatics tools. Applications in enhancing fine mapping and cloning of plant disease resistance genes are also discussed.

Secondary siRNAs from Medicago NB-LRRs modulated via miRNA-target interactions and their abundances.

Small RNAs are a class of non-coding RNAs that are of great importance in gene expression regulatory networks. Different families of small RNAs are generated via distinct biogenesis pathways. One such family specific to plants is that of phased, secondary siRNAs (phasiRNAs); these require RDR6, DCL4, and (typically) a microRNA (miRNA) trigger for their biogenesis. Protein-encoding genes are an important source of phasi-RNAs. The model legume Medicago truncatula generates phasiRNAs from many PHAS loci, and we aimed to investigate their biogenesis and mechanism by which miRNAs trigger these molecules. We modulated miRNA abundances in transgenic tissues showing that the abundance of phasiRNAs correlates with the levels of both miRNA triggers and the target, precursor transcripts. We identified sets of phasiRNAs or PHAS loci that predominantly and substantially increase in response to miRNA overexpression. In the process of validating targets from miRNA overexpression tissues, we found that in the miRNA-mRNA target pairing, the 3' terminal nucleotide (the 22nd position), but not the 10th position, is important for phasiRNA production. Mutating the single 3' terminal nucleotide dramatically diminishes phasiRNA production. Ectopic expression of Medicago NB-LRR-targeting miRNAs in Arabidopsis showed that only a few NB-LRRs are capable of phasiRNA production; our data indicate that this might be due to target inaccessibility determined by sequences flanking target sites. Our results suggest that target accessibility is an important component in miRNA-target interactions that could be utilized in target prediction, and the evolution of mRNA sequences flanking miRNA-target sites may be impacted.

An integrated workflow for DNA methylation analysis.

The analysis of cytosine methylation provides a new way to assess and describe epigenetic regulation at a whole-genome level in many eukaryotes. DNA methylation has a demonstrated role in the genome stability and protection, regulation of gene expression and many other aspects of genome function and maintenance. BS-seq is a relatively unbiased method for profiling the DNA methylation, with a resolution capable of measuring methylation at individual cytosines. Here we describe, as an example, a workflow to handle DNA methylation analysis, from BS-seq library preparation to the data visualization. We describe some applications for the analysis and interpretation of these data. Our laboratory provides public access to plant DNA methylation data via visualization tools available at our "Next-Gen Sequence" websites (http://mpss.udel.edu), along with small RNA, RNA-seq and other data types.

The p14ARF alternate reading frame protein enhances DNA binding of topoisomerase I by interacting with the serine 506-phosphorylated core domain.

In addition to its well-characterized function as a tumor suppressor, p14ARF (ARF) is a positive regulator of topoisomerase I (topo I), a central enzyme in DNA metabolism and a target for cancer therapy. We previously showed that topo I hyperphosphorylation, a cancer-associated event mediated by elevated levels of the protein kinase CK2, increases topo I activity and the cellular sensitivity to topo I-targeted drugs. Topo I hyperphosphorylation also increases its interaction with ARF. Because the ARF-topo I interaction could be highly relevant to DNA metabolism and cancer treatment, we identified the regions of topo I involved in ARF binding and characterized the effects of ARF binding on topo I function. Using a series of topo I deletion constructs, we found that ARF interacted with the topo I core domain, which encompasses most of the catalytic and DNA-interacting residues. ARF binding increased the DNA relaxation activity of hyperphosphorylated topo I by enhancing its association with DNA, but did not affect the topo I catalytic rate. In cells, ARF promoted the chromatin association of hyperphosphorylated, but not basal phosphorylated, topo I, and increased topo I-mediated DNA nicking under conditions of oxidative stress. The aberrant nicking was found to correlate with increased formation of DNA double-strand breaks, which are precursors of many genome destabilizing events. The results suggest that the convergent actions of oxidative stress and elevated CK2 and ARF levels, which are common features of cancer cells, lead to a dysregulation of topo I that may contribute both to the cellular response to topo I-targeted drugs and to genome instability.