Tracking bacterial pollution at a marine wastewater outfall site - A case study from Norway.

Coastal marine environments are increasingly affected by anthropogenic impacts, such as the release of sewage at outfall sites and agricultural run-off. Fecal pollution introduced to the sea through these activities poses risks of spreading microbial diseases and disseminating antibiotic resistant bacteria and their genes. The study area of this research, Bore beach, is situated between two such point sources, an outfall site where treated sewage is released 1 km off the coast and a stream that carries run-off from an agricultural area to the northern end of the beach. In order to investigate whether and to what extent fecal contamination from the sewage outfall reached the beach, we used microbial source tracking, based on whole community analysis. Samples were collected from sea water at varying distances from the sewage outfall site and along the beach, as well as from the sewage effluent and the stream. Amplicon sequencing of 16S rRNA genes from all the collected samples was carried out at two time points (June and September). In addition, the seawater at the sewage outfall site and the sewage effluent were subject to shotgun metagenomics. To estimate the contribution of the sewage effluent and the stream to the microbial communities at Bore beach, we employed SourceTracker2, a program that uses a Bayesian algorithm to perform such quantification. The SourceTracker2 results suggested that the sewage effluent is likely to spread fecal contamination towards the beach to a greater extent than anticipated based on the prevailing sea current. The estimated mixing proportions of sewage at the near-beach site (P4) were 0.22 and 0.035% in June and September, respectively. This was somewhat below that stream's contribution in June (0.028%) and 10-fold higher than the stream's contribution in September (0.004%). Our analysis identified a sewage signal in all the tested seawater samples.

Application of GCMS-pyrolysis to estimate the levels of microplastics in a drinking water supply system.

Communities value water and aquatic environments for a many diverse reasons. Ensuring safe drinking water is prioritized on the political agenda with a dedicated focus on safe and affordable drinking water under the 6th of the UN sustainable development goals. The occurrence of micron sized plastic fragments has been confirmed even in very remote areas. In the present study we analysed drinking water of a medium-sized Norwegian urban area for the presence of microplastics ≥ 1 µm. A modular filtering sampling devices was developed allowing a sequential in-situ enzymatic and mild oxidizing driven sample preparation prior to pyrolysis gas chromatography-mass spectrometry sample's analysis (pyr-GCMS). Samples were taken at different stages of the drinking water supply chain. The total amount of polymers per sites ranged from 6.1 to 93.1 µg/m3. Higher levels were detected in the raw water, but significant reduction rates ranging from 43% to 100% depending on the polymer type were scored after the water treatment processes. Polyethylene, polyamide, and polyester were the most frequently detected polymer types. Overall, the levels of MPs in the raw water influence the occurrence and polymer type occurrence and distribution is the drinking water supply net. This study contributes to the emerging field of plastics pollution in drinking water supply systems by providing effective methods helping with future routine monitoring of this source of human plastic uptake.

Computational identification of piRNA targets on mouse mRNAs.

MOTIVATION: PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are highly abundant in the germline. One important role of piRNAs is to defend genome integrity by guiding PIWI proteins to silence transposable elements (TEs), which have a high potential to cause deleterious effects on their host. The mechanism of piRNA-mediated post-transcriptional silencing was also observed to affect mRNAs, suggesting that piRNAs might play a broad role in gene expression regulation. However, there has been no systematic report with regard to how many protein-coding genes might be targeted and regulated by piRNAs. RESULTS: We trained a support vector machine classifier based on a combination of Miwi CLIP-Seq-derived features and position-derived features to predict the potential targets of piRNAs on mRNAs in the mouse. Reanalysis of a published microarray dataset suggested that the expression level of the 2587 protein-coding genes predicted as piRNA targets showed significant upregulation as a whole after abolishing the slicer activity of Miwi, supporting the conclusion that they are subject to piRNA-mediated regulation. AVAILABILITY AND IMPLEMENTATION: A web version of the method called pirnaPre as well as our results for browse is available at http://www.regulatoryrna.org/software/piRNA/piRNA_target_mRNA/index.php CONTACT: crs@sun5.ibp.ac.cn or heshunmin@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Functional diversity of CTCFs is encoded in their binding motifs.

