COMPSRA: a COMprehensive Platform for Small RNA-Seq data Analysis

A Holistic Analysis of Alzheimer's Disease-Associated lncRNA Communities Reveals Enhanced lncRNA-miRNA-RBP Regulatory Triad Formation Within Functionally Segregated Clusters

Recent studies on the regulatory networks implicated in Alzheimer's disease (AD) evince long non-coding RNAs (lncRNAs) as crucial regulatory players, albeit a poor understanding of the mechanism. Analyzing differential gene expression in the RNA-seq data from the post-mortem AD brain hippocampus, we categorized a list of AD-dysregulated lncRNA transcripts into functionally similar communities based on their k-mer profiles. Using machine-learning-based algorithms, their subcellular localizations were mapped. We further explored the functional relevance of each community through AD-dysregulated miRNA, RNA-binding protein (RBP) interactors, and pathway enrichment analyses. Further investigation of the miRNA-lncRNA and RBP-lncRNA networks from each community revealed the top RBPs, miRNAs, and lncRNAs for each cluster. The experimental validation community yielded ELAVL4 and miR-16-5p as the predominant RBP and miRNA, respectively. Five lncRNAs emerged as the top-ranking candidates from the RBP/miRNA-lncRNA networks. Further analyses of these networks revealed the presence of multiple regulatory triads where the RBP-lncRNA interactions could be augmented by the enhanced miRNA-lncRNA interactions. Our results advance the understanding of the mechanism of lncRNA-mediated AD regulation through their interacting partners and demonstrate how these functionally segregated but overlapping regulatory networks can modulate the disease holistically.

piRNA associates with immune diseases

PIWI-interacting RNA (piRNA) is the most abundant small non-coding RNA in animal cells, typically 26-31 nucleotides in length and it binds with PIWI proteins, a subfamily of Argonaute proteins. Initially discovered in germ cells, piRNA is well known for its role in silencing transposons and maintaining genome integrity. However, piRNA is also present in somatic cells as well as in extracellular vesicles and exosomes. While piRNA has been extensively studied in various diseases, particular cancer, its function in immune diseases remains unclear. In this review, we summarize current research on piRNA in immune diseases. We first introduce the basic characteristics, biogenesis and functions of piRNA. Then, we review the association of piRNA with different types of immune diseases, including autoimmune diseases, immunodeficiency diseases, infectious diseases, and other immune-related diseases. piRNA is considered a promising biomarker for diseases, highlighting the need for further research into its potential mechanisms in disease pathogenesis.

piRNAs in the human retina and retinal pigment epithelium reveal a potential role in intracellular trafficking and oxidative stress

Long considered active only in the germline, the PIWI/piRNA pathway is now known to play a significant role in somatic cells, especially neurons. In this study, piRNAs were profiled in the human retina and retinal pigment epithelium (RPE). Furthermore, RNA immunoprecipitation with HIWI2 (PIWIL4) in ARPE19 cells yielded 261 piRNAs, and the expression of selective piRNAs in donor eyes was assessed by qRT-PCR. Intriguingly, computational analysis revealed complete and partial seed sequence similarity between piR-hsa-26131 and the sensory organ specific miR-183/96/182 cluster. Furthermore, the expression of retina-enriched piR-hsa-26131 was positively correlated with miR-182 in HIWI2-silenced Y79 cells. In addition, the lnc-ZNF169 sequence matched with two miRNAs of the let-7 family, and piRNAs, piR-hsa-11361 and piR-hsa-11360, which could modulate the regulatory network of retinal differentiation. Interestingly, we annotated four enriched motifs among the piRNAs and found that the piRNAs containing CACAATG and CTCATCAKYG motifs were snoRNA-derived piRNAs, which are significantly associated with developmental functions. However, piRNAs consisting of ACCACTANACCAC and AKCACGYTCSC motifs were mainly tRNA-derived fragments linked to stress response and sensory perception. Additionally, co-expression network analysis revealed cell cycle control, intracellular transport and stress response as the important biological functions regulated by piRNAs in the retina. Moreover, loss of piRNAs in HIWI2 knockdown ARPE19 confirmed altered expression of targets implicated in intracellular transport, circadian clock, and retinal degeneration. Moreover, piRNAs were dysregulated under oxidative stress conditions, indicating their potential role in retinal pathology. Therefore, we postulate that piRNAs, miRNAs, and lncRNAs might have a functional interplay during retinal development and functions to regulate retinal homeostasis.

