Epigenomic and Transcriptomic Changes During Human RPE EMT in a Stem Cell Model of Epiretinal Membrane Pathogenesis and Prevention by Nicotinamide

Epithelial to mesenchymal transition (EMT) is a biological process involved in tissue morphogenesis and disease that causes dramatic changes in cell morphology, migration, proliferation, and gene expression. The retinal pigment epithelium (RPE), which supports the neural retina, can undergo EMT, producing fibrous epiretinal membranes (ERMs) associated with vision-impairing clinical conditions, such as macular pucker and proliferative vitreoretinopathy (PVR). We found that co-treatment with TGF-β and TNF-α (TNT) accelerates EMT in adult human RPE stem cell-derived RPE cell cultures. We captured the global epigenomic and transcriptional changes elicited by TNT treatment of RPE and identified putative active enhancers associated with actively transcribed genes, including a set of upregulated transcription factors that are candidate regulators. We found that the vitamin B derivative nicotinamide downregulates these key transcriptional changes, and inhibits and partially reverses RPE EMT, revealing potential therapeutic routes to benefit patients with ERM, macular pucker and PVR.

The Planktonic Core Microbiome and Core Functions in the Cattle Rumen by Next Generation Sequencing

The cow rumen harbors a great variety of diverse microbes, which form a complex, organized community. Understanding the behavior of this multifarious network is crucial in improving ruminant nutrient use efficiency. The aim of this study was to expand our knowledge by examining 10 Holstein dairy cow rumen fluid fraction whole metagenome and transcriptome datasets. DNA and mRNA sequence data, generated by Ion Torrent, was subjected to quality control and filtering before analysis for core elements. The taxonomic core microbiome consisted of 48 genera belonging to Bacteria (47) and Archaea (1). The genus Prevotella predominated the planktonic core community. Core functional groups were identified using co-occurrence analysis and resulted in 587 genes, from which 62 could be assigned to metabolic functions. Although this was a minimal functional core, it revealed key enzymes participating in various metabolic processes. A diverse and rich collection of enzymes involved in carbohydrate metabolism and other functions were identified. Transcripts coding for enzymes active in methanogenesis made up 1% of the core functions. The genera associated with the core enzyme functions were also identified. Linking genera to functions showed that the main metabolic pathways are primarily provided by Bacteria and several genera may serve as a “back-up” team for the central functions. The key actors in most essential metabolic routes belong to the genus Prevotella. Confirming earlier studies, the genus Methanobrevibacter carries out the overwhelming majority of rumen methanogenesis and therefore methane emission mitigation seems conceivable via targeting the hydrogenotrophic methanogenesis.

The Whole Transcriptome Involved in Denervated Muscle Atrophy Following Peripheral Nerve Injury

Peripheral nerve injury (PNI) usually leads to progressive muscle atrophy and poor functional recovery. Previous studies have demonstrated that non-coding ribonucleic acid (ncRNA) is a key regulator of muscle atrophy and beneficial for the treatment of PNI. We aimed to analyze the whole transcriptome involved in denervated muscle atrophy after PNI. Animal models of sciatic nerve injury were assessed at 0 (control group), 1, 2, 4, and 8 weeks after injury. The expression patterns in the whole transcriptome in the gastrocnemius muscle were profiled using RNA sequencing at each time point and compared to that obtained in the control group. Six-hundred and sixty-four long non-coding RNAs, 671 microRNAs, 236 circular RNAs, and 12,768 messenger RNAs (mRNAs) were differentially expressed (DE) after injury. Changes in some of the DE ncRNAs and mRNAs were validated using quantitative polymerase chain reaction. Gene Ontology and Kyoko Encyclopedia of Genes and Genomes analysis revealed the potential functions of and relationships among the DE ncRNAs and mRNAs. To our knowledge, this is the first study to expound the whole transcriptome involved in denervated muscle atrophy, and provides a theoretical basis for further research targeting ncRNAs.

Integrative analysis of gut microbiome and host transcriptomes reveals associations between treatment outcomes and immunotherapy-induced colitis

Immune checkpoint inhibitors (ICIs) are widely used to treat various malignancies. Although the gut microbiome is known to influence the efficacy of ICIs on epithelial tumors, the functional interactions between gut taxa and colonic mucosa remain poorly understood. Here we performed transcriptomic profiling and 16S rRNA sequencing to investigate the relationships between mucosal gene expression and microbial composition with ICI responses and gastrointestinal immune-related adverse events (GI irAEs). In responders, genes related to DNA repair and cell cycle signatures were enriched in responders whereas signatures related to innate immune response, NFAT and IFN-γ signaling pathways were enriched in nonresponders. Gut microbial composition revealed an association between moderate GI irAE and favorable response to ICI therapy. Favorable therapeutic responses to ICI and GI irAE treatments were associated with taxa classified as Enterobacteriaceae and were related to ribonucleoprotein complex biogenesis, cytokine-mediated signaling pathway, tRNA metabolic process, and ribonucleoprotein complex assembly in the colon. These findings open new perspectives for improving the efficacy and safety of cancer immunotherapy.

