Geographic variation of mutagenic exposures in kidney cancer genomes.

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

Discovery of high-confidence human protein-coding genes and exons by whole-genome PhyloCSF helps elucidate 118 GWAS loci.

The most widely appreciated role of DNA is to encode protein, yet the exact portion of the human genome that is translated remains to be ascertained. We previously developed PhyloCSF, a widely used tool to identify evolutionary signatures of protein-coding regions using multispecies genome alignments. Here, we present the first whole-genome PhyloCSF prediction tracks for human, mouse, chicken, fly, worm, and mosquito. We develop a workflow that uses machine learning to predict novel conserved protein-coding regions and efficiently guide their manual curation. We analyze more than 1000 high-scoring human PhyloCSF regions and confidently add 144 conserved protein-coding genes to the GENCODE gene set, as well as additional coding regions within 236 previously annotated protein-coding genes, and 169 pseudogenes, most of them disabled after primates diverged. The majority of these represent new discoveries, including 70 previously undetected protein-coding genes. The novel coding genes are additionally supported by single-nucleotide variant evidence indicative of continued purifying selection in the human lineage, coding-exon splicing evidence from new GENCODE transcripts using next-generation transcriptomic data sets, and mass spectrometry evidence of translation for several new genes. Our discoveries required simultaneous comparative annotation of other vertebrate genomes, which we show is essential to remove spurious ORFs and to distinguish coding from pseudogene regions. Our new coding regions help elucidate disease-associated regions by revealing that 118 GWAS variants previously thought to be noncoding are in fact protein altering. Altogether, our PhyloCSF data sets and algorithms will help researchers seeking to interpret these genomes, while our new annotations present exciting loci for further experimental characterization.

Clinical features and predictive biomarkers for bladder cancer in patients with type 2 diabetes presenting with haematuria.

AIMS: To identify clinical features and protein biomarkers associated with bladder cancer (BC) in individuals with type 2 diabetes mellitus presenting with haematuria. MATERIALS AND METHODS: Data collected from the Haematuria Biomarker (HaBio) study was used in this analysis. A matched sub-cohort of patients with type 2 diabetes and patients without diabetes was created based on age, sex, and BC diagnosis, using approximately a 1:2 fixed ratio. Randox Biochip Array Technology and ELISA were applied for measurement of 66 candidate serum and urine protein biomarkers. Hazard ratios and 95% confidence intervals were estimated by chi-squared and Wilcoxon rank sum test for clinical features and candidate protein biomarkers. Diagnostic protein biomarker models were identified using Lasso-based binominal regression analysis. RESULTS: There was no difference in BC grade, stage, and severity between individuals with type 2 diabetes and matched controls. Incidence of chronic kidney disease (CKD) was significantly higher in patients with type 2 diabetes (p = 0.008), and CKD was significantly associated with BC in patients with type 2 diabetes (p = 0.032). A biomarker model, incorporating two serum (monocyte chemoattractant protein 1 and vascular endothelial growth factor) and three urine (interleukin 6, cytokeratin 18, and cytokeratin 8) proteins, predicted incidence of BC with an Area Under the Curve (AUC) of 0.84 in individuals with type 2 diabetes. In people without diabetes, the AUC was 0.66. CONCLUSIONS: We demonstrate the potential clinical utility of a biomarker panel, which includes proteins related to BC pathogenesis and type 2 diabetes, for monitoring risk of BC in patients with type 2 diabetes. Earlier urology referral of patients with type 2 diabetes will improve outcomes for these patients. TRIAL REGISTRATION: http://www.isrctn.com/ISRCTN25823942.

A mutation in the putative CRP binding site of the dctA promoter of Salmonella enterica serovar Typhimurium enables growth with low orotate concentrations.