Novel Blood Biomarkers for Response Prediction and Monitoring of Stereotactic Ablative Radiotherapy and Immunotherapy in Metastatic Oligoprogressive Lung Cancer

Up to 80% of patients under immune checkpoint inhibitors (ICI) face resistance. In this context, stereotactic ablative radiotherapy (SABR) can induce an immune or abscopal response. However, its molecular determinants remain unknown. We present early results of a translational study assessing biomarkers of response to combined ICI and SABR (I-SABR) in liquid biopsy from oligoprogressive patients in a prospective observational multicenter study. Cohort A includes metastatic patients in oligoprogression to ICI maintaining the same ICI due to clinical benefit and who receive concomitant SABR. B is a comparative group of oligometastatic patients receiving only SABR. Blood samples are extracted at baseline (T1), after the first (T2) and last (T3) fraction, two months post-SABR (T4) and at further progression (TP). Response is evaluated by iRECIST and defined by the objective response rate (ORR)-complete and partial responses. We assess peripheral blood mononuclear cells (PBMCs), circulating cell-free DNA (cfDNA) and small RNA from extracellular vesicles. Twenty-seven patients could be analyzed (cohort A: n = 19; B: n = 8). Most were males with non-small cell lung cancer and one progressing lesion. With a median follow-up of 6 months, the last ORR was 63% (26% complete and 37% partial response). A decrease in cfDNA from T2 to T3 correlated with a good response. At T2, CD8+PD1+ and CD8+PDL1+ cells were increased in non-responders and responders, respectively. At T2, 27 microRNAs were differentially expressed. These are potential biomarkers of response to I-SABR in oligoprogressive disease.

Circulating miRNAs associate with historical childhood asthma hospitalization in different serum vitamin D groups

BACKGROUND

Vitamin D may help to alleviate asthma exacerbation because of its anti-inflammation effect, but the evidence is inconsistent in childhood asthma. MiRNAs are important mediators in asthma pathogenesis and also excellent non-invasive biomarkers. We hypothesized that circulating miRNAs are associated with asthma exacerbation and modified by vitamin D levels.

METHODS

We sequenced baseline serum miRNAs from 461 participants in the Childhood Asthma Management Program (CAMP). Logistic regression was used to associate miRNA expression with asthma exacerbation through interaction analysis first and then stratified by vitamin D insufficient and sufficient groups. Microarray from lymphoblastoid B-cells (LCLs) treated by vitamin D or sham of 43 subjects in CAMP were used for validation in vitro. The function of miRNAs was associated with gene modules by weighted gene co-expression network analysis (WGCNA).

RESULTS

We identified eleven miRNAs associated with asthma exacerbation with vitamin D effect modification. Of which, five were significant in vitamin D insufficient group and nine were significant in vitamin D sufficient group. Six miRNAs, including hsa-miR-143-3p, hsa-miR-192-5p, hsa-miR-151a-5p, hsa-miR-24-3p, hsa-miR-22-3p and hsa-miR-451a were significantly associated with gene modules of immune-related functions, implying miRNAs may mediate vitamin D effect on asthma exacerbation through immune pathways. In addition, hsa-miR-143-3p and hsa-miR-451a are potential predictors of childhood asthma exacerbation at different vitamin D levels.

CONCLUSIONS

miRNAs are potential mediators of asthma exacerbation and their effects are directly impacted by vitamin D levels.

A novel piwi-interacting RNA associates with type 2-high asthma phenotypes

BACKGROUND

Piwi-interacting RNAs (piRNAs), comprising the largest noncoding RNA group, regulate transcriptional processes. Whether piRNAs are associated with type 2 (T2)-high asthma is unknown.

OBJECTIVE

We sought to investigate the association between piRNAs and T2-high asthma in childhood asthma.

METHODS