Divergence in Morris Water Maze-Based Cognitive Performance under Chronic Stress Is Associated with the Hippocampal Whole Transcriptomic Modification in Mice

Individual susceptibility determines the magnitude of stress effects on cognitive function. The hippocampus, a brain region of memory consolidation, is vulnerable to stressful environments, and the impact of stress on hippocampus may determine individual variability in cognitive performance. Therefore, the purpose of this study was to define the relationship between the divergence in spatial memory performance under chronically unpredictable stress and an associated transcriptomic alternation in hippocampus, the brain region of spatial memory consolidation. Multiple strains of BXD (B6 × D2) recombinant inbred mice went through a 4-week chronic variable stress (CVS) paradigm, and the Morris water maze (MWM) test was conducted during the last week of CVS to assess hippocampal-dependent spatial memory performance and grouped animals into low and high performing groups based on the cognitive performance. Using hippocampal whole transcriptome RNA-sequencing data, differential expression, PANTHER analysis, WGCNA, Ingenuity's upstream regulator analysis in the Ingenuity Pathway Analysis® and phenotype association analysis were conducted. Our data identified multiple genes and pathways that were significantly associated with chronic stress-associated cognitive modification and the divergence in hippocampal dependent memory performance under chronic stress. Biological pathways associated with memory performance following chronic stress included metabolism, neurotransmitter and receptor regulation, immune response and cellular process. The Ingenuity's upstream regulator analysis identified 247 upstream transcriptional regulators from 16 different molecule types. Transcripts predictive of cognitive performance under high stress included genes that are associated with a high occurrence of Alzheimer's and cognitive impairments (e.g., Ncl, Eno1, Scn9a, Slc19a3, Ncstn, Fos, Eif4h, Copa, etc.). Our results show that the variable effects of chronic stress on the hippocampal transcriptome are related to the ability to complete the MWM task and that the modulations of specific pathways are indicative of hippocampal dependent memory performance. Thus, the divergence in spatial memory performance following chronic stress is related to the unique pattern of gene expression within the hippocampus.

Deciphering the RNA landscape by RNAome sequencing

Current RNA expression profiling methods rely on enrichment steps for specific RNA classes, thereby not detecting all RNA species in an unperturbed manner. We report strand-specific RNAome sequencing that determines expression of small and large RNAs from rRNA-depleted total RNA in a single sequence run. Since current analysis pipelines cannot reliably analyze small and large RNAs simultaneously, we developed TRAP, Total Rna Analysis Pipeline, a robust interface that is also compatible with existing RNA sequencing protocols. RNAome sequencing quantitatively preserved all RNA classes, allowing cross-class comparisons that facilitates the identification of relationships between different RNA classes. We demonstrate the strength of RNAome sequencing in mouse embryonic stem cells treated with cisplatin. MicroRNA and mRNA expression in RNAome sequencing significantly correlated between replicates and was in concordance with both existing RNA sequencing methods and gene expression arrays generated from the same samples. Moreover, RNAome sequencing also detected additional RNA classes such as enhancer RNAs, anti-sense RNAs, novel RNA species and numerous differentially expressed RNAs undetectable by other methods. At the level of complete RNA classes, RNAome sequencing also identified a specific global repression of the microRNA and microRNA isoform classes after cisplatin treatment whereas all other classes such as mRNAs were unchanged. These characteristics of RNAome sequencing will significantly improve expression analysis as well as studies on RNA biology not covered by existing methods.

Comparative transcriptome profiling approach to glean virulence and immunomodulation-related genes of Fasciola hepatica

BACKGROUND

Fasciola hepatica causes chronic liver disease, fasciolosis, leading to significant losses in the livestock economy and concerns for human health in many countries. The identification of F. hepatica genes involved in the parasite's virulence through modulation of host immune system is utmost important to comprehend evasion mechanisms of the parasite and develop more effective strategies against fasciolosis. In this study, to identify the parasite's putative virulence genes which are associated with host immunomodulation, we explored whole transcriptome of an adult F. hepatica using current transcriptome profiling approaches integrated with detailed in silico analyses. In brief, the comparison of the parasite transcripts with the specialised public databases containing sequence data of non-parasitic organisms (Dugesiidae species and Caenorhabditis elegans) or of numerous pathogens and investigation of the sequences in terms of nucleotide evolution (directional selection) and cytokine signaling relation were conducted.

RESULTS

NGS of the whole transcriptome resulted in 19,534,766 sequence reads, yielding a total of 40,260 transcripts (N₅₀ = 522 bp). A number of the parasite transcripts (n = 1,671) were predicted to be virulence-related on the basis of the exclusive homology with the pathogen-associated data, positive selection or relationship with cytokine signaling. Of these, a group of the virulence-related genes (n = 62), not previously described, were found likely to be associated with immunomodulation based on in silico functional categorisation, showing significant sequence similarities with various immune receptors (i.e. MHC I class, TGF-β receptor, toll/interleukin-1 receptor, T-cell receptor, TNF receptor, and IL-18 receptor accessory protein), cytokines (i.e. TGF-β, interleukin-4/interleukin-13 and TNF-α), cluster of differentiations (e.g. CD48 and CD147) or molecules associated with other immunomodulatory mechanisms (such as regulation of macrophage activation). Some of the genes (n = 5) appeared to be under positive selection (Ka/Ks > 1), imitating proteins associated with cytokine signaling (through sequence homologies with thrombospondin type 1, toll/interleukin-1 receptor, TGF-β receptor and CD147).

