The Shitou goose, the largest meat-type goose breed, is an ideal model for offering insights into enhancing meat production efficiency through understanding its genetic regulation of muscle development. Here, through whole-transcriptomic analysis of embryonic leg muscles, we identified 847 differentially expressed genes (DEG), 244 differentially expressed lncRNAs (DEL), 37 differentially expressed circRNAs (DEC), and 84 differentially expressed miRNAs (DEM). Gene ontology (GO) analysis highlighted the significant enrichment of differentially expressed RNAs in muscle structure development, actin filament-based processes, and the actin cytoskeleton pathway. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified pathways associated with the FoxO signaling pathway, AMPK signaling pathway, Wnt signaling pathway and calcium signaling pathway. Furthermore, we utilized Miranda, TargetScan, and miRDB to identify regulatory networks that involve interactions between lncRNA-mRNA, circRNA-mRNA, miRNA-mRNA, lncRNA-miRNA-mRNA, and circRNA-miRNA-mRNA, which regulated the growth and development of skeletal muscle. Notably, differentially expressed genes within the ceRNA network were most significantly enriched in the regulation of actin cytoskeletal organization. Additionally, a lncRNA/circRNA-miRNA-mRNA ceRNA network related to muscle growth and development was constructed based on protein-protein interaction (PPI) analysis and hub genes selection using Cytoscape. This further elucidated the regulatory roles of noncoding RNAs (ncRNA) in the formation of muscle fibers in Shitou goose. In summary, this study provides a valuable transcriptional regulatory network for goose muscle development laying the groundwork for further exploration of the molecular regulatory mechanisms underlying the excellent meat production performance of Shitou goose.

Coloration is a crucial trait that allows species to adapt and survive in different environments. Wild boars exhibit alternating black (dark) and yellow (light) longitudinal stripes on their back during their infancy (juvenile stripes), and as adults, they transform into uniform wild-type coat color. Aiming to record the procedure of juvenile stripes disappearing, piglets (WD) with juvenile stripes were produced by crossing a wild boar with Duroc sows, and photos of their coat color were taken from 20 d to 220 d. The pigments in the hairs from the black and yellow stripes were determined. Furthermore, the differentially expressed genes between the black and yellow stripes were investigated in 5 WD with the age of 30 d using whole-transcriptome sequencing to explore the genetic mechanism of the juvenile stripes. The juvenile stripes started to disappear at about 70 d, and stripes were not distinguished with the naked eye at about 160 d; that is, the juvenile stripe completely disappeared. A hotspot of a differentially expressing (DE) region was found on chromosome 13, containing/covering 2 of 13 DE genes and 8 of 10 DE lncRNAs in this region. A network among ZIC4, ssc-miR-532-3p, and ENSSSCG00000056225 might regulate the formation of juvenile stripes. Altogether, this study provides new insights into spatiotemporal coat color pattern.

Most viral infections can be self-limited, with no requirement for medical intervention. However, the same viruses can cause severe diseases in patients with compromised immunity due to single-gene diseases, acquired immune deficiency syndrome, or hematologic malignancies or those receiving immunosuppressive drugs. Occasionally, these immunocompromised patients harbor >1 infectious agent, requiring several concomitant diagnostic tests. We have developed, to our knowledge, a previously unreported whole-transcriptome sequencing-based pipeline that allows virome profiling, quantitation, and expression pattern analysis of 926 distinct viruses by sequencing of RNA isolated from a single lesional skin biopsy. This pipeline can also explore host genetics if there is a Mendelian predisposition to infection. We applied this pipeline to 6 Iranian patients with viral-induced skin lesions associated with immune deficiency secondary to HIV, human T-lymphotropic virus 1, chronic lymphocytic leukemia, and post transplant immunosuppression. In 5 cases, definitive human papillomavirus infections were identified, some caused by multiple viral types. In addition to human papillomavirus, coinfection with other viruses (Merkle cell polyomavirus, cytomegalovirus, and human herpesvirus 4) was detected in some lesions. In 1 case, whole-transcriptome sequencing validated the clinical diagnosis of adult T-cell leukemia/lymphoma in a patient with an initial diagnosis of mycosis fungoides/Sézary syndrome. These findings attest to the power of whole-transcriptome sequencing in profiling the cutaneous virome in the context of compromised immunity.

