Search

Simple Summary: The emergence of next-generation sequencing (NGS) signified a revolution in biology research by capturing the significance of DNA loci and RNA molecules on a genome-wide scale. The complexity and volume of NGS data highlight the need for robust and user-friendly computational tools. The framework presented here, aPEAch, is an automated pipeline for end-to-end analysis of DNA and RNA sequencing assays, including small RNA sequencing. Implemented in Python, it allows users to customize the analysis and results, handling single or multiple replicates in batches, while also automating advanced unsupervised learning analyses. aPEAch offers quality control reports, fragment size distribution plots, and intermediate files supporting reproducibility and interoperability, along with publication-ready visualizations. With the advent of next-generation sequencing (NGS), experimental techniques that capture the biological significance of DNA loci or RNA molecules have emerged as fundamental tools for studying the epigenome and transcriptional regulation on a genome-wide scale. The volume of the generated data and the underlying complexity regarding their analysis highlight the need for robust and easy-to-use computational analytic methods that can streamline the process and provide valuable biological insights. Our solution, aPEAch, is an automated pipeline that facilitates the end-to-end analysis of both DNA- and RNA-sequencing assays, including small RNA sequencing, from assessing the quality of the input sample files to answering meaningful biological questions by exploiting the rich information embedded in biological data. Our method is implemented in Python, based on a modular approach that enables users to choose the path and extent of the analysis and the representations of the results. The pipeline can process samples with single or multiple replicates in batches, allowing the ease of use and reproducibility of the analysis across all samples. aPEAch provides a variety of sample metrics such as quality control reports, fragment size distribution plots, and all intermediate output files, enabling the pipeline to be re-executed with different parameters or algorithms, along with the publication-ready visualization of the results. Furthermore, aPEAch seamlessly incorporates advanced unsupervised learning analyses by automating clustering optimization and visualization, thus providing invaluable insight into the underlying biological mechanisms. [ABSTRACT FROM AUTHOR]

Anthocyanins and proanthocyanidins are considered to be essential secondary metabolites in grapes and are used to regulate metabolic processes, while miRNAs are involved in their synthesis of anthocyanins and proanthocyanidins to regulate metabolic processes. The present research work was carried out to investigate the underlying regulatory mechanism of target genes in the grape cultivars 'Italia' and 'Benitaka'. miRNA and transnscriptomic sequencing technology were employed to characterize both the profiles of miRNAs and the transcripts of grape peels at 10 and 11 weeks post flowering (10 wpf and 11 wpf). The results revealed that the expression level of vvi-miR828a in 'Italia' at 10 and 11 wpf was significantly higher than that in 'Benitaka'. miRNA-seq analysis predicted MYBPA1 to be the target gene of vvi-miR828a. In transcriptome analysis, the expression level of the VvMYBPA1 gene in 'Benitaka' was significantly higher than that in 'Italia'; in addition, the TPM values (expression levels) of VvMYBPA1 and miR828a also showed an evident negative correlation. The determination of the proanthocyanidin (PA) content in 'Italia' and 'Benitaka' peels at 11 wpf demonstrated that the PA content of 'Benitaka' was significantly higher than that of 'Italia'. The outcomes of RT-qRCR analysis exhibited that the expression levels of the VdPAL, VdCHS, VdCHI, VdDFR, VdMYB5b, VdANR, and VdMYBPA1 genes related anthocyanin and proanthocyanidin pathways were reduced, while the expression levels of all of the above genes were increased after the transient expression of the VvMYBPA1 vector into grape leaves. The results of the transient overexpression experiment of vvi-miR828a before the veraison period of strawberry fruits showed that vvi-miR828a can significantly slow down the coloration of strawberries. The vvi-miR828a negatively regulates the accumulation of proanthocyanidins in grape fruits by inhibiting the expression of VvMYBPA1. [ABSTRACT FROM AUTHOR]

Obesity is one of the main causes of chronic kidney disease; however, the precise molecular mechanisms leading to the onset of kidney injury and dysfunction in obesity-associated nephropathy remain unclear. The present study aimed to unveil the kidney microRNA (miRNA) expression profile in a model of obesity-induced kidney disease in C57BL/6J mice using next-generation sequencing (NGS) analysis. High-fat diet (HFD)-induced obesity led to notable structural alterations in tubular and glomerular regions of the kidney, increased renal expression of proinflammatory and profibrotic genes, as well as an elevated renal expression of genes involved in cellular lipid metabolism. The miRNA sequencing analysis identified a set of nine miRNAs differentially expressed in the kidney upon HFD feeding, with miR-5099, miR-551b-3p, miR-223-3p, miR-146a-3p and miR-21a-3p showing the most significant differential expression between standard diet (STD) and HFD mice. A validation analysis showed that the expression levels of miR-5099, miR-551b-3p and miR-146a-3p were consistent with NGS results, while Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses revealed that these three validated miRNAs modulated target genes involved in metabolic and adipocytokine pathways, fatty acid and lipid metabolism, and inflammatory, senescence and profibrotic pathways. Our results suggest that differentially expressed miRNAs play pivotal roles in the intricate pathophysiology of obesity-associated kidney disease and could potentially create novel treatment strategies to counteract the deleterious effects of obesity on kidney function. [ABSTRACT FROM AUTHOR]

Breast milk, an indispensable source of immunological and nutrient components, is essential for the growth and development of newborn mammals. MicroRNAs (miRNAs) are present in various tissues and body fluids and are selectively packaged inside exosomes, a type of membrane vesicle. Milk exosomes have potential regulatory effects on the growth, development, and immunity of newborn piglets. To explore the differences in milk exosomes related to the breed and milk type, we isolated exosomes from colostrum and mature milk from domestic Bamei pigs and foreign Landrace pigs by using density gradient centrifugation and then characterized them by transmission electron microscopy (TEM) and nanoparticle tracking analysis (NTA). Furthermore, the profiles and functions of miRNAs in the two types of pig milk exosomes were investigated using miRNA-seq and bioinformatics analysis. We identified a total of 1081 known and 2311 novel miRNAs in pig milk exosomes from Bamei and Landrace pigs. These differentially expressed miRNAs (DE-miRNAs) are closely associated with processes such as cell signaling, cell physiology, and immune system development. Functional enrichment analysis showed that DE-miRNA target genes were significantly enriched in endocytosis, the T cell receptor signaling pathway, and the Th17 cell differentiation signaling pathway. The exosomal miRNAs in both the colostrum and mature milk of the two pig species showed significant differences. Based on related signaling pathways, we found that the colostrum of local pig breeds contained more immune-system-development-related miRNAs. This study provides new insights into the possible function of milk exosomal miRNAs in the development of the piglet immune system. [ABSTRACT FROM AUTHOR]

Background: Vitamin D is an immunomodulator, and its effects have been linked to many diseases, including the pathogenesis of cancer. However, the effect of vitamin D supplementation on the regulation of gene expression of the lungs is not fully understood. This study aims to determine the effect of the increased dose of cholecalciferol and a combination of cholecalciferol + calcidiol, as well as the replacement of cholecalciferol with calcidiol, on the miRNA profile of healthy swine lungs. Methods and results: The swine were long-term (88 days) supplemented with a standard dose (2000IU/kg) of cholecalciferol and calcidiol, the increased dose (3000 IU/kg) of cholecalciferol, and the cholecalciferol + calcidiol combination: grower: 3000 IU/Kg of vitamin D (67% of cholecalciferol and 33% of calcidiol), finisher 2500 IU/Kg of vitamin D (60% of cholecalciferol and 40% of calcidiol). Swine lung tissue was used for Next Generation Sequencing (NGS) of miRNA. Long-term supplementation with the cholecalciferol + calcidiol combination caused significant changes in the miRNA profile. They embraced altered levels of the expression of miR-150, miR-193, miR-145, miR-574, miR-340, miR-381, miR-148 and miR-96 (q-value < 0.05). In contrast, raising the dose of cholecalciferol only changed the expression of miR-215, and the total replacement of cholecalciferol with calcidiol did not significantly affect the miRNAome profile. Conclusions: The functional analysis of differentially expressed miRNAs suggests that the use of the increased dose of the cholecalciferol + calcidiol combination may affect tumorigenesis processes through, inter alia, modulation of gene regulation of the TGF- β pathway and pathways related to metabolism and synthesis of glycan. [ABSTRACT FROM AUTHOR]

Background: Vitamin D is an immunomodulator, and its effects have been linked to many diseases, including the pathogenesis of cancer. However, the effect of vitamin D supplementation on the regulation of gene expression of the lungs is not fully understood. This study aims to determine the effect of the increased dose of cholecalciferol and a combination of cholecalciferol + calcidiol, as well as the replacement of cholecalciferol with calcidiol, on the miRNA profile of healthy swine lungs. Methods and results: The swine were long-term (88 days) supplemented with a standard dose (2000IU/kg) of cholecalciferol and calcidiol, the increased dose (3000 IU/kg) of cholecalciferol, and the cholecalciferol + calcidiol combination: grower: 3000 IU/Kg of vitamin D (67% of cholecalciferol and 33% of calcidiol), finisher 2500 IU/Kg of vitamin D (60% of cholecalciferol and 40% of calcidiol). Swine lung tissue was used for Next Generation Sequencing (NGS) of miRNA. Long-term supplementation with the cholecalciferol + calcidiol combination caused significant changes in the miRNA profile. They embraced altered levels of the expression of miR-150, miR-193, miR-145, miR-574, miR-340, miR-381, miR-148 and miR-96 (q-value < 0.05). In contrast, raising the dose of cholecalciferol only changed the expression of miR-215, and the total replacement of cholecalciferol with calcidiol did not significantly affect the miRNAome profile. Conclusions: The functional analysis of differentially expressed miRNAs suggests that the use of the increased dose of the cholecalciferol + calcidiol combination may affect tumorigenesis processes through, inter alia, modulation of gene regulation of the TGF- β pathway and pathways related to metabolism and synthesis of glycan. [ABSTRACT FROM AUTHOR]