CONCLUSIONS

With a comparative transcriptome profiling approach, we have identified a number of potential immunomodulator genes of F. hepatica (n = 62), which are firstly described here, could be employed for the development of better strategies (including RNAi) in the battle against both zoonotically and economically important disease, fasciolosis.

[Serotonergic genes in the development of anxiety/depression-like state and pathology of aggressive behavior in male mice: RNA-seq data]

In course of daily agonistic interactions, mice tend to stratify into those with chronic social defeats and those that repeatedly display aggression, which lead to the development of mixed anxiety/depression-like state and the pathology of aggressive behavior, respectively. Using the data of whole transcriptome analysis (RNA-seq), the changes in the expression of serotonergic genes involved in the synthesis, inactivation, and reception of serotonin, as well as of the Creb1 (transcription factor) gene and the Bdnf (brain-derived neurotrophic factor) gene were detected in the striatum (STR), ventral tegmental area (VTA), midbrain raphe nuclei (MRN), hypothalamus (HYP), and hippocampus (HIP) of defeated and aggressive male mice. In mice of both groups, the Tph2, Ddc, Slc6a4, Htr2a, Htr3a, Htr5b, Slc18a2, and Bdnf genes were downregulated in the MRN and the Tph2, Ddc, and Slc6a4 genes were upregulated in the VTA. These changes were more significant in defeated mice. The Htr5b gene has first been shown to be involved in mechanisms of depression and pathology of aggressive behavior. In the defeated mice, the expression levels of the Htr4 and Aldh1b1 genes were increased in the MRN, and expression levels of the Maob, Htr4, Htr1a, and Slc18a2 genes were increased in the VTA, while the expression level of the Htr3a gene was decreased. In the HYP of aggressive mice the Maoa, Htr2a, Htr2c, and Creb1 genes were downregulated and the Htr6 gene was upregulated. In the defeated mice, the Maoa and Creb1 genes were downregulated and the Htr6 and Aldh1b1 genes were upregulated in the HYP. In the STR, the Htr1a gene was downregulated and the Htr7 and Bdnf genes were upregulated. The Htr1b gene was upregulated in the HIP. The coexpression of dopaminergic and serotonergic genes in the MRN and VTA in the control of pathological behaviors is discussed. Thus, the complex pattern of differential expression of serotonergic genes in brain regions developing under repeated agonistic interactions in mice in dependence on behavioral pathology have been observed.

Cryo-Treatment Enhances the Embryogenicity of Mature Somatic Embryos via the lncRNA-miRNA-mRNA Network in White Spruce

In conifers, somatic embryogenesis is uniquely initiated from immature embryos in a narrow time window, which is considerably hindered by the difficulty to induce embryogenic tissue (ET) from other tissues, including mature somatic embryos. In this study, the embryogenic ability of newly induced ET and DNA methylation levels was detected, and whole-transcriptome sequencing analyses were carried out. The results showed that ultra-low temperature treatment significantly enhanced ET induction from mature somatic embryos, with the induction rate from 0.4% to 15.5%, but the underlying mechanisms remain unclear. The newly induced ET showed higher capability in generating mature embryos than the original ET. DNA methylation levels fluctuated during the ET induction process. Here, WGCNA analysis revealed that OPT4, TIP1-1, Chi I, GASA5, GST, LAX3, WRKY7, MYBS3, LRR-RLK, PBL7, and WIN1 genes are involved in stress response and auxin signal transduction. Through co-expression analysis, lncRNAs MSTRG.505746.1, MSTRG.1070680.1, and MSTRG.33602.1 might bind to pre-novel_miR_339 to promote the expression of WRKY7 genes for stress response; LAX3 could be protected by lncRNAs MSTRG.1070680.1 and MSTRG.33602.1 via serving as sponges for novel_miR_495 to initiate auxin signal transduction; lncRNAs MSTRG.505746.1, MSTRG.1070680.1, and MSTRG.33602.1 might serve as sponges for novel_miR_527 to enhance the expression of Chi I for early somatic embryo development. This study provides new insight into the area of stress-enhanced early somatic embryogenesis in conifers, which is also attributable to practical applications.

Inhibition of MicroRNA 6937 Delays Photoreceptor and Vision Loss in a Mouse Model of Retinitis Pigmentosa

