Digital measurement of SARS-CoV-2 transmission risk from 7 million contacts.

How likely is it to become infected by SARS-CoV-2 after being exposed? Almost everyone wondered about this question during the COVID-19 pandemic. Contact-tracing apps1,2 recorded measurements of proximity3 and duration between nearby smartphones. Contacts-individuals exposed to confirmed cases-were notified according to public health policies such as the 2 m, 15 min guideline4,5, despite limited evidence supporting this threshold. Here we analysed 7 million contacts notified by the National Health Service COVID-19 app6,7 in England and Wales to infer how app measurements translated to actual transmissions. Empirical metrics and statistical modelling showed a strong relation between app-computed risk scores and actual transmission probability. Longer exposures at greater distances had risk similar to that of shorter exposures at closer distances. The probability of transmission confirmed by a reported positive test increased initially linearly with duration of exposure (1.1% per hour) and continued increasing over several days. Whereas most exposures were short (median 0.7 h, interquartile range 0.4-1.6), transmissions typically resulted from exposures lasting between 1 h and several days (median 6 h, interquartile range 1.4-28). Households accounted for about 6% of contacts but 40% of transmissions. With sufficient preparation, privacy-preserving yet precise analyses of risk that would inform public health measures, based on digital contact tracing, could be performed within weeks of the emergence of a new pathogen.

Integrated Analysis of Patient Networks and Plasmid Genomes to Investigate a Regional, Multispecies Outbreak of Carbapenemase-Producing Enterobacterales Carrying Both blaIMP and mcr-9 Genes.

BACKGROUND: Carbapenemase-producing Enterobacterales (CPE) are challenging in healthcare, with resistance to multiple classes of antibiotics. This study describes the emergence of imipenemase (IMP)-encoding CPE among diverse Enterobacterales species between 2016 and 2019 across a London regional network. METHODS: We performed a network analysis of patient pathways, using electronic health records, to identify contacts between IMP-encoding CPE-positive patients. Genomes of IMP-encoding CPE isolates were overlaid with patient contacts to imply potential transmission events. RESULTS: Genomic analysis of 84 Enterobacterales isolates revealed diverse species (predominantly Klebsiella spp, Enterobacter spp, and Escherichia coli); 86% (72 of 84) harbored an IncHI2 plasmid carrying blaIMP and colistin resistance gene mcr-9 (68 of 72). Phylogenetic analysis of IncHI2 plasmids identified 3 lineages showing significant association with patient contacts and movements between 4 hospital sites and across medical specialties, which was missed in initial investigations. CONCLUSIONS: Combined, our patient network and plasmid analyses demonstrate an interspecies, plasmid-mediated outbreak of blaIMPCPE, which remained unidentified during standard investigations. With DNA sequencing and multimodal data incorporation, the outbreak investigation approach proposed here provides a framework for real-time identification of key factors causing pathogen spread. Plasmid-level outbreak analysis reveals that resistance spread may be wider than suspected, allowing more interventions to stop transmission within hospital networks.SummaryThis was an investigation, using integrated pathway networks and genomics methods, of the emergence of imipenemase-encoding carbapenemase-producing Enterobacterales among diverse Enterobacterales species between 2016 and 2019 in patients across a London regional hospital network, which was missed on routine investigations.

Bayesian Inference of Clonal Expansions in a Dated Phylogeny.

Microbial population genetics models often assume that all lineages are constrained by the same population size dynamics over time. However, many neutral and selective events can invalidate this assumption and can contribute to the clonal expansion of a specific lineage relative to the rest of the population. Such differential phylodynamic properties between lineages result in asymmetries and imbalances in phylogenetic trees that are sometimes described informally but which are difficult to analyze formally. To this end, we developed a model of how clonal expansions occur and affect the branching patterns of a phylogeny. We show how the parameters of this model can be inferred from a given dated phylogeny using Bayesian statistics, which allows us to assess the probability that one or more clonal expansion events occurred. For each putative clonal expansion event, we estimate its date of emergence and subsequent phylodynamic trajectory, including its long-term evolutionary potential which is important to determine how much effort should be placed on specific control measures. We demonstrate the applicability of our methodology on simulated and real data sets. Inference under our clonal expansion model can reveal important features in the evolution and epidemiology of infectious disease pathogens. [Clonal expansion; genomic epidemiology; microbial population genomics; phylodynamics.].

