Search

Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development. Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in semen in other animals. Intriguingly we find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination. We conclude that small RNAs play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of origin-effects on gene expression and reproductive physiology.

Background High-throughput single-cell RNA-seq (scRNA-seq) is a powerful technology for studying gene expression variability in single cells; however, standard analysis approaches only consider the overall expression of each gene, masking additional heterogeneity that could exist through cell type-specific expression of alternative mRNA transcripts. Results Here we show that differential transcript usage (DTU) can be readily detected in data-sets generated from commonly used polyA-captured nanodroplet scRNA-seq technology. Our computational pipeline, Sierra, detects and quantifies polyadenylation sites in scRNA-seq data-sets, which are utilised to evaluate DTU between single-cell populations. We validate our approach by comparing cardiac cell populations derived from scRNA-seq to bulk RNA-seq of matched populations obtained through fluorescent activated cell sorting (FACS), demonstrating that we detect a significant overlap in cell-type DTU between the congruous data-sets. We further illustrate the utility of our method by detecting alternative transcript usage in human peripheral blood mononuclear cells (PBMCs), 3’UTR shortening in activated and proliferating cardiac fibroblasts from injured mouse hearts and finally by building an initial atlas of cell type-specific transcript usage across 12 mouse tissues. Conclusions We anticipate that Sierra will enable new avenues of transcriptional complexity and regulation to be explored in single-cell transcriptomic experiments. Sierra is available at https://github.com/VCCRI/Sierra .

Concerted examination of multiple collections of single cell RNA-Seq (scRNA-Seq) data promises further biological insights that cannot be uncovered with individual datasets. However, such integrative analyses are challenging and require sophisticated methodologies. To enable effective interrogation of multiple scRNA-Seq datasets, we have developed a novel algorithm, named scMerge, that removes unwanted variation by combining stably expressed genes and utilizing pseudo-replicates across datasets. Analysis of large collections of publicly available datasets demonstrates that scMerge performs well in multiple scenarios and enhances biological discovery, including inferring cell developmental trajectories.