Approaches for sRNA Analysis of Human RNA-Seq Data: Comparison, Benchmarking.

Expression analysis of small noncoding RNA (sRNA), including microRNA, piwi-interacting RNA, small rRNA-derived RNA, and tRNA-derived small RNA, is a novel and quickly developing field. Despite a range of proposed approaches, selecting and adapting a particular pipeline for transcriptomic analysis of sRNA remains a challenge. This paper focuses on the identification of the optimal pipeline configurations for each step of human sRNA analysis, including reads trimming, filtering, mapping, transcript abundance quantification and differential expression analysis. Based on our study, we suggest the following parameters for the analysis of human sRNA in relation to categorical analyses with two groups of biosamples: (1) trimming with the lower length bound = 15 and the upper length bound = Read length - 40% Adapter length; (2) mapping on a reference genome with bowtie aligner with one mismatch allowed (-v 1 parameter); (3) filtering by mean threshold > 5; (4) analyzing differential expression with DESeq2 with adjusted p-value < 0.05 or limma with p-value < 0.05 if there is very little signal and few transcripts.

ITAS: Integrated Transcript Annotation for Small RNA.

Transcriptomics analysis of various small RNA (sRNA) biotypes is a new and rapidly developing field. Annotations for microRNAs, tRNAs, piRNAs and rRNAs contain information on transcript sequences and loci that is vital for downstream analyses. Several databases have been established to provide this type of data for specific RNA biotypes. However, these sources often contain data in different formats, which makes the bulk analysis of several sRNA biotypes in a single pipeline challenging. Information on some transcripts may be incomplete or conflicting with other entries. To overcome these challenges, we introduce ITAS, or Integrated Transcript Annotation for Small RNA, a filtered, corrected and integrated transcript annotation containing information on several types of small RNAs, including tRNA-derived small RNA, for several species (Homo sapiens, Rattus norvegicus, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans). ITAS is presented in a format applicable for the vast majority of bioinformatic transcriptomics analysis, and it was tested in several case studies for human-derived data against existing alternative databases.

Optimization of small RNA extraction and comparative study of NGS library preparation from low count sperm samples.

Recent studies demonstrate that sperm epigenome is a vehicle that conveys paternal experiences to offspring phenotype. That evidence triggers interest of both experimental and epidemiological studies of epigenetic markers in sperm. Since samples are often unique in epidemiological studies, a careful and efficient use of the material is a critical requirement. The goal of this study was to provide optimization of methods for the isolation of small RNAs from spermatozoa and library preparation for sequencing. A total 67 fractionated sperm samples from the Russian Children's Study biobank prospectively collected at 18-20 years of age were used to isolate small RNAs with median (IQR) input total sperm count 17.0 (7.4-35.9) million. Twenty-four pairs of libraries were prepared using the NEBNext and NEXTFlex kits, 19 libraries using NEBNext and 6 using NEXTFlex. All libraries were sequenced on NextSeq 500, and the results were evaluated as a function of the number of small non-coding RNA (sncRNA) detected, quality parameters of sequencing libraries, as well as technical features of sample preparation. Although the same amount of miRNA input was used for NEBNext and NEXTFlex libraries, the concentration of DNA in NEBNext libraries was significantly higher in comparison with NEXTFlex libraries. In high input (sperm count >28 million and more than 25 ng miRNA in library) NEXTFlex Small RNA-Seq kit detected more microRNAs. In low input, the NEBNext proved more effective. The tricks and traps to protocol optimization are presented, including an efficient and effector gel-based system for the removal of sequencing library adaptors.

Different genome-wide transcriptome responses of Nocardioides simplex VKM Ac-2033D to phytosterol and cortisone 21-acetate.

