mRNA-lncRNA gene expression signature in HPV-associated neoplasia and cervical cancer.

Cervical cancer is one of the most frequent cancers in women and is associated with human papillomavirus (HPV) in 70 % of cases. Cervical cancer occurs because of progression of low-differentiated cervical intraepithelial neoplasia through grade 2 and 3 lesions. Along with the protein-coding genes, long noncoding RNAs (lncRNAs) play an important role in the development of malignant cell transformation. Although human papillomavirus is widespread, there is currently no well-characterized transcriptomic signature to predict whether this tumor will develop in the presence of HPV-associated neoplastic changes in the cervical epithelium. Changes in gene activity in tumors reflect the biological diversity of cellular phenotype and physiological functions and can be an important diagnostic marker. We performed comparative transcriptome analysis using open RNA sequencing data to assess differentially expressed genes between normal tissue, neoplastic epithelium, and cervical cancer. Raw data were preprocessed using the Galaxy platform. Batch effect correction, identification of differentially expressed genes, and gene set enrichment analysis (GSEA) were performed using R programming language packages. Subcellular localization of lncRNA was analyzed using Locate-R and iLoc-LncRNA 2.0 web services. 1,572 differentially expressed genes (DEGs) were recorded in the "cancer vs. control" comparison, and 1,260 DEGs were recorded in the "cancer vs. neoplasia" comparison. Only two genes were observed to be differentially expressed in the "neoplasia vs. control" comparison. The search for common genes among the most strongly differentially expressed genes among all comparison groups resulted in the identification of an expression signature consisting of the CCL20, CDKN2A, CTCFL, piR-55219, TRH, SLC27A6 and EPHA5 genes. The transcription level of the CCL20 and CDKN2A genes becomes increased at the stage of neoplastic epithelial changes and stays so in cervical cancer. Validation on an independent microarray dataset showed that the differential expression patterns of the CDKN2A and SLC27A6 genes were conserved in the respective gene expression comparisons between groups.

Cross-talk of MLST and transcriptome unveiling antibiotic resistance mechanism of carbapenem resistance Acinetobacter baumannii clinical strains isolated in Guiyang, China.

Introduction: Acinetobacter baumannii (A. baumannii) is an important opportunistic pathogen causing nosocomial infection in the clinic. The occurrence rate of antibiotic resistance is increasing year by year, resulting in a highly serious situation of bacterial resistance. Methods: To better understand the local epidemiology of multidrug-resistant A. baumannii, an investigation was conducted on the antibiotic resistance of different types of A. baumannii and its relationship with the genes of A. baumannii. Furthermore, the molecular mechanism underlying antibiotic resistance in A. baumannii was investigated through transcriptome analysis. Results: These results showed that a total of 9 STs were detected. It was found that 99% of the strains isolated in the hospital belonged to the same STs, and the clone complex CC208 was widely distributed in various departments and all kinds of samples. Furthermore, these A. baumannii strains showed high resistance to ertapenem, biapenem, meropenem, and imipenem, among which the resistance to ertapenem was the strongest. The detection rate of bla OXA-51 gene in these carbapenem resistance A. baumannii (CRAB) reached 100%; Additionally, the transcriptome results showed that the resistance genes were up-regulated in resistance strains, and these genes involved in biofilm formation, efflux pumps, peptidoglycan biosynthesis, and chaperonin synthesis. Discussion: These results suggest that the CC208 STs were the main clonal complex, and showed high carbapenem antibiotic resistance. All these resistant strains were distributed in various departments, but most of them were distributed in intensive care units (ICU). The bla OXA-23 was the main antibiotic resistance genotype; In summary, the epidemic trend of clinical A. baumannii in Guiyang, China was analyzed from the molecular level, and the resistance mechanism of A. baumannii to carbapenem antibiotics was analyzed with transcriptome, which provided a theoretical basis for better control of A. baumannii.

Transcriptome Analysis of the Striatum of Electroacupuncture-treated Naïve and Ischemic Stroke Mice.

