Discovering microRNAs from deep sequencing data using miRDeep

Climate challenges for fish larvae: Interactive multi-stressor effects impair acclimation potential of Atlantic herring larvae

Fish early life stages are particularly vulnerable and heavily affected by changing environmental factors. The interactive effects of multiple climate change-related stressors on fish larvae remain, however, largely underexplored. As rising temperatures can increase the abundance and virulence of bacteria, we investigated the combination of a spring heat wave and bacterial exposure on the development of Atlantic herring larvae (Clupea harengus). Eggs and larvae of Western Baltic Spring-spawners were reared at a normal and high temperature ramp and exposed to Vibrio alginolyticus and V. anguillarum, respectively. Subsequently, mRNA and miRNA transcriptomes, microbiota composition, growth and survival were assessed. Both high temperature and V. alginolyticus exposure induced a major downregulation of gene expression likely impeding larval cell proliferation. In contrast, interactive effects of elevated temperature and V. alginolyticus resulted in minimal gene expression changes, indicating an impaired plastic response, which may cause cellular damage reducing survival in later larval stages. The heat wave alone or in combination with V. alginolyticus induced a notable shift in miRNA expression leading to the down- but also upregulation of predicted target genes. Moreover, both increased temperature and the Vibrio exposures significantly altered the larval microbiota composition, with warming reducing microbial richness and diversity. The outcomes of this study highlight the high sensitivity of herring early life stages towards multiple climate change-related stressors. Our results indicate that interactive effects of rapidly changing environmental factors may exceed the larval stress threshold impairing essential acclimation responses, which may contribute to the ongoing recruitment decline of Western Baltic Spring-Spawning herring.

Genome-Wide Transcriptional Response of Avocado to Fusarium sp. Infection

The avocado crop is relevant for its economic importance and because of its unique evolutionary history. However, there is a lack of information regarding the molecular processes during the defense response against fungal pathogens. Therefore, using a genome-wide approach in this work, we investigated the transcriptional response of the Mexican horticultural race of avocado (Persea americana var. drymifolia), including miRNAs profile and their possible targets. For that, we established an avocado-Fusarium hydroponic pathosystem and studied the response for 21 days. To guarantee robustness in the analysis, first, we improved the avocado genome assembly available for this variety, resulting in 822.49 Mbp in length with 36,200 gene models. Then, using an RNA-seq approach, we identified 13,778 genes differentially expressed in response to the Fusarium infection. According to their expression profile across time, these genes can be clustered into six groups, each associated with specific biological processes. Regarding non-coding RNAs, 8 of the 57 mature miRNAs identified in the avocado genome are responsive to infection caused by Fusarium, and the analysis revealed a total of 569 target genes whose transcript could be post-transcriptionally regulated. This study represents the first research in avocados to comprehensively explore the role of miRNAs in orchestrating defense responses against Fusarium spp. Also, this work provides valuable data about the genes involved in the intricate response of the avocado during fungal infection.

Mechanisms underlying TRPV4-mediated regulation of miR-146a expression

Persistent inflammation is a major contributor in the development of various inflammatory diseases like atherosclerosis. Our study investigates how transient receptor potential vanilloid 4 (TRPV4), a mechanosensitive ion channel, interacts with microRNA-146a (miR-146a), within the context of inflammation and atherosclerosis. Micro-RNAs play a critical role in controlling gene expression, and miR-146a is notable for its anti-inflammatory actions. TRPV4 is activated by diverse soluble and mechanical stimuli, and often associated with inflammatory responses in various diseases. Here, we find that TRPV4 negatively regulates miR-146a expression in macrophages, especially following stimulation by lipopolysaccharides or alterations in matrix stiffness. We show that in atherosclerosis, a condition characterized by matrix stiffening, TRPV4 decreases miR-146a expression in aortic tissue macrophages. We find that TRPV4's impact on miR-146a is independent of activation of NFκB, Stat1, P38, and AKT, but is rather mediated through a mechanism involving histone deacetylation instead of DNA methylation at the miR-146a promoter site. Furthermore, we show that N-terminal residues 1 to 130 in TRPV4 is essential in suppression of miR-146a expression in LPS-stimulated macrophages. Altogether, this study identifies a regulatory mechanism of miR-146a expression by TRPV4 which may open new potential therapeutic strategies for managing inflammatory diseases.

A glimpse into the world of microRNAs and their putative roles in hard ticks

Ticks are important blood feeding ectoparasites that transmit pathogens to wildlife, domestic animals, and humans. Hard ticks can feed for several days to weeks, nevertheless they often go undetected. This phenomenon can be explained by a tick's ability to release analgesics, immunosuppressives, anticoagulants, and vasodilators within their saliva. Several studies have identified extracellular vesicles (EVs) as carriers of some of these effector molecules. Further, EVs, and their contents, enhance pathogen transmission, modulate immune responses, and delay wound healing. EVs are double lipid-membrane vesicles that transport intracellular cargo, including microRNAs (miRNAs) to recipient cells. miRNAs are involved in regulating gene expression post-transcriptionally. Interestingly, tick-derived miRNAs have been shown to enhance pathogen transmission and affect vital biological processes such as oviposition, blood digestion, and molting. miRNAs have been found within tick salivary EVs. This review focuses on current knowledge of miRNA loading into EVs and homologies reported in ticks. We also describe findings in tick miRNA profiles, including miRNAs packed within tick salivary EVs. Although no functional studies have been done to investigate the role of EV-derived miRNAs in tick feeding, we discuss the functional characterization of miRNAs in tick biology and pathogen transmission. Lastly, we propose the possible uses of tick miRNAs to develop management tools for tick control and to prevent pathogen transmission. The identification and functional characterization of conserved and tick-specific salivary miRNAs targeting important molecular and immunological pathways within the host could lead to the discovery of new therapeutics for the treatment of tick-borne and non-tick-borne human diseases.

HASTY-mediated miRNA dynamics modulate nitrogen starvation-induced leaf senescence in Arabidopsis

Nitrogen (N) deficiency responses are essential for plant survival and reproduction. Here, via an expression genome-wide association study (eGWAS), we reveal a mechanism that regulates microRNA (miRNA) dynamics necessary for N deficiency responses in Arabidopsis. Differential expression levels of three NAC transcription factor (TF) genes involved in leaf N deficiency responses among Arabidopsis accessions are most significantly associated with polymorphisms in HASTY (HST), which encodes an importin/exportin family protein responsible for the generation of mature miRNAs. HST acts as a negative regulator of N deficiency-induced leaf senescence, and the disruption and overexpression of HST differently modifies miRNA dynamics in response to N deficiency, altering levels of miRNAs targeting transcripts. Interestingly, N deficiency prevents the interaction of HST with HST-interacting proteins, DCL1 and RAN1, and some miRNAs. This suggests that HST-mediated regulation of miRNA dynamics collectively controls regulations mediated by multiple N deficiency response-associated NAC TFs, thereby being central to the N deficiency response network.

microRNA profile of endometrial cancer from Indian patients-identification of potential biomarkers for prognosis

Endometrial cancer is one of the major cancers in women throughout the world. If diagnosed early, these cancers are treatable and the prognosis is usually good. However, one major problem in treating endometrial cancer is accurate diagnosis and staging. Till date, the choice method for diagnosis and staging is histopathology. Although there are few molecular markers identified, they are not always sufficient in making accurate diagnosis and deciding on therapeutic strategy. As a result, very often patients are under treated or over treated. In this study, our group has profiled microRNAs from Indian patients using NGS-based approach. We have identified 212 differentially expressed microRNAs in endometrial cancer. Among these, there are 17 novel miRNAs. Since this data represents only Indian cohort and also lacks survival data, validation across other populations is necessary before being considered as biomarkers. As one approach towards this, these microRNAs have also been compared to data from TCGA, which represent other populations and also correlated to relevance in overall survival. Using in-silico approaches, mRNA targets of the miRNAs have been predicted. After comparing with TCGA, we have identified 16 miRNA-mRNA pairs which could be potential prognostic biomarkers for endometrial cancer. This is the first miRNA profiling report from Indian cohort and one of the very few studies which have identified potential biomarkers of prognosis in endometrial cancer.

Identification and Prediction of Differentially Expressed MicroRNAs Associated with Detoxification Pathways in Larvae of Spodoptera frugiperda