Inherited retinal dystrophies (IRDs) are a group of rare retinal conditions, including retinitis pigmentosa (RP), caused by monogenic mutations in 1 out of more than 250 genes. Despite recent advancements in gene therapy, there is still a lack of an effective treatment for this group of retinal conditions. MicroRNAs (miRNAs) are a class of highly conserved small non-coding RNAs that inhibit gene expression. Control of miRNAs-mediated protein expression has been described as a widely used mechanism for post-transcriptional regulation in many physiological and pathological processes in different organs, including the retina. Our main purpose was to test the hypothesis that modulation of a group of miRNAs can protect photoreceptor cells from death in the rd10 mouse model of retinitis pigmentosa. For this, we incorporated modulators of three miRNAs in adeno-associated viruses (AAVs), which were administered through sub-retinal injections. The results obtained indicate that inhibition of the miR-6937-5p slows down the visual deterioration of rd10 mice, reflected by an increased electroretinogram (ERG) wave response under scotopic conditions and significant preservation of the outer nuclear layer thickness. This work contributes to broadening our knowledge on the molecular mechanisms underlying retinitis pigmentosa and supports the development of novel therapeutic approaches for RP based on miRNA modulation.

Tissue Microbiome Associated With Human Diseases by Whole Transcriptome Sequencing and 16S Metagenomics

In recent years, a substantial number of tissue microbiome studies have been published, mainly due to the recent improvements in the minimization of microbial contamination during whole transcriptome analysis. Another reason for this trend is due to the capability of next-generation sequencing (NGS) to detect microbiome composition even in low biomass samples. Several recent studies demonstrate a significant role for the tissue microbiome in the development and progression of cancer and other diseases. For example, the increase of the abundance of Proteobacteria in tumor tissues of the breast has been revealed by gene expression analysis. The link between human papillomavirus infection and cervical cancer has been known for some time, but the relationship between the microbiome and breast cancer (BC) is more novel. There are also recent attempts to investigate the possible link between the brain microbiome and the cognitive dysfunction caused by neurological diseases. Such studies pointing to the role of the brain microbiome in Huntington’s disease (HD) and Alzheimer’s disease (AD) suggest that microbial colonization is a risk factor. In this review, we aim to summarize the studies that associate the tissue microbiome, rather than gut microbiome, with cancer and other diseases using whole-transcriptome analysis, along with 16S rRNA analysis. After providing several case studies for each relationship, we will discuss the potential role of transcriptome analysis on the broader portrayal of the pathophysiology of the breast, brain, and vaginal microbiome.

Integrated analysis of circRNA, lncRNA, miRNA and mRNA to reveal the ceRNA regulatory network of postnatal skeletal muscle development in Ningxiang pig

Introduction: The development of skeletal muscle is regulated by regulatory factors of genes and non-coding RNAs (ncRNAs). Methods: The objective of this study was to understand the transformation of muscle fiber type in the longissimus dorsi muscle of male Ningxiang pigs at four different growth stages (30, 90, 150, and 210 days after birth, n = 3) by histological analysis and whole transcriptome sequencing. Additionally, the study investigated the expression patterns of various RNAs involved in muscle fiber transformation and constructed a regulatory network for competing endogenous RNA (ceRNA) that includes circular RNA (circRNA)/long non-coding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA). Results: Histomorphology analysis showed that the diameter of muscle fiber reached its maximum at 150 days after birth. The slow muscle fiber transformation showed a pattern of initial decrease followed by an increase. 29,963 circRNAs, 2,683 lncRNAs, 986 miRNAs and 22,411 mRNAs with expression level ≥0 were identified by whole transcriptome sequencing. Furthermore, 642 differentially expressed circRNAs (DEc), 505 differentially expressed lncRNAs (DEl), 316 differentially expressed miRNAs (DEmi) and 6,090 differentially expressed mRNAs (DEm) were identified by differential expression analysis. Functions of differentially expressed mRNA were identified by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). GO enrichment analysis indicates that 40 known genes and 6 new genes are associated with skeletal muscle development. Additionally, KEGG analysis shows that these genes regulate skeletal muscle development via MAPK, FoxO, Hedgehog, PI3K-Akt, Notch, VEGF and other signaling pathways. Through protein-protein interaction (PPI) and transcription factor prediction (TFP), the action mode of skeletal muscle-related genes was explored. PPI analysis showed that there were stable interactions among 19 proteins, meanwhile, TFP analysis predicted 22 transcription factors such as HMG20B, MYF6, MYOD1 and MYOG, and 12 of the 19 interacting proteins were transcription factors. The regulatory network of ceRNA related to skeletal muscle development was constructed based on the correlation of various RNA expression levels and the targeted binding characteristics with miRNA. The regulatory network included 31 DEms, 59 miRNAs, 667 circRNAs and 224 lncRNAs. conclusion: Overall, the study revealed the role of ceRNA regulatory network in the transformation of skeletal muscle fiber types in Ningxiang pigs, which contributes to the understanding of ceRNA regulatory network in Ningxiang pigs during the skeletal muscle development period.

Whole-Transcriptome Sequencing Reveals a ceRNA Regulatory Network Associated with the Process of Periodic Albinism under Low Temperature in Baiye No. 1 (Camellia sinensis)

