An RNA-Seq Protocol for Differential Expression Analysis

Exploring the role of CBLB in acute myocardial infarction: transcriptomic, microbiomic, and metabolomic analyses

BACKGROUND

Specific alterations in gut microbiota and metabolites have been linked to AMI, with CBLB potentially playing an essential role. However, the precise interactions remain understudied, creating a significant gap in our understanding. This study aims to address this by exploring these interactions in CBLB-intervened AMI mice using transcriptome sequencing, 16 S rDNA, and non-targeted metabolite analysis.

METHODS

To probe the therapeutic potential and mechanistic underpinnings of CBLB overexpression in AMI, we utilized an integrative multi-omics strategy encompassing transcriptomics, metabolomics, and 16s rDNA sequencing. We selected these particular methods as they facilitate a holistic comprehension of the intricate interplay between the host and its microbiota, and the potential effects on the host's metabolic and gene expression profiles. The uniqueness of our investigation stems from utilizing a multi-omics approach to illuminate the role of CBLB in AMI, an approach yet unreported to the best of our knowledge. Our experimental protocol encompassed transfection of CBLB lentivirus-packaged vectors into 293T cells, followed by subsequent intervention in AMI mice. Subsequently, we conducted pathological staining, fecal 16s rDNA sequencing, and serum non-targeted metabolome sequencing. We applied differential expression analysis to discern differentially expressed genes (DEGs), differential metabolites, and differential microbiota. We performed protein-protein interaction analysis to identify core genes, and conducted correlation studies to clarify the relationships amongst these core genes, paramount metabolites, and key microbiota.

RESULTS

Following the intervention of CBLB in AMI, we observed a significant decrease in inflammatory cell infiltration and collagen fiber formation in the infarcted region of mice hearts. We identified key changes in microbiota, metabolites, and DEGs that were associated with this intervention. The findings revealed that CBLB has a significant correlation with DEGs, differential metabolites and microbiota, respectively. This suggests it could play a pivotal role in the regulation of AMI.

CONCLUSION

This study confirmed the potential of differentially expressed genes, metabolites, and microbiota in AMI regulation post-CBLB intervention. Our findings lay groundwork for future exploration of CBLB's role in AMI, suggesting potential therapeutic applications and novel research directions in AMI treatment strategies.

Evolution of major flowering pathway integrators in Orchidaceae

The Orchidaceae is a mega-diverse plant family with ca. 29,000 species with a large variety of life forms that can colonize transitory habitats. Despite this diversity, little is known about their flowering integrators in response to specific environmental factors. During the reproductive transition in flowering plants a vegetative apical meristem (SAM) transforms into an inflorescence meristem (IM) that forms bracts and flowers. In model grasses, like rice, a flowering genetic regulatory network (FGRN) controlling reproductive transitions has been identified, but little is known in the Orchidaceae. In order to analyze the players of the FRGN in orchids, we performed comprehensive phylogenetic analyses of CONSTANS-like/CONSTANS-like 4 (COL/COL4), FLOWERING LOCUS D (FD), FLOWERING LOCUS C/FRUITFULL (FLC/FUL) and SUPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) gene lineages. In addition to PEBP and AGL24/SVP genes previously analyzed, here we identify an increase of orchid homologs belonging to COL4, and FUL gene lineages in comparison with other monocots, including grasses, due to orchid-specific gene lineage duplications. Contrariwise, local duplications in Orchidaceae are less frequent in the COL, FD and SOC1 gene lineages, which points to a retention of key functions under strong purifying selection in essential signaling factors. We also identified changes in the protein sequences after such duplications, variation in the evolutionary rates of resulting paralogous clades and targeted expression of isolated homologs in different orchids. Interestingly, vernalization-response genes like VERNALIZATION1 (VRN1) and FLOWERING LOCUS C (FLC) are completely lacking in orchids, or alternatively are reduced in number, as is the case of VERNALIZATION2/GHD7 (VRN2). Our findings point to non-canonical factors sensing temperature changes in orchids during reproductive transition. Expression data of key factors gathered from Elleanthus auratiacus, a terrestrial orchid in high Andean mountains allow us to characterize which copies are actually active during flowering. Altogether, our data lays down a comprehensive framework to assess gene function of a restricted number of homologs identified more likely playing key roles during the flowering transition, and the changes of the FGRN in neotropical orchids in comparison with temperate grasses.

Transcriptomic analysis reveals zinc-mediated virulence and pathogenicity in multidrug-resistant Acinetobacter baumannii