Choroid plexus tumors (CPTs) are intraventricular tumors derived from the choroid plexus epithelium and occur frequently in children. The aim of this study was to investigate the genomic and epigenomic characteristics of CPT and identify the differences between choroid plexus papilloma (CPP) and choroid plexus carcinoma (CPC). We conducted multiomics analyses of 20 CPT patients including CPP and CPC. Multiomics analysis included whole-genome sequencing, whole-transcriptome sequencing, and methylation sequencing. Mutually exclusive TP53 and EPHA7 point mutations, coupled with the amplification of chromosome 1, were exclusively identified in CPC. In contrast, amplification of chromosome 9 was specific to CPP. Differential gene expression analysis uncovered a significant overexpression of genes related to cell cycle regulation and epithelial-mesenchymal transition pathways in CPC compared to CPP. Overexpression of genes associated with tumor metastasis and progression was observed in the CPC subgroup with leptomeningeal dissemination. Furthermore, methylation profiling unveiled hypomethylation in major repeat regions, including long interspersed nuclear elements, short interspersed nuclear elements, long terminal repeats, and retrotransposons in CPC compared to CPP, implying that the loss of epigenetic silencing of transposable elements may play a role in tumorigenesis of CPC. Finally, the differential expression of AK1, regulated by both genomic and epigenomic factors, emerged as a potential contributing factor to the histological difference of CPP against CPC. Our results suggest pronounced genomic and epigenomic disparities between CPP and CPC, providing insights into the pathogenesis of CPT at the molecular level.

The economic losses incurred due to reduced muscle pigmentation highlight the crucial role of melanin-based coloration in the meat of black-bone chickens. Melanogenesis in the breast muscle of black-bone chickens is currently poorly understood in terms of molecular mechanisms. This study employed whole-transcriptome sequencing to analyze black and white breast muscle samples from black-bone chickens, leading to the identification of 367 differentially expressed (DE) mRNAs, 48 DElncRNAs, 104 DEcircRNAs, and 112 DEmiRNAs involved in melanin deposition. Based on these findings, a competitive endogenous RNA (ceRNA) network was developed to better understand the complex mechanisms of melanin deposition. Furthermore, our analysis revealed key DEmRNAs (TYR, DCT, EDNRB, MLPH and OCA2) regulated by DEmiRNAs (gga-miR-140-5p, gga-miR-1682, gga-miR-3529, gga-miR-499-3p, novel-m0012-3p, gga-miR-200b-5p, gga-miR-203a, gga-miR-6651-5p, gga-miR-7455-3p, gga-miR-31-5p, miR-140-x, miR-455-x, novel-m0065-3p, gga-miR-29b-1-5p, miR-455-y, novel-m0085-3p, and gga-miR-196-1-3p). These DEmiRNAs competitively interacted with DElncRNAs including MSTRG.2609.2, MSTRG.4185.1, LOC112530666, LOC112533366, LOC771030, LOC107054724, LOC121107411, LOC100859072, LOC101750037, LOC121108550, LOC121109224, LOC121110876, and LOC101749016, as well as DEcircRNAs, such as novel_circ_000158, novel_circ_000623, novel_001518, and novel_circ_003596. The findings from this study provide insight into the mechanisms that regulate lncRNA, circRNA, miRNA, and mRNA expression in chicken melanin deposition.