The young shoots of the tea plant Baiye No. 1 display an albino phenotype in the early spring under low environmental temperatures, and the leaves re-green like those of common tea cultivars during the warm season. Periodic albinism is precisely regulated by a complex gene network that leads to metabolic differences and enhances the nutritional value of tea leaves. Here, we identified messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) to construct competing endogenous RNA (ceRNA) regulatory networks. We performed whole-transcriptome sequencing of 12 samples from four periods (Bud, leaves not expanded; Alb, albino leaves; Med, re-greening leaves; and Gre, green leaves) and identified a total of 6325 differentially expressed mRNAs (DEmRNAs), 667 differentially expressed miRNAs (DEmiRNAs), 1702 differentially expressed lncRNAs (DElncRNAs), and 122 differentially expressed circRNAs (DEcircRNAs). Furthermore, we constructed ceRNA networks on the basis of co-differential expression analyses which comprised 112, 35, 38, and 15 DEmRNAs, DEmiRNAs, DElncRNAs, and DEcircRNAs, respectively. Based on the regulatory networks, we identified important genes and their interactions with lncRNAs, circRNAs, and miRNAs during periodic albinism, including the ceRNA regulatory network centered on miR5021x, the GAMYB-miR159-lncRNA regulatory network, and the NAC035-miR319x-circRNA regulatory network. These regulatory networks might be involved in the response to cold stress, photosynthesis, chlorophyll synthesis, amino acid synthesis, and flavonoid accumulation. Our findings provide novel insights into ceRNA regulatory mechanisms involved in Baiye No. 1 during periodic albinism and will aid future studies of the molecular mechanisms underlying albinism mutants.

Identification of Smad3-related transcriptomes in type-2 diabetic nephropathy by whole transcriptome RNA sequencing

Smad3 deficiency prevents the development of type 2 diabetic nephropathy; however, the underlying molecular mechanisms remain unknown. In this study, we aimed to identify Smad3‐related genes involved in the pathogenesis of diabetic kidney disease. High‐throughput RNA sequencing was performed to profile the whole transcriptome in the diabetic kidney of Smad3 WT‐db/db, Smad3 KO‐db/db, Smad3^+/− db/db and their littermate control db/m mice at 20 weeks. The gene ontology, pathways and alternative splicing of differentially expressed protein‐coding genes and long non‐coding RNAs related to Smad3 in diabetic kidney were analysed. Compared to Smad3 WT‐db/db mice, Smad3 KO‐db/db mice exhibited an alteration of genes associated with RNA splicing and metabolism, whereas heterozygosity deletion of Smad3 (Smad3^+/− db/db mice) significantly altered genes related to cell division and cell cycle. Notably, three protein‐coding genes (Upk1b, Psca and Gdf15) and two lncRNAs (NONMMUG023520.2 and NONMMUG032975.2) were identified to be Smad3‐dependent and to be associated with the development of diabetic nephropathy. By using whole transcriptome RNA sequencing, we identified novel Smad3 transcripts related to the development of diabetic nephropathy. Thus, targeting these transcripts may represent a novel and effective therapy for diabetic nephropathy.

Effect of apigenin on whole transcriptome profile of TNFα-activated MDA-MB-468 triple negative breast cancer cells

The lack of hormone receptors in triple negative breast cancer (TNBC) is associated with the inefficacy of anti-estrogen chemotherapies, leaving fewer options for patient treatment and higher mortality rates. Additionally, as with numerous types of inflammatory breast cancer, infiltration of tumor associated macrophages and other leukocyte sub-populations within the tumor inevitably lead to aggressive, chemo-resistant, metastatic and invasive types of cancer which escape immune surveillance. These processes are orchestrated by the release of potent cytokines, including TNFα, IL-6 and CCL2 from the stroma, tumor and immune cells within the tumor microenvironment. The present study evaluated apigenin modulating effects on the pro-inflammatory activating action of TNFα in TNBC MDA-MB-468 cells, derived from an African American woman. Initially, cell viability was determined to establish an optimal sub-lethal dose of TNFα and apigenin in MDA-MB-468 cells. Subsequently, various treatments effects were evaluated using whole transcriptomic analysis of mRNA and long intergenic non-coding RNA with Affymetrix HuGene-2.1-st human microarrays. Gene level differential expression analysis was conducted on 48,226 genes where TNFα caused significant upregulation of 53 transcripts and downregulation of 11 transcripts. The largest upward differential shift was for CCL2 [+61.86 fold change (FC); false discovery rate (FDR), P<0.0001]; which was down regulated by apigenin (to +10.71 FC vs. Control; FDR P-value <0.001), equivalent to an 83% reduction. Several TNFα deferentially upregulated transcripts were reduced by apigenin, including CXCL10, C3, PGLYRP4, IL22RA2, KMO, IL7R, ROS1, CFB, IKBKe, SLITRK6 (a checkpoint target) and MMP13. Confirmation of CCL2 experimentally induced transcript alterations was corroborated at the protein level by ELISA assays. The high level of CCL2 transcript in the cell line was comparable to that in our previous studies in MDA-MB-231 cells. The differential effects of TNFα were corroborated by ELISA, where the data revealed a >10-fold higher releasing rate of CCL2 in MDA-MB-468 cells compared with in MDA-MB-231 cells, both of which were attenuated by apigenin. The data obtained in the present study demonstrated a high level of CCL2 in MDA-MB-468 cells and a possible therapeutic role for apigenin in downregulating TNFα-mediated processes in these TNBC cells.