Abstract Background Tan and Hu sheep are well-known local breeds in China, producing lamb fur with unique ornamental and practical values highly appreciated by consumers worldwide. Fur quality is optimal at one month of age and gradually declines with time. Despite active research on its genetic mechanism using transcriptomic and whole genome bisulfite sequencing analysis, the main effective gene locus has not been found, and its regulatory mechanism is still unclear, which limits the breeding and improvement of fur traits. Results Scapular skin samples from newborn (1-month old) and adult (24-month old) Tan sheep were utilized for small ribonucleic acid (RNA) sequencing Principal Component Analysis (PCA) showed that the newborn and adult groups were completely separated. Differential expression analysis of micro-RNAs (miRNAs) identified 32 up-regulated miRNAs and 48 down-regulated miRNAs in the newborn groups. All up-regulated miRNAs were located in the imprinted. Dlk1-Gtl2 locus on chromosome 18, whereas all down-regulated miRNAs were distributed across the sheep chromosomes, without a clear pattern of positional consistency. Further, by systematically analyzing the target genes and signaling pathways of all 32 up-regulated miRNAs, we found that the PI3K-AKT signaling pathway has the potential to be targeted and regulated by most of the miRNAs in the Dlk1-Gtl2 region. In addition, we also re-analyzed miRNA sequencing data from public databases on Hu lambs (full sibling Hu lambs with high- and low-quality fur characteristics). Again, it was found that most of the up-regulated miRNAs in lambs with high-quality fur were also located in the Dlk1-Gtl2 region, whereas this patter was not present for down-regulated miRNAs. Conclusion Sequencing of miRNAs in conjunction with public databases was employed to identify miRNAs within the imprinted Dlk1-Gtl2 region on chromosome 18, suggesting their potential roles as epigenetic regulators of fur traits. Small RNAs located at the Dlk1-Gtl2 locus were identified as having the potential to systematically regulate the PI3K-AKT signaling pathway, thereby indicating the relevance of the Dlk1-Gtl2/PI3K-AKT axis in the context of fur traits. Selection of parental specific expressed imprinted genes in the process of conserving and exploiting lamb fur traits should be emphasized.

With the advent of next-generation sequencing (NGS), experimental techniques that capture the biological significance of DNA loci or RNA molecules have emerged as fundamental tools for studying the epigenome and transcriptional regulation on a genome-wide scale. The volume of the generated data and the underlying complexity regarding their analysis highlight the need for robust and easy-to-use computational analytic methods that can streamline the process and provide valuable biological insights. Our solution, aPEAch, is an automated pipeline that facilitates the end-to-end analysis of both DNA- and RNA-sequencing assays, including small RNA sequencing, from assessing the quality of the input sample files to answering meaningful biological questions by exploiting the rich information embedded in biological data. Our method is implemented in Python, based on a modular approach that enables users to choose the path and extent of the analysis and the representations of the results. The pipeline can process samples with single or multiple replicates in batches, allowing the ease of use and reproducibility of the analysis across all samples. aPEAch provides a variety of sample metrics such as quality control reports, fragment size distribution plots, and all intermediate output files, enabling the pipeline to be re-executed with different parameters or algorithms, along with the publication-ready visualization of the results. Furthermore, aPEAch seamlessly incorporates advanced unsupervised learning analyses by automating clustering optimization and visualization, thus providing invaluable insight into the underlying biological mechanisms.

Anthocyanins and proanthocyanidins are considered to be essential secondary metabolites in grapes and are used to regulate metabolic processes, while miRNAs are involved in their synthesis of anthocyanins and proanthocyanidins to regulate metabolic processes. The present research work was carried out to investigate the underlying regulatory mechanism of target genes in the grape cultivars ‘Italia’ and ‘Benitaka’. miRNA and transnscriptomic sequencing technology were employed to characterize both the profiles of miRNAs and the transcripts of grape peels at 10 and 11 weeks post flowering (10 wpf and 11 wpf). The results revealed that the expression level of vvi-miR828a in ‘Italia’ at 10 and 11 wpf was significantly higher than that in ‘Benitaka’. miRNA-seq analysis predicted MYBPA1 to be the target gene of vvi-miR828a. In transcriptome analysis, the expression level of the VvMYBPA1 gene in ‘Benitaka’ was significantly higher than that in ‘Italia’; in addition, the TPM values (expression levels) of VvMYBPA1 and miR828a also showed an evident negative correlation. The determination of the proanthocyanidin (PA) content in ‘Italia’ and ‘Benitaka’ peels at 11 wpf demonstrated that the PA content of ‘Benitaka’ was significantly higher than that of ‘Italia’. The outcomes of RT-qRCR analysis exhibited that the expression levels of the VdPAL, VdCHS, VdCHI, VdDFR, VdMYB5b, VdANR, and VdMYBPA1 genes related anthocyanin and proanthocyanidin pathways were reduced, while the expression levels of all of the above genes were increased after the transient expression of the VvMYBPA1 vector into grape leaves. The results of the transient overexpression experiment of vvi-miR828a before the veraison period of strawberry fruits showed that vvi-miR828a can significantly slow down the coloration of strawberries. The vvi-miR828a negatively regulates the accumulation of proanthocyanidins in grape fruits by inhibiting the expression of VvMYBPA1.

Testis size is important for identifying breeding animals with adequate sperm production. The aim of this study was to survey the expression profile of mRNA and miRNA in testis tissue from rams carrying different FecB genotypes, including the wild-type and heterozygous genotypes in Tibetan sheep. Comparative transcriptome profiles for ovine testes were established for wild-type and heterozygote Tibetan sheep by next-generation sequencing. RNA-seq results identified 3,910 (2,034 up- and 1,876 downregulated) differentially expressed (DE) genes and 243 (158 up- and 85 downregulated) DE microRNAs (miRNAs) in wild-type vs heterozygote sheep, respectively. Combined analysis of mRNA-seq and miRNA-seq revealed that 20 miRNAs interacted with 48 true DE target genes in wild-type testes compared to heterozygous genotype testes. These results provide evidence for a functional series of genes operating in Tibetan sheep testis. In addition, quantitative real-time PCR analysis showed that the expression trends of randomly selected DE genes in testis tissues from different genotypes were consistent with high-throughput sequencing results.

Background: Tan and Hu sheep are well-known local breeds in China, producing lamb fur with unique ornamental and practical values highly appreciated by consumers worldwide. Fur quality is optimal at one month of age and gradually declines with time. Despite active research on its genetic mechanism using transcriptomic and whole genome bisulfite sequencing analysis, the main effective gene locus has not been found, and its regulatory mechanism is still unclear, which limits the breeding and improvement of fur traits. Results: Scapular skin samples from newborn (1-month old) and adult (24-month old) Tan sheep were utilized for small ribonucleic acid (RNA) sequencing Principal Component Analysis (PCA) showed that the newborn and adult groups were completely separated. Differential expression analysis of micro-RNAs (miRNAs) identified 32 up-regulated miRNAs and 48 down-regulated miRNAs in the newborn groups. All up-regulated miRNAs were located in the imprinted. Dlk1-Gtl2 locus on chromosome 18, whereas all down-regulated miRNAs were distributed across the sheep chromosomes, without a clear pattern of positional consistency. Further, by systematically analyzing the target genes and signaling pathways of all 32 up-regulated miRNAs, we found that the PI3K-AKT signaling pathway has the potential to be targeted and regulated by most of the miRNAs in the Dlk1-Gtl2 region. In addition, we also re-analyzed miRNA sequencing data from public databases on Hu lambs (full sibling Hu lambs with high- and low-quality fur characteristics). Again, it was found that most of the up-regulated miRNAs in lambs with high-quality fur were also located in the Dlk1-Gtl2 region, whereas this patter was not present for down-regulated miRNAs. Conclusion: Sequencing of miRNAs in conjunction with public databases was employed to identify miRNAs within the imprinted Dlk1-Gtl2 region on chromosome 18, suggesting their potential roles as epigenetic regulators of fur traits. Small RNAs located at the Dlk1-Gtl2 locus were identified as having the potential to systematically regulate the PI3K-AKT signaling pathway, thereby indicating the relevance of the Dlk1-Gtl2/PI3K-AKT axis in the context of fur traits. Selection of parental specific expressed imprinted genes in the process of conserving and exploiting lamb fur traits should be emphasized. [ABSTRACT FROM AUTHOR]