Acinetobacter baumannii is a gram-negative bacterium well known for its multidrug resistance and connection to nosocomial infections under ESKAPE pathogens. This opportunistic pathogen is ubiquitously associated with nosocomial infections, posing significant threats within healthcare environments. Its critical clinical symptoms, namely, meningitis, urinary tract infections, bloodstream infections, ventilator-associated pneumonia, and pneumonia, catalyze the imperative demand for innovative therapeutic interventions. The proposed research focuses on delineating the role of Zinc, a crucial metallo-binding protein and micronutrient integral to bacterial metabolism and virulence, to enhance understanding of the pathogenicity of A. baumannii. RNA sequencing and subsequent DESeq2 analytical methods were used to identify differential gene expressions influenced by zinc exposure. Exploiting the STRING database for functional enrichment analysis has demonstrated the complex molecular mechanisms underlying the enhancement of pathogenicity prompted by Zinc. Moreover, hub genes like gltB, ribD, AIL77834.1, sdhB, nuoI, acsA_1, acoC, accA, accD were predicted using the cytohubba tool in Cytoscape. This investigation underscores the pivotal role of Zinc in the virulence of A. baumannii elucidates the underlying molecular pathways responsible for its pathogenicity. The research further accentuates the need for innovative therapeutic strategies to combat A. baumannii infections, particularly those induced by multidrug-resistant strains.

A single-cell characterised signature integrating heterogeneity and microenvironment of lung adenocarcinoma for prognostic stratification

BACKGROUND

The high heterogeneity of tumour and the complexity of tumour microenvironment (TME) greatly impacted the tumour development and the prognosis of cancer in the era of immunotherapy. In this study, we aimed to portray the single cell-characterised landscape of lung adenocarcinoma (LUAD), and develop an integrated signature incorporating both tumour heterogeneity and TME for prognosis stratification.

METHODS

Single-cell tagged reverse transcription sequencing (STRT-seq) was performed on tumour tissues and matched normal tissues from 14 patients with LUAD for immune landscape depiction and candidate key genes selection for signature construction. Kaplan-Meier survival analyses and in-vitro cell experiments were conducted to confirm the gene functions. The transcriptomic profile of 1949 patients from 11 independent cohorts including nine public datasets and two in-house cohorts were obtained for validation.

FINDINGS

We selected 11 key genes closely related to cell-to-cell interaction, tumour development, T cell phenotype transformation, and Ma/Mo cell distribution, including HLA-DPB1, FAM83A, ITGB4, OAS1, FHL2, S100P, FSCN1, SFTPD, SPP1, DBH-AS1, CST3, and established an integrated 11-gene signature, stratifying patients to High-Score or Low-Score group for better or worse prognosis. Moreover, the prognostically-predictive potency of the signature was validated by 11 independent cohorts, and the immunotherapeutic predictive potency was also validated by our in-house cohort treated by immunotherapy. Additionally, the in-vitro cell experiments and drug sensitivity prediction further confirmed the gene function and generalizability of this signature across the entire RNA profile spectrum.

INTERPRETATION

This single cell-characterised 11-gene signature might offer insights for prognosis stratification and potential guidance for treatment selection.

FUNDING

Support for the study was provided by National key research and development project (2022YFC2505004, 2022YFC2505000 to Z.W. and J.W.), Beijing Natural Science Foundation (7242114 to J.X.), National Natural Science Foundation of China of China (82102886 to J.X., 81871889 and 82072586 to Z.W.), Beijing Nova Program (20220484119 to J.X.), NSFC general program (82272796 to J.W.), NSFC special program (82241229 to J.W.), CAMS Innovation Fund for Medical Sciences (2021-1-I2M-012, 2022-I2M-1-009 to Z.W. and J.W.), Beijing Natural Science Foundation (7212084 to Z.W.), CAMS Key lab of translational research on lung cancer (2018PT31035 to J.W.), Aiyou Foundation (KY201701 to J.W.). Medical Oncology Key Foundation of Cancer Hospital Chinese Academy of Medical Sciences (CICAMS-MOCP2022003 to J.X.).

Betaine attenuates age-related suppression in autophagy via Mettl21c/p97/VCP axis to delay muscle loss

Age-related impairment of autophagy accelerates muscle loss and lead to sarcopenia. Betaine can delay muscle loss as a dietary methyl donor via increasing S-adenosyl-L-methionine (SAM, a crucial metabolite for autophagy regulation) in methionion cycle. However, whether betaine can regulate autophagy level to attenuate degeneration in aging muscle remains unclear. Herein, male C57BL/6J young mice (YOU, 2-month-old), old mice (OLD, 15-month-old), and 2%-betaine-treated old mice (BET, 15-month-old) were employed and raised for 12 weeks. All mice underwent body composition examination and grip strength test before being sacrificed. Betaine alleviated age-related decline in muscle mass and strength. Meanwhile, betaine preserved the expression autophagy markers (Atg5, Atg7, LC3-II, and Beclin1) both at transcriptional and translational level during the aging process. RNA-sequencing results generated from mice gastrocnemius muscle found Mettl21c, a SAM-dependent autophagy-regulating methyltransferase, was significantly higher expressed in BET and YOU group. Results were further validated by qPCR and western bloting. In vitro, C2C12 cells with or without Mettl21c RNA interference were treated different concentration of betaine (0 mM, 10 mM) under methionine-starved condition. Compared with control group, betaine upregulated autophagy markers expression and autophagy flux. By increasing the SAM level, betaine facilitated trimethylation of p97 (Mettl21c downstream effector) into valosin-containing protein (VCP). Increased VCP promoted autophagic turnover of cellular components, ATP production, and cell differentiation. Knock-down of Metthl21c dismissed improvements mentioned above. Collectively, betaine could enhance aged skeletal muscle autophagy level via Mettl21c/p97/VCP axis to delay muscle loss.

