Direct RNA Nanopore Sequencing of Pseudomonas aeruginosa Clone C Transcriptomes.

The transcriptomes of Pseudomonas aeruginosa clone C isolates NN2 and SG17M during the mid-exponential and early stationary phases of planktonic growth were evaluated by direct RNA sequencing on the nanopore platform and compared with established short-read cDNA sequencing on the Illumina platform. Fifty to ninety percent of the sense RNAs turned out to be rRNA molecules, followed by similar proportions of mRNA transcripts and noncoding RNAs. The two platforms detected similar proportions of uncharged tRNAs and 29 yet-undescribed antisense tRNAs. For example, the rarest arginine codon was paired with the most abundant tRNAArg, and the tRNAArg gene is missing for the most frequent arginine codon. More than 90% of the antisense RNA molecules were complementary to a coding sequence. The antisense RNAs were evenly distributed in the genomes. Direct RNA sequencing identified more than 4,000 distinct nonoverlapping antisense RNAs during exponential and stationary growth. Besides highly expressed small antisense RNAs less than 200 bases in size, a population of longer antisense RNAs was sequenced that covered a broad range (a few hundred to thousands of bases) and could be complementary to a contig of several genes. In summary, direct RNA sequencing identified yet-undescribed RNA molecules and an unexpected composition of the pools of tRNAs and sense and antisense RNAs. IMPORTANCE Genome-wide gene expression of bacteria is commonly studied by high-throughput sequencing of size-selected cDNA fragment libraries of reverse-transcribed RNA preparations. However, the depletion of rRNAs, enzymatic reverse transcription, and the fragmentation, size selection, and amplification during library preparation lead to inevitable losses of information about the initial composition of the RNA pool. We demonstrate that direct RNA sequencing on the Nanopore platform can overcome these limitations. Nanopore sequencing of total RNA yielded novel insights into the Pseudomonas aeruginosa transcriptome that-if replicated in other species-will change our view of the bacterial RNA world. The discovery of sense-antisense pairs of transfer-messenger RNA (tmRNA), tRNAs, and mRNAs indicates a further and unknown level of gene regulation in bacteria.

Wochenende - modular and flexible alignment-based shotgun metagenome analysis.

BACKGROUND: Shotgun metagenome analysis provides a robust and verifiable method for comprehensive microbiome analysis of fungal, viral, archaeal and bacterial taxonomy, particularly with regard to visualization of read mapping location, normalization options, growth dynamics and functional gene repertoires. Current read classification tools use non-standard output formats, or do not fully show information on mapping location. As reference datasets are not perfect, portrayal of mapping information is critical for judging results effectively. RESULTS: Our alignment-based pipeline, Wochenende, incorporates flexible quality control, trimming, mapping, various filters and normalization. Results are completely transparent and filters can be adjusted by the user. We observe stringent filtering of mismatches and use of mapping quality sharply reduces the number of false positives. Further modules allow genomic visualization and the calculation of growth rates, as well as integration and subsequent plotting of pipeline results as heatmaps or heat trees. Our novel normalization approach additionally allows calculation of absolute abundance profiles by comparison with reads assigned to the human host genome. CONCLUSION: Wochenende has the ability to find and filter alignments to all kingdoms of life using both short and long reads, and requires only good quality reference genomes. Wochenende automatically combines multiple available modules ranging from quality control and normalization to taxonomic visualization. Wochenende is available at https://github.com/MHH-RCUG/nf_wochenende .

Reprogramming enriches for somatic cell clones with small-scale mutations in cancer-associated genes.