Background: Macrophages play an important role in regulating the course of spinal tuberculosis within the bone microenvironment. This study aimed to investigate the differential expression of miRNA in macrophage-derived exosomes within the tuberculosis-infected bone microenvironment, to identify miRNAs that hold potential as diagnostic markers and therapeutic targets. Methods: We established study cohorts for spinal tuberculosis, collected bone marrow blood samples, isolated macrophage exosomes, and performed exosome miRNA sequencing. A miRNA-mRNA co-expression network was constructed using WGCNA analysis. Gene GO analysis and KEGG pathway enrichment analysis were performed using KOBAS software. Target miRNAs were selected based on fold change, P-value, and false discovery rate, and their validation was carried out using qRT-PCR and ROC curve studies. Subsequently, we constructed a target gene network for these miRNAs and performed KEGG pathway enrichment analysis to explore the potential signaling mechanisms involved in regulating the disease course of spinal tuberculosis. Results: Our findings revealed that macrophages from the tuberculosis-infected bone microenvironment exhibited an M1 phenotype. The successful extraction of exosomes from macrophage supernatants was confirmed through electron microscopy, particle size analysis, and protein blot analysis. Exosome miRNAseq demonstrated that 28 miRNAs were up-regulated, while 34 miRNAs were down-regulated in individuals with spinal tuberculosis. GO analysis and KEGG pathway enrichment analysis indicated that the differentially expressed miRNAs were involved in various biological processes, cell components, molecular functions, and signaling pathways, which collectively contribute to the regulation of the disease course of spinal tuberculosis. Notably, miRNA-125b-5p was successfully selected based on fold change, p-value, and false discovery rate. qRT-PCR validation further confirmed the significant up-regulation of miRNA-125b-5p in spinal tuberculosis. The ROC curve revealed that miR125b-5p is a potential diagnostic biomarker for spinal tuberculosis. Moreover, construction of the miRNA-125b-5p target gene network and subsequent KEGG enrichment analysis highlighted the importance of MAPK, TNF, Ras, Rap1, and the PI3K-Akt signaling pathways in the regulation of the disease course of spinal tuberculosis. Conclusion: Our study demonstrates differential expression of miRNAs in macrophage-derived exosomes in the tuberculosis-infected bone microenvironment. Specifically, MiRNA-125b-5p is significantly up-regulated in spinal tuberculosis and shows potential as a diagnostic biomarker for spinal tuberculosis. [ABSTRACT FROM AUTHOR]

Abstract Background Despite significant advancements in multiple myeloma (MM) therapy, the highly heterogenous treatment response hinders reliable prognosis and tailored therapeutics. Herein, we have studied the clinical utility of miRNAs in ameliorating patients’ management. Methods miRNA-seq was performed in bone marrow CD138+ plasma cells (PCs) of 24 MM and smoldering MM (sMM) patients to analyze miRNAs profile. CD138+ and circulating miR-25 levels were quantified using in house RT-qPCR assays in our screening MM/sMM cohort (CD138+ plasma cells n = 167; subcohort of MM peripheral plasma samples n = 69). Two external datasets (Kryukov et al. cohort n = 149; MMRF CoMMpass study n = 760) served as institutional-independent validation cohorts. Patients’ mortality and disease progression were assessed as clinical endpoints. Internal validation was performed by bootstrap analysis. Clinical benefit was estimated by decision curve analysis. Results miRNA-seq highlighted miR-25 of CD138+ plasma cells to be upregulated in MM vs. sMM, R-ISS II/III vs. R-ISS I, and in progressed compared to progression-free patients. The analysis of our screening cohort highlighted that CD138+ miR-25 levels were correlated with short-term progression (HR = 2.729; p = 0.009) and poor survival (HR = 4.581; p = 0.004) of the patients; which was confirmed by Kryukov et al. cohort (HR = 1.878; p = 0.005) and MMRF CoMMpass study (HR = 1.414; p = 0.039) validation cohorts. Moreover, multivariate miR-25-fitted models contributed to superior risk-stratification and clinical benefit in MM prognostication. Finally, elevated miR-25 circulating levels were correlated with poor survival of MM patients (HR = 5.435; p = 0.021), serving as a potent non-invasive molecular prognostic tool. Conclusions Our study identified miR-25 overexpression as a powerful independent predictor of poor treatment outcome and post-treatment progression, aiding towards modern non-invasive disease prognosis and personalized treatment decisions.

Recently, a novel poplar mosaic disease caused by bean common mosaic virus (BCMV) was investigated in Populus alba var. pyramidalis in China. Symptom characteristics, physiological performance of the host, histopathology, genome sequences and vectors, and gene regulation at the transcriptional and posttranscriptional levels were analyzed and RT-qPCR (quantitative reverse transcription PCR) validation of expression was performed in our experiments. In this work, the mechanisms by which the BCMV pathogen impacts physiological performance and the molecular mechanisms of the poplar response to viral infection were reported. The results showed that BCMV infection decreased the chlorophyll content, inhibited the net photosynthesis rate (Pn) and stomatal conductance (Gs), and significantly changed chlorophyll fluorescence parameters in diseased leaves. Transcriptome analysis revealed that the expression of the majority of DEGs (differentially expressed genes) involved in the flavonoid biosynthesis pathway was promoted, but the expression of all or almost all DEGs associated with photosynthesis-antenna proteins and the photosynthesis pathway was inhibited in poplar leaves, suggesting that BCMV infection increased the accumulation of flavonoids but decreased photosynthesis in hosts. Gene set enrichment analysis (GSEA) illustrated that viral infection promoted the expression of genes involved in the defense response or plant-pathogen interaction. MicroRNA-seq analysis illustrated that 10 miRNA families were upregulated while 6 families were downregulated in diseased poplar leaves; moreover, miR156, the largest family with the most miRNA members and target genes, was only differentially upregulated in longperiod disease (LD) poplar leaves. Integrated transcriptome and miRNA-seq analyses revealed 29 and 145 candidate miRNA-target gene pairs; however, only 17 and 76 pairs, accounting for 2.2% and 3.2% of all DEGs, were authentically negatively regulated in short-period disease (SD) and LD leaves, respectively. Interestingly, 4 miR156/SPL (squamosa promoter-binding-like protein) miRNA -target gene pairs were identified in LD leaves: the miR156 molecules were upregulated, but SPL genes were downregulated. In conclusion, BCMV infection significantly changed transcriptional and posttranscriptional gene expression in poplar leaves, inhibited photosynthesis, increased the accumulation of flavonoids, induced systematic mosaic symptoms, and decreased physiological performance in diseased poplar leaves. This study elucidated the fine-tuned regulation of poplar gene expression by BCMV; moreover, the results also suggested that miR156/SPL modules played important roles in the virus response and development of viral systematic symptoms in plant virus disease. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

MicroRNA transcriptomes from fresh tumors and the adjacent normal tissues were profiled in 10 Korean patients diagnosed with lung adenocarcinoma using a next-generation sequencing (NGS) technique called miRNA-seq. The sequencing quality was assessed using FastQC, and low-quality or adapter-contaminated portions of the reads were removed using Trim Galore. Quality-assured reads were analyzed using miRDeep2 and Bowtie. The abundance of known miRNAs was estimated using the reads per million (RPM) normalization method. Subsequently, using DESeq2 and Wx, we identified differentially expressed miRNAs and potential miRNA biomarkers for lung adenocarcinoma tissues compared to adjacent normal tissues, respectively. We defined reliable miRNA biomarkers for lung adenocarcinoma as those detected by both methods. The miRNA-seq data are available in the Gene Expression Omnibus (GEO) database under accession number GSE196633, and all processed data can be accessed via the Mendeley data website. Dataset:  https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE196633 and https://data.mendeley.com/datasets/vp977psjcb/2. Dataset License: CC0 [ABSTRACT FROM AUTHOR]