Objectives: Electroacupuncture (EA) has been demonstrated to aid stroke recovery. However, few investigations have focused on identifying the potent molecular targets of EA by comparing EA stimulation between naïve and disease models. Therefore, this study was undertaken to identify the potent molecular therapeutic mechanisms underlying EA stimulation in ischemic stroke through a comparison of mRNA sequencing data obtained from EA-treated naïve control and ischemic stroke mouse models. Methods: Using both naïve control and middle cerebral artery occlusion (MCAO) mouse models, EA stimulation was administered at two acupoints, Baihui (GV20) and Dazhui (GV14), at a frequency of 2 Hz. Comprehensive assessments were conducted, including behavioral evaluations, RNA sequencing to identify differentially expressed genes (DEGs), functional enrichment analysis, protein-protein interaction (PPI) network analysis, and quantitative real-time PCR. Results: EA stimulation ameliorated the ischemic insult-induced motor dysfunction in mice with ischemic stroke. Comparative analysis between control vs. MCAO, control vs. control + EA, and MCAO vs. MCAO + EA revealed 4,407, 101, and 82 DEGs, respectively. Of these, 30, 7, and 1 were common across the respective groups. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed upregulated DEGs associated with the regulation of inflammatory immune response in the MCAO vs. MCAO + EA comparison. Conversely, downregulated DEGs in the control vs. control + EA comparison were linked to neuronal development. PPI analysis revealed major clustering related to the regulation of cytokines, such as Cxcl9, Pcp2, Ccl11, and Cxcl13, in the common DEGs of MCAO vs. MCAO + EA, with Esp8l1 identified as the only common downregulated DEG in both EA-treated naïve and ischemic models. Conclusion: These findings underscore the diverse potent mechanisms of EA stimulation between naïve and ischemic stroke mice, albeit with few overlaps. However, the potent mechanisms underlying EA treatment in ischemic stroke models were associated with the regulation of inflammatory processes involving cytokines.

Transcriptomic and proteomic elucidation of Z chromosome dosage compensation in Helicoverpa armigera.

Transcriptomic data have been used to study sex chromosome dosage compensation (SCDC) in approximately 10 Lepidoptera ZW species, yielding a consensus compensation pattern of Z ≈ ZZ < AA . $$ \approx \mathrm{ZZ}<\mathrm{AA}. $$ It remains unclear whether this compensation pattern holds when examining more Lepidoptera ZW species and/or using proteomic data to analyse SCDC. Here we combined transcriptomic and proteomic data as well as transcriptional level of six individual Z genes to reveal the SCDC pattern in Helicoverpa armigera, a polyphagous lepidopteran pest of economic importance. Transcriptomic analysis showed that the Z chromosome expression of H. armigera was balanced between male and female but substantially reduced relative to autosome expression, exhibiting an SCDC pattern of Z ≈ ZZ < AA $$ \approx \mathrm{ZZ}<\mathrm{AA} $$ . When using H. amigera midgut proteomic data, the SCDC pattern of this species changed from Z ≈ ZZ < AA $$ \approx \mathrm{ZZ}<\mathrm{AA} $$ at transcriptomic level to Z = ZZ = AA at the proteomic level. RT-qPCR analysis of transcript abundance of six Z genes found that compensation for each Z gene could vary from no compensation to overcompensation, depending on the individual genes and tissues tested. These results demonstrate for the first time the existence of a translational compensation mechanism, which is operating in addition to a translational mechanism, such as has been reported in other lepidopteran species. And the transcriptional compensation mechanism functions to accomplish Z chromosome dosage balance between the sexes (M = F on the Z chromosome), whereas the translation compensation mechanism operates to achieve dosage compensation between Z chromosome and autosome (Z = AA).

Adjuvant Effects of a CC Chemokine for Enhancing the Efficacy of an Inactivated Streptococcus agalactiae Vaccine in Nile Tilapia (Oreochromis niloticus).

In this study, the ability of a CC chemokine (On-CC1) adjuvant to enhance the efficacy of a formalin-killed Streptococcus agalactiae vaccine (WC) in inducing immune responses against S. agalactiae in Nile tilapia was investigated through immune-related gene expression analysis, enzyme-linked immunosorbent assay (ELISA), transcriptome sequencing, and challenge tests. Significantly higher S. agalactiae-specific IgM levels were detected in fish in the WC+CC group than in the WC alone or control groups at 8 days postvaccination (dpv). The WC vaccine group exhibited increased specific IgM levels at 15 dpv, comparable to those of the WC+CC group, with sustained higher levels observed in the latter group at 29 dpv and after challenge with S. agalactiae for 14 days. Immune-related gene expression analysis revealed upregulation of all target genes in the control group compared to those in the vaccinated groups, with notable differences between the WC and WC+CC groups at various time intervals. Additionally, transcriptome analysis revealed differential gene expression profiles between the vaccinated (24 and 96 hpv) and control groups, with notable upregulation of immune-related genes in the vaccinated fish. Differential gene expression (DGE) analysis revealed significant upregulation of immunoglobulin and other immune-related genes in the control group compared to those in the vaccinated groups (24 and 96 hpv), with distinct patterns observed between the WC and WC+CC vaccine groups. Finally, challenge with a virulent strain of S. agalactiae resulted in significantly higher survival rates for fish in the WC and WC+CC groups compared to fish in the control group, with a notable increase in survival observed in fish in the WC+CC group.