Spodoptera frugiperda poses a severe threat to crops, causing substantial economic losses. The increased use of chemical pesticides has led to resistance in S. frugiperda populations. Micro ribonucleic acids (MicroRNAs or miRNAs) are pivotal in insect growth and development. This study aims to identify miRNAs across different developmental stages of S. frugiperda to explore differential expression and predict target gene functions. High-throughput sequencing of miRNAs was conducted on eggs, 3rd instar larvae, pupae, and adults. Bioinformatics analyses identified differentially expressed miRNAs specifically in larvae, with candidate miRNAs screened to predict target genes, particularly those involved in detoxification pathways. A total of 184 known miRNAs and 209 novel miRNAs were identified across stages. Comparative analysis revealed 54, 15, and 18 miRNAs differentially expressed in larvae, compared to egg, pupa, and adult stages, respectively. Eight miRNAs showed significant differential expression across stages, validated by quantitative reverse transcription PCR (qRT-PCR). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses predicted target genes' functions, identifying eight differentially expressed miRNAs targeting 10 gene families associated with detoxification metabolism, including P450s, glutathione S-transferase (GSTs), ATP-binding cassette (ABC) transporters, and sodium channels. These findings elucidate the species-specific miRNA profiles and regulatory mechanisms of detoxification-related genes in S. frugiperda larvae, offering insights and strategies for effectively managing this pest.

Primer exchange reaction assisted CRISPR/Cas9 cleavage for detection of dual microRNAs with electrochemistry method

An electrochemical sensor assisted by primer exchange reaction (PER) and CRISPR/Cas9 system (PER-CRISPR/Cas9-E) was established for the sensitive detection of dual microRNAs (miRNAs). Two PER hairpin (HP) were designed to produce a lot of extended PER products, which could hybridize with two kinds of hairpin probes modified on the electrode and initiate the cleavage of two CRISPR/Cas9 systems guided by single guide RNAs (sgRNAs) with different recognition sequences. The decrease of the two electrochemical redox signals indicated the presence of dual-target miRNAs. With the robustness and high specificity of PER amplification and CRISPR/Cas9 cleavage system, simultaneous detection of two targets was achieved and the detection limits for miRNA-21 and miRNA-155 were 0.43 fM and 0.12 fM, respectively. The developed biosensor has the advantages of low cost, easy operation, and in-situ detection, providing a promising platform for point-of-care detection of multiple miRNAs.

Identification of Differentially Expressed Genes and microRNAs in the Gray and White Feather Follicles of Shitou Geese

The Shitou goose, a highly recognized indigenous breed with gray plumage originating from Chaozhou Raoping in Guangdong Province, China, is renowned for being the largest goose species in the country. Notably, during the pure breeding process of Shitou geese, approximately 2% of the offspring in each generation unexpectedly exhibited white plumage. To better understand the mechanisms underlying white plumage color formation in Shitou geese, we conducted a comparative transcriptome analysis between white and gray feather follicles, aiming to identify key genes and microRNAs that potentially regulate white plumage coloration in this unique goose breed. Our results revealed a number of pigmentation genes, encompassing TYR, TYRP1, EDNRB2, MLANA, SOX10, SLC45A2, GPR143, TRPM1, OCA2, ASIP, KIT, and SLC24A5, which were significantly down-regulated in the white feather follicles of Shitou geese. Among these genes, EDNRB2 and KIT emerged as the most promising candidate genes for white plumage coloration in Shitou geese. Additionally, our analysis also uncovered 46 differentially expressed miRNAs. Of these, miR-144-y may play crucial roles in the regulation of feather pigmentation. Furthermore, the expression of novel-m0086-5p, miR-489-y, miR-223-x, miR-7565-z, and miR-3535-z exhibits a significant negative correlation with the expression of pigmentation genes including TYRP1, EDNRB2, MLANA, SOX10, TRPM1, and KIT, suggesting these miRNAs may indirectly regulate the expression of these genes, thereby influencing feather color. Our findings provide valuable insights into the genetic mechanisms underlying white plumage coloration in Shitou geese and contribute to the broader understanding of avian genetics and coloration research.

A novel microRNA miR-4433a-3p as a potential diagnostic biomarker for lung adenocarcinoma

Background

Lung adenocarcinoma is one of the leading causes of cancer-related deaths because of the lack of early specific clinical indicators. MicroRNAs (miRNAs) have become the focus in lung cancer diagnosis. Further studies are required to explore miRNA expression in the serum of lung adenocarcinoma patients and their correlation with therapy and analyse specific messenger RNA targets to improve the specificity and sensitivity of early diagnosis.

Methods

The Toray 3D-Gene miRNA array was used to compare the expression levels of various miRNAs in the sera of patients with lung adenocarcinoma and healthy volunteers. Highly expressed miRNAs were selected for further analysis. To verify the screening results, serum and pleural fluid samples were analysed using qRT-PCR. Serum levels of the miRNAs and their correlation with the clinical information of patients with lung adenocarcinoma were analysed. The functions of miRNAs were further analysed using the Kyoto Encyclopedia of Gene and Genomes and Gene Ontology databases.

Results

Microarray analysis identified 60 and 50 miRNAs with upregulated and downregulated expressions, respectively, in the serum of patients with lung adenocarcinoma compared to those in healthy individuals. Using qRT-qPCR to detection of miRNAs expression in the serum or pleural effusion of patients with early and advanced lung adenocarcinoma, we found that miR-4433a-3p could be used as a diagnostic marker and therapeutic evaluation indicator for lung adenocarcinoma. Serum of miR-4433a-3p levels significantly correlated with the clinical stage. miR-4433a-3p may be more suitable than other tumour markers for the early diagnosis and evaluation of therapeutic effects in lung adenocarcinoma. miR-4433a-3p may affect tumour growth and metastasis by acting on target genes (PIK3CD, UBE2J2, ICMT, PRDM16 and others) and regulating tumour-related signalling pathways (MAPK signal pathway, Ras signalling pathway and others).

Conclusion

miR-4433a-3p may serve as a biomarker for the early diagnosis of lung adenocarcinoma and monitoring of therapeutic effects.

A time course analysis through diapause reveals dynamic temporal patterns of microRNAs associated with endocrine regulation in the butterfly Pieris napi

Organisms inhabiting highly seasonal environments must cope with a wide range of environmentally induced challenges. Many seasonal challenges require extensive physiological modification to survive. In winter, to survive extreme cold and limited resources, insects commonly enter diapause, which is an endogenously derived dormant state associated with minimized cellular processes and low energetic expenditure. Due to the high degree of complexity involved in diapause, substantial cellular regulation is required, of which our understanding primarily derives from the transcriptome via messenger RNA expression dynamics. Here we aim to advance our understanding of diapause by investigating microRNA (miRNA) expression in diapausing and direct developing pupae of the butterfly Pieris napi. We identified coordinated patterns of miRNA expression throughout diapause in both head and abdomen tissues of pupae, and via miRNA target identification, found several expression patterns to be enriched for relevant diapause-related physiological processes. We also identified two candidate miRNAs, miR-14-5p and miR-2a-3p, that are likely involved in diapause progression through their activity in the ecdysone pathway, a critical regulator of diapause termination. miR-14-5p targets phantom, a gene in the ecdysone synthesis pathway, and is upregulated early in diapause. miR-2a-3p has been found to be expressed in response to ecdysone, and is upregulated during diapause termination. Together, the expression patterns of these two miRNAs match our current understanding of the timing of hormonal regulation of diapause in P. napi and provide interesting candidates to further explore the mechanistic role of microRNAs in diapause regulation.

Mechanisms underlying TRPV4-mediated regulation of miR-146a expression

Persistent inflammation is a major contributor in the development of various inflammatory diseases like atherosclerosis. Our study investigates how transient receptor potential vanilloid 4 (TRPV4), a mechanosensitive ion channel, interacts with microRNA-146a (miR-146a), within the context of inflammation and atherosclerosis. Micro-RNAs play a critical role in controlling gene expression, and miR-146a is notable for its anti-inflammatory actions. TRPV4 is activated by diverse soluble and mechanical stimuli, and often associated with inflammatory responses in various diseases. Here, we find that TRPV4 negatively regulates miR-146a expression in macrophages, especially following stimulation by lipopolysaccharides or alterations in matrix stiffness. We show that in atherosclerosis, a condition characterized by matrix stiffening, TRPV4 decreases miR-146a expression in aortic tissue macrophages. We find that TRPV4's impact on miR-146a is independent of activation of NFκB, Stat1, P38, and AKT, but is rather mediated through a mechanism involving histone deacetylation instead of DNA methylation at the miR-146a promoter site. Furthermore, we show that N-terminal residues 1 to 130 in TRPV4 is essential in suppression of miR-146a expression in LPS-stimulated macrophages. Altogether, this study identifies a regulatory mechanism of miR-146a expression by TRPV4 which may open new potential therapeutic strategies for managing inflammatory diseases.

Human proximal tubular epithelial cell-derived small extracellular vesicles mediate synchronized tubular ferroptosis in hypoxic kidney injury