Exosomes cargo tumour-characterized biomolecules secreted from cancer cells and play a pivotal role in tumorigenesis and cancer progression, thus providing their potential for non-invasive cancer monitoring. Since cancer cell-derived exosomes are often mixed with those from healthy cells in liquid biopsy of tumour patients, accurately measuring the purity of tumour cell-derived exosomes is not only critical for the early detection but also essential for unbiased identification of diagnosis biomarkers. Here, we propose 'ExosomePurity', a tumour purity deconvolution model to estimate tumour purity in serum exosomes of cancer patients based on microribonucleic acid (miRNA)-Seq data. We first identify the differently expressed miRNAs as signature to distinguish cancer cell- from healthy cell-derived exosomes. Then, the deconvolution model was developed to estimate the proportions of cancer exosomes and normal exosomes in serum. The purity predicted by the model shows high correlation with actual purity in simulated data and actual data. Moreover, the model is robust under the different levels of noise background. The tumour purity was also used to correct differential expressed gene analysis. ExosomePurity empowers the research community to study non-invasive early diagnosis and to track cancer progression in cancers more efficiently. It is implemented in R and is freely available from GitHub (https://github.com/WangHYLab/ExosomePurity). [ABSTRACT FROM AUTHOR]

Recently, we reported that forkhead box A2 (FOXA2) is required for the development of human pancreatic α- and β-cells. However, whether miRNAs play a role in regulating pancreatic genes during pancreatic development in the absence of FOXA2 expression is largely unknown. Here, we aimed to capture the dysregulated miRNAs and to identify their pancreatic-specific gene targets in pancreatic progenitors (PPs) derived from wild-type induced pluripotent stem cells (WT-iPSCs) and from iPSCs lacking FOXA2 (FOXA2–/–iPSCs). To identify differentially expressed miRNAs (DEmiRs), and genes (DEGs), two different FOXA2–/–iPSC lines were differentiated into PPs. FOXA2–/– PPs showed a significant reduction in the expression of the main PP transcription factors (TFs) in comparison to WT-PPs. RNA sequencing analysis demonstrated significant reduction in the mRNA expression of genes involved in the development and function of exocrine and endocrine pancreas. Furthermore, miRNA profiling identified 107 downregulated and 111 upregulated DEmiRs in FOXA2–/– PPs compared to WT-PPs. Target prediction analysis between DEmiRs and DEGs identified 92 upregulated miRNAs, predicted to target 1498 downregulated genes in FOXA2–/– PPs. Several important pancreatic TFs essential for pancreatic development were targeted by multiple DEmiRs. Selected DEmiRs and DEGs were further validated using RT-qPCR. Our findings revealed that FOXA2 expression is crucial for pancreatic development through regulating the expression of pancreatic endocrine and exocrine genes targeted by a set of miRNAs at the pancreatic progenitor stage. These data provide novel insights of the effect of FOXA2 deficiency on miRNA-mRNA regulatory networks controlling pancreatic development and differentiation. [ABSTRACT FROM AUTHOR]

Objective Adipocyte differentiation is regulated by a variety of functional genes and non-coding RNAs. However, the role of miRNAs in lipid deposition of goat white adipose tissue is still unclear. Therefore, this study revealed the miRNA expression profile in goat subcutaneous adipocytes by sRNA-seq. Methods The miRNA expressed in goat subcutaneous preadipocytes and the mature adipocytes were sequenced by sRNA-seq. The differentially expressed miRNAs (DEm) were screened and gene ontology (GO) and Kyoto encyclopedia for genes and genomes (KEGG) analyses were performed. Gain-of-function and loss-of-function combined with oil red O staining, Bodipy staining, and quantitative reverse-transcription polymerase chain reaction (qPCR) were utilized to determine the effect of miR-133a-3p on adipocyte differentiation. Results A total of 218 DEm were screened out. The target genes of these DEm were significantly enriched in GO items such as biological regulation and in KEGG terms such as FAK signaling pathway and MAPK signaling pathway. qPCR verified that the expression trend of miRNA was consistent with miRNA-seq. The gain-of-function or loss-of-function of miR-133a-3p showed that it promoted or inhibited the accumulation of lipid droplets, and CCAAT enhancer binding protein α (C/EBPα) and C/EBPβ were extremely significantly up-regulated or down-regulated respectively (p

Background: Despite significant advancements in multiple myeloma (MM) therapy, the highly heterogenous treatment response hinders reliable prognosis and tailored therapeutics. Herein, we have studied the clinical utility of miRNAs in ameliorating patients' management. Methods: miRNA-seq was performed in bone marrow CD138+ plasma cells (PCs) of 24 MM and smoldering MM (sMM) patients to analyze miRNAs profile. CD138+ and circulating miR-25 levels were quantified using in house RT-qPCR assays in our screening MM/sMM cohort (CD138+ plasma cells n = 167; subcohort of MM peripheral plasma samples n = 69). Two external datasets (Kryukov et al. cohort n = 149; MMRF CoMMpass study n = 760) served as institutional-independent validation cohorts. Patients' mortality and disease progression were assessed as clinical endpoints. Internal validation was performed by bootstrap analysis. Clinical benefit was estimated by decision curve analysis. Results: miRNA-seq highlighted miR-25 of CD138+ plasma cells to be upregulated in MM vs. sMM, R-ISS II/III vs. R-ISS I, and in progressed compared to progression-free patients. The analysis of our screening cohort highlighted that CD138+ miR-25 levels were correlated with short-term progression (HR = 2.729; p = 0.009) and poor survival (HR = 4.581; p = 0.004) of the patients; which was confirmed by Kryukov et al. cohort (HR = 1.878; p = 0.005) and MMRF CoMMpass study (HR = 1.414; p = 0.039) validation cohorts. Moreover, multivariate miR-25-fitted models contributed to superior risk-stratification and clinical benefit in MM prognostication. Finally, elevated miR-25 circulating levels were correlated with poor survival of MM patients (HR = 5.435; p = 0.021), serving as a potent non-invasive molecular prognostic tool. Conclusions: Our study identified miR-25 overexpression as a powerful independent predictor of poor treatment outcome and post-treatment progression, aiding towards modern non-invasive disease prognosis and personalized treatment decisions. [ABSTRACT FROM AUTHOR]

Astrocytes are central nervous system (CNS)-restricted glial cells involved in synaptic function and CNS blood flow regulation. Astrocyte extracellular vesicles (EVs) participate in neuronal regulation. EVs carry RNAs, either surface-bound or luminal, which can be transferred to recipient cells. We characterized the secreted EVs and RNA cargo of human astrocytes derived from an adult brain. EVs were isolated by serial centrifugation and characterized with nanoparticle tracking analysis (NTA), Exoview, and immuno-transmission electron microscopy (TEM). RNA from cells, EVs, and proteinase K/RNase-treated EVs was analyzed by miRNA-seq. Human adult astrocyte EVs ranged in sizes from 50 to 200 nm, with CD81 as the main tetraspanin marker and larger EVs positive for integrin β1. Comparison of the RNA between the cells and EVs identified RNA preferentially secreted in the EVs. In the case of miRNAs, enrichment analysis of their mRNA targets indicates that they are good candidates for mediating EV effects on recipient cells. The most abundant cellular miRNAs were also abundant in EVs, and the majority of their mRNA targets were found to be downregulated in mRNA-seq data, but the enrichment analysis lacked neuronal specificity. Proteinase K/RNase treatment of EV-enriched preparations identified RNAs secreted independently of EVs. Comparing the distribution of cellular and secreted RNA identifies the RNAs involved in intercellular communication via EVs. [ABSTRACT FROM AUTHOR]

High-throughput sequencing of small RNA molecules such as microRNAs (miRNAs) has become a widely used approach for studying gene expression and regulation. However, analyzing miRNA-Seq data can be challenging because it requires multiple steps, from quality control and preprocessing to differential expression and pathway-enrichment analyses, with many tools and databases available for each step. Furthermore, reproducibility of the analysis pipeline is crucial to ensure that the results are accurate and reliable. Here, we present myBrain-Seq, a comprehensive and reproducible pipeline for analyzing miRNA-Seq data that incorporates miRNA-specific solutions at each step of the analysis. The pipeline was designed to be flexible and user-friendly, allowing researchers with different levels of expertise to perform the analysis in a standardized and reproducible manner, using the most common and widely used tools for each step. In this work, we describe the implementation of myBrain-Seq and demonstrate its capacity to consistently and reproducibly identify differentially expressed miRNAs and enriched pathways by applying it to a real case study in which we compared schizophrenia patients who responded to medication with treatment-resistant schizophrenia patients to obtain a 16-miRNA treatment-resistant schizophrenia profile. [ABSTRACT FROM AUTHOR]

Introduction: Atrial fibrillation (AF) is the most common cardiac arrhythmia in the general population, and stroke is the most severe complication of AF. Exosomal miRNAs have been reported to be candidates as biomarkers for cardiovascular diseases, including AF and stroke. This study aimed to identify differentially expressed miRNAs (DEMs) in serum exosomes of AF and AF-associated ischemic stroke (AF-IS) patients and evaluate their potential in distinguishing AF and AF-IS patients. Methods: Serum exosomes were isolated from 8 healthy individuals with sinus rhythm (SR controls), 8 AF patients, and 8 AF-IS patients. miRNA-seq was performed to identify DEMs, and qRT-PCR analysis was performed to confirm the sequencing results. A support vector machine (SVM) model was developed using Python to distinguish AF and AF-IS patients. Results: 68 and 86 DEMs were identified in serum exosomes of AF patients compared to AF-IS patients and SR controls, respectively. Levels of miR-641 and miR-30e-5p were found significantly higher in AF-IS patients. The SVM model achieved an accuracy of 100%, with an area under curve of 1. Conclusions: The results indicated that miRNA expression profiles of serum exosomes in AF patients were distinct from those in AF-IS patients, and based on the distinction, AF and AF-IS patients can be distinguished. [ABSTRACT FROM AUTHOR]

17α-Methyltestosterone (17MT), a synthetic organic compound commonly found in sewage waters, can affect reproduction in aquatic animals, such as tilapia and yellow catfish. In the present study, male Gobiocypris rarus were exposed to 25, 50, and 100 ng/L of 17α-methyltestosterone (17MT) for 7 days. We first analyzed miRNA- and RNA-seq results to determine miRNA-target gene pairs and then developed miRNA-mRNA interactive networks after 17MT administration. Total weights, total lengths, and body lengths were not significantly different between the test groups and control groups. The paraffin slice method was applied to testes of G. rarus in the MT exposure and control groups. We found that there were more mature sperm (S) and fewer secondary spermatocytes (SSs) and spermatogonia (SGs) in the testes of control groups. As 17MT concentration increased, fewer and fewer mature sperm (S) were observed in the testes of male G. rarus. The results showed that FSH, 11-KT, and E2 were significantly higher in individuals exposed to 25 ng/L 17MT compared with the control groups. VTG, FSH, LH, 11-KT, and E2 were significantly lower in the 50 ng/L 17MT exposure groups compared to the control groups. VTG, FSH, LH, 11-KT, E2, and T were significantly lower in the groups exposed to 100 ng/L 17MT. High-throughput sequencing revealed 73,449 unigenes, 1205 known mature miRNAs, and 939 novel miRNAs in the gonads of G. rarus. With miRNA-seq, 49 (MT25-M vs. Con-M), 66 (MT50-M vs. Con-M), and 49 (MT100-M vs. Con-M) DEMs were identified in the treatment groups. Five mature miRNAs (miR-122-x, miR-574-x, miR-430-y, lin-4-x, and miR-7-y), as well as seven differentially expressed genes (soat2, inhbb, ihhb, gatm, faxdc2, ebp, and cyp1a1), which may be associated with testicular development, metabolism, apoptosis, and disease response, were assayed using qRT-PCR. Furthermore, miR-122-x (related to lipid metabolism), miR-430-y (embryonic development), lin-4-x (apoptosis), and miR-7-y (disease) were differentially expressed in the testes of 17MT-exposed G. rarus. This study highlights the role of miRNA-mRNA pairs in the regulation of testicular development and immune response to disease and will facilitate future studies on the miRNA-RNA-associated regulation of teleost reproduction. [ABSTRACT FROM AUTHOR]

The complete replacement of fresh baits (FB) by formulated artificial feeds (AF) after domestication in mandarin fish (Siniperca chuatsi) has matured in recent years, while the underlying molecular events of long-term AF diet feeding after diet shifts remained unclear. Therefore, mandarin fish, which had been successfully domesticated, were fed with the FB or AF for 180 days to explore the molecular mechanisms underlying the differences due to shifted food habit. Compared to the FB group, slight intestinal histological damage was found in the AF group of mandarin fish. An integrated analysis of miRNA-mRNA sequencing and their interaction networks was conducted between the FB and AF groups. The result showed that a total of 209 differentially expressed mRNAs (DEmRNAs) and 2 differentially expressed miRNAs (DEmiRNAs) with significance were identified. The DEmRNAs were enriched to 39 pathways including tight junction, steroid biosynthesis and protein processing in ER. Oni-Let-7a and aca-Let-7i were the two significantly identified DEmiRNAs, whose target prediction was enriched with GO terms and KEGG pathways including TGFβ signaling pathway, steroid biosynthesis, and retinol metabolism. In general, our research was the first exploration to characterize the interacted miRNA-mRNA sequencing of mandarin fish fed the FB and AF diets, which could reveal the mechanisms of long-term adaptation to new foods in organisms who subjected to a food conversion in their early life. Data Availability Statement: The dataset generated by deep-sequencing platforms was deposited in the NCBI Sequence Read Archive (SRA, http://www.ncbi.nlm.nih.gov/Traces/sra). The NCBI Sequence Read Archive accession number is PRJNA897748.

Breast milk, an indispensable source of immunological and nutrient components, is essential for the growth and development of newborn mammals. MicroRNAs (miRNAs) are present in various tissues and body fluids and are selectively packaged inside exosomes, a type of membrane vesicle. Milk exosomes have potential regulatory effects on the growth, development, and immunity of newborn piglets. To explore the differences in milk exosomes related to the breed and milk type, we isolated exosomes from colostrum and mature milk from domestic Bamei pigs and foreign Landrace pigs by using density gradient centrifugation and then characterized them by transmission electron microscopy (TEM) and nanoparticle tracking analysis (NTA). Furthermore, the profiles and functions of miRNAs in the two types of pig milk exosomes were investigated using miRNA-seq and bioinformatics analysis. We identified a total of 1081 known and 2311 novel miRNAs in pig milk exosomes from Bamei and Landrace pigs. These differentially expressed miRNAs (DE-miRNAs) are closely associated with processes such as cell signaling, cell physiology, and immune system development. Functional enrichment analysis showed that DE-miRNA target genes were significantly enriched in endocytosis, the T cell receptor signaling pathway, and the Th17 cell differentiation signaling pathway. The exosomal miRNAs in both the colostrum and mature milk of the two pig species showed significant differences. Based on related signaling pathways, we found that the colostrum of local pig breeds contained more immune-system-development-related miRNAs. This study provides new insights into the possible function of milk exosomal miRNAs in the development of the piglet immune system.

Florfenicol (FLR), a broad-spectrum antibacterial agent, is commonly used in the breeding and cultivation of aquatic animals. microRNAs (miRNAs) control gene expression by mediating post-transcriptional mRNA repression. However, not much has been reported regarding their functions when crustaceans encounter FLR stress. Herein we observed histological alterations and performed small RNA sequencing to characterize differentially expressed miRNAs in the hepatopancreas of Portunus trituberculatus exposed to FLR (80 mg/kg) stress for 3 h and 12 h. Our results indicated irregularities in the tubular structure of hepatopancreas, separation between epithelium and myoepithelial layer, ruptured epithelial cells, and tubule degeneration. Thirty-five (19 up- and 16 downregulated) and 31 (17 up- and 14 downregulated) differentially expressed miRNAs were detected between Group-C and Group-FLR-3 h and between Group-C and Group-FLR-12 h, respectively. Dual luciferase reporter assay was performed for predicting target genes of miR-34 and miR-263b. We found that mir-263b directly inhibited cytochrome C expression at the cell level in vitro. Functional analyses of predicted target genes of differentially expressed miRNAs revealed that FLR stress induced xenobiotic metabolism, immune toxicity, and cell apoptosis, as well as influenced diverse metabolic processes. To our knowledge, this is the first report on the response of P. trituberculatus to FLR stress. Our findings provide evidence that miRNAs participate in cell apoptosis and immune regulation in crabs.

Intensive shrimp farming is often threatened by microcystins Hepatopancreas is the primary target organ of MCs in shrimp. To investigate the response of hepatopancreas to acute MC-LR exposure, the expression profiles of RNA-seq and miRNA-seq in the hepatopancreas of L. vannamei were determined, and data integration analysis was performed at 72 h after MC-LR injection. The expression of 5 DEGs and three DEMs were detected by Quantitative PCR (qPCR). The results showed that the cumulative mortality rate of shrimp in MC-LR treatment group was 41.1%. A total of 1229 differentially expressed genes (844 up- and 385 downregulated) and 86 differentially expressed miRNAs (40 up- and 46 downregulated) were identified after MC-LR exposure. Functional analysis indicated that DEGs is mainly involved in the oxidative activity process in molecular functional categories, and proteasome was the most enriched KEGG pathway for mRNAs profile. According to the functional annotation of target genes of DEMs, protein binding was the most important term in the GO category, and protein processing in endoplasmic reticulum (ER) was the most enriched KEGG pathway. The regulatory network of miRNAs and DEGs involved in the pathway related to protein degradation in endoplasmic reticulum was constructed, and miR-181-5p regulated many genes in this pathway. The results of qPCR showed that there were significant differences in the expression of five DEGs and three DEMs, which might play an important role in the toxicity and hepatopancreas detoxification of MC-LR in shrimp. The results revealed that MC-LR exposure affected the degradation pathway of misfolded protein in ER of L. vannamei hepatopancreas, and miR-181-5p might play an important role in the effect of MC-LR on the degradation pathway of misfolded protein. [ABSTRACT FROM AUTHOR]

Testis size is important for identifying breeding animals with adequate sperm production. The aim of this study was to survey the expression profile of mRNA and miRNA in testis tissue from rams carrying different FecB genotypes, including the wild-type and heterozygous genotypes in Tibetan sheep. Comparative transcriptome profiles for ovine testes were established for wild-type and heterozygote Tibetan sheep by next-generation sequencing. RNA-seq results identified 3,910 (2,034 up- and 1,876 downregulated) differentially expressed (DE) genes and 243 (158 up- and 85 downregulated) DE microRNAs (miRNAs) in wild-type vs heterozygote sheep, respectively. Combined analysis of mRNA-seq and miRNA-seq revealed that 20 miRNAs interacted with 48 true DE target genes in wildtype testes compared to heterozygous genotype testes. These results provide evidence for a functional series of genes operating in Tibetan sheep testis. In addition, quantitative real-time PCR analysis showed that the expression trends of randomly selected DE genes in testis tissues from different genotypes were consistent with high-throughput sequencing results. [ABSTRACT FROM AUTHOR]

• A total of 51 known miRNAs and 95 novel miRNAs were identified from the nine small RNA libraries. • 6 down-regulated DE miRNAs responded to the hypo-salinity stress. • 1 up-regulated and 5 down-regulated DE miRNAs responded to the hyper-salinity stress. • 931 and 768 potential target genes were respectively predicted from the hypo- and hyper-salinity stress. • Genes involved in the vesicle-mediated transport and metal ion binding responded to the hypo-salinity stress. • Genes related to the signal transduction and metabolic process responded to the hyper-salinity stress. The Hong Kong oyster, Crassostrea hongkongensis , is a significant bivalve species with economic importance. It primarily inhabits the estuarine intertidal zones in southern China, making it susceptible to salinity fluctuations. Consequently, investigating the molecular mechanisms governing salinity regulation in C. hongkongensis is essential. In this study, we conducted miRNA-seq on C. hongkongensis to compare miRNA expression differences under varying salinities (5‰, 25‰, and 35‰). The miRNA sequencing revealed 51 known miRNAs and 95 novel miRNAs across nine small RNA libraries (S5, S25, and S35). Among these miRNAs, we identified 6 down-regulated differentially expressed (DE) miRNAs in response to hypo-salinity stress (5‰), while 1 up-regulated DE miRNA and 5 down-regulated DE miRNAs were associated with hyper-salinity stress (35‰). Additionally, we predicted 931 and 768 potential target genes for hypo- and hyper-salinity stress, respectively. Functional gene annotation indicated that the target genes under hypo-salinity stress were linked to vesicle-mediated transport and metal ion binding. Conversely, those under hyper-salinity stress were primarily involved in signal transduction and metabolic processes. These findings have provided insights into the regulatory role of miRNAs, their potential target genes and associated pathways in oyster hypo- and hyper-salinity stress, which establish a foundation for future studies on the roles of miRNAs in salinity acclimation mechanisms in C. hongkongensis. [ABSTRACT FROM AUTHOR]

ABSTRACT: The regulation of skeletal muscle growth and development in chicken is complex. MicroRNAs (miRNAs) have been found to play an important role in the process, and more research is needed to further understand the regulatory mechanism of miRNAs. In this study, leg muscles of Jinghai yellow chickens at 300 d with low body weight (slow-growing group) and high body weight (fast-growing group) were collected for miRNA sequencing (miRNA-seq) and Bioinformatics analysis revealed 12 differentially expressed miRNAs (DEMs) between the two groups. We predicted 150 target genes for the DEMs, and GO and KEGG pathway analysis showed the target genes of miR-24-3p and novel_miR_133 were most enriched in the terms related to growth and development. Moreover, networks of DEMs and target genes showed that miR-24-3p and novel_miR_133 were the 2 core miRNAs. Hence, miR-24-3p was selected for further functional exploration in chicken primary myoblasts (CPMs) with molecular biology technologies including qPCR, cell counting kit-8 (CCK-8), 5-ethynyl-2’-deoxyuridine (EdU) and immunofluorescence. When proliferating CPMs were transfected with miR-24-3p mimic, the expression of cyclin dependent kinase inhibitor 1A (P21) was up-regulated and both CCK-8 and EdU assays showed that the proliferation of CPMs was inhibited. However, when the inhibitor was transfected into the proliferating CPMs, the opposite results were found. In differentiated CPMs, transfection with miR-24-3p mimic resulted in up regulation of MYOD, MYOG and MYHC after 48 h. Myotube areas also increased significantly compared to the mimic negative control (NC) group. When treated with inhibitor, differentiation CPMs produced the opposite effects. Overall, we revealed 2 miRNAs (novel_miR_133 and miR-24-3p) significantly related with growth and development and further proved that miR-24-3p could suppress the proliferation and promote differentiation of CPMs. The results would facilitate understanding the effects of miRNAs on the growth and development of chickens at the post-transcriptional level and could also have an important guiding role in yellow-feathered chicken breeding.

Myometrial contraction is essential for successful delivery. Recent studies have highlighted the vital roles of tissue-derived exosomes in disease diagnostic, prognostic, and therapeutic applications; however, the characteristics of uterine myometrium-derived exosomes are unclear. Here, we successfully isolated exosomes from myometrial tissues, human myometrial smooth muscle cells (HMSMCs), and human umbilical vein endothelial cells (HUVECs), then performed quantitative liquid chromatography-tandem mass spectrometry and miRNA sequencing to investigate the cargo of the exosomes. Fifty-two proteins and five miRNAs were differentially expressed (DE) in term non-labor and term labor myometrium-derived exosomes. Among them, seven proteins (SERPINE1, THBS1, MGAT1, VIM, FGB, FGG, and VWF) were differentially expressed both in the myometrial exosomes and tissues, three miRNAs (miR-363-3p, miR-203a-3p, and miR-205-5p) target 13 DE genes. The top three miRNA derived from HMSMCs (miR-125b-1-3p, miR-337-5p, and miR-503-5p) and HUVECs (miR-663a, miR-4463, and miR-3622a-5p) were identified. Two proteins, GJA1 and SLC39A14, exist in female blood exosomes and are highly expressed in HMSMCs exosomes, are also upregulated in the laboring myometrium, which verified increased in laboring blood samples, might be novel potential biomarkers for myometrial activation. The proteomic and miRNA profile of exosomes derived from laboring myometrium revealed some molecules in the exosomes that affect the intercellular communication and the function of the myometrium. [ABSTRACT FROM AUTHOR]

Following radical surgery, patients may suffer a relapse. It is important to identify such patients so that therapy tactics can be modified appropriately. Existing stratification schemes do not display the probability of recurrence with enough precision since locally advanced prostate cancer (PCa) is classified as high-risk but is not ranked in greater detail. Between 40 and 50% of PCa cases belong to the TMPRSS2-ERG subtype that is a sufficiently homogeneous group for high-precision prognostic marker search to be possible. This study includes two independent cohorts and is based on high throughput sequencing and qPCR data. As a result, we have been able to suggest a perspective-trained model involving a deep neural network based on both qPCR data for mRNA and miRNA and clinicopathological criteria that can be used for recurrence risk forecasts in patients with TMPRSS2-ERG-positive, locally advanced PCa (the model uses ALDH3A2 + ODF2 + QSOX2 + hsa-miR-503-5p + ISUP + pT, with an AUC = 0.944). In addition to the prognostic model's use of identified differentially expressed genes and miRNAs, miRNA–target pairs were found that correlate with the prognosis and can be presented as an interactome network. [ABSTRACT FROM AUTHOR]

Pre-analytical factors have a significant influence on circulating microRNA (miRNA) profiling. The aim of this study was a comprehensive assessment of the impact of the anticoagulant type in blood collection tubes on circulating plasma miRNA profiles using small RNA sequencing. Blood from ten healthy participants (five males and five females from 25 to 40 years old) was taken in collection tubes with four different anticoagulants: acid citrate dextrose (ACD-B), sodium citrate, citrate-theophylline-adenosine-dipyridamole (CTAD) and dipotassium-ethylenediaminetetraacetic acid (K2 EDTA). Platelet-free plasma samples were obtained by double centrifugation. EDTA plasma samples had elevated levels of hemolysis compared to samples obtained using other anticoagulants. Small RNA was extracted from plasma samples and small RNA sequencing was performed on the Illumina NextSeq 500 system. A total of 30 samples had been successfully sequenced starting from ~1 M reads mapped to miRNAs, allowing us to analyze their diversity and isoform content. The principal component analysis showed that the EDTA samples have distinct circulating plasma miRNA profiles compared to samples obtained using other anticoagulants. We selected 50 miRNA species that were differentially expressed between the sample groups based on the type of anticoagulant. We found that the EDTA samples had elevated levels of miRNAs which are abundant in red blood cells (RBC) and associated with hemolysis, while the levels of some platelet-specific miRNAs in these samples were lowered. The ratio between RBC-derived and platelet-derived miRNAs differed between the EDTA samples and other sample groups, which was validated by quantitative PCR. This study provides full plasma miRNA profiles of 10 healthy adults, compares them with previous studies and shows that the profile of circulating miRNAs in the EDTA plasma samples is altered primarily due to an increased level of hemolysis. [ABSTRACT FROM AUTHOR]

We aimed to identify miRNAs and pathways specifically deregulated in adolescent and young adult (AYA) T-ALL patients. Small RNA-seq showed no major differences between AYA and pediatric T-ALL, but it revealed downregulation of miR-143-3p in T-ALL patients. Prediction algorithms identified several known and putative oncogenes targeted by this miRNA, including KRAS, FGF1, and FGF9. Pathway analysis indicated signaling pathways related to cell growth and proliferation, including FGFR signaling and PI3K-AKT signaling, with the majority of genes overrepresented in these pathways being predicted targets of hsa-miR-143-3p. By luciferase reporter assays, we validated direct interactions of this miRNA with KRAS, FGF1 and FGF9. In cell proliferation assays, we showed reduction of cell growth upon miR-143-3p overexpression in two T-ALL cell lines. Our study is the first description of the miRNA transcriptome in AYA T-ALL patients and the first report on tumor suppressor potential of miR-143-3p in T-ALL. Downregulation of this miRNA in T-ALL patients might contribute to enhanced growth and viability of leukemic cells. We also discuss the potential role of miR-143-3p in FGFR signaling. Although this requires more extensive validation, it might be an interesting direction, since FGFR inhibition proved promising in preclinical studies in various cancers. [ABSTRACT FROM AUTHOR]

Vibrio parahaemolyticus, as one of the main pathogens of marine vibriosis, has brought huge losses to aquaculture. However, the interaction mechanism between V. parahaemolyticus and Epinephelus coioides remains unclear. Moreover, there is a lack of comprehensive multi-omics analysis of the immune response of grouper spleen to V. parahaemolyticus. Herein, E. coioides was artificially injected with V. parahaemolyticus, and it was found that the mortality was 16.7% in the early stage of infection, and accompanied by obvious histopathological lesions in the spleen. Furthermore, 1586 differentially expressed genes were screened by mRNA-seq. KEGG analysis showed that genes were significantly enriched in immune-related pathways, Acute-phase immune response, Apoptosis, Complement system and Cytokine-cytokine receptor interaction. As for miRNA-seq analysis, a total of 55 significantly different miRNAs were identified. Further functional annotation analysis indicated that the target genes of differentially expressed miRNAs were enriched in three important pathways (Phosphatidylinositol signaling system, Lysosome and Focal adhesions). Through mRNA-miRNA integrated analysis, 1427 significant miRNA-mRNA pairs were obtained and "p53 signaling pathway", "Intestinal immune network for IgA production" were considered as two crucial pathways. Finally, miR-144-y, miR-497-x, novel-m0459-5p, miR-7133-y, miR-378-y, novel-m0440-5p and novel-m0084-3p may be as key miRNAs to regulate immune signaling pathways via the miRNA-mRNA interaction network. The above results suggest that the mRNA-miRNA integrated analysis not only sheds new light on the molecular mechanisms underlying the interaction between host and V. parahaemolyticus but also provides valuable and new insights into resistance to vibrio infection. [ABSTRACT FROM AUTHOR]

Diazinon (DZN), a persistent organophosphate insecticide, has been associated with neurotoxic effects, particularly in the hippocampus. However, the specific molecular mechanisms of DZN-induced hippocampal toxicity remain unknown. In this study, we analyzed the mRNA and miRNA expression patterns of HT22 cells following exposure to DZN (125 μM), and the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted subsequently. The integration of miRNA sequencing (miRNA-seq) and mRNA sequencing (mRNA-seq) data identified 33 differentially expressed miRNAs (DEMIs, 15 up-regulated and 18 down-regulated) and 271 differentially expressed mRNAs (DEMs, 69 up-regulated and 202 down-regulated) targeted by the DEMIs. Moreover, the 3 most central mRNAs (ITGAV, FN1, and EGFR) and 7 associated miRNAs (mmu-miR-700-5p, mmu-miR-26a-2-3p, mmu-miR-452-3p, mmu-miR-25-3p, mmu-miR-582-5p, mmu-miR-467a-5p, and mmu-miR-467b-5p) were screened and validated using quantitative real-time PCR. Furthermore, the GO analysis revealed that the identified DEMs were enriched in biological adhesion extracellular matrix, and growth factor binding, while the KEGG analysis suggested that the enriched DEMs were involved in ECM-receptor interaction, mTOR signaling pathway, MAPK signaling pathway, and AMPK signaling pathway. Our results may aid in elucidating the underlying mechanisms associated with DZN-induced hippocampal toxicity and provide valuable insights into the pathogenesis of neurotoxicity triggered by other organophosphorus pesticides., Competing Interests: Declaration of competing interest The authors declare no conflicts of interest!, (Copyright © 2024 Elsevier Ltd. All rights reserved.)

MicroRNA transcriptomes from fresh tumors and the adjacent normal tissues were profiled in 10 Korean patients diagnosed with lung adenocarcinoma using a next-generation sequencing (NGS) technique called miRNA-seq. The sequencing quality was assessed using FastQC, and low-quality or adapter-contaminated portions of the reads were removed using Trim Galore. Quality-assured reads were analyzed using miRDeep2 and Bowtie. The abundance of known miRNAs was estimated using the reads per million (RPM) normalization method. Subsequently, using DESeq2 and Wx, we identified differentially expressed miRNAs and potential miRNA biomarkers for lung adenocarcinoma tissues compared to adjacent normal tissues, respectively. We defined reliable miRNA biomarkers for lung adenocarcinoma as those detected by both methods. The miRNA-seq data are available in the Gene Expression Omnibus (GEO) database under accession number GSE196633, and all processed data can be accessed via the Mendeley data website.

Purpose: Hypoxia plays an essential role in the progression of hepatocellular carcinoma (HCC), whereas hypoxia inducible factor-1 (HIF-1) is the key transcription factor allowing HCC to survive hypoxia. The aim of this study was to define the essential mRNAs and miRNAs regulated by HIF1A and dissect their functions, interactions, and tumor-infiltrating immune cells in HCC. Methods: A human HCC cell line HepG2 was used as a cell model of HCC. The CRISPR/Cas9 system was used to knock out HIF1A in HepG2 cells, and RNA sequencing was utilized to characterize differentially expressed mRNAs and miRNAs in the HIF1A -knockout HepG2 cells; the identified candidates were then analyzed by GO annotation and KEGG pathway enrichment to study their function and establish a PPI network. Quantitative (q) PCR was used to verify if there were significant differences in the expression of mRNAs, and the association of the selected mRNAs expression with immune cell infiltration levels was further analyzed using The Cancer Genome Atlas (TCGA) pan-cancer data. Results: Using RNA-sequencing, we discovered that there were 1535 mRNAs differentially expressed (adjusted p < 0.05, |fold change|>1.5) in the HIF1A -knockout HepG2 cells, among which there were 644 mRNAs upregulated and 891 mRNAs downregulated. GO annotation and KEGG pathway enrichment showed that these mRNAs were involved in glycolysis/gluconeogenesis, PI3K-Akt signaling pathways, and HIF-1 signaling pathways. In addition, we found that there were 309 miRNAs differentially expressed (adjusted p < 0.05, |fold change|>1.5) in the HIF1A-knockout HepG2 cells, of which there were 213 miRNAs upregulated and 96 miRNAs downregulated. Our further analyses uncovered that these miRNA putative targets were involved in the hippo signaling pathway, axon guidance, and tight junction. Moreover, the construction and analysis of the PPI network showed that OASL , IL6 , and TAF1 were recognized as hub genes with the highest connectivity degrees. Importantly, in the HIF1A -knockout HepG2 cells, our qRT-PCR data confirmed the selected mRNA changes revealed by RNA-sequencing, and with TCGA pan-cancer data, we revealed that the expressional levels of these three genes, LUM , SCOC , and CCL2 , were associated with immune cell infiltration levels. Conclusion: The identified potential key network of mRNAs and miRNAs regulated by HIF1A in the HCC cells suggests a key role of HIF1A in the tumorigenesis of HCC. [ABSTRACT FROM AUTHOR]

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder, with mutations in hundreds of genes contributing to its risk. Herein, we studied lymphoblastoid cell lines (LCLs) from children diagnosed with autistic disorder (n = 10) and controls (n = 7) using RNA and miRNA sequencing profiles. The sequencing analysis identified 1700 genes and 102 miRNAs differentially expressed between the ASD and control LCLs (p ≤ 0.05). The top upregulated genes were GABRA4, AUTS2, and IL27, and the top upregulated miRNAs were hsa-miR-6813-3p, hsa-miR-221-5p, and hsa-miR-21-5p. The RT-qPCR analysis confirmed the sequencing results for randomly selected candidates: AUTS2, FMR1, PTEN, hsa-miR-15a-5p, hsa-miR-92a-3p, and hsa-miR-125b-5p. The functional enrichment analysis showed pathways involved in ASD control proliferation of neuronal cells, cell death of immune cells, epilepsy or neurodevelopmental disorders, WNT and PTEN signaling, apoptosis, and cancer. The integration of mRNA and miRNA sequencing profiles by miRWalk2.0 identified correlated changes in miRNAs and their targets' expression. The integration analysis found significantly dysregulated miRNA–gene pairs in ASD. Overall, these findings suggest that mRNA and miRNA expression profiles in ASD are greatly altered in LCLs and reveal numerous miRNA–gene interactions that regulate critical pathways involved in the proliferation of neuronal cells, cell death of immune cells, and neuronal development. [ABSTRACT FROM AUTHOR]

Infection with Salmonella Typhimurium (S. Typhimurium) is a common cause of food-borne zoonosis leading to acute gastroenteritis in humans and pigs, causing economic losses to producers and farmers, and generating a food security risk. In a previous study, we demonstrated that S. Typhimurium infection produces a severe transcriptional activation of inflammatory processes in ileum. However, little is known regarding how microRNAs regulate this response during infection. Here, small RNA sequencing was used to identify 28 miRNAs differentially expressed (DE) in ileum of S. Typhimurium-infected pigs, which potentially regulate 14 target genes involved in immune system processes such as regulation of cytokine production, monocyte chemotaxis, or cellular response to interferon gamma. Using in vitro functional and gain/loss of function (mimics/CRISPR-Cas system) approaches, we show that porcine miR-194a-5p (homologous to human miR-194-5p) regulates TLR4 gene expression, an important molecule involved in pathogen virulence, recognition and activation of innate immunity in Salmonella infection. [ABSTRACT FROM AUTHOR]

Throughout the world, numerous individuals are infected with Toxoplasma gondii, which may improve immunity against cancer. Furthermore, microRNAs (miRs) may be differentially expressed in the host upon infection with T. gondii. In the present study, RNA-sequencing analysis and reverse transcription-quantitative PCR revealed that miR-429-3p, miR-145a-5p, miR-211-5p, miR-31-3p and miR-135a-5p were determined to be downregulated, while miR-21a-3p, miR-135b-5p, miR-210-5p and miR-146-3p were upregulated in mice post-infection with T. gondii. Antitumor genes [TNF receptor superfamily member 11b, large tumor suppressor kinase (Lats)2 and Lats1] were identified as targets of miR-429-3p, miR-145a-5p, miR-211-5p, miR-31-3p and miR-135a-5p with a luciferase reporter assay. In addition, the protein levels of Lats2 and Lats1 were detected to be higher in T. gondii-infected mice than in the control group. Therefore, these results provide favorable evidence for the suppression of cancer upon T. gondii infection and may give novel ideas for the treatment of tumors. [ABSTRACT FROM AUTHOR]

Introduction: In previous work, we have shown the synergistic effects of combinations of vitamins A, D2, and D3 in reducing the proliferation of HCT-116 colon cancer cells. This combination also induced apoptosis and altered gene expression patterns as determined by transcriptomic profiling. Aims: The aims of this study were to investigate the effects of vitamins A and D in HCT-116 Crl:NU(NCr)-Foxn1nu mouse xenographs. Also, to determine potential mechanisms of action using miRNA-seq and correlated these data with results from RNA-seq. Methods: HCT-116 colon cancer cells were cultured and used for the xenograph study. Crl:NU(NCr)-Male Foxn1nu mice were injected subcutaneously with HCT-116-cells (1 X 106 in 100 µl 50% Matrigel and Sterile PBS), and after 7 days, were treated with vitamin A and D in feed for 21 days. The mice were imaged on day 28, sacrificed and the tumors were excised and measured. RNA was isolated from the HCT-116 cells and tumors, and RNA-seq and miRNA-seq were performed. esults: There was a concentration-dependent reduction in the HCT-116 colon cancer cell viability after treatment with vitamin A and D combinations. In Crl:NU(NCr)-Foxn1nu mice injected with HCT-116 colon tumor cells, treatment with vitamin A (25,000 IU) and vitamin D (4,000 IU) or vitamin A (35,000 IU) and vitamin D (5,000 IU) for 21 days significantly reduced tumor growth by ~38% and ~58%, respectively (p < 0.001). In the HCT-116 cells and excised tumors, treatment upregulated the expression of Bim, Bax, p53, and PTEN, and downregulated AXIN2, ID2 and DDX20 mRNAs, all well-known pro-apoptotic proteins, tumor suppressors, and molecules involved in the epithelial-mesenchymal transition. Ingenuity Pathway Analysis of miRNA-seq paired with RNA-seq showed correlations between miRNAs expression and expression of these genes. For example, miRNA-mRNA correlation pairing showed that treatment downregulated miR-30c-3p and miR125b-3p which was associated with upregulated TP53. Upregulation of the tumor suppressor PTEN was associated with the downregulation of eleven miRNAs, including oncogenic miR17 and miR21. Conclusion: The results suggest that the combinations of vitamins A and D reduced the HCT-116 tumor burden in mice and altered the expression of miRNAs directly associated with genes in the apoptosis, tumor suppression, and epithelialmesenchymal transition pathways. The data supports the hypothesis that vitamin A and D combinations impact multiple cancer signaling pathways and thus may be more effective in the treatment and prevention of cancer, as well as reducing cancer metastasis. [ABSTRACT FROM AUTHOR]

Purpose: Hypoxia plays an essential role in the progression of hepatocellular carcinoma (HCC), whereas hypoxia inducible factor-1 (HIF-1) is the key transcription factor allowing HCC to survive hypoxia. The aim of this study was to define the essential mRNAs and miRNAs regulated by HIF1A and dissect their functions, interactions, and tumor-infiltrating immune cells in HCC.Methods: A human HCC cell line HepG2 was used as a cell model of HCC. The CRISPR/Cas9 system was used to knock out HIF1A in HepG2 cells, and RNA sequencing was utilized to characterize differentially expressed mRNAs and miRNAs in the HIF1A-knockout HepG2 cells; the identified candidates were then analyzed by GO annotation and KEGG pathway enrichment to study their function and establish a PPI network. Quantitative (q) PCR was used to verify if there were significant differences in the expression of mRNAs, and the association of the selected mRNAs expression with immune cell infiltration levels was further analyzed using The Cancer Genome Atlas (TCGA) pan-cancer data.Results: Using RNA-sequencing, we discovered that there were 1535 mRNAs differentially expressed (adjusted p < 0.05, |fold change|>1.5) in the HIF1A-knockout HepG2 cells, among which there were 644 mRNAs upregulated and 891 mRNAs downregulated. GO annotation and KEGG pathway enrichment showed that these mRNAs were involved in glycolysis/gluconeogenesis, PI3K-Akt signaling pathways, and HIF-1 signaling pathways. In addition, we found that there were 309 miRNAs differentially expressed (adjusted p < 0.05, |fold change|>1.5) in the HIF1A-knockout HepG2 cells, of which there were 213 miRNAs upregulated and 96 miRNAs downregulated. Our further analyses uncovered that these miRNA putative targets were involved in the hippo signaling pathway, axon guidance, and tight junction. Moreover, the construction and analysis of the PPI network showed that OASL, IL6, and TAF1 were recognized as hub genes with the highest connectivity degrees. Importantly, in the HIF1A-knockout HepG2 cells, our qRT-PCR data confirmed the selected mRNA changes revealed by RNA-sequencing, and with TCGA pan-cancer data, we revealed that the expressional levels of these three genes, LUM, SCOC, and CCL2, were associated with immune cell infiltration levels.Conclusion: The identified potential key network of mRNAs and miRNAs regulated by HIF1A in the HCC cells suggests a key role of HIF1A in the tumorigenesis of HCC.

• Inhibition of JMJD3 and UTX methyltransferase activity with the GSK-J4 treatment modulate miR-27a-3p and miR-423-5p expression levels. • miR-27a-3p and miR-423-5p are involved in proliferation of differentiated liver cells induced by GSK-J4. • miR-27a-3p and miR-423-5p do not influence retro -differentiation of differentiated liver cells induced by GSK-J4. • miR-27a-3p and miR-423-5p contribute to the modulation of hepatic cell plasticity. Epigenetic mechanisms, including DNA methylation, histone modifications, and chromatin remodeling, are highly involved in the regulation of hepatocyte viability, proliferation, and plasticity. We have previously demonstrated that repression of H3K27 methylation in differentiated hepatic HepaRG cells by treatment with GSK-J4, an inhibitor of JMJD3 and UTX H3K27 demethylase activity, changed their phenotype, inducing differentiated hepatocytes to proliferate. In addition to the epigenetic enzymatic role in the regulation of the retro -differentiation process, emerging evidence indicate that microRNAs (miRNAs) are involved in controlling hepatocyte proliferation during liver regeneration. Hence, the aim of this work is to investigate the impact of H3K27 methylation on miRNAs expression profile and its role in the regulation of the differentiation status of human hepatic progenitors HepaRG cells. A miRNA-sequencing was carried out in differentiated HepaRG cells treated or not with GSK-J4. Target searching and Gene Ontology analysis were performed to identify the molecular processes modulated by differentially expressed miRNAs. The biological functions of selected miRNAs was further investigated by transfection of miRNAs inhibitors or mimics in differentiated HepaRG cells followed by qPCR analysis, albumin ELISA assay, CD49a FACS analysis and EdU staining. We identified 12 miRNAs modulated by GSK-J4; among these, miR-27a-3p and miR- 423-5p influenced the expression of several proliferation genes in differentiated HepaRG cells. MiR-27a-3p overexpression increased the number of hepatic cells reentering proliferation. Interestingly, both miR-27a-3p and miR-423-5p did not affect the expression levels of genes involved in the differentiation of progenitors HepaRG cells. Modulation of H3K27me3 methylation in differentiated HepaRG cells, by GSK-J4 treatment, influenced miRNA' s expression profile pushing liver cells towards a proliferating phenotype. We demonstrated the involvement of miR-27a-3p in reinducing proliferation of differentiated hepatocytes suggesting a potential role in liver plasticity. [ABSTRACT FROM AUTHOR]

High-throughput sequencing of small RNA molecules such as microRNAs (miRNAs) has become a widely used approach for studying gene expression and regulation. However, analyzing miRNA-Seq data can be challenging because it requires multiple steps, from quality control and preprocessing to differential expression and pathway-enrichment analyses, with many tools and databases available for each step. Furthermore, reproducibility of the analysis pipeline is crucial to ensure that the results are accurate and reliable. Here, we present myBrain-Seq, a comprehensive and reproducible pipeline for analyzing miRNA-Seq data that incorporates miRNA-specific solutions at each step of the analysis. The pipeline was designed to be flexible and user-friendly, allowing researchers with different levels of expertise to perform the analysis in a standardized and reproducible manner, using the most common and widely used tools for each step. In this work, we describe the implementation of myBrain-Seq and demonstrate its capacity to consistently and reproducibly identify differentially expressed miRNAs and enriched pathways by applying it to a real case study in which we compared schizophrenia patients who responded to medication with treatment-resistant schizophrenia patients to obtain a 16-miRNA treatment-resistant schizophrenia profile.

17α-Methyltestosterone (17MT), a synthetic organic compound commonly found in sewage waters, can affect reproduction in aquatic animals, such as tilapia and yellow catfish. In the present study, male Gobiocypris rarus were exposed to 25, 50, and 100 ng/L of 17α-methyltestosterone (17MT) for 7 days. We first analyzed miRNA- and RNA-seq results to determine miRNA-target gene pairs and then developed miRNA-mRNA interactive networks after 17MT administration. Total weights, total lengths, and body lengths were not significantly different between the test groups and control groups. The paraffin slice method was applied to testes of G. rarus in the MT exposure and control groups. We found that there were more mature sperm (S) and fewer secondary spermatocytes (SSs) and spermatogonia (SGs) in the testes of control groups. As 17MT concentration increased, fewer and fewer mature sperm (S) were observed in the testes of male G. rarus. The results showed that FSH, 11-KT, and E2 were significantly higher in individuals exposed to 25 ng/L 17MT compared with the control groups. VTG, FSH, LH, 11-KT, and E2 were significantly lower in the 50 ng/L 17MT exposure groups compared to the control groups. VTG, FSH, LH, 11-KT, E2, and T were significantly lower in the groups exposed to 100 ng/L 17MT. High-throughput sequencing revealed 73,449 unigenes, 1205 known mature miRNAs, and 939 novel miRNAs in the gonads of G. rarus. With miRNA-seq, 49 (MT25-M vs. Con-M), 66 (MT50-M vs. Con-M), and 49 (MT100-M vs. Con-M) DEMs were identified in the treatment groups. Five mature miRNAs (miR-122-x, miR-574-x, miR-430-y, lin-4-x, and miR-7-y), as well as seven differentially expressed genes (soat2, inhbb, ihhb, gatm, faxdc2, ebp, and cyp1a1), which may be associated with testicular development, metabolism, apoptosis, and disease response, were assayed using qRT-PCR. Furthermore, miR-122-x (related to lipid metabolism), miR-430-y (embryonic development), lin-4-x (apoptosis), and miR-7-y (disease) were differentially expressed in the testes of 17MT-exposed G. rarus. This study highlights the role of miRNA-mRNA pairs in the regulation of testicular development and immune response to disease and will facilitate future studies on the miRNA-RNA-associated regulation of teleost reproduction.

By analyzing the sequencing results, we found that PKM2 was significantly enriched in the KEGG pathway associated with muscle development and was screened as a hub gene in the fast-growing group; hence, we speculate that PKM2 may be associated with the growth and development of Bian chickens. Keywords: miRNA-seq; skeletal muscle; Bian chicken; growth and development EN miRNA-seq skeletal muscle Bian chicken growth and development N.PAG N.PAG 15 04/26/22 20220415 NES 220415 1. Apart from research on genes related to Bian chicken growth, there are few studies on miRNA functions in Bian chicken growth. miRNA is a noncoding sequence with a length of approximately 22 nt, which inhibits the expression of target genes by binding to the 3'UTR region of mRNA [[5]]. In order to further explore the key DEMs and predicted target genes, we performed an interaction network analysis of genes enriched in KEGG pathways associated with muscle growth and development. [Extracted from the article]