Gastric cancer remains a highly prevalent malignancy worldwide with its molecular features poorly understood. To gain full insight into its genetic landscape, we performed whole-transcriptome sequencing on human tumors and adjacent non-tumors to predict the function of microRNA, long coding RNA, circular RNA, and mRNA, as well as estimate their correlation with gastric cancer characteristics through construction of ceRNA, WGCNA and PPI network. Functional enrichment analysis annotated nucleic acid binding, enzyme activity and binding related to differentially expressed miRNAs (dif-miRNAs); energy binding and enzyme binding related to dif-lncRNAs; protein binding and enzyme activity related to dif-circRNAs; protein digestion and absorption related to dif-mRNAs. The expression of key miR-135a-5p, lncRNAs-MSTRG.48856.1, ENST00000569981, MSTRG.22826.1, ENST00000564492, circRNAs-CCSER2, FNDC3B, CORO1C, FAM214A were validated by real-time PCR. The ceRNA network filtered 14 miRNAs, 30 lncRNAs, and 6 mRNA in the lncRNA-ceRNA axis and 8 miRNAs, 9 circRNAs, and 3 mRNA in the circRNA-ceRNA axis. Genes involved in ceRNA were annotated to be closely related to tumor material synthesis and metabolism. The WGCNA network filtered gene clusters related to TNM traits and extracted the hub genes CLDN10, CD177, newGene_35523, newGene_51201, CEACAM7, and newGene_46634. These genes were associated with cell proliferation, metabolism, and enzyme activity regulation. The PPI network analyzed the stable interaction relationships of the hub genes. Our research provides a valuable resource for understanding the molecular mechanisms of gastric cancer from the perspective of tumor metabolism.

Angiosarcomas (ASs) are rare malignant vascular entities that can affect several regions in our body, including the heart. Cardiac ASs comprise 25-40% of cardiac sarcomas and can cause death within months of diagnosis. Thus, our aim was to identify potential differences and/or similarities between cardiac and extra-cardiac ASs to enhance targeted therapies and, consequently, patients\' prognosis. Whole-transcriptome analysis of three cardiac and eleven extra-cardiac non-cutaneous samples was performed to investigate differential gene expression and mutational events between the two groups. The gene signature of cardiac and extra-cardiac non-cutaneous ASs was also compared to that of cutaneous angiosarcomas (n = 9). H/N/K-RAS and TP53 alterations were more recurrent in extra-cardiac ASs, while POTE-gene family overexpression was peculiar to cardiac ASs. Additionally, in vitro functional analyses showed that POTEH upregulation conferred a growth advantage to recipient cells, partly supporting the cardiac AS aggressive phenotype and patients\' scarce survival rate. These features should be considered when investigating alternative treatments.

Duck meat is pivotal in providing high-quality protein for human nutrition, underscoring the importance of studying duck myogenesis. The regulatory mechanisms governing duck myogenesis involve both coding and non-coding RNAs, yet their specific expression patterns and molecular mechanisms remain elusive. To address this knowledge gap, we performed expression profiling analyses of mRNAs, lncRNAs, circRNAs, and miRNAs involved in duck myogenesis using whole-transcriptome RNA-seq. Our analysis identified 1733 differentially expressed (DE)-mRNAs, 1116 DE-lncRNAs, 54 DE-circRNAs, and 174 DE-miRNAs when comparing myoblasts and myotubes. A GO analysis highlighted the enrichment of DE molecules in the extracellular region, protein binding, and exocyst. A KEGG analysis pinpointed pathways related to ferroptosis, PPAR signaling, nitrogen metabolism, cell cycle, cardiac muscle contraction, glycerolipid metabolism, and actin cytoskeleton. A total of 51 trans-acting lncRNAs, including ENSAPLT00020002101 and ENSAPLT00020012069, were predicted to participate in regulating myoblast proliferation and differentiation. Based on the ceRNAs, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA ceRNA networks involving five miRNAs (miR-129-5p, miR-133a-5p, miR-22-3p, miR-27b-3p, and let-7b-5p) that are relevant to myogenesis. Furthermore, the GO and KEGG analyses of the DE-mRNAs within the ceRNA network underscored the significant enrichment of the glycerolipid metabolism pathway. We identified five different DE-mRNAs, specifically ENSAPLG00020001677, ENSAPLG00020002183, ENSAPLG00020005019, ENSAPLG00020010497, and ENSAPLG00020017682, as potential target genes that are crucial for myogenesis in the context of glycerolipid metabolism. These five mRNAs are integral to ceRNA networks, with miR-107_R-2 and miR-1260 emerging as key regulators. In summary, this study provides a valuable resource elucidating the intricate interplay of mRNA-lncRNA-circRNA-miRNA in duck myogenesis, shedding light on the molecular mechanisms that govern this critical biological process.