Hypoxia is the key pathobiological trigger of tubular oxidative stress and cell death that drives the transition of acute kidney injury (AKI) to chronic kidney disease (CKD). The mitochondrial-rich proximal tubular epithelial cells (PTEC) are uniquely sensitive to hypoxia and thus, are pivotal in propagating the sustained tubular loss of AKI-to-CKD transition. Here, we examined the role of PTEC-derived small extracellular vesicles (sEV) in propagating the 'wave of tubular death'. Ex vivo patient-derived PTEC were cultured under normoxia (21 % O2) and hypoxia (1 % O2) on Transwell inserts for isolation and analysis of sEV secreted from apical versus basolateral PTEC surfaces. Increased numbers of sEV were secreted from the apical surface of hypoxic PTEC compared with normoxic PTEC. No differences in basolateral sEV numbers were observed between culture conditions. Biological pathway analysis of hypoxic-apical sEV cargo identified distinct miRNAs linked with cellular injury pathways. In functional assays, hypoxic-apical sEV selectively induced ferroptotic cell death (↓glutathione peroxidase-4, ↑lipid peroxidation) in autologous PTEC compared with normoxic-apical sEV. The addition of ferroptosis inhibitors, ferrostatin-1 and baicalein, attenuated PTEC ferroptosis. RNAse A pretreatment of hypoxic-apical sEV also abrogated PTEC ferroptosis, demonstrating a role for sEV RNA in ferroptotic 'wave of death' signalling. In line with these in vitro findings, in situ immunolabelling of diagnostic kidney biopsies from AKI patients with clinical progression to CKD (AKI-to-CKD transition) showed evidence of ferroptosis propagation (increased numbers of ACSL4+ PTEC), while urine-derived sEV (usEV) from these 'AKI-to-CKD transition' patients triggered PTEC ferroptosis (↑lipid peroxidation) in functional studies. Our data establish PTEC-derived apical sEV and their intravesicular RNA as mediators of tubular lipid peroxidation and ferroptosis in hypoxic kidney injury. This concept of how tubular pathology is propagated from the initiating insult into a 'wave of death' provides novel therapeutic check-points for targeting AKI-to-CKD transition.

Distinctive physical properties of DNA shared by RNA polymerase II gene promoters and 5'-flanking regions of tRNA genes

Numerous noncoding (nc)RNAs have been identified. Similar to the transcription of protein-coding (mRNA) genes, long noncoding (lnc)RNA genes and most of micro (mi)RNA genes are transcribed by RNA polymerase II (Pol II). In the transcription of mRNA genes, core promoters play an indispensable role; they support the assembly of the preinitiation complex (PIC). However, the structural and/or physical properties of the core promoters of lncRNA and miRNA genes remain largely unexplored, in contrast with those of mRNA genes. Using the core promoters of human genes, we analyzed the repertoire and population ratios of residing core promoter elements (CPEs) and calculated the following five DNA physical properties (DPPs): duplex DNA free energy, base stacking energy, protein-induced deformability, rigidity and stabilizing energy of Z-DNA. Here, we show that their CPE and DPP profiles are similar to those of mRNA gene promoters. Importantly, the core promoters of these three classes of genes have two highly distinctive sites in their DPP profiles around the TSS and position -27. Similar characteristics in DPPs are also found in the 5'-flanking regions of tRNA genes, indicating their common essential roles in transcription initiation over the kingdom of RNA polymerases.

Identification of novel microRNAs in the embryonic mouse brain using deep sequencing

Many advances in small RNA-seq technology and bioinformatics pipelines have been made recently, permitting the discovery of novel miRNAs in the embryonic day 15.5 (E15.5) mouse brain. We aimed to improve miRNA discovery in this tissue to expand our knowledge of the regulatory networks that underpin normal neurodevelopment, find new candidates for neurodevelopmental disorder aetiology, and deepen our understanding of non-coding RNA evolution. A high-quality small RNA-seq dataset of 458 M reads was generated. An unbiased miRNA discovery pipeline identified fifty putative novel miRNAs, six of which were selected for further validation. A combination of conservation analysis and target functional prediction was used to determine the authenticity of novel miRNA candidates. These findings demonstrate that miRNAs remain to be discovered, particularly if they have the features of other small RNA species.

Tissue-specific profiling of age-dependent miRNAomic changes in Caenorhabditis elegans

Ageing exhibits common and distinct features in various tissues, making it critical to decipher the tissue-specific ageing mechanisms. MiRNAs are essential regulators in ageing and are recently highlighted as a class of intercellular messengers. However, little is known about the tissue-specific transcriptomic changes of miRNAs during ageing. C. elegans is a well-established model organism in ageing research. Here, we profile the age-dependent miRNAomic changes in five isolated worm tissues. Besides the diverse ageing-regulated miRNA expression across tissues, we discover numerous miRNAs in the tissues without their transcription. We further profile miRNAs in the extracellular vesicles and find that worm miRNAs undergo inter-tissue trafficking via these vesicles in an age-dependent manner. Using these datasets, we uncover the interaction between body wall muscle-derived mir-1 and DAF-16/FOXO in the intestine, suggesting mir-1 as a messenger in inter-tissue signalling. Taken together, we systematically investigate worm miRNAs in the somatic tissues and extracellular vesicles during ageing, providing a valuable resource to study tissue-autonomous and nonautonomous functions of miRNAs in ageing.

Discovery of circulating miRNAs as biomarkers of chronic Chagas heart disease via a small RNA-Seq approach

Chagas disease affects approximately 7 million people worldwide in Latin America and is a neglected tropical disease. Twenty to thirty percent of chronically infected patients develop chronic Chagas cardiomyopathy decades after acute infection. Identifying biomarkers of Chagas disease progression is necessary to develop better therapeutic and preventive strategies. Circulating microRNAs are increasingly reliable biomarkers of disease and therapeutic targets. To identify new circulating microRNAs for Chagas disease, we performed exploratory small RNA sequencing from the plasma of patients and performed de novo miRNA prediction, identifying potential new microRNAs. The levels of the new microRNAs temporarily named miR-Contig-1519 and miR-Contig-3244 and microRNAs that are biomarkers for nonchagasic cardiomyopathies, such as miR-148a-3p and miR-224-5p, were validated by quantitative reverse transcription. We found a specific circulating microRNA signature defined by low miR-Contig-3244, miR-Contig-1519, and miR-148a-3 levels but high miR-224-5p levels for patients with chronic Chagas disease. Finally, we predicted in silico that these altered circulating microRNAs could affect the expression of target genes involved in different cellular pathways and biological processes, which we will explore in the future.

Analysis of extracellular vesicle microRNA profiles reveals distinct blood and lymphatic endothelial cell origins

Extracellular vesicles (EVs) are crucial mediators of cell-to-cell communication in physiological and pathological conditions. Specifically, EVs released from the vasculature into blood were found to be quantitatively and qualitatively different in diseases compared to healthy states. However, our understanding of EVs derived from the lymphatic system is still scarce. In this study, we compared the mRNA and microRNA (miRNA) expression in blood vascular (BEC) and lymphatic (LEC) endothelial cells. After characterization of the EVs by fluorescence-triggered flow cytometry, nanoparticle tracking analysis and cryo-transmission electron microscopy (cryo-TEM) we utilized small RNA-sequencing to characterize miRNA signatures in the EVs and identify cell-type specific miRNAs in BEC and LEC. We found miRNAs specifically enriched in BEC and LEC on the cellular as well as the extracellular vesicle level. Our data provide a solid basis for further functional in vitro and in vivo studies addressing the role of EVs in the blood and lymphatic vasculature.

Resequencing of the TMF-1 (TATA Element Modulatory Factor) regulated protein (TRNP1) gene in domestic and wild canids

BACKGROUND

Cortical folding is related to the functional organization of the brain. The TMF-1 regulated protein (TRNP1) regulates the expansion and folding of the mammalian cerebral cortex, a process that may have been accelerated by the domestication of dogs. The objectives of this study were to sequence the TRNP1 gene in dogs and related canid species, provide evidence of its expression in dog brain and compare the genetic variation within dogs and across the Canidae. The gene was located in silico to dog chromosome 2. The sequence was experimentally confirmed by amplifying and sequencing the TRNP1 exonic and promoter regions in 72 canids (36 purebred dogs, 20 Gy wolves and wolf-dog hybrids, 10 coyotes, 5 red foxes and 1 Gy fox).

RESULTS

A partial TRNP1 transcript was isolated from several regions in the dog brain. Thirty genetic polymorphisms were found in the Canis sp. with 17 common to both dogs and wolves, and only one unique to dogs. Seven polymorphisms were observed only in coyotes. An additional 9 variants were seen in red foxes. Dogs were the least genetically diverse. Several polymorphisms in the promoter and 3'untranslated region were predicted to alter TRNP1 function by interfering with the binding of transcriptional repressors and miRNAs expressed in neural precursors. A c.259_264 deletion variant that encodes a polyalanine expansion was polymorphic in all species studied except for dogs. A stretch of 15 nucleotides that is found in other mammalian sequences (corresponding to 5 amino acids located between Pro58 and Ala59 in the putative dog protein) was absent from the TRNP1 sequences of all 5 canid species sequenced. Both of these aforementioned coding sequence variations were predicted to affect the formation of alpha helices in the disordered region of the TRNP1 protein.

