Global distribution of mcr gene variants in 214,095 metagenomic samples

Since the initial discovery of a mobilized colistin resistance gene (mcr-1), several other variants have been reported, some of which might have circulated a while before being discovered. Metagenomic data provides an opportunity to re-analyze available older data to understand the evolutionary history of recently discovered antimicrobial resistance genes (ARGs). Here, we present a large-scale metagenomic study of 442 Tbp of sequencing reads from 214,095 samples to identify the host and geographical distribution and genomic context of nine mcr gene variants (mcr-1 to mcr-9). Our results show that the dissemination of each variant is not uniform. Instead, the source and location play a role in the spread. Despite the very diverse distribution, the genomic background of the mcr genes remains unchanged as the same mobile genetic elements and plasmid replicons occur. This work emphasizes the importance of sharing genomic data for surveillance of ARGs in our fight against antimicrobial resistance.

Premature Termination Codon in 5′ Region of Desmoplakin and Plakoglobin Genes May Escape Nonsense-Mediated Decay through the Reinitiation of Translation

Arrhythmogenic cardiomyopathy is a heritable heart disease associated with desmosomal mutations, especially premature termination codon (PTC) variants. It is known that PTC triggers the nonsense-mediated decay (NMD) mechanism. It is also accepted that PTC in the last exon escapes NMD; however, the mechanisms involving NMD escaping in 5′-PTC, such as reinitiation of translation, are less known. The main objective of the present study is to evaluate the likelihood that desmosomal genes carrying 5′-PTC will trigger reinitiation. HL1 cell lines were edited by CRISPR/Cas9 to generate isogenic clones carrying 5′-PTC for each of the five desmosomal genes. The genomic context of the ATG in-frame in the 5′ region of desmosomal genes was evaluated by in silico predictions. The expression levels of the edited genes were assessed by Western blot and real-time PCR. Our results indicate that the 5′-PTC in PKP2, DSG2 and DSC2 acts as a null allele with no expression, whereas in the DSP and JUP gene, N-truncated protein is expressed. In concordance with this, the genomic context of the 5′-region of DSP and JUP presents an ATG in-frame with an optimal context for the reinitiation of translation. Thus, 5′-PTC triggers NMD in the PKP2, DSG2* and DSC2 genes, whereas it may escape NMD through the reinitiation of the translation in DSP and JUP genes, with no major effects on ACM-related gene expression.

MvPPT: a highly efficient and sensitive pathogenicity prediction tool for missense variants

Next generation sequencing technologies both boost the discovery of variants in the human genome and exacerbate the challenges of pathogenic variant identification. In this study, we developed mvPPT (Pathogenicity Prediction Tool for missense variants), a highly sensitive and accurate missense variant classifier based on gradient boosting. MvPPT adopts high-confidence training sets with a wide spectrum of variant profiles, and extracts three categories of features, including scores from existing prediction tools, allele, amino acid and genotype frequencies, and genomic context. Compared with established predictors, mvPPT achieved superior performance in all test sets, regardless of data source. In addition, our study also provides guidance for training set and feature selection strategies, as well as reveals highly relevant features, which may further provide biological insights of variant pathogenicity.

Genomic diversity of antimicrobial resistance in non-typhoidal Salmonella in Victoria, Australia

Non-typhoidal Salmonella (NTS) is the second most common cause of foodborne bacterial gastroenteritis in Australia with antimicrobial resistance (AMR) increasing in recent years. Whole-genome sequencing (WGS) provides opportunities for in silico detection of AMR determinants. The objectives of this study were two-fold: (1) establish the utility of WGS analyses for inferring phenotypic resistance in NTS, and (2) explore clinically relevant genotypic AMR profiles to third generation cephalosporins (3GC) in NTS lineages. The concordance of 2490 NTS isolates with matched WGS and phenotypic susceptibility data against 13 clinically relevant antimicrobials was explored. In silico serovar prediction and typing was performed on assembled reads and interrogated for known AMR determinants. The surrounding genomic context, plasmid determinants and co-occurring AMR patterns were further investigated for multidrug resistant serovars harbouring bla CMY-2, bla CTX-M-55 or bla CTX-M-65. Our data demonstrated a high correlation between WGS and phenotypic susceptibility testing. Phenotypic-genotypic concordance was observed between 2440/2490 (98.0 %) isolates, with overall sensitivity and specificity rates >98 % and positive and negative predictive values >97 %. The most common AMR determinants were bla TEM-1, sul2, tet(A), strA-strB and floR. Phenotypic resistance to cefotaxime and azithromycin was low and observed in 6.2 % (151/2486) and 0.9 % (16/1834) of the isolates, respectively. Several multi-drug resistant NTS lineages were resistant to 3GC due to different genetic mechanisms including bla CMY-2, bla CTX-M-55 or bla CTX-M-65. This study shows WGS can enhance existing AMR surveillance in NTS datasets routinely produced in public health laboratories to identify emerging AMR in NTS. These approaches will be critical for developing capacity to detect emerging public health threats such as resistance to 3GC.

Genome context influences evolutionary flexibility of nearly identical type III effectors in two phytopathogenic Pseudomonads

