Transcriptome assembly and functional annotation are essential for understanding gene expression and biological function. Nevertheless, many existing tools lack the flexibility to integrate both short- and long-read sequencing data or fail to provide a complete, customizable workflow for transcriptome analysis. Here, we present TrAnnoScope, a comprehensive transcriptome analysis pipeline that provides a complete, customizable workflow capable of efficiently processing and integrating short- and long-read sequencing data to generate high-quality, full-length transcripts with detailed functional annotation. The pipeline encompasses steps from quality control to functional annotation, employing a range of tools and established databases, such as SwissProt, Pfam, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes. As a case study, TrAnnoScope was applied to RNA-Seq data from zebra finch brain, ovary, and testis tissues. TrAnnoScope produced comprehensive transcriptome, demonstrating strong alignment with the reference genome (99.63%) and capturing a significant proportion of nearly complete protein sequences (92.7%). The functional annotation process yielded extensive matches to known protein databases and successfully assigned relevant functional terms to majority of the transcripts. As such, TrAnnoScope successfully integrates multiple sequencing technologies to generate comprehensive transcriptomes with minimal user input. Its modular design, flexibility, and ease of use make it a valuable tool for researchers working with complex datasets, particularly for non-model organisms.

Circular RNAs are a special type of RNAs which recently attracted a lot of research interest in studying its formation and function. RNA binding proteins (RBPs) that bind circRNAs are important in these processes but are relatively less studied. CLIP-Seq technology has been invented and applied to profile RBP-RNA interactions on the genome-wide scale. While mRNAs are usually the focus of CLIP-Seq experiments, RBP-circRNA interactions could also be identified through specialized analysis of CLIP-Seq datasets. However, many technical difficulties are involved in this process, such as the usually short read length of CLIP-Seq reads. In this study, we created a pipeline called Clirc specialized for profiling circRNAs in CLIP-Seq data and analyzing the characteristics of RBP- circRNAs interactions. In conclusion, this is one of the first few studies to investigate circRNAs and their binding partners through repurposing CLIP-Seq datasets to our knowledge, and we hope our work will become a valuable resource for future studies into the biogenesis and function of circRNAs. Clirc software is available at https://github.com/Minzhe/Clirc

Long non-coding RNAs (lncRNAs) have been found to be involved in many biological processes, including the regulation of cell differentiation, but a complete characterization of lncRNA is still lacking. Additionally, there is evidence that lncRNAs interact with ribosomes, raising questions about their functions in cells. Here, we used a developmentally staged protocol to induce cardiogenic commitment of hESCs and then investigated the differential association of lncRNAs with polysomes. Our results identified lncRNAs in both the ribosome-free and polysome-bound fractions during cardiogenesis and showed a very well-defined temporal lncRNA association with polysomes. Clustering of lncRNAs was performed according to the gene expression patterns during the five timepoints analyzed. In addition, differential lncRNA recruitment to polysomes was observed when comparing the differentially expressed lncRNAs in the ribosome-free and polysome-bound fractions or when calculating the polysome-bound vs ribosome-free ratio. The association of lncRNAs with polysomes could represent an additional cytoplasmic role of lncRNAs, e.g., in translational regulation of mRNA expression.

A definitive transcriptome atlas for the non-coding expressed elements of pathogenic mycobacteria does not exist. Incomplete lists of non-coding transcripts can be obtained for some of the reference genomes (e.g. Mycobacterium tuberculosis H37Rv) but to what extent these transcripts have homologues in closely related species or even strains is not clear. This has implications for the analysis of transcriptomic data; non-coding parts of the transcriptome are often ignored in the absence of formal, reliable annotation. Here, we review the state of our knowledge of non-coding RNAs in pathogenic mycobacteria, emphasising the disparities in the information included in commonly used databases. We then proceed to review ways of combining computational solutions for predicting the non- coding transcriptome with experiments that can help refine and confirm these predictions.