Salmonella enterica and Escherichia coli use the inner membrane transporter DctA to import the pyrimidine biosynthetic pathway intermediate orotate from the environment. To study the regulation of dctA expression, we used an S. enterica serovar Typhimurium pyrimidine auxotroph to select a mutant that could grow in an otherwise nonpermissive culture medium containing glucose and a low concentration of orotate. Whole genome sequencing revealed a point mutation upstream of dctA in the putative cyclic AMP receptor protein (CRP) binding site. The C→T transition converted the least favourable base to the most favourable base for CRP-DNA affinity. A dctA::lux transcriptional fusion confirmed that the mutant dctA promoter gained responsiveness to CRP even in the presence of glucose. Moreover, dctA expression was higher in the mutant than the wild type in the presence of alternative carbon sources that activate CRP.

Redefining the H-NS protein family: a diversity of specialized core and accessory forms exhibit hierarchical transcriptional network integration.

H-NS is a nucleoid structuring protein and global repressor of virulence and horizontally-acquired genes in bacteria. H-NS can interact with itself or with homologous proteins, but protein family diversity and regulatory network overlap remain poorly defined. Here, we present a comprehensive phylogenetic analysis that revealed deep-branching clades, dispelling the presumption that H-NS is the progenitor of varied molecular backups. Each clade is composed exclusively of either chromosome-encoded or plasmid-encoded proteins. On chromosomes, stpA and newly discovered hlpP are core genes in specific genera, whereas hfp and newly discovered hlpC are sporadically distributed. Six clades of H-NS plasmid proteins (Hpp) exhibit ancient and dedicated associations with plasmids, including three clades with fidelity for plasmid incompatibility groups H, F or X. A proliferation of H-NS homologs in Erwiniaceae includes the first observation of potentially co-dependent H-NS forms. Conversely, the observed diversification of oligomerization domains may facilitate stable co-existence of divergent homologs in a genome. Transcriptomic and proteomic analysis in Salmonella revealed regulatory crosstalk and hierarchical control of H-NS homologs. We also discovered that H-NS is both a repressor and activator of Salmonella Pathogenicity Island 1 gene expression, and both regulatory modes are restored by Sfh (HppH) in the absence of H-NS.

Observation of an Intermediate to H2 Binding in a Metal-Organic Framework.

Coordinatively unsaturated metal sites within certain zeolites and metal-organic frameworks can strongly adsorb a wide array of substrates. While many classical examples involve electron-poor metal cations that interact with adsorbates largely through physical interactions, unsaturated electron-rich metal centers housed within porous frameworks can often chemisorb guests amenable to redox activity or covalent bond formation. Despite the promise that materials bearing such sites hold in addressing myriad challenges in gas separations and storage, very few studies have directly interrogated mechanisms of chemisorption at open metal sites within porous frameworks. Here, we show that nondissociative chemisorption of H2 at the trigonal pyramidal Cu+ sites in the metal-organic framework CuI-MFU-4l occurs via the intermediacy of a metastable physisorbed precursor species. In situ powder neutron diffraction experiments enable crystallographic characterization of this intermediate, the first time that this has been accomplished for any material. Evidence for a precursor intermediate is also afforded from temperature-programmed desorption and density functional theory calculations. The activation barrier separating the precursor species from the chemisorbed state is shown to correlate with a change in the Cu+ coordination environment that enhances π-backbonding with H2. Ultimately, these findings demonstrate that adsorption at framework metal sites does not always follow a concerted pathway and underscore the importance of probing kinetics in the design of next-generation adsorbents.

Shiga toxin sub-type 2a increases the efficiency of Escherichia coli O157 transmission between animals and restricts epithelial regeneration in bovine enteroids.

