Towards a public analysis database for LHC new physics searches using MadAnalysis 5

We present the implementation, in the MadAnalysis 5 framework, of several ATLAS and CMS searches for supersymmetry in data recorded during the first run of the LHC. We provide extensive details on the validation of our implementations and propose to create a public analysis database within this framework.Comment: 20 pages, 15 figures, 5 recast codes; version accepted by EPJC (Dec 22, 2014) including a new section with guidelines for the experimental collaborations as well as for potential contributors to the PAD; complementary information can be found at http://madanalysis.irmp.ucl.ac.be/wiki/PhysicsAnalysisDatabas

Statistical and agent-based modelling of the transmissibility of different SARS-CoV-2 variants in England and impact of different interventions

The English SARS-CoV-2 epidemic has been affected by the emergence of new viral variants such as B.1.177, Alpha and Delta, and changing restrictions. We used statistical models and the agent-based model Covasim, in June 2021, to estimate B.1.177 to be 20% more transmissible than the wild type, Alpha to be 50-80% more transmissible than B.1.177 and Delta to be 65-90% more transmissible than Alpha. Using these estimates in Covasim (calibrated 1 September 2020 to 20 June 2021), in June 2021, we found that due to the high transmissibility of Delta, resurgence in infections driven by the Delta variant would not be prevented, but would be strongly reduced by delaying the relaxation of restrictions by one month and with continued vaccination. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'

Natural Islands for a 125 GeV Higgs in the scale-invariant NMSSM

We study whether a 125 GeV standard model-like Higgs boson can be accommodated within the scale-invariant NMSSM in a way that is natural in all respects, i.e., not only is the stop mass and hence its loop contribution to Higgs mass of natural size, but we do not allow significant tuning of NMSSM parameters as well. We pursue as much as possible an analytic approach which gives clear insights on various ways to accommodate such a Higgs mass, while conducting complementary numerical analyses. We consider both scenarios with singlet-like state being heavier and lighter than SM-like Higgs. With A-terms being small, we find for the NMSSM to be perturbative up to GUT scale, it is not possible to get 125 GeV Higgs mass, which is true even if we tune parameters of NMSSM. If we allow some of the couplings to become non-perturbative below the GUT scale, then the non-tuned option implies that the singlet self-coupling, kappa, is larger than the singlet-Higgs coupling, lambda, which itself is order 1. This leads to a Landau pole for these couplings close to the weak scale, in particular below ~10^4 TeV. In both the perturbative and non-perturbative NMSSM, allowing large A_lambda, A_kappa gives "more room" to accommodate a 125 GeV Higgs, but a tuning of these A-terms may be needed. In our analysis we also conduct a careful study of the constraints on the parameter space from requiring global stability of the desired vacuum fitting a 125 GeV Higgs, which is complementary to existing literature. In particular, as the singlet-Higgs coupling lambda increases, vacuum stability becomes more serious of an issue.Comment: 34 pages, 4 figures, references added, minor corrections to text and figures, version to be published in JHE

Les Houches 2011: Physics at TeV Colliders New Physics Working Group Report

We present the activities of the "New Physics" working group for the "Physics at TeV Colliders" workshop (Les Houches, France, 30 May-17 June, 2011). Our report includes new agreements on formats for interfaces between computational tools, new tool developments, important signatures for searches at the LHC, recommendations for presentation of LHC search results, as well as additional phenomenological studies.Comment: 243 pages, report of the Les Houches 2011 New Physics Group; fix three figure

Horizontal DNA transfer mechanisms of bacteria as weapons of intragenomic conflict