T proliferating cells derived autophagy signature associated with prognosis and immunotherapy resistance in a pan-cancer analysis

Despite autophagy modulating tumor immunity in the tumor microenvironment (TME), the immunotherapeutic efficacy and potential mechanism of autophagy signature was not explicit. We manually curated an autophagy gene set and defined a pan-cancer autophagy signature by comparing malignant tissues and normal tissues in The Cancer Genome Atlas (TCGA) cohort. The pan-cancer autophagy signature was derived from T proliferating cells as demonstrated in multiple single-cell RNA sequencing (scRNA-seq) datasets. The pan-cancer autophagy signature could influence the cell-cell interactions in the TME and predict the responsiveness of immune checkpoint inhibitors (ICIs) in the metastatic renal cell carcinoma, non-small cell lung cancer, bladder cancer, and melanoma cohorts. Metabolism inactivation accompanied with dysregulation of autophagy was investigated with transcriptomic and proteomic data. The immunotherapeutic predictive role and mechanism regulation of the autophagy signature was validated in an in-house cohort. Our study provides valuable insights into the mechanisms of ICI resistance.

Asperuloside regulates the proliferation, apoptosis, and differentiation of chronic myeloid leukemia cell line K562 through the RAS/MEK/ERK pathway

Context

Chronic myeloid leukemia (CML) is a malignant hematopoietic stem cell disease caused by excessive proliferation and abnormal differentiation of hematopoietic stem cells. Asperuloside (ASP) is considered to have good biological activity and may be a good anti-CML drug.

Objective

This study aimed to explore the effects and possible mechanisms of ASP on the biological behavior of K562 cells based on RNA-seq.

Materials and methods

The IC50 of ASP in K562 cells was calculated by the concentration-effect curve. Cell viability, apoptosis, and differentiation were detected by CCK8, flow cytometry, benzidine staining, and WB analysis, respectively. Further, RNA-seq was used to analyze the possible mechanism of ASP regulating K562 cells.

Results

ASP significantly inhibited the proliferation, and promoted apoptosis and differentiation of K562 cells. A total of 117 differentially expressed genes were screened by RNA-seq, mainly involved in the RAS/MEK/ERK pathway. PD98059 was used to inhibit the RAS/MEK/ERK pathway in K562 cells, and results confirmed that PD98059 could not only inhibit the RAS/MEK/ERK pathway, but also inhibit the regulation of ASP on the proliferation and differentiation of K562 cells.

Conclusion

ASP inhibited the proliferation, promoted apoptosis and differentiation of K562 cells by regulating the RAS/MEK/ERK pathway, and played a good anti-CML role.

Multi-omics indicators of long-term survival benefits after immune checkpoint inhibitor therapy

Molecular indicators of long-term survival (LTS) in response to immune-checkpoint inhibitor (ICI) treatment have the potential to provide both mechanistic and therapeutic insights. In this study, we construct predictive models of LTS following ICI therapy based on data from 158 clinical trials involving 21,023 patients of 25 cancer types with available 1-year overall survival (OS) rates. We present evidence for the use of 1-year OS rate as a surrogate for LTS. Based on these and corresponding TCGA multi-omics data, total neoantigen, metabolism score, CD8+ T cell, and MHC_score were identified as predictive biomarkers. These were integrated into a Gaussian process regression model that estimates "long-term survival predictive score of immunotherapy" (iLSPS). We found that iLSPS outperformed the predictive capabilities of individual biomarkers and successfully predicted LTS of patient groups with melanoma and lung cancer. Our study explores the feasibility of modeling LTS based on multi-omics indicators and machine-learning methods.

Transcriptome-based analysis of the toxic effects of aluminum chloride exposure on spermatocytes