The bounded coalescent model: Conditioning a genealogy on a minimum root date.

The coalescent model represents how individuals sampled from a population may have originated from a last common ancestor. The bounded coalescent model is obtained by conditioning the coalescent model such that the last common ancestor must have existed after a certain date. This conditioned model arises in a variety of applications, such as speciation, horizontal gene transfer or transmission analysis, and yet the bounded coalescent model has not been previously analysed in detail. Here we describe a new algorithm to simulate from this model directly, without resorting to rejection sampling. We show that this direct simulation algorithm is more computationally efficient than the rejection sampling approach. We also show how to calculate the probability of the last common ancestor occurring after a given date, which is required to compute the probability density of realisations under the bounded coalescent model. Our results are applicable in both the isochronous (when all samples have the same date) and heterochronous (where samples can have different dates) settings. We explore the effect of setting a bound on the date of the last common ancestor, and show that it affects a number of properties of the resulting phylogenies. All our methods are implemented in a new R package called BoundedCoalescent which is freely available online.

PRDM9 Diversity at Fine Geographical Scale Reveals Contrasting Evolutionary Patterns and Functional Constraints in Natural Populations of House Mice.

One of the major challenges in evolutionary biology is the identification of the genetic basis of postzygotic reproductive isolation. Given its pivotal role in this process, here we explore the drivers that may account for the evolutionary dynamics of the PRDM9 gene between continental and island systems of chromosomal variation in house mice. Using a data set of nearly 400 wild-caught mice of Robertsonian systems, we identify the extent of PRDM9 diversity in natural house mouse populations, determine the phylogeography of PRDM9 at a local and global scale based on a new measure of pairwise genetic divergence, and analyze selective constraints. We find 57 newly described PRDM9 variants, this diversity being especially high on Madeira Island, a result that is contrary to the expectations of reduced variation for island populations. Our analysis suggest that the PRDM9 allelic variability observed in Madeira mice might be influenced by the presence of distinct chromosomal fusions resulting from a complex pattern of introgression or multiple colonization events onto the island. Importantly, we detect a significant reduction in the proportion of PRDM9 heterozygotes in Robertsonian mice, which showed a high degree of similarity in the amino acids responsible for protein-DNA binding. Our results suggest that despite the rapid evolution of PRDM9 and the variability detected in natural populations, functional constraints could facilitate the accumulation of allelic combinations that maintain recombination hotspot symmetry. We anticipate that our study will provide the basis for examining the role of different PRDM9 genetic backgrounds in reproductive isolation in natural populations.

Re-emergence of methicillin susceptibility in a resistant lineage of Staphylococcus aureus.

Objectives: MRSA is a leading cause of hospital-associated infection. Acquired resistance is encoded by the mecA gene or its homologue mecC , but little is known about the evolutionary dynamics involved in gain and loss of resistance. The objective of this study was to obtain an expanded understanding of Staphylococcus aureus methicillin resistance microevolution in vivo , by focusing on a single lineage. Methods: We compared the whole-genome sequences of 231 isolates from a single epidemic lineage [clonal complex 30 (CC30) and spa -type t018] of S. aureus that caused an epidemic in the UK. Results: We show that resistance to methicillin in this single lineage was gained on at least two separate occasions, one of which led to a clonal expansion around 1995 presumably caused by a selective advantage. Resistance was, however, subsequently lost in vivo by nine strains isolated between 2008 and 2012. We describe the genetic mechanisms involved in this loss of resistance and the imperfect relationship between genotypic and phenotypic resistance. Conclusions: The recent re-emergence of methicillin susceptibility in this epidemic lineage suggests a significant fitness cost of resistance and reduced selective advantage following the introduction in the mid-2000s of MRSA hospital control measures throughout the UK.

Overcoming challenges in the economic evaluation of interventions to optimise antibiotic use.

Bacteria are becoming increasingly resistant to antibiotics, reducing our ability to treat infections and threatening to undermine modern health care. Optimising antibiotic use is a key element in tackling the problem. Traditional economic evaluation methods do not capture many of the benefits from improved antibiotic use and the potential impact on resistance. Not capturing these benefits is a major obstacle to optimising antibiotic use, as it fails to incentivise the development and use of interventions to optimise the use of antibiotics and preserve their effectiveness (stewardship interventions). Estimates of the benefits of improving antibiotic use involve considerable uncertainty as they depend on the evolution of resistance and associated health outcomes and costs. Here we discuss how economic evaluation methods might be adapted, in the face of such uncertainties. We propose a threshold-based approach that estimates the minimum resistance-related costs that would need to be averted by an intervention to make it cost-effective. If it is probable that without the intervention costs will exceed the threshold then the intervention should be deemed cost-effective.