CONCLUSIONS

Potentially functionally important polymorphisms in the TRNP1 gene are found within and across various Canis species as well as the red fox, and unique differences in protein structure have evolved and been conserved in the Canidae compared to all other mammalian species.

The microRNA-mediated gene regulatory network in the hippocampus and hypothalamus of the aging mouse

Aging leads to time-dependent functional decline of all major organs. In particular, the aging brain is prone to cognitive decline and several neurodegenerative diseases. Various studies have attempted to understand the aging process and underlying molecular mechanisms by monitoring changes in gene expression in the aging mouse brain using high-throughput sequencing techniques. However, the effect of microRNA (miRNA) on the post-transcriptional regulation of gene expression has not yet been comprehensively investigated. In this study, we performed global analysis of mRNA and miRNA expression simultaneously in the hypothalamus and hippocampus of young and aged mice. We identified aging-dependent differentially expressed genes, most of which were specific either to the hypothalamus or hippocampus. However, genes related to immune response-related pathways were enriched in upregulated differentially expressed genes, whereas genes related to metabolism-related pathways were enriched in downregulated differentially expressed genes in both regions of the aging brain. Furthermore, we identified many differentially expressed miRNAs, including three that were upregulated and three that were downregulated in both the hypothalamus and hippocampus. The two downregulated miRNAs, miR-322-3p, miR-542-3p, and the upregulated protein-encoding coding gene C4b form a regulatory network involved in complement and coagulation cascade pathways in the hypothalamus and hippocampus of the aging brain. These results advance our understanding of the miRNA-mediated gene regulatory network and its influence on signaling pathways in the hypothalamus and hippocampus of the aging mouse brain.

Glass Nanopipette-Based Plasmonic SERS Platform for Single-Cell MicroRNA-21 Sensing during Apoptosis

As one of the most widely distributed microRNAs, microRNA-21 (miRNA-21) significantly regulates target genes' expression levels and participates in many cellular and intercellular activities, and its abnormal expression is always related to some diseases, especially cancer. Hence, detecting miRNA-21, as a biomarker, at the single-cell level helps us to reveal cell heterogeneity and expression level variation during the state change of cells. In this study, we constructed a gold nanoparticles nanomembrane (AuNPs-NM)-modified plasmonic glass nanopipette (P-nanopipette) surface-enhanced Raman scattering (SERS) sensing platform to sensitively detect content variation of the intracellular miRNA-21 during the electrostimulus (ES)-induced apoptosis process. The cytoplasm-located miRNA-21 was first extracted by using the extraction DNA (HP1)-modified P-nanopipette through a hybridization chain reaction (HCR). The nanopipette was then incubated with a labeling DNA (HP2) and reporter 4-MBA-modified Raman tag. The Raman signal (collected from the tip area near the orifice within 1 μm) showed a good response to the content variation of intracellular miRNA-21 under ES, and the proposed single-cell SERS detection platform provides a simple way to study intracellular substance change and evaluate cancer treatment outcomes.

Substances in the mandibular glands mediate queen effects on larval development and colony organization in an annual bumble bee

Social organization is commonly dynamic, with extreme examples in annual social insects, but little is known about the underlying signals and mechanisms. Bumble bee larvae with close contact to a queen do not differentiate into gynes, pupate at an earlier age, and are commonly smaller than siblings that do not contact a queen. We combined detailed observations, proteomics, microRNA transcriptomics, and gland removal surgery to study the regulation of brood development and division of labor in the annual social bumble bee Bombus terrestris. We found that regurgitates fed to larvae by queens and workers differ in their protein and microRNA composition. The proteome of the regurgitate overlaps significantly with that of the mandibular (MG) and hypopharyngeal glands (HPG), suggesting that these exocrine glands are sources of regurgitate proteins. The proteome of the MG and HPG, but not the salivary glands, differs between queens and workers, with caste-specificity preserved for the MG and regurgitate proteomes. Queens subjected to surgical removal of the MG showed normal behavior, brood care, and weight gain, but failed to shorten larval development. These findings suggest that substances in the queen MG are fed to larvae and influence their developmental program. We suggest that when workers emerge and contribute to larval feeding, they dilute the effects of the queen substances, until she can no longer manipulate the development of all larvae. Longer developmental duration may allow female larvae to differentiate into gynes rather than to workers, mediating the colony transition from the ergonomic to the reproductive phase.

Integrated Transcriptome Analysis of miRNAs and mRNAs in the Skeletal Muscle of Wuranke Sheep

MicroRNAs (miRNAs) are regarded as important regulators in skeletal muscle development. To reveal the regulatory roles of miRNAs and their target mRNAs underlying the skeletal muscle development of Wuranke sheep, we investigated the miRNA and mRNA expression profiles in the biceps femoris of these sheep at the fetal (3 months of gestation) and 3- and 15-month-old postnatal stages. Consequently, a total of 1195 miRNAs and 24,959 genes were identified. Furthermore, 474, 461, and 54 differentially expressed miRNAs (DEMs) and 6783, 7407, and 78 differentially expressed genes (DEGs) were detected among three comparative groups. Functional analysis demonstrated that the target mRNAs of the DEMs were enriched in multiple pathways related to muscle development. Moreover, the interactions among several predicted miRNA-mRNA pairs (oar-miR-133-HDAC1, oar-miR-1185-5p-MYH1/HADHA/OXCT1, and PC-5p-3703_578-INSR/ACTG1) that potentially affect skeletal muscle development were verified using dual-luciferase reporter assays. In this study, we identified the miRNA and mRNA differences in the skeletal muscle of Wuranke sheep at different developmental stages and revealed that a series of candidate miRNA-mRNA pairs may act as modulators of muscle development. These results will contribute to future studies on the function of miRNAs and their target mRNAs during skeletal muscle development in Wuranke sheep.

The Identification and Comparative Analysis of Non-Coding RNAs in Spores and Mycelia of Penicillium expansum

Penicillium expansum is the most popular post-harvest pathogen and causes blue mold disease in pome fruit and leads to significant economic losses worldwide every year. However, the fundamental regulation mechanisms of growth in P. expansum are unclear. Recently, non-coding RNAs (ncRNAs) have attracted more attention due to critical roles in normalizing gene expression and maintaining cellular genotypes in organisms. However, the research related to ncRNAs in P. expansum have not been reported. Therefore, to provide an overview of ncRNAs on composition, distribution, expression changes, and potential targets in the growth process, a comparative transcriptomic analysis was performed on spores and mycelia of P. expansum in the present study. A total of 2595 novel mRNAs, 3362 long non-coding RNAs (lncRNAs), 10 novel microRNAs (miRNAs), 86 novel small interfering RNAs (siRNAs), and 11,238 circular RNAs (circRNAs) were predicted and quantified. Of these, 1482 novel mRNAs, 5987 known mRNAs, 2047 lncRNAs, 40 miRNAs, 38 novel siRNAs, and 9235 circRNAs were differentially expressed (DE) in response to the different development stages. Afterward, the involved functions and pathways of DE RNAs were revealed via Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database enrichment analysis. The interaction networks between mRNAs, lncRNAs, and miRNAs were also predicted based on their correlation coefficient of expression profiles. Among them, it was found that miR168 family members may play important roles in fungal growth due to their central location in the network. These findings will contribute to a better understanding on regulation machinery at the RNA level on fungal growth and provide a theoretical basis to develop novel control strategies against P. expansum.

Genome-wide transcriptomics and microRNAomics analyses uncover multi-faceted mechanisms to cope with copper stress in ancient macrobenthos amphioxus

The mechanisms underlying the toxicity of environmental stress are unclear for marine macrobenthos. Copper/Cu has posed the most serious threats to amphioxus, an ancient and model benthic cephalochordate. Herein, a dynamic change in the physiological parameters (GR, SOD, ATP, and MDA) was detected with ROS accumulation in Branchiostoma belcheri exposed to 0.3 mg·L-1 Cu. Transcriptomes and microRNAomes of B. belcheri were generated to investigate the molecular mechanisms by which this amphioxus copes with Cu exposure. Time-specific genes identified at different time points after exposure were involved in the stimulus and immune response, detoxification and ionic homeostasis, aging and the nervous system, sequentially, with prolongation of exposure time, forming a dynamic process of molecular response to Cu stress. In total, 57 differentially expressed miRNAs were identified under Cu stress. Transcriptomics-miRNAomics analyses indicate that these miRNAs targeted genes associated with many key biological processes such as xenobiotics degradation, oxidative stress, and energy metabolism. The constructed miRNA-mRNA-pathway network uncovered a broad post-transcriptional regulatory mechanism in B. belcheri to cope with Cu stress. Overall, this integrated analyses show that enhanced defense response, accelerated ROS elimination, and repressed ATP production constitute a comprehensive strategy to cope with Cu toxicity in the ancient macrobenthos.

