MapMi: automated mapping of microRNA loci.

BACKGROUND: A large effort to discover microRNAs (miRNAs) has been under way. Currently miRBase is their primary repository, providing annotations of primary sequences, precursors and probable genomic loci. In many cases miRNAs are identical or very similar between related (or in some cases more distant) species. However, miRBase focuses on those species for which miRNAs have been directly confirmed. Secondly, specific miRNAs or their loci are sometimes not annotated even in well-covered species. We sought to address this problem by developing a computational system for automated mapping of miRNAs within and across species. Given the sequence of a known miRNA in one species it is relatively straightforward to determine likely loci of that miRNA in other species. Our primary goal is not the discovery of novel miRNAs but the mapping of validated miRNAs in one species to their most likely orthologues in other species. RESULTS: We present MapMi, a computational system for automated miRNA mapping across and within species. This method has a sensitivity of 92.20% and a specificity of 97.73%. Using the latest release (v14) of miRBase, we obtained 10,944 unannotated potential miRNAs when MapMi was applied to all 21 species in Ensembl Metazoa release 2 and 46 species from Ensembl release 55. CONCLUSIONS: The pipeline and an associated web-server for mapping miRNAs are freely available on http://www.ebi.ac.uk/enright-srv/MapMi/. In addition precomputed miRNA mappings of miRBase miRNAs across a large number of species are provided

Large-scale analysis of microRNA evolution.

BACKGROUND: In animals, microRNAs (miRNA) are important genetic regulators. Animal miRNAs appear to have expanded in conjunction with an escalation in complexity during early bilaterian evolution. Their small size and high-degree of similarity makes them challenging for phylogenetic approaches. Furthermore, genomic locations encoding miRNAs are not clearly defined in many species. A number of studies have looked at the evolution of individual miRNA families. However, we currently lack resources for large-scale analysis of miRNA evolution. RESULTS: We addressed some of these issues in order to analyse the evolution of miRNAs. We perform syntenic and phylogenetic analysis for miRNAs from 80 animal species. We present synteny maps, phylogenies and functional data for miRNAs across these species. These data represent the basis of our analyses and also act as a resource for the community. CONCLUSIONS: We use these data to explore the distribution of miRNAs across phylogenetic space, characterise their birth and death, and examine functional relationships between miRNAs and other genes. These data confirm a number of previously reported findings on a larger scale and also offer novel insights into the evolution of the miRNA repertoire in animals, and it's genomic organization

Generation of a Cancer Testis Antigen Mcherry Reporter Hct116 Colorectal Carcinoma Cell Line

In the context of cancer immunotherapy, agents that target the immune system to cancer cells need to fulfil two criteria: 1) that they are only expressed on the desired target cell and 2) that they can elicit a potent immunological response. Cancer Testis Antigens are a large disparate family of factors ordinarily expressed in the germ-line but aberrantly expressed across multiple types of cancer. The ability to enforce their expression on tumour cells is an attractive strategy that could render such cells potent targets of the immune system, but very little is known about their regulation. We describe the generation of an mCherry reporter cell line using HCT116 colorectal carcinoma cells that we anticipate will be useful for screen-based approaches to identify novel regulators of CTA expression. Discoveries arising from their use could in future be exploited to enhance tumour cell immunogenicity and improve cancer immuno-therapy

GenomeChronicler: The Personal Genome Project UK Genomic Report Generator Pipeline