To guide analysts to select the right tool and parameters in differential gene expression analysis of single-cell RNA sequencing (scRNA-seq) data, we developed a novel simulator that recapitulates the data characteristics of real scRNA-seq datasets while accounting for all the relevant sources of variation in a multi-subject, multi-condition scRNA-seq experiment: the cell-to-cell variation within a subject, the variation across subjects, the variability across cell types, the mean/variance relationship of gene expression across genes, library size effects, group effects, and covariate effects. By applying it to benchmark 12 differential gene expression analysis methods (including cell-level and pseudo-bulk methods) on simulated multi-condition, multi-subject data of the 10x Genomics platform, we demonstrated that methods originating from the negative binomial mixed model such as glmmTMB and NEBULA-HL outperformed other methods. Utilizing NEBULA-HL in a statistical analysis pipeline (https://github.com/interactivereport/scRNAseq_DE) for single cell analysis will enable scientists to better understand cell-type specific transcriptomic response to disease or treatment effects and to discover new drug targets. Further, application to two real datasets showed the outperformance of our differential expression (DE) pipeline, with unified findings of differentially expressed genes (DEG) and a pseudo-time trajectory transcriptomic result. In the end, we made recommendations of filtering strategies of cells and genes based on simulation results to achieve optimal experimental goals.

The exogenous application of phenolic compounds is increasingly recognized as a valuable strategy for promoting growth and mitigating the adverse effects of abiotic stress. However, the biostimulant effect under optimal conditions have not been thoroughly explored. In this study, we investigated the impact of foliar application of flavonoids, specifically Crop Bio Life (CBL), on tomato plants grown under controlled conditions. Our study focused on determining growth parameters, such as cell size, and assessing the concentration of hormones. Principal component analysis (PCA) from all physiological variables were determined. Additionally, we utilized high-throughput mRNA sequencing technology and bioinformatic methodologies to robustly analyze the transcriptomes of tomato leaves regulated by flavonoids. The findings revealed that CBL primarily influenced cells enlargement, leading to increased growth. Furthermore, CBL-treated plants exhibited higher concentrations of the hormone zeatin but lower concentrations of IAA. Moreover, RNAseq analysis indicated that CBL-treated plants required increased mineral transport and water uptake, as evidenced by gene expression patterns. Genes related to pathways such as fatty acid degradation, phenylpropanoid biosynthesis, and ABC transporters showed regulatory mechanisms governing internal flavonoid biosynthesis, transport, and tissue concentration, ultimately resulting in higher flavonoid concentrations in tomato leaves.

Plant life-long organogenesis involves sequential, time and tissue specific expression of developmental genes. This requires activities of Polycomb Group (PcG) and trithorax Group complexes, respectively responsible for repressive Histone 3 trimethylation at lysine 27 (H3K27me3) and activation-related H3K4me3. However, the genome-wide dynamics in histone modifications that occur during developmental processes have remained elusive. Here, we report the distributions of H3K27me3 and H3K4me3 along with transcriptional changes, in a developmental series including Arabidopsis leaf and three stages of flower development. We found that chromatin mark levels are highly dynamic over the time series on nearly half of all Arabidopsis genes. Moreover, during early flower morphogenesis, changes in H3K4me3 prime over changes in H3K27me3 and quantitatively correlate with transcription changes, while H3K27me3 changes occur after prolonged expression changes. Notably, early activation of PcG target genes is dominated by increases in H3K4me3 while H3K27me3 remains present at the locus. Our results reveal H3K4me3 as greater predictor over H3K27me3 for transcription dynamics, unveil unexpected chromatin mechanisms at gene activation and underline the relevance of tissue-specific temporal epigenomics.

Prostate cancer (Pca) is a highly heterogeneous disease and the second more common tumor in males. Molecular and genetic profiles have been used to identify subtypes and guide therapeutic intervention. However, roughly 26% of primary Pca are driven by unknown molecular lesions. We use Principal Component Analysis (PCA) and custom RNAseq-data normalization to identify a gene expression signature which segregates primary PRAD from normal tissues. This Core-Expression Signature (PRAD-CES) includes 33 genes and accounts for 39% of data complexity along the PC1-cancer axis. The PRAD-CES is populated by protein-coding (AMACR, TP63, HPN) and RNA-genes (PCA3, ARLN1) sparsely found in previous studies, validated/predicted biomarkers (HOXC6, TDRD1, DLX1), and/or cancer drivers (PCA3, ARLN1, PCAT-14). Of note, the PRAD-CES also comprises six over-expressed LncRNAs without previous Pca association, four of them potentially modulating driver’s genes TMPRSS2, PRUNE2 and AMACR. Overall, our PCA capture 57% of data complexity within PC1-3. GO enrichment and correlation analysis involving major clinical features (i.e., Gleason Score, AR Score, TMPRSS2-ERG fusion and Tumor Cellularity) suggest that PC2 and PC3 gene signatures might describe more aggressive and inflammation-prone transitional forms of PRAD. Of note, surfaced genes may entail novel prognostic biomarkers and molecular alterations to intervene. Particularly, our work uncovered RNA genes with appealing implications on Pca biology and progression.

Mulberries is a type of fruit tree that is rich in nutrients and possesses both medicinal and dietary properties, deeply loved by consumers. The fragrance is a significant factor in the sensory quality of mulberries, playing a crucial role in consumer purchasing behavior. 'Baichang' mulberry is a unique variety of mulberry that emits a strong, creamy-like fragrance, which has enormous potential market value. However, the key metabolites and mechanisms responsible for the formation of this fragrance are still unclear. In this study, we utilized GC-MS technology and RNA-seq to examine the Volatile Organic Compounds (VOCs) and gene expression in 'Baichang' (Morus macroura) mulberries and 'Hongchang' (Morus laevigata) mulberries, the latter of which do not have a distinct fragrance. The results showed that 661 volatile metabolites were detected in both mulberry varieties. Most of these were found in higher relative amounts in the 'Baichang' mulberry variety, with 47 types being present only in the 'Baichang' mulberries, suggesting a rich composition of VOCs in the 'Baichang' mulberry. 312 metabolites were identified as differential metabolites; combining the sensory annotation information and relative content of VOCs, 70 candidate VOCs were annotated as related to the formation of 'sweet, waxy, coconut, and oily' sensations. They could potentially be the key metabolites contributing to the unique aroma of 'Baichang' mulberries. The relative odor activity value (rOAV) analysis revealed that 2(3H)-Furanone, 5-hexyldihydro- was the key VOC contributing to the formation of the special aroma of 'Baichang' mulberries. By RNA-seq, 19 key structural genes related to ester synthesis were identified, including AAT , ACX , ALDH, CYP, EHL, FAD, HPL, and LOX. The high expression of these structural genes in the "Baichang" mulberry fruit is the main reason for its unique creamy aroma. This study provides insights into the aroma differences among different kinds of mulberrie.

This study explored the effects of pBD114 (porcine β-defensin 114) on inflammatory response. The results in RAW264.7 cells showed that 100 μg/mL of pBD114 significantly increased the concentrations of TNFα and IL-10, and 10 μg/mL of pBD114 significantly decreased the concen-tration of TNFα in LPS-induced RAW264.7 cell cultures at 10 ng/mL. These results suggest that pBD114 can exhibit pro-inflammatory or anti-inflammatory activities under certain conditions. RNA-seq analysis was performed to gain further insight into the effects of pBD114 on the in-flammatory response. Among pBD114 promoting RAW264.7 pro-inflammatory response, pBD114 significantly up-regulated 1170 genes and down-regulated 724 genes. KEGG enrichment found that the DEGs were significantly enriched in immune- and signal-transduction-related signaling pathways. PPI and KDA analyses revealed that Bcl10 and Bcl3 were the key genes. In addition, pBD114 significantly up-regulated 12 genes and down-regulated 38 genes in the an-ti-inflammatory response. KEGG enrichment analysis revealed that the DEGs were mainly en-riched in the "Cytokine-cytokine receptor interaction" signaling pathway, and PPI and KDA analyses showed that Stat1 and Csf2 were the key genes. The results of qRT-PCR verified those of RNA-seq. In vivo mice tests also confirmed the pro- or anti-inflammatory activity of pBD114. Although the inflammatory response is a rapid and complex physiological reaction to noxious stimuli, this study found that pBD114 plays an essential role in the inflammatory response mainly by acting on genes related to immunity, signal transduction, signaling molecules, and interactions.

Cervical cancer presents a significant global health concern with high-risk human papilloma-viruses (HPVs) being identified as the cause of this cancer. Although current treatment methods for cervical cancer can eliminate lesions, preventing metastatic spread and minimizing tissue damage remains a major challenge. Therefore, the development of a safer and innovative thera-peutic approach is of utmost importance. Natural products like Fig latex, derived from the Ficus Carica tree, have demonstrated promising anti-cancer properties when tested on cervical cancer cell lines. However, the specific mechanisms by which Fig latex exerts its effects are still un-known. In this study, we conducted RNA-seq analysis to explore how Fig latex may counteract carcinogenesis in HPV-positive cervical cancer cell lines, namely CaSki (HPV type 16-positive) and HeLa (HPV type 18-positive). Our data from this investigation indicates that Fig latex influ-ences the expression of genes associated with the development and progression of cervical can-cer, including pathways related to " Nonsense-Mediated Decay (NMD) "cell cycle regulation" and Transcriptional Regulation by TP53. This selective impact of Fig latex on cancer-related pathways suggests a potential novel therapeutic approach for HPV-related cervical cancer.

Alternative splicing can produce transcripts that affect cancer development and shows potential for cancer diagnosis and treatment. However, intron retention (IR), a type of alternative splicing, has been less systematically studied in cancer biology research. Here, we generated a pan-cancer IR landscape for more than 10,000 samples across 33 cancer types from The Cancer Genome Atlas (TCGA). We characterized differentially retained introns between tumor and normal samples and identified retained introns associated with survival. We discovered 988 differentially retained introns in 14 cancers, some of which demonstrated diagnostic potential in multiple cancer types. We also inferred a large number of prognosis-related introns in 33 cancer types, and the associated genes included well-known cancer hallmarks such as angiogenesis, metastasis, and DNA mutations. Notably, we discovered a novel intron retention event inside 5′UTR of STN1 that is associated with the survival of lung cancer patients. The retained intron reduces translation efficiency by producing upstream open reading frames (uORFs) and thereby inhibits colony formation and cell migration of lung cancer cells. Besides, the IR-based prognostic model achieved good stratification on certain cancers, as illustrated in acute myeloid leukemia. Taken together, we performed a comprehensive IR survey at a pan-cancer level, and the results implied that IR has the potential to be diagnostic and prognostic cancer biomarkers, as well as new drug targets.

The transcriptomic analysis of the dehydration rate of mature rice seeds was conducted to explore candidate genes related to the dehydration rate and provide a theoretical basis for breeding and utilization. In this study, we selected Baghlani Nangarhar, an extremely fast dehydration material, and Saturn, a slow dehydration material, based on the results of the previous studies of screening of 165 germplasm materials for dehydration rate phenotypes. Fast dehydration experiment on these two types of seeds was conducted. Four comparative groups were set up under control and dehydration conditions. The differentially expressed genes (DEGs) were quantified by transcriptome sequencing and quantitative real-time PCR (qRT-PCR). GO and KEGG analyses were carried out. The results showed that in Baghlani Nangarhar, 53 DEGs were screened, of which 33 were up-regulated and 20 were down-regulated. Similarly, in Saturn, 25 DEGs were screened, of which 19 were up-regulated and 6 down-regulated. The results of the GO analysis showed that the sites of action of the differentially expressed genes enriched in the fast dehydration modes were concentrated in the cytoplasm, internal components of the membrane and nucleosomes. They play regulatory roles in catalysis, binding, translocation, transcription, protein folding, degradation and replication. They are involved in adaptive responses to adverse external environments such as reactive oxygen species and high temperature. The KEGG analysis showed that the main metabolic pathways enriched are protein processing in the endoplasmic reticulum, amino acid biosynthesis, and oxidative phosphorylation. The key differentially expressed genes and the most important metabolic pathways in the rapidly and slowly dehydrated materials obtained in this study were protein processing in the endoplasmic reticulum and oxidative phosphorylation metabolism, which are presumed to have important regulatory roles in stress/defense, energy metabolism, protein synthesis/folding, and signal transduction during dehydration and drying of mature seeds. The results of this study could potentially provide a valuable reference for further research on the genes and metabolic pathways related to the dehydration rate of mature rice seeds and provide theoretical guidance for the selection and breeding of new rice germplasm that can be rapidly dehydrated at the mature stage.

Caseous lymphadenitis (CL) is a chronic contagious disease that affects small ruminants and is characterized by the formation of pyogranulomas in lymph nodes and other organs. However, the pathogenesis of this disease and the response of the host genome to infection are not yet fully understood. This study aimed to investigate the whole blood transcriptome and evaluate differential gene expression during the later stages of CL in naturally infected ewes. The study included diseased, serologically positive (EP), exposed, serologically negative (EN) ewes from the same infected flock and healthy ewes (CN) from a different flock. RNA sequencing was performed using the Illumina NextSeq system, and differential gene expression was estimated using DESeq2 and Edge R approaches. The analysis identified 191 annotated differentially expressed genes (DEGs) in the EP group (102 upregulated and 89 downregulated) and 256 DEGs in the EN group (106 upregulated and 150 downregulated). Infection influenced numerous immunoregulatory interactions between lymphoid and nonlymphoid cells in both EP and EN ewes. Immune DEGs were preferentially assigned to antigen presentation through the MHC complex, T-lymphocyte mediated immunity, and extracellular matrix interactions. Furthermore, the EP group showed altered regulation of cytokine and chemokine signalling and activation and recombination of B-cell receptors. Conversely, NF-kappa B signalling, apoptosis, and stress response were the main processes influenced in the EN group. In addition, statistically significant enrichment of the essential immune pathways of Binding and Uptake of Ligands by Scavenger Receptors in EP and p53 signalling in the EN group was found. In conclusion, this study provides new insights into the disease course and host‒pathogen interaction in naturally CL-infected sheep by investigating the blood transcriptome.

Brown-marbled grouper is one of the most important mariculture species in China, which is used as an important crossbreeding parent in grouper industry. Enhancing growth rates is a key target in fish breeding, and gaining insight into the underlying mechanisms responsible for growth differences among individuals can aid in the improvement of grouper growth rates. However, the mechanism behind this difference in growth in this fish is unclear. The difference of transcriptome profiles of muscle tissue between fast- and slow-growing brown-marbled grouper was analyzed by RNA-Seq. 77 significantly up-regulated genes and 92 significantly down-regulated genes were identified in the growth extreme groups. The up-regulated of ghr and tnni2 and the down-regulated of stc2 led to the growth advantages of brown-marbled grouper. The differently expressed genes (DEGs) were used for GO and KEGG enrichment analysis. The results of GO enrichments indicated that the significantly upregulated genes in the fast-growing group were involved in protein folding, actin cytoskeleton, myosin complex, etc. The results of KEGG enrichments indicated that the significantly upregulated genes in the fast-growing group were involved in glycolysis/ gluconeogenesis, adipocytokine signaling pathway, MAPK signaling pathway, carbon metabolism, PI3K-Akt signaling pathway, etc. To analysis the difference gene sets between fast- and slow-growing group, the RNA-seq data were used by gene set enrichment analysis (GSEA). The results showed that the PI3K/AKT/mTOR pathway was up-regulated in the fast-growing group. The up-regulated of this pathway could lead to higher nutrient absorption efficiency and lead to muscle growth in the fast-growing group. These results contribute to understanding of the molecular mechanisms of fast growth and regulative pathways regulating growth in brown-marbled grouper.

Research Highlights: This study identified the cell cycle genes in birch that likely play important roles during plant growth and development. This analysis provides a basis for understanding the regulatory mechanism of various cell cycles in Betula pendula. Background and Objectives: The cell cycle factors not only influence cell cycle progression together, but also regulate accretion, division and differentiation of cells, and then regulate growth and development of plant. In this study, we identified the putative cell cycle genes in B. pendula genome, based on the annotated cell cycle genes in A. thaliana. It could serve as a foundation for further functional studies. Materials and Methods: The transcript abundance was determined for all the cell cycle genes in xylem, root, leaf and flower tissues using RNA-seq technology. Results: We identified 59 cell cycle gene models in the genome of B. pendula, 17 highly expression genes among them. These genes were BpCDKA.1, BpCDKB1.1, BpCDKB2.1, BpCKS1.2, BpCYCB1.1, BpCYCB1.2, BpCYCB2.1, BpCYCD3.1, BpCYCD3.5, BpDEL1, BpDpa2, BpE2Fa, BpE2Fb, BpKRP1, BpKRP2, BpRb1 and BpWEE1. Conclusions: We identified 17 core cell cycle genes in the genome of birch by combining phylogenetic analysis and tissue specific expression data.

RNA sequencing (RNA-Seq) is a complicated protocol, both in the laboratory in generation of data and at the computer in analysis of results. Several decisions during RNA-Seq library construction have important implications for analysis, most notably strandedness during complementary DNA (cDNA) library construction. Here we clarify bioinformatic decisions related to strandedness in both alignment of DNA sequencing reads to reference genomes and subsequent determination of transcript abundance.

Astrocytes perform a wide variety of essential functions defining normal operation of the nervous system, and are active contributors to the pathogenesis of neurodegenerative disorders such as Alzheimer among others. Recent data provide compelling evidence that distinct reactive astrocyte states are associated with specific stages of Alzheimer´s disease. The advent of transcriptomics technologies enables rapid progress in the characterisation of such pathological astrocyte states. In this review, we provide an overview of the origin, main functions, molecular and morphological features of astrocytes in physiological as well as pathological conditions related to Alzheimer´s disease. We will also explore the main roles of astrocytes in the pathogenesis of Alzheimer´s disease and summarize main transcriptional changes and altered molecular pathways observed in astrocytes during the course of the disease.

Comparative metatranscriptomics of bacterial and yeast communities of two milk kefir beverages (MKAA1 and MKAA2) obtained by fermentation with two different frozen stocks of the kefir grain CIDCA AGK1, and differing in rheological features and organic acid amounts production was carried out to figure out the relation between differences in physicochemical and rheological properties of kefir and the microbial active profile associated with each product. The dominance of lactic acid bacteria and yeast and a marginal amount of acetic acid bacteria marked the microbiome. The bacterial families Lactobacillaceae and Streptococcaceae account for almost all bacterial gene transcripts, with Lactobacillus helveticus, L. kefiranofaciens, L. gallinarum and Lactococcus lactis being most frequent in the microbiome of MKAA1 beverage, and L. kefiranofaciens, Lc. lactis and Leuconostoc mesenteroides being most prevalent in the MKAA2. Dipodascaceae and Saccharomycetaceae were the leading yeast families represented by Yarrowia lipolytica, Saccharomyces unisporus, and Kluyveromyces marxianus. The MKAA1 and MKAA2 shared >75% KEGG Ortologs (KO) in bacteria and yeast libraries. Considerable decreases in total expressed genes (KEGG Ortolog) assigned to Lactobacillus helveticus and L. gallinarum might be related to the variations in rheological features of the beverages, probably by compromising the interrelations with L. kefiranofaciens which might explain the variations in rheological features of the beverages.

The Jeryak, a hybrid offspring of Gannan yak and Jersey cattle, exhibits evident hybrid advantages over the Gannan yak in terms of production performance and other various aspects. The small non-coding RNAs known as miRNAs exert a significant regulatory influence on gene expression post-transcriptionally. However, the regulatory mechanism of miRNA associated with muscle development in Jeryak remains elusive. To elucidate the regulatory role of miRNAs in orches-trating the skeletal muscle development in Jeryak, we selected dorsal longissimus muscle tissues from Gannan yak and Jeryak for transcriptome sequencing analysis. A total of 230 (DE)miRNAs and 1819 differentially expressed mRNAs (DEMs) were identified in the longest dorsal muscle of Gannan yak and Jeryak. The functional enrichment analysis revealed a significantly enrichment of target genes from differentially expressed (DE)miRNAs and DEMs in signaling pathways associ-ated with muscle growth, such as the Ras signaling pathway and the MAPK signaling pathway. The network of interactions between miRNA and mRNA suggested that some differentially ex-pressed miRNAs, including miR-2478-z, miR-339-x, novel-m0036-3p and novel-m0037-3p, played a pivotal role in facilitating muscle development. These findings help us to deepen our under-standing of the hybrid dominance of Jeryaks and provide a theoretical basis for further research on the regulatory mechanisms of miRNAs associated with Jeryak muscle growth and development.

The formation and development of potato tissues and organs is a complex process regulated by a variety of genes and environmental factors. However, the regulatory mechanisms underlying the growth and development are still unclear. In this study, we used autotetraploid potato JC14 as experimental subject to analyze the transcriptome of root, stem and leaf at seedling, tuber formation and tuber expansion stages to explore the spatio-temporal expression pattern of genes and genetic development characteristics. The results identified thousands of differentially expressed genes and KEGG pathway enrichment analysis showed that these genes were mainly involved in defense response and carbohydrate metabolism pathways. A total of 12 co-expressed Gene modules were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and 4 modules were screened out with the highest correlation with potato stem developmental traits. Core genes in the network were further investigated and functionally annotated by computing the connectivity of genes within the module. The results unveiled number of hub genes in stems at different developmental stages, including carbohydrate metabolism related genes, the defense response related genes, and transcription factors. These findings provide important leads for further understanding of the molecular regulation and genetic mechanisms of potato tissue development.

Iron is the most abundant micronutrient in plant mitochondria and it has a crucial role in biochemical reactions involving electron transfer. It has been described in Oryza sativa that Mitochondrial Iron Transporter (MIT) is an essential gene and that knockdown mutant rice plants have a decreased amount of iron in mitochondria, strongly suggesting that OsMIT is involved in mitochondrial iron uptake. In Arabidopsis thaliana, two genes encode MIT homologues. In this study, we analyzed different AtMIT1 and AtMIT2 mutant alleles, confirming that individually AtMIT1 nor AtMIT2 genes are essential. When we generated crosses between Atmit1 and Atmit2 alleles we were able to isolate homozygous double mutant plants. Interestingly, homozygous double mutant plants were obtained only when mutant alleles of Atmit2 with the T-DNA insertion in the intron region were used for crossings, and in these cases a correctly spliced AtMIT2 mRNA was generated, although at a low level. Atmit1 Atmit2 double homozygous mutant plants, which were knockout for AtMIT1 and knockdown for AtMIT2, were grown and chacterised in iron sufficient conditions. Pleiotropic developmental defects were observed including abnormal seeds, increased number of cotyledons, slow growth rate, pinoid stems, defects in flower structures and reduced seed set. We observed a possible phenomenon of T-DNA suppression in the next generation of Atmit1 Atmit2 double homozygous mutant plants, correlating with an increased splicing of the AtMIT2 intron containing the T-DNA. Molecular analysis of gene expression markers for mitochondrial and oxidative stress showed that Atmit1 Atmit2 double homozygous mutant plants express a degree of mitochondrial perturbation. A RNA-Seq study was performed and we could identify more than 760 genes differentially expressed in Atmit1 Atmit2, including genes involved in iron transport, coumarin metabolism, and hormones metabolism, transport and signaling. Our data suggest that some of the phenotypes observed in Atmit1 Atmit2 double homozygous mutant plants are mediated by defects in auxin homeostasis.

The mosaic disease in maize is caused by Sugarcane mosaic virus (SCMV), a member of the Potyviridae family. The best strategy to cope with viral infections is the use of disease-resistant maize lines. To better understand the resistance response to SCMV, we analyzed differentially expressed genes among a resistant line (CI-RL1), a susceptible line (B73), and the F1 progeny from a cross between both lines using RNA-Seq data. We also analyzed transcript expression pattern clustering to allocate previously reported resistance candidate genes. GO enrichment analysis of biological processes highlighted a strong regulation in ROS detoxification in both the susceptible and resistant lines. The enrichment of cellular components led to the identification of an integral component of the plasma membrane in the RL line. Transcript expression patterns provide evidence of the importance of host translation in virus response, showing the diverse and complex behavior of eIF4E homologs and the presence of eleven eEF1α factors in maize. In addition, we identified two genes putatively implied in long-distance movement: ZmPiezo and ZmPVIP1. Finally, we propose an ABC transporter to be associated with viral resistance.

The major threats to the sustainable supply of forest tree products are adverse climate, pests and diseases. Climate change, exemplified by increased drought, poses a unique threat to global forest health. This is attributed to the unpredictable behavior of forest pathosystems, which can favor fungal pathogens over the host under persistent drought stress conditions in the future. Currently, the effects of drought on tree resistance against pathogens are hypothetical, thus research is needed to identify these correlations. Norway spruce (Picea abies) is one of the most economically important tree species in Europe, and is considered highly vulnerable to changes in climate. Dedicated experiments to investigate how disturbances will affect the Norway spruce - Heterobasidion sp. pathosystem are important, in order to develop different strategies to limit the spread of H. annosum s.l. under the predicted climate change. Here, we report a transcriptional study to compare Norway spruce gene expressions to evaluate the effects of water availability and the infection of Heterobasidion parviporum. We performed inoculation studies of three-year-old saplings in a greenhouse (purchased from a nursery). Norway spruce saplings were treated in either high (+) or low (-) water groups: high water group received double the water amount than the low water group. RNA was extracted and sequenced. Similarly, we quantified gene expression levels of candidate genes in biotic stress and jasmonic acid (JA) signaling pathways using qRT-PCR, through which we discovered a unique preferential defense response of H. parviporum-infected Norway spruce under drought stress at the molecular level. Disturbances related to water availability, especially low water conditions can have negative effects on the tree host and benefit the infection ability of the pathogens in the host. From our RNA-seq analysis, 114 differentially expressed gene regions were identified between high (+) and low (-) water groups under pathogen attack. None of these gene pathways were identified to be differentially expressed from both non-treated and mock-control treatments between high (+) and low (-) water groups. Finally, only four genes were found to be associated with drought in all treatments.

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from liver, muscle, hypothalamus, pituitary and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.

Removal of the proliferation component of gene expression by PCNA adjustment has been addressed in numerous survival prediction studies for breast cancer and all cancers in the TCGA. These studies indicate that widespread co-regulation of proliferation upwardly biases survival prediction when gene selection is performed on a genome-wide basis. In addition, removal of the correlative effects of proliferation does not reduce the random bias associated with survival prediction using random gene selection. Since most cancers become addicted to DNA repair as a result of forced cellular replication, increased oxidation, and repair deficiencies from oncogenic loss or genetic polymorphisms, we pursued an investigation to remove the proliferation component of expression in DNA repair genes to determine survival prediction. This translational hypothesis-driven focus on DNA repair genes is directly amenable to finding new sets of DNA repair genes that could potentially be studied for inhibition therapy. Overall survival (OS) prediction was evaluated in 18 cancers by using normalized RNA-Seq data for 126 DNA repair genes with expression available in TCGA. Transformations for normality and adjustments for age at diagnosis, stage, and PCNA metagene expression were performed for all DNA repair genes. We also analyzed genomic event rates (GER) for somatic mutations, deletions, and amplification in driver genes and DNA repair genes. After performing empirical p-value testing with use of randomly selected gene sets, it was observed that OS could be predicted significantly by sets of DNA repair genes for 61% (11/18) of the cancers. Interestingly, PARP1 was not a significant predictor of survival for any of the 11 cancers. Results from cluster analysis of GERs indicates that the most opportunistic cancers for inhibition therapy may be AML, colorectal, and renal papillary, because of potentially less confounding due to lower GERs for mutations, deletions, and amplifications in DNA repair genes. However, the most opportunistic cancer for inhibition therapy is likely to be AML, since it showed the lowest GERs for mutations, deletions, and amplifications in DNA repair genes. In conclusion, our hypothesis-driven focus to target DNA repair gene expression adjusted for the PCNA metagene as a means of predicting OS in various cancers resulted in statistically significant sets of genes.

Peanut (Arachis hypogaea) is a crop plant with high economic value, but the epigenetic regulation of its growth and development has only rarely been studied. The peanut histone deacetylase 1 gene (AhHDA1) has been isolated and is known to be ABA- and drought-responsive. In this paper, we investigate the role of AhHDA1 in more detail, focussing on the effect of altered AhHDA1 expression in hairy roots at both the phenotypic and transcriptional levels. Agrobacterium rhizogenes-mediated transformation of A. hypogaea hairy roots was used to analyse how overexpression or RNA interference of AhHDA1 affects this tissue. In both types of transgenic hairy root, RNA sequencing was adopted to identify genes that were differentially expressed, and these genes were assigned to specific metabolic pathways. AhHDA1-overexpressing hairy roots were growth-retarded after 20 d in vitro cultivation, and superoxide anions and hydrogen peroxide accumulated to a greater extent than in control or RNAi groups. Overexpression of AhHDA1 is likely to accelerate flux through various secondary synthetic metabolic pathways in hairy roots, as well as reduce photosynthesis and oxidative phosphorylation. Genes encoding the critical enzymes caffeoyl-CoA O-methyltransferase (Araip.XGB85) and caffeic acid 3-O-methyltransferase (Araip.Z3XZX) in the phenylpropanoid biosynthesis pathway, chalcone synthase (Araip.B8TJ0) and polyketide reductase (Araip.MKZ27) in the flavonoid biosynthesis pathway, and hydroxyisoflavanone synthase (Araip.0P3RJ) and isoflavone 2'-hydroxylase (Araip.S5EJ7) in the isoflavonoid biosynthesis pathway were significantly upregulated by AhHDA1 overexpression, while their expression in AhHDA1-RNAi and control hairy roots remained at a lower level or was unchanged. Our results suggest that alteration of secondary metabolism activities is related to overexpression of AhHDA1, which is mainly reflected in phenylpropanoid, flavonoid and flavonoid biosynthesis. Future studies will focus on the function of AhHDA1 interacting proteins and their action on cell growth and stress responses.

Recently, selection in pigs has been focused on improving the lean meat content in carcasses; this focus has been most evident in breeds constituting a paternal component in breeding. Such sire-breeds are used to improve the meat quantity of cross-breed pig lines. However, even in one breed, a significant variation in the meatiness level can be observed. In the present study, the comprehensive analysis of genes and microRNA expression profiles in porcine muscle tissue was applied to identify the genetic background of meat content. The comparison was performed between whole gene expression and miRNA profiles of muscle tissue collected from two sire-line pig breeds (Piertain, Hampshire). The RNA-seq approach allowed the identification of 627 and 416 differentially expressed genes (DEGs) between pig groups differing in terms of loin weight between Pietrain and Hampshire breeds, respectively. The comparison of miRNA profiles showed differential expression of 57 microRNAs for Hampshire and 34 miRNAs for Pietrain pigs. Next, 43 genes and 18 miRNAs were selected as differentially expressed in both breeds and potentially related to muscle development. According to Gene Ontology analysis, identified DEGs and microRNAs were involved in the regulation of the cell cycle, fatty acid biosynthesis and regulation of the actin cytoskeleton. The most deregulated pathways dependent on muscle mass were the Hippo signalling pathway connected with the TGF-beta signalling pathway and controlling organ size via the regulation of ubiquitin-mediated proteolysis, cell proliferation and apoptosis. The identified target genes were also involved in pathways such as the FoxO signalling pathway, signalling pathways regulating pluripotency of stem cells and the PI3K-Akt signalling pathway. The obtained results indicate molecular mechanisms controlling porcine muscle growth and development. Identified genes (SOX2, SIRT1, KLF4, PAX6 and genes belonging to the transforming growth factor beta superfamily) could be considered candidate genes for determining muscle mass in pigs.

We developed and applied a method for characterizing bacterial promoters genome-wide by in vitro transcription coupled to transcriptome sequencing specific for native 5’-ends of transcripts. This method called ROSE (Run-Off transcription/RNA-SEquencing), only requires chromosomal DNA, ribonucleotides, RNA polymerase (RNAP) core enzyme, and a specific sigma factor, recognizing the corresponding promoters, which have to be analyzed. ROSE was performed on E. coli K-12 MG1655 genomic DNA using E. coli RNAP holoenzyme (including σ70) and yielded 3,226 transcription start sites, 2,167 of which were also identified in in vivo studies, and 598 were new. Many new promoters not yet identified by in vivo experiments might be repressed under the tested conditions. Complementary in vivo experiments with E. coli K-12 strain BW25113 and isogenic transcription factor gene knockout mutants of fis, fur, and hns were used to test this hypothesis. Comparative transcriptome analysis demonstrated that ROSE could identify bona fide promoters that were apparently repressed in vivo. In this sense, ROSE is well-suited as a bottom-up approach for characterizing transcriptional networks in bacteria and ideally complementary to top-down in vivo transcriptome studies.

Expression analysis of small noncoding RNA (sRNA), including microRNA, piwi-interacting RNA, small rRNA-derived RNA, and tRNA-derived small RNA, is a novel and quickly developing field. Despite a range of proposed approaches, selecting and adapting a particular pipeline for transcriptomic analysis of sRNA remains a challenge. This paper focuses on the identification of the optimal pipeline configurations for each step of human sRNA analysis, including reads trimming, filtering, mapping, transcript abundance quantification and differential expression analysis. Based on our study, we suggest the following parameters for analysis of human sRNA in relation to categorical analyses with two groups of biosamples: (1) trimming with the lower length bound = 15 and the upper length bound = \(Read\ length - 40\% Adapter\ length\); (2) mapping on a reference genome with bowtie aligner with one mismatch allowed (-v 1 parameter); (3) filtering by mean threshold > 5; and (4) analyzing differential expression with DESeq2 with adjusted p-value < 0.05 or limma with p-value < 0.05 if there is very little signal and few transcripts.

Natural killer cells are associated with the pathogenesis of ulcerative colitis (UC), but their precise contributions remain unclear. The present study sought to investigate the diagnostic value of activated NK-associated genes (ANAGs) in UC. Bulk RNA-seq and scRNA-seq datasets were obtained from the Gene Expression Omnibus (GEO) and Single Cell Portal (SCP) databases. In the bulk RNA-seq, 92 differentially expressed genes (DEGs) were screened out by the “Batch correction” and “Robust rank aggregation” (RRA) methods. The immune infiltration landscape was estimated by single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT, which revealed a higher abundance of activated NK cells in noninflamed UC tissues. 54 DEGs correlated with activated NK cells were identified as ANAGs. Protein-protein interaction (PPI) analysis and least absolute shrinkage and selection operator (LASSO) regression were utilized to screen out 4 key ANAGs (SELP, TIMP1, MMP7, and ABCG2) and establish an activated NK-associated gene score (ANAG score). The ANAG score demonstrated excellent diagnostic value and was validated in three external datasets. The expression of the 4 key ANAGs was validated in UC patients and healthy controls (HC) samples. Through scRNA-seq data analysis, higher expression levels of SELP, TIMP1, MMP7, and ABCG2 were observed in post-capillary venules, inflammatory fibroblasts, enterocytes, and immature enterocytes. The cell scores based on the ANAGs showed enrichment in endothelial cells and fibroblasts. In conclusion, we established and validated an ANAG score with the ability to precisely diagnose UC. The 4 key ANAGs have the potential to serve as therapeutic targets in UC.

It is established that BCG vaccination results in the development of both a specific immune response to mycobacterial infections and a nonspecific (heterologous) immune response, designated as trained immunity (TRIM), to other pathogens. We hypothesized that local BCG immunization may induce early immune response of bone marrow and spleen innate immunity cells. The early transcriptomic response of various populations of innate immune cells, including monocytes, neutrophils, and natural killer (NK) cells, to BCG vaccination was examined. To this end, C57Bl/6J mice were subcutaneously immunized with 10⁶ CFU of BCG. Three days following BCG administration, the three cell populations were collected from the control and BCG-vaccinated groups by FACS. All cell populations obtained were utilized for the preparation and sequencing of RNA-seq libraries. The analysis of FACS data revealed an increase in the proportion of splenic NK cells and monocytes 3 days post-vaccination. Transcriptomic analysis revealed deregulation of genes associated with regulation of immune response (according to Gene Ontology terms) in NK cells, monocytes, and unsorted bone marrow. Two NK cell-specific immune ligands (Tnfsf14 and S100a8) and two bone marrow-specific immune receptors (C5ar1 and Csf2rb) were identified among differentially expressed genes. No alterations were identified in neutrophils in either the percentage or at the transcriptomic level. Thus, in this study we demonstrated that BCG vaccination provides the early increase in the proportion of murine bone marrow and spleen immune cell populations, as well as the transcriptomic alterations in monocytes, NK cells, and non-sorted bone marrow. This early innate immune response may be beneficial for enhancing TRIM.

Gestational Diabetes Mellitus (GDMs) is worldly wide pregnancy complication. Gestational Diabetes can significantly impact the fetus development. However, the effects of high glucose on embryological development post-fertilization are yet to be researched. Danio rerio embryos are great model to study the embryonic development. In this study, the effects on embryological (morphological and genetic) development were examined in the presence of a high glucose environment that mimic the developing fetus in pregnant women with GMD. Fertilized Zebrafish embryos were treated with normal media and high glucose for 5 days from 3 hours post-fertilization(hpf) to 96hpf respectively as control and experimental groups. Morphological changes are recorded with Microscope images. Hatch rate and heart rate are compared between groups in set up time points. RNA-Seq is performed to examine the gene changes in experimental group. Glucose delayed the zebrafish embryo development by slowing the hatch rate about 24 hours. The brain, heart, and tail started showing morphology smaller in the glucose group compared to the control group at 24 hpf. Heart rate is faster in the glucose group compared to the control group on day 2 and 3 with a statistically significant difference. Among zebrafish whole genome, the significantly changed genes are 556 genes up-regulated and 1118 genes down-regulated respectively in high glucose group. The metabolic and Wnt pathways are altered under high glucose condition. These conditions contribute to significant physiological differences that may provide insight on functionality post-embryological development.

Background: We investigated the tumor immune response in gastric cancer patients receiving third line nivolumab monotherapy to identify immune-related biomarkers for better patient se-lection. Methods: Nineteen patients (10 males, median age 67 years) who received nivolumab as a third or later line therapy were enrolled. We analyzed the tumor immune response in durable clinical benefit (DCB) and non-DCB patients. Pre-treatment and early-on-treatment tumor transcriptomes were examined, and gene expression profiles, immunograms, and T cell receptor (TCR) repertoire were analyzed. Results: DCB was observed in 15.8% of patients, with comparable secondary endpoints (ORR, OS, PFS) to previous trials. Individual immunograms showed no significant changes before and ear-ly-on-treatment, nor consistent alterations among DCB cases. The intratumoral immune response was suppressed by previous treatments in most third or later line nivolumab recipients. TCR repertoire analysis revealed newly emerged clonotypes in on-treatment tumors, but clonal re-placement did not impact efficacy. High T cells/Tregs ratios and a low UV radiation response gene signature were linked to DCB and treatment response. Conclusions: This study emphasizes the tumor immune response's importance in nivolumab ef-ficacy for gastric cancer. High T cells/Tregs ratios and specific gene expression signatures show promise as potential biomarkers for treatment response. Larger cohort validation is crucial to optimize immune checkpoint inhibitors in gastric cancer treatment.

The achievement of reproductive competence in male mammals is dependent on the testis. Goat testis’ development and spermatogenesis involve physiological events with high complexity. In the current work, 6 testes were respectively collected from immature, sexually mature and physically mature Qianbei Ma goats (1, 6 and 12 months old, respectively). RNA-Seq was carried out to reveal changes in testis mRNA expression levels in Qianbei Ma goats at various developmental stages, and gene expression profiling at different ages was established. Totally 18 libraries were established for screening genes and pathways associated with testis development and spermatogenic processes. Totally 9,724 upregulated and 4,153 downregulated genes were identified between immature (I) and sexually mature (S) testes; 7 upregulated and 3 downregulated genes were detected between sexually mature (S) and physically mature (P) testes, and approximately 4% of genes were alternately spliced between the I and S groups. Selected genes were verified by qRT-PCR, in agreement with sequencing data indicating their reliability. Those genes have critical functions in various developmental stages of goat testicular development and spermatogenesis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out to evaluate differentially expressed genes (DEGs). GO analysis suggested DEGs were involved in “reproduction process”, “channel activity” and “cell periphery part” between I and S, and in “ion transport process”, “channel activity” and “transporter complex part” between S and P. KEGG analysis indicated that pathways such as “glycerolipid metabolism”, “steroid hormone biosynthesis” and “MAPK signaling pathway” may be involved in testis development and spermatogenesis. Genes including IGF1, TGFB1, TGFBR1 and EGFR may regulate the development of the testis from immature to sexually mature, which may be key candidate genes for the development of goat testis. These findings provide novel insights into goat testicular development and spermatogenesis.

Rose (Rosa spp.), especially R. hybrida, is one of the most popular ornamental plants in the world and the third largest cut flower crop in Taiwan. Rose mosaic disease (RMD), showing mosaic, line patterns and ringspots on leaves, is a common rose disease caused by complex infection of various viruses. Due to pests and diseases, the rose planting area in Taiwan has been decreasing since 2008; however, no rose virus disease has been reported in the past five decades. In the spring of 2020, rose samples showing RMD-like symptoms were observed at an organic farm in Chiayi, central Taiwan. The virome in the farm was analyzed by RNA-seq. Rose genomic sequences were filtered from the obtained reads. The remaining reads were de novo assembled to generate 294 contigs, 50 of which were annotated as viral sequences corresponding to 10 viruses. Through reverse transcription-polymerase chain reaction validation, a total of seven viruses were detected, including six known rose viruses, namely apple mosaic virus, prunus necrotic ringspot virus, rose partitivirus, apple stem grooving virus, rose spring dwarf-associated virus and rose cryptic virus 1, and a novel ilarvirus. After completing the whole genome sequencing and sequence analysis, the unknown ilarvirus was demonstrated as a new species, tentatively named rose ilarvirus 2. This is the first report of the rose viruses in Taiwan.

The objective of this study was to identify key transcription factors involved in lipid metabolism and immune response related to the differentially expressed genes (DEG) from the liver samples of 35 pig model for metabolic diseases fed diets containing either 1.5 or 3.0% soybean oil (SOY1.5 or SOY3.0). A total of 281 DEG between SOY1.5 and SOY3.0 diets (log2fold-change ≥ 1 or ≤ −1; FDR-corrected p-value &lt; 0.1) were identified, in which 129 were down-regulated and 152 were up-regulated in SOY1.5 group. The functional annotation analysis detected transcription factors linked to lipid homeostasis and immune response, such as RXRA, EGFR, and SREBP2 precursor. These findings demonstrated that key transcription factors related to lipid metabolism could be modulated by dietary inclusion of soybean oil. It could contribute to nutrigenomics research field that aims to elucidate dietary interventions in animal and human health, as well as to drive the food technology and science.

Background: Jerusalem artichoke (Helianthus tuberosus L.) is tolerant to salinity stress and has high economic value. The salt tolerance mechanisms of Jerusalem artichoke are still unclear. Especially in the early stage of Jerusalem artichoke exposure to salt stress, the plant physiology, biochemistry and gene transcription are likely to undergo large changes. Elucidating these changes may be of great significance to understanding the salt tolerance mechanisms of it. Results: We obtained high-quality transcriptome from leaves and roots of Jerusalem artichoke exposed to salinity (300 mM NaCl) for 0 h, 6 h, 12 h, 24 h and 48 h, with 150,129 unigenes and 9023 DEGs (Differentially Expressed Genes). The RNA-seq data were clustered into time-dependent groups (nine clusters each in leaves and roots); gene functions were distributed evenly among the groups convergence. KEGG enrichment analysis showed the genes related to plant hormone signal transduction were enriched in almost all treatment comparisons. Under salt stress, genes belongs to PYL (abscisic acid receptor PYR / PYL family), PP2C (Type 2C protein phosphatases), GH3 (Gretchen Hagen3), ETR (ethylene receptor), EIN2/3 (ethylene-insensitive protein 2/3), JAZ (Genes such as jasmonate ZIM-domain gene) and MYC2 (Transcription factor MYC2) had extremely similar expression patterns. The results of qPCR of 12 randomly selected genes confirmed the accuracy of RNA-seq. Conclusions: Under the impact of high salinity (300mM) environment, Jerusalem artichoke in the seedling stage was difficult to survive for a long time, and the phenotype was severe in the short term. Based on the expression of genes on the time scale, we found that the distribution of gene functions in time is relatively even. Upregulation of the phytohormone signal transduction had a crucial role in the response of Jerusalem artichoke seedlings to salt stress, the genes of abscisic acid, auxin, ethylene, and jasmonic acid had the most obvious change pattern.

Single cell RNA-seq (scRNA-seq) profiles conceal temporal and spatial tissue developmental information. De novo reconstruction of single cell temporal trajectory has been fairly addressed, but reverse engineering single cell 3D spatial tissue localization is hitherto landmark based, and de novo spatial reconstruction is a compelling computational open problem. Here we show that a new algorithm - named D-CE - for coalescent embedding of single cell transcriptomic networks can address this open problem. We rely merely on the spatial information encoded in the expression patterns of developmental signal transcription factor (DST) genes, and we find that D-CE of cell-cell association DST-transcriptomic networks reliably reconstructs the Geo-seq or single cell samples’ 3D spatial tissue distribution. Comparison to the novoSpaRC and CSOmap (recent and only available de novo 3D spatial reconstruction methods) on 16 datasets and 681 reconstructions, reveals a significantly distinctive superior performance of D-CE.

Background: DSCAM-AS1 is a cancer-related long noncoding RNA with higher expression levels in Luminal A, B and HER2-positive Breast Cancer (BC), where its expression is strongly dependent on Estrogen Receptor Alpha (ERα). Methods: To decipher its function, DSCAM-AS1 expression was measured by qRT-PCR in tissue samples from 93 BC patients in addition to a meta-analysis of 30 gene expression datasets, together with the evaluation of its association with clinical data. By computational analyses of our RNA-Seq in MCF-7 cells, we investigated the DSCAM-AS1 knock-down effects at both gene and isoform levels. Results: We confirmed DSCAM-AS1 overexpression in high grade Luminal A, B and HER2+ BCs and found a significant correlation with disease relapse. 908 genes were regulated by DSCAM-AS1-silencing, primarily involved in cell cycle and inflammatory response. Noteworthy, the analysis of alternative splicing and isoform regulation revealed 2,085 splicing events regulated by DSCAM-AS1, enriched in differential polyadenylation sites and 3’UTR shortening events. Finally, the DSCAM-AS1-interacting splicing factor hnRNPL was predicted as the most enriched RBP for exon skipping and 3’UTR events. Conclusion: The relevance of DSCAM-AS1 overexpression in BC is confirmed by clinical data and further enhanced by its possible involvement in the regulation of RNA processing, which is emerging as one of the most important dysfunctions in cancer.

Powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt) is one of many severe diseases that threaten bread wheat (Triticum aestivum L.) yield and quality worldwide. The discovery and deployment of powdery mildew resistance genes (Pm) can prevent this disease epidemic in wheat. In a previous study, we transferred the powdery mildew resistance gene Pm57 from Aegilops searsii into common wheat and cytogenetically mapped the gene in a chromosome region with the fraction length (FL) 0.75-0.87, which represents 12% of 2Ss#1 segment on the long arm of chromosome 2Ss#1. In this study, we performed RNA-Seq on three infected and mock-infected wheat-Ae. searsii 2Ss#1 introgression lines with Bgt-isolates inoculation at 0, 12, 24, and 48 hours after inoculation. Then we designed 79 molecular markers based on transcriptome sequences and physically mapped them to Ae. searsii chromosome 2Ss#1- in seven intervals. We used these markers to identify 46 wheat-Ae. searsii 2Ss#1 recombinants induced by ph1b, a deletion mutant of pairing homoelogous (Ph) genes. Analysis of the 46 ph1b-induced 2Ss#1L recombinants with different Bgt-responses using 28 2Ss#1L-specific molecular markers in the interval FL0.72-0.87 where Pm57 is located, and the flanking intervals, we physically mapped Pm57 gene on the long arm of 2Ss#1 in a 5.13 Mb genomic region, which was flanked by markers X67593 (773.72 Mb) and X62492 (778.85 Mb). By comparative synteny analysis of the corresponding region on chromosome 2B in Chinese spring (T. aestivum L.) with other model species we identified ten genes that are putative plant defense-related (R) genes which includes six coiled-coil nucleotide-binding site-leucine-rich repeat (CNL), three nucleotide-binding site-leucine-rich repeat (NL) and a leucine-rich receptor-like repeat (RLP) encoding proteins. This study will lay a foundation for further cloning of Pm57, and benefit the understanding of interactions between resistance genes of wheat and powdery mildew pathogens.