Lower cold tolerance of tropical Porites lutea is possibly detrimental to its migration to relatively high latitude refuges in the South China Sea

As high temperature stress due to climate change threatens tropical corals, cooler areas at relatively high latitudes may be potential refuges. Tolerance to low temperatures is critical in determining whether corals can successfully migrate to higher latitudes. However, the physiological and molecular adaptations that protect corals from low temperature stress are unclear. In this study, scleractinian Porites lutea samples from the tropical Xisha Islands (XS) and subtropical Daya Bay (DY) in the South China Sea were subjected to a reduction in ambient temperature from 26 to 12°C. Differences in physiological changes and gene expression were analysed. P. lutea from both XS and DY exhibited physiological bleaching under low temperature stress, and the Symbiodiniaceae density, Fv/Fm, and chlorophyll-α content were significantly reduced. Symbiosome antioxidative stress and metabolic enzyme activity first increased and then decreased. RNA-seq analysis showed that the host responded to low temperature stress by activating immune, apoptotic, and autophagic pathways and reducing metabolic levels. Nevertheless, Symbiodiniaceae lacked the physiological regulatory capacity to adapt to low temperatures. The lower cold tolerance of XS tropical P. lutea may attribute to lower oxidative stress resistance, lower photosynthetic capacity, worse energy supply, and higher susceptibility to bacterial and viral infections and diseases in XS corals. The difference in cold tolerance may result from genetic differences between the geographic populations and is possibly detrimental to the migration of tropical coral to relatively high latitude refuges. This study provides a theoretical basis for anthropogenically assisted coral migration as a response to global change.

Molecular Changes Induced in Melanoma by Cell Culturing in 3D Alginate Hydrogels

Simple Summary

The research field of 3D cell cultivation in hydrogels is continuously growing. To be able to analyze the reaction of melanoma cells to 3D cultivation in alginate hydrogel on a molecular level, whole transcriptome sequencing was performed. Intriguingly, we could not only unravel differences between the gene regulation in 2D and 3D cultures but could also correlate the culture switch to the physiological process of tumor plasticity based on the observed patterns. Thereby, the role of EGR1 in controlling tumor plasticity and progression in melanoma was revealed. We conclude that the combination of cell culture models using biomaterials and whole transcriptome analysis leads to a deeper molecular understanding of cancer cells, herewith defining new therapeutic targets.

Abstract

Alginate hydrogels have been used as a biomaterial for 3D culturing for several years. Here, gene expression patterns in melanoma cells cultivated in 3D alginate are compared to 2D cultures. It is well-known that 2D cell culture is not resembling the complex in vivo situation well. However, the use of very intricate 3D models does not allow performing high-throughput screening and analysis is highly complex. 3D cell culture strategies in hydrogels will better mimic the in vivo situation while they maintain feasibility for large-scale analysis. As alginate is an easy-to-use material and due to its favorable properties, it is commonly applied as a bioink component in the growing field of cell encapsulation and biofabrication. Yet, only a little information about the transcriptome in 3D cultures in hydrogels like alginate is available. In this study, changes in the transcriptome based on RNA-Seq data by cultivating melanoma cells in 3D alginate are analyzed and reveal marked changes compared to cells cultured on usual 2D tissue culture plastic. Deregulated genes represent valuable cues to signaling pathways and molecules affected by the culture method. Using this as a model system for tumor cell plasticity and heterogeneity, EGR1 is determined to play an important role in melanoma progression.

Emerging role of a systems biology approach to elucidate factors of reduced penetrance: transcriptional changes in THAP1-linked dystonia as an example

Pathogenic variants in THAP1 can cause dystonia with a penetrance of about 50 %. The underlying mechanisms are unknown and can be considered as means of endogenous disease protection. Since THAP1 encodes a transcription factor, drivers of this variability putatively act at the transcriptome level. Several transcriptome studies tried to elucidate THAP1 function in diverse cellular and mouse models, including mutation carrier-derived cells and iPSC-derived neurons, unveiling various differentially expressed genes and affected pathways. These include nervous system development, dopamine signalling, myelination, or cell-cell adhesion. A network diffusion analysis revealed mRNA splicing, mitochondria, DNA repair, and metabolism as significant pathways that may represent potential targets for therapeutic interventions.

Systematic Identification and Comparison of the Expressed Profiles of lncRNAs, miRNAs, circRNAs, and mRNAs with Associated Co-Expression Networks in Pigs with Low and High Intramuscular Fat