Horizontal DNA transfer (HDT) is a pervasive mechanism of diversification in many microbial species, but its primary evolutionary role remains controversial. Much recent research has emphasised the adaptive benefit of acquiring novel DNA, but here we argue instead that intragenomic conflict provides a coherent framework for understanding the evolutionary origins of HDT. To test this hypothesis, we developed a mathematical model of a clonally descended bacterial population undergoing HDT through transmission of mobile genetic elements (MGEs) and genetic transformation. Including the known bias of transformation toward the acquisition of shorter alleles into the model suggested it could be an effective means of counteracting the spread of MGEs. Both constitutive and transient competence for transformation were found to provide an effective defence against parasitic MGEs; transient competence could also be effective at permitting the selective spread of MGEs conferring a benefit on their host bacterium. The coordination of transient competence with cell-cell killing, observed in multiple species, was found to result in synergistic blocking of MGE transmission through releasing genomic DNA for homologous recombination while simultaneously reducing horizontal MGE spread by lowering the local cell density. To evaluate the feasibility of the functions suggested by the modelling analysis, we analysed genomic data from longitudinal sampling of individuals carrying Streptococcus pneumoniae. This revealed the frequent within-host coexistence of clonally descended cells that differed in their MGE infection status, a necessary condition for the proposed mechanism to operate. Additionally, we found multiple examples of MGEs inhibiting transformation through integrative disruption of genes encoding the competence machinery across many species, providing evidence of an ongoing "arms race." Reduced rates of transformation have also been observed in cells infected by MGEs that reduce the concentration of extracellular DNA through secretion of DNases. Simulations predicted that either mechanism of limiting transformation would benefit individual MGEs, but also that this tactic's effectiveness was limited by competition with other MGEs coinfecting the same cell. A further observed behaviour we hypothesised to reduce elimination by transformation was MGE activation when cells become competent. Our model predicted that this response was effective at counteracting transformation independently of competing MGEs. Therefore, this framework is able to explain both common properties of MGEs, and the seemingly paradoxical bacterial behaviours of transformation and cell-cell killing within clonally related populations, as the consequences of intragenomic conflict between self-replicating chromosomes and parasitic MGEs. The antagonistic nature of the different mechanisms of HDT over short timescales means their contribution to bacterial evolution is likely to be substantially greater than previously appreciated

Easy and accurate reconstruction of whole HIV genomes from short-read sequence data

Abstract Next-generation sequencing has yet to be widely adopted for HIV. The difficulty of accurately reconstructing the consensus sequence of a quasispecies from reads (short fragments of DNA) in the presence of rapid between- and within-host evolution may have presented a barrier. In particular, mapping (aligning) reads to a reference sequence leads to biased loss of information; this bias can distort epidemiological and evolutionary conclusions. De novo assembly avoids this bias by effectively aligning the reads to themselves, producing a set of sequences called contigs. However contigs provide only a partial summary of the reads, misassembly may result in their having an incorrect structure, and no information is available at parts of the genome where contigs could not be assembled. To address these problems we developed the tool shiver to preprocess reads for quality and contamination, then map them to a reference tailored to the sample using corrected contigs supplemented with existing reference sequences. Run with two commands per sample, it can easily be used for large heterogeneous data sets. We use shiver to reconstruct the consensus sequence and minority variant information from paired-end short-read data produced with the Illumina platform, for 65 existing publicly available samples and 50 new samples. We show the systematic superiority of mapping to shiver ’s constructed reference over mapping the same reads to the standard reference HXB2: an average of 29 bases per sample are called differently, of which 98.5% are supported by higher coverage. We also provide a practical guide to working with imperfect contigs

HIV-1 full-genome phylogenetics of generalized epidemics in sub-Saharan Africa: impact of missing nucleotide characters in next-generation sequences.

To characterize HIV-1 transmission dynamics in regions where the burden of HIV-1 is greatest, the 'Phylogenetics and Networks for Generalised HIV Epidemics in Africa' consortium (PANGEA-HIV) is sequencing full-genome viral isolates from across sub-Saharan Africa. We report the first 3,985 PANGEA-HIV consensus sequences from four cohort sites (Rakai Community Cohort Study, n=2,833; MRC/UVRI Uganda, n=701; Mochudi Prevention Project, n=359; Africa Health Research Institute Resistance Cohort, n=92). Next-generation sequencing success rates varied: more than 80% of the viral genome from the gag to the nef genes could be determined for all sequences from South Africa, 75% of sequences from Mochudi, 60% of sequences from MRC/UVRI Uganda, and 22% of sequences from Rakai. Partial sequencing failure was primarily associated with low viral load, increased for amplicons closer to the 3' end of the genome, was not associated with subtype diversity except HIV-1 subtype D, and remained significantly associated with sampling location after controlling for other factors. We assessed the impact of the missing data patterns in PANGEA-HIV sequences on phylogeny reconstruction in simulations. We found a threshold in terms of taxon sampling below which the patchy distribution of missing characters in next-generation sequences has an excess negative impact on the accuracy of HIV-1 phylogeny reconstruction, which is attributable to tree reconstruction artifacts that accumulate when branches in viral trees are long. The large number of PANGEA-HIV sequences provides unprecedented opportunities for evaluating HIV-1 transmission dynamics across sub-Saharan Africa and identifying prevention opportunities. Molecular epidemiological analyses of these data must proceed cautiously because sequence sampling remains below the identified threshold and a considerable negative impact of missing characters on phylogeny reconstruction is expected