Specific Escherichia coli isolates lysogenised with prophages that express Shiga toxin (Stx) can be a threat to human health, with cattle being an important natural reservoir. In many countries the most severe pathology is associated with enterohaemorrhagic E. coli (EHEC) serogroups that express Stx subtype 2a. In the United Kingdom, phage type (PT) 21/28 O157 strains have emerged as the predominant cause of life-threatening EHEC infections and this phage type commonly encodes both Stx2a and Stx2c toxin types. PT21/28 is also epidemiologically linked to super-shedding (>103 cfu/g of faeces) which is significant for inter-animal transmission and human infection as demonstrated using modelling studies. We demonstrate that Stx2a is the main toxin produced by stx2a+/stx2c+ PT21/28 strains induced with mitomycin C and this is associated with more rapid induction of gene expression from the Stx2a-encoding prophage compared to that from the Stx2c-encoding prophage. Bacterial supernatants containing either Stx2a and/or Stx2c were demonstrated to restrict growth of bovine gastrointestinal organoids with no restriction when toxin production was not induced or prevented by mutation. Isogenic strains that differed in their capacity to produce Stx2a were selected for experimental oral colonisation of calves to assess the significance of Stx2a for both super-shedding and transmission between animals. Restoration of Stx2a expression in a PT21/28 background significantly increased animal-to-animal transmission and the number of sentinel animals that became super-shedders. We propose that while both Stx2a and Stx2c can restrict regeneration of the epithelium, it is the relatively rapid and higher levels of Stx2a induction, compared to Stx2c, that have contributed to the successful emergence of Stx2a+ E. coli isolates in cattle in the last 40 years. We propose a model in which Stx2a enhances E. coli O157 colonisation of in-contact animals by restricting regeneration and turnover of the colonised gastrointestinal epithelium.

Site Specific Relationships between COVID-19 Cases and SARS-CoV-2 Viral Load in Wastewater Treatment Plant Influent.

Wastewater based epidemiology (WBE) has become an important tool during the COVID-19 pandemic, however the relationship between SARS-CoV-2 RNA in wastewater treatment plant influent (WWTP) and cases in the community is not well-defined. We report here the development of a national WBE program across 28 WWTPs serving 50% of the population of Scotland, including large conurbations, as well as low-density rural and remote island communities. For each WWTP catchment area, we quantified spatial and temporal relationships between SARS-CoV-2 RNA in wastewater and COVID-19 cases. Daily WWTP SARS-CoV-2 influent viral RNA load, calculated using daily influent flow rates, had the strongest correlation (ρ > 0.9) with COVID-19 cases within a catchment. As the incidence of COVID-19 cases within a community increased, a linear relationship emerged between cases and influent viral RNA load. There were significant differences between WWTPs in their capacity to predict case numbers based on influent viral RNA load, with the limit of detection ranging from 25 cases for larger plants to a single case in smaller plants. SARS-CoV-2 viral RNA load can be used to predict the number of cases detected in the WWTP catchment area, with a clear statistically significant relationship observed above site-specific case thresholds.

Expert curation of the human and mouse olfactory receptor gene repertoires identifies conserved coding regions split across two exons.

BACKGROUND: Olfactory receptor (OR) genes are the largest multi-gene family in the mammalian genome, with 874 in human and 1483 loci in mouse (including pseudogenes). The expansion of the OR gene repertoire has occurred through numerous duplication events followed by diversification, resulting in a large number of highly similar paralogous genes. These characteristics have made the annotation of the complete OR gene repertoire a complex task. Most OR genes have been predicted in silico and are typically annotated as intronless coding sequences. RESULTS: Here we have developed an expert curation pipeline to analyse and annotate every OR gene in the human and mouse reference genomes. By combining evidence from structural features, evolutionary conservation and experimental data, we have unified the annotation of these gene families, and have systematically determined the protein-coding potential of each locus. We have defined the non-coding regions of many OR genes, enabling us to generate full-length transcript models. We found that 13 human and 41 mouse OR loci have coding sequences that are split across two exons. These split OR genes are conserved across mammals, and are expressed at the same level as protein-coding OR genes with an intronless coding region. Our findings challenge the long-standing and widespread notion that the coding region of a vertebrate OR gene is contained within a single exon. CONCLUSIONS: This work provides the most comprehensive curation effort of the human and mouse OR gene repertoires to date. The complete annotation has been integrated into the GENCODE reference gene set, for immediate availability to the research community.