scPLAN: a hierarchical computational framework for single transcriptomics data annotation, integration and cell-type label refinement.

MOTIVATION: In the past decade, single-cell RNA sequencing (scRNA-seq) has emerged as a pivotal method for transcriptomic profiling in biomedical research. Precise cell-type identification is crucial for subsequent analysis of single-cell data. And the integration and refinement of annotated data are essential for building comprehensive databases. However, prevailing annotation techniques often overlook the hierarchical organization of cell types, resulting in inconsistent annotations. Meanwhile, most existing integration approaches fail to integrate datasets with different annotation depths and none of them can enhance the labels of outdated data with lower annotation resolutions using more intricately annotated datasets or novel biological findings. RESULTS: Here, we introduce scPLAN, a hierarchical computational framework designed for scRNA-seq data analysis. scPLAN excels in annotating unlabeled scRNA-seq data using a reference dataset structured along a hierarchical cell-type tree. It identifies potential novel cell types in a systematic, layer-by-layer manner. Additionally, scPLAN effectively integrates annotated scRNA-seq datasets with varying levels of annotation depth, ensuring consistent refinement of cell-type labels across datasets with lower resolutions. Through extensive annotation and novel cell detection experiments, scPLAN has demonstrated its efficacy. Two case studies have been conducted to showcase how scPLAN integrates datasets with diverse cell-type label resolutions and refine their cell-type labels. AVAILABILITY: https://github.com/michaelGuo1204/scPLAN.

A novel plant growth regulator B2 mediates drought resistance by regulating reactive oxygen species, phytohormone signaling, phenylpropanoid biosynthesis, and starch metabolism pathways in Carex breviculmis.

Drought is one of the most common environmental stressors that severely threatens plant growth, development, and productivity. B2 (2,4-dichloroformamide cyclopropane acid), a novel plant growth regulator, plays an essential role in drought adaptation, significantly enhancing the tolerance of Carex breviculmis seedlings. Its beneficial effects include improved ornamental value, sustained chlorophyll content, increased leaf dry weight, elevated relative water content, and enhanced root activity under drought conditions. B2 also directly scavenges hydrogen peroxide and superoxide anion contents while indirectly enhancing the activities of antioxidant enzymes (superoxide dismutase, peroxidase, catalase, and ascorbate peroxidase) to detoxify reactive oxygen species (ROS) oxidative damage. Transcriptome analysis demonstrated that B2 activates drought-responsive transcription factors (AP2/ERF-ERF, WRKY, and mTERF), leading to significant upregulation of genes associated with phenylpropanoid biosynthesis (HCT, POD, and COMT). Additionally, these transcription factors were found to suppress the degradation of starch. B2 regulates phytohormone signaling related-genes, leading to an increase in abscisic acid contents in drought-stressed plants. Collectively, these findings offer new insights into the intricate mechanisms underlying C. breviculmis' resistance to drought damage, highlighting the potential application of B2 for future turfgrass establishment and management with enhanced drought tolerance.

Protective efficacy and immune responses of largemouth bass (Micropterus salmoides) immunized with an inactivated vaccine against the viral hemorrhagic septicemia virus genotype IVa.

Viral hemorrhagic septicemia virus (VHSV) poses a significant threat to the aquaculture industry, prompting the need for effective preventive measures. Here, we developed an inactivated VHSV and revealed the molecular mechanisms underlying the host's protective response against VHSV. The vaccine was created by treating VHSV with 0.05% formalin at 16°C for 48 h, which was determined to be the most effective inactivation method. Compared with nonvaccinated fish, vaccinated fish exhibited a remarkable increase in survival rate (99%) and elevated levels of serum neutralizing antibodies, indicating strong immunization. To investigate the gene changes induced by vaccination, RNA sequencing was performed on spleen samples from control and vaccinated fish 14 days after vaccination. The analysis revealed 893 differentially expressed genes (DEGs), with notable upregulation of immune-related genes such as annexin A1a, coxsackievirus and adenovirus receptor homolog, V-set domain-containing T-cell activation inhibitor 1-like, and heat shock protein 90 alpha class A member 1 tandem duplicate 2, indicating a vigorous innate immune response. Furthermore, KEGG enrichment analysis highlighted significant enrichment of DEGs in processes related to antigen processing and presentation, necroptosis, and viral carcinogenesis. GO enrichment analysis further revealed enrichment of DEGs related to the regulation of type I interferon (IFN) production, type I IFN production, and negative regulation of viral processes. Moreover, protein-protein interaction network analysis identified central hub genes, including IRF3 and HSP90AA1.2, suggesting their crucial roles in coordinating the immune response elicited by the vaccine. These findings not only confirm the effectiveness of our vaccine formulation but also offer valuable insights into the underlying immunological mechanisms, which can be valuable for future vaccine development and disease management in the aquaculture industry.