Aluminum chloride (AlCl3) exposure is pervasive in our daily lives. Numerous studies have demonstrated that exposure to AlCl3 can lead to male reproductive toxicity. However, the precise mechanism of action remains unclear. The objective of this study is to investigate the mechanism of aluminum-induced toxicity by analyzing the alterations in the global transcriptome gene profile of mouse spermatocytes (GC-2spd cells) exposed to AlCl3. GC-2spd cells were exposed to concentrations of 0, 1, 2, and 4 mM AlCl3, and high-throughput mRNA-seq was performed to investigate the changes in the transcriptome after exposure to 4 mM AlCl3. Our findings indicate that exposure to AlCl3 led to an increase in oxidative stress, disrupted glutathione metabolism, reduced cell viability, and altered gene expression in mouse spermatocytes. Gene enrichment analysis revealed that the differentially expressed genes (DEGs) were associated with various biological functions such as mitochondrial inner membrane, response to oxidative stress. Furthermore, these DEGs were found to be enriched in pathways including proteasome, glutathione metabolism, oxidative phosphorylation, and Hif-1 signaling pathway. Real-time PCR and western blot were employed to validate the expression alterations of pivotal genes, and the outcomes exhibited concordance with the mRNA-seq findings. This study provides a theoretical basis for revealing the potential mechanism of male reproductive toxicity caused by aluminum exposure.

Toxicity identification evaluation for hydraulic fracturing flowback and produced water during shale gas exploitation in China: Evidence from tissue residues and gene expression

Hydraulic fracturing flowback and produced water (HF-FPW) from shale gas extraction processes is a highly complex medium with potential threats to the environment. Current research on ecological risks of FPW in China is limited, and the link between major components of FPW and their toxicological effects on freshwater organisms is largely unknown. By integrating chemical and biological analyses, toxicity identification evaluation (TIE) was used to reveal causality between toxicity and contaminants, potentially disentangling the complex toxicological nature of FPW. Here, FPW from different shale gas wells, treated FPW effluent, and a leachate from HF sludge were collected from southwest China, and TIE was applied to obtain a comprehensive toxicity evaluation in freshwater organisms. Our results showed that FPW from the same geographic zone could cause significantly different toxicity. Salinity, solid phase particulates, and organic contaminants were identified as the main contributors to the toxicity of FPW. In addition to water chemistry, internal alkanes, PAHs, and HF additives (e.g., biocides and surfactants) were quantified in exposed embryonic fish by target and non-target tissue analyses. The treated FPW failed to mitigate the toxicity associated with organic contaminants. Transcriptomic results illustrated that organic compounds induced toxicity pathways in FPW-exposed embryonic zebrafish. Similar zebrafish gene ontologies were affected between treated and untreated FPW, again confirming that sewage treatment did not effectively remove organic chemicals from FPW. Thus, zebrafish transcriptome analyses revealed organic toxicant-induced adverse outcome pathways and served as evidence for TIE confirmation in complex mixtures under data-poor scenarios.

4-methylumbelliferone (4-MU) enhances drought tolerance of apple by regulating rhizosphere microbial diversity and root architecture

The dwarfing rootstocks-mediated high-density apple orchard is becoming the main practice management. Currently, dwarfing rootstocks are widely used worldwide, but their shallow root system and drought sensitivity necessitate high irrigation requirements. Here, the root transcriptome and metabolome of dwarfing (M9-T337, a drought-sensitive rootstock) and vigorous rootstocks (Malus sieversii, a drought-tolerant species, is commonly used as a rootstock) showed that a coumarin derivative, 4-Methylumbelliferon (4-MU), was found to accumulate significantly in the roots of vigorous rootstock under drought condition. When exogenous 4-MU was applied to the roots of dwarfing rootstock under drought treatment, the plants displayed increased root biomass, higher root-to-shoot ratio, greater photosynthesis, and elevated water use efficiency. In addition, diversity and structure analysis of the rhizosphere soil microbial community demonstrated that 4-MU treatment increased the relative abundance of putatively beneficial bacteria and fungi. Of these, Pseudomonas, Bacillus, Streptomyces, and Chryseolinea bacterial strains and Acremonium, Trichoderma, and Phoma fungal strains known for root growth, or systemic resistance against drought stress, were significantly accumulated in the roots of dwarfing rootstock after 4-MU treatment under drought stress condition. Taken together, we identified a promising compound-4-MU, as a useful tool, to strengthen the drought tolerance of apple dwarfing rootstock.

Regulation of capsule spine formation in castor

Castor (Ricinus communis L.) is a dicotyledonous oilseed crop that can have either spineless or spiny capsules. Spines are protuberant structures that differ from thorns or prickles. The developmental regulatory mechanisms governing spine formation in castor or other plants have remained largely unknown. Herein, using map-based cloning in 2 independent F2 populations, F2-LYY5/DL01 and F2-LYY9/DL01, we identified the RcMYB106 (myb domain protein 106) transcription factor as a key regulator of capsule spine development in castor. Haplotype analyses demonstrated that either a 4,353-bp deletion in the promoter or a single nucleotide polymorphism leading to a premature stop codon in the RcMYB106 gene could cause the spineless capsule phenotype in castor. Results of our experiments indicated that RcMYB106 might target the downstream gene RcWIN1 (WAX INDUCER1), which encodes an ethylene response factor known to be involved in trichome formation in Arabidopsis (Arabidopsis thaliana) to control capsule spine development in castor. This hypothesis, however, remains to be further tested. Nevertheless, our study reveals a potential molecular regulatory mechanism underlying the spine capsule trait in a nonmodel plant species.