Integrative Conjugative Elements (ICEs) are replicons that can insert and excise from chromosomal locations in a site specific manner, can conjugate across strains, and which often carry a variety of genes useful for bacterial growth and survival under specific conditions. Although ICEs have been identified and vetted within certain clades of the agricultural pathogen Pseudomonas syringae, the impact of ICE carriage and transfer across the entire P. syringae species complex remains underexplored. Here we identify and vet an ICE (PmaICE-DQ) from P. syringae pv. maculicola ES4326, a strain commonly used for laboratory virulence experiments, demonstrate that this element can excise and conjugate across strains, and contains loci encoding multiple type III effector proteins. Moreover, genome context suggests that another ICE (PmaICE-AOAB) is highly similar in comparison with and found immediately adjacent to PmaICE-DQ within the chromosome of strain ES4326, and also contains multiple type III effectors. Lastly, we present passage data from in planta experiments that suggests that genomic plasticity associated with ICEs may enable strains to more rapidly lose type III effectors that trigger R-gene mediated resistance in comparison to strains where nearly isogenic effectors are not present in ICEs. Taken together, our study sheds light on a set of ICE elements from P. syringae pv. maculicola ES4326 and highlights how genomic context may lead to different evolutionary dynamics for shared virulence genes between strains.

Alu insertion variants alter gene transcript levels

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC. These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

A Multigraph model of the 3D genome

Spatial organisation of the genome has a fundamental effect on its biological functions. Chromatin-chromatin interactions and 3D spatial structures are involved in transcriptional regulation and have a decisive role in DNA replication and repair. To understand how individual genes and their regulatory elements function within the larger genomic context, and how the genome reacts as a whole to environmental stimuli, the linear sequence information needs to be interpreted in 3-dimensional space. While recent advances in chromatin conformation capture technologies including Hi-C, considerably advanced our understanding of the genomes, defining the DNA, as it is organized in the cell nucleus is still a challenging task. 3D genome modelling needs to reflect the DNA as a flexible polymer, which can wind up to the fraction of its total length and greatly unwind and stretch to implement a multitude of functions. Here we propose a novel approach to model genomes as a multigraph based on Hi-C contact data. Multigraph-based 3D genome modelling of barley and rice revealed the well-known Rabl and Rosetta chromatin organizations, respectively, as well as other higher order structures. Our results shows that the well-established toolset of Graph theory is highly valuable in modelling large genomes in 3D.

Whole genome sequencing facilitates intragenic variant interpretation following modifier screening in C. elegans

Background Intragenic modifiers (in-phase, second-site variants) are known to have dramatic effects on clinical outcomes, affecting disease attributes such as severity or age of onset. However, despite their clinical importance, the focus of many genetic screens in model systems is on the discovery of extragenic variants, with many labs still relying upon more traditional methods to identify modifiers. However, traditional methods such as PCR and Sanger sequencing can be time-intensive and do not permit a thorough understanding of the intragenic modifier effects in the context of non-isogenic genomic backgrounds. Results Here, we apply high throughput approaches to identify and understand intragenic modifiers using Caenorhabditis elegans. Specifically, we applied whole genome sequencing (WGS) to a mutagen-induced forward genetic screen to identify intragenic suppressors of a temperature-sensitive zyg-1(it25) allele in C. elegans. ZYG-1 is a polo kinase that is important for centriole function and cell divisions, and mutations that truncate its human orthologue, PLK4, have been associated with microcephaly. Combining WGS and CRISPR/Cas9, we rapidly identify intragenic modifiers, show that these variants are distributed non-randomly throughout zyg-1 and that genomic context plays an important role on phenotypic outcomes. Conclusions Ultimately, our work shows that WGS facilitates high-throughput identification of intragenic modifiers in clinically relevant genes by reducing hands-on research time and overall costs and by allowing thorough understanding of the intragenic phenotypic effects in the context of different genetic backgrounds.

Structural basis of terephthalate recognition by solute binding protein TphC

Biological degradation of Polyethylene terephthalate (PET) plastic and assimilation of the corresponding monomers ethylene glycol and terephthalate (TPA) into central metabolism offers an attractive route for bio-based molecular recycling and bioremediation applications. A key step is the cellular uptake of the non-permeable TPA into bacterial cells which has been shown to be dependent upon the presence of the key tphC gene. However, little is known from a biochemical and structural perspective about the encoded solute binding protein, TphC. Here, we report the biochemical and structural characterisation of TphC in both open and TPA-bound closed conformations. This analysis demonstrates the narrow ligand specificity of TphC towards aromatic para-substituted dicarboxylates, such as TPA and closely related analogues. Further phylogenetic and genomic context analysis of the tph genes reveals homologous operons as a genetic resource for future biotechnological and metabolic engineering efforts towards circular plastic bio-economy solutions.

Systematic analysis of intrinsic enhancer-promoter compatibility in the mouse genome.

Gene expression is in part controlled by cis-regulatory elements (CREs) such as enhancers and repressive elements. Anecdotal evidence has indicated that a CRE and a promoter need to be biochemically compatible for promoter regulation to occur, but this compatibility has remained poorly characterised in mammalian cells. We used high-throughput combinatorial reporter assays to test thousands of CRE - promoter pairs from three Mb-sized genomic regions in mouse cells. This revealed that CREs vary substantially in their promoter compatibility, ranging from striking specificity for a single promoter to quantitative differences in activation across a broad set of promoters. More than half of the tested CREs exhibit significant promoter selectivity. Housekeeping promoters tend to have similar CRE preferences, but other promoters exhibit a wide diversity of compatibilities. Higher-order TF motif combinations may account for compatibility. CRE-promoter selectivity does not correlate with looping interactions in the native genomic context, suggesting that chromatin folding and compatibility are two orthogonal mechanisms that confer specificity to gene regulation.