RNA sequencing studies have identified hundreds of non‐coding RNAs in bacteria, including regulatory small RNA (sRNA). However, our understanding of sRNA function has lagged behind their identification due to a lack of tools for the high‐throughput analysis of RNA–RNA interactions in bacteria. Here we demonstrate that in vivo sRNA–mRNA duplexes can be recovered using UV‐crosslinking, ligation and sequencing of hybrids (CLASH). Many sRNAs recruit the endoribonuclease, RNase E, to facilitate processing of mRNAs. We were able to recover base‐paired sRNA–mRNA duplexes in association with RNase E, allowing proximity‐dependent ligation and sequencing of cognate sRNA–mRNA pairs as chimeric reads. We verified that this approach captures bona fide sRNA–mRNA interactions. Clustering analyses identified novel sRNA seed regions and sets of potentially co‐regulated target mRNAs. We identified multiple mRNA targets for the pathotype‐specific sRNA Esr41, which was shown to regulate colicin sensitivity and iron transport in E. coli. Numerous sRNA interactions were also identified with non‐coding RNAs, including sRNAs and tRNAs, demonstrating the high complexity of the sRNA interactome.

BackgroundThe cane toad (Rhinella marina formerly Bufo marinus) is a species native to Central and South America that has spread across many regions of the globe. Cane toads are known for their rapid adaptation and deleterious impacts on native fauna in invaded regions. However, despite an iconic status, there are major gaps in our understanding of cane toad genetics. The availability of a genome would help to close these gaps and accelerate cane toad research.FindingsWe report a draft genome assembly for R. marina, the first of its kind for the Bufonidae family. We used a combination of long-read Pacific Biosciences RS II and short-read Illumina HiSeq X sequencing to generate 359.5 Gb of raw sequence data. The final hybrid assembly of 31,392 scaffolds was 2.55 Gb in length with a scaffold N50 of 168 kb. BUSCO analysis revealed that the assembly included full length or partial fragments of 90.6% of tetrapod universal single-copy orthologs (n = 3950), illustrating that the gene-containing regions have been well assembled. Annotation predicted 25,846 protein coding genes with similarity to known proteins in Swiss-Prot. Repeat sequences were estimated to account for 63.9% of the assembly.ConclusionsThe R. marina draft genome assembly will be an invaluable resource that can be used to further probe the biology of this invasive species. Future analysis of the genome will provide insights into cane toad evolution and enrich our understanding of their interplay with the ecosystem at large.

Nanopore sequencing depends on the FAST5 file format, which does not allow efficient parallel analysis. Here we introduce SLOW5, an alternative format engineered for efficient parallelization and acceleration of nanopore data analysis. Using the example of DNA methylation profiling of a human genome, analysis runtime is reduced from more than two weeks to approximately 10.5 h on a typical high-performance computer. SLOW5 is approximately 25% smaller than FAST5 and delivers consistent improvements on different computer architectures.

BackgroundAlthough high risk HPVs are associated with an increased risk of prostate cancer it is not known if they have a causal role. The purpose of this study is to investigate the potential role of human papilloma viruses (HPVs) in prostate cancer. The aims are (i) to investigate the presence and confirm the identity of high risk HPVs in benign prostate tissues prior to the development of HPV positive prostate cancer in the same patients, and (ii) to determine if HPVs are biologically active.MethodsWe used polymerase chain reaction (PCR) to identify HPVs in specimens from 52 Australian men with benign prostate biopsies who 1 to 10 years later developed prostate cancer. Immunohistochemistry (IHC) was used to assess the expression of HPV E7 oncoproteins, cytokeratin and prostate specific antigen (PSA).We used RNASeq data from The Cancer Genome Atlas (TCGA) to identify possible HPV RNA sequences in prostate cancer.ResultsHPV screening using standard PCR was conducted on 28 of the 52 sets of benign and later prostate cancers. HPV L1 genes were identified in 13 (46%) benign and 8 (29%) of 28 later prostate cancers in the same patients. HPV E7 genes were identified in 23 (82%) benign and 19 (68%) of 28 subsequent prostate cancers in the same patients. The same HPV types were present in both the benign and subsequent prostate cancers in 9 sets of specimens. HPV type 16 was identified in 15% of benign and 3% of prostate cancers. HPV type 18 was identified in 26% of benign and 16% of prostate cancers. Small numbers of HPV types 45, 47, 76 and 115 were also identified.High confidence RNA-Seq evidence for high risk HPV types 16 and 18 was identified in 12 (2%) of the 502 TCGA prostate cancer transcriptomes.High risk HPV E7 oncoprotein was positively expressed in 23 (82%) of 28 benign prostate specimens but only in 8 (29%) of 28 of the later prostate cancer specimens. This difference is statistically significant (p = 0.001). Prostate specific antigen (PSA) was more highly expressed in 26 (50%) of 52 prostate cancer specimens as compared to prior benign prostate specimens in the same patients.ConclusionsHigh risk HPVs are present in benign prostate tissues prior to the development of HPV positive prostate cancer. There is a significantly higher expression of HPV E7 oncoproteins in benign prostate tissues as compared to late prostate cancer that subsequently developed in the same patients. This observation suggests that HPV oncogenic activity is an early phenomenon in a majority of prostate oncogenesis. TCGA RNA-Seq data suggests that HPV is biologically active in some prostate tumour samples.Electronic supplementary materialThe online version of this article (doi:10.1186/s13027-017-0157-2) contains supplementary material, which is available to authorized users.

BackgroundMECP2, the gene mutated in the majority of Rett syndrome cases, is a transcriptional regulator that can activate or repress transcription. Although the transcription regulatory function of MECP2 has been known for over a decade, it remains unclear how transcriptional dysregulation leads to the neurodevelopmental disorder. Notably, little convergence was previously observed between the genes abnormally expressed in the brain of Rett syndrome mouse models and those identified in human studies.MethodsHere we carried out a comprehensive transcriptome analysis of human brain tissue from Rett syndrome brain using both RNA-seq and microarrays.ResultsWe identified over two hundred differentially expressed genes, and identified the complement C1Q complex genes (C1QA, C1QB and C1QC) as a point of convergence between gene expression changes in human and mouse Rett syndrome brain.ConclusionsThe results of our study support a role for alterations in the expression level of C1Q complex genes in RTT pathogenesis.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2746-7) contains supplementary material, which is available to authorized users.