Transcriptome Profiling Reveals Distinct Alterations in B-Cell Signature and Dysregulation of Peripheral B-Cell Subsets in Sickle Cell Anemia Patients.

INTRODUCTION: Sickle cell anemia (SCA) is characterized by immune system activation and heightened susceptibility to infections. We hypothesized that SCA patients exhibit transcriptional alterations in B-cell-related genes, impacting their peripheral B-cell compartment and leading to dysregulated humoral immunity and increased infection susceptibility. OBJECTIVE: Our objective was to conduct an in-silico analysis of whole blood transcriptomes from SCA patients and healthy controls obtained from public repositories. We aimed to identify alterations in the adaptive immune system and validate these findings in our own SCA patient cohort. METHODS: Bioinformatic analyses unveiled significant transcriptional alterations in B-cell signatures, developmental pathways, and signaling pathways. These results were validated in peripheral blood mononuclear cells from our SCA patient cohort and controls using real-time PCR and flow cytometry. RESULTS: Ninety genes exhibited differential expression, with 70 upregulated and 20 downregulated. Dysregulation in the B-cell compartment of SCA patients was evident, characterized by increased frequencies of immature and naive B-cells and decreased percentages of memory B-cell subsets compared to healthy controls. CONCLUSION: Our findings highlight previously unexplored transcriptional and quantitative alterations in peripheral B-cells among SCA patients. Understanding these changes sheds light on the mechanisms contributing to the heightened infection risk in this population. Future studies should delve deeper into these molecular changes to develop targeted interventions and therapeutic strategies aimed at mitigating infection susceptibility in individuals with SCA.

Exploring the role of CD8+ T cells in clear renal cell carcinoma metastasis.

Clear cell renal cell carcinoma (ccRCC) accounts for approximately 75-80% of all patients with renal cell carcinoma. Despite its prevalence, little is known regarding the key components involved in ccRCC metastasis. In this study, scRNA-seq analysis was employed to classify CD8+ T cells into four sub-clusters based on their genetic profiles and immunofluorescence experiments were used to validate two key clusters. Through gene set enrichment analysis, these newly identified sub-clusters were found to exhibit distinct biological characteristics. Notably, TYMP, TOP2A, CHI3L2, CDKN3, CENPM, and RZH2 were highly expressed in these sub-clusters, indicating a correlation with poor prognosis. Among these sub-clusters, CD8+ T cells (MT-ND4) were identified as potentially playing a critical role in mediating ccRCC metastasis. These results contribute to our understanding of CD8+ T cell heterogeneity in ccRCC and shed light on the mechanisms underlying the loss of immune response against cancer.

Odorant receptor co-receptors affect expression of tuning receptors in Drosophila.