Catalytic Hairpin Assembly-Enhanced Graphene Transistor for Ultrasensitive miRNA Detection

MicroRNAs (miRNAs) have emerged as powerful biomarkers for disease diagnosis and screening. Traditional miRNA analytical techniques are inadequate for point-of-care testing due to their reliance on specialized expertise and instruments. Graphene field-effect transistors (GFETs) offer the prospect of simple and label-free diagnostics. Herein, a GFET biosensor based on tetrahedral DNA nanostructure (TDN)-assisted catalytic hairpin assembly (CHA) reaction (TCHA) has been fabricated and applied to the sensitive and specific detection of miRNA-21. TDN structures are assembled to construct the biosensing interface, facilitating CHA reaction by providing free space and preventing unwanted entanglements, aggregation, and adsorption of probes on the graphene channel. Owing to synergistic effects of TDN-assisted in situ nucleic acid amplification on the sensing surface, as well as inherent signal sensitization of GFETs, the biosensor exhibits ultrasensitive detection of miRNA-21 down to 5.67 × 10-19 M, approximately three orders of magnitude lower than that normally achieved by graphene transistors with channel functionalization of single-stranded DNA probes. In addition, the biosensor demonstrates excellent analytical performance regarding selectivity, stability, and reproducibility. Furthermore, the practicability of the biosensor is verified by analyzing targets in a complex serum environment and cell lysates, showing tremendous potential in bioanalysis and clinical diagnosis.

Accurate microRNA annotation of animal genomes using trained covariance models of curated microRNA complements in MirMachine

The annotation of microRNAs depends on the availability of transcriptomics data and expert knowledge. This has led to a gap between the availability of novel genomes and high-quality microRNA complements. Using >16,000 microRNAs from the manually curated microRNA gene database MirGeneDB, we generated trained covariance models for all conserved microRNA families. These models are available in our tool MirMachine, which annotates conserved microRNAs within genomes. We successfully applied MirMachine to a range of animal species, including those with large genomes and genome duplications and extinct species, where small RNA sequencing is hard to achieve. We further describe a microRNA score of expected microRNAs that can be used to assess the completeness of genome assemblies. MirMachine closes a long-persisting gap in the microRNA field by facilitating automated genome annotation pipelines and deeper studies into the evolution of genome regulation, even in extinct organisms.

Transcriptomic analyses of ovarian clear-cell carcinoma with concurrent endometriosis

Introduction

Endometriosis, a benign inflammatory disease whereby endometrial-like tissue grows outside the uterus, is a risk factor for endometriosis-associated ovarian cancers. In particular, ovarian endometriomas, cystic lesions of deeply invasive endometriosis, are considered the precursor lesion for ovarian clear-cell carcinoma (OCCC).

Methods

To explore this transcriptomic landscape, OCCC from women with pathology-proven concurrent endometriosis (n = 4) were compared to benign endometriomas (n = 4) by bulk RNA and small-RNA sequencing.

Results

Analysis of protein-coding genes identified 2449 upregulated and 3131 downregulated protein-coding genes (DESeq2, P< 0.05, log2 fold-change > |1|) in OCCC with concurrent endometriosis compared to endometriomas. Gene set enrichment analysis showed upregulation of pathways involved in cell cycle regulation and DNA replication and downregulation of pathways involved in cytokine receptor signaling and matrisome. Comparison of pathway activation scores between the clinical samples and publicly-available datasets for OCCC cell lines revealed significant molecular similarities between OCCC with concurrent endometriosis and OVTOKO, OVISE, RMG1, OVMANA, TOV21G, IGROV1, and JHOC5 cell lines. Analysis of miRNAs revealed 64 upregulated and 61 downregulated mature miRNA molecules (DESeq2, P< 0.05, log2 fold-change > |1|). MiR-10a-5p represented over 21% of the miRNA molecules in OCCC with endometriosis and was significantly upregulated (NGS: log2fold change = 4.37, P = 2.43e-18; QPCR: 8.1-fold change, P< 0.05). Correlation between miR-10a expression level in OCCC cell lines and IC50 (50% inhibitory concentration) of carboplatin in vitro revealed a positive correlation (R2 = 0.93). MiR-10a overexpression in vitro resulted in a significant decrease in proliferation (n = 6; P< 0.05) compared to transfection with a non-targeting control miRNA. Similarly, the cell-cycle analysis revealed a significant shift in cells from S and G2 to G1 (n = 6; P< 0.0001). Bioinformatic analysis predicted that miR-10a-5p target genes that were downregulated in OCCC with endometriosis were involved in receptor signaling pathways, proliferation, and cell cycle progression. MiR-10a overexpression in vitro was correlated with decreased expression of predicted miR-10a target genes critical for proliferation, cell-cycle regulation, and cell survival including [SERPINE1 (3-fold downregulated; P< 0.05), CDK6 (2.4-fold downregulated; P< 0.05), and RAP2A (2-3-fold downregulated; P< 0.05)].

Discussion

These studies in OCCC suggest that miR-10a-5p is an impactful, potentially oncogenic molecule, which warrants further studies.

MiR-196b-3p and miR-450b-3p are key regulators of adipogenesis in porcine intramuscular and subcutaneous adipocytes

BACKGROUND

As components of white adipose tissue, porcine intramuscular (IM) and subcutaneous (SC) adipocytes undergo similar differentiation and adipogenesis processes. However, the adipogenic capacity of IM adipocytes is weaker than that of SC adipocytes. Identifying key regulators underlying this difference between IM and SC adipocytes will benefit pig breeding.

RESULTS

In this study, we used BGISEQ-500 sequencing technology to analyze the expression of small RNAs in primary cultured IM and SC adipocytes on day 8 after adipogenic induction, and found 32-fold higher miR-196b-3p expression, as well as 8-fold lower miR-450b-3p expression in IM adipocytes than in SC adipocytes. Functional studies revealed that miR-196b-3p inhibits adipogenesis by targeting CD47 via the AMPK signaling pathway, and its effect was attenuated by the specific p-AMPKα activator AICAR. We also found that miR-450b-3p promotes adipogenesis by targeting SIRT1 via the Wnt/β-catenin signaling pathway, and its effect was weakened by the Wnt/β-catenin signaling activator LiCl.

CONCLUSIONS

Our findings suggest that miR-196b-3p and miR-450b-3p are novel key regulatory factors that play opposite roles in porcine adipogenesis, helping us decipher the regulatory differences between porcine IM and SC fat deposition.

Characterizations of microRNAs involved in the molecular mechanisms underlying the therapeutic effects of noni (Morinda citrifolia L.) fruit juice on hyperuricemia in mice

Background

Hyperuricemia is generally defined as the high level of serum uric acid and is well known as an important risk factor for the development of various medical disorders. However, the medicinal treatment of hyperuricemia is frequently associated with multiple side-effects.

Methods

The therapeutic effect of noni (Morinda citrifolia L.) fruit juice on hyperuricemia and the underlying molecular mechanisms were investigated in mouse model of hyperuricemia induced by potassium oxonate using biochemical and high-throughput RNA sequencing analyses.

Results

The levels of serum uric acid (UA) and xanthine oxidase (XOD) in mice treated with noni fruit juice were significantly decreased, suggesting that the noni fruit juice could alleviate hyperuricemia by inhibiting the XOD activity and reducing the level of serum UA. The contents of both serum creatinine and blood urine nitrogen of the noni fruit juice group were significantly lower than those of the model group, suggesting that noni fruit juice promoted the excretion of UA without causing deleterious effect on the renal functions in mice. The differentially expressed microRNAs involved in the pathogenesis of hyperuricemia in mice were identified by RNA sequencing with their target genes further annotated based on both Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases to explore the metabolic pathways and molecular mechanisms underlying the therapeutic effect on hyperuricemia by noni fruit juice.

Conclusion

Our study provided strong experimental evidence to support the further investigations of the potential application of noni fruit juice in the treatment of hyperuricemia.

Intracellular metabolomics and microRNAomics unveil new insight into the regulatory network for potential biocontrol mechanism of stress-tolerant Tricho-fusants interacting with phytopathogen Sclerotium rolfsii Sacc