BACKGROUND: The CCCTC-binding factor (CTCF) has diverse regulatory functions. However, the definitive characteristics of the CTCF binding motif required for its functional diversity still remains elusive. RESULTS: Here, we describe a new motif discovery workflow by which we have identified three CTCF binding motif variations with highly divergent functionalities. Supported by transcriptomic, epigenomic and chromatin-interactomic data, we show that the functional diversity of the CTCF binding motifs is strongly associated with their GC content, CpG dinucleotide coverage and relative DNA methylation level at the 12th position of the motifs. Further analysis suggested that the co-localization of cohesin, the key factor in cohesion of sister chromatids, is negatively correlated with the CpG coverage and the relative DNA methylation level at the 12th position. Finally, we present evidences for a hypothetical model in which chromatin interactions between promoters and distal regulatory regions are likely mediated by CTCFs binding to sequences with high CpG. CONCLUSION: These results demonstrate the existence of definitive CTCF binding motifs corresponding to CTCF's diverse functions, and that the functional diversity of the motifs is strongly associated with genetic and epigenetic features at the 12th position of the motifs.

RNAcentral: an international database of ncRNA sequences.

The field of non-coding RNA biology has been hampered by the lack of availability of a comprehensive, up-to-date collection of accessioned RNA sequences. Here we present the first release of RNAcentral, a database that collates and integrates information from an international consortium of established RNA sequence databases. The initial release contains over 8.1 million sequences, including representatives of all major functional classes. A web portal (http://rnacentral.org) provides free access to data, search functionality, cross-references, source code and an integrated genome browser for selected species.

piRBase: a web resource assisting piRNA functional study.

piRNAs are a class of small RNAs that is most abundantly expressed in the animal germ line. Presently, substantial research is going on to reveal the functions of piRNAs in the epigenetic and post-transcriptional regulation of transposons and genes. A piRNA database for collection, annotation and structuring of these data will be a valuable contribution to the field, and we have therefore developed the piRBase platform which integrates various piRNA-related high-throughput data. piRBase has the largest collection of piRNAs among existing databases, and contains at present 77 million piRNA sequences from nine organisms. Repeat-derived and gene-derived piRNAs, which possibly participate in the regulation of the corresponding elements, have been given particular attention. Furthermore, epigenetic data and reported piRNA targets were also collected. To our knowledge, this is the first piRNA database that systematically integrates epigenetic and post-transcriptional regulation data to support piRNA functional analysis. We believe that piRBase will contribute to a better understanding of the piRNA functions. Database URL: http://www.regulatoryrna.org/database/piRNA/

Analysis of the p53/CEP-1 regulated non-coding transcriptome in C. elegans by an NSR-seq strategy.

In recent years, large numbers of non-coding RNAs (ncRNAs) have been identified in C. elegans but their functions are still not well studied. In C. elegans, CEP-1 is the sole homolog of the p53 family of genes. In order to obtain transcription profiles of ncRNAs regulated by CEP-1 under normal and UV stressed conditions, we applied the 'not-so-random' hexamers priming strategy to RNA sequencing in C. elegans, This NSR-seq strategy efficiently depleted rRNA transcripts from the samples and showed high technical replicability. We identified more than 1,000 ncRNAs whose apparent expression was repressed by CEP-1, while around 200 were activated. Around 40% of the CEP-1 activated ncRNAs promoters contain a putative CEP-1-binding site. CEP-1 regulated ncRNAs were frequently clustered and concentrated on the X chromosome. These results indicate that numerous ncRNAs are involved in CEP-1 transcriptional network and that these are especially enriched on the X chromosome in C. elegans.

Evolutionary relationships between miRNA genes and their activity.

BACKGROUND: The emergence of vertebrates is characterized by a strong increase in miRNA families. MicroRNAs interact broadly with many transcripts, and the evolution of such a system is intriguing. However, evolutionary questions concerning the origin of miRNA genes and their subsequent evolution remain unexplained. RESULTS: In order to systematically understand the evolutionary relationship between miRNAs gene and their function, we classified human known miRNAs into eight groups based on their evolutionary ages estimated by maximum parsimony method. New miRNA genes with new functional sequences accumulated more dynamically in vertebrates than that observed in Drosophila. Different levels of evolutionary selection were observed over miRNA gene sequences with different time of origin. Most genic miRNAs differ from their host genes in time of origin, there is no particular relationship between the age of a miRNA and the age of its host genes, genic miRNAs are mostly younger than the corresponding host genes. MicroRNAs originated over different time-scales are often predicted/verified to target the same or overlapping sets of genes, opening the possibility of substantial functional redundancy among miRNAs of different ages. Higher degree of tissue specificity and lower expression level was found in young miRNAs. CONCLUSIONS: Our data showed that compared with protein coding genes, miRNA genes are more dynamic in terms of emergence and decay. Evolution patterns are quite different between miRNAs of different ages. MicroRNAs activity is under tight control with well-regulated expression increased and targeting decreased over time. Our work calls attention to the study of miRNA activity with a consideration of their origin time.

The human microbiome: a hot spot of microbial horizontal gene transfer.

The human body harbors numerous microbes, and here exists a close relationship between microbes and human health. The Human Microbiome Project has generated whole genome sequences of several hundred human microbes. In this study, we identified horizontal gene transfer (HGT) events in human microbes and tried to elucidate the relationships between the gene-transferring microbes. A total of 13,514 high confidence HGT genes were identified in 308 human microbes. The horizontally transferred genes were enriched for Gene Ontology terms pertaining to catalytic functions and metabolic processes. Construction of an HGT event network suggested that the human microbes could be divided into specific communities which only partly overlap their distribution in human body. Our research suggests that human microbiome may facilitate frequent horizontal gene transfer among bacteria in human body. Awareness of HGT in human microbiome may aid our understanding of the relationship between the human microbiome and human health.

Distinct MicroRNA Subcellular Size and Expression Patterns in Human Cancer Cells.

Introduction. Small noncoding RNAs have important regulatory functions in different cell pathways. It is believed that most of them mainly play role in gene post-transcriptional regulation in the cytoplasm. Recent evidence suggests miRNA and siRNA activity in the nucleus. Here, we show distinct genome-wide sub-cellular localization distribution profiles of small noncoding RNAs in human breast cancer cells. Methods. We separated breast cancer cell nuclei from cytoplasm, and identified small RNA sequences using a high-throughput sequencing platform. To determine the relationship between miRNA sub-cellular distribution and cancer progression, we used microarray analysis to examine the miRNA expression levels in nucleus and cytoplasm of three human cell lines, one normal breast cell line and two breast cancer cell lines. Logistic regression and SVM were used for further analysis. Results. The sub-cellular distribution of small noncoding RNAs shows that numerous miRNAs and their isoforms (isomiR) not only locate to the cytoplasm but also appeare in the nucleus. Subsequent microarray analyses indicated that the miRNA nuclear-cytoplasmic-ratio is a significant characteristic of different cancer cell lines. Conclusions. Our results indicate that the sub-cellular distribution is important for miRNA function, and that the characterization of the small RNAs sub-cellular localizome may contribute to cancer research and diagnosis.

Integrated sequence-structure motifs suffice to identify microRNA precursors.

BACKGROUND: Upwards of 1200 miRNA loci have hitherto been annotated in the human genome. The specific features defining a miRNA precursor and deciding its recognition and subsequent processing are not yet exhaustively described and miRNA loci can thus not be computationally identified with sufficient confidence. RESULTS: We rendered pre-miRNA and non-pre-miRNA hairpins as strings of integrated sequence-structure information, and used the software Teiresias to identify sequence-structure motifs (ss-motifs) of variable length in these data sets. Using only ss-motifs as features in a Support Vector Machine (SVM) algorithm for pre-miRNA identification achieved 99.2% specificity and 97.6% sensitivity on a human test data set, which is comparable to previously published algorithms employing combinations of sequence-structure and additional features. Further analysis of the ss-motif information contents revealed strongly significant deviations from those of the respective training sets, revealing important potential clues as to how the sequence and structural information of RNA hairpins are utilized by the miRNA processing apparatus. CONCLUSION: Integrated sequence-structure motifs of variable length apparently capture nearly all information required to distinguish miRNA precursors from other stem-loop structures.

A differential sequencing-based analysis of the C. elegans noncoding transcriptome.

Noncoding RNAs are increasingly being recognized as important players in eukaryote biology. However, despite major efforts in mapping the Caenorhabditis elegans transcriptome over the last couple of years, nonpolyadenylated and intermediate-size noncoding RNAs (is-ncRNAs) are still incompletely explored. We have combined an enzymatic approach with full-length RNA-Seq of is-ncRNAs in C. elegans. A total of 473 novel is-ncRNAs has been identified, of which a substantial fraction was associated with transcription factor binding sites and developmentally regulated expression patterns. Analysis of sequence and secondary structure permitted classification of more than 200 is-ncRNAs into several known RNA classes, while another 33 is-ncRNAs were identified as belonging to two previously uncharacterized groups of is-ncRNAs. Three of the unclassified is-ncRNAs contain the 5' Alu domain common to SRP RNAs and specifically bound with the SRP9/14 heterodimer in vitro. One of these (inc394) showed 65% sequence identity with the human, neuron-specific BC200 RNA. Structure-based clustering analysis and in vitro binding experiments supported the notion that the nematode stem-bulge RNAs (sbRNAs) are homologs (or functional analogs) of the Y RNAs. Moreover, analysis of the differential libraries showed that some mature snoRNAs undergo secondary 5' cap modification after processing of the primary transcript, thus suggesting the existence of a wider range of functional RNAs arising from processed and modified fragments of primary transcripts.

NONCODE v3.0: integrative annotation of long noncoding RNAs.

Facilitated by the rapid progress of high-throughput sequencing technology, a large number of long noncoding RNAs (lncRNAs) have been identified in mammalian transcriptomes over the past few years. LncRNAs have been shown to play key roles in various biological processes such as imprinting control, circuitry controlling pluripotency and differentiation, immune responses and chromosome dynamics. Notably, a growing number of lncRNAs have been implicated in disease etiology. With the increasing number of published lncRNA studies, the experimental data on lncRNAs (e.g. expression profiles, molecular features and biological functions) have accumulated rapidly. In order to enable a systematic compilation and integration of this information, we have updated the NONCODE database (http://www.noncode.org) to version 3.0 to include the first integrated collection of expression and functional lncRNA data obtained from re-annotated microarray studies in a single database. NONCODE has a user-friendly interface with a variety of search or browse options, a local Genome Browser for visualization and a BLAST server for sequence-alignment search. In addition, NONCODE provides a platform for the ongoing collation of ncRNAs reported in the literature. All data in NONCODE are open to users, and can be downloaded through the website or obtained through the SOAP API and DAS services.

ncFANs: a web server for functional annotation of long non-coding RNAs.

Recent interest in the non-coding transcriptome has resulted in the identification of large numbers of long non-coding RNAs (lncRNAs) in mammalian genomes, most of which have not been functionally characterized. Computational exploration of the potential functions of these lncRNAs will therefore facilitate further work in this field of research. We have developed a practical and user-friendly web interface called ncFANs (non-coding RNA Function ANnotation server), which is the first web service for functional annotation of human and mouse lncRNAs. On the basis of the re-annotated Affymetrix microarray data, ncFANs provides two alternative strategies for lncRNA functional annotation: one utilizing three aspects of a coding-non-coding gene co-expression (CNC) network, the other identifying condition-related differentially expressed lncRNAs. ncFANs introduces a highly efficient way of re-using the abundant pre-existing microarray data. The present version of ncFANs includes re-annotated CDF files for 10 human and mouse Affymetrix microarrays, and the server will be continuously updated with more re-annotated microarray platforms and lncRNA data. ncFANs is freely accessible at http://www.ebiomed.org/ncFANs/ or http://www.noncode.org/ncFANs/.

The Caenorhabditis elegans intermediate-size transcriptome shows high degree of stage-specific expression.

Earlier studies have revealed a substantial amount of transcriptional activity occurring outside annotated protein-coding genes of the Caenorhabditis elegans genome. One important fraction of this transcriptional activity relates to intermediate-size (70-500 nt) transcripts (is-ncRNAs) of mostly unknown function. Profiling the expression of this segment of the transcriptome on a tiling array through the C. elegans life cycle identified 5866 hitherto unannotated transcripts. The novel loci were distributed across intronic and intergenic space, with some enrichment toward protein-coding gene termini. The majority of the putative is-ncRNAs showed either stage-specific expression, or distinct developmental variation in their expression levels. More than 200 loci showed male-specific expression, and conserved loci were significantly enriched on the X chromosome, both observations strongly suggesting involvement of is-ncRNAs in sex-specific functions. Half of the novel loci were conserved in other nematodes, and numerous loci showed significant conservational correlations to nearby coding genes. Assuming functional roles for most of the novel loci, the data imply a nematode is-ncRNA tool kit of considerable size and variety.

Large-scale prediction of long non-coding RNA functions in a coding-non-coding gene co-expression network.

Although accumulating evidence has provided insight into the various functions of long-non-coding RNAs (lncRNAs), the exact functions of the majority of such transcripts are still unknown. Here, we report the first computational annotation of lncRNA functions based on public microarray expression profiles. A coding-non-coding gene co-expression (CNC) network was constructed from re-annotated Affymetrix Mouse Genome Array data. Probable functions for altogether 340 lncRNAs were predicted based on topological or other network characteristics, such as module sharing, association with network hubs and combinations of co-expression and genomic adjacency. The functions annotated to the lncRNAs mainly involve organ or tissue development (e.g. neuron, eye and muscle development), cellular transport (e.g. neuronal transport and sodium ion, acid or lipid transport) or metabolic processes (e.g. involving macromolecules, phosphocreatine and tyrosine).

Systematic identification and characterization of chicken (Gallus gallus) ncRNAs.

Recent studies have demonstrated that non-coding RNAs (ncRNAs) play important roles during development and evolution. Chicken, the first genome-sequenced non-mammalian amniote, possesses unique features for developmental and evolutionary studies. However, apart from microRNAs, information on chicken ncRNAs has mainly been obtained from computational predictions without experimental validation. In the present study, we performed a systematic identification of intermediate size ncRNAs (50-500 nt) by ncRNA library construction and identified 125 chicken ncRNAs. Importantly, through the bioinformatics and expression analysis, we found the chicken ncRNAs has several novel features: (i) comparative genomic analysis against 18 sequenced vertebrate genomes revealed that the majority of the newly identified ncRNA candidates is not conserved and most are potentially bird/chicken specific, suggesting that ncRNAs play roles in lineage/species specification during evolution. (ii) The expression pattern analysis of intronic snoRNAs and their host genes suggested the coordinated expression between snoRNAs and their host genes. (iii) Several spatio-temporal specific expression patterns suggest involvement of ncRNAs in tissue development. Together, these findings provide new clues for future functional study of ncRNAs during development and evolution.

Genomic regions with distinct genomic distance conservation in vertebrate genomes.

BACKGROUND: A number of vertebrate highly conserved elements (HCEs) have been detected and their genomic interval distances have been reported to be more conserved than protein coding genes among mammalian genomes. A characteristic of the human - non-mammalian comparisons is a bimodal distribution of relative distance difference of conserved consecutive HCE pairs; and it is difficult to attribute such profile to a random assortment. We therefore undertook an analysis of the human genomic regions confined by consecutive HCE pairs common to eight genomes (human, mouse, rat, chicken, frog, zebrafish, tetradon and fugu). RESULTS: Among HCE pairs, we found that some consistently preserve highly conserved interval distance among genomes while others have relatively low distance conservation. Using a partition method, we detected two groups of inter-HCE regions (IHRs) with distinct distance conservation pattern in vertebrate genomes: IHR1s that are bordered by HCE pairs with relative small distance variation, and IHR2s with larger distance difference values. Compared to random background, annotated repeat sequences are significantly less frequent in IHR1s than IHR2s, which reflects a correlation between repeat sequences and the length expansion of IHRs. Both groups of IHRs are unexpectedly enriched in human indel (i.e. insertion and deletion) polymorphism-variations than random background. The correlation between the percentage of conserved sequence and human IHR length was stronger for IHR1 than IHR2. Both groups of IHRs are significantly enriched for CpG islands. CONCLUSION: The data suggest that subsets of HCE pairs may undergo different evolutionary paths in light of their genomic distance conservation, and that sets of genomic regions pertain to HCEs, as well as the region in which HCEs reside, should be treated as integrated domains.

Genome-scale identification of Caenorhabditis elegans regulatory elements by tiling-array mapping of DNase I hypersensitive sites.

BACKGROUND: A major goal of post-genomics research is the integrated analysis of genes, regulatory elements and the chromatin architecture on a genome-wide scale. Mapping DNase I hypersensitive sites within the nuclear chromatin is a powerful and well-established method of identifying regulatory element candidates. RESULTS: Here, we report the first genome-wide analysis of DNase I hypersensitive sites (DHSs) in Caenorhabditis elegans. The data was obtained by hybridizing DNase I-treated and end-captured material from young adult worms to a high-resolution tiling microarray. The data show that C. elegans DHSs were significantly enriched within intergenic regions located 2 kb upstream and downstream of coding genes, and also that a considerable fraction of all DHSs mapped to intergenic positions distant to annotated coding genes. Annotated transcribed loci were generally depleted in DHSs relative to intergenic regions, but DHSs were nonetheless enriched in coding exons and UTRs, whereas introns were significantly depleted in DHSs. Many DHSs appeared to be associated with annotated non-coding RNAs and recently detected transcripts of unknown function. It has been reported that nematode highly conserved non-coding elements were associated with cis-regulatory elements, and we also found that DHSs, particularly distal intergenic DHSs, were significantly enriched in regions that were conserved between the C. elegans and C. briggsae genomes. CONCLUSION: We describe the first genome-wide analysis of C. elegans DHSs, and show that the distribution of DHSs is strongly associated with functional elements in the genome.

Transcriptional inhibiton of Hoxd4 expression by miRNA-10a in human breast cancer cells.

BACKGROUND: Small noncoding RNAs (ncRNAs), including short interfering RNAs (siRNAs) and microRNAs (miRNAs), can silence genes at the transcriptional, post-transcriptional or translational level 12. RESULTS: Here, we show that microRNA-10a (miR-10a) targets a homologous DNA region in the promoter region of the hoxd4 gene and represses its expression at the transcriptional level. Mutational analysis of the miR-10a sequence revealed that the 3' end of the miRNA sequence is the most critical element for the silencing effect. MicroRNA-10a-induced transcriptional gene inhibition requires the presence of Dicer and Argonautes 1 and 3, and it is related to promoter associated noncoding RNAs. Bisulfite sequencing analysis showed that the reduced hoxd4 expression was accompanied by de novo DNA methylation at the hoxd4 promoter. We further demonstrated that trimethylation of histone 3 lysine 27 (H3K27me3) is involved in the miR-10a-induced hoxd4 transcriptional gene silence. CONCLUSION: In conclusion, our results demonstrate that miR-10a can regulate human gene expression in a transcriptional manner, and indicate that endogenous small noncoding RNA-induced control of transcription may be a potential system for expressional regulation in human breast cancer cells.