Search

The skin microbiota of atopic dermatitis (AD) patients is characterized by increased Staphylococcus aureus colonization, which exacerbates disease symptoms and has been linked to reduced bacterial diversity. Skin bacterial communities in AD patients have mostly been described at family and genus levels, while species-level characterization has been limited. In this study, we investigated the role of the bacteria belonging to the Staphylococcus genus using targeted sequencing of the tuf gene with genus-specific primers. We compared staphylococcal communities on lesional and non-lesional skin of AD patients, as well as AD patients with healthy controls, and determined the absolute abundance of bacteria present at each site. We observed that the staphylococcal community, bacterial alpha diversity, and bacterial densities were similar on lesional and non-lesional skin, whereas AD severity was associated with significant changes in staphylococcal composition. Increased S. aureus, Staphylococcus capitis, and Staphylococcus lugdunensis abundances were correlated with increased severity. Conversely, Staphylococcus hominis abundance was negatively correlated with severity. Furthermore, S. hominis relative abundance was reduced on AD skin compared to healthy skin. In conclusion, various staphylococcal species appear to be important for skin health.

Background: The pathophysiology in atopic dermatitis (AD) is not fully understood, but immune dysfunction, skin barrier defects, and alterations of the skin microbiota are thought to play important roles. AD skin is frequently colonized with Staphylococcus aureus (S. aureus) and microbial diversity on lesional skin (LS) is reduced compared to on healthy skin. Treatment with narrow-band ultraviolet B (nb-UVB) leads to clinical improvement of the eczema and reduced abundance of S. aureus. However, in-depth knowledge of the temporal dynamics of the skin microbiota in AD in response to nb-UVB treatment is lacking and could provide important clues to decipher whether the microbial changes are primary drivers of the disease, or secondary to the inflammatory process. Objectives: To map the temporal shifts in the microbiota of the skin, nose, and throat in adult AD patients after nb-UVB treatment. Methods: Skin swabs were taken from lesional AD skin (n = 16) before and after 3 treatments of nb-UVB, and after 6–8 weeks of full-body treatment. We also obtained samples from non-lesional skin (NLS) and from the nose and throat. All samples were characterized by 16S rRNA gene sequencing. Results: We observed shifts towards higher diversity in the microbiota of lesional AD skin after 6–8 weeks of treatment, while the microbiota of NLS and of the nose/throat remained unchanged. After only 3 treatments with nb-UVB, there were no significant changes in the microbiota. Conclusion: Nb-UVB induces changes in the skin microbiota towards higher diversity, but the microbiota of the nose and throat are not altered.

Since the late 1990s, changes in the epidemiology of methicillin-resistant Staphylococcus aureus (MRSA) were recognized with the emergence of community-associated MRSA (CA-MRSA). CA-MRSA belonging to clonal complex 152 (CC152), carrying the small staphylococcal cassette chromosome mec (SCCmec) type V and encoding the Panton-Valentine leukocidin (PVL), has been observed in Europe. The aim of this study was to investigate its origin, evolution, and dissemination. Whole-genome sequencing was performed on a global collection of 149 CC152 isolates spanning 20 years (93 methicillin-susceptible S. aureus [MSSA] and 56 MRSA isolates). Core genome phylogeny, Bayesian inference, in silico resistance analyses, and genomic characterization were applied. Phylogenetic analysis revealed two major distinct clades, one dominated by MSSA and the other populated only by MRSA. The MSSA isolates were predominately from sub-Saharan Africa, whereas MRSA was almost exclusively from Europe. The European MRSA isolates all harbored an SCCmec type V (5C2&5) element, whereas other SCCmec elements were sporadically detected in MRSA from the otherwise MSSA-dominated clade, including SCCmec types IV (2B), V (5C2), and XIII (9A). In total, 93% of the studied CC152 isolates were PVL positive. Bayesian coalescent inference suggests an emergence of the European CC152-MRSA in the 1990s, while the CC152 lineage dates back to the 1970s. The CA-MRSA CC152 clone mimics the European CC80 CA-MRSA lineage by its emergence from a PVL-positive MSSA ancestor from North Africa or Europe. The CC152 lineage has acquired SCCmec several times, but acquisition of SCCmec type V (5C2&5) seems associated with expansion of MRSA CC152 in Europe.IMPORTANCE Understanding the evolution of CA-MRSA is important in light of the increasing importance of this reservoir in the dissemination of MRSA. Here, we highlight the story of the CA-MRSA CC152 lineage using whole-genome sequencing on an international collection of CC152. We show that the evolution of this lineage is novel and that antibiotic usage may have the potential to select for the phage-encoded Panton-Valentine leukocidin. The diversity of the strains correlated highly to geography, with higher level of resistance observed among the European MRSA isolates. The mobility of the SCCmec element is mandatory for the emergence of novel MRSA lineages, and we show here distinct acquisitions, one of which is linked to the successful clone found throughout Europe today.

Staphylococcus aureus is a colonizing opportunistic pathogen and a leading cause of bloodstream infection with high morbidity and mortality. S. aureus carriage frequency is reportedly between 20 and 40 % among healthy adults, with S. aureus colonization considered to be a risk factor for S. aureus bacteraemia. It is unknown whether a genetic component of the bacterium is associated with S. aureus bacteraemia in comparison to nasal carriage strains. Previous association studies primarily focusing on the clinical outcome of an S. aureus infection have produced conflicting results, often limited by study design challenged by sample collections and the clonal diversity of S. aureus . To date, no study has investigated whether genomic features separate nasal carriage isolates from S. aureus bacteraemia isolates within a single clonal lineage. Here we have investigated whether genomic features, including single-nucleotide polymorphisms (SNPs), genes, or kmers, distinguish S. aureus nasal carriage isolates from bacteraemia isolates that all belong to the same clonal lineage [clonal complex 45 (CC45)] using whole-genome sequencing (WGS) and a genome-wide association (GWA) approach. From CC45, 100 isolates (50 bacteraemia and 50 nasal carriage, geographically and temporally matched) from Denmark were whole-genome sequenced and subjected to GWA analyses involving gene copy number variation, SNPs, gene content, kmers and gene combinations, while correcting for lineage effects. No statistically significant association involving SNPs, specific genes, gene variants, gene copy number variation, or a combination of genes was identified that could distinguish bacteraemia isolates from nasal carriage isolates. The presented results suggest that all S. aureus nasal CC45 isolates carry the potential to cause invasive disease, as no core or accessory genome content or variations were statistically associated with invasiveness.