Drivers of epidemic dynamics in real time from daily digital COVID-19 measurements.

Understanding the drivers of respiratory pathogen spread is challenging, particularly in a timely manner during an ongoing epidemic. In this work, we present insights that we obtained using daily data from the National Health Service COVID-19 app for England and Wales and that we shared with health authorities in almost real time. Our indicator of the reproduction number R(t) was available days earlier than other estimates, with an innovative capability to decompose R(t) into contact rates and probabilities of infection. When Omicron arrived, the main epidemic driver switched from contacts to transmissibility. We separated contacts and transmissions by day of exposure and setting and found pronounced variability over days of the week and during Christmas holidays and events. For example, during the Euro football tournament in 2021, days with England matches showed sharp spikes in exposures and transmissibility. Digital contact-tracing technologies can help control epidemics not only by directly preventing transmissions but also by enabling rapid analysis at scale and with unprecedented resolution.

Natural selection drives the accumulation of amino acid tandem repeats in human proteins.

Amino acid tandem repeats are found in a large number of eukaryotic proteins. They are often encoded by trinucleotide repeats and exhibit high intra- and interspecies size variability due to the high mutation rate associated with replication slippage. The extent to which natural selection is important in shaping amino acid repeat evolution is a matter of debate. On one hand, their high frequency may simply reflect their high probability of expansion by slippage, and they could essentially evolve in a neutral manner. On the other hand, there is experimental evidence that changes in repeat size can influence protein-protein interactions, transcriptional activity, or protein subcellular localization, indicating that repeats could be functionally relevant and thus shaped by selection. To gauge the relative contribution of neutral and selective forces in amino acid repeat evolution, we have performed a comparative analysis of amino acid repeat conservation in a large set of orthologous proteins from 12 vertebrate species. As a neutral model of repeat evolution we have used sequences with the same DNA triplet composition as the coding sequences--and thus expected to be subject to the same mutational forces--but located in syntenic noncoding genomic regions. The results strongly indicate that selection has played a more important role than previously suspected in amino acid tandem repeat evolution, by increasing the repeat retention rate and by modulating repeat size. The data obtained in this study have allowed us to identify a set of 92 repeats that are postulated to play important functional roles due to their strong selective signature, including five cases with direct experimental evidence.

Origin, Genetic Variation and Molecular Epidemiology of SARS-CoV-2 Strains Circulating in Sardinia (Italy) during the First and Second COVID-19 Epidemic Waves.

Understanding how geography and human mobility shape the patterns and spread of infectious diseases such as COVID-19 is key to control future epidemics. An interesting example is provided by the second wave of the COVID-19 epidemic in Europe, which was facilitated by the intense movement of tourists around the Mediterranean coast in summer 2020. The Italian island of Sardinia is a major tourist destination and is widely believed to be the origin of the second Italian wave. In this study, we characterize the genetic variation among SARS-CoV-2 strains circulating in northern Sardinia during the first and second Italian waves using both Illumina and Oxford Nanopore Technologies Next Generation Sequencing methods. Most viruses were placed into a single clade, implying that despite substantial virus inflow, most outbreaks did not spread widely. The second epidemic wave on the island was actually driven by local transmission of a single B.1.177 subclade. Phylogeographic analyses further suggest that those viral strains circulating on the island were not a relevant source for the second epidemic wave in Italy. This result, however, does not rule out the possibility of intense mixing and transmission of the virus among tourists as a major contributor to the second Italian wave.

Genome-wide analysis of histidine repeats reveals their role in the localization of human proteins to the nuclear speckles compartment.