OBJECTIVE: To explore the key genes involved in the occurrence and development of glioblastoma (GBM) by analyzing whole-transcriptome sequencing and biologic data from GBM and normal cerebral cortex tissues and to search for important noncoding RNA (ncRNA) molecular markers based on the competitive endogenous RNA (ceRNA) network.
METHODS: Ten GBM and normal cerebral cortex tissues were collected for full transcriptome sequencing, screened for differentially expressed (DE) mRNAs, miRNAs, lncRNAs, and circRNAs, and subjected to bioinformatic analysis. We constructed a Protein-Protein Interaction (PPI) network and a circRNA/lncRNA-miRNA-mRNA regulatory network and identified them using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Finally, The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) databases were used to validate and conduct a survival analysis of the target genes.
RESULTS: A total of 5341 DEmRNAs, 259 DEmiRNAs, 3122 DElncRNAs, and 2135 DEcircRNAs were identified. Enrichment analysis showed that target genes regulated by DEmiRNA, DElncRNA, and DEcircRNA were closely related to chemical synaptic transmission and ion transmembrane transport. A PPI network analysis screened 10 hub genes that directly participate in tumor cell mitosis regulation. In addition, the ceRNA composite network showed that hsa-miR-296-5p and hsa-miR-874-5p were the central nodes of the network, and the reliability of relevant key molecules was successfully verified through RT-qPCR identification and the TCGA database. The CGGA database survival analysis produced 8 DEmRNAs closely related to GBM patient survival prognosis.
CONCLUSIONS: This study revealed the important regulatory functions and molecular mechanisms of ncRNA molecules and identified hsa-miR-296-5p and hsa-miR-874-5p as key molecules in the ceRNA network. They may play an important role in GBM pathogenesis, treatment, and prognosis.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is one of the world\'s leading causes of death and a major chronic respiratory disease. Aerobic exercise, the cornerstone of pulmonary rehabilitation, improves prognosis of COPD patients; however, few studies have comprehensively examined the changes in RNA transcript levels and the crosstalk between various transcripts in this context. This study identified the expression of RNA transcripts in COPD patients who engaged in aerobic exercise training for 12 weeks, and further constructions of the possible RNAs networks were made.
METHODS: Peripheral blood samples for all four COPD patients who benefited from 12 weeks of PR were collected pre- and post-aerobic exercises and evaluated for the expression of mRNA, miRNA, lncRNA, and circRNA with high-throughput RNA sequencing followed by GEO date validation. In addition, enrichment analyses were conducted on different expressed mRNAs. LncRNA-mRNA and circRNA-mRNA coexpression networks, as well as lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA competing expression networks (ceRNAs) in COPD were constructed.
RESULTS: We identified and analyzed the differentially expressed mRNAs and noncoding RNAs in the peripheral blood of COPD patients\' post-exercise. Eighty-six mRNAs, 570 lncRNAs, 8 miRNAs, and 2087 circRNAs were differentially expressed. Direct function enrichment analysis and Gene Set Variation Analysis showed that differentially expressed RNAs(DE-RNAs) correlated with several critical biological processes such as chemotaxis, DNA replication, anti-infection humoral response, oxidative phosphorylation, and immunometabolism, which might affect the progression of COPD. Some DE-RNAs were validated by Geo databases and RT-PCR, and the results were highly correlated with RNA sequencing. We constructed ceRNA networks of DE-RNAs in COPD.
CONCLUSIONS: The systematic understanding of the impact of aerobic exercise on COPD was achieved using transcriptomic profiling. This research offers a number of potential candidates for clarifying the regulatory mechanisms that exercise has on COPD, which could ultimately help in understanding the pathophysiology of COPD.