E. coli ST131 is an important extraintestinal pathogenic lineage. A signature characteristic of ST131 is its ability to asymptomatically colonize the gastrointestinal tract and then opportunistically cause extraintestinal infections, such as cystitis, pyelonephritis, and urosepsis. In this study, we identified an ST131 H27 sublineage that has acquired the enteroaggregative diarrheagenic phenotype, spread across multiple continents, and caused multiple outbreaks of community-acquired ESBL-associated bloodstream infections in Denmark. The strain’s ability to both cause diarrhea and innocuously colonize the human gastrointestinal tract may facilitate its dissemination and establishment in the community., Escherichia coli sequence type 131 (ST131) is a major cause of urinary and bloodstream infections. Its association with extended-spectrum β-lactamases (ESBLs) significantly complicates treatment. Its best-described component is the rapidly expanding H30Rx clade, containing allele 30 of the type 1 fimbrial adhesin gene fimH. This lineage appears to have emerged in the United States and spread around the world in part due to the acquisition of the ESBL-encoding blaCTX-M-15 gene and resistance to fluoroquinolones. However, non-H30 ST131 sublineages with other acquired CTX-M-type resistance genes are also emerging. Based on whole-genome analyses, we describe here the presence of an (fimH) H27 E. coli ST131 sublineage that has recently caused an outbreak of community-acquired bacteremia and recurrent urinary tract infections (UTIs) in Denmark. This sublineage has acquired both a virulence plasmid (pAA) that defines the enteroaggregative E. coli (EAEC) diarrheagenic pathotype and multiple genes associated with extraintestinal E. coli (ExPEC); combined, these traits have made this particular ST131 sublineage successful at colonizing its human host and causing recurrent UTI. Moreover, using a historic World Health Organization (WHO) E. coli collection and publicly available genome sequences, we identified a global H27 EAEC ST131 sublineage that dates back as far as 1998. Most H27 EAEC ST131 isolates harbor pAA or pAA-like plasmids, and our analysis strongly implies a single ancestral acquisition among these isolates. These findings illustrate both the profound plasticity of this important pathogenic E. coli ST131 H27 sublineage and genetic acquisitions of EAEC-specific virulence traits that likely confer an enhanced ability to cause intestinal colonization.

The human vagina harbor a rich microbiota. The optimal state is dominated by lactobacilli that help to maintain health and prevent various diseases. However, the microbiota may rapidly change to a polymicrobial state that has been linked to a number of diseases. In the present study, the temporal changes of the vaginal microbiota in patients treated for sexually transmitted diseases or bacterial vaginosis (BV) and in untreated controls were studied for 26 days. The patients included 52 women treated with azithromycin, tetracyclines or moxifloxacin for present or suspected infection withChlamydia trachomatisorMycoplasma genitalium. Women with concurrent BV were also treated with metronidazole. The controls were 10 healthy women of matching age. The microbiota was analyzed by 16S rRNA gene deep sequencing, specific qPCRs and microscopy. There was generally good correlation between Nugent score and community state type (CST) and qPCR confirmed the sequencing results. By sequencing, more than 600 different taxa were found, but only 33 constituted more than 1 parts per thousand of the sequences. In both patients and controls the microbiota could be divided into three different community state types, CST-I, CST-III and CST-IV. Without metronidazole, the microbiota remained relatively stable regarding CST although changes were seen during menstrual periods. Administration of metronidazole changed the microbiota from CST-IV to CST-III in approximately 50% of the treated patients. In contrast, the CST was generally unaffected by azithromycin or tetracyclines. In 30% of the BV patients,Gardnerella vaginaliswas not eradicated by metronidazole. The majority of women colonized withUreaplasma parvumremained positive after azithromycin whileU.urealyticumwas eradicated.

Staphylococcus aureus is a colonizing opportunistic pathogen and a leading cause of bloodstream infection with high morbidity and mortality. S. aureus carriage frequency is reportedly between 20 and 40 % among healthy adults, with S. aureus colonization considered to be a risk factor for S. aureus bacteraemia. It is unknown whether a genetic component of the bacterium is associated with S. aureus bacteraemia in comparison to nasal carriage strains. Previous association studies primarily focusing on the clinical outcome of an S. aureus infection have produced conflicting results, often limited by study design challenged by sample collections and the clonal diversity of S. aureus . To date, no study has investigated whether genomic features separate nasal carriage isolates from S. aureus bacteraemia isolates within a single clonal lineage. Here we have investigated whether genomic features, including single-nucleotide polymorphisms (SNPs), genes, or kmers, distinguish S. aureus nasal carriage isolates from bacteraemia isolates that all belong to the same clonal lineage [clonal complex 45 (CC45)] using whole-genome sequencing (WGS) and a genome-wide association (GWA) approach. From CC45, 100 isolates (50 bacteraemia and 50 nasal carriage, geographically and temporally matched) from Denmark were whole-genome sequenced and subjected to GWA analyses involving gene copy number variation, SNPs, gene content, kmers and gene combinations, while correcting for lineage effects. No statistically significant association involving SNPs, specific genes, gene variants, gene copy number variation, or a combination of genes was identified that could distinguish bacteraemia isolates from nasal carriage isolates. The presented results suggest that all S. aureus nasal CC45 isolates carry the potential to cause invasive disease, as no core or accessory genome content or variations were statistically associated with invasiveness.