Single amino acid repeats are prevalent in eukaryote organisms, although the role of many such sequences is still poorly understood. We have performed a comprehensive analysis of the proteins containing homopolymeric histidine tracts in the human genome and identified 86 human proteins that contain stretches of five or more histidines. Most of them are endowed with DNA- and RNA-related functions, and, in addition, there is an overrepresentation of proteins expressed in the brain and/or nervous system development. An analysis of their subcellular localization shows that 15 of the 22 nuclear proteins identified accumulate in the nuclear subcompartment known as nuclear speckles. This localization is lost when the histidine repeat is deleted, and significantly, closely related paralogous proteins without histidine repeats also fail to localize to nuclear speckles. Hence, the histidine tract appears to be directly involved in targeting proteins to this compartment. The removal of DNA-binding domains or treatment with RNA polymerase II inhibitors induces the re-localization of several polyhistidine-containing proteins from the nucleoplasm to nuclear speckles. These findings highlight the dynamic relationship between sites of transcription and nuclear speckles. Therefore, we define the histidine repeats as a novel targeting signal for nuclear speckles, and we suggest that these repeats are a way of generating evolutionary diversification in gene duplicates. These data contribute to our better understanding of the physiological role of single amino acid repeats in proteins.

Hospital outbreak of carbapenem-resistant Enterobacterales associated with a blaOXA-48 plasmid carried mostly by Escherichia coli ST399.

A hospital outbreak of carbapenem-resistant Enterobacterales was detected by routine surveillance. Whole genome sequencing and subsequent analysis revealed a conserved promiscuous blaOXA-48 carrying plasmid as the defining factor within this outbreak. Four different species of Enterobacterales were involved in the outbreak. Escherichia coli ST399 accounted for 35 of all the 55 isolates. Comparative genomics analysis using publicly available E. coli ST399 genomes showed that the outbreak E. coli ST399 isolates formed a unique clade. We developed a mathematical model of pOXA-48-like plasmid transmission between host lineages and used it to estimate its conjugation rate, giving a lower bound of 0.23 conjugation events per lineage per year. Our analysis suggests that co-evolution between the pOXA-48-like plasmid and E. coli ST399 could have played a role in the outbreak. This is the first study to report carbapenem-resistant E. coli ST399 carrying blaOXA-48 as the main cause of a plasmid-borne outbreak within a hospital setting. Our findings suggest complementary roles for both plasmid conjugation and clonal expansion in the emergence of this outbreak.

Lessons for preparedness and reasons for concern from the early COVID-19 epidemic in Iran.

INTRODUCTION: Many countries with an early outbreak of SARS-CoV-2 struggled to gauge the size and start date of the epidemic mainly due to limited testing capacities and a large proportion of undetected asymptomatic and mild infections. Iran was among the first countries with a major outbreak outside China. METHODS: We constructed a globally representative sample of 802 genomes, including 46 samples from patients inside or with a travel history to Iran. We then performed a phylogenetic analysis to identify clades related to samples from Iran and estimated the start of the epidemic and early doubling times in cases. We leveraged air travel data from 36 exported cases of COVID-19 to estimate the point-prevalence and the basic reproductive number across the country. We also analysed the province-level all-cause mortality data during winter and spring 2020 to estimate under-reporting of COVID-19-related deaths. Finally, we use this information in an SEIR model to reconstruct the early outbreak dynamics and assess the effectiveness of intervention measures in Iran. RESULTS: By identifying the most basal clade that contained genomes from Iran, our phylogenetic analysis showed that the age of the root is placed on 2019-12-21 (95 % HPD: 2019-09-07 - 2020-02-14). This date coincides with our estimated epidemic start date on 2019-12-25 (95 %CI: 2019-12-11 - 2020-02-24) based air travel data from exported cases with an early doubling time of 4.0 (95 %CI: 1.4-6.7) days in cases. Our analysis of all-cause mortality showed 21.9 (95 % CI: 16.7-27.2) thousand excess deaths by the end of summer. Our model forecasted the second epidemic peak and suggested that by 2020-08-31 a total of 15.0 (95 %CI: 4.9-25.0) million individuals recovered from the disease across the country. CONCLUSION: These findings have profound implications for assessing the stage of the epidemic in Iran despite significant levels of under-reporting. Moreover, the results shed light on the dynamics of SARS-CoV-2 transmissions in Iran and central Asia. They also suggest that in the absence of border screening, there is a high risk of introduction from travellers from areas with active outbreaks. Finally, they show both that well-informed epidemic models are able to forecast episodes of resurgence following a relaxation of interventions, and that NPIs are key to controlling ongoing epidemics.

Database of epidemic trends and control measures during the first wave of COVID-19 in mainland China.