BACKGROUND: Bacterial degradation/transformation of steroids is widely investigated to create biotechnologically relevant strains for industrial application. The strain of Nocardioides simplex VKM Ac-2033D is well known mainly for its superior 3-ketosteroid Δ1-dehydrogenase activity towards various 3-oxosteroids and other important reactions of sterol degradation. However, its biocatalytic capacities and the molecular fundamentals of its activity towards natural sterols and synthetic steroids were not fully understood. In this study, a comparative investigation of the genome-wide transcriptome profiling of the N. simplex VKM Ac-2033D grown on phytosterol, or in the presence of cortisone 21-acetate was performed with RNA-seq. RESULTS: Although the gene patterns induced by phytosterol generally resemble the gene sets involved in phytosterol degradation pathways in mycolic acid rich actinobacteria such as Mycolicibacterium, Mycobacterium and Rhodococcus species, the differences in gene organization and previously unreported genes with high expression level were revealed. Transcription of the genes related to KstR- and KstR2-regulons was mainly enhanced in response to phytosterol, and the role in steroid catabolism is predicted for some dozens of the genes in N. simplex. New transcription factors binding motifs and new candidate transcription regulators of steroid catabolism were predicted in N. simplex. Unlike phytosterol, cortisone 21-acetate does not provide induction of the genes with predicted KstR and KstR2 sites. Superior 3-ketosteroid-Δ1-dehydrogenase activity of N. simplex VKM Ac-2033D is due to the kstDs redundancy in the genome, with the highest expression level of the gene KR76_27125 orthologous to kstD2, in response to cortisone 21-acetate. The substrate spectrum of N. simplex 3-ketosteroid-Δ1-dehydrogenase was expanded in this study with progesterone and its 17α-hydroxylated and 11α,17α-dihydroxylated derivatives, that effectively were 1(2)-dehydrogenated in vivo by the whole cells of the N. simplex VKM Ac-2033D. CONCLUSION: The results contribute to the knowledge of biocatalytic features and diversity of steroid modification capabilities of actinobacteria, defining targets for further bioengineering manipulations with the purpose of expansion of their biotechnological applications.

Aging-induced changes in sperm DNA methylation are modified by low dose of perinatal flame retardants.

Aims: Paternal age is increasing in developed countries. Understanding of aging-related epigenetic changes in sperm is needed as well as factors that modify such changes. Materials & methods: Young pubertal and mature rats were exposed perinatally to vehicle or environmental xenobiotic 2,2',4,4'-tetrabromodiphenyl ether. Epididymal sperm was reduced representation bisulfite sequenced. Differentially methylated regions (DMRs) were identified via MethPipe. Results: In control animals, 5319 age-dependent DMRs were identified. Age-related DMRs were enriched for embryonic development. In exposed rats, DNA methylation was higher in young and lower in mature animals then in controls. Conclusions: Sperm methylome undergoes significant age-dependent changes, which may represent a causal link between paternal age and offspring phenotype. Environmental xenobiotics can interfere with the natural process of epigenetic aging.

Aging Induces Profound Changes in sncRNA in Rat Sperm and These Changes Are Modified by Perinatal Exposure to Environmental Flame Retardant.

Advanced paternal age at fertilization is a risk factor for multiple disorders in offspring and may be linked to age-related epigenetic changes in the father's sperm. An understanding of aging-related epigenetic changes in sperm and environmental factors that modify such changes is needed. Here, we characterize changes in sperm small non-coding RNA (sncRNA) between young pubertal and mature rats. We also analyze the modification of these changes by exposure to environmental xenobiotic 2,2',4,4'-tetrabromodiphenyl ether (BDE-47). sncRNA libraries prepared from epididymal spermatozoa were sequenced and analyzed using DESeq 2. The distribution of small RNA fractions changed with age, with fractions mapping to rRNA and lncRNA decreasing and fractions mapping to tRNA and miRNA increasing. In total, 249 miRNA, 908 piRNA and 227 tRNA-derived RNA were differentially expressed (twofold change, false discovery rate (FDR) p ≤ 0.05) between age groups in control animals. Differentially expressed miRNA and piRNA were enriched for protein-coding targets involved in development and metabolism, while piRNA were enriched for long terminal repeat (LTR) targets. BDE-47 accelerated age-dependent changes in sncRNA in younger animals, decelerated these changes in older animals and increased the variance in expression of all sncRNA. Our results indicate that the natural aging process has profound effects on sperm sncRNA profiles and this effect may be modified by environmental exposure.

Genome-Wide Transcriptome Profiling Provides Insight on Cholesterol and Lithocholate Degradation Mechanisms in Nocardioides simplex VKM Ac-2033D.

Steroid microbial degradation plays a significant ecological role for biomass decomposition and removal/detoxification of steroid pollutants. In this study, the initial steps of cholesterol degradation and lithocholate bioconversion by a strain with enhanced 3-ketosteroid dehydrogenase (3-KSD) activity, Nocardioides simplex VKM Ac-2033D, were studied. Biochemical, transcriptomic, and bioinformatic approaches were used. Among the intermediates of sterol sidechain oxidation cholest-5-en-26-oic acid and 3-oxo-cholesta-1,4-dien-26-oic acid were identified as those that have not been earlier reported for N. simplex and related species. The transcriptomic approach revealed candidate genes of cholesterol and lithocholic acid (LCA) catabolism by the strain. A separate set of genes combined in cluster and additional 3-ketosteroid Δ1-dehydrogenase and 3-ketosteroid 9α-hydroxylases that might be involved in LCA catabolism were predicted. Bioinformatic calculations based on transcriptomic data showed the existence of a previously unknown transcription factor, which regulates cholate catabolism gene orthologs. The results contribute to the knowledge on diversity of steroid catabolism regulation in actinobacteria and might be used at the engineering of microbial catalysts for ecological and industrial biotechnology.