A retrospective study of Staphylococcus aureus isolates from orthopaedic patients treated between 2000 and 2017 at Akershus University Hospital, Norway was performed using a genome-wide association approach. The aim was to characterize and investigate molecular characteristics unique to S. aureus isolates from HHA associated prosthetic joint infections and potentially explain the HHA patients’ elevated 1-year mortality compared to a non-HHA group. The comparison group consisted of patients with non-HHA lower-extremity implant-related S. aureus infections. S. aureus isolates from diagnostic patient samples were whole-genome sequenced. Univariate and multivariate analyses were performed to detect group-associated genetic signatures. A total of 62 HHA patients and 73 non-HHA patients were included. Median age (81 years vs. 74 years; p p S. aureus can cause invasive infections given the opportunity.

Since its emergence in the early 2000s, livestock-associated methicillin-resistant Staphylococcus aureus clonal complex 398 (LA-MRSA CC398) has led to an increasing number of human infections in Denmark and other European countries with industrial pig production. LA-MRSA CC398 is primarily associated with skin infections among pig farm workers but is also increasingly recognized as a cause of life-threatening disease among elderly and immunocompromised people. Pig farm workers may serve as vehicles for the spread of LA-MRSA CC398 and other farm-origin bacteria between farms and into the general population. Yet, little is known about the bacterial community dynamics in pig farm workers and other persons with long- and short-term exposure to the pig farm environment. To gain insight into this, we investigated the nasal microbiomes in pig farm workers during a workweek on four LA-MRSA CC398-positive pig farms, as well as in short-term visitors two hours before, immediately after, and 48 hours after a 1-hour visit to another LA-MRSA CC398-positive pig farm. S. aureus and LA-MRSA CC398 carriage was quantified by means of culture, and the composition of the bacterial communities was investigated through sequencing of the 16S rRNA gene. Pig farm workers often carried LA-MRSA CC398 and other bacteria from the pig farm environment, both at work and at home, although at lower levels at home. In contrast, short-term visitors were subject to a less dramatic and rapidly reversible change in the nasal bacterial community composition. These results suggest that pig farm workers may be an important source of LA-MRSA CC398 and perhaps other pathogens of human and veterinary relevance.

Today, next-generation whole-genome sequencing (WGS) is increasingly used to determine the genetic relationships of bacteria on a nearly whole-genome level for infection control purposes and molecular surveillance. Here, we conducted a multicenter ring trial comprising five laboratories to determine the reproducibility and accuracy of WGS-based typing. The participating laboratories sequenced 20 blind-coded Staphylococcus aureus DNA samples using 250-bp paired-end chemistry for library preparation in a single sequencing run on an Illumina MiSeq sequencer. The run acceptance criteria were sequencing outputs >5.6 Gb and Q30 read quality scores of >75%. Subsequently, spa typing, multilocus sequence typing (MLST), ribosomal MLST, and core genome MLST (cgMLST) were performed by the participants. Moreover, discrepancies in cgMLST target sequences in comparisons with the included and also published sequence of the quality control strain ATCC 25923 were resolved using Sanger sequencing. All five laboratories fulfilled the run acceptance criteria in a single sequencing run without any repetition. Of the 400 total possible typing results, 394 of the reported spa types, sequence types (STs), ribosomal STs (rSTs), and cgMLST cluster types were correct and identical among all laboratories; only six typing results were missing. An analysis of cgMLST allelic profiles corroborated this high reproducibility; only 3 of 183,927 (0.0016%) cgMLST allele calls were wrong. Sanger sequencing confirmed all 12 discrepancies of the ring trial results in comparison with the published sequence of ATCC 25923. In summary, this ring trial demonstrated the high reproducibility and accuracy of current next-generation sequencing-based bacterial typing for molecular surveillance when done with nearly completely locked-down methods.

More than 55 million mink skins were produced globally in 2017. As a consequence, a large number of people are employed in mink production worldwide. In Denmark, farmed mink were found to constitute a reservoir of methicillin-resistant Staphylococcus aureus (MRSA) clonal complex (CC) 398 and 6000 mink farm workers in Denmark are potentially exposed to LA-MRSA CC398. The study aim was to elucidate the source of LA-MRSA CC398 in mink farms and to investigate possible transmission to humans. In total, 161 LA-MRSA CC398 isolates from mink (n = 65), mink feed (n = 16) and humans (n = 80) with reported contact to mink, were whole-genome sequenced and compared to 183 LA-MRSA CC398 isolates from Danish pigs and an international collection of 89 S. aureus CC398 isolates. Most of the mink-associated isolates clustered within the predominant LA-MRSA CC398 lineages circulating in the Danish pig production, supporting that pigs are a source of LA-MRSA CC398 in mink feed, mink, and mink farmers.

APEC causes a range of infections in poultry, collectively called colibacillosis, and is the leading cause of mortality and is associated with major economic significance in the poultry industry. A growing number of studies have suggested APEC as an external reservoir of human ExPEC, including UPEC, which is a reservoir. ExPEC belonging to ST95 is considered one of the most important pathogens in both poultry and humans. This study is the first in-depth whole-genome-based comparison of ST95 E. coli which investigates both the core genomes as well as the accessory genomes of avian and human ExPEC. We demonstrated that multiple lineages of ExPEC belonging to ST95 exist, of which the majority may cause infection in humans, while only part of the ST95 cluster seem to be avian pathogenic. These findings further support the idea that urinary tract infections may be a zoonotic infection., Avian-pathogenic Escherichia coli (APEC) is a subgroup of extraintestinal pathogenic E. coli (ExPEC) presumed to be zoonotic and to represent an external reservoir for extraintestinal infections in humans, including uropathogenic E. coli (UPEC) causing urinary tract infections. Comparative genomics has previously been applied to investigate whether APEC and human ExPEC are distinct entities. Even so, whole-genome-based studies are limited, and large-scale comparisons focused on single sequence types (STs) are not available yet. In this study, comparative genomic analysis was performed on 323 APEC and human ExPEC genomes belonging to sequence type 95 (ST95) to investigate whether APEC and human ExPEC are distinct entities. Our study showed that APEC of ST95 did not constitute a unique ExPEC branch and was genetically diverse. A large genetic overlap between APEC and certain human ExPEC was observed, with APEC located on multiple branches together with closely related human ExPEC, including nearly identical APEC and human ExPEC. These results illustrate that certain ExPEC clones may indeed have the potential to cause infection in both poultry and humans. Previously described ExPEC-associated genes were found to be encoded on ColV plasmids. These virulence-associated plasmids seem to be crucial for ExPEC strains to cause avian colibacillosis and are strongly associated with strains of the mixed APEC/human ExPEC clusters. The phylogenetic analysis revealed two distinct branches consisting of exclusively closely related human ExPEC which did not carry the virulence-associated plasmids, emphasizing a lower avian virulence potential of human ExPEC in relation to an avian host. IMPORTANCE APEC causes a range of infections in poultry, collectively called colibacillosis, and is the leading cause of mortality and is associated with major economic significance in the poultry industry. A growing number of studies have suggested APEC as an external reservoir of human ExPEC, including UPEC, which is a reservoir. ExPEC belonging to ST95 is considered one of the most important pathogens in both poultry and humans. This study is the first in-depth whole-genome-based comparison of ST95 E. coli which investigates both the core genomes as well as the accessory genomes of avian and human ExPEC. We demonstrated that multiple lineages of ExPEC belonging to ST95 exist, of which the majority may cause infection in humans, while only part of the ST95 cluster seem to be avian pathogenic. These findings further support the idea that urinary tract infections may be a zoonotic infection.

BACKGROUND: MicroRNAs are small noncoding RNAs involved in the post-transcriptional regulation of protein synthesis. Extracellular microRNAs are accessible in a stable form in biofluids.OBJECTIVES: The aim was to identify individual microRNAs and/or subsets of microRNAs in CSF with biomarker potential and thus identify specific putative pathophysiological pathways.METHODS: In a two-step exploratory study design of PD, MSA, PSP, and controls, we initially profiled CSF microRNAs in a pilot cohort (n = 40) by screening for 372 microRNAs. Subsequently, we attempted to validate findings in an independent study cohort in CSF (n = 118) and ethylenediaminetetraacetic acid plasma (n = 114). This study cohort encompassed 46 microRNAs, of which 26 were singled out from the pilot cohort, and an additional 20 microRNAs were added based on previous publications. The most accurate diagnostic microRNA classifiers were identified in a multivariable logistic regression model adjusted for age and sex.RESULTS: A set of three microRNAs in CSF discriminated PD and MSA from controls with good diagnostic accuracy by receiver operating characteristics curve evaluation. The microRNAs were for PD versus controls: miR-7-5p, miR-331-5p, and miR-145-5p (area under the curve = 0.88) and MSA versus controls: miR-7-5p, miR-34c-3p, and miR-let-7b-5p (area under the curve = 0.87). The classifier that best distinguished MSA and PD consisted of two microRNAs: miR-9-3p and miR-106b-5p (area under the curve = 0.73). A single microRNA, miR-106b-5p, provided the best discrimination between PD and PSP (area under the curve = 0.85) in the CSF.CONCLUSIONS: Levels of specific trios of CSF-microRNAs discriminate well between α-synucleinopathies (PD and MSA) and controls. The results need to be validated in larger, independent cohorts. © 2018 International Parkinson and Movement Disorder Society.

Recurrent urinary tract infection (rUTI) remains a major problem for many women and therefore the pursuit for genomic and phenotypic traits which could define rUTI has been ongoing. The present study applied a genomic approach to investigate recurrent urinary tract infections by comparative analyses of recurrent and non-recurrent Escherichia coli isolates from general practice. From whole-genome sequencing data, phylogenetic clustering and genomic traits were studied on a collection of isolates which caused recurrent infection compared to non-recurrent isolates. In addition, genomic variation between the 1st and following infection was studied on a subset of the isolates. Evidence of limited adaptation between the recurrent infections based on single nucleotide polymorphism analyses with a range of 0–13 non-synonymous single nucleotide polymorphisms (SNPs) between the paired isolates. This included an overrepresentation of SNPs in metabolism genes. We identified several genes which were more common in rUTI isolates, including nine fimbrial genes, however, not significantly after false-discovery rate. Finally, the results show that recurrent isolates of the present dataset are not distinctive by variation in the core genome, and thus, did not cluster distinct from non-rUTI isolates in a SNP phylogeny.

OBJECTIVES We unexpectedly identified MRSA isolates carrying mecC (mecC-MRSA) from a Danish swine farm located in eastern Zealand. The objective of the present study was to investigate the origin of these isolates and their genetic relatedness to other mecC-MRSA isolates from Zealand. METHODS WGS was used to infer the phylogenetic relationship between 19 identified mecC-MRSA isolates from the swine farm and 34 additional epidemiologically unrelated human isolates from the same geographical region of Denmark. Variations in the accessory genome were investigated by bioinformatics tools, and antibiotic susceptibility profiles were assessed by MIC determination. RESULTS mecC-MRSA was isolated from a domestic swine farm, but not from cattle reared at the same farm. Phylogenetic analysis revealed that all mecC-MRSA isolates from both farm animals and workers formed a separate cluster, whereas human isolates from the same municipality belonged to a closely related cluster. Analysis of the accessory genome supported this relationship. CONCLUSIONS To the best of our knowledge, this is the first report of mecC-MRSA isolated from domestic swine. The investigation strongly indicates that transmission of mecC-MRSA has taken place on the swine farm between the farmers and swine. The close clustering of farm isolates and isolates from the same municipality suggests a local transmission of mecC-MRSA.

Genes & Development 28: 2498–2517 (2014) While revisiting next-generation sequencing (NGS) data analyses from the above-mentioned article for another project, we discovered four unintentional errors. After correction of the errors, we reapplied the reported methods and found the originally stated conclusions to be unchanged. However, several figure panels and reported numbers were altered. Corrected Figures 1​1​​​​–7 are shown below, and corrected Supplemental Figures S2, S3, S5, S6, and S7 and Supplemental Tables S2, S3, and S4 are accessible on the journal site online. Figure 1. Figure 2. Figure 3. Figure 4. Figure 5. Figure 6. Figure 7. The following errors were found and are corrected here (a flow chart illustrating the errors and how figures/tables were affected is also available in the Revised Supplemental Material): (1) Due to a programming error during the processing of “5′ end-seq” data, reads spanning exon–exon junctions were erroneously omitted. This effectively led to a subsampling of the data, leaving out roughly a quarter of all “5′ end-seq” reads. In this Corrigendum, we now present analyses based on all of the reads. The increased amount of data influences the calculation of nonsense-mediated RNA decay (NMD) thresholds as well as most of the downstream computational analyses (see Box 1 in the flow chart). This affected Figure 1, E and G; Figure 2, C and D; Figure 3, A and B; Figure 4, A and B; Figure 5, A, B, and D; Figure 6B; Supplemental Figure S2, A, C, E, and H; Supplemental Figure S3, A, D, and G; Supplemental Figure S5, A–H; Supplemental Figure S7D; Supplemental Table S2; Supplemental Table S3; and Supplemental Table S4. (2) The plots displaying NGS data mapping to specific genes or splice junctions were based on an alternative mapping procedure (RUM) (Grant et al. 2011) rather than the one reported in the Materials and Methods section (TopHat version 2). The two procedures gave very similar mapping results. In the corrected Figures 1, E–G, and ​and77 (below) as well as revised Supplemental Figures 2, A–E, and 3, A, B, D, E, and G (Revised Supplemental Material), we now present plots based on the reported mapping procedure (see Box 2 in the flow chart). (3) In the top panel of Figure 6C and middle panel of Supplemental Figure S7D, we compared the number of transcript isoforms produced from snoRNA host genes with similarly expressed normal protein-coding genes. Although we reported the use of de novo transcripts based on our own RNA sequencing (RNA-seq) data from HEK293 cells, in fact, Gencode version 17 annotated transcripts were used. In the corrected Figure 6 (below) and Supplemental Figure S7 (Revised Supplemental Material), we now present analyses based on de novo assembled transcripts as stated in the corresponding legends (see Box 3 in the flow chart). We sincerely apologize for any inconvenience caused by this negligence and stress once more that all original conclusions are left intact after the conducted modifications of the article. The text changes reflecting the corrected data are noted below: p. 2499, column 2, paragraph 3: By combining global identification of nonsense RNAs and their corresponding decay intermediates, we identified primary NMD-responsive isoforms from up to 15% of all expressed genes. p. 2500, column 2, paragraph 3: To investigate whether this observation was supported by our global data, we characterized NMD-specific endocleavage events in a reference set of annotated NMD substrates, consisting of 801 transcripts with a total of 1544 potential endocleavage sites (“NMD reference set” derived from annotation) (see Materials and Methods; Supplemental Table S2; Harrow et al. 2012), by mapping their positions relative to the annotated termination codons. p. 2503, legend to Figure 2D: (D) Box plots illustrating the distribution of RNA-seq-based expression levels of the transcripts corresponding to the 5′ end-seq signals plotted in C. (***) P-value ≤ 0.001. p. 2505, column 2, paragraph 1: Finally, the NMD-sensitive transcripts identified independently by 5′ end-seq and RNA-seq (Supplemental Table S3) were compared, which yielded 3223 NMD-sensitive transcripts arising from 1563 assembled genes (Fig. 4B, right panel; Supplemental Table S4). p. 2505, column 2, paragraph 1: When applying a less stringent (“relaxed”) cutoff for the initial identification of peaks in the 5′ end-seq data (noncorrected P-value ≤ 0.005 based on negative binomial fitting), we could detect 9060 transcripts corresponding to 4245 genes in the overlap with the RNA-seq data (Supplemental Fig. S5D; Supplemental Table S4). p. 2505, column 2, paragraph 1: Using the present data, we confirmed that 25 of 53 SRSF protein-coding genes (Supplemental Table S5; total number of SRSF genes based on Long and Caceres 2009) could be detected in our most stringent NMD substrate set based on the combined analysis of RNA-seq and 5′ end-seq. p. 2505, column 2, paragraph 1: Even more genes could be included by using the “relaxed” criteria (40 out of 53) or taking the maximum union of independent RNA-seq and 5′ end-seq approaches (51 out of 53) (Supplemental Fig. S5G). p. 2506, column 1, paragraph 3: We found that 59 (34%) and 104 (60%) of these genes were responsive to NMD based on the “stringent” and “relaxed” NMD substrate set, respectively (Supplemental Table S4). p. 2506, column 1, paragraph 3: Even more snoRNA hosts are potential targets of NMD, since the independent 5′ end-seq and RNA-seq procedures combined include up to 166 (96%) of these genes (Supplemental Fig. S5H; Supplemental Table S3). p. 2507, legend to Figure 5A: Three asterisks indicate that snoRNA host genes are significantly enriched compared with protein-coding genes (Fisher two-sided test, P = 1.5 × 10−13). p. 2508, column 2, paragraph 1: RNAs containing the latter splice junction constitute ∼60% of the total amount of exon 2-containing RNA under control conditions (Fig. 7, histograms next to the schematics of the four isoforms show the percent of the total amount of exon 2-containing RNA within samples). p. 2509, column 1, paragraph 1: However, upon inhibition of NMD, this variant, although being somewhat NMD-sensitive, only constitutes ∼25% of the total RNA. p. 2509, column 1, paragraph 1: Conversely, the isoforms encoding only SNORD49A comprise ∼60% of the total RNA when NMD is inhibited as opposed to ∼30% under normal conditions (Fig. 7, i,ii). p. 2509, column 1, paragraph 1: Similarly, the relative level of the transcript variants giving rise to both SNORD49A and SNORD49B increase from ∼5% to ∼10% upon inhibition of NMD (Fig. 7, iii). p. 2509, column 2, paragraph 3: This strategy revealed nonsense RNA isoforms produced from 6% and 15% of expressed gene loci, respectively. p. 2513, column 1, paragraph 3: This gave us a total set of 1321 transcripts corresponding to 916 genes. p. 2513, column 1, paragraph 3: This gave us a total set of 256 transcripts derived from 134 genes. p. 2513, column 2, paragraph 4: SMG6/XRN1 > 1.56 × CTRL and UPF1/XRN1 > 1.56 × CTRL; SMG6/XRN1 > 1.23 × XRN1, and UPF1/XRN1 > 1.23 × XRN1. p. 2514, column 1, paragraph 6: NMD-specific endocleavage site: XRN1 > 1.36 × SMG6/XRN1 and XRN1 > 1.36 × UPF1/XRN1; NMD-specific decapping site: SMG6/XRN1 > 1.93 × XRN1 and UPF1/XRN1 > 1.93 × XRN1. p. 2516, column 1 The reference below was also mistakenly omitted from the published article. Grant GR, Farkas MH, Pizarro AD, Lahens NF, Schug J, Brunk BP, Stoeckert CJ, Hogenesch JB, Pierce EA. 2011. Comparative analysis of RNA-seq alignment algorithms and the RNA-seq unified mapper (RUM). Bioinformatics 27: 2518–2528. Supplemental Material, p. 12, legend to Supplemental Figure S5E: (E) As Fig. 5A, but based on “relaxed” criteria for identification of NMD-responsive genes. *** Indicates that snoRNA host genes are significantly enriched compared to protein-coding genes (Fisher two-sided test, P = 8.9 × 10−12). doi: 10.1101/gad.281881.116

Escherichia colisequence type 131 (ST131) is a major cause of urinary and bloodstream infections and its association with extended-spectrum β-lactamases (ESBL) significantly complicates treatment. Most notorious is its rapidly expandingH30-Rx clade (named for containing allele 30 of the type-1 fimbrial adhesin genefimHand extensive antimicrobial resistance), which appears to have emerged in the United States due in part due to the acquisition of the ESBL-encodingblaCTX-M-15gene and resistance to fluoroquinolones. However, non-H30 ST131 lineages with acquired CTX-M-type resistance genes also are emerging. Based on whole-genome analyses, we describe here the presence of an (fimH)H27E. coliST131 lineage that currently is causing an outbreak of community-acquired bacteremia and recurrent urinary tract infections (UTIs) in Denmark. This lineage has acquired both a virulence plasmid (pAA) that defines the enteroaggregativeE. coli(EAEC) diarrheagenic pathotype and multiple genes associated with extraintestinalE. coli(ExPEC) that combined has made this particular ST131 lineage highly successful at colonizing its human host and cause recurrent UTI. Moreover, using a historic World Health OrganizationE. colicollection and publically available genome sequences, we identify a globalH27 EAEC ST131 lineage dating back as far as 1998. MostH27 EAEC ST131 isolates harbor pAA or pAA-like plasmids, which analysis strongly imply was caused by a single ancestral acquisition. These findings illustrate the profound plasticity of this important pathogenicE. coli H27 lineage in general, and the genetic acquisitions of EAEC-specific virulence traits that likely confer an enhanced ability to cause intestinal colonization.ImportanceTheE. coliST131 lineage is a notorious extraintestinal pathogen. A signature characteristic of ST131 is its ability to asymptomatically colonize the gastrointestinal tract and then opportunistically cause extraintestinal infections, such as cystitis, pyelonephritis and urosepsis. In this study, we report a novel ST131 sublineage that has acquired the enteroaggregative diarrheagenic phenotype, spread across multiple continents and has been associated with multiple outbreaks of community-acquired bloodstream infections in Denmark. The strain’s ability to both cause diarrhea and colonize the human gastrointestinal tract may facilitate its dissemination and establishment in the community, whereas the strain’s clonal nature may facilitate targeted control strategies, such as vaccination.

Inflammatory bowel disease (IBD) is a chronic intestinal disorder, with two main types: Crohn’s disease (CD) and ulcerative colitis (UC), whose molecular pathology is not well understood. The majority of IBD-associated SNPs are located in non-coding regions and are hard to characterize since regulatory regions in IBD are not known. Here we profile transcription start sites (TSSs) and enhancers in the descending colon of 94 IBD patients and controls. IBD-upregulated promoters and enhancers are highly enriched for IBD-associated SNPs and are bound by the same transcription factors. IBD-specific TSSs are associated to genes with roles in both inflammatory cascades and gut epithelia while TSSs distinguishing UC and CD are associated to gut epithelia functions. We find that as few as 35 TSSs can distinguish active CD, UC, and controls with 85% accuracy in an independent cohort. Our data constitute a foundation for understanding the molecular pathology, gene regulation, and genetics of IBD., Many SNPs associated with inflammatory bowel disease are located in non-coding genomic regions. Here, the authors perform CAGE-sequencing on descending colon biopsies of Crohn’s disease and ulcerative colitis patients to map transcription start sites and enhancer activity for analysis of regulatory regions.

Uncaging promoter and enhancer dynamics In order to understand cellular differentiation, it is important to understand the timing of the regulation of gene expression. Arner et al. used cap analysis of gene expression (CAGE) to analyze gene enhancer and promoter activities in a number of human and mouse cell types. The RNA of enhancers was transcribed first, followed by that of transcription factors, and finally by genes that are not transcription factors. Science , this issue p. 1010

Most Staphylococcus aureus isolates can cause invasive disease given the right circumstances, but it is unknown if some isolates are more likely to cause severe infections than others. S. aureus bloodstream isolates from 120 patients with definite infective endocarditis and 121 with S. aureus bacteraemia without infective endocarditis underwent whole-genome sequencing. Genome-wide association analysis was performed using a variety of bioinformatics approaches including SNP analysis, accessory genome analysis and k-mer based analysis. Core and accessory genome analyses found no association with either of the two clinical groups. In this study, the genome sequences of S. aureus bloodstream isolates did not discriminate between bacteraemia and infective endocarditis. Based on our study and the current literature, it is not convincing that a specific S. aureus genotype is clearly associated to infective endocarditis in patients with S. aureus bacteraemia.

Objectives: To describe the changing epidemiology of vancomycin-resistant Enterococcus faecium and Enterococcus faecalis in clinical samples in Denmark 2005-15 according to species and van type, and, furthermore, to investigate the genetic relatedness of the clinical E. faecium isolates from 2015.Methods: During 2005-14, all clinical VRE isolates were tested for the presence of vanA/B/C genes by PCR. In 2015, all clinical VRE isolates were whole-genome sequenced. From the WGS data, the presence of van genes and MLST STs were extracted in silico . Core-genome MLST (cgMLST) analysis was performed for the vancomycin-resistant E. faecium isolates.Results: During 2005-15, 1043 vanA E. faecium , 25 vanB E. faecium , 4 vanA E. faecalis and 28 vanB E. faecalis were detected. The number of VRE was 200 isolates/year in 2013-15. In 2015, 368 vanA E. faecium and 1 vanB E. faecium were detected along with 1 vanA E. faecalis and 1 vanB E. faecalis . cgMLST subdivided the 368 vanA E. faecium isolates into 33 cluster types (CTs), whereas the vanB E. faecium isolate belonged to a different CT. ST203-CT859 was most prevalent (51%), followed by ST80-CT14 (22%), ST117-CT24 (6%), ST80-CT866 (4%) and ST80-CT860 (2%). Comparison with the cgMLST.org database, previous studies and personal communications with neighbouring countries revealed that the novel cluster ST203-CT859 emerged in December 2014 and spread to the south of Sweden and the Faroe Islands during 2015.Conclusions: VRE increased in Denmark during 2005-15 due to the emergence of several vanA E. faecium clones.

The faecal flora is a common reservoir for urinary tract infection (UTI), and Escherichia coli (E. coli) is frequently found in this reservoir without causing extraintestinal infection. We investigated these E. coli reservoirs by whole-genome sequencing a large collection of E. coli from healthy controls (faecal), who had never previously had UTI, and from UTI patients (faecal and urinary) sampled from the same geographical area. We compared MLST types, phylogenetic relationship, accessory genome content and FimH type between patient and control faecal isolates as well as between UTI and faecal-only isolates, respectively. Comparison of the accessory genome of UTI isolates to faecal isolates revealed 35 gene families which were significantly more prevalent in the UTI isolates compared to the faecal isolates, although none of these were unique to one of the two groups. Of these 35, 22 belonged to a genomic island and three putatively belonged to a type VI secretion system (T6SS). MLST types and SNP phylogeny indicated no clustering of the UTI or faecal E. coli from patients distinct from the control faecal isolates, although there was an overrepresentation of UTI isolates belonging to clonal lineages CC73 and CC12. One combination of mutations in FimH, N70S/S78N, was significantly associated to UTI, while phylogenetic analysis of FimH and fimH identified no signs of distinct adaptation of UTI isolates compared to faecal-only isolates not causing UTI. In summary, the results showed that (i) healthy women who had never previously had UTI carried faecal E. coli which were overall closely related to UTI and faecal isolates from UTI patients; (ii) UTI isolates do not cluster separately from faecal-only isolates based on SNP analysis; and (iii) 22 gene families of a genomic island, putative T6SS proteins as well as specific metabolism and virulence associated proteins were significantly more common in UTI isolates compared to faecal-only isolates and (iv) evolution of fimH for these isolates was not linked to the clinical source of the isolates, apart from the mutation combination N70S/S78N, which was correlated to UTI isolates of phylogroup B2. Combined, these findings illustrate that faecal and UTI isolates, as well as faecal-only and faecal-UTI isolates, are closely related and can only be distinguished, if at all, by their accessory genome.

Background Escherichia coli infections known as colibacillosis constitute a considerable challenge to poultry farmers worldwide, in terms of decreased animal welfare and production economy. Colibacillosis is caused by avian pathogenic E. coli (APEC). APEC strains are extraintestinal pathogenic E. coli and have in general been characterized as being a genetically diverse population. In the Nordic countries, poultry farmers depend on import of Swedish broiler breeders which are part of a breeding pyramid. During 2014 to 2016, an increased occurrence of colibacillosis on Nordic broiler chicken farms was reported. The aim of this study was to investigate the genetic diversity among E. coli isolates collected on poultry farms with colibacillosis issues, using whole genome sequencing. Methods Hundred and fourteen bacterial isolates from both broilers and broiler breeders were whole genome sequenced. The majority of isolates were collected from poultry with colibacillosis on Nordic farms. Subsequently, comparative genomic analyses were carried out. This included in silico typing (sero- and multi-locus sequence typing), identification of virulence and resistance genes and phylogenetic analyses based on single nucleotide polymorphisms. Results In general, the characterized poultry isolates constituted a genetically diverse population. However, the phylogenetic analyses revealed a major clade of 47 closely related ST117 O78:H4 isolates. The isolates in this clade were collected from broiler chickens and breeders with colibacillosis in multiple Nordic countries. They clustered together with a human ST117 isolate and all carried virulence genes that previously have been associated with human uropathogenic E. coli. Conclusions The investigation revealed a lineage of ST117 O78:H4 isolates collected in different Nordic countries from diseased broilers and breeders. The data indicate that the closely related ST117 O78:H4 strains have been transferred vertically through the broiler breeding pyramid into distantly located farms across the Nordic countries. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3415-6) contains supplementary material, which is available to authorized users.