Transcriptional Dysregulations of Seven Non-Differentially Expressed Genes as Biomarkers of Metastatic Colon Cancer

Background: Colon cancer (CC) is common, and the mortality rate greatly increases as the disease progresses to the metastatic stage. Early detection of metastatic colon cancer (mCC) is crucial for reducing the mortality rate. Most previous studies have focused on the top-ranked differentially expressed transcriptomic biomarkers between mCC and primary CC while ignoring non-differentially expressed genes. Results: This study proposed that the complicated inter-feature correlations could be quantitatively formulated as a complementary transcriptomic view. We used a regression model to formulate the correlation between the expression levels of a messenger RNA (mRNA) and its regulatory transcription factors (TFs). The change between the predicted and real expression levels of a query mRNA was defined as the mqTrans value in the given sample, reflecting transcription regulatory changes compared with the model-training samples. A dark biomarker in mCC is defined as an mRNA gene that is non-differentially expressed in mCC but demonstrates mqTrans values significantly associated with mCC. This study detected seven dark biomarkers using 805 samples from three independent datasets. Evidence from the literature supports the role of some of these dark biomarkers. Conclusions: This study presented a complementary high-dimensional analysis procedure for transcriptome-based biomarker investigations with a case study on mCC.

A prognostic and immunotherapeutic predictive model based on the cell-originated characterization of tumor microenvironment in lung adenocarcinoma

Tumor microenvironment (TME) plays a crucial role in predicting prognosis and response to therapy in lung cancer. Our study established a prognostic and immunotherapeutic predictive model, the tumor immune cell score (TICS), by differentiating cell origins in lung adenocarcinoma (LUAD) based on the transcriptomic data of 2,510 patients in 14 independent cohorts, including 12 public datasets and two in-house cohorts. The high TICS was associated with prolonged overall survival (OS), especially in the early-stage LUAD. For the advanced-stage LUAD, high TICS predicted a superior OS in patients who were treated with immunotherapy instead of chemotherapy or TKI. The result suggested that TICS could serve as an indicator for the prognostic stratification management of patients in the early-stage LUAD, and as a potential guide for therapeutic decision-marking in the advanced-stage LUAD. Our findings provided an insight into prognosis stratification and potential guidance for treatment strategy selection.

Whole transcriptome expression profiles in kidney samples from rats with hyperuricaemic nephropathy

Hyperuricaemic nephropathy (HN) is a common clinical complication of hyperuricaemia (HUA) and poses a huge threat to human health. Hence, we aimed to prospectively investigate the dysregulated genes, pathways and networks involved in HN by performing whole transcriptome sequencing using RNA sequencing. Six kidney samples from HN group (n = 3) and a control group (n = 3) were obtained to conduct RNA sequencing. To disclose the relevant signalling pathways, we conducted the analysis of differentially expressed genes (DEGs), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. A competitive endogenous RNA (ceRNA) network was established to reveal the interactions between lncRNAs, circRNAs, mRNAs and miRNAs and investigate the potential mechanisms of HN. Ultimately, 2250 mRNAs, 306 lncRNAs, 5 circRNAs, and 70 miRNAs were determined to be significantly differentially expressed in the HN group relative to the control group. We further authenticated 8 differentially expressed (DE)-ncRNAs by quantitative real-time polymerase chain reaction, and these findings were in accordance with the sequencing results. The analysis results evidently showed that these DE-ncRNAs were significantly enriched in pathways related to inflammatory reaction. In conclusion, HUA may generate abnormal gene expression changes and regulate signalling pathways in kidney samples. Potentially related genes and pathways involved in HN were identified.

Omics technologies in allergy and asthma research: An EAACI position paper

Allergic diseases and asthma are heterogenous chronic inflammatory conditions with several distinct complex endotypes. Both environmental and genetic factors can influence the development and progression of allergy. Complex pathogenetic pathways observed in allergic disorders present a challenge in patient management and successful targeted treatment strategies. The increasing availability of high-throughput omics technologies, such as genomics, epigenomics, transcriptomics, proteomics, and metabolomics allows studying biochemical systems and pathophysiological processes underlying allergic responses. Additionally, omics techniques present clinical applicability by functional identification and validation of biomarkers. Therefore, finding molecules or patterns characteristic for distinct immune-inflammatory endotypes, can subsequently influence its development, progression, and treatment. There is a great potential to further increase the effectiveness of single omics approaches by integrating them with other omics, and nonomics data. Systems biology aims to simultaneously and longitudinally understand multiple layers of a complex and multifactorial disease, such as allergy, or asthma by integrating several, separated data sets and generating a complete molecular profile of the condition. With the use of sophisticated biostatistics and machine learning techniques, these approaches provide in-depth insight into individual biological systems and will allow efficient and customized healthcare approaches, called precision medicine. In this EAACI Position Paper, the Task Force "Omics technologies in allergic research" broadly reviewed current advances and applicability of omics techniques in allergic diseases and asthma research, with a focus on methodology and data analysis, aiming to provide researchers (basic and clinical) with a desk reference in the field. The potential of omics strategies in understanding disease pathophysiology and key tools to reach unmet needs in allergy precision medicine, such as successful patients' stratification, accurate disease prognosis, and prediction of treatment efficacy and successful prevention measures are highlighted.