The present study employed microRNA (miRNA) sequencing and metabolome profiling of Trichoderma parental strains and fusants during normal growth and interaction with the phytopathogen Sclerotium rolfsii Sacc. In-vitro antagonism indicated that abiotic stress-tolerant Tricho-fusant FU21 was examined as a potent biocontroller with mycoparasitic action after 10 days. During interaction with the test pathogen, the most abundant uprising intracellular metabolite was recognized as l-proline, which corresponds to held-down l-alanine, associated with arginine and proline metabolism, biosynthesis of secondary metabolites, and nitrogen metabolism linked to predicted genes controlled by miRNAs viz., cel-miR-8210-3p, hsa-miR-3613-5p, and mml-miR-7174-3p. The miRNAs- mml-miR-320c and mmu-miR-6980-5p were found to be associated with phenylpropanoid biosynthesis, transcription factors, and signal transduction pathways, respectively, and were ascertained downregulated in potent FU21_IB compared with FU21_CB. The amino benzoate degradation and T cell receptor signaling pathways were regulated by miRNAs cel-miR-8210 and tca-miR-3824 as stress tolerance mechanisms of FU21. The intracellular metabolites l-proline, maleic acid, d-fructose, Myo-inositol, arabinitol, d-xylose, mannitol, and butane were significantly elevated as potential biocontrol and stress-tolerant constituents associated with miRNA regulatory pathways in potent FU21_IB. A network analysis between regulatory miRNA-predicted genes and intracellular metabolomics acknowledged possible biocontrol pathways/mechanisms in potent FU21_IB to restrain phytopathogen.

Construction of ceRNA network based on RNA-seq for identifying prognostic lncRNA biomarkers in Perthes disease

Introduction: Legg-Calvé-Perthes disease or Perthes disease is a condition that occurs in children aged 2 to 15 years, and is characterized by osteonecrosis of the femoral head, which results in physical limitations. Despite ongoing research, the pathogenesis and molecular mechanisms underlying the development of Perthes disease remain unclear. In order to obtain further insights, the expression patterns of long non-coding RNAs (lncRNAs), miRNAs, and mRNAs in a rabbit model of Perthes disease were analyzed in this study by transcriptome sequencing. Methods and results: The results of RNA-seq analyses revealed that 77 lncRNAs, 239 miRNAs, and 1027 mRNAs were differentially expressed in the rabbit model. This finding suggested that multiple genetic pathways are involved in the development of Perthes disease. A weighted gene co-expression network analysis (WGCNA) network was subsequently constructed using the differentially expressed mRNAs (DEmRNAs), and network analysis revealed that the genes associated with angiogenesis and platelet activation were downregulated, which was consistent with the findings of Perthes disease. A competing endogenous RNA (ceRNA) network was additionally constructed using 29 differentially expressed lncRNAs (including HIF3A and LOC103350994), 28 differentially expressed miRNAs (including ocu-miR-574-5p and ocu-miR-324-3p), and 76 DEmRNAs (including ALOX12 and PTGER2). Disscusion: The results obtained herein provide novel perspectives regarding the pathogenesis and molecular mechanisms underlying the development of Perthes disease. The findings of this study can pave the way for the development of effective therapeutic strategies for Perthes disease in future.

Application of annotation-agnostic RNA sequencing data analysis tools for biomarker discovery in liquid biopsy

RNA sequencing analysis is an important field in the study of extracellular vesicles (EVs), as these particles contain a variety of RNA species that may have diagnostic, prognostic and predictive value. Many of the bioinformatics tools currently used to analyze EV cargo rely on third-party annotations. Recently, analysis of unannotated expressed RNAs has become of interest, since these may provide complementary information to traditional annotated biomarkers or may help refine biological signatures used in machine learning by including unknown regions. Here we perform a comparative analysis of annotation-free and classical read-summarization tools for the analysis of RNA sequencing data generated for EVs isolated from persons with amyotrophic lateral sclerosis (ALS) and healthy donors. Differential expression analysis and digital-droplet PCR validation of unannotated RNAs also confirmed their existence and demonstrates the usefulness of including such potential biomarkers in transcriptome analysis. We show that find-then-annotate methods perform similarly to standard tools for the analysis of known features, and can also identify unannotated expressed RNAs, two of which were validated as overexpressed in ALS samples. We demonstrate that these tools can therefore be used for a stand-alone analysis or easily integrated into current workflows and may be useful for re-analysis as annotations can be integrated post hoc.

Identification and comparative analysis of miRNA transcriptomes after allograft and xenograft transplantation in Pinctada fucata martensii

Effective immune regulation after transplantation during pearl production is crucial for the cultivation of high-quality pearls. MicroRNAs (miRNAs) play an important role in a variety of physiological processes. To understand the regulatory rules of miRNAs after transplantation in Pinctada funcata martensii, we constructed 13 miRNA transcriptomes, including the control group (Con), allograft (Al), and xenograft (Xe) transplantation at six time points (6, 12, and 24 h and 3, 6, and 12 days), in which the xenografted mantle tissue was from Pinctada maxima. We identified 159 differentially expressed miRNAs (DEMs) and found that these DEMs showed high expression at 12 h, 24 h, and 3 days after transplantation. A total of 130 DEMs, such as Let-7, were present in the Al and Xe groups; miR-34 and 16 other DEMs were specifically present in the Al group; miR-216b and 13 other DEMs were specifically present in the Xe group. Compared with the Con group, the target genes of DEMs in the Al group were significantly enriched in protein complex, cytoskeleton, and macromolecular complex, and the Xe group was significantly enriched in ribonucleoside metabolic process, nucleoside binding, and cell division. Compared with the Al group, the target genes in the Xe group were significantly enriched in response to DNA damage stimulation. Overall, multiple pathways associated with cellular activity were enriched in higher numbers of genes in the Xe group than in the Al group. These findings enriched the information on immune regulatory mechanisms at the expression level of miRNAs in P. f. martensii after transplantation.

Symbiosis preservation: Putative regulation of fatty acyl-CoA reductase by miR-31a within the symbiont harboring bacteriome through tsetse evolution

Tsetse flies are the sole vectors of African trypanosomes. In addition to trypanosomes, tsetse harbor obligate Wigglesworthia glossinidia bacteria that are essential to tsetse biology. The absence of Wigglesworthia results in fly sterility, thus offering promise for population control strategies. Here, microRNA (miRNAs) and mRNA expression are characterized and compared between the exclusive Wigglesworthia-containing bacteriome and adjacent aposymbiotic tissue in females of two evolutionarily distant tsetse species (Glossina brevipalpis and G. morsitans). A total of 193 miRNAs were expressed in either species, with 188 of these expressed in both species, 166 of these were novel to Glossinidae, and 41 miRNAs exhibited comparable expression levels between species. Within bacteriomes, 83 homologous mRNAs demonstrated differential expression between G. morsitans aposymbiotic and bacteriome tissues, with 21 of these having conserved interspecific expression. A large proportion of these differentially expressed genes are involved in amino acid metabolism and transport, symbolizing the essential nutritional role of the symbiosis. Further bioinformatic analyses identified a sole conserved miRNA::mRNA interaction (miR-31a::fatty acyl-CoA reductase) within bacteriomes likely catalyzing the reduction of fatty acids to alcohols which comprise components of esters and lipids involved in structural maintenance. The Glossina fatty acyl-CoA reductase gene family is characterized here through phylogenetic analyses to further understand its evolutionary diversification and the functional roles of members. Further research to characterize the nature of the miR-31a::fatty acyl-CoA reductase interaction may find novel contributions to the symbiosis to be exploited for vector control.

Exploring conserved and novel MicroRNA-like small RNAs from stress tolerant Trichoderma fusants and parental strains during interaction with fungal phytopathogen Sclerotium rolfsii Sacc

The study investigated potential microRNA-like small RNAs (milRNAs) from multi-stress-tolerant Tricho-fusants and parental strains (P1- Trichoderma virens NBAIITvs12 and P2- Trichoderma koningii MTCC796) for antagonistic activity during interaction with phytopathogen Sclerotium rolfsii. The Trichoderma was cultured in-vitro, with and without antagonism, against the pathogen and total RNA was extracted followed by small RNA library construction and sequencing. The milRNAs were identified by mapping high-quality unique reads against a reference genome. The milRNAs were recognized higher in antagonist Trichoderma during interaction with test pathogen compared to normal growth. The novel milRNAs candidates were found to vary during interaction with the pathogen and normal growth. The gene ontology and functional analysis illustrated that a total of 5828 potential targeted genes were recognized for 93 milRNAs of potent Fu21_IB and 3053 genes for 62 milRNAs of least fusant Fu28_IL. Functional annotation of milRNA-predicted genes integrating KEGG pathways indicates new insights into regulatory mechanisms, by interfering with milRNAs, associated with signal transduction, amino sugar metabolism, benzoate degradation, amino acid metabolism, and steroid and alkaloid metabolism for potential biocontrol of stress-tolerant Tricho-fusant FU21 during interaction with S. rolfsii. The present investigation is the first report of conserved and novel milRNAs from Tricho-fusants and parental strains interacting with S. rolfsii.

Expression analysis and targets prediction of microRNAs in OGD/R treated astrocyte-derived exosomes by smallRNA sequencing