Cellular therapies based on induced pluripotent stem cells (iPSCs) come out of age and an increasing number of clinical trials applying iPSC-based transplants are ongoing or in preparation. Recent studies, however, demonstrated a high number of small-scale mutations in iPSCs. Although the mutational load in iPSCs seems to be largely derived from their parental cells, it is still unknown whether reprogramming may enrich for individual mutations that could lead to loss of functionality and tumor formation from iPSC derivatives. 30 hiPSC lines were analyzed by whole exome sequencing. High accuracy amplicon sequencing showed that all analyzed small-scale variants pre-existed in their parental cells and that individual mutations present in small subpopulations of parental cells become enriched among hiPSC clones during reprogramming. Among those, putatively actionable driver mutations affect genes related to cell-cycle control, cell death, and pluripotency and may confer a selective advantage during reprogramming. Finally, a short hairpin RNA (shRNA)-based experimental approach was applied to provide additional evidence for the individual impact of such genes on the reprogramming efficiency. In conclusion, we show that enriched mutations in curated onco- and tumor suppressor genes may account for an increased tumor risk and impact the clinical value of patient-derived hiPSCs.

Genome-Wide Methylation Mapping Using Nanopore Sequencing Technology Identifies Novel Tumor Suppressor Genes in Hepatocellular Carcinoma.

Downregulation of multiple tumor suppressor genes (TSGs) plays an important role in cancer formation. Recent evidence has accumulated that cancer progression involves genome-wide alteration of epigenetic modifications, which may cause downregulation of the tumor suppressor gene. Using hepatocellular carcinoma (HCC) as a system, we mapped 5-methylcytosine signal at a genome-wide scale using nanopore sequencing technology to identify novel TSGs. Integration of methylation data with gene transcription profile of regenerated liver and primary HCCs allowed us to identify 10 potential tumor suppressor gene candidates. Subsequent validation led us to focus on functionally characterizing one candidate-glucokinase (GCK). We show here that overexpression of GCK inhibits the proliferation of HCC cells via induction of intracellular lactate accumulation and subsequently causes energy crisis due to NAD+ depletion. This suggests GCK functions as a tumor suppressor gene and may be involved in HCC development. In conclusion, these data provide valuable clues for further investigations of the process of tumorigenesis in human cancer.

Long-Term Microevolution of Pseudomonas aeruginosa Differs between Mildly and Severely Affected Cystic Fibrosis Lungs.

Chronic airway infections with Pseudomonas aeruginosa determine morbidity in most individuals with cystic fibrosis (CF). P. aeruginosa may persist for decades in CF lungs, which provides a rare opportunity to study the long-term within-host evolution of a bacterial airway pathogen. In this work, we sought to resolve the genetic adaptation of P. aeruginosa in CF lungs from the onset of colonization until the patient's death or permanent replacement by another P. aeruginosa clone. We followed the microevolution of the first persisting P. aeruginosa clone by whole-genome sequencing of serial isolates from highly divergent clinical courses of airway infection, i.e., a fatal outcome because of respiratory insufficiency within less than 15 years, or a rather normal daily life 25-35 years after acquisition of P. aeruginosa. Nonneutral mutations predominantly emerged in P. aeruginosa genes relevant for protection against and communication with signals from the lung environment, i.e., antibiotic resistance, cell wall components, and two-component systems. Drastic and loss-of-function mutations preferentially happened during the severe courses of infection, and the bacterial lineages of the mild courses more proficiently incorporated extra metabolic genes into their accessory genome. P. aeruginosa followed different evolutionary paths depending on whether the bacterium had taken up residence in a patient with CF and normal or already compromised lung function.

Functional characterization of the mouse Serpina1 paralog DOM-7.

The generation of authentic mouse-models for human α1-antitrypsin (A1AT)-deficiency is difficult due to the high complexity of the mouse Serpina1 gene locus. Depending on the exact mouse strain, three to five paralogs are expressed, with different proteinase inhibitory properties. Nowadays with CRISPR-technology, genome editing of complex genomic loci is feasible and could be employed for the generation of A1AT-deficiency mouse models. In preparation of a CRISPR/Cas9-based genome-engineering approach we identified cDNA clones with a functional CDS for the Serpina1-paralog DOM-7. Here, we show that DOM-7 functionally inhibits neutrophil elastase (ELANE) and chymotrypsin, and therefore needs to be considered when aiming at the generation of A1AT-deficient models.

