tsRNAsearch: A pipeline for the identification of tRNA and ncRNA fragments from small RNA-sequencing data

MOTIVATION

tRNAs were originally considered uni-functional RNA molecules involved in the delivery of amino acids to growing peptide chains on the ribosome. More recently, the liberation of tRNA fragments from tRNAs via specific enzyme cleavage has been characterized. Detection of tRNA fragments in sequencing data is difficult due to tRNA sequence redundancy and the short length of both tRNAs and their fragments.

RESULTS

Here we introduce tsRNAsearch, a Nextflow pipeline for the identification of differentially abundant tRNA fragments and other non-coding RNAs from small RNA-sequencing data. tsRNAsearch is intended for use when comparing two groups of datasets, such as control and treatment groups. tsRNAsearch comparatively searches for tRNAs and ncRNAs with irregular read distribution profiles (a proxy for RNA cleavage) using a combined score made up of four novel methods and a differential expression analysis, and reports the top ranked results in simple PDF and TEXT files. In this study, we used publicly available small RNA-seq data to replicate the identification of tsRNAs from chronic hepatitis-infected liver tissue data. In addition, we applied tsRNAsearch to pancreatic ductal adenocarcinoma (PDAC) and matched healthy pancreatic tissue small RNA-sequencing data. Our results support the identification of miR135b from the original study as a potential biomarker of PDAC and identify other potentially stronger miRNA biomarkers of PDAC.

AVAILABILITY

https://github.com/GiantSpaceRobot/tsRNAsearch.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Polymorphic centromere locations in the pathogenic yeast Candida parapsilosis

Centromeres pose an evolutionary paradox: strongly conserved in function but rapidly changing in sequence and structure. However, in the absence of damage, centromere locations are usually conserved within a species. We report here that isolates of the pathogenic yeast species Candida parapsilosis show within-species polymorphism for the location of centromeres on two of its eight chromosomes. Its old centromeres have an inverted-repeat (IR) structure, whereas its new centromeres have no obvious structural features but are located within 30 kb of the old site. Centromeres can therefore move naturally from one chromosomal site to another, apparently spontaneously and in the absence of any significant changes in DNA sequence. Our observations are consistent with a model in which all centromeres are genetically determined, such as by the presence of short or long IRs or by the ability to form cruciforms. We also find that centromeres have been hotspots for genomic rearrangements in the C. parapsilosis clade.

TPP riboswitch-dependent regulation of an ancient thiamin transporter in Candida

Riboswitches are non-coding RNA molecules that regulate gene expression by binding to specific ligands. They are primarily found in bacteria. However, one riboswitch type, the thiamin pyrophosphate (TPP) riboswitch, has also been described in some plants, marine protists and fungi. We find that riboswitches are widespread in the budding yeasts (Saccharomycotina), and they are most common in homologs of DUR31, originally described as a spermidine transporter. We show that DUR31 (an ortholog of N. crassa gene NCU01977) encodes a thiamin transporter in Candida species. Using an RFP/riboswitch expression system, we show that the functional elements of the riboswitch are contained within the native intron of DUR31 from Candida parapsilosis, and that the riboswitch regulates splicing in a thiamin-dependent manner when RFP is constitutively expressed. The DUR31 gene has been lost from Saccharomyces, and may have been displaced by an alternative thiamin transporter. TPP riboswitches are also present in other putative transporters in yeasts and filamentous fungi. However, they are rare in thiamin biosynthesis genes THI4 and THI5 in the Saccharomycotina, and have been lost from all genes in the sequenced species in the family Saccharomycetaceae, including S. cerevisiae.

Thiamin, or Vitamin B1, is an essential requirement in all living organisms because it is a co-factor for many enzymes in metabolism. Unlike animals, many yeasts can synthesize thiamin, or they can import it from the environment. Expression of thiamin biosynthesis genes and of thiamin transporters is strictly regulated in response to the presence of thiamin. In many filamentous fungi, expression of thiamin biosynthesis genes is regulated by TPP riboswitches, RNA regulatory elements that are located within messenger RNA. TPP riboswitches are rare in yeasts. However, we find that TPP riboswitches are conserved in an ancient thiamin transporter, found in filamentous fungi, yeasts and other related organisms. There appears to be a high turnover of thiamin transporters in fungi, and there has been a gradual loss of TPP riboswitches in yeasts.

Identification of Non-Coding RNAs in the Candida parapsilosis Species Group

The Candida CTG clade is a monophyletic group of fungal species that translates CTG as serine, and includes the pathogens Candida albicans and Candida parapsilosis. Research has typically focused on identifying protein-coding genes in these species. Here, we use bioinformatic and experimental approaches to annotate known classes of non-coding RNAs in three CTG-clade species, Candida parapsilosis, Candida orthopsilosis and Lodderomyces elongisporus. We also update the annotation of ncRNAs in the C. albicans genome. The majority of ncRNAs identified were snoRNAs. Approximately 50% of snoRNAs (including most of the C/D box class) are encoded in introns. Most are within mono- and polycistronic transcripts with no protein coding potential. Five polycistronic clusters of snoRNAs are highly conserved in fungi. In polycistronic regions, splicing occurs via the classical pathway, as well as by nested and recursive splicing. We identified spliceosomal small nuclear RNAs, the telomerase RNA component, signal recognition particle, RNase P RNA component and the related RNase MRP RNA component in all three genomes. Stem loop IV of the U2 spliceosomal RNA and the associated binding proteins were lost from the ancestor of C. parapsilosis and C. orthopsilosis, following the divergence from L. elongisporus. The RNA component of the MRP is longer in C. parapsilosis, C. orthopsilosis and L. elongisporus than in S. cerevisiae, but is substantially shorter than in C. albicans.

Effect of stopper processing conditions on moisture content and ramifications for lyophilized products: comparison of "low" and "high" moisture uptake stoppers

The purpose of this study was to evaluate the effect of processing and storage on the moisture content of two commercially available, 13-mm lyophilization stoppers designated as low moisture (LM) and high moisture (HM) uptake stoppers. The stopper moisture studies included the effect of steam sterilization time, drying time and temperature, equilibrium moisture content, lyophilization and moisture transfer from stopper to a model-lactose lyophilized cake. Results indicated that both stoppers absorbed significant amounts of moisture during sterilization and that the HM stopper absorbed significantly more water than the LM stopper. LM and HM stoppers required approximately 2 and 8 h drying at 105 degrees C, respectively, to achieve final moisture content of not more than 0.5 mg/stopper. Following drying, stopper moisture levels equilibrated rapidly to ambient storage conditions. The apparent equilibrium moisture level was approximately 7 times higher in the HM versus LM stopper. Freeze-drying had minimal effect on the moisture content of dried stoppers. Finally, moisture transfer from the stopper to the lyophilized product is dependent on the initial stopper water content and storage temperature. To better quantify the ramifications of stopper moisture, projections of moisture uptake over the shelf life of a drug product were calculated based on the product-contact surface area of stoppers. Attention to stopper storage conditions prior to use, in addition to processing steps, are necessary to minimize stability issues especially in low-fill, mass lyophilized products.