Astrocytes activate and crosstalk with neurons influencing inflammatory responses following ischemic stroke. The distribution, abundance, and activity of microRNAs in astrocytes-derived exosomes after ischemic stroke remains largely unknown. In this study, exosomes were extracted from primary cultured mouse astrocytes via ultracentrifugation, and exposed to oxygen glucose deprivation/re‑oxygenation injury to mimic experimental ischemic stroke. SmallRNAs from astrocyte-derived exosomes were sequenced, and differentially expressed microRNAs were randomly selected and verified by stem-loop real time quantitative polymerase chain reaction. We found that 176 microRNAs, including 148 known and 28 novel microRNAs, were differentially expressed in astrocyte-derived exosomes following oxygen glucose deprivation/re‑oxygenation injury. In gene ontology enrichment, Kyoto encyclopedia of genes and genomes pathway analyses, and microRNA target gene prediction analyses, these alteration in microRNAs were associated to a broad spectrum of physiological functions including signaling transduction, neuroprotection and stress responses. Our findings warrant further investigating of these differentially expressed microRNAs in human diseases particularly ischemic stroke.

Circulating miRNAs are associated with frailty and ST-elevation myocardial infarction pathways

BACKGROUND

Frailty and ST-Elevation Myocardial Infarction (STEMI) share similar molecular pathways. Specific biomarkers, such as microRNAs (miRNAs), may provide insights into the molecular mechanisms that cause the relationship between frailty and STEMI.

OBJECTIVE

Our aim was to identify and compare circulating miRNA levels between frail and non-frail older adults following STEMI and comprehend the regulatory miRNA-gene networks and pathways involved in this condition.

METHODS

This exploratory study is a subanalysis of a larger observational study. In this study, we selected patients ≥ 65 years old, following STEMI, with pre-frail/frail (n=5) and non-frail (n=4) phenotype evaluated using the Clinical Frailty Scale and serum circulating miRNA levels were analyzed.

RESULTS

Pre-frail/frail patients had greater serum levels of 53 miRNAs, compared with non-frail patients. Notably, miR-103a-3p, miR-598-3p, and miR-130a-3p were the top three significantly deregulated miRNAs predicted to modulate gene expression associated with aging. Additional computational analyses showed 7,420 predicted miRNA gene targets, which were regulated by at least two of the 53 identified miRNAs. Pathway enrichment analysis showed that axon guidance and MAPK signaling were among pathways regulated by miRNA target genes.

CONCLUSIONS

These novel findings suggest a correlation between the identified miRNAs, target genes, and pathways in pre-frail and frail patients with myocardial infarction.

OsNAR2.1 induced endogenous nitrogen concentration variation affects transcriptional expression of miRNAs in rice

The studies of rice nitrogen concentration on the expression of miRNA so far are mostly limited to the exogenous nitrogen, leaving the effect of endogenous nitrogen largely unexplored. OsNAR2.1 is a high-affinity nitrate transporter partner protein which plays a central role in nitrate absorption and translocation in rice. The expression of OsNAR2.1 could influence the concentration of the endogenous nitrogen in rice. We showed that the expression and production of miRNA in rice can be influenced by manipulating the endogenous nitrogen concentration via OsNAR2.1 transgenic lines. The small RNA content, particularly 24 nucleotides small RNA, expressed differently in two transgenic rice lines (nitrogen efficient line with overexpression of OsNAR2.1 (Ov199), nitrogen-inefficient line with knockdown OsNAR2.1 by RNAi (RNAi)) compared to the wild-type (NP). Comparative hierarchical clustering expression pattern analysis revealed that the expression profiles of mature miRNA in both transgenic lines were different from NP. Several previously unidentified miRNAs were identified to be differentially expressed under different nitrogen concentrations, namely miR1874, miR5150, chr3-36147, chr4-27017 and chr5-21745. In conclusion, our findings suggest that the level of endogenous nitrogen concentration variation by overexpression or knockdown OsNAR2.1 could mediate the expression pattern and intensity of miRNA in rice, which is of high potential to be used in molecular breeding to improve the rice responses towards nitrogen utilization.

Identification and analysis of miRNAs differentially expressed in male and female Trichosanthes kirilowii maxim

Trichosanthes kirilowii Maxim. (TK) is a dioecious plant in the Cucurbitaceae family of which different sexes have separate medicinal uses. We used Illumina high-throughput sequencing technology to sequence miRNAs from male and female flower buds of TK. We performed bioinformatics analysis, miRNA identification, and target gene prediction on the data obtained from sequencing, and association analysis was performed in combination with the results of a previous transcriptome sequencing study. As a result, there were 80 differentially expressed miRNAs (DESs) between the female and male plants (48 upregulated and 32 downregulated in female plants). Moreover, 27 novel miRNAs in DESs were predicted to have 282 target genes, and 51 known miRNAs were predicted to have 3418 target genes. By establishing a regulatory network between miRNAs and target genes, 12 core genes were screened, including 7 miRNAs and 5 target genes. Among them, tkmiR157a-5p, tkmiR156c, tkmiR156_2, and tkmiR156k_2 jointly target the regulation of tkSPL18 and tkSPL13B. These two target genes are specifically expressed in male and female plants, respectively, and are involved in the biosynthesis process of BR, which is closely related to the sex differentiation process of TK. The identification of these miRNAs will provide a reference for the analysis of the sex differentiation mechanism of TK.

A comprehensive survey of C. elegans argonaute proteins reveals organism-wide gene regulatory networks and functions

Argonaute (AGO) proteins associate with small RNAs to direct their effector function on complementary transcripts. The nematode Caenorhabditis elegans contains an expanded family of 19 functional AGO proteins, many of which have not been fully characterized. In this work, we systematically analyzed every C. elegans AGO using CRISPR-Cas9 genome editing to introduce GFP::3xFLAG tags. We have characterized the expression patterns of each AGO throughout development, identified small RNA binding complements, and determined the effects of ago loss on small RNA populations and developmental phenotypes. Our analysis indicates stratification of subsets of AGOs into distinct regulatory modules, and integration of our data led us to uncover novel stress-induced fertility and pathogen response phenotypes due to ago loss.

Acute green tea intake attenuates circulating microRNA expression induced by a high-fat, high-saturated meal in obese women: A randomized crossover study

The objective of this study was to assess whether acute green tea (GT) supplementation attenuates inflammatory and oxidative stress biomarkers induced by high-fat, high-saturated (HFHS) meals in obese women, and to assess its ability to modulate circulating microRNA (miRNA) expression. This was a randomized, double-blind, crossover study. The study included obese women over 18 years old who had no comorbidities. In the first moment, patients were instructed to take 2 capsules of placebo or GT (738 mg) at 10:00 p.m. and to fast overnight. The next morning, a blood sample was collected, and an HFHS meal was offered to the patients. Another blood sample was collected 5 hours after the meal. In the second moment, patients who received placebo in the first moment now received the GT and vice-versa. Serum inflammatory and oxidative stress biomarkers were measured, and circulating levels of miRNA were evaluated. Fifteen women with mean age of 35.5±9.9 years were included in the final analysis. There was no difference regarding inflammatory and oxidative stress biomarkers. However, patients who consumed GT had lower circulating expression of 62 miRNAs compared with patients who did not consume GT. Predictive analysis of target genes showed 1,757 targets regulated by the 62 miRNAs. Notably, 5 miRNAs (miR-1297, miR-192-5p, miR-373-3p, miR-595 and miR-1266-5p) regulate genes associated with TGF-beta, CARM1, RSK, and BMP pathways. Our study showed that GT inhibited the expression of miRNAs induced by HFHS meal intake. These results shed light on the mechanisms involved in the beneficial effects of GT ingestion.

MicroRNAs as clinical tools for diagnosis, prognosis, and therapy in prostate cancer

Prostate cancer (PCa) is one of the most commonly diagnosed cancers among men worldwide. Despite the presence of accumulated clinical strategies for PCa management, limited prognostic/sensitive biomarkers are available to follow up on disease occurrence and progression. MicroRNAs (miRNAs) are small non-coding RNAs that control gene expression through post-transcriptional regulation of their complementary target messenger RNA (mRNA). MiRNAs modulate fundamental biological processes and play crucial roles in the pathology of various diseases, including PCa. Multiple evidence proved an aberrant miRNA expression profile in PCa, which is actively involved in the carcinogenic process. The robust and pleiotropic impact of miRNAs on PCa suggests them as potential candidates to help more understand the molecular landscape of the disease, which is likely to provide tools for early diagnosis and prognosis as well as additional therapeutic strategies to manage prostate tumors. Here, we emphasize the most consistently reported dysregulated miRNAs and highlight the contribution of their altered downstream targets with PCa hallmarks. Also, we report the potential effectiveness of using miRNAs as diagnostic/prognostic biomarkers in PCa and the high-throughput profiling technologies that are being used in their detection. Another key aspect to be discussed in this review is the promising implication of miRNAs molecules as therapeutic tools and targets for fighting PCa.