Quantifying HIV transmission flow between high-prevalence hotspots and surrounding communities: a population-based study in Rakai, Uganda.

BACKGROUND: International and global organisations advocate targeting interventions to areas of high HIV prevalence (ie, hotspots). To better understand the potential benefits of geo-targeted control, we assessed the extent to which HIV hotspots along Lake Victoria sustain transmission in neighbouring populations in south-central Uganda. METHODS: We did a population-based survey in Rakai, Uganda, using data from the Rakai Community Cohort Study. The study surveyed all individuals aged 15-49 years in four high-prevalence Lake Victoria fishing communities and 36 neighbouring inland communities. Viral RNA was deep sequenced from participants infected with HIV who were antiretroviral therapy-naive during the observation period. Phylogenetic analysis was used to infer partial HIV transmission networks, including direction of transmission. Reconstructed networks were interpreted through data for current residence and migration history. HIV transmission flows within and between high-prevalence and low-prevalence areas were quantified adjusting for incomplete sampling of the population. FINDINGS: Between Aug 10, 2011, and Jan 30, 2015, data were collected for the Rakai Community Cohort Study. 25 882 individuals participated, including an estimated 75·7% of the lakeside population and 16·2% of the inland population in the Rakai region of Uganda. 5142 participants were HIV-positive (2703 [13·7%] in inland and 2439 [40·1%] in fishing communities). 3878 (75·4%) people who were HIV-positive did not report antiretroviral therapy use, of whom 2652 (68·4%) had virus deep-sequenced at sufficient quality for phylogenetic analysis. 446 transmission networks were reconstructed, including 293 linked pairs with inferred direction of transmission. Adjusting for incomplete sampling, an estimated 5·7% (95% credibility interval 4·4-7·3) of transmissions occurred within lakeside areas, 89·2% (86·0-91·8) within inland areas, 1·3% (0·6-2·6) from lakeside to inland areas, and 3·7% (2·3-5·8) from inland to lakeside areas. INTERPRETATION: Cross-community HIV transmissions between Lake Victoria hotspots and surrounding inland populations are infrequent and when they occur, virus more commonly flows into rather than out of hotspots. This result suggests that targeted interventions to these hotspots will not alone control the epidemic in inland populations, where most transmissions occur. Thus, geographical targeting of high prevalence areas might not be effective for broader epidemic control depending on underlying epidemic dynamics. FUNDING: The Bill & Melinda Gates Foundation, the National Institute of Allergy and Infectious Diseases, the National Institute of Mental Health, the National Institute of Child Health and Development, the Division of Intramural Research of the National Institute for Allergy and Infectious Diseases, the World Bank, the Doris Duke Charitable Foundation, the Johns Hopkins University Center for AIDS Research, and the President's Emergency Plan for AIDS Relief through the Centers for Disease Control and Prevention

Statistical and agent-based modelling of the transmissibility of different SARS-CoV-2 variants in England and impact of different interventions

The English SARS-CoV-2 epidemic has been affected by the emergence of new viral variants such as B.1.177, Alpha and Delta, and changing restrictions. We used statistical models and the agent-based model Covasim, in June 2021, to estimate B.1.177 to be 20% more transmissible than the wild type, Alpha to be 50-80% more transmissible than B.1.177 and Delta to be 65-90% more transmissible than Alpha. Using these estimates in Covasim (calibrated 1 September 2020 to 20 June 2021), in June 2021, we found that due to the high transmissibility of Delta, resurgence in infections driven by the Delta variant would not be prevented, but would be strongly reduced by delaying the relaxation of restrictions by one month and with continued vaccination. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'