Insects detect odorants using two large families of heteromeric receptors, the Odorant Receptors (ORs) and Ionotropic Receptors (IRs). Most OR and IR genes encode odorant-binding "tuning" subunits, whereas four (Orco, Ir8a, Ir25a, and Ir76b) encode co-receptor subunits required for receptor function. Olfactory neurons are thought to degenerate in the absence of Orco in ants and bees, and limited data suggest this may happen to some olfactory neurons in Drosophila fruit flies as well. Here, we thoroughly examined the role of co-receptors on olfactory neuron survival in Drosophila. Leveraging knowledge that olfactory neuron classes are defined by the expression of different tuning receptors, we used tuning receptor expression in antennal transcriptomes as a proxy for the survival of distinct olfactory neuron classes. Consistent with olfactory neuron degeneration, expression of many OR-family tuning receptors is decreased in Orco mutants relative to controls, and transcript loss is progressive with age. The effects of Orco are highly receptor-dependent, with expression of some receptor transcripts nearly eliminated and others unaffected. Surprisingly, further studies revealed that olfactory neuron classes with reduced tuning receptor expression generally survive in Orco mutant flies. Furthermore, there is little apoptosis or neuronal loss in the antenna of these flies. We went on to investigate the effects of IR family co-receptor mutants using similar approaches and found that expression of IR tuning receptors is decreased in the absence of Ir8a and Ir25a, but not Ir76b. As in Orco mutants, Ir8a-dependent olfactory neurons mostly endure despite near-absent expression of associated tuning receptors. Finally, we used differential expression analysis to identify other antennal genes whose expression is changed in IR and OR co-receptor mutants. Taken together, our data indicate that odorant co-receptors are necessary for maintaining expression of many tuning receptors at the mRNA level. Further, most Drosophila olfactory neurons persist in OR and IR co-receptor mutants, suggesting that the impact of co-receptors on neuronal survival may vary across insect species.

Cost-Effective Whole Transcriptome Sequencing Landscape and Diagnostic Potential Biomarkers in Active Tuberculosis.

Tuberculosis (TB) is a prevalent and severe infectious disease that poses a significant threat to human health. However, it is frequently disregarded as there are not enough quick and accurate ways to diagnose tuberculosis. Here, we develop a strategy for tuberculosis detection to address the challenges, including an experimental strategy, namely, Double Adapter Directional Capture sequencing (DADCSeq), an easily operated and low-cost whole transcriptome sequencing method, and a computational method to identify hub differentially expressed genes as well as the diagnosis of TB based on whole transcriptome data using DADCSeq on peripheral blood mononuclear cells (PBMCs) from active TB and latent TB or healthy control. Applying our approach to create a robust and stable TB multi-mRNA risk probability model (TBMMRP) that can accurately distinguish active and latent TB patients, including active TB and healthy controls in clinical cohorts, this diagnostic biomarker was successfully validated by several independent cross-platform cohorts with favorable performance in differentiating active TB from latent TB or active TB from healthy controls and further demonstrated superior or similar diagnostic accuracy compared to previous diagnostic markers. Overall, we develop a low-cost and effective strategy for tuberculosis diagnosis; as the clinical cohort increases, we can expand to different disease kinds and learn new features through our disease diagnosis strategy.

Effect of static magnetic field on gene expression of human umbilical cord mesenchymal stem cells by transcriptome analysis.

PURPOSE: Static magnetic fields (SMFs) induce various biological reactions and have been applied in the biological therapy of diseases, especially in combination with mesenchymal stem cells (MSCs) and tissue engineering. However, the underlying influence of SMFs on MSCs gene expression remains largely unclear. In this study, we aim to investigate the effects of SMFs on gene expression of human MSCs. MATERIALS AND METHODS: We exposed human MSCs to two different intensities (0.35 â€‹T and 1.0 â€‹T) of SMFs and observed the effects of SMFs on cell morphology. Subsequently, RNA-sequencing was performed to explore the gene expression changes. RESULTS: Compared with control group cells, no significant differences in cell morphology were observed under a phase contrast inverted microscope, but the transcriptome of SMF-exposed MSCs were significantly changed in both 0.35 â€‹T and 1.0 â€‹T groups and the differential expressed genes are involved in multiple pathways, such as ubiquitin mediated proteolysis, TNF signaling pathway, NF-kappa B signaling pathway, TGF-beta signaling pathway, metabolic pathways, and apoptosis, which regulate the biological functions of MSCs. CONCLUSIONS: SMFs stimulation could affect the gene expression of human MSCs and the biological effects vary by the different intensities of SMFs. These data offer the molecular foundation for future application of SMFs in stem cell technology as well as tissue engineering medicine.

Comparative analysis of testicular fusion in Spodoptera litura (cutworm) and Bombyx mori (silkworm): Histological and transcriptomic insights.