Epigenetic mechanisms, including histone post-translational modifications (PTMs), play a critical role in regulating pain perception and the pathophysiology of burn injury. However, the epigenetic regulation and molecular mechanisms underlying burn injury-induced pain remain insufficiently explored. Spinal dynorphinergic (Pdyn) neurons contribute to heat hyperalgesia induced by severe scalding-type burn injury through p-S10H3-dependent signaling. Beyond p-S10H3, burn injury may impact various other histone H3 PTMs. Double immunofluorescent staining and histone H3 protein analyses demonstrated significant hypermethylation at H3K4me1 and H3K4me3 sites and hyperphosphorylation at S10H3 within the spinal cord. By analyzing Pdyn neurons in the spinal dorsal horn, we found evidence of chromatin activation with a significant elevation in p-S10H3 immunoreactivity. We used RNA-seq analysis to compare the effects of burn injury and forma-lin-induced inflammatory pain on spinal cord transcriptomic profiles. We identified 98 DEGs for burn injury and 86 DEGs for formalin-induced inflammatory pain. A limited number of shared differentially expressed genes (DEGs) suggest distinct central pain processing mechanisms between burn injury and formalin models. KEGG pathway analysis supported this divergence, with burn injury activating Wnt signaling. This study enhances our understanding of burn injury mechanisms and uncovers converging and diverging pathways in pain models with different origins.