Objective Clostridium perfringens causes gastrointestinal diseases in both humans and domestic animals. Type A strains expressing the NetB toxin are the main cause of necrotic enteritis (NE) in chickens, which has remarkable impact on animal welfare and production economy in the international poultry industry. Three pathogenicity loci NELoc-1, -2 and -3 and a collagen adhesion gene cnaA have been found to be associated with NE in chickens, whereas the presence of these has not been investigated in diseased turkeys. The purpose was to investigate the virulence associated genome content and the genetic relationship among 30 C. perfringens isolates from both healthy and NE infected chickens and turkeys, applying whole-genome sequencing. Results NELoc-1, -3, netB and cnaA were significantly associated with NE isolates from chickens, whereas only NELoc-2 was commonly observed in both diseased turkeys and chickens. A putative collagen adhesion gene that encodes a von Willebrand Factor (vWF) domain was identified in all diseased turkeys and designated as cnaD. The phylogenetic analysis based on single nucleotide polymorphisms showed that the isolates generally were not closely related. These results indicate that virulence factors and pathogenicity loci associated with NE in chickens are not important to the same extent in diseased turkeys except for NELoc-2. A putative collagen adhesion gene which potentially could be of importance in regard to the NE pathogenesis in turkeys was identified and need to be further investigated. Thus, the pathogenesis of NE in turkeys appears to be different from that of broiler chickens. Electronic supplementary material The online version of this article (doi:10.1186/s13104-017-2594-9) contains supplementary material, which is available to authorized users.