Spodoptera litura commonly known as the cutworm, is among the most destructive lepidopteran pests affecting over 120 plants species. The powerful destructive nature of this lepidopteran is attributable to its high reproductive capacity. The testicular fusion that occurs during metamorphosis from larvae to pupa in S.litura positively influences the reproductive success of the offspring. In contrast, Bombyx mori, the silkworm, retains separate testes throughout its life and does not undergo this fusion process. Microscopic examination reveals that during testicular fusion in S.litura, the peritoneal sheath becomes thinner and more translucent, whereas in B.mori, the analogous region thickens. The outer basement membrane in S.litura exhibits fractures, discontinuity, and uneven thickness accompanied by a significant presence of cellular secretions, large cell size, increased vesicles, liquid droplets, and a proliferation of rough endoplasmic reticulum and mitochondria. In contrast, the testicular peritoneal sheath of B.mori at comparable developmental stage exhibits minimal change. Comparative transcriptomic analysis of the testicular peritoneal sheath reveals a substantial difference in gene expression between the two species. The disparity in differential expressed genes (DEGs) is linked to an enrichment of numerous transcription factors, intracellular signaling pathways involving Ca2+ and GTPase, as well as intracellular protein transport and signaling pathways. Meanwhile, structural proteins including actin, chitin-binding proteins, membrane protein fractions, cell adhesion, extracellular matrix proteins are predominantly identified. Moreover, the study highlights the enrichment of endopeptidases, serine proteases, proteolytic enzymes and matrix metalloproteins, which may play a role in the degradation of the outer membrane. Five transcription factors-Slforkhead, Slproline, Slcyclic, Slsilk, and SlD-ETS were identified, and their expression pattern were confirmed by qRT-PCR. they are candidates for participating in the regulation of testicular fusion. Our findings underscore significant morphological and trancriptomic variation in the testicular peritoneal sheath of S.litura compared to the silkworm, with substantial changes at the transcriptomic level coinciding with testicular fusion. The research provides valuable clues for understanding the complex mechanisms underlying this unique phenomenon in insects.

Transcriptome analysis of the Japanese eel (Anguilla japonica) during larval metamorphosis.

BACKGROUND: Anguillid eels spend their larval period as leptocephalus larvae that have a unique and specialized body form with leaf-like and transparent features, and they undergo drastic metamorphosis to juvenile glass eels. Less is known about the transition of leptocephali to the glass eel stage, because it is difficult to catch the metamorphosing larvae in the open ocean. However, recent advances in rearing techniques for the Japanese eel have made it possible to study the larval metamorphosis of anguillid eels. In the present study, we investigated the dynamics of gene expression during the metamorphosis of Japanese eel leptocephali using RNA sequencing. RESULTS: During metamorphosis, Japanese eels were classified into 7 developmental stages according to their morphological characteristics, and RNA sequencing was used to collect gene expression data from each stage. A total of 354.8 million clean reads were generated from the body and 365.5 million from the head, after the processing of raw reads. For filtering of genes that characterize developmental stages, a classification model created by a Random Forest algorithm was built. Using the importance of explanatory variables feature obtained from the created model, we identified 46 genes selected in the body and 169 genes selected in the head that were defined as the "most characteristic genes" during eel metamorphosis. Next, network analysis and subsequently gene clustering were conducted using the most characteristic genes and their correlated genes, and then 6 clusters in the body and 5 clusters in the head were constructed. Then, the characteristics of the clusters were revealed by Gene Ontology (GO) enrichment analysis. The expression patterns and GO terms of each stage were consistent with previous observations and experiments during the larval metamorphosis of the Japanese eel. CONCLUSION: Genome and transcriptome resources have been generated for metamorphosing Japanese eels. Genes that characterized metamorphosis of the Japanese eel were identified through statistical modeling by a Random Forest algorithm. The functions of these genes were consistent with previous observations and experiments during the metamorphosis of anguillid eels.

Integrated transcriptomic and metabolomic analysis provides insight into the pollen development of CMS-D1 rice.

BACKGROUND: Cytoplasmic male sterility (CMS) has greatly improved the utilization of heterosis in crops due to the absence of functional male gametophyte. The newly developed sporophytic D1 type CMS (CMS-D1) rice exhibits unique characteristics compared to the well-known sporophytic CMS-WA line, making it a valuable resource for rice breeding. RESULTS: In this research, a novel CMS-D1 line named Xingye A (XYA) was established, characterized by small, transparent, and shriveled anthers. Histological and terminal deoxynucleotidyl transferase-mediated dUTP nick-end labeling (TUNEL) assays conducted on anthers from XYA and its maintainer line XYB revealed that male sterility in XYA is a result of delayed degradation of tapetal cells and abnormal programmed cell death (PCD) of microspores. Transcriptome analysis of young panicles revealed that differentially expressed genes (DEGs) in XYA, compared to XYB, were significantly enriched in processes related to chromatin structure and nucleosomes during the microspore mother cell (MMC) stage. Conversely, processes associated with sporopollenin biosynthesis, pollen exine formation, chitinase activity, and pollen wall assembly were enriched during the meiosis stage. Metabolome analysis identified 176 specific differentially accumulated metabolites (DAMs) during the meiosis stage, enriched in pathways such as α-linoleic acid metabolism, flavone and flavonol biosynthesis, and linolenic acid metabolism. Integration of transcriptomic and metabolomic data underscored the jasmonic acid (JA) biosynthesis pathway was significant enriched in XYA during the meiosis stage compared to XYB. Furthermore, levels of JA, MeJA, OPC4, OPDA, and JA-Ile were all higher in XYA than in XYB at the meiosis stage. CONCLUSIONS: These findings emphasize the involvement of the JA biosynthetic pathway in pollen development in the CMS-D1 line, providing a foundation for further exploration of the molecular mechanisms involved in CMS-D1 sterility.