Mechanisms involved in the adaptation of Escherichia coli O157:H7 to the host intestinal microenvironment.

Host adaptation of pathogens may increase intra- and interspecies transmission. We showed previously that the passage of a clinically isolated enterohemorrhagic Escherichia coli (EHEC) O157 strain (125/99) through the gastrointestinal tract of mice increases its pathogenicity in the same host. In this work, we aimed to elucidate the underlying mechanism(s) involved in the patho-adaptation of the stool-recovered (125RR) strain. We assessed the global transcription profile by microarray and found almost 100 differentially expressed genes in 125RR strain compared with 125/99 strain. We detected an overexpression of Type Three Secretion System (TTSS) proteins at the mRNA and protein levels and demonstrated increased adhesion to epithelial cell lines for the 125RR strain. Additional key attributes of the 125RR strain were: increased motility on semisolid agar, which correlated with an increased fliC mRNA level; reduced Stx2 production at the mRNA and protein levels; increased survival at pH 2.5, as determined by acid resistance assays. We tested whether the overexpression of the LEE-encoded regulator (ler) in trans in the 125/99 strain could recreate the increased pathogenicity observed in the 125RR strain. As anticipated ler overexpression led to increased expression of TTSS proteins and bacterial adhesion to epithelial cells in vitro but also increased mortality and intestinal colonization in vivo. We conclude that this host-adaptation process required changes in several mechanisms that improved EHEC O157 fitness in the new host. The research highlights some of the bacterial mechanisms required for horizontal transmission of these zoonotic pathogens between their animal and human populations.

Design Evaluation of a Novel Multicompartment Unloader Knee Brace.

Knee osteoarthritis (OA) is a significant problem in the aging population, causing pain, impaired mobility, and decreased quality of life. Conservative treatment methods are necessary to reduce rapidly increasing rates of knee joint surgery. Recommended strategies include weight loss and knee bracing to unload knee joint forces. Although weight loss can be beneficial for joint unloading, knee OA patients often find it difficult to lose weight or begin exercise due to knee pain, and not all patients are overweight. Unicompartment offloader knee braces can redistribute joint forces away from one tibiofemoral (TF) compartment; however, <5% of patients have unicompartmental tibiofemoral osteoarthritis (TFOA), while patients with isolated patellofemoral or multicompartmental OA are much more common. By absorbing body weight (BW) and assisting the knee extension moment using a spring-loaded hinge, sufficiently powerful knee-extension-assist (KEA) braces could be useful for unloading the whole knee. This paper (1) describes the design of a spring-loaded tricompartment unloader (TCU) knee brace intended to provide unloading in all three compartments of the knee while weight-bearing, (2) measures and compares the force output of the TCU against the only published and commercially available KEA brace, and (3) calculates the static unloading capacity of each device. The TCU and KEA braces delivered maximum assistive moments equivalent to reducing BW by approximately 45 and 6 lbs, respectively. The paper concludes that sufficiently powerful spring-loaded knee braces show promise in a new class of multicompartment unloader knee orthoses, capable of providing a clinically meaningful unloading effect across all three knee compartments.

Application of a New Multiplexed Array for Rapid, Sensitive, Simultaneous and Quantitative Assessment of Spliced and Unspliced XBP1.

