Comparison of alternative approaches for analysing multi-level RNA-seq data

RNA sequencing (RNA-seq) is widely used for RNA quantification in the environmental, biological and medical sciences. It enables the description of genome-wide patterns of expression and the identification of regulatory interactions and networks. The aim of RNA-seq data analyses is to achieve rigorous quantification of genes/transcripts to allow a reliable prediction of differential expression (DE), despite variation in levels of noise and inherent biases in sequencing data. This can be especially challenging for datasets in which gene expression differences are subtle, as in the behavioural transcriptomics test dataset from D. melanogaster that we used here. We investigated the power of existing approaches for quality checking mRNA-seq data and explored additional, quantitative quality checks. To accommodate nested, multi-level experimental designs, we incorporated sample layout into our analyses. We employed a subsampling without replacement-based normalization and an identification of DE that accounted for the hierarchy and amplitude of effect sizes within samples, then evaluated the resulting differential expression call in comparison to existing approaches. In a final step to test for broader applicability, we applied our approaches to a published set of H. sapiens mRNA-seq samples, The dataset-tailored methods improved sample comparability and delivered a robust prediction of subtle gene expression changes. The proposed approaches have the potential to improve key steps in the analysis of RNA-seq data by incorporating the structure and characteristics of biological experiments

Small RNA analysis in Sindbis virus infected human HEK293 cells

In contrast to the defence mechanism of RNA interference (RNAi) in plants and invertebrates, its role in the innate response to virus infection of mammals is a matter of debate. Since RNAi has a well-established role in controlling infection of the alphavirus Sindbis virus (SINV) in insects, we have used this virus to investigate the role of RNAi in SINV infection of human cells

Small RNA populations revealed by blocking rRNA fragments in Drosophila melanogaster reproductive tissues

RNA interference (RNAi) is a complex and highly conserved regulatory mechanism mediated via small RNAs (sRNAs). Recent technical advances in high throughput sequencing have enabled an increasingly detailed analysis of sRNA abundances and profiles in specific body parts and tissues. This enables investigations of the localized roles of microRNAs (miRNAs) and small interfering RNAs (siRNAs). However, variation in the proportions of non-coding RNAs in the samples being compared can hinder these analyses. Specific tissues may vary significantly in the proportions of fragments of longer non-coding RNAs (such as ribosomal RNA or transfer RNA) present, potentially reflecting tissue-specific differences in biological functions. For example, in Drosophila, some tissues contain a highly abundant 30nt rRNA fragment (the 2S rRNA) as well as abundant 5’ and 3’ terminal rRNA fragments. These can pose difficulties for the construction of sRNA libraries as they can swamp the sequencing space and obscure sRNA abundances. Here we addressed this problem and present a modified “rRNA blocking” protocol for the construction of high-definition (HD) adapter sRNA libraries, in D. melanogaster reproductive tissues. The results showed that 2S rRNAs targeted by blocking oligos were reduced from >80% to < 0.01% total reads. In addition, the use of multiple rRNA blocking oligos to bind the most abundant rRNA fragments allowed us to reveal the underlying sRNA populations at increased resolution. Side-by-side comparisons of sequencing libraries of blocked and non-blocked samples revealed that rRNA blocking did not change the miRNA populations present, but instead enhanced their abundances. We suggest that this rRNA blocking procedure offers the potential to improve the in-depth analysis of differentially expressed sRNAs within and across different tissues

Viral transduction of primary human lymphoma B cells reveals mechanisms of NOTCH-mediated immune escape

Hotspot mutations in the PEST-domain of NOTCH1 and NOTCH2 are recurrently identified in B cell malignancies. To address how NOTCH-mutations contribute to a dismal prognosis, we have generated isogenic primary human tumor cells from patients with Chronic Lymphocytic Leukemia (CLL) and Mantle Cell Lymphoma (MCL), differing only in their expression of the intracellular domain (ICD) of NOTCH1 or NOTCH2. Our data demonstrate that both NOTCH-paralogs facilitate immune-escape of malignant B cells by up-regulating PD-L1, partly dependent on autocrine interferon-? signaling. In addition, NOTCH-activation causes silencing of the entire HLA-class II locus via epigenetic regulation of the transcriptional co-activator CIITA. Notably, while NOTCH1 and NOTCH2 govern similar transcriptional programs, disease-specific differences in their expression levels can favor paralog-specific selection. Importantly, NOTCH-ICD also strongly down-regulates the expression of CD19, possibly limiting the effectiveness of immune-therapies. These NOTCH-mediated immune escape mechanisms are associated with the expansion of exhausted CD8+ T cells in vivo.© 2022. The Author(s)

Small RNA Profile in Moso Bamboo Root and Leaf Obtained by High Definition Adapters

Moso bamboo (Phyllostachy heterocycla cv. pubescens L.) is an economically important fast-growing tree. In order to gain better understanding of gene expression regulation in this important species we used next generation sequencing to profile small RNAs in leaf and roots of young seedlings. Since standard kits to produce cDNA of small RNAs are biased for certain small RNAs, we used High Definition adapters that reduce ligation bias. We identified and experimentally validated five new microRNAs and a few other small non-coding RNAs that were not microRNAs. The biological implication of microRNA expression levels and targets of microRNAs are discussed

Viral transduction of primary human lymphoma B cells reveals mechanisms of NOTCH-mediated immune escape