Transcriptomic analysis reveals the regulatory mechanisms of messenger RNA (mRNA) and long non-coding RNA (lncRNA) in response to waterlogging stress in rye (Secale cereale L.).

BACKGROUND: Waterlogging stress (WS) negatively impacts crop growth and productivity, making it important to understand crop resistance processes and discover useful WS resistance genes. In this study, rye cultivars and wild rye species were subjected to 12-day WS treatment, and the cultivar Secale cereale L. Imperil showed higher tolerance. Whole transcriptome sequencing was performed on this cultivar to identify differentially expressed (DE) messenger RNAs (DE-mRNAs) and long non-coding RNAs (DE-lncRNAs) involved in WS response. RESULTS: Among the 6 species, Secale cereale L. Imperil showed higher tolerance than wild rye species against WS. The cultivar effectively mitigated oxidative stress, and regulated hydrogen peroxide and superoxide anion. A total of 728 DE-mRNAs and 60 DE-lncRNAs were discovered. Among these, 318 DE-mRNAs and 32 DE-lncRNAs were upregulated, and 410 DE-mRNAs and 28 DE-lncRNAs were downregulated. GO enrichment analysis discovered metabolic processes, cellular processes, and single-organism processes as enriched biological processes (BP). For cellular components (CC), the enriched terms were membrane, membrane part, cell, and cell part. Enriched molecular functions (MF) terms were catalytic activity, binding, and transporter activity. LncRNA and mRNA regulatory processes were mainly related to MAPK signaling pathway-plant, plant hormone signal transduction, phenylpropanoid biosynthesis, anthocyanin biosynthesis, glutathione metabolism, ubiquitin-mediated proteolysis, ABC transporter, Cytochrome b6/f complex, secondary metabolite biosynthesis, and carotenoid biosynthesis pathways. The signalling of ethylene-related pathways was not mainly dependent on AP2/ERF and WRKY transcription factors (TF), but on other factors. Photosynthetic activity was active, and carotenoid levels increased in rye under WS. Sphingolipids, the cytochrome b6/f complex, and glutamate are involved in rye WS response. Sucrose transportation was not significantly inhibited, and sucrose breakdown occurs in rye under WS. CONCLUSIONS: This study investigated the expression levels and regulatory functions of mRNAs and lncRNAs in 12-day waterlogged rye seedlings. The findings shed light on the genes that play a significant role in rye ability to withstand WS. The findings from this study will serve as a foundation for further investigations into the mRNA and lncRNA WS responses in rye.

The nuclear and mitochondrial genome assemblies of Tetragonisca angustula (Apidae: Meliponini), a tiny yet remarkable pollinator in the Neotropics.

BACKGROUND: The field of bee genomics has considerably advanced in recent years, however, the most diverse group of honey producers on the planet, the stingless bees, are still largely neglected. In fact, only eleven of the ~ 600 described stingless bee species have been sequenced, and only three using a long-read (LR) sequencing technology. Here, we sequenced the nuclear and mitochondrial genomes of the most common, widespread and broadly reared stingless bee in Brazil and other neotropical countries-Tetragonisca angustula (popularly known in Brazil as jataí). RESULTS: A total of 48.01 Gb of DNA data were generated, including 2.31 Gb of Pacific Bioscience HiFi reads and 45.70 Gb of Illumina short reads (SRs). Our preferred assembly comprised 683 contigs encompassing 284.49 Mb, 62.84 Mb of which (22.09%) corresponded to 445,793 repetitive elements. N50, L50 and complete BUSCOs reached 1.02 Mb, 91 contigs and 97.1%, respectively. We predicted that the genome of T. angustula comprises 17,459 protein-coding genes and 4,108 non-coding RNAs. The mitogenome consisted of 17,410 bp, and all 37 genes were found to be on the positive strand, an unusual feature among bees. A phylogenomic analysis of 26 hymenopteran species revealed that six odorant receptor orthogroups of T. angustula were found to be experiencing rapid evolution, four of them undergoing significant contractions. CONCLUSIONS: Here, we provided the first nuclear and mitochondrial genome assemblies for the ecologically and economically important T. angustula, the fourth stingless bee species to be sequenced with LR technology thus far. We demonstrated that even relatively small amounts of LR data in combination with sufficient SR data can yield high-quality genome assemblies for bees.