Complete genome assembly data of paenibacillus sp. RUD330, a hypothetical symbiont of euglena gracilis.

An unknown bacterial strain was detected in the cytostome of Euglena gracilis and on the cell surface of Euglena gracilis using transmission electron microscopy. To identify the unknown bacterium and its function, we performed isolation experiments. Here we present the genome sequence of the isolate that was determined to be Paenibacillus sp. The genome of the bacterium was sequenced four times using Illumina technology with pair-end reads, Illumina technology with mate pair reads (inserts 3-4 and 6-8 Kb), and Nanopore technology with long reads (tens of thousands of nucleotides). Assemblies based on Illumina reads including mate-pair reads could not resolve issues caused by long tandem copies of rRNA, other tandem repeats, and extremely GC-rich regions (90-100%). Only long Nanopore reads resolved those gaps and made it possible to complete the entire genome; moreover, we found one plasmid. The length of the genome is 5.56 Mbp, and the average GC content is 59%. The genome of Paenibacillus sp. RUD330 included 8 copies of all the rRNA genes (23S; 16S; 5S), the length of the plasmid was 8.3 Kb. We hope that our genome assembly and the methods used can help other investigators in the assembly of complex genomes. Our reliable assembly could be a good basis for further physiological and genetic engineering studies of similar strains.

Rat liver epigenome programing by perinatal exposure to 2,2',4'4'-tetrabromodiphenyl ether.

Perinatal exposures to polybrominated diphenyl ethers permanently reprogram liver metabolism and induce a nonalcoholic fatty liver disease-like phenotype and insulin resistance in rodents. Aim: To test if these changes are associated with altered liver epigenome. Materials & methods: Expression of small RNA and changes in DNA methylation in livers of adult rats were analyzed following perinatal exposure to 2,2',4,4'-tetrabromodiphenyl ether, the polybrominated diphenyl ether congener most prevalent in human tissues. Results: We identified 33 differentially methylated DNA regions and 15 differentially expressed miRNAs. These changes were enriched for terms related to lipid and carbohydrate metabolism, insulin signaling, Type-2 diabetes and nonalcoholic fatty liver disease. Conclusion: Changes in the liver epigenome are a likely candidate mechanism of long-term maintenance of an aberrant metabolic phenotype.

Genome-wide response on phytosterol in 9-hydroxyandrostenedione-producing strain of Mycobacterium sp. VKM Ac-1817D.

BACKGROUND: Aerobic side chain degradation of phytosterols by actinobacteria is the basis for the industrial production of androstane steroids which are the starting materials for the synthesis of steroid hormones. A native strain of Mycobacterium sp. VKM Ac-1817D effectively produces 9α-hydroxyandrost-4-ene-3,17-dione (9-OH-AD) from phytosterol, but also is capable of slow steroid core degradation. However, the set of the genes with products that are involved in phytosterol oxidation, their organisation and regulation remain poorly understood. RESULTS: High-throughput sequencing of the global transcriptomes of the Mycobacterium sp. VKM Ac-1817D cultures grown with or without phytosterol was carried out. In the presence of phytosterol, the expression of 260 genes including those related to steroid catabolism pathways significantly increased. Two of the five genes encoding the oxygenase unit of 3-ketosteroid-9α-hydroxylase (kshA) were highly up-regulated in response to phytosterol (55- and 25-fold, respectively) as well as one of the two genes encoding its reductase subunit (kshB) (40-fold). Only one of the five putative genes encoding 3-ketosteroid-∆1-dehydrogenase (KstD_1) was up-regulated in the presence of phytosterol (61-fold), but several substitutions in the conservative positions of its product were revealed. Among the genes over-expressed in the presence of phytosterol, several dozen genes did not possess binding sites for the known regulatory factors of steroid catabolism. In the promoter regions of these genes, a regularly occurring palindromic motif was revealed. The orthologue of TetR-family transcription regulator gene Rv0767c of M. tuberculosis was identified in Mycobacterium sp. VKM Ac-1817D as G155_05115. CONCLUSIONS: High expression levels of the genes related to the sterol side chain degradation and steroid 9α-hydroxylation in combination with possible defects in KstD_1 may contribute to effective 9α-hydroxyandrost-4-ene-3,17-dione accumulation from phytosterol provided by this biotechnologically relevant strain. The TetR-family transcription regulator gene G155_05115 presumably associated with the regulation of steroid catabolism. The results are of significance for the improvement of biocatalytic features of the microbial strains for the steroid industry.