OBJECTIVES: In this data collation study, we aimed to provide a comprehensive database describing the epidemic trends and responses during the first wave of coronavirus disease 2019 (COVID-19) throughout the main provinces in China. METHODS: From mid-January to March 2020, we extracted publicly available data regarding the spread and control of COVID-19 from 31 provincial health authorities and major media outlets in mainland China. Based on these data, we conducted descriptive analyses of the epidemic in the six most-affected provinces. RESULTS: School closures, travel restrictions, community-level lockdown, and contact tracing were introduced concurrently around late January but subsequent epidemic trends differed among provinces. Compared with Hubei, the other five most-affected provinces reported a lower crude case fatality ratio and proportion of critical and severe hospitalised cases. From March 2020, as the local transmission of COVID-19 declined, switching the focus of measures to the testing and quarantine of inbound travellers may have helped to sustain the control of the epidemic. CONCLUSIONS: Aggregated indicators of case notifications and severity distributions are essential for monitoring an epidemic. A publicly available database containing these indicators and information regarding control measures is a useful resource for further research and policy planning in response to the COVID-19 epidemic.

Patterns of RNA Editing in Newcastle Disease Virus Infections.

The expression of accessory non-structural proteins V and W in Newcastle disease virus (NDV) infections depends on RNA editing. These proteins are derived from frameshifts of the sequence coding for the P protein via co-transcriptional insertion of one or two guanines in the mRNA. However, a larger number of guanines can be inserted with lower frequencies. We analysed data from deep RNA sequencing of samples from in vitro and in vivo NDV infections to uncover the patterns of mRNA editing in NDV. The distribution of insertions is well described by a simple Markov model of polymerase stuttering, providing strong quantitative confirmation of the molecular process hypothesised by Kolakofsky and collaborators three decades ago. Our results suggest that the probability that the NDV polymerase would stutter is about 0.45 initially, and 0.3 for further subsequent insertions. The latter probability is approximately independent of the number of previous insertions, the host cell, and viral strain. However, in LaSota infections, we also observe deviations from the predicted V/W ratio of about 3:1 according to this model, which could be attributed to deviations from this stuttering model or to further mechanisms downregulating the abundance of W protein.

Specialization of small non-conjugative plasmids in Escherichia coli according to their family types.

We undertook a comprehensive comparative analysis of a collection of 30 small (<25 kb) non-conjugative Escherichia coli plasmids previously classified by the gene sharing approach into 10 families, as well as plasmids found in the National Center for Biotechnology Information (NCBI) nucleotide database sharing similar genomic sequences. In total, 302 mobilizable (belonging to 2 MOBrep and 5 MOBRNA families) and 106 non-transferable/relaxase-negative (belonging to three ReLRNA families) plasmids were explored. The most striking feature was the specialization of the plasmid family types that was not related to their transmission mode and replication system. We observed a range of host strain specificity, from narrow E. coli host specificity to broad host range specificity, including a wide spectrum of Enterobacteriaceae. We found a wide variety of toxin/antitoxin systems and colicin operons in the plasmids, whose numbers and types varied according to the plasmid family type. The plasmids carried genes conferring resistance spanning almost all of the antibiotic classes, from those to which resistance developed early, such as sulphonamides, to those for which resistance has only developed recently, such as colistin. However, the prevalence of the resistance genes varied greatly according to the family type, ranging from 0 to 100 %. The evolutionary history of the plasmids based on the family type core genes showed variability within family nucleotide divergences in the range of E. coli chromosomal housekeeping genes, indicating long-term co-evolution between plasmids and host strains. In rare cases, a low evolutionary divergence suggested the massive spread of an epidemic plasmid. Overall, the importance of these small non-conjugative plasmids in bacterial adaptation varied greatly according to the type of family they belonged to, with each plasmid family having specific hosts and genetic traits.

Extended-spectrum ß-lactamase-encoding genes are spreading on a wide range of Escherichia coli plasmids existing prior to the use of third-generation cephalosporins.