In recent years, there has been a significant increase in whole genome sequencing data of individual genomes produced by research projects as well as direct to consumer service providers. While many of these sources provide their users with an interpretation of the data, there is a lack of free, open tools for generating reports exploring the data in an easy to understand manner. GenomeChronicler was developed as part of the Personal Genome Project UK (PGP-UK) to address this need. PGP-UK provides genomic, transcriptomic, epigenomic and self-reported phenotypic data under an open-access model with full ethical approval. As a result, the reports generated by GenomeChronicler are intended for research purposes only and include information relating to potentially beneficial and potentially harmful variants, but without clinical curation. GenomeChronicler can be used with data from whole genome or whole exome sequencing, producing a genome report containing information on variant statistics, ancestry and known associated phenotypic traits. Example reports are available from the PGP-UK data page (personalgenomes.org.uk/data). The objective of this method is to leverage existing resources to find known phenotypes associated with the genotypes detected in each sample. The provided trait data is based primarily upon information available in SNPedia, but also collates data from ClinVar, GETevidence, and gnomAD to provide additional details on potential health implications, presence of genotype in other PGP participants and population frequency of each genotype. The analysis can be run in a self-contained environment without requiring internet access, making it a good choice for cases where privacy is essential or desired: any third party project can embed GenomeChronicler within their off-line safe-haven environments. GenomeChronicler can be run for one sample at a time, or in parallel making use of the Nextflow workflow manager. The source code is available from GitHub (https://github.com/PGP-UK/GenomeChronicler), container recipes are available for Docker and Singularity, as well as a pre-built container from SingularityHub (https://singularity-hub.org/collections/3664) enabling easy deployment in a variety of settings. Users without access to computational resources to run GenomeChronicler can access the software from the Lifebit CloudOS platform (https://lifebit.ai/cloudos) enabling the production of reports and variant calls from raw sequencing data in a scalable fashion

Identification of Recurrent Mutations in the microRNA-Binding Sites of B-Cell Lymphoma-Associated Genes in Follicular Lymphoma

Follicular lymphoma (FL) is a common indolent B-cell lymphoma that can transform into the more aggressive transformed FL (tFL). However, the molecular process driving this transformation is uncertain. In this work, we aimed to identify microRNA (miRNA)-binding sites recurrently mutated in follicular lymphoma patients, as well as in transformed FL patients. Using whole-genome sequencing data from FL tumors, we discovered 544 mutations located in bioinformatically predicted microRNA-binding sites. We then studied these specific regions using targeted sequencing in a cohort of 55 FL patients, found 16 recurrent mutations, and identified a further 69 variants. After filtering for QC, we identified 21 genes with mutated miRNA-binding sites that were also enriched for B-cell-associated genes by Gene Ontology. Over 40% of mutations identified in these genes were present exclusively in tFL patients. We validated the predicted miRNA-binding sites of five of the genes by luciferase assay and demonstrated that the identified mutations in BCL2 and EZH2 genes impaired the binding efficiency of miR-5008 and miR-144 and regulated the endogenous levels of messenger RNA (mRNA)

Genome-wide saturation mutagenesis of Burkholderia pseudomallei K96243 predicts essential genes and novel targets for antimicrobial development

This is the final version of the article. Available from the publisher via the DOI in this record.Burkholderia pseudomallei is the causative agent of melioidosis, an often fatal infectious disease for which there is no vaccine. B. pseudomallei is listed as a tier 1 select agent, and as current therapeutic options are limited due to its natural resistance to most antibiotics, the development of new antimicrobial therapies is imperative. To identify drug targets and better understand the complex B. pseudomallei genome, we sought a genome-wide approach to identify lethal gene targets. As B. pseudomallei has an unusually large genome spread over two chromosomes, an extensive screen was required to achieve a comprehensive analysis. Here we describe transposon-directed insertion site sequencing (TraDIS) of a library of over 10(6) transposon insertion mutants, which provides the level of genome saturation required to identify essential genes. Using this technique, we have identified a set of 505 genes that are predicted to be essential in B. pseudomallei K96243. To validate our screen, three genes predicted to be essential, pyrH, accA, and sodB, and a gene predicted to be nonessential, bpss0370, were independently investigated through the generation of conditional mutants. The conditional mutants confirmed the TraDIS predictions, showing that we have generated a list of genes predicted to be essential and demonstrating that this technique can be used to analyze complex genomes and thus be more widely applied. IMPORTANCE: Burkholderia pseudomallei is a lethal human pathogen that is considered a potential bioterrorism threat and has limited treatment options due to an unusually high natural resistance to most antibiotics. We have identified a set of genes that are required for bacterial growth and thus are excellent candidates against which to develop potential novel antibiotics. To validate our approach, we constructed four mutants in which gene expression can be turned on and off conditionally to confirm that these genes are required for the bacteria to survive.This work was partially funded by the Defence Science and Technology Laboratories (DSTL)

Post-vaccination COVID-19: A case-control study and genomic analysis of 119 breakthrough infections in partially vaccinated individuals

BACKGROUND: Post-vaccination infections challenge the control of the COVID-19 pandemic. METHODS: We matched 119 cases of post-vaccination SARS-CoV-2 infection with BNT162b2 mRNA, or ChAdOx1 nCOV-19, to 476 unvaccinated patients with COVID-19 (Sept 2020-March 2021), according to age and sex. Differences in 60-day all-cause mortality, hospital admission, and hospital length of stay were evaluated. Phylogenetic, single nucleotide polymorphism (SNP) and minority variant allele (MVA) full genome sequencing analysis was performed. RESULTS: 116/119 cases developed COVID-19 post first vaccination dose (median 14 days, IQR 9 - 24 days). Overall, 13/119 (10∙9%) cases and 158/476 (33∙2%) controls died (p<0.001), corresponding to 4∙5 number needed to treat (NNT). Multivariably, vaccination was associated with 69∙3% (95%CI 45∙8 - 82∙6) relative risk (RR) reduction in mortality. Similar results were seen in subgroup analysis for patients with infection onset ≥14 days after first vaccination (RR reduction 65∙1%, 95%CI 27∙2 - 83∙2, NNT 4∙5), and across vaccine subgroups (BNT162b2: RR reduction 66%, 95%CI 34∙9 - 82∙2, NNT 4∙7, ChAdOx1: RR reduction 78∙4%, 95%CI 30∙4 - 93∙3, NNT 4∙1). Hospital admissions (OR 0∙80, 95%CI 0∙51 - 1∙28), and length of stay (-1∙89 days, 95%CI -4∙57 - 0∙78) were lower for cases, while Ct values were higher (30∙8 versus 28∙8, p = 0.053). B.1.1.7 was the predominant lineage in cases (100/108, 92.6%) and controls (341/446, 76.5%). Genomic analysis identified one post-vaccination case harboring the E484K vaccine escape mutation (B.1.525 lineage). CONCLUSIONS: Previous vaccination reduces mortality when B.1.1.7 is the predominant lineage. No significant lineage-specific genomic changes during phylogenetic, SNP and MVA analysis were detected

A blood pressure-associated variant of the SLC39A8 gene influences cellular cadmium accumulation and toxicity.

Genome-wide association studies have revealed a relationship between inter-individual variation in blood pressure and the single nucleotide polymorphism rs13107325 in the SLC39A8 gene. This gene encodes the ZIP8 protein which co-transports divalent metal cations, including heavy metal cadmium, the accumulation of which has been associated with increased blood pressure. The polymorphism results in two variants of ZIP8 with either an alanine (Ala) or a threonine (Thr) at residue 391. We investigated the functional impact of this variant on protein conformation, cadmium transport, activation of signalling pathways and cell viability in relation to blood pressure regulation. Following incubation with cadmium, higher intracellular cadmium was detected in cultured human embryonic kidney cells (HEK293) expressing heterologous ZIP8-Ala391, compared with HEK293 cells expressing heterologous ZIP8-Thr391. This Ala391-associated cadmium accumulation also increased the phosphorylation of the signal transduction molecule ERK2, activation of the transcription factor NFκB, and reduced cell viability. Similarly, vascular endothelial cells with the Ala/Ala genotype had higher intracellular cadmium concentration and lower cell viability than their Ala/Thr counterpart following cadmium exposure. These results indicate that the ZIP8 Ala391-to-Thr391 substitution has an effect on intracellular cadmium accumulation and cell toxicity, providing a potential mechanistic explanation for the association of this genetic variant with blood pressure

Identification of Recurrent Mutations in the microRNA-Binding Sites of B-Cell Lymphoma-Associated Genes in Follicular Lymphoma.

Follicular lymphoma (FL) is a common indolent B-cell lymphoma that can transform into the more aggressive transformed FL (tFL). However, the molecular process driving this transformation is uncertain. In this work, we aimed to identify microRNA (miRNA)-binding sites recurrently mutated in follicular lymphoma patients, as well as in transformed FL patients. Using whole-genome sequencing data from FL tumors, we discovered 544 mutations located in bioinformatically predicted microRNA-binding sites. We then studied these specific regions using targeted sequencing in a cohort of 55 FL patients, found 16 recurrent mutations, and identified a further 69 variants. After filtering for QC, we identified 21 genes with mutated miRNA-binding sites that were also enriched for B-cell-associated genes by Gene Ontology. Over 40% of mutations identified in these genes were present exclusively in tFL patients. We validated the predicted miRNA-binding sites of five of the genes by luciferase assay and demonstrated that the identified mutations in BCL2 and EZH2 genes impaired the binding efficiency of miR-5008 and miR-144 and regulated the endogenous levels of messenger RNA (mRNA)