Comparative Transcriptome Analysis of Two Sugarcane Cultivars in Response to Paclobutrazol Treatment

Sugarcane is an important crop across the globe, and the rapid multiplication of excellent cultivars is an important object of the sugarcane industry. As one of the plant growth regulators, paclobutrazol (PBZ) has been frequently used in the tissue culture of sugarcane seedlings. However, little is known about the molecular mechanisms of response to PBZ in this crop. Here, we performed a comparative transcriptome analysis between sensitive (LC05-136) and non-sensitive (GGZ001) sugarcane cultivars treated by PBZ at three time points (0 d, 10 d, and 30 d) using RNA sequencing (RNA-Seq). The results showed that approximately 70.36 Mb of clean data for each sample were generated and assembled into 239,212 unigenes. A total of 6108 and 4404 differentially expressed genes (DEGs) were identified within the sensitive and non-sensitive sugarcane cultivars, respectively. Among them, DEGs in LC05-136 were most significantly enriched in the photosynthesis and valine, leucine and isoleucine degradation pathways, while in GGZ001, DEGs associated with ion channels and plant-pathogen interaction were mainly observed. Notably, many interesting genes, including those encoding putative regulators, key components of photosynthesis, amino acids degradation and glutamatergic synapse, were identified, revealing their importance in the response of sugarcane to PBZ. Furthermore, the expressions of sixteen selected DEGs were tested by quantitative reverse transcription PCR (RT-qPCR), confirming the reliability of the RNA-seq data used in this study. These results provide valuable information regarding the transcriptome changes in sugarcane treated by PBZ and provide an insight into understanding the molecular mechanisms underlying the resistance to PBZ in sugarcane.

Integration of RNA-seq and ATAC-seq identifies muscle-regulated hub genes in cattle

As the main product of livestock, muscle itself plays an irreplaceable role in maintaining animal body movement and regulating metabolism. Therefore, it is of great significance to explore its growth, development and regeneration to improve the meat yield and quality of livestock. In this study, we attempted to use RNA-seq and ATAC-seq techniques to identify differentially expressed genes (DEGs) specifically expressed in bovine skeletal muscle as potential candidates for studying the regulatory mechanisms of muscle development. Microarray data from 8 tissue samples were selected from the GEO database for analysis. First, we obtained gene modules related to each tissue through WGCNA analysis. Through Gene Ontology (GO) functional annotation, the module of lightyellow (ME_lightyellow) was closely related to muscle development, and 213 hub genes were screened as follow-up research targets. Further, the difference analysis showed that, except for PREB, all other candidate hub genes were up-regulated (muscle group vs. other-group). ATAC-seq analysis showed that muscle-specific accessible chromatin regions were mainly located in promoter of genes related to muscle structure development (GO:0061061), muscle cell development (GO:0055001) and muscle system process (GO:0003012), which were involved in cAMP, CGMP-PKG, MAPK, and other signaling pathways. Next, we integrated the results of RNA-seq and ATAC-seq analysis, and 54 of the 212 candidate hub genes were identified as key regulatory genes in skeletal muscle development. Finally, through motif analysis, 22 of the 54 key genes were found to be potential target genes of transcription factor MEF2C. Including CAPN3, ACTN2, MB, MYOM3, SRL, CKM, ALPK3, MAP3K20, UBE2G1, NEURL2, CAND2, DOT1L, HRC, MAMSTR, FSD2, LRRC2, LSMEM1, SLC29A2, FHL3, KLHL41, ATXN7L2, and PDRG1. This provides a potential reference for studying the molecular mechanism of skeletal muscle development in mammals.

Buffalo long non-coding RNA gene11007 promotes myoblasts proliferation

Buffalo meat is of good quality because it is lean and tender, and could bring significant cardiovascular benefits. The underlying difference in muscle development and meat quality is a complex and precisely orchestrated process which has been demonstrated to be regulated by long non-coding RNAs (lncRNAs). However, the regulatory role of lncRNAs in the growth and development of buffalo skeletal muscle is still unclear. In this study, the Ribo-Zero RNA-Seq method was used to explore the lncRNA expression profiles of buffalo myoblasts during the proliferation and differentiation phases. A specific set of 9,978 lncRNAs was found. By comparing the expression profiles of lncRNAs, it was found that there were 1,576 differentially expressed lncRNAs (DELs) during buffalo myoblast differentiation. Twelve DELs were chosen and subsequently verified in eight different buffalo tissues during fetal and adult stages by using qPCR. Gene11007 was found to be one of the most down-regulated lncRNAs during buffalo myoblasts differentiation and it was subsequently characterized. EdU, CCK-8, qPCR and western blotting assays showed that gene11007 promoted the proliferation of buffalo myoblasts but it had no effect on cell differentiation. Our research may enrich the genome annotations of buffalo and provide a new molecular target for the in-depth understanding of the regulation of lncRNAs in skeletal muscle.