Between July and November 2014, 15 community-acquired cases of Legionnaires´ disease (LD), including four with Legionella pneumophila serogroup 1 sequence type (ST) 82, were diagnosed in Northern Zealand, Denmark. An outbreak was suspected. No ST82 isolates were found in environmental samples and no external source was established. Four putative-outbreak ST82 isolates were retrospectively subjected to whole genome sequencing (WGS) followed by phylogenetic analyses with epidemiologically unrelated ST82 sequences. The four putative-outbreak ST82 sequences fell into two clades, the two clades were separated by ca 1,700 single nt polymorphisms (SNP)s when recombination regions were included but only by 12 to 21 SNPs when these were removed. A single putative-outbreak ST82 isolate sequence segregated in the first clade. The other three clustered in the second clade, where all included sequences had L. pneumophila evolution. On the other hand, strains belonging to the same ST can have only few SNP differences despite being sampled over both large timespans and geographic distances. These are two important factors to consider in outbreak investigations.

The Caco-2 cell line is one of the most important in vitro models for enterocytes, and is used to study drug absorption and disease, including inflammatory bowel disease and cancer. In order to use the model optimally, it is necessary to map its functional entities. In this study, we have generated genome-wide maps of active transcription start sites (TSSs), and active enhancers in Caco-2 cells with or without tumour necrosis factor (TNF)-α stimulation to mimic an inflammatory state. We found 520 promoters that significantly changed their usage level upon TNF-α stimulation; of these, 52% are not annotated. A subset of these has the potential to confer change in protein function due to protein domain exclusion. Moreover, we locate 890 transcribed enhancer candidates, where ∼50% are changing in usage after TNF-α stimulation. These enhancers share motif enrichments with similarly responding gene promoters. As a case example, we characterize an enhancer regulating the laminin-5 γ2-chain (LAMC2) gene by nuclear factor (NF)-κB binding. This report is the first to present comprehensive TSS and enhancer maps over Caco-2 cells, and highlights many novel inflammation-specific promoters and enhancers.