Peripubertal serum dioxin concentrations and subsequent sperm methylome profiles of young Russian adults.

BACKGROUND: The association of exposure to endocrine disrupting chemicals in the peripubertal period with subsequent sperm DNA methylation is unknown. OBJECTIVE: We examined the association of peripubertal serum 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) concentrations with whole-genome bisulfite sequencing (WGBS) of sperm collected in young adulthood. METHODS: The Russian Children's Study is a prospective cohort of 516 boys who were enrolled at 8-9 years of age and provided semen samples at 18-19 years of age. WGBS of sperm was conducted to identify differentially methylated regions (DMR) between highest (n = 4) and lowest (n = 4) peripubertal TCDD groups. RESULTS: We found 52 DMRs that distinguished lowest and highest peripubertal serum TCDD concentrations. One of the top scoring networks, "Cellular Assembly and Organization, Cellular Function and Maintenance, Carbohydrate Metabolism", identified estrogen receptor alpha as its central regulator. CONCLUSION: Findings from our limited sample size suggest that peripubertal environmental exposures are associated with sperm DNA methylation in young adults.

The complete genome of the oil emulsifying strain Thalassolituus oleivorans K-188 from the Barents Sea.

Gammaproteobacterium Thalassolituus oleivorans plays an important role in oil degradation in sea water through emulsifying crude oil and alkanes at low temperatures in polar sea environment. Here we report the complete genome sequence of K-188 strain (VKPM B-9394) isolated in the Barents Sea and compare it with other known Thalassolituus oleivorans strains. The Thalassolituus strains are differed in orthologs number of the genes of alkane degradation, transport proteins, genes of sugar utilization, endonucleases, signaling proteins, transcriptional regulators and presence of CRISPR/Cas locus. Also only the genome of K-188 contains the 3-hydroxyalkanoate synthetase.

Genome Sequencing of Steroid-Producing Bacteria with Illumina Technology.

Illumina technology is widely used for bacterial whole-genome sequencing due to its simplicity, cheapness, reliability, and abundant software for manipulation with raw data. Illumina technology belongs to a second generation of whole genome sequencing that yields great amount of short reads for genome regions. Genomic DNA is fragmented to short pieces. DNA fragments are amplified for signal increasing, and are read using sequencing-by-synthesis. Millions of short reads up to 100-300 bp in length are assembled in continuous sequences. Mate-pair technology allows resolving a long repeat.Here, we describe the principles of standard and mate-pair library preparation from DNA samples, library quality control, sequencing with MiSeq instrument and following data bioinformatics treatment. Software for genome assembly and completion are listed that assemble, map, annotate, visualize, edit and allow doing other manipulations with genomic sequences. The whole genomes sequencing of the steroid-producing Actinobacteria using these protocols is exemplified.

Effect of methyl-ß-cyclodextrin on gene expression in microbial conversion of phytosterol.

Modified ß-cyclodextrins are widely used for the enhancement of microbial conversions of lipophilic compounds such as steroids. Multiple mechanisms of cyclodextrin-mediated enhancement of phytosterol bioconversion by mycobacteria had previously been shown to include steroid solubilization, alterations in the cell wall permeability for both steroids and nutrients, facilitation of protein leaking, and activity suppression of some steroid-transforming enzymes.In this work, we studied whether cyclodextrins might affect expression of the genes involved in the steroid catabolic pathway. Phytosterol bioconversion with 9α-hydroxy-androst-4-ene-3,17-dione accumulation by Mycobacterium sp. VKM Ac-1817D in the presence of methylated ß-cyclodextrin (MCD) was investigated. RNA sequencing of the whole transcriptomes in different combinations of phytosterol and MCD showed a similar expression level of the steroid catabolism genes related to the KstR-regulon and was responsible for side chain and initial steps of steroid core oxidation; whereas, induction levels of the genes related to the KstR2-regulon were attenuated in the presence of MCD in this strain. The data were attenuated with quantitative real-time PCR.The results contribute to the understanding of cyclodextrin effects on microbial steroid conversion and provide a basis for the use of cyclodextrins as expression enhancers for studies of sterol catabolism in actinobacteria.

