Identification of knowledge gaps in whole-genome sequence analysis of multi-resistant thermotolerant Campylobacter spp.

BACKGROUND: Campylobacter spp. is the most frequent cause of bacterial food-borne gastroenteritis and a high priority antibiotic resistant bacterium according to the World Health Organization (WHO). European monitoring of thermotolerant Campylobacter spp. does not reflect the global burden of resistances already circulating within the bacterial population worldwide. METHODS: We systematically compared whole genome sequencing with comprehensive phenotypic antimicrobial susceptibility, analyzing 494 thermotolerant Campylobacter poultry isolates from Vietnam and Germany. Any discrepancy was checked by repeating the wet lab and improving the dry lab part. Selected isolates were additionally analyzed via long-read Oxford Nanopore technology, leading to closed chromosomes and plasmids. RESULTS: Overall, 22 different resistance genes and gene variants (e. g. erm(B), aph(3')-IIIa, aph(2'')-If, catA, lnu(C), blaOXA, sat4) and point mutations in three distinct genes (gyrA, 23S rRNA, rpsL) associated with AMR were present in the Campylobacter isolates. Two AMR genes were missing in the database and one falsely associated with resistance. Bioinformatic analysis based on short-read data partly failed to identify tet(O) and aadE, when the genes were present as duplicate or homologous gene variants. Intriguingly, isolates also contained different determinants, redundantly conferring resistance to chloramphenicol, gentamicin, kanamycin, lincomycin and streptomycin. We found a novel tet(W) in tetracycline sensitive strains, harboring point mutations. Furthermore, analysis based on assemblies from short-read data was impaired to identify full length phase variable aad9, due to variations of the poly-C tract within the gene. The genetic determinant responsible for gentamicin resistance of one isolate from Germany could not be identified. GyrT86I, presenting the main determinant for (fluoro-)quinolone resistance led to a rare atypical phenotype of ciprofloxacin resistance but nalidixic acid sensitivity. Long-read sequencing predicted AMR genes were mainly located on the chromosome, and rarely on plasmids. Predictions from long- and short-read sequencing, respectively, often differed. AMR genes were often organized in multidrug resistance islands (MDRI) and partially located in proximity to transposase genes, suggesting main mobilization of resistance determinants is via natural transformation and transposition in Campylobacter. CONCLUSIONS: The results of this study suggest that there is frequent resistance gene duplication, mosaicism, and mutation leading to gene variation and truncation in Campylobacter strains that have not been reported in previous studies and are missing from databases. Furthermore, there is a need for deciphering yet unknown resistance mechanisms and resistance spread in thermotolerant Campylobacter spp. that may pose a challenge to global food safety.

Multiplex Real-Time PCR for the Detection of Tetracycline, Ciprofloxacin, and Erythromycin Resistance Determinants from Human and Foodborne Campylobacter jejuni and Campylobacter coli.

Campylobacter jejuni and Campylobacter coli are the predominant thermophilic species responsible for foodborne gastroenteritis worldwide. Elevated resistance to certain antibiotics was observed due to antimicrobial therapy in farm animals and humans, while reduced antimicrobial usage partially reduced antibiotic resistance. Monitoring the antimicrobial resistance demonstrated a substantial fraction of multi-resistant isolates, indicating the necessity of reliable tools for their detection. In this study, resistance determinants in 129 German and 21 Vietnamese isolates were selected to establish a novel multiplex real-time PCR (qPCR), facilitating the simultaneous detection of four resistance determinants. These comprised tet(O) gene variants associated with tetracycline resistance, point mutations GyrA_T86I and GyrA_T86V associated with ciprofloxacin resistance, and the erm(B) gene together with the point mutation A2075G in the 23S rRNA gene, associated with erythromycin resistance. Moreover, the performance of the qPCR assay was evaluated by comparing the results of qPCR to phenotypic antimicrobial resistance profiles, obtained with standardized EUCAMP3 microdilution panel, which showed 100% similarity (inclusivity and exclusivity). Variation in measurement methods, including qPCR machines and master mixes showed robustness, essential for laboratories. The assay can be used for the rapid detection of resistance determinants, and is beneficial for monitoring the spread of antibiotic resistance in C. jejuni and C. coli.

Genetic Characterization of Listeria from Food of Non-Animal Origin Products and from Producing and Processing Companies in Bavaria, Germany.

Reported cases of listeriosis from food of non-animal origin (FNAO) are increasing. In order to assess the risk of exposure to Listeria monocytogenes from FNAO, the genetic characterization of the pathogen in FNAO products and in primary production and processing plants needs to be investigated. For this, 123 samples of fresh and frozen soft fruit and 407 samples of 39 plants in Bavaria, Germany that produce and process FNAO were investigated for Listeria contamination. As a result, 64 Listeria spp. isolates were detected using ISO 11290-1:2017. Environmental swabs and water and food samples were investigated. L. seeligeri (36/64, 56.25%) was the most frequently identified species, followed by L. monocytogenes (8/64, 12.50%), L. innocua (8/64, 12.50%), L. ivanovii (6/64, 9.38%), L. newyorkensis (5/64, 7.81%), and L. grayi (1/64, 1.56%). Those isolates were subsequently sequenced by whole-genome sequencing and subjected to pangenome analysis to retrieve data on the genotype, serotype, antimicrobial resistance (AMR), and virulence markers. Eight out of sixty-four Listeria spp. isolates were identified as L. monocytogenes. The serogroup analysis detected that 62.5% of the L. monocytogenes isolates belonged to serogroup IIa (1/2a and 3a) and 37.5% to serogroup IVb (4b, 4d, and 4e). Furthermore, the MLST (multilocus sequence typing) analysis of the eight detected L. monocytogenes isolates identified seven different sequence types (STs) and clonal complexes (CCs), i.e., ST1/CC1, ST2/CC2, ST6/CC6, ST7/CC7, ST21/CC21, ST504/CC475, and ST1413/CC739. The core genome MLST analysis also showed high allelic differences and suggests plant-specific isolates. Regarding the AMR, we detected phenotypic resistance against benzylpenicillin, fosfomycin, and moxifloxacin in all eight L. monocytogenes isolates. Moreover, virulence factors, such as prfA, hly, plcA, plcB, hpt, actA, inlA, inlB, and mpl, were identified in pathogenic and nonpathogenic Listeria species. The significance of L. monocytogenes in FNAO is growing and should receive increasing levels of attention.

Whole-Genome Sequence Comparisons of Listeria monocytogenes Isolated from Meat and Fish Reveal High Inter- and Intra-Sample Diversity.

Interpretation of whole-genome sequencing (WGS) data for foodborne outbreak investigations is complex, as the genetic diversity within processing plants and transmission events need to be considered. In this study, we analyzed 92 food-associated Listeria monocytogenes isolates by WGS-based methods. We aimed to examine the genetic diversity within meat and fish production chains and to assess the applicability of suggested thresholds for clustering of potentially related isolates. Therefore, meat-associated isolates originating from the same samples or processing plants as well as fish-associated isolates were analyzed as distinct sets. In silico serogrouping, multilocus sequence typing (MLST), core genome MLST (cgMLST), and pangenome analysis were combined with screenings for prophages and genetic traits. Isolates of the same subtypes (cgMLST types (CTs) or MLST sequence types (STs)) were additionally compared by SNP calling. This revealed the occurrence of more than one CT within all three investigated plants and within two samples. Analysis of the fish set resulted in predominant assignment of isolates from pangasius catfish and salmon to ST2 and ST121, respectively, potentially indicating persistence within the respective production chains. The approach not only allowed the detection of distinct subtypes but also the determination of differences between closely related isolates, which need to be considered when interpreting WGS data for surveillance.

Oestrogen receptor ß regulates epigenetic patterns at specific genomic loci through interaction with thymine DNA glycosylase.

BACKGROUND: DNA methylation is one way to encode epigenetic information and plays a crucial role in regulating gene expression during embryonic development. DNA methylation marks are established by the DNA methyltransferases and, recently, a mechanism for active DNA demethylation has emerged involving the ten-eleven translocator proteins and thymine DNA glycosylase (TDG). However, so far it is not clear how these enzymes are recruited to, and regulate DNA methylation at, specific genomic loci. A number of studies imply that sequence-specific transcription factors are involved in targeting DNA methylation and demethylation processes. Oestrogen receptor beta (ERß) is a ligand-inducible transcription factor regulating gene expression in response to the female sex hormone oestrogen. Previously, we found that ERß deficiency results in changes in DNA methylation patterns at two gene promoters, implicating an involvement of ERß in DNA methylation. In this study, we set out to explore this involvement on a genome-wide level, and to investigate the underlying mechanisms of this function. RESULTS: Using reduced representation bisulfite sequencing, we compared genome-wide DNA methylation in mouse embryonic fibroblasts derived from wildtype and ERß knock-out mice, and identified around 8000 differentially methylated positions (DMPs). Validation and further characterisation of selected DMPs showed that differences in methylation correlated with changes in expression of the nearest gene. Additionally, re-introduction of ERß into the knock-out cells could reverse hypermethylation and reactivate expression of some of the genes. We also show that ERß is recruited to regions around hypermethylated DMPs. Finally, we demonstrate here that ERß interacts with TDG and that TDG binds ERß-dependently to hypermethylated DMPs. CONCLUSION: We provide evidence that ERß plays a role in regulating DNA methylation at specific genomic loci, likely as the result of its interaction with TDG at these regions. Our findings imply a novel function of ERß, beyond direct transcriptional control, in regulating DNA methylation at target genes. Further, they shed light on the question how DNA methylation is regulated at specific genomic loci by supporting a concept in which sequence-specific transcription factors can target factors that regulate DNA methylation patterns.

Whole genome and transcriptome analyses of environmental antibiotic sensitive and multi-resistant Pseudomonas aeruginosa isolates exposed to waste water and tap water.

The fitness of sensitive and resistant Pseudomonas aeruginosa in different aquatic environments depends on genetic capacities and transcriptional regulation. Therefore, an antibiotic-sensitive isolate PA30 and a multi-resistant isolate PA49 originating from waste waters were compared via whole genome and transcriptome Illumina sequencing after exposure to municipal waste water and tap water. A number of different genomic islands (e.g. PAGIs, PAPIs) were identified in the two environmental isolates beside the highly conserved core genome. Exposure to tap water and waste water exhibited similar transcriptional impacts on several gene clusters (antibiotic and metal resistance, genetic mobile elements, efflux pumps) in both environmental P. aeruginosa isolates. The MexCD-OprJ efflux pump was overexpressed in PA49 in response to waste water. The expression of resistance genes, genetic mobile elements in PA49 was independent from the water matrix. Consistently, the antibiotic sensitive strain PA30 did not show any difference in expression of the intrinsic resistance determinants and genetic mobile elements. Thus, the exposure of both isolates to polluted waste water and oligotrophic tap water resulted in similar expression profiles of mentioned genes. However, changes in environmental milieus resulted in rather unspecific transcriptional responses than selected and stimuli-specific gene regulation.

Leveraging cross-species transcription factor binding site patterns: from diabetes risk loci to disease mechanisms.

Genome-wide association studies have revealed numerous risk loci associated with diverse diseases. However, identification of disease-causing variants within association loci remains a major challenge. Divergence in gene expression due to cis-regulatory variants in noncoding regions is central to disease susceptibility. We show that integrative computational analysis of phylogenetic conservation with a complexity assessment of co-occurring transcription factor binding sites (TFBS) can identify cis-regulatory variants and elucidate their mechanistic role in disease. Analysis of established type 2 diabetes risk loci revealed a striking clustering of distinct homeobox TFBS. We identified the PRRX1 homeobox factor as a repressor of PPARG2 expression in adipose cells and demonstrate its adverse effect on lipid metabolism and systemic insulin sensitivity, dependent on the rs4684847 risk allele that triggers PRRX1 binding. Thus, cross-species conservation analysis at the level of co-occurring TFBS provides a valuable contribution to the translation of genetic association signals to disease-related molecular mechanisms.

Detecting and visualizing gene fusions.

In recent years, gene fusions have gained significant recognition as biomarkers. They can assist treatment decisions, are seldom found in normal tissue and are detectable through Next-generation sequencing (NGS) of the transcriptome (RNA-seq). To transform the data provided by the sequencer into robust gene fusion detection several analysis steps are needed. Usually the first step is to map the sequenced transcript fragments (RNA-seq) to a reference genome. One standard application of this approach is to estimate expression and detect variants within known genes, e.g. SNPs and indels. In case of gene fusions, however, completely novel gene structures have to be detected. Here, we describe the detection of such gene fusion events based on our comprehensive transcript annotation (ElDorado). To demonstrate the utility of our approach, we extract gene fusion candidates from eight breast cancer cell lines, which we compare to experimentally verified gene fusions. We discuss several gene fusion events, like BCAS3-BCAS4 that was only detected in the breast cancer cell line MCF7. As supporting evidence we show that gene fusions occur more frequently in copy number enriched regions (CNV analysis). In addition, we present the Transcriptome Viewer (TViewer) a tool that allows to interactively visualize gene fusions. Finally, we support detected gene fusions through literature mining based annotations and network analyses. In conclusion, we present a platform that allows detecting gene fusions and supporting them through literature knowledge as well as rich visualization capabilities. This enables scientists to better understand molecular processes, biological functions and disease associations, which will ultimately lead to better biomedical knowledge for the development of biomarkers for diagnostics and therapies.

RNA-sequencing as useful screening tool in the combat against the misuse of anabolic agents.

The abuse of anabolic substances in animal husbandry is forbidden within the EU and well controlled by detecting substance residues in different matrices. The application of newly designed drugs or substance cocktails represents big problems. Therefore developing sensitive test methods is important. The analysis of physiological changes caused by the use of anabolic agents on the molecular level, for example, by quantifying gene expression response, is a new approach to develop such screening methods. A novel technology for holistic gene expression analysis is RNA sequencing. In this study, the potential of this high-throughput method for the identification of biomarkers was evaluated. The effect of trenbolone acetate plus estradiol on gene expression in liver from Nguni heifers was analyzed with RNA sequencing. The expression of 40 selected candidate genes was verified via RT-qPCR, whereby 20 of these genes were significantly regulated. To extract the intended information from these regulated genes, biostatistical tools for pattern recognition were applied and resulted in a clear separation of the treatment groups. Those candidate genes could be verified in boars and in calves treated with anabolic substances. These results show the potential of RNA sequencing to screen for biomarker candidates to detect the abuse of anabolics. The verification of these biomarkers in boars and calves leads to the assumption that gene expression biomarkers are independent of breed or even species and that biomarkers, identified in farm animals could also act as potential biomarker candidates to detect the abuse of anabolic substances in human sports.

SFTA2--a novel secretory peptide highly expressed in the lung--is modulated by lipopolysaccharide but not hyperoxia.

Tissue-specific transcripts are likely to be of importance for the corresponding organ. While attempting to define the specific transcriptome of the human lung, we identified the transcript of a yet uncharacterized protein, SFTA2. In silico analyses, biochemical methods, fluorescence imaging and animal challenge experiments were employed to characterize SFTA2. Human SFTA2 is located on Chr. 6p21.33, a disease-susceptibility locus for diffuse panbronchiolitis. RT-PCR verified the abundance of SFTA2-specific transcripts in human and mouse lung. SFTA2 is synthesized as a hydrophilic precursor releasing a 59 amino acid mature peptide after cleavage of an N-terminal secretory signal. SFTA2 has no recognizable homology to other proteins while orthologues are present in all mammals. SFTA2 is a glycosylated protein and specifically expressed in nonciliated bronchiolar epithelium and type II pneumocytes. In accordance with other hydrophilic surfactant proteins, SFTA2 did not colocalize with lamellar bodies but colocalized with golgin97 and clathrin-labelled vesicles, suggesting a classical secretory pathway for its expression and secretion. In the mouse lung, Sfta2 was significantly downregulated after induction of an inflammatory reaction by intratracheal lipopolysaccharides paralleling surfactant proteins B and C but not D. Hyperoxia, however, did not alter SFTA2 mRNA levels. We have characterized SFTA2 and present it as a novel unique secretory peptide highly expressed in the lung.

Estrogen induces repression of the breast cancer and salivary gland expression gene in an estrogen receptor alpha-dependent manner.

The focus of this study is on the expression and regulation of the estrogen-regulated breast cancer and salivary gland expression (BASE) gene that may function as a breast cancer marker. In MCF7 cells, BASE is repressed by estrogen in an estrogen receptor alpha (ER alpha)-dependent manner. Promoter analysis of the BASE gene led to the identification of a 2-kb upstream enhancer that harbors binding sites for ER alpha and FoxA1. The recruitment of both ER alpha and FoxA1 to this region was shown by chromatin immunoprecipitation analysis. Furthermore, mutation studies and knockdown experiments show a clear separation between gene expression mediated by FoxA1 and ER alpha-dependent gene regulation. Additionally, we provide information on BASE expression in human breast tumor samples.

E2-mediated cathepsin D (CTSD) activation involves looping of distal enhancer elements.

Estrogen receptor alpha (ERalpha) is a ligand dependent transcription factor that regulates the expression of target genes through interacting with cis-acting estrogen response elements (EREs). However, only a minority of ERalpha binding sites are located within the proximal promoter regions of responsive genes. Here we report the characterization of an ERE located 9kbp upstream of the TSS of the cathepsin D gene (CTSD) that up-regulates CTSD expression upon estrogen stimulation in MCF-7 cells. Using ChIP, we show recruitment of ERalpha and phosphorylated PolII at the CTSD distal enhancer region. Moreover, we determine the kinetics of transient CpG methylation on the promoter region of CTSD and for the first time, at a distal enhancer element. We show that ERalpha is crucial for long-distance regulation of CTSD expression involving a looping mechanism.

Yeast adapt to near-freezing temperatures by STRE/Msn2,4-dependent induction of trehalose synthesis and certain molecular chaperones.

Virtually nothing is known about the biochemical adaptations in eukaryotic cells that may enhance survival at low temperatures or upon freezing. Here we demonstrate an adaptive response in yeast that is activated below 10 degrees C and increases tolerance to low temperatures and freezing. This response involves a dramatic accumulation of the chemical chaperone trehalose and induction of trehalose-synthesizing enzymes (Tps1, Tps2) and certain heat shock proteins (Hsp104, Hsp42, Hsp12, Ssa4). mRNAs for these proteins increase dramatically below 10 degrees C and even at 0 degrees C. Their expression requires Msn2,4 transcription factors but also involves marked mRNA stabilization. Upon return to 30 degrees C, TPS1, TPS2, and HSP104 mRNAs, trehalose levels and tolerance to freezing fall dramatically within minutes. Mutants lacking trehalose or Msn2,4 die more rapidly at 0 degrees C and upon freezing. Thus, below 10 degrees C, yeast show an adaptive response that sustains viability at low or freezing temperatures, which are commonly encountered in natural environments and laboratory refrigerators.

Genome-wide identification of high-affinity estrogen response elements in human and mouse.

Although estrogen receptors (ERs) recognize 15-bp palindromic estrogen response elements (EREs) with maximal affinity in vitro, few near-consensus sequences have been characterized in estrogen target genes. Here we report the design of a genome-wide screen for high-affinity EREs and the identification of approximately 70000 motifs in the human and mouse genomes. EREs are enriched in regions proximal to the transcriptional start sites, and approximately 1% of elements appear conserved in the flanking regions (-10 kb to +5 kb) of orthologous human and mouse genes. Conserved and nonconserved elements were also found, often in multiple occurrences, in more than 230 estrogen-stimulated human genes previously identified from expression studies. In genes containing known EREs, we also identified additional distal elements, sometimes with higher in vitro binding affinity and/or better conservation between the species considered. Chromatin immunoprecipitation experiments in breast cancer cell lines indicate that most novel elements present in responsive genes bind ERalpha in vivo, including some EREs located up to approximately 10 kb from transcriptional start sites. Our results demonstrate that near-consensus EREs occur frequently in both genomes and that whereas chromatin structure likely modulates access to binding sites, far upstream elements can be evolutionarily conserved and bind ERs in vivo.