A 24-hour light and dark cycle-dependent rhythmicity pervades physiological processes in virtually all living organisms including humans. These regular oscillations are caused by external cues to endogenous, independent biological time-keeping systems (clocks). The rhythm is reflected by gene expression that varies in a circadian and specific fashion in different organs and tissues and is regulated largely by dynamic epigenetic and post-transcriptional mechanisms. This leads to well-documented oscillations of specific electrolytes, hormones, metabolites, and plasma proteins in blood samples. An emerging, important class of gene regulators is short single-stranded RNA (micro-RNA, miRNA) that interferes post-transcriptionally with gene expression and thus may play a role in the circadian variation of gene expression. MiRNAs are promising biomarkers by virtue of their disease-specific tissue expression and because of their presence as stable entities in the circulation. However, no studies have addressed the putative circadian rhythmicity of circulating, cell-free miRNAs. This question is important both for using miRNAs as biological markers and for clues to miRNA function in the regulation of circadian gene expression. Here, we investigate 92 miRNAs in plasma samples from 24 young male, healthy volunteers repeatedly sampled 9 times during a 24-hour stay in a regulated environment. We demonstrate that a third (26/79) of the measurable plasma miRNAs (using RT-qPCR on a microfluidic system) exhibit a rhythmic behavior and are distributed in two main phase patterns. Some of these miRNAs weakly target known clock genes and many have strong targets in intracellular MAPK signaling pathways. These novel findings highlight the importance of considering bio-oscillations in miRNA biomarker studies and suggest the further study of a set of specific circulating miRNAs in the regulation and functioning of biological clocks.