To understand the evolutionary dynamics of extended-spectrum ß-lactamase (ESBL)-encoding genes in Escherichia coli, we undertook a comparative genomic analysis of 116 whole plasmid sequences of human or animal origin isolated over a period spanning before and after the use of third-generation cephalosporins (3GCs) using a gene-sharing network approach. The plasmids included 82 conjugative, 22 mobilizable and 9 non-transferable plasmids and 3 P-like bacteriophages. ESBL-encoding genes were found on 64 conjugative, 6 mobilizable, 2 non-transferable plasmids and 2 P1-like bacteriophages, indicating that these last three types of mobile elements also play a role, albeit modest, in the diffusion of the ESBLs. The network analysis showed that the plasmids clustered according to their genome backbone type, but not by origin or period of isolation or by antibiotic-resistance type, including type of ESBL-encoding gene. There was no association between the type of plasmid and the phylogenetic history of the parental strains. Finer scale analysis of the more abundant clusters IncF and IncI1 showed that ESBL-encoding plasmids and plasmids isolated before the use of 3GCs had the same diversity and phylogenetic history, and that acquisition of ESBL-encoding genes had occurred during multiple independent events. Moreover, the blaCTX-M-15 gene, unlike other CTX-M genes, was inserted at a hot spot in a blaTEM-1-Tn2 transposon. These findings showed that ESBL-encoding genes have arrived on wide range of pre-existing plasmids and that the successful spread of blaCTX-M-15 seems to be favoured by the presence of well-adapted IncF plasmids that carry a Tn2-blaTEM-1 transposon.

Decomposing the Site Frequency Spectrum: The Impact of Tree Topology on Neutrality Tests.

We investigate the dependence of the site frequency spectrum on the topological structure of genealogical trees. We show that basic population genetic statistics, for instance, estimators of θ or neutrality tests such as Tajima's D, can be decomposed into components of waiting times between coalescent events and of tree topology. Our results clarify the relative impact of the two components on these statistics. We provide a rigorous interpretation of positive or negative values of an important class of neutrality tests in terms of the underlying tree shape. In particular, we show that values of Tajima's D and Fay and Wu's H depend in a direct way on a peculiar measure of tree balance, which is mostly determined by the root balance of the tree. We present a new test for selection in the same class as Fay and Wu's H and discuss its interpretation and power. Finally, we determine the trees corresponding to extreme expected values of these neutrality tests and present formulas for these extreme values as a function of sample size and number of segregating sites.

The conserved nhaAR operon is drastically divergent between B2 and non-B2 Escherichia coli and is involved in extra-intestinal virulence.

The Escherichia coli species is divided in phylogenetic groups that differ in their virulence and commensal distribution. Strains belonging to the B2 group are involved in extra-intestinal pathologies but also appear to be more prevalent as commensals among human occidental populations. To investigate the genetic specificities of B2 sub-group, we used 128 sequenced genomes and identified genes of the core genome that showed marked difference between B2 and non-B2 genomes. We focused on the gene and its surrounding region with the strongest divergence between B2 and non-B2, the antiporter gene nhaA. This gene is part of the nhaAR operon, which is in the core genome but flanked by mobile regions, and is involved in growth at high pH and high sodium concentrations. Consistently, we found that a panel of non-B2 strains grew faster than B2 at high pH and high sodium concentrations. However, we could not identify differences in expression of the nhaAR operon using fluorescence reporter plasmids. Furthermore, the operon deletion had no differential impact between B2 and non-B2 strains, and did not result in a fitness modification in a murine model of gut colonization. Nevertheless, sequence analysis and experiments in a murine model of septicemia revealed that recombination in nhaA among B2 strains was observed in strains with low virulence. Finally, nhaA and nhaAR operon deletions drastically decreased virulence in one B2 strain. This effect of nhaAR deletion appeared to be stronger than deletion of all pathogenicity islands. Thus, a population genetic approach allowed us to identify an operon in the core genome without strong effect in commensalism but with an important role in extra-intestinal virulence, a landmark of the B2 strains.

Detecting correlations among functional-sequence motifs.

Sequence motifs are words of nucleotides in DNA with biological functions, e.g., gene regulation. Identification of such words proceeds through rejection of Markov models on the expected motif frequency along the genome. Additional biological information can be extracted from the correlation structure among patterns of motif occurrences. In this paper a log-linear multivariate intensity Poisson model is estimated via expectation maximization on a set of motifs along the genome of E. coli K12. The proposed approach allows for excitatory as well as inhibitory interactions among motifs and between motifs and other genomic features like gene occurrences. Our findings confirm previous stylized facts about such types of interactions and shed new light on genome-maintenance functions of some particular motifs. We expect these methods to be applicable to a wider set of genomic features.