Impact of Molecular Genetics on Outcome in Myelofibrosis Patients after Allogeneic Stem Cell Transplantation.

Molecular genetics may influence outcome for patients with myelofibrosis. To determine the impact of molecular genetics on outcome after allogeneic stem cell transplantation, we screened 169 patients with primary myelofibrosis (n = 110), post-essential thrombocythemia/polycythemia vera myelofibrosis (n = 46), and myelofibrosis in transformation (n = 13) for mutations in 16 frequently mutated genes. The most frequent mutation was JAK2V617F (n = 101), followed by ASXL1 (n = 49), calreticulin (n = 34), SRSF2 (n = 16), TET2 (n = 10), U2AF1 (n = 11), EZH2 (n = 7), MPL (n = 6), IDH2 (n = 5), IDH1 (n = 4), and CBL (n = 1). The cumulative incidence of nonrelapse mortality (NRM) at 1 year was 21% and of relapse at 5 years 25%. The 5-year rates progression-free (PFS) and overall survival (OS) were and 56%, respectively. In a multivariate analysis CALR mutation was an independent factor for lower NRM (HR, .415; P = .05), improved PFS (HR, .393; P = .01), and OS (HR, .448; P = .03). ASXL1 and IDH2 mutations were independent risk factors for lower PFS (HR, 1.53 [P = .008], and HR, 5.451 [P = .002], respectively), whereas no impact was observed for "triple negative" patients. Molecular genetics, especially CALR, IDH2, and ASXL1 mutations, may thus be useful to predict outcome independently from known clinical risk factors after allogeneic stem cell transplantation for myelofibrosis.

Molecular Epidemiology of Mutations in Antimicrobial Resistance Loci of Pseudomonas aeruginosa Isolates from Airways of Cystic Fibrosis Patients.

The chronic airway infections with Pseudomonas aeruginosa in people with cystic fibrosis (CF) are treated with aerosolized antibiotics, oral fluoroquinolones, and/or intravenous combination therapy with aminoglycosides and ß-lactam antibiotics. An international strain collection of 361 P. aeruginosa isolates from 258 CF patients seen at 30 CF clinics was examined for mutations in 17 antimicrobial susceptibility and resistance loci that had been identified as hot spots of mutation by genome sequencing of serial isolates from a single CF clinic. Combinatorial amplicon sequencing of pooled PCR products identified 1,112 sequence variants that were not present in the genomes of representative strains of the 20 most common clones of the global P. aeruginosa population. A high frequency of singular coding variants was seen in spuE, mexA, gyrA, rpoB, fusA1, mexZ, mexY, oprD, ampD, parR, parS, and envZ (amgS), reflecting the pressure upon P. aeruginosa in lungs of CF patients to generate novel protein variants. The proportion of nonneutral amino acid exchanges was high. Of the 17 loci, mexA, mexZ, and pagL were most frequently affected by independent stop mutations. Private and de novo mutations seem to play a pivotal role in the response of P. aeruginosa populations to the antimicrobial load and the individual CF host.

Mutations in the cohesin complex in acute myeloid leukemia: clinical and prognostic implications.

Mutations in the cohesin complex are novel, genetic lesions in acute myeloid leukemia (AML) that are not well characterized. In this study, we analyzed the frequency, clinical, and prognostic implications of mutations in STAG1, STAG2, SMC1A, SMC3, and RAD21, all members of the cohesin complex, in a cohort of 389 uniformly treated AML patients by next generation sequencing. We identified a total of 23 patients (5.9%) with somatic mutations in 1 of the cohesin genes. All gene mutations were mutually exclusive, and STAG1 (1.8%), STAG2 (1.3%), and SMC3 (1.3%) were most frequently mutated. Patients with any cohesin complex mutation had lower BAALC expression levels. We found a strong association between mutations affecting the cohesin complex and NPM1. Mutated allele frequencies were similar between NPM1 and cohesin gene mutations. Overall survival (OS), relapse-free survival (RFS), and complete remission rates (CR) were not influenced by the presence of cohesin mutations (OS: hazard ratio [HR] 0.98; 95% confidence interval [CI], 0.56-1.72 [P = .94]; RFS: HR 0.7; 95% CI, 0.36-1.38 [P = .3]; CR: mutated 83% vs wild-type 76% [P = .45]). The cohesin complex presents a novel pathway affected by recurrent mutations in AML. This study is registered at www.clinicaltrials.gov as #NCT00209833.