Antisense oligonucleotides (ASOs) have emerged as a promising therapeutic approach for a wide variety of diseases, offering the potential to modulate gene expression with a high degree of specificity. These small, synthetic nucleic acid-like molecules provide unique advantages over traditional pharmacological treatments, including the ability to target previously "undruggable" genes. Despite this promise, several biological barriers severely limit their clinical efficacy Upon administration, ASOs primarily enter cells through endocytosis, becoming trapped inside membrane-bound vesicles known as endosomes. Studies estimate that only 1-2% of ASOs successfully escape endosomal compartments to reach the cytosol, and ultimately the nucleus, where they bind target mRNA and exert their therapeutic effect. Endosomal entrapment and inefficient nuclear localization are therefore critical bottlenecks in the therapeutic application of ASOs. This review explores the current understanding of endosomal escape, and nuclear transport of ASOs. Furthermore, we examine strategies aimed at overcoming these challenges, including the use of endosomal escape agents, peptide-ASO conjugates, non-viral delivery vehicles, and nuclear localization signals. By improving both endosomal escape and nuclear localization, significant advances in ASO-based therapeutics can be realized, ultimately enhancing their potency and expanding their clinical utility

Through the integration of genomic information, transcriptome sequencing data, and bioinformatics methods, we conducted a comprehensive identification of the ALDH gene family in melon. We explored the impact of this gene family on melon growth, development, and their expression patterns in various tissues and under different stress conditions. Our study discovered a total of 17 ALDH genes spread across chromosomes 1, 2, 3, 4, 5, 7, 8, 11, and 12 in the melon genome. Through a phylogenetic analysis, these genes were classified into 10 distinct subfamilies. Notably, genes within the same subfamily exhibited consistent gene structures and conserved motifs. Our study discovered a pair of fragment repeats within the melon ALDH gene. Furthermore, there was a noticeable collinearity relationship between the melon's ALDH gene and that of Arabidopsis (12 times), and rice (3 times). Transcriptome data reanalysis revealed that some ALDH genes consistently expressed highly across all tissues and developmental stages, while others were tissue or stage-specific. We analysed the ALDH gene's expression patterns under six stress types, namely salt, chilling, waterlogged, powdery mildew, fusarium wilt, and gummy stem blight. The results showed differential expression of MELO3C024328 and MELO3C004430 under all stress conditions, signifying their crucial roles in melon growth and stress response. RT-qPCR analysis further corroborated these findings. This study paves the way for future genetic improvements in melon molecular breeding.

Cotton seeds as the main by-product of cotton are not only an important raw material for edible oil and feed but also a source of biofuel. The quality of cotton seeds directly affects cotton planting and is closely related to the yield and fiber quality. However, the molecular mechanism governing cotton seed size remains largely unexplored. This study investigates the regulatory mechanisms of cotton seed size by focusing on two cotton genotypes, N10 and N12, which exhibit notable phenotypic variations across multiple environments. Developing seeds were sampled at various stages (5, 20, 30, and 35 DAP) and subjected to RNA-seq. The temporal pattern clustering and WGCNA on differentially expressed genes identified 413 candidate genes, including these related to sugar metabolism that were significantly enriched in transcriptional regulation. A genetic transformation experiment indicated that the overexpression of the GhUXS5 gene encoding UDP-glucuronate decarboxylase 5 significantly increased seed size, suggesting an important role of GhUXS5 in regulating cotton seed size. This discovery provides crucial insights into the molecular mechanisms controlling cotton seed size, helping to unravel the complex regulatory network and offering new strategies and targets for cotton breeding to enhance the economic value of cotton seeds and overall cotton yield.

Streptococcus agalactiae is a highly contagious Gram- positive bacterium that causes mastitis, has a high infectivity for mammary epithelial cells, and is challenging to treat. However, the molecular interactions between it and mammary epithelial cells remains poorly understood. This study analyzed differential gene expression in mammary epithelial cells with varying levels of S. agalactiae infection using UID- Dual Seq and bioinformatics tools. This study identified 211 differentially expressed mRNAs (DEmRNAs) and 452 differentially expressed lncRNAs (DElncRNAs) in host cells, primarily enriched in anti-inflammatory responses, immune responses, and cancer-related processes. Additionally, 854 pathogen differentially expressed mRNAs (pDEmRNAs) were identified, mainly enriched in protein metabolism, gene expression, and biosynthesis processes. Mammary epithelial cells activate pathways, such as the ERK1/2 pathway to produce, to produce reactive oxygen species (ROS) to eliminate bacteria. However, the bacteria disrupt the host's innate immune mechanisms by interfering with the alternative splicing processes of mammary epithelial cells. Specifically, the bacterial genes of tsf, prfB, and infC can interfere with lncRNAs targeting RUNX1 and BCL2L11 in mammary epithelial cells, affecting the alternative splicing of target genes and altering normal molecular regulation.

With the rapid development of gene therapy technology in recent years, its abuse as a method of sports doping in athletic has become a concern. However, there is still room for improvement in gene doping testing methods, and a robust animal model needs to be developed. Therefore, the purposes of this study were to establish a model of gene doping using recombinant adeno-associated virus vector-9 including the human erythropoietin gene (rAAV9-hEPO), and to establish a relevant testing method. First, it was attempted to establish the model using rAAV9-hEPO on mice. The results showed a significant increase in erythrocyte volume accompanied by an increase in spleen weight, confirming the validity of the model. Next, we attempted to detect proof of gene doping by targeting DNA and RNA. Direct proof of gene doping was detected using a TaqMan-qPCR assay with certain primers/probes. In addition, some indirect proof was identified in RNAs through the combination of a TB Green qPCR assay with RNA sequencing. Taken together, these results could provide the foundation for an effective test for gene doping in human athletes in the future.

The addition of 50% vitreous humor to the media bathing lens explants results in the induction of both fiber cell differentiation and a significant immune/inflammatory response. Although Fgfr loss blocks differentiation in lens epithelial explants, this loss is partially rescued by deleting Pten. To investigate the functions of the Fgfrs and Pten during lens fiber cell differentiation, we utilized a lens epithelial explant system and conducted RNA sequencing on vitreous humor-exposed explants lacking Fgfrs, or Pten or both Fgfrs and Pten. We found that Fgfr loss impairs both vitreous-induced differentiation and inflammation while the additional loss of Pten restores these responses. Furthermore, transcriptomic analysis suggested that PDGFR-signaling in FGFR-deficient explants is required to mediate the rescue of vitreous-induced fiber differentiation in explants lacking both Fgfrs and Pten. The blockage of β-crystallin induction in explants lacking both Fgfrs and Pten in the presence of a PDGFR-inhibitor supports this hypothesis. Our findings demonstrate that a wide array of genes associated with fiber cell differentiation are downstream of FGFR-signaling and that the vitreous-induced immune responses also depend on FGFR-signaling. Our data also demonstrates that many of the vitreous-induced gene expression changes Fgfr-deficient explants are rescued in explants lacking both Fgfrs and Pten.