Comparative analysis of plastid genomes of non-photosynthetic Ericaceae and their photosynthetic relatives.

Although plastid genomes of flowering plants are typically highly conserved regarding their size, gene content and order, there are some exceptions. Ericaceae, a large and diverse family of flowering plants, warrants special attention within the context of plastid genome evolution because it includes both non-photosynthetic and photosynthetic species with rearranged plastomes and putative losses of "essential" genes. We characterized plastid genomes of three species of Ericaceae, non-photosynthetic Monotropa uniflora and Hypopitys monotropa and photosynthetic Pyrola rotundifolia, using high-throughput sequencing. As expected for non-photosynthetic plants, M. uniflora and H. monotropa have small plastid genomes (46 kb and 35 kb, respectively) lacking genes related to photosynthesis, whereas P. rotundifolia has a larger genome (169 kb) with a gene set similar to other photosynthetic plants. The examined genomes contain an unusually high number of repeats and translocations. Comparative analysis of the expanded set of Ericaceae plastomes suggests that the genes clpP and accD that are present in the plastid genomes of almost all plants have not been lost in this family (as was previously thought) but rather persist in these genomes in unusual forms. Also we found a new gene in P. rotundifolia that emerged as a result of duplication of rps4 gene.

Genome-wide bioinformatics analysis of steroid metabolism-associated genes in Nocardioides simplex VKM Ac-2033D.

Actinobacteria comprise diverse groups of bacteria capable of full degradation, or modification of different steroid compounds. Steroid catabolism has been characterized best for the representatives of suborder Corynebacterineae, such as Mycobacteria, Rhodococcus and Gordonia, with high content of mycolic acids in the cell envelope, while it is poorly understood for other steroid-transforming actinobacteria, such as representatives of Nocardioides genus belonging to suborder Propionibacterineae. Nocardioides simplex VKM Ac-2033D is an important biotechnological strain which is known for its ability to introduce ∆(1)-double bond in various 1(2)-saturated 3-ketosteroids, and perform convertion of 3ß-hydroxy-5-ene steroids to 3-oxo-4-ene steroids, hydrolysis of acetylated steroids, reduction of carbonyl groups at C-17 and C-20 of androstanes and pregnanes, respectively. The strain is also capable of utilizing cholesterol and phytosterol as carbon and energy sources. In this study, a comprehensive bioinformatics genome-wide screening was carried out to predict genes related to steroid metabolism in this organism, their clustering and possible regulation. The predicted operon structure and number of candidate gene copies paralogs have been estimated. Binding sites of steroid catabolism regulators KstR and KstR2 specified for N. simplex VKM Ac-2033D have been calculated de novo. Most of the candidate genes grouped within three main clusters, one of the predicted clusters having no analogs in other actinobacteria studied so far. The results offer a base for further functional studies, expand the understanding of steroid catabolism by actinobacteria, and will contribute to modifying of metabolic pathways in order to generate effective biocatalysts capable of producing valuable bioactive steroids.

Comparative analysis reveals similarities between cultured submandibular salivary gland cells and liver progenitor cells.

Mouse submandibular salivary gland cells and liver progenitor cells from long-term in vitro cultures with a high proliferation potential were side-by-side compared by methods of immunocytochemistry, quantitative real-time PCR, flow cytometry, and transcriptome analysis. The two cell types were found to be similar in expressing cell markers such as EpCAM, CD29, c-Kit, Sca-1, and c-Met. In addition, both cell types expressed cytokeratins 8, 18, and 19, alpha-fetoprotein, and (weakly) albumin. Unlike the liver cells, however, the salivary gland cells in culture showed high-level expression of cytokeratin 14 and CD49f, which was indicative of their origin from salivary gland ducts. Quantitative real-time PCR and deep-sequencing transcriptome analysis revealed similarities in the expression pattern of transcription factors between the two cell types. In this respect, however, the cultured salivary gland cells proved to be closer to exocrine cells of the pancreas than to the liver progenitor cells. Thus, ductal cells of postnatal submandibular salivary glands in culture show phenotypic convergence with progenitor cells of endodermal origin, suggesting that these glands may serve as a potential cell source for cellular therapy of hepatic and pancreatic disorders. The results of this study provide a deeper insight into the molecular features of salivary gland cells and may help optimize procedures for stimulating their differentiation in a specified direction.