Hotspot mutations in the PEST-domain of NOTCH1 and NOTCH2 are recurrently identified in B cell malignancies. To address how NOTCH-mutations contribute to a dismal prognosis, we have generated isogenic primary human tumor cells from patients with Chronic Lymphocytic Leukemia (CLL) and Mantle Cell Lymphoma (MCL), differing only in their expression of the intracellular domain (ICD) of NOTCH1 or NOTCH2. Our data demonstrate that both NOTCH-paralogs facilitate immune-escape of malignant B cells by up-regulating PD-L1, partly dependent on autocrine interferon-γ signaling. In addition, NOTCH-activation causes silencing of the entire HLA-class II locus via epigenetic regulation of the transcriptional co-activator CIITA. Notably, while NOTCH1 and NOTCH2 govern similar transcriptional programs, disease-specific differences in their expression levels can favor paralog-specific selection. Importantly, NOTCH-ICD also strongly down-regulates the expression of CD19, possibly limiting the effectiveness of immune-therapies. These NOTCH-mediated immune escape mechanisms are associated with the expansion of exhausted CD8+ T cells in vivo

Mitochondrial complex I activity in microglia sustains neuroinflammation

Sustained smouldering, or low-grade activation, of myeloid cells is a common hallmark of several chronic neurological diseases, including multiple sclerosis1. Distinct metabolic and mitochondrial features guide the activation and the diverse functional states of myeloid cells2. However, how these metabolic features act to perpetuate inflammation of the central nervous system is unclear. Here, using a multiomics approach, we identify a molecular signature that sustains the activation of microglia through mitochondrial complex I activity driving reverse electron transport and the production of reactive oxygen species. Mechanistically, blocking complex I in pro-inflammatory microglia protects the central nervous system against neurotoxic damage and improves functional outcomes in an animal disease model in vivo. Complex I activity in microglia is a potential therapeutic target to foster neuroprotection in chronic inflammatory disorders of the central nervous system3

Divergence in transcriptional and regulatory responses to mating in male and female fruitflies

Mating induces extensive physiological, biochemical and behavioural changes in female animals of many taxa. In contrast, the overall phenotypic and transcriptomic consequences of mating for males, hence how they might differ from those of females, are poorly described. Post mating responses in each sex are rapidly initiated, predicting the existence of regulatory mechanisms in addition to transcriptional responses involving de novo gene expression. That post mating responses appear different for each sex also predicts that the genome-wide signatures of mating should show evidence of sex-specific specialisation. In this study, we used high resolution RNA sequencing to provide the first direct comparisons of the transcriptomic responses of male and female Drosophila to mating, and the first comparison of mating-responsive miRNAs in both sexes in any species. As predicted, the results revealed the existence of sex- and body part-specific mRNA and miRNA expression profiles. More genes were differentially expressed in the female head-thorax than the abdomen following mating, whereas the opposite was true in males. Indeed, the transcriptional profile of male head-thorax tissue was largely unaffected by mating, and no differentially expressed genes were detected at the most stringent significance threshold. A subset of ribosomal genes in females were differentially expressed in both body parts, but in opposite directions, consistent with the existence of body part-specific resource allocation switching. Novel, mating-responsive miRNAs in each sex were also identified, and a miRNA-mRNA interactions analysis revealed putative targets among mating-responsive genes. We show that the structure of genome-wide responses by each sex to mating is strongly divergent, and provide new insights into how shared genomes can achieve characteristic distinctiveness

Cloning and Characterization of Maize miRNAs Involved in Responses to Nitrogen Deficiency

Although recent studies indicated that miRNAs regulate plant adaptive responses to nutrient deprivation, the functional significance of miRNAs in adaptive responses to nitrogen (N) limitation remains to be explored. To elucidate the molecular biology underlying N sensing/signaling in maize, we constructed four small RNA libraries and one degradome from maize seedlings exposed to N deficiency. We discovered a total of 99 absolutely new loci belonging to 47 miRNA families by small RNA deep sequencing and degradome sequencing, as well as 9 new loci were the paralogs of previously reported miR169, miR171, and miR398, significantly expanding the reported 150 high confidence genes within 26 miRNA families in maize. Bioinformatic and subsequent small RNA northern blot analysis identified eight miRNA families (five conserved and three newly identified) differentially expressed under the N-deficient condition. Predicted and degradome-validated targets of the newly identified miRNAs suggest their involvement in a broad range of cellular responses and metabolic processes. Because maize is not only an important crop but is also a genetic model for basic biological research, our research contributes to the understanding of the regulatory roles of miRNAs in plant adaption to N-deficiency stress

Identification and characterization of microRNAs expressed in the African malaria vector Anopheles funestus life stages using high throughput sequencing

Background: Over the past several years, thousands of microRNAs (miRNAs) have been identified in the genomes of various insects through cloning and sequencing or even by computational prediction. However, the number of miRNAs identified in anopheline species is low and little is known about their role. The mosquito Anopheles funestus is one of the dominant malaria vectors in Africa, which infects and kills millions of people every year. Therefore, small RNA molecules isolated from the four life stages (eggs, larvae, pupae and unfed adult females) of An. funestus were sequenced using next generation sequencing technology. Results: High throughput sequencing of four replicates in combination with computational analysis identified 107 mature miRNA sequences expressed in the An. funestus mosquito. These include 20 novel miRNAs without sequence identity in any organism and eight miRNAs not previously reported in the Anopheles genus but are known in non-anopheles mosquitoes. Finally, the changes in the expression of miRNAs during the mosquito development were determined and the analysis showed that many miRNAs have stage-specific expression, and are co-transcribed and co-regulated during development. Conclusions: This study presents the first direct experimental evidence of miRNAs in An. funestus and the first profiling study of miRNA associated with the maturation in this mosquito. Overall, the results indicate that miRNAs play important roles during the growth and development. Silencing such molecules in a specific life stage could decrease the vector population and therefore interrupt malaria transmission.IS