Trichothecene mycotoxin T-2 toxin is an important environmental pollutant and poses a global threat to human and animal health. T-2 toxin could induce nephrotoxicity, however, the precise molecular mechanism remains unclear. In this study, the mice were intraperitoneally administrated with a single dose of 2 mg/kg T-2 toxin. The kidney function and ultrastructural observation were also assessed and found T-2 toxin exposure could cause kidney damage. Transmission electron microscopy showed a significant swelling and vacuolization of the mitochondria of renal cell. A total of 1122, 58, and 391 differential expression genes (DEGs) that were predominantly expressed in kidney tissues at 1 d, 3 d, and 7 d after T-2 toxin exposure, respectively. Early transcriptomic changes at 1 d after exposure were down-regulation of DEGs involved in the cell cycle, p53 signaling pathway, and cellular senescence, while up-regulation of DEGs referred to the ribosome pathway. The temporal variations in gene expression pattern in T-2 toxin exposure in mice kidney were presented and cellular metabolism was disturbed at the exposure recovery period with 7 d. In conclusion, this study, for the first time, provided a comprehensive comparative transcriptomic analysis of gene regulation in T-2 toxin exposure-induced nephrotoxicity at different temporal periods and explored the nephrotoxicity mechanism of T-2 toxin at the mRNA level.

The shoot apical meristem (SAM) gives rise to the aerial structure of plants by producing lateral organs and secondary meristems. The SAM is responsible for plant developmental patterns, thus determining plant morphology and, consequently, many agronomic traits such as the number and size of fruits and flowers and kernel yield. Our current understanding of SAM morphology and regulation is based on studies conducted mainly on some angiosperms, including economically important crops such as maize (Zea mays) and rice (Oryza sativa), and the model species Arabidopsis (Arabidopsis thaliana). However, studies in other plant species from the gymnosperm class are scant, making difficult comparative analyses that help us understand SAM regulation in diverse plant species. This limitation prevents deciphering the mechanisms by which evolution gave rise to the multiple plant structures within the plant kingdom, and determined the conserved mechanisms involved in SAM maintenance and operation. This review aims to integrate and analyze the current knowledge of SAM evolution by combining the morphological and molecular information recently reported from the plant kingdom.

Because rice is native to tropical and subtropical regions, it is generally more sensitive to cold stress compared to other cereals. Low temperature stress has become one of the most important conditions that affects the growth, development, and yield of rice. Plant NAC (NAM, ATAF1/2, and CUC) transcription factors are known to play significant regulatory roles in the stress response. In our research, we found that OsNAC050 (LOC_Os03g60080） is mainly expressed in leaves, and low temperature can further up-regulate its expression. OsNAC050 mutants created using CRISPR-Cas9 gene editing technology showed significantly enhanced tolerance to low temperature treatment. Detection of enzyme activities related to the redox pathway also showed that the mutants had stronger viability under low temperature stress. Comparative transcriptome analysis showed that photosynthesis and soluble sugar metabolism were significantly affected in the osnac050 mutant lines, suggesting that OsNAC050 may participate in the above molecular pathways in the response to low temperature stress. The results will enhance our understanding of molecular mechanisms underlying the responses to cold stress in rice and can provide new strategies for engineering cold-tolerance in high-yielding rice varieties.

The emergence of drug resistant pathogenic bacteria is increasingly challenging conventional antibiotics. Plant derived flavonoids are considered as potential alternatives to antibiotics due to their antimicrobial properties. However, the mechanisms by which flavonoids modulate pathogenic microorganisms’ growth are not fully understood. In our previous studies we found that Rutin, a kind of flavonoids showed inhibition against Klebsiella pneumoniae strains. In order to better understand its inhibitory mechanism on strain growth, we used both Rutin and Luteolin as treatments to incubate K. pneumoniae ATCC700603 strain. After incubation for 4 hours, the homogenous and differential effects of two different flavonoids on the growth of K. pneumoniae were evaluated. The integrated metabolomic and transcriptomic analysis was performed. The expression levels of 5,483 genes and the contents of 882 metabolites were measured. The differentially expressed genes (DEGs) and metabolites (DEMs) were screened and analyzed. The correlation between DEGs and DEMs were also studied. Our results showed that Rutin could inhibit the strain growth by changing metabolic pathways and ABC transporters pathways. Our study also revealed FU841_RS17580 and FU841_RS19145 as functional genes that played vital roles in the strain growth.

Systems genetics is key for integrating a large number of variants associated with diseases. Vitamin K (VK) is one of the scarcely studied conditions in lieu of ascertaining either the differentially expressed genes (DEGs) or variants in an individual subpopulation of diseased phenotypes associated with VK, viz. myocardial infarction, renal failure, prostate cancer, thrombosis, thrombocytopenia, coagulation related diseases to name a few. In this work, we have screened characteristic DEGs common to three VK-related diseases, viz. myocardial infarction, renal failure and prostate cancer and asked whether or not any DEGs in addition to pathogenic variants are common to these conditions. We attempt to bridge the gap in finding characteristic biomarkers and discuss the role of long noncoding RNAs (lncRNAs) in the biogenesis of VK deficiencies.

Single-cell RNA sequencing data facilitates investigation of cell heterogeneity and subpopulations as well as differentially abundant states however modern single-cell RNA sequencing datasets are growing in size and complexity requiring advances in the bioinformatic methods that analyze them. Many methods exist for each step of analysis including read alignment, normalization, quality control, batch effect correction, imputation and dimensionality reduction. With so many options to choose from at each step of the analysis, benchmarking and a synthesis of the literature on the methods available is necessary to inform biological researchers on the most optimal workflow for their data. Here, recent key methods of analysis are highlighted with a focus on methods that facilitate identification of cell subpopulations and differentially abundant cell states. With a constantly expanding toolset for each step in single-cell RNA sequencing dataset analysis, biological researchers should stay informed to utilize the most applicable methods for their own analyses.

Microglia activation causes neuroinflammation, which is a hallmark of neurodegenerative disorders, brain injury, and aging. Ladostigil, a bifunctional reagent with antioxidant and anti-inflammatory properties, reduced microglial activation and enhanced brain functioning in elderly rats. In this study, we studied SH-SY5Y, a human neuroblastoma cell line, and tested viability in the presence of hydrogen peroxide and Sin1 (3-morpholinosydnonimine), which generates reactive oxygen and nitrogen species (ROS/RNS). Both stressors caused significant apoptosis and necrotic cell death that was attenuated by ladostigil. Our results from RNA-seq experiments show that long non-coding RNAs (lncRNAs) account for 30% of all transcripts in SH-SY5Y cells treated with Sin1 for 24 hours. Altogether, we identify 94 differently expressed lncRNAs in the presence of Sin1, including MALAT1, a highly expressed lncRNA with anti-inflammatory and anti-apoptotic functions. Additional activities of Sin-1 upregulated lncRNAs include redox homeostasis (e.g., MIAT, GABPB1-AS1), energy metabolism (HAND2-AS1), and neurodegeneration (e.g., MIAT, GABPB1-AS1, NEAT1). Four lncRNAs implicated as enhancers were significantly upregulated in cells exposed to Sin1 and ladostigil. Finally, we show that H2O2 and Sin1 increased the expression of DJ-1, a redox sensor and modulator of Nrf2 (nuclear factor erythroid 2–related factor 2). Nrf2 (NFE2L2 gene) is a major transcription factor regulating antioxidant genes. In the presence of ladostigil, DJ-1 expression is restored to its baseline. The mechanisms governing SH-SY5Y cell survival and homeostasis are highlighted by the beneficial role of ladostigil in the crosstalk involving Nrf2, antioxidant transcription factor DJ-1, and lncRNAs. Stress-dependent induction of lncRNAs represents an underappreciated regulatory level that contributes to cellular homeostasis and the capacity of SH-SY5Y to cope with oxidative stress.

Intratumoral heterogeneity associates with more aggressive disease progression and worse patient outcomes. Understanding the reasons enabling the emergence of such heterogeneity remains incomplete, which restricts our ability to manage it from a therapeutic perspective. Technological advancements such as high-throughput molecular imaging, single-cell omics and spatial transcriptomics now allow recording the patterns of spatiotemporal heterogeneity in a longitudinal manner, thus offering insights into the multi-scale dynamics of its evolution. Here, we review latest technological trends and biological insights from molecular diagnostics as well as spatial transcriptomics, both of which have witnessed a burgeoning growth in recent past in terms of mapping heterogeneity within tumor cell types as well as stromal constitution. We also discuss ongoing challenges, indicating possible ways to integrate insights across these methods to have a systems-level spatiotemporal map of heterogeneity in each tumor, and a more systematic investigation of implications of heterogeneity for the patient outcomes.

As a major infectious disease in chickens, Newcastle disease causes considerable economic losses in the poultry industry, especially in developing countries where there is limited access to effective vaccination. Therefore, enhancing resistance to the virus in commercial chickens through breeding is a promising way to promote poultry production. In this study, we investigated gene expression changes at 2 and 6 dpi after infection at day21 with a lentogenic Newcastle disease virus in a commercial egg-laying chicken hybrid using RNA sequencing analysis. By comparing NDV challenged and nonchallenged groups, 526 differentially expressed genes (DEGs) (FDR &lt; 0.05) were identified at 2 dpi, and only 36 at 6 dpi. For the DEGs at 2 dpi, IPA analysis predicted inhibition of multiple signaling pathways in response to NDV that regulate immune cell development and activity, neurogenesis and angiogenesis. Upregulation of Interferon Induced Protein with Tetratricopeptide Repeats 5 (IFIT5) in response to NDV was consistent between the current and most previous studies. Sprouty RTK Signaling Antagonist 1 (SPRY1), a DEG in the current study is located in a significant QTL associated with virus load at 6 dpi in the same population. These identified pathways and DEGs provide potential targets to further study breeding strategy to enhance NDV resistance in chickens.

An outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) occurred in China towards the end of 2019, and has spread rapidly ever since. Previous studies showed that some virus could affect the reproductive system and cause long-term complications. Recent studies exploring the source of SARS-CoV-2 using genomic sequencing have revealed that SARS-CoV-2 enters the host cells via the angiotensin-converting enzyme II (ACE2), the receptor that recognizes SARS-CoV. To investigate the expression of ACE2 and to explore the potential risk of infection in the reproductive system, we performed a thorough bioinformatic analysis on data from public databases involving RNA expression, protein expression, and single-cell RNA expression studies. The analyzed data showed high levels of ACE2 mRNA and protein expression in the testis and spermatids and equal levels of ACE2 expression in the uterus and lung. Comprehensive single-cell analysis identified ACE2 expression in the lung, testis, spermatids, and uterus. In conclusion, this study revealed the potential risk associated with the SARS-CoV-2 infection in the reproductive system and predicted that long-term complications might have a significant impact on the prevention and management of COVID-19, the disease caused upon infection with SARS-CoV-2.

PARK7, also known as DJ-1, plays a critical role in protecting cells by functioning as a sensitive oxidation sensor and modulator of antioxidants. DJ-1 acts to maintain mitochondrial function and regulate transcription in response to different stressors. In this study, we showed that cell lines vary based on their antioxidation potential at basal conditions. The transcriptome of HEK293 cells was tested following knockdown (KD) of DJ-1 using siRNAs which reduced the DJ-1 transcripts to only 12% of the original level. We compared the expression levels of 14k protein-coding transcripts, and 4.2k non-coding RNAs relative to control cells treated with non-specific siRNAs. Among the coding genes, ~200 upregulated differentially expressed genes (DEGs) signified a coordinated antiviral innate immune response. Most genes were associated with regulation of type 1 interferons (IFN) and induction of inflammatory cytokines. About a quarter of these genes were also induced in cells treated by non-specific siRNAs that were used as a negative control. Beyond the antiviral response, 114 genes were specific to KD of DJ-1 with enrichment in RNA metabolism and mitochondrial functions. A smaller set of downregulated genes (58 genes) were associated with dysregulation in membrane structure, cell viability, and mitophagy. We propose that KD of DJ-1 diminish its protective potency against oxidative stress, rendering the cells labile and responsive to dsRNA signal by activation of a large number of genes, many of which drive apoptosis, cell death, and inflammatory signatures. The KD of DJ-1 highlights its crucial role in regulating genes associated with antiviral responses, RNA metabolism, and mitochondrial functions, apparently through alteration in STAT activity and downstream signaling. Given that DJ-1 is highly expressed in metastatic cancers, targeting DJ-1 could be a promising therapeutic strategy where manipulation of DJ-1 levels may reduce cancer cell viability and enhance the efficacy of cancer treatments.

The interaction between native and invasive plant species plays a crucial role in understanding invasion mechanisms and developing control strategies. This study investigated the differences in photosynthetic and fluorescence parameters, gene expression profiles, and rhizosphere soil microbial composition between the native species Ficus tikoua Bur. and the invasive species Alternanthera philoxeroides. For the first time, it was demonstrated that that F. tikoua Bur. may enhance photoinhibition ability and the expression of genes involved in carbon metabolism, protein processing in the endoplasmic reticulum, and plant-pathogen interactions to maintain growth. Additionally, it recruits beneficial microbes, e.g. Actinomarina (Bacteria), Nitrosotenuis (Archaea) and Laccaria (fungi), to suppress pathogens, thereby preserving soil DMSP and carbon cycling pathways to resist invasion by A. philoxeroides. These findings provide novel insights into the mechanisms underlying plant invasion and suggest potential strategies for invasive plant management. Further research is needed to investigate this relationship and determine the mechanisms underlying the potential resistance of F. tikoua Bur. to plant invasion.

Intramuscular fat content (IFC) and meat color are vital indicators affecting pork quality. A significant positive correlation between IFC and redness meat color (a* value) indicates that these two traits are likely to be regulated by shared molecular pathways. To identify candidate genes, hub genes, and signaling pathways that regulate these two traits, we measured the IFC and a* value in 147 hybrid pigs and selected individuls with extreme phenotypes for transcriptome analysis. The results revealed 485 and 394 overlapping differentially expressed genes (DEGs) using DESeq2, limma, and edgeR packages affecting IFC and a* value, respectively. Weighted gene co-expression network analysis (WGCNA) identified four modules significantly correlated with IFC and a* value. Moreover, we integrated functional enrichment analysis results based on DEGs, GSEA, and WGCNA conditions to identify candidate genes, and identified 47 and 53 candidate genes affecting IFC and a* value, respectively. The protein-protein interaction (PPI) network analysis of candidate genes showed that 6 and 14 hub genes affecting IFC and a* value, respectively. These genes mainly participated in various pathways related to lipid metabolism and redox reactions. Notably, four crucial hub genes (MYC, SOX9, CEBPB, and PPAGRC1A) were shared for these two traits. After functional annotation of these four hub genes, we hypothesized that the SOX9/CEBPB/PPARGC1A axis could co-regulate lipid metabolism and myoglobin redox response. Further research on these hub genes, especially the SOX9/CEBPB/PPARGC1A axis, will help to understand the molecular mechanism of the co-regulation of IFC and a* value, which will provide the theoretical basis for improving pork quality.