Liver transcriptome analysis reveals biological pathways and transcription factors in response to high ammonia exposure

Aim: Ammonia is a toxic gas that not only causes environmental pollution, but also is harmful to human health after inhalation. Liver is an important detoxification organ that can convert external or metabolized toxic substances into nontoxic substances. However, the toxic effects of ammonia exposure on livers have not been well studied.Method: In this study, pigs were used as an animal model and were exposed to 80 ppm ammonia (8 h during 12 days), and then, RNA-seq were conducted to explore the key genes in response to high ammonia exposure in livers.Result: Gene set enrichment analysis (GSEA) showed that the genes associated with hypoxia, inflammatory response, and apoptosis were up-regulated in the ammonia group, but the genes associated with DNA replication, linoleic acid metabolism, and glycolysis were down-regulated. Totally, 556 differentially expressed genes (DEGs) including 54 genes that encode the transcription factors (TFs) were identified between the exposure and control groups. GO and KEGG pathway analysis suggested that these DEGs were involved in inflammatory response, oxidative stress, apoptosis, immune, and cell cycle. Furthermore, the TF-target interaction analysis showed that FOS, HIF-1α, JUNB, ATF3, REL, and KLF4 were important TFs in regulating the hepatic gene expression in response to high ammonia exposure.Conclusion: Altogether, our findings not only presented a comprehensive mRNA transcriptome profile of liver after high ammonia exposure, but also found some key genes and TFs that could be used to investigate the toxicity mechanism of high ammonia on livers.

Systems neuroimmunology: a review of multiomics methodologies to characterize neuroimmunological interactions in spinal and cranial diseases

Neuroimmunology plays a critical role in our understanding of the pathophysiological processes that underlie a variety of diseases treated by neurosurgeons, including degenerative disc disease (DDD), glioblastoma (GBM), aneurysmal subarachnoid hemorrhage (aSAH), and others. Compared with traditional methods in neuroimmunology, which study one pathway or gene at a time, emerging multiomics methodologies allow for holistic interrogation of multiple immune-signaling pathways to test hypotheses and the effects of therapeutics at a systems level. In this review, the authors summarize key concepts for gathering and analyzing multiomics data so that neurosurgeons can contribute to the emerging field of systems neuroimmunology. Additionally, they describe 3 use cases, based on original research published by their group and others, that utilize transcriptomic, metabolomic, and proteomic analyses to study immune-signaling pathways in DDD, aSAH, and GBM. Through these use cases, techniques for performing machine learning and network-based analyses to generate new clinical insights from multiomics data are shared. The authors hope that neurosurgeons might use this review as a summary of common tools and principles in systems immunology to better engage in creating the immunotherapies of tomorrow.

CBX3 is a Prognostic Biomarker Correlated with ATR Activation and Immune Infiltration in Head and Neck Squamous Cell Carcinoma

Background

Chromobox protein homolog (CBX) family members play important roles in the progression and prognosis of many cancers. However, their functional role in head and neck squamous cell carcinoma (HNSCC) remains largely unknown.

Methods

In this study, we analyzed the expression and functions of CBX family members using The Cancer Genome Atlas data. Most CBX family members were found to be differentially expressed in various tumors, including HNSCC, compared to normal tissues. Multivariate Cox regression analysis showed that CBX3 expression is an independent prognostic factor for HNSCC patients. A nomogram based on CBX3 expression was constructed for use as a diagnostic indicator for HNSCC patients. We also used qPCR to validate the expression of CBX3.

Results

Gene set enrichment analysis suggested that CBX3 participates in ataxia-telangiectasia mutated and Rad3-related protein kinase (ATR) activation and tumor progression. Analysis of immune infiltration indicated that CBX3 expression is negatively correlated with mast cells, DCs, immature DCs, and neutrophils.

Conclusion

Our findings show that high CBX3 expression predicts poor prognosis in HNSCC and that CBX3 may act as an oncoprotein by activating ATR and affecting immune infiltration.

The Mesorhizobium huakuii transcriptional regulator AbiEi plays a critical role in nodulation and is important for bacterial stress response

Background

Bacterial abortive infection (Abi) systems are type IV toxin–antitoxin (TA) system, which could elicit programmed cell death and constitute a native survival strategy of pathogenic bacteria under various stress conditions. However, no rhizobial AbiE family TA system has been reported so far. Here, a M. huakuii AbiE TA system was identified and characterized.