Intramuscular fat (IMF) content is a complex trait that affects meat quality and determines pork quality. In order to explore the potential mechanisms that affect the intramuscular fat content of pigs, a Large white × Min pigs F2 resource populations were constructed, then whole-transcriptome profile analysis was carried out for five low-IMF and five high-IMF F2 individuals. In total, 218 messenger RNA (mRNAs), 213 long non-coding RNAs (lncRNAs), 18 microRNAs (miRNAs), and 59 circular RNAs (circRNAs) were found to be differentially expressed in the longissimus dorsi muscle. Gene ontology analysis and Kyoto Encyclopedia of Genes and Genomes annotations revealed that these differentially expressed (DE) genes or potential target genes (PTGs) of DE regulatory RNAs (lncRNAs, miRNAs, and circRNAs) are mainly involved in cell differentiation, fatty acid synthesis, system development, muscle fiber development, and regulating lipid metabolism. In total, 274 PTGs were found to be differentially expressed between low- and high-IMF pigs, which indicated that some DE regulatory RNAs may contribute to the deposition/metabolism of IMF by regulating their PTGs. In addition, we analyzed the quantitative trait loci (QTLs) of DE RNAs co-located in high- and low-IMF groups. A total of 97 DE regulatory RNAs could be found located in the QTLs related to IMF. Co-expression networks among different types of RNA and competing endogenous RNA (ceRNA) regulatory networks were also constructed, and some genes involved in type I diabetes mellitus were found to play an important role in the complex molecular process of intramuscular fat deposition. This study identified and analyzed some differential RNAs, regulatory RNAs, and PTGs related to IMF, and provided new insights into the study of IMF formation at the level of the genome-wide landscape.

Differential age-related transcriptomic analysis of ovarian granulosa cells in Kazakh horses

Introduction

The Kazakh horse, renowned for its excellence as a breed, exhibits distinctive reproductive traits characterized by early maturity and seasonal estrus. While normal reproductive function is crucial for ensuring the breeding and expansion of the Kazakh horse population, a noteworthy decline in reproductive capabilities is observed after reaching 14 years of age.

Methods

In this study, ovarian granulosa cells (GCs) were meticulously collected from Kazakh horses aged 1, 2, 7, and above 15 years old (excluding 15 years old) for whole transcriptome sequencing.

Results

The analysis identified and selected differentially expressed mRNAs, lncRNAs, miRNAs, and circRNAs for each age group, followed by a thorough examination through GO enrichment analysis. The study uncovered significant variations in the expression profiles of mRNAs, lncRNAs, miRNAs, and circRNAs within GCs at different stages of maturity. Notably, eca-miR-486-3p and miR-486-y exhibited the highest degree of connectivity. Subsequent GO, KEGG, PPI, and ceRNA network analyses elucidated that the differentially expressed target genes actively participate in signaling pathways associated with cell proliferation, apoptosis, and hormonal regulation. These pathways include but are not limited to the MAPK signaling pathway, Hippo signaling pathway, Wnt signaling pathway, Calcium signaling pathway, Aldosterone synthesis and secretion, Cellular senescence, and NF-kappa B signaling pathway-essentially encompassing signal transduction pathways crucial to reproductive processes.

Discussion

This research significantly contributes to unraveling the molecular mechanisms governing follicular development in Kazakh horses. It establishes and preliminarily validates a differential regulatory network involving lncRNA-miRNA-mRNA, intricately associated with processes such as cell proliferation, differentiation, and apoptosis and integral to the developmental intricacies of stromal follicles. The findings of this study provide a solid theoretical foundation for delving deeper into the realm of reproductive aging in Kazakh mares, presenting itself as a pivotal regulatory pathway in the context of horse ovarian development.

Whole blood transcriptome profile identifies motor neurone disease RNA biomarker signatures

Blood-based biomarkers for motor neuron disease are needed for better diagnosis, progression prediction, and clinical trial monitoring. We used whole blood-derived total RNA and performed whole transcriptome analysis to compare the gene expression profiles in (motor neurone disease) MND patients to the control subjects. We compared 42 MND patients to 42 aged and sex-matched healthy controls and described the whole transcriptome profile characteristic for MND. In addition to the formal differential analysis, we performed functional annotation of the genomics data and identified the molecular pathways that are differentially regulated in MND patients. We identified 12,972 genes differentially expressed in the blood of MND patients compared to age and sex-matched controls. Functional genomic annotation identified activation of the pathways related to neurodegeneration, RNA transcription, RNA splicing and extracellular matrix reorganisation. Blood-based whole transcriptomic analysis can reliably differentiate MND patients from controls and can provide useful information for the clinical management of the disease and clinical trials.

Cadmium Exposure Alters Rhizospheric Microbial Community and Transcriptional Expression of Vetiver Grass

Vetiver grass (Chrysopogon zizanioides L.) has been used to remediate cadmium (Cd)-contaminated soil, while there have been few studies on the influence of Cd exposure on the rhizospheric microbial community and transcriptional expression of C. zizanioides. In this study, we investigated the response of the rhizospheric microbial community and transcriptional expression of C. zizanioides in 20 mg/kg Cd-contaminated soil. The results showed that Cd levels in the roots and shoots of C. zizanioides reached 250.80 and 73.40 mg/kg, respectively. The Cd exposure changed the rhizospheric bacterial community, resulting in the significant enrichment of Sphingomonas, Lysobacter, and Gemmatimonadetes in Cd-contaminated soil. In addition, 880 and 3,419 differentially expressed genes were identified in the plant roots and shoots, respectively, in response to Cd stress. Among these, the overexpressed genes associated with redox homeostasis, glutathione (GSH) metabolism, cell wall biosynthesis, and transmembrane transport pathways were found to participate in Cd detoxification in C. zizanioides. These findings could be useful for understanding the selective variation of the rhizospheric microbial community and the detoxification mechanisms of C. zizanioides in Cd phytoremediation.