Injury to the central nervous system (CNS), in most vertebrate animals, results in permanent damage and lack of function, due to their limited regenerative capacities. In contrast, echinoderms can fully regenerate their radial nerve cord (RNC) following transection, with little or no scarring. Investigators have associated the regenerative capacity of some organisms with the stress response and inflammation produced by the injury. Here we explore the gene activation profile of the stressed holothurian CNS. To do this, we performed RNA sequencing on isolated RNC explants submitted to the stress of transection and enzyme dissection and compared them to explants kept in culture for 3 days following dissection. We describe stress-associated genes, including members of heat-shock families, ubiquitin-related pathways, transposons, and apoptosis that were differentially expressed. Surprisingly, the stress response does not induce apoptosis in this system. Other genes associated with stress in other animal models, such as hero proteins and those associated with the integrated stress response, were not found to be differentially expressed either. Our results provide a new viewpoint on the stress response in the nervous system of an organism with an amazing regenerative capacity. This is the first step to deciphering the molecular processes that allow echinoderms to undergo fully functional CNS regeneration while also providing a comparative view for students of the stress response in other organisms.

Numerous proteomic and transcriptomic studies have been carried out to better understand the current multi-variant SARS-CoV-2 virus mechanisms of action and effects. However, they are mostly centered on mRNAs and proteins. The effect of the virus on human post-transcriptional regulatory agents such as microRNAs (miRNAs) involved in the regulation of 60% of human gene activity remains poorly explored. Similar to what we have previously done with other viruses such as Ebola and HIV, in this study we investigated the miRNA profile of lung epithelial cells following infection with SARS-CoV-2. At the 24 and 72 hours post-infection, SARS-CoV-2 did not drastically alter the miRNome. About 90% of the miRNAs remained non-differentially expressed. The results revealed that miR-1246, miR-1290 and miR-4728-5p were the most upregulated over time. miR-196b-5p and miR-196a-5p were the most downregulated at 24 h while at 72 h, miR-3924, miR-30e-5p and miR-145-3p showed the highest level of downregulation. In the top significantly enriched KEGG pathways of genes targeted by differentially expressed miRNAs we found, among others, MAPK, RAS, P13K-Akt and renin secretion signaling pathways. By RT-qPCR, we also showed that SARS-CoV-2 may regulate several predicted host mRNA targets involved in the entry of the virus into host cells (ACE2, TMPRSS2, ADAM17 and FURIN), in renin–angiotensin system (RAS) (Renin, Angiotensinogen, ACE), innate immune response (IL-6, IFN1β, CXCL10, SOCS4) and fundamental cellular processes (AKT, NOTCH, WNT). Finally, we demonstrated by dual luciferase assay a direct interaction between miR-1246 and ACE-2 mRNA. This study highlights the modulatory role of miRNAs in the pathogenesis of SARS-CoV-2.

Inducible genetically defined mouse models of cancer uniquely facilitate the investigation of early events in cancer progression, however there are valid concerns about the ability of such models to faithfully recapitulate human disease. We developed an inducible mouse model of progressive lung adenocarcinoma (LuAd) that combines sporadic activation of oncogenic KRas^G12D with modest overexpression of c-MYC (KM model). Histological examination revealed a highly reproducible transition from adenoma to locally invasive adenocarcinoma within 6 weeks of oncogene activation. Laser-capture microdissection coupled with RNA-SEQ was employed to determine transcriptional changes associated with tumour progression. Upregulated genes were triaged for relevance to human LuAd using datasets from Oncomine and cBioportal. Selected genes were validated by RNAi screening in human lung cancer cell lines and examined for association with lung cancer patient overall survival using KMplot.com. Depletion of progression-associated genes resulted in pronounced viability and/or cell migration defects in human lung cancer cells. Progression-associated genes moreover exhibited strong associations with overall survival, specifically in human lung adenocarcinoma, but not in squamous cell carcinoma. The KM mouse model faithfully recapitulates key molecular events in human lung cancer and is a useful tool for mechanistic interrogation of LuAd progression.

Plants grown under fluctuating light impact plant developments compared with those grown under non-fluctuating light conditions. However, our knowledge on the underlying regulatory mechanisms is still quite limited, particularly from the transcriptional perspective. In order to investigate the influence of different light intensity distributions on tomato plant development, we designed three fluctuating light intensity distributions with the non-fluctuating light intensity as control and compared the transcriptional differences after five weeks of treatment. We found plant height and aerial/root weight were significantly reduced under all fluctuating light treatments. Transcriptome analysis revealed that the number of up and down regulated genes had a distinct distribution pattern between different treatments and control. The largest difference between the numbers of down and up regulated genes was found between treatment 1 and 3, reaching to a total of 416 genes. The number and type of the top 20 enriched pathways differed between treatments and control. The largest number of genes enriched was involved in the biosynthesis of secondary metabolites. These results provide insights into the transcriptional regulations of tomato under different light intensity distributions.

B chromosomes (B) are supernumerary elements found in many taxonomic groups. Most B chromosomes are rich in heterochromatin and composed of abundant repetitive sequences, especially transposable elements (TEs). Bs origin is generally linked to the A chromosome complement (A). The first report of a B chromosome in African cichlids was on Astatotilapia latifasciata, which can harbor 0, 1 or 2 B chromosomes. Classical cytogenetics studies found high TE content on the species B chromosome. In this study, we aim to understand TE composition and expression on A. latifasciata genome and its relation to the B chromosome. We use bioinformatics analysis to explore TEs genome organization and also their composition on the B chromosome. Bioinformatics findings were validated by fluorescent in situ hybridization (FISH) and real-time PCR (qPCR). A. latifasciata has a TE content similar to other cichlid fishes and several expanded elements on its B chromosome. With RNA sequencing data (RNA-seq) we showed that all major TE classes are transcribed in brain, muscle and male/female gonads. The evaluation of TE expression between B- and B+ individuals showed that few elements have differential expression among groups and expanded B elements were not highly transcribed. Putative silencing mechanisms may the acting on the B chromosome of A. latifasciata to prevent adverse consequences of repeat transcription and mobilization in the genome.

Solanesol is a noncyclic terpene alcohol composed of nine isoprene units and it mainly accumulates in solanaceous plants, especially tobacco (Nicotiana tabacum L.). Here, RNA-seq analyses of tobacco leaves, stems, and roots were used to identify solanesol biosynthesis genes. Six 1-deoxy-d-xylulose 5-phosphate synthase, two 1-deoxy-d-xylulose 5-phosphate reductoisomerase, two 2-C-methyl-d-erythritol 4-phosphate cytidylyltransferase, four 4-diphosphocytidyl-2-C-methyl-d-erythritol kinase, two 2-C-methyl-d-erythritol 2,4-cyclodiphosphate synthase, four 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase, two 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase, six isopentenyl diphosphate isomerase, and two solanesyl diphosphate synthase (SPS) genes were identified to be involved in solanesol biosynthesis. Furthermore, the two N. tabacum SPS (NtSPS1 and NtSPS2), which had two conserved aspartate-rich DDxxD domains, were highly homologous with SPS enzymes from other solanaceous plant species. In addition, the solanesol contents of three organs, and leaves from four growing stages, corresponded with the distribution of chlorophyll. Our findings provide a comprehensive evaluation of the correlation between the expression of different biosynthetic genes and the accumulation of solanesol in tobacco.

Apple red flesh coloration is a result of the performance of biochemical pathway involved in the biosynthesis of anthocyanins and anthocyanidins. Based on apple genome analysis, high number of regulatory genes, mainly transcription factors such as MYB, which are components of regulatory complex MYB-bHLH-WD40, and several structural genes (PAL, 4CL, CHS, CHI, F3H, DFR, ANS, UFGT) involved in anthocyanin biosynthesis have been identified so far. In this study, we investigated previously unrecognized genes involved in the constitution of red-flesh apple phenotype. These genes could be concerned as molecular markers for early selection of new apple cultivars. Based on comparative transcriptome analysis of apples with different fruit flesh coloration, we successfully identified and characterized ten genes from plant hormone transduction pathway of: auxin (GH3); cytokinins (B-ARR); gibberellins (DELLA); abscisic acid (SnRK2 and ABF); brassinosteroids (BRI1, BZR1 and TCH4); jasmonic acid (MYC2); and salicylic acid (NPR1). In presented studies an analysis of their expression profiles was performed in immature and ripe fruits of red-fleshed cultivars. Generally, we have uncovered, the set of the novel genes, mediating regulation of the abscisic acid, salicylic acid, cytokinins and jasmonic acid signaling and preliminarily described their role in anthocyanin biosynthesis, accumulation and degradation. Their expression was not characterized in apple fruits so far. Presented study, allowed to underline the relationship of genes from hormone signal transduction pathway with UFGT gene, directly responsible for anthocyanin color transformation as well as confirmed their crucial role of plant hormone regulation and anthocyanin accumulation during apple fruit ripening.

Legume-rhizobial symbiosis leads to to the formation of root nodules, where rhizobia reside and develop into bacteroids to reduce nitrogen into ammonium for plant growth, which leaves an opening question as how plant immunity is attenuated in nodules in the presence of large number of foreign bacteria. In Medicago truncatula, mutation in NAD1 (Nodules with Activated Defense 1) only produces necrotic nodules with overactivated immunity, indicating NAD1 is an indispensable component required for suppressing immunity in nodules. In this study, a dual RNA-seq transcriptomic technology was performed to extensively analyze gene expressions in nad1-1 mutant nodules. We indenifited 89 DEGs in symbiotic nitrogen fixation and 89 DEGs in immunity in Medicago truncatula at 10 dpi. Simultaneously, we indenifited 27 DEGs in fix and nif genes in Sinorhizobium meliloti. Then we identified 56 DEGs related to stress, including ROS, NO, and NCR, most of which were upregulated in Sinorhizobium meliloti. Our analyses of nitrogen fixation-defective plant nad1-1 mutants with overactivated defense indicate that host use plant immunity to control massive bacterial colonization during early and late symbiotic stages. Our findings provide insight into the function of NAD1 in improvement or inhibition of plant immune response to maintain symbiosis during rhizobial endosymbiosis.

Gene expression data is key for the functional annotation of single nucleotide polymorphisms (SNPs) identified in genome-wide association studies (GWAS). Expression and splicing quantitative trait loci (e/sQTLs) in normal colon tissue, such as those from the University of Barcelona and University of Virginia RNA sequencing project (BarcUVa-Seq) and the Genotype-Tissue Expression project (GTEx), are required to gain biological insight of colon-related diseases risk loci. Moreover, transcriptome-wide association studies (TWAS) rely on reference gene expression imputation panels in the tissue of interest to nominate susceptibility genes. Also, it is of high interest to study the relationships between genes in a network framework. For facilitating these analyses, we have updated and expanded the scope of the Colon Transcriptome Explorer (CoTrEx) to the version 2.0. This web-based resource provides exhaustive visualization and analysis of transcriptome-wide gene expression profiles of normal colon tissue from BarcUVa-Seq and GTEx. In addition to the integration of new datasets, CoTrEx 2.0 provides additional e/sQTLs sets, as well as gene expression prediction models and regulatory and co-expression networks. It is freely available at https://barcuvaseq.org/cotrex/. Overall, it is of high interest for researchers aiming to investigate the genetic susceptibility to colon-related complex traits and diseases.

Hungateiclostridium thermocellum ATCC 27405 is a promising bacterium with a robust ability to degrade lignocellulosic biomass complexes, including crystalline cellulose components, through a multienzyme cellulosomal system. In contrast, it exhibits poor growth on simple monosaccharides such as fructose and glucose. This phenomenon raises many important questions concerning its glycolytic pathways and sugar transport systems. Until now, the detailed mechanisms of H. thermocellum adaptation to growth on monosaccharides have been poorly explored. In this study, adaptive laboratory evolution was applied to train the bacterium on monosaccharides, and genome resequencing was used to detect the genes that had mutated during adaptation. RNA-seq data of the 1st-generation culture growing on either fructose or glucose revealed that several glycolytic genes in the EMP pathway were expressed at lower levels in these cells than in cellobiose-grown cells. After 8 generations of culture on fructose and glucose, the evolved H. thermocellum strains grew faster and yielded greater biomass than the nonadapted strains. Genomic screening also revealed several mutation events in the genomes of the evolved strains, especially in genes responsible for sugar transport and central carbon metabolism. Consequently, these genes could be applied as targets for further metabolic engineering to improve this bacterium for bioindustrial usage.

COVID-19 was an inflammatory pneumonia caused by a respiratory infection with a coronavirus named SARS-CoV-2. Severely ill COVID-19 patients could die of acute pulmonary and systemic inflammation. Thus, there was a theory that long COVID-19 shared many similarities with systemic autoimmune diseases. There was also a theory that the SARS-CoV-2 virus could infect many cell- and tissue-types distributed throughout multiple organ systems of the body. I aim to review the above theories using unbiased high-throughput datasets cataloguing gene expression. Studies reviewed used the technology of next generation sequencing, either as bulk tissue RNA sequencing, or as single-cell RNAs sequencing; ingenious methods for signal deconvolution were then used to identify individual cell types. Datasets reviewed suggested that severe COVID-19 induced expression of genes associated with pro-inflammatory signaling, tissue development or remodeling, scarring, or cell death. The datasets also provided evidence for several inflammatory syndromes associated with COVID-19, such as neuronal COVID, acute respiratory disease syndrome, vascular inflammation, or multisystem inflammatory syndrome. Indeed, inflammatory signal during severe COVID-19 was reported in many cell types, for example in structural cells, in infiltrating immune cells, in the vascular system, or in circulating cells in the blood. There was also frequently a change in proportions of immune cell-types in the blood, epithelial cells in the lungs, or among infiltrating immune cells.

Carbamate insecticides have been used for over four decades to control brown planthopper, Nilaparvata lugens, but resistance has been reported in many countries, including Korea. The test results on resistance to fenobucarb showed that the LC₅₀ values were 3.08 for the susceptible strain, 10.61 for the 2015 population, and 73.98 for the 2019 population. Compared to the susceptible strain, the 2015 and 2019 populations exhibited resistance levels 3.4 and 24.2 times higher, respectively. To elucidate the reason for the varying levels of resistance to fenobucarb in these populations, mutations in the ace1 gene, the target site of carbamate, were investigated, but no previously reported mutations were confirmed. Through RNAseq analysis focusing on the expression of detoxification enzyme genes as an alternative resistance mechanism, it was found that the carboxylesterase gene Nl-EST1 was overexpressed 2.4 times in the 2015 population and 4.7 times in the 2019 population compared to the susceptible strain. This indicates a strong correlation between the level of resistance development in each population and the expression level of Nl-EST1. Previously, Nl-EST1 was reported in an organophosphorus insecticide-resistant strain of Sri Lanka 2000. Thus, Nl-EST1 is crucial for developing resistance to organophosphorus and carbamate insecticides. Resistance-related genes such as Nl-EST1 could serve as expression markers for resistance diagnosis, and it can apply to integrated resistance management of N. lugens.