THOC5 controls 3'end-processing of immediate early genes via interaction with polyadenylation specific factor 100 (CPSF100).

Transcription of immediate early genes (IEGs) in response to extrinsic and intrinsic signals is tightly regulated at multiple stages. It is known that untranslated regions of the RNA can play a role in these processes. Here we show that THOC5, a member of the TREX (transcription/export) complex, plays a role in expression of only a subset of constitutively active genes, however transcriptome analysis reveals that more than 90% of IEG were not induced by serum in THOC5 depleted cells. Furthermore, THOC5 depletion does not influence the expression of the most rapidly induced IEGs, e.g. Fos and Jun. One group of THOC5 target genes, including Id1, Id3 and Wnt11 transcripts, were not released from chromatin in THOC5 depleted cells. Genes in another group, including Myc and Smad7 transcripts, were released with shortening of 3'UTR by alternative cleavage, and were spliced but export was impaired in THOC5 depleted cells. By interactome analysis using THOC5 as bait, we show that upon stimulation with serum THOC5 forms a complex with polyadenylation-specific factor 100 (CPSF100). THOC5 is required for recruitment of CPSF100 to 3'UTR of THOC5 target genes. These data suggest the presence of a novel mechanism for the control of IEG response by THOC5 via 3'end-processing.

Interclonal gradient of virulence in the Pseudomonas aeruginosa pangenome from disease and environment.

The population genomics of Pseudomonas aeruginosa was analysed by genome sequencing of representative strains of the 15 most frequent clonal complexes in the P. aeruginosa population and of the five most common clones from the environment of which so far no isolate from a human infection has been detected. Gene annotation identified 5892-7187 open reading frame (ORFs; median 6381 ORFs) in the 20 6.4-7.4 Mbp large genomes. The P. aeruginosa pangenome consists of a conserved core of at least 4000 genes, a combinatorial accessory genome of a further 10 000 genes and 30 000 or more rare genes that are present in only a few strains or clonal complexes. Whole genome comparisons of single nucleotide polymorphism synteny indicated unrestricted gene flow between clonal complexes by recombination. Using standardized acute lettuce, Galleria mellonella and murine airway infection models the full spectrum of possible host responses to P. aeruginosa was observed with the 20 strains ranging from unimpaired health following infection to 100% lethality. Genome comparisons indicate that the differential genetic repertoire of clones maintains a habitat-independent gradient of virulence in the P. aeruginosa population.

Sputum metagenomics of people with bronchiectasis.

Background: The microbiota in the sputum of people with bronchiectasis has repeatedly been investigated in cohorts of different geographic origin, but so far has not been studied to the species level in comparison to control populations including healthy adults and smokers without lung disease. Methods: The microbial metagenome from sputa of 101 European Bronchiectasis Registry (EMBARC) study participants was examined by using whole-genome shotgun sequencing. Results: Our analysis of the metagenome of people with bronchiectasis revealed four clusters characterised by a predominance of Haemophilus influenzae, Pseudomonas aeruginosa or polymicrobial communities with varying compositions of nonpathogenic commensals and opportunistic pathogens. The metagenomes of the severely affected patients showed individual profiles characterised by low alpha diversity. Importantly, nearly 50% of patients with severe disease were grouped in a cluster characterised by commensals. Comparisons with the sputum metagenomes of healthy smokers and healthy nonsmokers revealed a gradient of depletion of taxa in bronchiectasis, most often Neisseria subflava, Fusobacterium periodonticum and Eubacterium sulci. Conclusion: The gradient of depletion of commensal taxa found in healthy airways is a key feature of bronchiectasis associated with disease severity.