Spatially Resolved in vivo CRISPR Screen Sequencing via Perturb-DBiT

Perturb-seq enabled the profiling of transcriptional effects of genetic perturbations in single cells but lacks the ability to examine the impact on tissue environments. We present Perturb-DBiT for simultaneous co-sequencing of spatial transcriptome and guide RNAs (gRNAs) on the same tissue section for in vivo CRISPR screen with genome-scale gRNA libraries, offering a comprehensive understanding of how genetic modifications affect cellular behavior and tissue architecture. This platform supports a variety of delivery vectors, gRNA library sizes, and tissue preparations, along with two distinct gRNA capture methods, making it adaptable to a wide range of experimental setups. In applying Perturb-DBiT, we conducted un-biased knockouts of tens of genes or at genome-wide scale across three cancer models. We mapped all gRNAs in individual colonies and corresponding transcriptomes in a human cancer metastatic colonization model, revealing clonal dynamics and cooperation. We also examined the effect of genetic perturbation on the tumor immune microenvironment in an immune-competent syngeneic model, uncovering differential and synergistic perturbations in promoting immune infiltration or suppression in tumors. Perturb-DBiT allows for simultaneously evaluating the impact of each knockout on tumor initiation, development, metastasis, histopathology, and immune landscape. Ultimately, it not only broadens the scope of genetic inquiry, but also lays the groundwork for developing targeted therapeutic strategies.

Whole Transcriptome-Based ceRNA Regulatory Network Analysis of Radiation-Induced Esophageal Epithelial Cell Injury

Introduction

Esophageal epithelial cells are essential for esophageal homeostasis and defense against harmful stimuli, but the mechanisms of radiation-induced injury in these cells are poorly understood. The competitive endogenous RNA (ceRNA) network, involved in various physiological processes and diseases, may also play a role in radiation-induced injury, although its mechanism remains unclear. This study aimed to investigate the effects of ionizing radiation on human esophageal epithelial cells and explore the role of the ceRNA network in this injury.

Methods

Cellular phenotype experiments assessed the effects of ionizing radiation on human esophageal epithelial cells. Whole transcriptome sequencing (lncRNA, circRNA, miRNA, and mRNA) was performed on cells exposed to 0, 2, and 4 Gy radiation. Differentially expressed RNAs (dd-DERs) were identified through differential expression analysis and dose-dependent screening. A ceRNA network was constructed using co-expression analysis and binding site prediction. Real-time quantitative PCR validated the expression levels of selected dd-DERs, and gene set enrichment analysis explored affected pathways.

Results

We identified 41 lncRNAs, 18 miRNAs, and 192 mRNAs as dose-dependent differentially expressed RNAs. A ceRNA network comprising 10 lncRNAs, 5 miRNAs, and 55 mRNAs was established. Real-time PCR confirmed the expression levels of 8 dd-DERs within the network. Gene set enrichment analysis showed that radiation disrupted channel activity, cell replication, repair, and immune response. Functional enrichment analysis revealed modulation of metabolic pathways, particularly involving UGT1A family members.

Discussion

This study established a ceRNA network related to radiation-induced esophageal epithelial cell injury, advancing our understanding of its pathophysiology. The ceRNA network may mediate injury through metabolic pathway modulation. Future work should focus on elucidating specific ceRNA interactions and exploring therapeutic potential for mitigating radiation-induced esophageal injury.

High-throughput sequencing reveals the existence of a novel fucose binding lectin gene from the whole gut and fat body tissues of the grub of darkling beetle, Zophobas morio

Fucose binding lectins are a taxonomically, evolutionarily, and ecologically relevant class of lectins that have been identified among a wide range of taxa. These lectins possess a distinctive ability to identify and bind fucose-containing glycans, thereby lending them the ability to function as immune cell mitogens, diagnostic markers and anti-cancer agents. Our preliminary analysis revealed the existence of a single D-fucose binding lectin from the grub serum of the darkling beetle, Z. morio. Here, whole transcriptome analysis using entire gut and fat body tissues of grub of Z. morio with Illumina NovaSeq6000 sequencing platform revealed the existence of a novel fucose binding lectin (ZmFBL) displaying significant similarity to Diaphorina citri fucolectin-3 and possessing fucolectin-tachylectin4-pentraxin domain. Molecular docking between the structure of the predicted FBL transcript from Z. morio and D-fucose demonstrated the lowest binding energy of -5.4 kcal/mol. The MD simulation and MM/PBSA analysis furnished insights into the binding stability of D-fucose and the ZmFBL. A similar transcript was also identified from the fat body transcriptome. A 2.49-fold increase in the ZmFBL expression was observed in the fat body than the whole gut as evidenced from the relative quantitation using RT-PCR.