Fusarium wilt (FW) caused by Fusarium oxysporum f. sp. Cucumerinum (Foc) of cucumber is a destructive soil-borne disease in cucumber (Cucumis sativus.L). There is limited knowledge on molecular mechanisms to FW resistance-mediated defense responses in cucumber. In this study, metabolome and transcriptome profiling were carried out in two FW, resistant (NR) and susceptible (NS), near isogenic lines (NILs) before and after Foc inoculation. The NILs showed consistent and stable FW resistances in multiple greenhouse and laboratory screening tests. Widely targeted metabolomic analysis identified differentially accumulated metabolites (DMAs) with significantly more accumulation in NR in response to Foc infection including many phenolic acid and flavonoid compounds from the flavonoid biosynthesis pathway. Transcriptome analysis identified differentially expressed genes (DEGs) between the NILs upon Foc inoculation including genes for secondary metabolite biosynthesis in or transcription factor genes regulating the flavonoid biosynthesis pathway. Joint analysis of the metabolomic and transcriptomic data identified DAMs and DEGs closely associated with biosynthesis of phenolic acid and flavonoid DAMs. The association of these compounds with NR-conferred FW resistance was exemplified by in vivo assays in which two phenolic acid compounds, bis (2-ethylhexyl) phthalate and diisooctyl phthalate, and the flavonoid compound gallocatechin 3-O-gallate were shown to have significant inhibitory effect on Foc growth. The antifungal effect of three compounds seems to be a novel finding. We conclude that phenolic acids and flavonoids play an important role in NR mediated FW resistance in cucumber.

The Pashmina and Barbari are two famous goat breeds found in the wide areas of the Indo-Pak region. Pashmina is famous for its long hair-fiber (Cashmere) production while Barbari is not-selected for this trait. So, the mRNA expression profiling in the skin samples of both breeds would be an attractive and judicious approach for detecting putative genes involved in this valued trait. Here, we performed differential gene expression analysis on publicly available RNA-Seq data from both breeds. Out of 44,617,994 filtered reads of Pashmina and 55,995,999 of Barbari which are 76.48% and 73.69% mapped to the ARS1 reference transcriptome assembly respectively. A pairwise comparison of both breeds resulted in 47,159 normalized expressed transcripts while 8,414 transcripts are differentially expressed above the significant threshold. Among these, 4,788 are upregulated in Pashmina while 3,626 transcripts are upregulated in Barbari. Fifty-nine transcripts harbor 57 genes including 32 LOC genes and 24 are annotated genes which were selected on the basis of TMM counts > 500. Genes with ectopic expressions other than uncharacterized and LOC symbol genes are Keratins (KRT) and Keratin Associated Proteins (KRTAPs), CystatinA&6, TCHH, SPRR4, PPIA, SLC25A4, S100A11, DMKN, LOR, ANXA2, PRR9 and SFN. All of these genes are likely to be involved in keratinocyte differentiation, sulfur matrix proteins, dermal papilla cells, hair follicles proliferation, hair curvature, wool fiber diameter, hair transition, hair shaft differentiation and its keratinization. These differentially expressed reported genes are critically valuable for enhancing the quality and quantity of the pashmina fiber and overall breed improvement. This study will also provide important information on hair follicle differentiation for further enrichment analyses and introducing this valued trait to other goat breeds as well.

Corynebacterium pseudotuberculosis is a Gram-positive bacterium that causes caseous lymphadenitis, a disease that predominantly affects sheep, goat, cattle, buffalo, and horses, but has also been recognized in other animals. This bacterium generates a severe economic impact on countries producing meat. Gene expression studies using RNA-seq is one of the most commonly used techniques to perform transcriptional experiments. Computational analysis on such data through reverse-engineering algorithms leads to a better understanding of the genome-wide complexity of gene interactomes, enabling the identification of genes having the most significant functions inferred by the activated stress response pathways. In this study, we identified the influential or causal genes from four RNA-seq data-sets from different stress conditions (high iron, low iron, acid, osmosis, and PH) in C. pseudotuberculosis, using a consensus-based network inference algorithm called miRsig and identified the causal genes in the network using the miRinfluence tool, which is based on the influence diffusion model. We found that over 50\% of the genes identified as influential have some essential cellular functions in the genomes. In the strains analyzed, most of the causal genes have crucial roles or participate in processes associated with response to extracellular stresses, pathogenicity, membrane components, and essential genes. This research brings new insight into the understanding of virulence and infection by C. pseudotuberculosis.

Single-cell RNA sequencing (scRNA-seq) technique has enabled detailed analysis of gene expression at the single cell level, enhancing the understanding of subtle mechanisms that underly pathologies and drug resistance. To derive such biological meaning from sequencing data in oncology, some critical processing must be performed, including identification of the tumor cells by markers and algorithms that infer copy number variations (CNVs). We compared the performance of sciCNV, InferCNV, CopyKAT and SCEVAN tools that identify tumour cells by inferring CNVs from scRNA-seq data. Sequencing data from Pancreatic Ductal Adenocarcinoma (PDAC) patients, adjacent and healthy tissues were analysed, and the predicted tumor cells were compared to those identified by well-assessed PDAC markers. Results from InferCNV, CopyKAT and SCEVAN overlapped by less than 30% with InferCNV showing the highest sensitivity (0.72) and SCEVAN the highest specificity (0.75). We show that the predictions are highly dependent on the sample and the software used, and that they return so many false positives hence are of little use in verifying or filtering predictions made via tumor biomarkers. We highlight how critical this processing can be, warn against the blind use of these software and point out the great need for more reliable algorithms.

One of the ways to regulate the sensitivity of human cells to the influenza virus is to knock out genes of the innate immune response. Promising targets for the knockout are genes of the interferon-inducible transmembrane protein (IFITM) family, in particular the IFITM3 gene, whose product limits the entry of a virus into the cell by blocking the fusion of the viral and endosomal membranes. In this study, by means of genome-editing system CRISPR/Cas9, monoclonal cell lines with an IFITM3 knockout were obtained based on WI-38 VA13 cells (human origin). It was found that such cell lines are more sensitive to infection by influenza A viruses of various subtypes. Nonetheless, this feature is not accompanied by an increased titer of newly formed viral particles in a culture medium.

Polo-Like Kinase-1 (PLK1) is a key regulator of the G2/M cell cycle and has been linked to treatment resistance in prostate cancer (PCa). This study aimed to explore the correlation between PLK1 expression, determined through RNA-seq analysis, and clinical outcomes in PCa patients. We analyzed tumor samples from 452 patients diagnosed with prostate adenocarcinoma between 2014 and 2020, sourced from nine institutions affiliated with the Oncology Research Information Exchange Network (ORIEN). PLK1 expression levels were determined using RNA-seq data and stratified into high (≥75th percentile) vs. low (<75th percentile) expression. Univariate and multivariate Cox proportional hazard models compared overall survival between patients with low vs. high PLK1 expression, adjusting for age, race, PSA, and stage. Of the 452 tumors, 241 (53.3%) had high PLK1 and 211 (46.7%) had low PLK1 expression. Median follow-up was 2.74 years for low PLK1 and 2.71 years for high PLK1. No significant difference in overall survival was observed between low and high PLK1 groups in univariate (p=0.86; HR 1.09) or multivariate (p=1.0; HR 0.99) analyses. This study found no association between RNA-seq-based PLK1 expression and prognosis in PCa. PLK1 expression is not a prognostic marker for prostate cancer.

Aedes aegypti transmits major arboviruses of public health importance, including dengue, chikungunya, Zika, and yellow fever. The use of insecticides represents the cornerstone of vector control, however insecticide resistance in Ae. aegypti has become widespread. Understanding the molecular basis of insecticide resistance in this species is crucial to design effective resistance management strategies. Here, we applied Illumina RNA-Seq to study the gene expression patterns associated with resistance to three widely used insecticides (malathion, alpha-cypermethrin and lambda-cyhalothrin) in Ae. aegypti populations from 2 sites (Manatí and Isabela) in Puerto Rico (PR). Cytochrome P450s were the most over-expressed detoxification genes across all resistant phenotypes. Some detoxification genes (CYP6Z7, CYP28A5, CYP9J2, CYP6Z6, CYP6BB2, CYP6M9, and two CYP9F2 orthologs) were commonly over-expressed in mosquitoes that survived exposure to all three insecticides (independent of geographical origin) while others including CYP6BY1 (malathion), GSTD1 (alpha-cypermethrin), CYP4H29 and GSTE6 (lambda-cyhalothrin) were uniquely overexpressed in mosquitoes that survived exposure to specific insecticides. The gene ontology (GO) terms associated with monooxygenase, iron binding, and passive transmembrane transporter activities were significantly enriched in 4 out of 6 resistant vs susceptible comparisons while serine protease activity was elevated in all insecticide resistant groups relative to the susceptible strain. Interestingly, cuticular-related protein genes (chinase and chitin) were predominantly downregulated, which was also confirmed in the functional enrichment analysis. This RNA-Seq analysis presents a detailed picture of the candidate detoxification genes and other pathways that are potentially associated with pyrethroid and organophosphate resistance in Ae. aegypti populations from PR. These results could inform development of novel molecular tools for detection of resistance-associated gene expression in this important arbovirus vector and guide the design and implementation of resistance management strategies.

Intratumor heterogeneity (ITH) is associated with therapeutic resistance and poor prognosis in cancer patients, and attributed to genetic, epigenetic, and microenvironmental factors. We developed a new computational platform, GATHER, for geostatistical modeling of single cell RNA-seq data to synthesize high-resolution and continuous gene expression landscapes of a given tumor sample. Such landscapes allow GATHER to map the enriched regions of pathways of interest in the tumor space and identify genes that have spatial differential expressions at locations representing specific phenotypic contexts using measures based on optimal transport. GATHER provides new applications of spatial entropy measures for quantification and objective characterization of ITH. It includes new tools for insightful visualization of spatial transcriptomic phenomena. We illustrate the capabilities of GATHER using real data from breast cancer tumor to study hallmarks of cancer in the phenotypic contexts defined by cancer associated fibroblasts.

Heat stress affects the granulosa cells (GCs) and ovarian follicular microenvironment, causing poor oocyte developmental competence and fertility. This study aimed to investigate the physical responses and global transcriptomic changes in bovine GCs to acute heat stress (43 ℃ for 2 h) in-vitro and gave essential insights into the general interaction at cell–stress nexus. Heat-stressed GCs exhibited transient proliferation senescence, resumed proliferation at 48 h post-stress. While post-stress immediate culture-media change had a relatively positive effect on proliferation resumption. Increased accumulation of reactive oxygen species and apoptosis was observed in heat stress group. In spite of the upregulation of pro-apoptotic and caspase executioner genes, antioxidants and anti-apoptotic genes were also upregulated in heat-stressed GCs. Progesterone and Estrogen hormones along with steroidogenic genes expression, declined significantly, in spite of the upregulation of genes involved in cholesterol synthesis. Out of 12385 differentially expressed genes (DEGs), 330 significant DEGs (75 upregulated, 225 downregulated) were subjected to KEGG functional pathway annotation, gene ontology enrichment, and STRING network analyses. Based on the manual query of DEGs, pathway and enrichment analyses, a vast interplay observed among all major signaling pathways strongly evidence the repression of cellular transcriptional and proliferation activity, averting the effects of heat stress through remodeling of cellular structural proteins and energetic-homeostasis. This study presents detailed responses of acute heat-stressed GCs at physical, transcriptional, and pathway levels and presents interesting insights into future studies regarding GCs adaptation and their interaction with oocyte and reproductive system at ovarian level.

Post-translational modifications (PTMs) of histone residues shape the landscape of gene expression by modulating the dynamic process of RNAPII transcription. The contribution of particular histone modifications to the definition of distinct RNAPII transcription stages remains poorly characterized in plants. Chromatin Immuno-precipitation combined with next-generation sequencing (ChIP-seq) resolves the genomic distribution of histone modifications. Here, we review histone PTM ChIP-seq data in Arabidopsis thaliana and find support for a Genomic Positioning System (GPS) that guides RNAPII transcription. We review the roles of histone PTM “readers”, “writers” and “erasers”, with a focus on the regulation of gene expression and biological functions in plants. The distinct functions of RNAPII transcription during the plant transcription cycle may in part rely on the characteristic histone PTMs profiles that distinguish transcription stages.

Injuries to the anterior cruciate ligament (ACL) often result in post-traumatic osteoarthritis (PTOA). To better understand the molecular mechanisms behind PTOA development following ACL injury, we profiled ACL injury-induced gene expression changes in knee joints of three mouse strains with varying susceptibility to OA: STR/ort (highly susceptible), C57BL/6 (moderately susceptible) and super-healer MRL/MpJ (not susceptible). Right knee joints of the mice were injured using a non-invasive tibial compression injury model that closely mimics ACL rupture in humans and global gene expression was quantified before and at 1-day, 1-week, and 2-weeks post-injury using RNA-seq. Following injury, STR/ort displayed severe cartilage degeneration while MRL/MpJ had little cartilage damage. Gene expression analysis suggested that prolonged inflammation and elevated catabolic activity in STR/ort injured joints, compared to the other two strains may be responsible for the severe PTOA phenotype observed in this strain. MRL/MpJ had the lowest expression values for several inflammatory cytokines and catabolic enzymes activated in response to ACL injury. Furthermore, we identified several genes highly expressed in MRL/MpJ compared to the other two strains including B4galnt2 and Tpsab1 which may contribute to enhanced healing in the MRL/MpJ. Overall, this study has increased our knowledge of early molecular changes associated with PTOA development.

In the face of rising global demand and unsustainable production methods, cultivated crustacean meat (CCM) is proposed as an alternative means to produce delicious lobster, shrimp and crab products. Cultivated meat requires starting stem cells that may vary in terms of potency, and propensity to proliferate or differentiate into myogenic (muscle-related) tissues. We suggest that regenerating crustacean limbs harbor a range of stem cells suitable for cultivated meat production and present a humane tissue source because they can be accessed non-lethally. To investigate stem cell activity in this model, we have conducted RNA-Seq analysis across six stages of claw regeneration in the crayfish Cherax quadricarinatus (four pre-molt and two post-molt stages), along with histology and real time quantitative PCR (qPCR). Differential expression analysis revealed expression patterns of significantly up- or down-regulated genes in 14 functional categories. Relative expression of select stem cell target genes (including endocrine, myogenic, cell cycle and pluripotency factors) highlighted two growth peaks during pre-molt regeneration and limited stem cell activity after the molt. Histology and qPCR results support the overall finding that pre-molt limb regeneration stages harbor more stem cell activity, likely making them a more prolific cell source for CCM development.

The role of desmoglein-3 (DSG3) in oncogenesis is unclear. This study aimed to uncover molecular mechanisms through comparative transcriptome analysis in oral cancer cells, defining potential key genes and associated biological processes related to DSG3 expression. Four mRNA libraries of oral squamous carcinoma H413 cell lines were sequenced and 599 candidate genes exhibited differential expression between DSG3-overexpressing and matched control lines, with 12 genes highly significantly differentially expressed, including 9 upregulated and 3 downregulated. Genes with known implications in cancer, such as MMP-13, KRT84, OLFM4, GJA1, AMOT, and ADAMTS1, were strongly linked to DSG3 overexpression. Gene ontology analysis indicated that the DSG3-associated candidate gene products participated in crucial cellular processes such as junction assembly, focal adhesion, extracellular matrix formation, intermediate filament organization, and keratinocyte differentiation. Validation of RNA-Seq was performed through qRT-PCR, Western blotting, and immunofluorescence analyses. Furthermore, transmission electron microscopy meticulously examined desmosome morphology, revealing slightly immature desmosome structure in DSG3-overexpressing cells compared to controls. No changes in desmosome frequency and diameter were observed between the two conditions. This study underscores intricate and multifaceted alterations associated with DSG3 in oral squamous carcinoma cells, implying a potential oncogenic role of this gene in biological processes that enable cell communication, motility and survival.