Results

A mutation in M. huakuii abiEi gene, encoding an adjacent GntR-type transcriptional regulator, was generated by homologous recombination. The abiEi mutant strain grew less well in rich TY medium, and displayed increased antioxidative capacity and enhanced gentamicin resistance, indicating the abiEi operon was negatively regulated by the antitoxin AbiEi in response to the oxidative stress and a particular antibiotic. The mRNA expression of abiEi gene was significantly up-regulated during Astragalus sinicus nodule development. The abiEi mutant was severely impaired in its competitive ability in rhizosphere colonization, and was defective in nodulation with 97% reduction in nitrogen-fixing capacity. The mutant infected nodule cells contained vacuolation and a small number of abnormal bacteroids with senescence character. RNA-seq experiment revealed it had 5 up-regulated and 111 down-regulated genes relative to wild type. Of these down-regulated genes, 21 are related to symbiosis nitrogen fixation and nitrogen mechanism, 16 are involved in the electron transport chain and antioxidant responses, and 12 belong to type VI secretion system (T6SS).

Conclusions

M. huakuii AbiEi behaves as a key transcriptional regulator mediating root nodule symbiosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02304-0.

Whole exome and transcriptome sequencing reveal clonal evolution and exhibit immune-related features in metastatic colorectal tumors

Liver is the most common site where metastatic lesions of colorectal cancer (CRC) arise. Although researches have shown mutations in driver genes, copy number variations (CNV) and alterations in relevant signaling pathways promoted the tumor evolution and immune escape during colorectal liver metastasis (CLM), the underlying mechanism remains largely elusive. Tumor and matched metastatic tissues were collected from 16 patients diagnosed with colorectal cancer and subjected to whole-exome sequencing (WES) and RNA sequencing (RNA-seq) for studying colorectal cancer clonal evolution and immune escape during CLM. Shared somatic mutations between primary and metastatic tissues with a commonly observed subclonal-clonal (S-C) changing pattern indicated a common clonal origin between two lesions. The recurrent mutations with S-C changing pattern included those in KRAS, SYNE1, CACNA1H, PCLO, FBXL2, and DNAH11. The main CNV events underwent clonal-clonal evolution (20q amplification (amp), 17p deletion (del), 18q del and 8p del), subclonal-clonal evolution (8q amp, 13q amp, 8p del) and metastasis-specific evolution (8q amp) during the process of CLM. In addition, we revealed a potential mechanism of tumor cell immune escape by analyzing human leukocytes antigens (HLA) related clonal neoantigens and immune cell components in CLM. Our study proposed a novel liver metastasis-related evolutionary process in colorectal cancer and emphasized the theory of neo-immune escape in colorectal liver metastasis.

The Potential Genes Mediate the Pathogenicity of Allogeneic CD4+T Cell in aGVHD Mouse Model

Acute graft-versus-host disease (aGVHD) remains a significant and severe complication of allogeneic hematopoietic stem cell transplantation (allo-HSCT). Due to the occurrence of aGVHD, allo-HSCT significantly increases the mortality rate compared with autologous hematopoietic stem cell transplantation (auto-HSCT). In this study, auto-HSCT and allo-HSCT aGVHD mouse models were built to detect the difference in CD4⁺ lymphocyte in different tissues based on ribonucleic acid sequencing (RNA-Seq) analysis. Clustering analysis, functional annotation, and pathway enrichment analysis were performed on differentially expressed genes (DEGs). The protein-protein interaction (PPI) network was used to find hub genes. CD4⁺T cells were activated by MLR and cytokine stimulation. Cells were sorted out by a flow cell sorter. The selected genes were verified by qRT-PCR, histology, and immunofluorescence staining. The GSE126518 GEO dataset was used to verify the hub genes. Enrichment analysis revealed four immune-related pathways that play an important role in aGVHD, including immunoregulatory interactions between a lymphoid and a nonlymphoid cell, chemokine receptors binding chemokines, cytokine and cytokine receptor interaction, and the chemokine signaling pathway. At the same time, with the PPI network, 11 novel hub genes that were most likely to participate in immunoregulation in aGVHD were identified, which were further validated by qRT-PCR and the GSE126518 dataset. Besides, the protein expression level of Cxcl7 was consistent with the sequencing results. In summary, this study revealed that immunoregulation-related DEGs and pathways played a vital role in the onset of aGVHD. These findings may provide some new clues for probing the pathogenesis and treatment of aGVHD.

Transcriptomics and Proteomics Methods for Xenopus Embryos and Tissues

The general field of quantitative biology has advanced significantly on the back of recent improvements in both sequencing technology and proteomics methods. The development of high-throughput, short-read sequencing has revolutionized RNA-based expression studies, while improvements in proteomics methods have enabled quantitative studies to attain better resolution. Here we introduce methods to undertake global analyses of gene expression through RNA and protein quantification in Xenopus embryos and tissues.