BACKGROUND: IRE1α-mediated unconventional splicing of XBP1 is emerging as a biomarker in several disease states and is indicative of activation of the unfolded protein response sensor IRE1. Splicing of XBP1 mRNA results in the translation of two distinct XBP1 protein isoforms (XBP1s and XBP1u) which, due to post-translational regulation, do not correlate with mRNA levels. As both XBP1 isoforms are implicated in pathogenic or disease progression mechanisms there is a need for a reliable, clinically applicable method to detect them. METHODS: A multiplexed isoform-specific XBP1 array utilising Biochip array technology (BAT™) was assessed for specificity and suitability when using cell protein lysates. The array was applied to RIPA protein lysates from several relevant pre-clinical models with an aim to quantify XBP1 isoforms in comparison with RT-PCR or immunoblot reference methods. RESULTS: A novel reliable, specific and sensitive XBP1 biochip was successfully utilised in pre-clinical research. Application of this biochip to detect XBP1 splicing at the protein level in relevant breast cancer models, under basal conditions as well as pharmacological inhibition and paclitaxel induction, confirmed the findings of previous studies. The biochip was also applied to non-adherent cells and used to quantify changes in the XBP1 isoforms upon activation of the NLRP3 inflammasome. CONCLUSIONS: The XBP1 biochip enables isoform specific quantification of protein level changes upon activation and inhibition of IRE1α RNase activity, using a routine clinical methodology. As such it provides a research tool and potential clinical tool with a quantified, simultaneous, rapid output that is not available from any other published method.

Development of a new biochip array for APOE4 classification from plasma samples using immunoassay-based methods.

BACKGROUND: Apolipoprotein E (APOE) is a key player in lipid transport and metabolism and exists in three common isoforms: APOE2, APOE3 and APOE4. The presence of the E4 allelic variant is recognized as a major genetic risk factor for dementia and other chronic (neuro)degenerative diseases. The availability of a validated assay for rapid and reliable APOE4 classification is therefore advantageous. METHODS: Biochip array technology (BAT) was successfully applied to identify directly the APOE4 status from plasma within 3 h, through simultaneous immunoassay-based detection of both specific APOE4 and total APOE levels. RESULTS: Samples (n=432) were first genotyped by polymerase chain reaction (PCR), and thereafter, using BAT, the corresponding plasma was identified as null, heterozygous or homozygous for the E4 allele by calculating the ratio of APOE4 to total APOE protein. Two centers based in Austria and Ireland correctly classified 170 and 262 samples, respectively, and achieved 100% sensitivity and specificity. CONCLUSIONS: This chemiluminescent biochip-based sandwich immunoarray provides a novel platform to detect rapidly and accurately an individual's APOE4 status directly from plasma. The E4 genotype of individuals has been shown previously to affect presymptomatic risk, prognosis and treatment response for a variety of diseases, including Alzheimer's disease. The biochip's potential for being incorporated in quantitative protein biomarker arrays capable of analyzing disease stages makes it a superior alternative to PCR-based APOE genotyping and may deliver additional protein-specific information on a variety of diseases in the future.

Insights into hominid evolution from the gorilla genome sequence.

Gorillas are humans' closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago. In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

Ensembl 2016.