We sequenced plasma samples from 462 subjects in the Childhood Asthma Management Program (CAMP) as the discovery cohort and 1165 subjects in the Genetics of Asthma in Costa Rica Study (GACRS) as a replication cohort. Sequencing reads were filtered first, and piRNA reads were annotated and normalized. Linear regression was used for the association analysis of piRNAs and peripheral blood eosinophil count, total serum IgE level, and long-term asthma exacerbation in children with asthma. Mediation analysis was performed to investigate the effect direction. We then ascertained if the circulating piRNAs were present in asthmatic airway epithelial cells in a Gene Expression Omnibus (GEO; www.ncbi.nlm.nih.gov/geo) public data set.

RESULTS

Fifteen piRNAs were significantly associated with eosinophil count in CAMP (P ≤ .05), and 3 were successfully replicated in GACRS. Eleven piRNAs were associated with total IgE in CAMP, and one of these was replicated in GACRS. All 22 significant piRNAs were identified in epithelial cells in vitro, and 6 of these were differentially expressed between subjects with asthma and healthy controls. Fourteen piRNAs were associated with long-term asthma exacerbation, and effect of piRNAs on long-term asthma exacerbation are mediated through eosinophil count and serum IgE level.

CONCLUSION

piRNAs are associated with peripheral blood eosinophils and total serum IgE in childhood asthma and may play important roles in T2-high asthma.

tiny-count: a counting tool for hierarchical classification and quantification of small RNA-seq reads with single-nucleotide precision

Summary

tiny-count is a highly flexible counting tool that allows for hierarchical classification and quantification of small RNA reads from high-throughput sequencing data. Selection rules can be used to filter reads by 5' nucleotide, length, position of alignments in relation to reference features, and by the number of mismatches to reference sequences. tiny-count can quantify reads aligned to a genome or directly to small RNA or transcript sequences. With tiny-count, users can quantify a single class of small RNAs or multiple classes in parallel. tiny-count can resolve distinct classes of small RNAs, for example, piRNAs and siRNAs, produced from the same locus. It can distinguish small RNA variants, such as miRNAs and isomiRs, with single-nucleotide precision. tRNA, rRNA, and other RNA fragments can also be quantified. tiny-count can be run alone or as part of tinyRNA, a workflow that provides a basic all-in-one command line-based solution for small RNA-seq data analysis, with documentation and statistics generated at each step for accurate and reproducible results.

Availability and implementation