Characterization and microRNA Expression Analysis of Serum-Derived Extracellular Vesicles in Severe Liver Injury from Chronic HBV Infection

BACKGROUND

Extracellular vesicle (EV) microRNAs have been documented in several studies to have significantly different expressions in hepatitis B virus (HBV)-related liver diseases, such as hepatocellular carcinoma (HCC). The current work aimed to observe the characteristics of EVs and EV miRNA expressions in patients with severe liver injury chronic hepatitis B (CHB) and patients with HBV-associated decompensated cirrhosis (DeCi).

METHODS

The characterization of the EVs in the serum was carried out for three different groups, namely, patients with severe liver injury-CHB, patients with DeCi, and healthy controls. EV miRNAs were analyzed using miRNA-seq and RT-qPCR arrays. Additionally, we assessed the predictive and observational values of the miRNAs with significant differential expressions in serum EVs.

RESULTS

Patients with severe liver injury-CHB had the highest EV concentrations when compared to the normal controls (NCs) and patients with DeCi (p < 0.001). The miRNA-seq of the NC and severe liver injury-CHB groups identified 268 differentially expressed miRNAs (|FC| > 2, p < 0.05). In this case, 15 miRNAs were verified using RT-qPCR, and it was found that novel-miR-172-5p and miR-1285-5p in the severe liver injury-CHB group showed marked downregulation in comparison to the NC group (p < 0.001). Furthermore, compared with the NC group, three EV miRNAs (novel-miR-172-5p, miR-1285-5p, and miR-335-5p) in the DeCi group showed various degrees of downregulated expression. However, when comparing the DeCi group with the severe liver injury-CHB group, only the expression of miR-335-5p in the DeCi group decreased significantly (p < 0.05). For the severe liver injury-CHB and DeCi groups, the addition of miR-335-5p improved the predictive accuracy of the serological levels, while miR-335-5p was significantly correlated with ALT, AST, AST/ALT, GGT, and AFP. Conclusions: The patients with severe liver injury-CHB had the highest number of EVs. The combination of novel-miR-172-5p and miR-1285-5p in serum EVs helped in predicting the progression of the NCs to severe liver injury-CHB, while the addition of EV miR-335-5p improved the serological accuracy of predicting the progression of severe liver injury-CHB to DeCi.

The hitchhikers' guide to RNA sequencing and functional analysis

DNA and RNA sequencing technologies have revolutionized biology and biomedical sciences, sequencing full genomes and transcriptomes at very high speeds and reasonably low costs. RNA sequencing (RNA-Seq) enables transcript identification and quantification, but once sequencing has concluded researchers can be easily overwhelmed with questions such as how to go from raw data to differential expression (DE), pathway analysis and interpretation. Several pipelines and procedures have been developed to this effect. Even though there is no unique way to perform RNA-Seq analysis, it usually follows these steps: 1) raw reads quality check, 2) alignment of reads to a reference genome, 3) aligned reads' summarization according to an annotation file, 4) DE analysis and 5) gene set analysis and/or functional enrichment analysis. Each step requires researchers to make decisions, and the wide variety of options and resulting large volumes of data often lead to interpretation challenges. There also seems to be insufficient guidance on how best to obtain relevant information and derive actionable knowledge from transcription experiments. In this paper, we explain RNA-Seq steps in detail and outline differences and similarities of different popular options, as well as advantages and disadvantages. We also discuss non-coding RNA analysis, multi-omics, meta-transcriptomics and the use of artificial intelligence methods complementing the arsenal of tools available to researchers. Lastly, we perform a complete analysis from raw reads to DE and functional enrichment analysis, visually illustrating how results are not absolute truths and how algorithmic decisions can greatly impact results and interpretation.

Exosomal microRNAs regulate the heat stress response in sea cucumber Apostichopus japonicus

Exosomes are small extracellular vesicles that contain nucleic acids such as microRNAs and may participate in important biological processes. We made the initial report of exosomes from sea cucumber Apostichopus japonicus, that were classically cup-shaped and had an average size of 74.65 nm, and identified specific exosome biomarkers (HSP70, TSG101, and CD9). We explored changes in the global expression of microRNAs in exosomes from the commercially important A. japonicus under normal conditions and heat-stressed conditions for 3 and 7 d. We found that heat stress increased exosome production and modified the expression profiles of the microRNAs that they contained. Novel_mir31, novel_mir132, novel_mir26, miR-92_1, and novel_mir27 were commonly found to be differentially expressed in three comparison groups, indicating their importance in the heat stress response. The microRNA expression levels were validated by qPCR. Function analysis of the target genes of these microRNAs indicated they were involved mainly in replication and repair in the initial response of A. japonicus to heat stress exposure. Conversely, during acclimation to the high temperature conditions, the target genes of the differentially expressed microRNAs were primarily involved in metabolism adjustments. Our results will contribute to a better understanding of the regulatory roles of exosomes in sea cucumber, and provide insights into the functions of sea cucumber exosome-shuttled microRNAs against environmental stresses exacerbated by global warming.

Identification and profiling of Sogatella furcifera microRNAs and their potential roles in regulating the developmental transitions of nymph-adult

Sogatella furcifera is one of the most serious insect pests that affect rice in Asia. One class of small RNAs (sRNAs; ~22 nt long) is miRNAs, which participate in various biological processes by regulating the expression of target genes in a spatiotemporal manner. However, the role of miRNAs in nymph-to-adult transition in S. furcifera remains unknown. In this study, we sequenced sRNA libraries of S. furcifera prepared from individuals at three different developmental stages (pre-moult, moulting and early adult). A total of 253 miRNAs (134 known and 119 novel) were identified, of which 12 were differentially expressed during the nymph-to-adult developmental transition. Moreover, Real time quantitative PCR (RT-qPCR) analysis revealed that all 12 miRNAs were differentially expressed among five different nymph tissues and 14 different developmental stages (first to fifth instar nymphs and 1-day-old adults). Injection of miR-2a-2 mimic/antagomir and miR-305-5p-1 mimic/antagomir into 1-day-old fifth instar nymphs significantly increased the mortality rate. In addition, a defective moulting phenotype was observed in nymphs injected with miR-2a-2 and miR-305-5p-1, suggesting that these miRNAs are involved in S. furcifera nymph-adult transition. In conclusion, these results reveal the function of critical miRNAs in S. furcifera nymph-adult transition, and also provide novel potential targets of insecticides for the long-term sustainable management of S. furcifera.

MicroRNAs are deeply linked to the emergence of the complex octopus brain

Soft-bodied cephalopods such as octopuses are exceptionally intelligent invertebrates with a highly complex nervous system that evolved independently from vertebrates. Because of elevated RNA editing in their nervous tissues, we hypothesized that RNA regulation may play a major role in the cognitive success of this group. We thus profiled messenger RNAs and small RNAs in three cephalopod species including 18 tissues of the Octopus vulgaris. We show that the major RNA innovation of soft-bodied cephalopods is an expansion of the microRNA (miRNA) gene repertoire. These evolutionarily novel miRNAs were primarily expressed in adult neuronal tissues and during the development and had conserved and thus likely functional target sites. The only comparable miRNA expansions happened, notably, in vertebrates. Thus, we propose that miRNAs are intimately linked to the evolution of complex animal brains.

BrumiR: A toolkit for de novo discovery of microRNAs from sRNA-seq data

MicroRNAs (miRNAs) are small noncoding RNAs that are key players in the regulation of gene expression. In the past decade, with the increasing accessibility of high-throughput sequencing technologies, different methods have been developed to identify miRNAs, most of which rely on preexisting reference genomes. However, when a reference genome is absent or is not of high quality, such identification becomes more difficult. In this context, we developed BrumiR, an algorithm that is able to discover miRNAs directly and exclusively from small RNA (sRNA) sequencing (sRNA-seq) data. We benchmarked BrumiR with datasets encompassing animal and plant species using real and simulated sRNA-seq experiments. The results demonstrate that BrumiR reaches the highest recall for miRNA discovery, while at the same time being much faster and more efficient than the state-of-the-art tools evaluated. The latter allows BrumiR to analyze a large number of sRNA-seq experiments, from plants or animal species. Moreover, BrumiR detects additional information regarding other expressed sequences (sRNAs, isomiRs, etc.), thus maximizing the biological insight gained from sRNA-seq experiments. Additionally, when a reference genome is available, BrumiR provides a new mapping tool (BrumiR2reference) that performs an a posteriori exhaustive search to identify the precursor sequences. Finally, we also provide a machine learning classifier based on a random forest model that evaluates the sequence-derived features to further refine the prediction obtained from the BrumiR-core. The code of BrumiR and all the algorithms that compose the BrumiR toolkit are freely available at https://github.com/camoragaq/BrumiR.