The Ensembl project (http://www.ensembl.org) is a system for genome annotation, analysis, storage and dissemination designed to facilitate the access of genomic annotation from chordates and key model organisms. It provides access to data from 87 species across our main and early access Pre! websites. This year we introduced three newly annotated species and released numerous updates across our supported species with a concentration on data for the latest genome assemblies of human, mouse, zebrafish and rat. We also provided two data updates for the previous human assembly, GRCh37, through a dedicated website (http://grch37.ensembl.org). Our tools, in particular the VEP, have been improved significantly through integration of additional third party data. REST is now capable of larger-scale analysis and our regulatory data BioMart can deliver faster results. The website is now capable of displaying long-range interactions such as those found in cis-regulated datasets. Finally we have launched a website optimized for mobile devices providing views of genes, variants and phenotypes. Our data is made available without restriction and all code is available from our GitHub organization site (http://github.com/Ensembl) under an Apache 2.0 license.

Alignathon: a competitive assessment of whole-genome alignment methods.

Multiple sequence alignments (MSAs) are a prerequisite for a wide variety of evolutionary analyses. Published assessments and benchmark data sets for protein and, to a lesser extent, global nucleotide MSAs are available, but less effort has been made to establish benchmarks in the more general problem of whole-genome alignment (WGA). Using the same model as the successful Assemblathon competitions, we organized a competitive evaluation in which teams submitted their alignments and then assessments were performed collectively after all the submissions were received. Three data sets were used: Two were simulated and based on primate and mammalian phylogenies, and one was comprised of 20 real fly genomes. In total, 35 submissions were assessed, submitted by 10 teams using 12 different alignment pipelines. We found agreement between independent simulation-based and statistical assessments, indicating that there are substantial accuracy differences between contemporary alignment tools. We saw considerable differences in the alignment quality of differently annotated regions and found that few tools aligned the duplications analyzed. We found that many tools worked well at shorter evolutionary distances, but fewer performed competitively at longer distances. We provide all data sets, submissions, and assessment programs for further study and provide, as a resource for future benchmarking, a convenient repository of code and data for reproducing the simulation assessments.

Hydrogen Uptake on Coordinatively Unsaturated Metal Sites in VSB-5: Strong Binding Affinity Leading to High-Temperature D2/H2 Selectivity.

We examine the adsorption of hydrogen and deuterium into the nanoporous nickel phosphate, VSB-5. On the basis of gas sorption analysis, VSB-5 exhibits one of the highest measured H2 heats of adsorption (HOA) for hydrogen (16 kJ/mol) yet reported. This high HOA is consistent with an unusually large red shift in the Q(1) and Q(0) hydrogen vibrational modes as measured with in situ infrared spectroscopy. The HOA for D2 is measured to be 2 kJ/mol higher than that for H2. "Ideal adsorbed solution theory" analysis of H2 and D2 isotherms provides selectivities above 4 for deuterium at 140 K, suggesting that VSB-5 is a promising adsorbent for pressure-swing adsorption-type separations of hydrogen isotopes.

Explaining the symmetry breaking observed in the endofullerenes H2@C60, HF@C60, and H2O@C60.

Symmetry breaking has been recently observed in the endofullerenes M@C60 (M = H2, HF, H2O), manifesting in the splittings of the three-fold degenerate ground states of the endohedral ortho-H2, ortho-H2O and the j = 1 level of HF. The nature of the interaction causing the symmetry breaking is established in this study. A fragment of the solid C60 is considered, comprised of the central C60 molecule surrounded by twelve nearest-neighbor (NN) C60 molecules. The fullerenes have either P (major) or H (minor) orientational orderings, and are assumed to be rigid with Ih symmetry. Only the central C60 is occupied by the guest molecule M, while the NN fullerenes are all empty. The key proposition of the study is that the electrostatic interactions between the charge densities on the NN C60 molecules and that on M inside the central C60 give rise to the symmetry breaking responsible for the measured level splittings. Using this model, the M@C60 level splittings of interest are calculated variationally and using perturbation theory, for both the P and H orientations. Those obtained for the dominant P orientation are in excellent agreement with the experimental results, with respect to the splitting magnitudes and patterns, for all three M@C60 systems considered, pointing strongly to the quadrupolar M-NN interactions as the main cause of the symmetry breaking. The level splittings calculated for the H orientation are about 30 times smaller than the ones in the P orientation.

A high-resolution map of human evolutionary constraint using 29 mammals.

The comparison of related genomes has emerged as a powerful lens for genome interpretation. Here we report the sequencing and comparative analysis of 29 eutherian genomes. We confirm that at least 5.5% of the human genome has undergone purifying selection, and locate constrained elements covering ∼4.2% of the genome. We use evolutionary signatures and comparisons with experimental data sets to suggest candidate functions for ∼60% of constrained bases. These elements reveal a small number of new coding exons, candidate stop codon readthrough events and over 10,000 regions of overlapping synonymous constraint within protein-coding exons. We find 220 candidate RNA structural families, and nearly a million elements overlapping potential promoter, enhancer and insulator regions. We report specific amino acid residues that have undergone positive selection, 280,000 non-coding elements exapted from mobile elements and more than 1,000 primate- and human-accelerated elements. Overlap with disease-associated variants indicates that our findings will be relevant for studies of human biology, health and disease.