Importance Skin microbiome correlates with disease severity for lesional and nonlesional skin, indicating a global influence of atopic dermatitis (AD). A relation between skin microbiome and filaggrin gene (FLG) mutations proposes a possible association between skin microbiome and host genetics. Objectives To assess skin and nasal microbiome diversity and composition in patients with AD and compare with healthy controls, and to investigate the microbiome in relation to disease severity and FLG mutations in patients with AD. Design, Setting, and Participants An observational case-control study of 45 adult healthy controls and 56 adult patients with AD was carried out from January 2015 to June 2015 in a tertiary referral center, Department of Dermatology, Bispebjerg Hospital, Denmark. Exposures Bacterial swabs were taken from patients with AD (lesional skin, nonlesional skin, and anterior nares) and from healthy controls (nonlesional skin and anterior nares). Eczema severity was assessed and FLG mutations noted. Bacterial DNA was extracted from swabs, and V3-V4 16S rDNA regions amplified with PCR. Samples were analyzed at Statens Serum Institut September 2015 to September 2016. Bioinformatics analyses of the microbiome were analyzed using R statistical software (version 3.3.1, R Foundation Inc). Main Outcomes and Measures Skin microbiomes were investigated using next-generation sequencing targeting 16S ribosomal RNA. Results Microbiome alpha diversity was lower in patients with AD compared with healthy controls in nonlesional skin (effect size, 0.710; 95% CI, 0.27-1.15; P = .002), lesional skin (effect size, 0.728; 95% CI, 0.35-1.33; P = .001), and nose (effect size, 1.111; 95% CI, 0.48-0.94; P P = .02) and nonlesional skin (effect size, 0.451; 95% CI, 0.04-2.44; P = .04) in patients with AD. Microbiome composition in AD nonlesional skin was linked to FLG mutations. Conclusions and Relevance An altered microbiome composition in patients with AD in nonlesional skin, lesional skin, as well as nose, suggests a global influence of AD. Microbiome composition in AD nonlesional skin is associated with FLG mutations, proposing a possible association between the skin microbiome and host genetics.