THOC5, a member of the mRNA export complex, contributes to processing of a subset of wingless/integrated (Wnt) target mRNAs and integrity of the gut epithelial barrier.

BACKGROUND: THO (Suppressors of the transcriptional defects of hpr1 delta by overexpression) complex 5 (THOC5), an mRNA export protein, is involved in the expression of only 1% of all genes. Using an interferon inducible knockout mouse system, we have previously shown that THOC5 is an essential element in the maintenance of hematopoietic stem cells and cytokine-mediated hematopoiesis in adult mice. Here we interrogate THOC5 function in cell differentiation beyond the hematopoietic system and study pathological changes caused by THOC5 deficiency. RESULTS: To examine whether THOC5 plays a role in general differentiation processes, we generated tamoxifen inducible THOC5 knockout mice. We show here that the depletion of THOC5 impaired not only hematopoietic differentiation, but also differentiation and self renewal of the gut epithelium. Depletion of the THOC5 gene did not cause pathological alterations in liver or kidney. We further show that THOC5 is indispensable for processing of mRNAs induced by Wnt (wingless/integrated) signaling which play key roles in epithelial cell differentiation/proliferation. A subset of Wnt target mRNAs, SRY-box containing gene 9 (Sox9), and achaete-scute complex homolog 2 (Ascl2), but not Fibronectin 1 (Fn1), were down-regulated in THOC5 knockout intestinal cells. The down-regulated Wnt target mRNAs were able to bind to THOC5. Furthermore, pathological alterations in the gastrointestinal tract induced translocation of intestinal bacteria and caused sepsis in mice. The bacteria translocation may cause Toll-like receptor activation. We identified one of the Toll-like receptor inducible genes, prostaglandin-endoperoxidase synthase 2 (Ptgs2 or COX2) transcript as THOC5 target mRNA. CONCLUSION: THOC5 is indispensable for processing of only a subset of mRNAs, but plays a key role in processing of mRNAs inducible by Wnt signals. Furthermore, THOC5 is dispensable for general mRNA export in terminally differentiated organs, indicating that multiple mRNA export pathways exist. These data imply that THOC5 may be a useful tool for studying intestinal stem cells, for modifying the differentiation processes and for cancer therapy.

Intraclonal diversity of the Pseudomonas aeruginosa cystic fibrosis airway isolates TBCF10839 and TBCF121838: distinct signatures of transcriptome, proteome, metabolome, adherence and pathogenicity despite an almost identical genome sequence.

Microevolution of closely related Pseudomonas aeruginosa was compared in the clone TB strains TBCF10839 and TBCF121838 which had been isolated from two unrelated individuals with cystic fibrosis who had acquired clone TB during a local outbreak. Compared with the strain PAO1 reference sequence the two clone TB genomes shared 23 155 nucleotide exchanges, 32 out-of-frame indels in the coding region and another repertoire of replacement and genomic islands such as PAGI-1, PAGI-2, PAGI-5, LESGI-1 and LES-prophage 4. Only TBCF121838 carried a genomic island known from Ralstonia pickettii. Six of the seven strain-specific sequence variations in the core genome were detected in genes affecting motility, biofilm formation or virulence, i.e. non-synonymous nucleotide substitutions in mexS, PA3729, PA5017, mifR, a frameshift mutation in pilF (TBCF121838) and an intragenic deletion in pilQ (TBCF10839). Despite their almost identical genome sequence the two strains differed strongly from each other in transcriptome and metabolome profiles, mucin adherence and phagocytosis assays. TBCF121838 was susceptible to killing by neutrophils, but TBCF10839 could grow in leucocytes. Microevolution in P. aeruginosa apparently can generate novel complex traits by few or even single mutations provided that predisposing mutational events had occurred before in the clonal lineage.