tiny-count and other tinyRNA tools are implemented in Python, C++, Cython, and R, and the workflow is coordinated with CWL. tiny-count and tinyRNA are free and open-source software distributed under the GPLv3 license. tiny-count can be installed via Bioconda (https://anaconda.org/bioconda/tiny-count) and both tiny-count and tinyRNA documentation and software downloads are available at https://github.com/MontgomeryLab/tinyRNA. Reference data, including genome and feature information, for certain species can be found at https://www.MontgomeryLab.org.

Omics technologies in allergy and asthma research: An EAACI position paper

Allergic diseases and asthma are heterogenous chronic inflammatory conditions with several distinct complex endotypes. Both environmental and genetic factors can influence the development and progression of allergy. Complex pathogenetic pathways observed in allergic disorders present a challenge in patient management and successful targeted treatment strategies. The increasing availability of high-throughput omics technologies, such as genomics, epigenomics, transcriptomics, proteomics, and metabolomics allows studying biochemical systems and pathophysiological processes underlying allergic responses. Additionally, omics techniques present clinical applicability by functional identification and validation of biomarkers. Therefore, finding molecules or patterns characteristic for distinct immune-inflammatory endotypes, can subsequently influence its development, progression, and treatment. There is a great potential to further increase the effectiveness of single omics approaches by integrating them with other omics, and nonomics data. Systems biology aims to simultaneously and longitudinally understand multiple layers of a complex and multifactorial disease, such as allergy, or asthma by integrating several, separated data sets and generating a complete molecular profile of the condition. With the use of sophisticated biostatistics and machine learning techniques, these approaches provide in-depth insight into individual biological systems and will allow efficient and customized healthcare approaches, called precision medicine. In this EAACI Position Paper, the Task Force "Omics technologies in allergic research" broadly reviewed current advances and applicability of omics techniques in allergic diseases and asthma research, with a focus on methodology and data analysis, aiming to provide researchers (basic and clinical) with a desk reference in the field. The potential of omics strategies in understanding disease pathophysiology and key tools to reach unmet needs in allergy precision medicine, such as successful patients' stratification, accurate disease prognosis, and prediction of treatment efficacy and successful prevention measures are highlighted.

EVAtool: an optimized reads assignment tool for small ncRNA quantification and its application in extracellular vesicle datasets

Extracellular vesicles (EVs) carrying various small non-coding RNAs (sncRNAs) play a vital roles in cell communication and diseases. Correct quantification of multiple sncRNA biotypes simultaneously in EVs is a challenge due to the short reads (<30 bp) could be mapped to multiple sncRNA types. To address this question, we developed an optimized reads assignment algorithm (ORAA) to dynamically map multi-mapping reads to the sncRNA type with a higher proportion. We integrated ORAA with reads processing steps into EVAtool Python-package (http://bioinfo.life.hust.edu.cn/EVAtool) to quantify sncRNAs, especially for sncRNA-seq from EV samples. EVAtool allows users to specify interested sncRNA types in advanced mode or use default seven sncRNAs (microRNA, small nucleolar RNA, PIWI-interacting RNAs, small nuclear RNA, ribosomal RNA, transfer RNA and Y RNA). To prove the utilities of EVAtool, we quantified the sncRNA expression profiles for 200 samples from cognitive decline and multiple sclerosis. We found that more than 20% of short reads on average were mapped to multiple sncRNA biotypes in multiple sclerosis. In cognitive decline, the proportion of Y RNA is significantly higher than other sncRNA types. EVAtool is a flexible and extensible tool that would benefit to mine potential biomarkers and functional molecules in EVs.

Blood miRNAs Are Linked to Frequent Asthma Exacerbations in Childhood Asthma and Adult COPD

MicroRNAs have been independently associated with asthma and COPD; however, it is unclear if microRNA associations will overlap when evaluating retrospective acute exacerbations. Objective: We hypothesized that peripheral blood microRNAs would be associated with retrospective acute asthma exacerbations in a pediatric asthma cohort and that such associations may also be relevant to acute COPD exacerbations. Methods: We conducted small-RNA sequencing on 374 whole-blood samples from children with asthma ages 6-14 years who participated in the Genetics of Asthma in Costa Rica Study (GACRS) and 450 current and former adult smokers with and without COPD who participated in the COPDGene study. Measurements and Main Results: After QC, we had 351 samples and 649 microRNAs for Differential Expression (DE) analysis between the frequent (n = 183) and no or infrequent exacerbation (n = 168) groups in GACRS. Fifteen upregulated miRs had odds ratios (OR) between 1.22 and 1.59 for a doubling of miR counts, while five downregulated miRs had ORs between 0.57 and 0.8. These were assessed for generalization in COPDGene, where three of the upregulated miRs (miR-532-3p, miR-296-5p, and miR-766-3p) and two of the downregulated miRs (miR-7-5p and miR-451b) replicated. Pathway enrichment analysis showed MAPK and PI3K-Akt signaling pathways were strongly enriched for target genes of DE miRNAs and miRNAs generalizing to COPD exacerbations, as well as infection response pathways to various pathogens. Conclusion: miRs (451b; 7-5p; 532-3p; 296-5p and 766-3p) associated with both childhood asthma and adult COPD exacerbations may play a vital role in airflow obstruction and exacerbations and point to shared genomic regulatory machinery underlying exacerbations in both diseases.

gga-miRNOME, a microRNA-sequencing dataset from chick embryonic tissues

MicroRNAs (miRNAs) are small non-coding RNA molecules, with sizes ranging from 18 to 25 nucleotides, which are key players in gene expression regulation. These molecules play an important role in fine-tuning early vertebrate embryo development. However, there are scarce publicly available miRNA datasets from non-mammal embryos, such as the chicken (Gallus gallus), which is a classical model system to study vertebrate embryogenesis. Here, we performed microRNA-sequencing to characterize the early stages of trunk and limb development in the chick embryo. For this, we profiled three chick embryonic tissues, namely, Undetermined Presomitic Mesoderm (PSM_U), Determined Presomitic Mesoderm (PSM_D) and Forelimb Distal Cyclic Domain (DCD). We identified 926 known miRNAs, and 1,141 novel candidate miRNAs, which nearly duplicates the number of Gallus gallus entries in the miRBase database. These data will greatly benefit the avian research community, particularly by highlighting new miRNAs potentially involved in the regulation of early vertebrate embryo development, that can be prioritized for further experimental testing.

Measurement Biases Distort Cell-Free DNA Fragmentation Profiles and Define the Sensitivity of Metagenomic Cell-Free DNA Sequencing Assays

BACKGROUND

Metagenomic sequencing of microbial cell-free DNA (cfDNA) in blood and urine is increasingly used as a tool for unbiased infection screening. The sensitivity of metagenomic cfDNA sequencing assays is determined by the efficiency by which the assay recovers microbial cfDNA vs host-specific cfDNA. We hypothesized that the choice of methods used for DNA isolation, DNA sequencing library preparation, and sequencing would affect the sensitivity of metagenomic cfDNA sequencing.

METHODS

We characterized the fragment length biases inherent to select DNA isolation and library preparation procedures and developed a model to correct for these biases. We analyzed 305 cfDNA sequencing data sets, including publicly available data sets and 124 newly generated data sets, to evaluate the dependence of the sensitivity of metagenomic cfDNA sequencing on pre-analytical variables.

RESULTS

Length bias correction of fragment length distributions measured from different experimental procedures revealed the ultrashort (<100 bp) nature of microbial-, mitochondrial-, and host-specific urinary cfDNA. The sensitivity of metagenomic sequencing assays to detect the clinically reported microorganism differed by more than 5-fold depending on the combination of DNA isolation and library preparation used.

CONCLUSIONS

Substantial gains in the sensitivity of microbial and other short fragment recovery can be achieved by easy-to-implement changes in the sample preparation protocol, which highlights the need for standardization in the liquid biopsy field.

Transposable elements and Piwi‑interacting RNAs in hemato‑oncology with a focus on myelodysplastic syndrome (Review)

Our current understanding of hematopoietic stem cell differentiation and the abnormalities that lead to leukemogenesis originates from the accumulation of knowledge regarding protein‑coding genes. However, the possible impact of transposable element (TE) mobilization and the expression of P‑element‑induced WImpy testis‑interacting RNAs (piRNAs) on leukemogenesis has been beyond the scope of scientific interest to date. The expression profiles of these molecules and their importance for human health have only been characterized recently due to the rapid progress of high‑throughput sequencing technology development. In the present review, current knowledge on the expression profile and function of TEs and piRNAs was summarized, with specific focus on their reported involvement in leukemogenesis and pathogenesis of myelodysplastic syndrome.

Circulating MicroRNA: Incident Asthma Prediction and Vitamin D Effect Modification

Of children with recurrent wheezing in early childhood, approximately half go on to develop asthma. MicroRNAs have been described as excellent non-invasive biomarkers due to their prognostic utility. We hypothesized that circulating microRNAs can predict incident asthma and that that prediction might be modified by vitamin D. We selected 75 participants with recurrent wheezing at 3 years old from the Vitamin D Antenatal Asthma Reduction Trial (VDAART). Plasma samples were collected at age 3 and sequenced for small RNA-Seq. The read counts were normalized and filtered by depth and coverage. Logistic regression was employed to associate miRNAs at age 3 with asthma status at age 5. While the overall effect of miRNA on asthma occurrence was weak, we identified 38 miRNAs with a significant interaction effect with vitamin D and 32 miRNAs with a significant main effect in the high vitamin D treatment group in VDAART. We validated the VDAART results in Project Viva for both the main effect and interaction effect. Meta-analysis was performed on both cohorts to obtain the combined effect and a logistic regression model was used to predict incident asthma at age 7 in Project Viva. Of the 23 overlapped miRNAs in the stratified and interaction analysis above, 9 miRNAs were replicated in Project Viva with strong effect size and remained in the meta-analysis of the two populations. The target genes of the 9 miRNAs were enriched for asthma-related Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways. Using logistic regression, microRNA hsa-miR-574-5p had a good prognostic ability for incident asthma prognosis with an area under the receiver operating characteristic (AUROC) of 0.83. In conclusion, miRNAs appear to be good biomarkers of incident asthma, but only when vitamin D level is considered.