Importance Although the pathogenesis of hidradenitis suppurativa (HS) remains enigmatic, several factors point to potential involvement of the cutaneous microbiome. Insight into the cutaneous microbiome in HS using next-generation sequencing may provide novel data on the microbiological diversity of the skin. Objective To investigate the follicular skin microbiome in patients with HS and in healthy controls. Design, Setting, and Participants This case-control study obtained punch biopsy specimens from patients with HS (lesional and nonlesional) and healthy controls between October 1, 2014, and August 1, 2016. Data were analyzed from March to November 2016. Patients with HS were recruited from the Department of Dermatology, Zealand University Hospital, Roskilde, Denmark. Biopsy specimens were analyzed at the Department of Microbiology and Infection Control, Statens Serum Institut, Copenhagen, Denmark. None of the participants received any antibiotics (systemic or topical therapy) within 1 month before the study. In patients with HS, biopsy specimens were obtained from lesional skin (axilla or groin) and nonlesional skin. Only nodules containing at least 1 visible hair follicle were biopsied. Biopsy specimens from healthy controls were obtained from the axilla only. Main Outcomes and Measures The different microbiomes were investigated using next-generation sequencing targeting 16S and 18S ribosomal RNA. Results The skin microbiome was characterized in 30 patients with HS (mean [SD] age, 46.9 [14.0] years; 19 [63% female]) and 24 healthy controls (mean [SD] age, 32.2 [12.0] years; 13 [54% female]). The next-generation sequencing data provided a previously unreported (to our knowledge) characterization of the skin microbiome in HS. The study demonstrated that the microbiome in HS differs significantly from that in healthy controls in lesional and nonlesional skin. Overall, the following 5 microbiome types were identified: Corynebacterium species (type I), Acinetobacter and Moraxella species (type II), Staphylococcus epidermidis (type III), Porphyromonas and Peptoniphilus species (type IV), and Propionibacterium acnes (type V). In lesional skin, microbiome types consisted predominantly of type I or type IV. Microbiome type IV was not detected in healthy controls. Several taxa, including Propionibacterium , showed a significantly higher relative abundance in healthy controls vs HS skin, indicating that Propionibacterium may be part of the pathogenesis in HS. Conclusions and Relevance The study findings suggest a link between a dysbiotic cutaneous microbiome and HS.

Eukaryotic RNAs with premature termination codons (PTCs) are eliminated by nonsense-mediated decay (NMD). While human nonsense RNA degradation can be initiated either by an endonucleolytic cleavage event near the PTC or through decapping, the individual contribution of these activities on endogenous substrates has remained unresolved. Here we used concurrent transcriptome-wide identification of NMD substrates and their 5′–3′ decay intermediates to establish that SMG6-catalyzed endonucleolysis widely initiates the degradation of human nonsense RNAs, whereas decapping is used to a lesser extent. We also show that a large proportion of genes hosting snoRNAs in their introns produce considerable amounts of NMD-sensitive splice variants, indicating that these RNAs are merely by-products of a primary snoRNA production process. Additionally, transcripts from genes encoding multiple snoRNAs often yield alternative transcript isoforms that allow for differential expression of individual coencoded snoRNAs. Based on our findings, we hypothesize that snoRNA host genes need to be highly transcribed to accommodate high levels of snoRNA production and that the expression of individual snoRNAs and their cognate spliced RNA can be uncoupled via alternative splicing and NMD.

Epigenetic information is available from contemporary organisms, but is difficult to track back in evolutionary time. Here, we show that genome-wide epigenetic information can be gathered directly from next-generation sequence reads of DNA isolated from ancient remains. Using the genome sequence data generated from hair shafts of a 4000-yr-old Paleo-Eskimo belonging to the Saqqaq culture, we generate the first ancient nucleosome map coupled with a genome-wide survey of cytosine methylation levels. The validity of both nucleosome map and methylation levels were confirmed by the recovery of the expected signals at promoter regions, exon/intron boundaries, and CTCF sites. The top-scoring nucleosome calls revealed distinct DNA positioning biases, attesting to nucleotide-level accuracy. The ancient methylation levels exhibited high conservation over time, clustering closely with modern hair tissues. Using ancient methylation information, we estimated the age at death of the Saqqaq individual and illustrate how epigenetic information can be used to infer ancient gene expression. Similar epigenetic signatures were found in other fossil material, such as 110,000- to 130,000-yr-old bones, supporting the contention that ancient epigenomic information can be reconstructed from a deep past. Our findings lay the foundation for extracting epigenomic information from ancient samples, allowing shifts in epialleles to be tracked through evolutionary time, as well as providing an original window into modern epigenomics.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

The diversity of OXA-48-like carbapenemases is continually expanding. In this study, we describe the dissemination and characteristics of a novel carbapenem-hydrolyzing class D β-lactamase (CHDL) named OXA-436. In total, six OXA-436-producing Enterobacteriaceae isolates, including Enterobacter asburiae ( n = 3), Citrobacter freundii ( n = 2), and Klebsiella pneumoniae ( n = 1), were identified in four patients in the period between September 2013 and April 2015. All three species of OXA-436-producing Enterobacteriaceae were found in one patient. The amino acid sequence of OXA-436 showed 90.4 to 92.8% identity to the amino acid sequences of other acquired OXA-48-like variants. Expression of OXA-436 in Escherichia coli and kinetic analysis of purified OXA-436 revealed an activity profile similar to that of OXA-48 and OXA-181, with activity against penicillins, including temocillin; limited or no activity against extended-spectrum cephalosporins; and activity against carbapenems. The bla OXA-436 gene was located on a conjugative ∼314-kb IncHI2/IncHI2A plasmid belonging to plasmid multilocus sequence typing sequence type 1 in a region surrounded by chromosomal genes previously identified to be adjacent to bla OXA genes in Shewanella spp. In conclusion, OXA-436 is a novel CHDL with functional properties similar to those of OXA-48-like CHDLs. The described geographical spread among different Enterobacteriaceae and the plasmid location of bla OXA-436 illustrate its potential for further dissemination.