Longitudinal development of the airway metagenome of preterm very low birth weight infants during the first two years of life.

Preterm birth is accompanied with many complications and requires severe therapeutic regimens at the neonatal intensive care unit. The influence of the above-mentioned factors on the premature-born infants' respiratory metagenome or more generally its maturation is unknown. We therefore applied shotgun metagenome sequencing of oropharyngeal swabs to analyze the airway metagenome development of 24 preterm infants from one week postpartum to 15 months of age. Beta diversity analysis revealed a distinct clustering of airway microbial communities from hospitalized preterms and samples after hospital discharge. At nine and 15 months of age, the preterm infants lost their hospital-acquired individual metagenome signatures towards a common taxonomic structure. However, ecological network analysis and Random Forest classification of cross-sectional data revealed that by this age the preterm infants did not succeed in establishing the uniform and stable bacterial community structures that are characteristic for healthy full-term infants.

A VirB4 ATPase of the mobile accessory genome orchestrates core genome-encoded features of physiology, metabolism, and virulence of Pseudomonas aeruginosa TBCF10839.

Pseudomonas aeruginosa TBCF10839 is a highly virulent strain that can persist and replicate in human neutrophils. Screening of a signature-tagged mutagenesis (STM) TBCF10839 transposon library in phagocytosis tests identified a mutant that carried the transposon in the VirB4 homolog 5PG21 of an integrative and conjugative element (ICE)-associated type IV secretion system of the pKLC102 subtype. 5P21 TBCF10839 insertion mutants were deficient in metabolic versatility, secretion, quorum sensing, and virulence. The mutants were efficiently killed in phagocytosis tests in vitro and were avirulent in an acute murine airway infection model in vivo. The inactivation of 5PG21 silenced the rhl, las, and pqs operons and the gene expression for the synthesis of hydrogen cyanide, the antimetabolite l-2-amino-4-methoxy-trans-3-butenoic acid, and the H2- and H3-type VI secretion systems and their associated effectors. The mutants were impaired in the utilization of carbon sources and stored compounds that are not funneled into intermediary metabolism. This showcase demonstrates that a single gene of the mobile accessory genome can become an essential element to operate the core genome-encoded features of metabolism and virulence.

THOC5 complexes with DDX5, DDX17, and CDK12 to regulate R loop structures and transcription elongation rate.

THOC5, a member of the THO complex, is essential for the 3'processing of some inducible genes, the export of a subset of mRNAs and stem cell survival. Here we show that THOC5 depletion results in altered 3'cleavage of >50% of mRNAs and changes in RNA polymerase II binding across genes. THOC5 is recruited close to high-density polymerase II sites, suggesting that THOC5 is involved in transcriptional elongation. Indeed, measurement of elongation rates in vivo demonstrated decreased rates in THOC5-depleted cells. Furthermore, THOC5 is preferentially recruited to its target genes in slow polymerase II cells compared with fast polymerase II cells. Importantly chromatin-associated THOC5 interacts with CDK12 (a modulator of transcription elongation) and RNA helicases DDX5, DDX17, and THOC6 only in slow polymerase II cells. The CDK12/THOC5 interaction promotes CDK12 recruitment to R-loops in a THOC6-dependent manner. These data demonstrate a novel function of THOC5 in transcription elongation.

The Cystic Fibrosis Upper and Lower Airway Metagenome.