Combining transcriptome and metabolome analysis to understand the response of sorghum to Melanaphis sacchari.

BACKGROUND: The sorghum aphid Melanaphis sacchari (Zehntner) (Homoptera: Aphididae) is an important insect in the late growth phase of sorghum (Sorghum bicolor L.). However, the mechanisms of sorghum response to aphid infestation are unclear. RESULTS: In this paper, the mechanisms of aphid resistance in different types of sorghum varieties were revealed by studying the epidermal cell structure and performing a transcriptome and metabolome association analysis of aphid-resistant and aphid-susceptible varieties. The epidermal cell results showed that the resistance of sorghum to aphids was positively correlated with epidermal cell regularity and negatively correlated with the intercellular space and leaf thickness. Transcriptome and metabolomic analyses showed that differentially expressed genes in the resistant variety HN16 and susceptible variety BTX623 were mainly enriched in the flavonoid biosynthesis pathway and differentially expressed metabolites were mainly related to isoflavonoid biosynthesis and flavonoid biosynthesis. The q-PCR results of key genes were consistent with the transcriptome expression results. Meanwhile, the metabolome test results showed that after aphidinfestation, naringenin and genistein were significantly upregulated in the aphid-resistant variety HN16 and aphid-susceptible variety BTX623 while luteolin was only significantly upregulated in BTX623. These results show that naringenin, genistein, and luteolin play important roles in plant resistance to aphid infestation. The results of exogenous spraying tests showed that a 1‰ concentration of naringenin and genistein is optimal for improving sorghum resistance to aphid feeding. CONCLUSIONS: In summary, the physical properties of the sorghum leaf structure related to aphid resistance were studied to provide a reference for the breeding of aphid-resistant varieties. The flavonoid biosynthesis pathway plays an important role in the response of sorghum aphids and represents an important basis for the biological control of these pests. The results of the spraying experiment provide insights for developing anti-aphid substances in the future.

Whole transcriptome analysis revealed the regulatory network and related pathways of non-coding RNA regulating ovarian atrophy in broody hens.

Poultry broodiness can cause ovarian atresia, which has a detrimental impact on egg production. Non-coding RNAs (ncRNAs) have become one of the most talked-about topics in life sciences because of the increasing evidence of their novel biological roles in regulatory systems. However, the molecular mechanisms of ncRNAs functions and processes in chicken ovarian development remain largely unknown. Whole-transcriptome RNA sequencing of the ovaries of broodiness and laying chickens was thus performed to identify the ncRNA regulatory mechanisms associated with ovarian atresia in chickens. Subsequent analysis revealed that the ovaries of laying chickens and those with broodiness had 40 differentially expressed MicroRNA (miRNAs) (15 up-regulated and 25 down-regulated), 379 differentially expressed Long Noncoding RNA (lncRNAs) (213 up-regulated and 166 down-regulated), and 129 differentially expressed circular RNA (circRNAs) (63 up-regulated and 66 down-regulated). The competing endogenous RNAs (ceRNA) network analysis further revealed the involvement of ECM-receptor interaction, AGE-RAGE signaling pathway, focal adhesion, cytokine-cytokine receptor interaction, inflammatory mediator regulation of TRP channels, renin secretion, gap junction, insulin secretion, serotonergic synapse, and IL-17 signaling pathways in broodiness. Upon further analysis, it became evident that THBS1 and MYLK are significant candidate genes implicated in the regulation of broodiness. The expression of these genes is linked to miR-155-x, miR-211-z, miR-1682-z, gga-miR-155, and gga-miR-1682, as well as to the competitive binding of novel_circ_014674 and MSTRG.3306.4. The findings of this study reveal the existence of a regulatory link between non-coding RNAs and their competing mRNAs, which provide a better comprehension of the ncRNA function and processes in chicken ovarian development.