During the infection by SARS-CoV-2, the virus is changing infected host cell into its own factory producing new viral particles. As infection progresses, infected cell undergoes many changes in various pathways. One of the events caused by changes is cytokine storm, which leads to the severe symptoms. In this study, we investigated transcriptomic changes caused by COVID-19 disease using RNA-seq data obtained from COVID-19-positive patients and COVID-19-negative donors. RNA-seq data were collected for the purpose of identification of potential biomarkers associated with a different course of the disease. Here, the first datasets of 96 samples were analyzed to validate the methods. The aim of this publication is to report pilot results. In search of potential biomarkers associated with different disease severity, we performed differential expression analysis of human transcriptome, focusing on COVID-19 positivity and symptom severity. Since we detected plenty of potential biomarkers, we performed KEGG enrichment analysis to get better view of altered pathways. Results show, that affected were pathways related to immune processes and response to infection, also multiple signaling pathways, while most of them were also reported to be influenced by SARS-CoV-2 infection in previous studies.

The generation of different types of defective viral genomes (DVG) is an unavoidable consequence of the error-prone replication of RNA viruses. In recent years, a particular class of DVGs, those containing long deletions or genome rearrangements, has gain interest due to their potential therapeutic and biotechnological applications. Identifying such DVGs in high-throughput sequencing data has become an interesting computational problem. Up to nowadays, several algorithms have been proposed, though all incur in false positives, a problem of practical interest if such DVGs have to be synthetized and tested in the laboratory. Here we develop a novel software, DVGfinder, that wraps the two most commonly used algorithms into a pipeline that predicts DVGs. Using a gradient boosting classifier machine learning algorithm, we evaluate the performance of DVGfinder compared to previous algorithms and found that it outcompetes their precision and sensitivity in simulated datasets. DVGfinder generates user-friendly output files in HTML format that can assist users to identify DVGs based on their associated probability of being true positives.

Precision medicine emphasizes fine-grained diagnostics, taking individual variability into account to enhance treatment effectiveness. Parkinson's Disease (PD) heterogeneity among individuals is a proof that disease subtypes exist, and assigning individuals to subgroups is necessary for a better understanding of disease mechanisms and designing precise treatment approaches. The purpose of this study was to identify PD subtypes using RNA-Seq data in a combined pipeline including unsupervised machine learning, bioinformatics, and network analysis. 210 post mortem brain RNA-Seq samples from PD (n = 115) and Normal Controls (NC, n = 95) were obtained with a systematic data retrieval following PRISMA statements and a fully data-driven clustering pipeline was performed to identify PD subtypes. Bioinformatics and Network analyses were performed to characterize the disease mechanisms of the identified PD subtypes and to identify target genes for drug repurposing. Two PD clusters were identified and 42 DEGs were found (p.adjusted ≤ 0.01). PD clusters had significantly different gene network structures (p &lt; 0.0001) and phenotype-specific disease mechanisms, highlighting the differential involvement of the Wnt/β-catenin pathway regulating adult neurogenesis. NEUROD1 was identified as a key regulator of gene networks and ISX9 and PD98059 were identified as NEUROD1-interacting compounds with disease-modifying potential, reducing the effects of dopaminergic neurodegeneration. This hybrid data analysis approach could enable precision medicine applications by providing insights for the identification and characterization of pathological subtypes. This workflow has proven useful on PD brain RNA-Seq, but its application to other neurodegenerative diseases is encouraged.

In a changing climate, extreme weather events such as heat waves will be more frequent and could affect grain weight and the quality of crops such as wheat, one of the most significant crops in terms of global food security. In this work, we characterized the response of Triticum turgidum spp. durum wheat to a short-term heat-stress (HS) treatment at transcriptomic and physiological levels during early grain filling in glasshouse experiments. We found a significant reduction in grain weight and size from HS treatment. Grain quality was also affected, showing a decrease in starch content in addition to increments in grain protein levels. Moreover, an RNA-seq analysis of durum wheat grains allowed us to identify 1590 differentially expressed genes related to photosynthesis, response to heat, and carbohydrate metabolic process. A gene regulatory network analysis of HS-responsive genes uncovered novel transcription factors (TFs) controlling the expression of genes involved in abiotic stress response and grain quality, such as a member of the DOF family predicted to regulate glycogen and starch biosynthetic processes in response to HS in grains. In summary, our results provide new insights into the extensive transcriptome reprogramming that occurs during short-term HS in durum wheat grains.

Cervical carcinogenesis is the leading cause of cancer related death in women and the role of high-risk human papillomavirus (HR-HPV) as a possible risk factor in the development of this cancer is well recognised. Despite the availability of multi-therapeutic approaches, there is still a major concern regarding the prevention of metastatic dissemination and excessive tissue injuries. Therefore, it is imperative to develop a safer and more efficient treatment modality. Ficus Carica, a natural plant, has shown potential therapeutic properties through its fruit latex when applied to HPV-positive cervical cancer cell lines. However, the mechanisms of action of Ficus Carica (fig) latex are not well understood. This study aims to provide a deeper insight into the biological activities of fig latex on human cervical cancer cell lines expressing high-risk HPV types 16 and 18. The data obtained from this study reveals that fig latex influences the expression of genes in-volved in "Class I MHC-mediated antigen presentation" as well as "Antigen processing: Ubiquitination and Proteasome degradation." These genes play a crucial role in host immune surveillance and the resolution of infection. Findings from this study suggest that fig latex may enhance T cell responses against oncogenic HPV, which could be beneficial for the clearance of early-stage cancer.

Multiple mRNA isoforms of the same gene are produced via alternative splicing, a biological mechanism that regulates protein diversity while maintaining genome size. Alternatively spliced mRNA isoforms of the same gene may sometimes have very similar sequence, but they can have significantly diverse effects on cellular function and regulation. The products of alternative splicing have important and diverse functional roles, such as response to environmental stress, regulation of gene expression, human heritable and plant diseases. The mRNA isoforms of the same gene, such as the apoptosis associated CASP3 gene, can have dramatically different functions. The shorter mRNA isoform product CASP3-S inhibits apoptosis, while the longer CASP3-L mRNA isoform promotes apoptosis. Despite the functional importance of mRNA isoforms, very little has been done to annotate their functions. The recent years have however seen the development of several computational methods aimed at predicting mRNA isoform level biological functions. These methods use a wide array of proteo-genomic data to develop machine learning-based mRNA isoform function prediction tools. In this review, we discuss the computational methods developed for predicting the biological function at the individual mRNA isoform level.

Prebiotic oligosaccharides are widely used as human and animal feed additives for their beneficial effects on the gut microbiota. However, there are limited data to assess the direct effect of such functional foods on the transcriptome of intestinal epithelial cells. The purpose of this study is to describe the differential transcriptomes and cellular pathways of colonic cells directly exposed to galacto-oligosaccharides (GOS) and fructo-oligosaccharides (FOS). We have examined the differential gene expression of polarized Caco-2 cells treated with GOS or FOS and their respective mock-treated cells using mRNA sequencing (RNA-seq). A total of 89 significant differentially expressed genes were identified between GOS and mock-treated groups. For FOS treatment, a reduced number of 12 significant genes were observed to be differentially expressed relative to the control group. KEGG and Gene Ontology functional analysis revealed that genes up-regulated in the presence of GOS were involved in digestion and absorption processes, fatty acids and steroids metabolism, potential antimicrobial proteins, energy-dependent and -independent transmembrane trafficking of solutes and amino acids. Using our data, we have established complementary non-prebiotic modes of action for these frequently used dietary fibers.

Combined awareness about the power and limitations of bioinformatics and molecular biology enables advanced research based on high-throughput data. Despite an increasing demand for scientists with a combined background in both fields, the education in dry lab and wet lab is often separated. This work describes an example of integrated education with focus on genomics and transcriptomics. Participants learn computational and molecular biology methods in the same practical course. Peer-review is applied as a teaching method to foster cooperative learning of students with heterogeneous backgrounds. Evaluation results indicate acceptance and appreciation of this approach.

Transposable elements (TEs) are repetitive elements which make up around 45% of the human genome. A class of TEs known as SINE-VNTR-Alu (SVA) demonstrate the capacity to mobilise throughout the genome, resulting in SVA polymorphisms for presence or absence within the population. Although studies have previously highlighted the involvement of TEs within neurodegenerative diseases, such as Parkinson’s disease and amyotrophic lateral sclerosis (ALS), however the exact mechanism has yet to be identified. In this study we used whole genome sequencing and RNA sequencing data of ALS patients and healthy controls from the New York Genome Center ALS Consortium, to elucidate the influence of reference SVA elements on gene expression genome-wide within central nervous system (CNS) tissues. To investigate this, we applied matrix expression quantitative trait loci analysis and demonstrated that reference SVA insertion polymorphisms can significantly modulate the expression of numerous genes, preferentially in the trans position, and in a tissue-specific manner. We also highlight that SVAs significantly regulate mitochondrial genes as well as genes within the HLA and MAPT loci, previously associated within neurodegenerative disease. In conclusion, this study continues to bring to light the effects of polymorphic SVAs on gene regulation and further highlights the importance of TEs within disease pathology.

Lysine is the first limiting essential amino acid in rice because it is present in the lowest quantity compared to all the other amino acids. Amino acids are the building block of proteins and play an essential role in maintaining the human body’s healthy functioning. Rice is a staple food for large proportion of the global population, thus increasing the lysine content in rice will improve its nutritional value. In this paper, we studied the lysine biosynthesis pathway in rice (Oryza Sativa) to identify the regulators of the lysine reporter gene LYSA (LOC_Os02g24354). Genetically intervening at the regulators has the potential to increase the overall lysine content in rice. We modeled the lysine biosynthesis pathway in rice seedlings under normal and saline (NaCl) stress conditions using Bayesian networks. We estimated the model parameters using experimental data and identified the gene DAPF(LOC_Os12g37960) as a positive regulator of the lysine reporter gene LYSA under both normal and saline stress conditions. Based on this analysis, we conclude that the gene DAPF is a potent candidate for genetic intervention. Upregulating DAPF using methods such as CRISPR-Cas9 has the potential to upregulate the lysine reporter gene LYSA and increase the overall lysine content in rice.

Placental growth factor (PlGF or PGF) is a member of the VEGF family, which is known to play a critical role in pathological angiogenesis, inflammation, and endothelial cell barrier function. However, the molecular mechanisms by which PlGF mediates its effects in non-proliferative diabetic retinopathy (DR) remain elusive. In this study, we performed transcriptome-wide profiling of differential gene expression for human retinal endothelial cells (HRECs) treated with PlGF antibody. The effect of antibody treatment on the samples was validated using trans-endothelial electric resistance (TEER), and western blot. A total of 3760 genes (1750 upregulated and 2010 downregulated) were found to be differentially expressed between the control and PlGF antibody treatment group. These differentially expressed genes (DEGs) were used for gene ontology and enrichment analysis to identify gene function, signal pathway, and interaction networks. The gene ontology results revealed that catalytic activity (GO:0003824) of molecular function, cell (GO:0005623) of the cellular component, and cellular process (GO:0009987) were among the most enriched biological processes. Pathways such as TGF-β, VEGF-VEGFR2, p53, apoptosis, pentose phosphate pathway, and ubiquitin-proteasome pathway, were among the most enriched, and TGF-β1 was identified as a primary upstream regulator. These data provide new insights into the underlying molecular mechanisms of PlGF in mediating biological functions, in relation to DR.

The Drosophila imaginal disc has been an excellent model for the study of developmental gene regulation. In particular, long non-coding RNAs (lncRNAs) have gained widespread attention in recent years due to their important role in gene regulation. Their specific spatiotemporal expressions further support their role in developmental processes and diseases. In this study, we explored the role of a novel lncRNA in Drosophila leg development by dissecting and dissociating w1118 third-instar larval third leg (L3) discs into single cells and single nuclei, and performing single-cell RNA-sequencing (scRNA-seq) and single-cell assays for transposase-accessible chromatin (scATAC-seq). Single-cell transcriptomics analysis of the L3 discs across three developmental timepoints revealed different cell types and identified lncRNA:CR33938 as a distal specific gene with high expression in late development. This was further validated by fluorescence in-situ hybridization (FISH). The scATAC-seq results reproduced the single-cell transcriptomics landscape and elucidated the distal cell functions at different timepoints. Furthermore, overexpression of lncRNA:CR33938 in the S2 cell line increased the expression of leg development genes, further confirming its important role in development.

High rates of phosphorus (P) currently being applied to soils for the production of vegetables in the Mekong Delta, Vietnam has led to the concern of its negative effect on the economics and the environment. This research presents a comprehensive study on the determination of P supplying capacity in this region of Vietnam to examine the possibility of reducing P fertilizer input. One hundred twenty (120) soil samples were collected to evaluate total P and Bray 1 available P in the soils. Phosphorus maximum sorption, degree of P saturation, P release, and the effect of P fertilizer on corn (Zea mays L.) yield in greenhouse and fields were also determined. Total P concentrations of 56.7% soil samples evaluated yielded high P concentrations (>560 mg P/kg), while 74.2% of the samples had high Bray 1 available P concentrations (>20 mg P/kg soil). Maximum P sorption ranged from 149 to 555 mg P/kg soil, respectively and has negative correlation to available P (r = - 0.63*), degrees of P saturation ranged from 0.63 to 5.46% correlated to available P (r = 0.98**) and maximum P release ranged from 1.2 to 61.9 mg P/kg soil, respectively correlated to available P (r = 0.96**). Corn grown in soils with available P concentrations >15 mg P/kg did not respond to P fertilizer in greenhouse or field experiments. We conclude that many farmers in this region can reduce P fertilizer input, thus increasing their profits and reducing negative environmental impacts associated with excess soil P for sustainable agriculture.

The snATAC + snRNA platform allows epigenomic profiling of open chromatin and gene expression with single-cell resolution. The most critical assay step is to isolate high-quality nuclei to proceed with droplet-base single nuclei isolation and barcoding. With the increasing popularity of multiomic profiling in various fields, there is a need for optimized and reliable nuclei isolation methods, mainly for human tissue samples. Herein we compared different nuclei isolation methods for cell suspensions, such as peripheral blood mononuclear cells (PBMC, n=18) and a solid tumor type, ovarian cancer (OC) (n=18), derived from debulking surgery. Nuclei morphology and sequencing output parameters were used to evaluate the quality of preparation. Our results show that NP-40 detergent-based nuclei isolation yields better sequencing results than collagenase tissue dissociation for OC, significantly impacting cell identification and analysis. Given the utility of applying such techniques to frozen samples, we also tested frozen preparation and digestion (n=6). A paired comparison between frozen and fresh samples validated the quality of both specimens. Finally, we demonstrate the reproducibility of scRNA and snATAC + snRNA platform, by comparing the gene expression profiling of PBMC. Our results highlight how the choice of nuclei isolation methods is critical for obtaining quality data in multiomic assays. It also shows that the measurement of expression between scRNA and snRNA is comparable and effective for cell type identification.