The microbial metagenome in cystic fibrosis (CF) airways was investigated by whole-genome shotgun sequencing of total DNA isolated from nasal lavage samples, oropharyngeal swabs, and induced sputum samples collected from 65 individuals with CF aged 7 to 50 years. Each patient harbored a personalized microbial metagenome unique in microbial load and composition, the exception being monocultures of the most common CF pathogens Staphylococcus aureus and Pseudomonas aeruginosa from patients with advanced lung disease. The sampling of the upper airways by nasal lavage uncovered the fungus Malassezia restricta and the bacterium Staphylococcus epidermidis as prominent species. Healthy and CF donors harbored qualitatively and quantitatively different spectra of commensal bacteria in their sputa, even in the absence of any typical CF pathogen. If P. aeruginosa, S. aureus, or Stenotrophomonas maltophilia belonged to the trio of the most abundant species in the CF sputum metagenome, common inhabitants of the respiratory tract of healthy subjects, i.e., Eubacterium sulci, Fusobacterium periodonticum, and Neisseria subflava, were present only in low numbers or not detectable. Random forest analysis identified the numerical ecological parameters of the bacterial community, such as Shannon and Simpson diversity, as the key parameters that globally distinguish sputum samples from CF and healthy donors. IMPORTANCE Cystic fibrosis (CF) is the most common life-limiting monogenetic disease in European populations and is caused by mutations in the CFTR gene. Chronic airway infections with opportunistic pathogens are the major morbidity that determines prognosis and quality of life in most people with CF. We examined the composition of the microbial communities of the oral cavity and upper and lower airways in CF patients across all age groups. From early on, the spectrum of commensals is different in health and CF. Later on, when the common CF pathogens take up residence in the lungs, we observed differential modes of depletion of the commensal microbiota in the presence of S. aureus, P. aeruginosa, S. maltophilia, or combinations thereof. It remains to be seen whether the implementation of lifelong CFTR (cystic fibrosis transmembrane conductance regulator) modulation will change the temporal evolution of the CF airway metagenome.

Comparison of methylation estimates obtained via MinION nanopore sequencing and sanger bisulfite sequencing in the TRPA1 promoter region.

BACKGROUND: Bisulfite sequencing has long been considered the gold standard for measuring DNA methylation at single CpG resolution. However, in recent years several new approaches like nanopore sequencing have been developed due to hints for a partial error-proneness of bisulfite sequencing. Since these errors were shown to be sequence-specific, we aimed to verify the methylation data of a particular region of the TRPA1 promoter from our previous studies obtained by bisulfite sequencing. METHODS: We compared methylation rates determined by direct bisulfite sequencing and nanopore sequencing following Cas9-mediated PCR-free enrichment. RESULTS: We could show that CpG methylation levels above 20% corroborate well with our previous data. Within the range between 0 and 20% methylation, however, Sanger sequencing data have to be interpreted cautiously, at least in the investigated region of interest (TRPA1 promotor region). CONCLUSION: Based on the investigation of the TRPA1- region as an example, the present work can help in choosing the right method out of the two current main approaches for methylation analysis for different individual settings regarding many factors like cohort size, costs and prerequisites that should be fulfilled for each method. All in all, both methods have their raison d'être. Furthermore, the present paper contains and illustrates some important basic information and explanation of how guide RNAs should be located for an optimal outcome in Cas9 mediated PCR free target enrichment.

Pseudomonas aeruginosa population structure revisited under environmental focus: impact of water quality and phage pressure.

Pseudomonas aeruginosa attracts research attention as a common opportunistic nosocomial pathogen causing severe health problems in humans. Nevertheless, its primary habitat is the natural environment. Here, we relate the genetic diversity of 381 environmental isolates from rivers in northern Germany to ecological factors such as river system, season of sampling and different levels of water quality. From representatives of 99 environmental clones, also in comparison with 91 clinical isolates, we determined motility phenotypes, virulence factors, biofilm formation, serotype and the resistance to seven environmental P.aeruginosa phages. The integration of genetic, ecological and phenotypic data showed (i) the presence of several extended clonal complexes (ecc) which are non-uniformly distributed across different water qualities, and (ii) a correlation of the hosts' serotype composition with susceptibility towards distinct groups of environmental phages. For at least one ecc (eccB), we assumed the ecophysiological differences on environmental water adaptation and phage resistance to be so distinct as to reinforce an environmentally driven cladogenic split from the remainder of P.aeruginosa. In summary, we conclude that the majority of the microevolutionary population dynamics of P.aeruginosa were shaped by the natural environment and not by the clinical habitat.